流场结构对PEM电解槽性能影响模拟_王华.pdf

第 卷第 期洁 净 煤 技 术.年 月 流场结构对 电解槽性能影响模拟王 华，马晓锋，何 勇，徐超群，朱燕群，王智化（浙江富兴海运有限公司，浙江 杭州；浙江大学 能源清洁利用国家重点实验室，浙江 杭州）移动阅读收稿日期：；责任编辑：白娅娜 ：基金项目：国家自然科学基金资助项目（）作者简介：王 华（），男，浙江杭州人，工程师，硕士。：通讯作者：王智化（），男，河南荥阳人，教授，博士。：引用格式：王华，马晓锋，何勇，等 流场结构对 电解槽性能影响模拟洁净煤技术，（）：，（）：摘 要：为了提高质子交换膜（）电解槽的性能，降低电压损耗，提升运行稳定性，在一定假设条件下，使用 软件建立了全尺寸多通道的质子交换膜电解槽三维模型，并将仿真计算结果与文献中同尺寸 电解槽试验测试结果进行对比验证，仿真结果与文献试验测试结果基本吻合。基于此电解槽模型研究了流道高度、堵块、脊宽度等对 电解槽性能的影响，从而优化流场结构。仿真结果显示，在研究设定的 电解槽尺寸下，最佳流道高度为 ；在流道顶部设置堵块可使扩散层内氧气质量分数下降约，膜电极平均温度下降 左右，电解槽的电解电势减少 左右；流场脊宽度由 减小至 时，平均氧气质量分数下降约，膜电极平均温度下降，电解槽电解电势下降 左右。优化后的流场结构有利于循环水带走扩散层内氧气，降低氧气气泡堵塞扩散层孔隙的可能性，减小传质阻力，增强传热过程，及时排出电解槽中多余热量，降低电解槽电解电势，提高电解槽性能和运行稳定性。关键词：质子交换膜电解槽；流场结构；电解电势；氧气质量分数；膜电极温度中图分类号：文献标志码：文章编号：（），（，；，）：（），：；王 华等：流场结构对 电解槽性能影响模拟 年第 期 引 言稳定持续的能源供给是人类社会发展的重要前提。据国际能源署（）预测，世界能源供应需求仍将继续增长，到 年，全球能源需求总量将达 亿 石油当量。而在双碳时代背景下，可再生清洁能源的开发和利用尤为重要。年 月 日，国务院发布了关于完整准确全面贯彻新发展理念做好碳达峰碳中和工作的意见，标志着我国坚定不移地走生态优先、绿色低碳的高质量发展道路，确保如期实现碳达峰、碳中和目标。氢能清洁、高效，且反应产物是水，无 排放，具有较高质量能量密度（），是未来最理想的能源形式之一，发展氢能经济是实现双碳目标的必由之路。氢气在自然界中并不存在，需通过输入其他一次能源生产。迄今为止，氢气主要通过天然气重整、煤气化、电解水等方式制备，然而化石燃料蒸汽重整产生的低纯度氢中还含有高浓度含碳物种，如、等。且蒸汽重整不会减轻对化石燃料的依赖，也不会减少污染物，对实现双碳目标作用甚微。氢能作为一种储能方式，可解决可再生能源难以持续供应的问题。随着可再生能源产业的大力发展和水电解制氢技术不断进步，水电解制氢必将成为未来制氢产业的发展趋势。国际能源署预测，到 年电解水制氢占比增至 左右。而质子交换膜（）电解水制氢技术具有电流密度高、运行负荷范围宽、动态响应速度快等优点，可直接耦合可再生能源实现绿氢生产，极具应用前景。经过 多年发展，水电解制氢技术已从实验室研究阶段过渡到产业化初始阶段。尽管如此，质子交换膜电解槽还存在成本高、耐久性不足等问题。如何降低 电解槽的电压损耗、维持电解槽长期稳定运行是质子交换膜电解水技术研究的首要任务。在高电流密度下，大量水在阳极催化层发生电解产生氧气，一段时间后电解槽的阳极气体扩散层内会积聚大量氧气，气体扩散层孔隙被氧气气泡堵塞，未反应的水无法穿过扩散层到达催化层发生连续反应，增加了电解槽内部的传质阻力，进而降低电解效率。由于氧气导热系数远小于水，若不及时排出产生的氧气，会影响电解槽散热，造成电解槽局部过热，严重时会导致质子交换膜发生破裂，存在安全隐患，。通过改善流场设计，有望强化传质，及时排出电解水过程产生的氧气，减少传质阻力，从而降低电压损耗。同时优化流场结构还可强化传热，带走热量，提升电解槽运行稳定性。质子交换膜电解水制氢过程涉及电化学反应、多组分流体流动，固体流体传热、多孔介质传热。但在试验测试中仅能通过电化学工作站监测电解槽运行电流和电压，通过在线氢气分析仪在线测量电解槽氢气产量，但无法监测电解槽运行时内部的物质分布、温度分布等。因此，通过对 水电解槽进行合理化建模，借助数值仿真软件实现电解槽内部物质分布和热量分布的仿真，对电解槽运行时内部传质传热情况进行分析，为优化流道结构设计提供理论依据。等研究了单通道 电解池的阳极流道深度对电解槽电解性能的影响。在低、中电流密度下，流道深度对电解槽性能的影响可忽略不计，但在高电流密度下影响显著。何旭等借助 方法模拟研究了单通道下纯水流速、电解电流密度对电解产生氧气在流动通道中运输特性的影响。增大纯水流速也有利于氧气泡脱离，减小气泡体积；较高的电解电流密度导致氧气气泡生成速率较高，增大多孔气体扩散层表面的气体覆盖率，阻碍电解水反应连续发生。以往针对质子交换膜电解槽流道结构研究大多围绕单通道 电解槽进行。笔者基于 仿真软件，建立一个全尺寸多通道 电解槽三维模型，研究流道高度、堵块、脊宽度等流场结构设计对 电解槽电解性能和内部传质传热现象的影响，通过比较相同工作电流密度下阳极气体扩散层内平均氧气质量分数、膜电极平均温度和电解电势，为优化电解槽流场结构设计提供理论依据。模型建立及验证.电解水制氢技术 电解槽装置的基本结构如图 所示，主要由膜电极、阴阳极扩散层、阴阳极流道以及阴阳极极板、端板组合而成。电解槽工作原理为：水进入阳极端板，经流道后穿过扩散层，在一定电压下，在阳极催化层发生析氧反应产生氧气和氢离子，氧气与未经反应的水一起排出；氢离子穿过质子交换膜，在阴极催化层得到电子发生析氢反应，产生氢气，穿过扩散层经流道排出。电解槽反应方程式为总反应：，（）阳极反应：，（）年第 期洁 净 煤 技 术第 卷阴极反应：。（）图 电解槽结构示意 .模型假设建立一个同时涉及电化学反应、两相流体流动及传热、多孔介质传热的全尺寸的多通道 电解槽三维模型（图），做以下假设：参与电解反应的纯水及产生的氧气和氢气均是不可压缩流体；只考虑阳极侧的纯水和氧气混合流体流动，只在阳极侧供水，阴极侧的氢气产生后直接排出；忽略模型界面间的接触电阻；认为质子交换膜只能透过氢离子，不存在氧气和氢气的交叉扩散现象；模型求解为带电流初始化的稳态结果。图 质子交换膜电解槽模型剖面 .模型建立.电化学模型 电解槽的电解电势 主要由 部分组成（式），分别为开路电压、活化过电势、欧姆过电势：。（）在不同温度和压力下开路电压有所差别，可通过能斯特方程计算得到：，（）（）（）（）|。（）其中，为一般气体常数，（）；为法拉第常数，；为电解槽运行时的温度，；（）为、的平衡压力，；为电子转移数；，（）可通过式（）计算：，（）（）。（）活化过电势分为阴阳两极的活化过电势（，和，），通过 方程计算活化过电势：，|（），|，（），|（），|，（）式中，为电化学活性比表面积，；、分别为阳极和阴极的电荷转移系数；，、，分别为阳极和阴极交换电流密度，。电子在电解槽内部转移时，会有欧姆损失，主要由质子交换膜、阴阳极催化剂层、阴阳极气体扩散层、阴阳极流场板的电阻产生。将这几部分组件假设成各向同性材料，且只考虑电解槽运行温度对质子交换膜电阻的影响。质子交换膜电导率 为|。（）欧姆过电势 计算公式为|，（）式中，、分别为质子交换膜和电解池固相组件的电导率，；为电流，；为电解质电势，；为电子电势，。.流体流动模型通过质量守恒定律和动量守恒定律方程对 电解槽内部的两相流体流动进行描述。假设流体为不可压缩，且包含了多孔介质流体流动，通过 方程描述电解池内部各组分的对流和扩散。.热量传递模型 电解槽的热量来源主要有水流带来的热量和电流做功产生的热量，热量散失途径主要为对流散热、辐射散热及电解产生的气体和未反应的水流出电解池带走的热量。热量分布主要通过能量守恒方程描述。.模型验证模型边界条件设置如下：初始状态下，流道内充满液态水，阳极和阴极出口均为大气压，在阳极流道板侧添加电势，因不考虑阴极侧的气体流动，阴极侧的扩散层接地。模型参数设定见表。王 华等：流场结构对 电解槽性能影响模拟 年第 期表 模型参数，参数数值阳极传递系数 阴极传递系数 扩散层电导率 （）扩散层孔隙率 质子交换膜热导率 （）氧气热导率（）（）氢气热导率（）（）扩散层热导率 （）电解槽的尺寸参数设定见表。表 尺寸参数 项目数值极板边长 极板厚度 流道宽度 流道高度 单脊宽度 质子交换膜厚度 气体扩散层厚度 对模型进行网格加密后计算，将仿真结果与单通道 电解槽模型及文献中同尺寸 电解槽数据进行比较，如图 所示。全尺寸模型计算的电化学性能与 等测得数据更接近，平均工作电流密度为 时相对误差仅，而单通道模型为，说明全尺寸多通道模型更完善，模型建立基本正确。图 全尺寸模型仿真结果与单通道模型、文献数据对比 流场结构优化液态水进入流道后，在流场板分配作用下穿过阳极扩散层后在催化层发生分解，电解产生的氧气穿过扩散层后与未参加反应的水一起流出流道后排出电解槽。基于建立的全尺寸多通道三维模型进行不同流场结构的仿真计算，研究流道高度、堵块、脊宽度（图）对 电解槽性能的影响。图 质子交换膜电解槽三维模型 .直行流道高度的影响控制入口流速相同，探究流道高度对 电解槽性能的影响。流道高度设定为、五种尺寸，相应的入口水流量设定为、，电解槽入口水流初始温度为 ，电解槽平均工作电流密度为 。电解槽模型仿真计算结果见表。表 不同流道高度的 电解槽模型仿真结果 流道高度扩散层平均氧气质量分数 膜电极平均温度 电势 仿真结果表明，在研究设定的 电解槽尺寸下，流道高度为 时，高电流密度下稳态运行时阳极气体扩散层内的平均氧气质量分数最小为，膜电极平均温度最低为 ，达 电流密度所需电解电势最低为 ，此时电解槽性能最佳，电压损耗最小。随流道高度增加，扩散层内平均氧气质量分数、膜电极平均温度、电解电势均增大。流道高度越小，越有利于及时排出电解产生的氧气，扩散层孔隙中氧气气泡越少则越有利于水穿过扩散层在催化层发生持续反应，同时水含量增大有利于水流带走多余热量，增加导热能力，提高电解槽性能和运行稳定性。膜电极平均温度与扩散层内平均氧气质量分数呈正相关，这主要是由于水的导热系数远大于氧气，扩散层内氧气质量分数越低，水含量越高，可明显增强电解槽散 年第 期洁 净 煤 技 术第 卷热，降低电解槽整体温度，从而提升电解槽运行稳定性。但流道高度减小至 时，电解电势、膜电极的平均温度、阳极扩散层的平均氧气质量分数均上升，这是由于流道高度过低时，水流量较小，水流能带走的氧气和热量有限。控制入口水流速度不变时，电解槽阳极流场的流道高度会影响电解槽内部氧气和水的物质传递过程及传热过程，从而影响电解槽运行时所需电解电势。适当减小流道高度，可增强传质传热过程，提升电解槽性能和稳定运行。该电解槽尺寸下最佳流道高度为 。.不同形状堵块的影响扩散层内氧气传质过程对电解槽性能有重要影响，增强电解槽内部流体扰动有助于水流带出氧气。在直行流道顶部设置堵块，可改变水在流道内的运动方向，促进水流向气体扩散层，如图 所示。图 扩散层内水和氧气流动示意 建立了 种不同流道结构的 电解槽模型，分别在流道顶部设置半圆形、三角形和矩形堵块，如图 所示。图 设置不同形状的堵块的直行流道 直行流道高度设置为最佳高度 ，堵块高度为 ，为流道高度的，入口水流量设为，入口水温为 ，膜电极平均电流密度为 。电解槽的模型仿真计算结果见表。由表 可知，最佳流道高度下，在直行流道顶部设置不同形状堵块后，扩散层内平均氧气质量分数最大下降约 个百分点，膜电极平均温度最多降低 ，电势下降 左右。证明在流道顶表 不同形状堵块的 电解槽模型仿真结果 堵块形状扩散层平均氧气质量分数 膜电极平均温度 电势 无堵块 半圆形 三角形 矩形 部设置堵块可增强流体扰动，改变水流流动方向，增加水流流向扩散层的速度，克服传质阻力，水流更易穿过扩散层到达催化层持续发生析氧反应；同时有利于带走热量并增强散热，降低膜电极温度，降低电解电势且提升电解槽稳定性。另外，流道顶部设置不同形状堵块对电解槽性能的影响差别不大，其中设置矩形堵块优化效果最佳。.脊宽度的影响适当增加流体扰动可一定程度上提高电解槽电化学性能，但需考虑到在流道末端增加扰动会导致原本已穿过扩散层进入流道的氧气再次随水流进入扩散层，影响氧气排出，降低电解槽性能。电解槽稳态运行时阳极气体扩散层内氧气质量分数分布热力图如图 所示，可知脊对应位置的扩散层内氧气质量分数明显高于平均值，特别是流场两侧边界脊对应的扩散层处，氧气质量分数高达，这主要是由于两侧边界处的扩散层不仅无法与水直接接触，且不同于中间处脊对应的扩散层，氧气可被两边流道内水流带走，所以两侧边界扩散层内积聚了大量氧气，导致电解水反应缓慢。图 脊宽 的电解槽扩散层内氧气质量分数分布热力图 因此在不改变流场流道水流总面积的情况下，王 华等：流场结构对 电解槽性能影响模拟 年第 期研究减小单个脊宽度对 电解槽性能的影响。不同宽度脊的直行流道结构如图 所示。图 不同宽度脊的直行流道 设定水流量为 ，入口水温度为，电解槽稳态运行时设定膜电极的平均电流密度为 ，流道深度为 ，建立 种不同宽度脊直行流道的 电解槽模型进行模拟，仿真计算结果见表。表 不同宽度脊的 电解槽模型仿真结果 脊宽度扩散层平均氧气质量分数 膜电极平均温度 电势 仿真结果表明，随单个脊宽度减小，阳极气体扩散层内平均质量分数、膜电极平均温度、电解槽电解电势均明显下降。电解槽脊宽度从 减小至 时，扩散层内平均氧气质量分数下降约 个百分点，膜电极平均温度下降 ，电解槽的电势减少 左右。说明减小单个脊的宽度，有利于提高 电解槽的电化学性能，达到相同电流密度所需电势更低。种不同宽度脊的 电解槽模型稳态运行时阳极气体扩散层内氧气质量分数分布的最大值和最小值见表。可知随脊宽度减小，扩散层中氧气质量分数最大值下降，由 降至 时，氧气质量分数最大值由 降至，降幅明显，氧气质量分数最小值变化不明显，说明减小单个脊宽度可促进脊对应处扩散层内氧气被水流带出，减少氧气在脊处扩散层内积聚的可能性，从而增强电解槽内部传质传热过程，降低电解电势，提升电解槽稳定性。电解槽脊宽度为 时，电解槽性能较好。尽管减小脊宽度有利于提高电解槽性能和稳定性，但减小脊宽度还应考虑流场结构的加工成本和技术可行性。表 扩散层中氧气质量分数极值 脊宽度氧气质量分数最大值 氧气质量分数最小值.综合影响基于上述结果，设置流道高度为 ，脊宽度为 ，在流道顶部设置矩形堵块，仿真结果显示，电 解 槽 扩 散 层 内 平 均 氧 气 质 量 分 数 为，膜电极平均温度为 ，电解槽的电解电势为 ，说明电解槽的电解性能进一步提升。结 论）在该研究设定的 电解槽模型尺寸下，最佳流道高度为 。）在流道顶部设置不同形状的堵块可提升电解槽性能，电解槽扩散层内的氧气质量分数下降约 个百分点，膜电极平均温度下降 左右，电解槽电解电势降低 左右，矩形堵块对电解槽的性能提升最大。）脊宽度从 减小至 时可使扩散层内平均氧气质量分数下降 个百分点，膜电极平均温度下降 ，电解槽的电势降低 。）结合最佳流道高度、堵块形状和脊宽度可使扩散层内平均氧气质量分数下降 个百分点，膜电极平均温度下降 ，电解电势降低 。参考文献（）：，（），（）：，（）：，：，（）：，：，年第 期洁 净 煤 技 术第 卷 ，（）：，（）：，（）：，：，（）：，（）：，（）：，（）：何旭，罗马吉，陈奔 质子交换膜电解池内氧气泡输运过程特性 江苏大学学报（自然科学版），（）：，（），（）：，（）：，（）：，：，（）：，（）：，（）：，（）：，：，：，：，：，：，（）：，：