2. 长安大学 能源与电气工程学院, 陕西 西安 710018
为了实现交通领域碳中和目标, 发展氢能源燃料电池汽车十分必要。近年来, 我国先后出台了相关政策, 有力促进了以质子交换膜燃料电池汽车的快速发展[1-2]。相比于纯电动汽车, 质子交换膜燃料电池汽车以氢为“燃料”, 具有续航足、加氢快、寿命长等优势, 近年来受到各国的青睐[3]。
PEMFC作为质子交换膜燃料电池汽车的核心动力部件, 具有能量转换效率高、排放零污染、耐久性强等优点[4-5]。工作温度与进气湿度对PEMFC性能产生重要影响, 不当的工作温度和进气湿度会造成“阳极拖干”和“阴极水淹”等问题[6]。
为解决上述问题, 国内外学者对工作温度和进气湿度开展了相关研究。樊林浩等[7]通过构建三维模型发现: 在阳极低加湿的工况下, 逆流进气比顺流进气更有利于提升电池性能。Xing等[8]研究发现适当提升阴极侧的工作温度, 会提高PEM的电导率、提升还原反应动力学速率, 从而提升电池性能。陆佳斌等[9]采用实验与仿真相结合的手段, 研究了阴极湿度对电池性能的影响,结果表明:增大电流密度、降低阴极湿度, 可提升电池性能; 此外电池反应产生的水会优先补充到PEM中, 当PEM的含水量达到工作阈值后, 反应生成的水会进入扩散层, 使扩散层内液态水饱和度增加。蒋杨等[10]构建了一维PEMFC模型, 发现采用阳极充分加湿或阴极低加湿的方式, 可有效降低“阳极拖干”和“阴极水淹”的风险, 改善电池内部受热情况。但是, 在实际运行中, PEMFC处于多参数动态变化过程。目前研究大多集中于分析单一参数变化对PEMFC性能的影响[11-16], 关于多参数变化对PEMFC性能的影响研究相对较少。
本文通过耦合分析工作温度与阴极湿度2个参数对不同电位下电池输出性能和水分布的影响, 得到在不同电压阶段下, 电池对温湿度的敏感性影响, 为电池的实际运行策略提供参考。
1 物理模型和几何参数 1.1 PEMFC物理模型本文以逆流形式[17]的单直通道PEMFC为研究对象, 探究工作温度与阴极湿度对电池性能的影响。采用Nafion112质子交换膜, 催化剂为铂金属, 催化层铂碳质量比为20%, 扩散层孔材料为碳。图 1为单体PEMFC的几何模型, 具体几何参数见表 1。
![]() |
图 1 PEMFC几何模型 |
组成部件 | 长/mm | 宽/mm | 高/mm |
质子交换膜 | 30 | 2 | 0.03 |
阴/阳极催化层 | 30 | 2 | 0.01 |
阴/阳极扩散层 | 30 | 2 | 0.20 |
阴/阳极流道 | 30 | 1 | 1.00 |
阴/阳极极板 | 30 | 2 | 1.50 |
为防止出现“阳极拖干”, 阳极进气湿度假定为100%[10]。在保证各项传输与反应合理的前提下, 对数值模型做出如下假设:
1) 工作压力为常压, 混合气体遵守理想气体状态方程[18];
2) 在稳态条件下运行;
3) 由于电池内部流速较低, 为层流流动[19];
4) 气体扩散层和催化剂层均为各向同性、均质、不可压缩的多孔材料[20];
5) PEM不允许气体通过[21];
6) 不考虑流道内液态水传输。
基于上述假设, 建立PEMFC的三维数学模型。
![]() |
(1) |
式中: ρ为混合气体密度, kg/m3; ε为多孔介质孔隙率; u为速度场, m/s; Sm为质量源项, kg/(m3·s)。
![]() |
(2) |
式中:μ为混合物的动力黏度, kg/(m·s); Su为动量源项, kg·m/s。
能量守恒方程[21]
![]() |
(3) |
式中:cp为气体比热容, J/(kg·K); T为温度, K; keff为有效导热系数, W/(m·K); SQ为能量源项, W/m3。
![]() |
(4) |
式中, 下标k为气体组分H2、O2、H2O和N2; ck为组分浓度, mol/m3; Dkeff为有效扩散系数; Sk为组分源项, mol/(m3·s)。
电荷守恒方程[8]
![]() |
(5) |
![]() |
(6) |
式中: σseff和σmeff为固相和膜相有效电导率, S/m; ϕs和ϕm分别为固相和膜相电势, V; is和im分别为电子流和质子流源项, A/m3。
电化学动力学守恒方程[19]
![]() |
(7) |
![]() |
(8) |
式中: a为单位体积内电催化表面积, m-1; iaref为参考电流密度, A/m3; β为传递系数; η为过电位, V; crefH2和crefO2为氢气和氧气参考浓度, mol/m3; s为液态水饱和度。
液态水传输方程[8]
![]() |
(9) |
式中: ρwl为液态水密度, kg/m3; Dc为毛细管扩散系数; krl和kgr分别为液态水和气态水相对渗透率, m2; μwl和μgw分别为液态水和气态水的动力黏度, kg/(m·s); Mw为水的分子质量, kg/mol; Swl为液态水向气态水及溶解水转化的源项, mol/(m3·s)。
溶解水传输方程[8]
![]() |
(10) |
![]() |
(11) |
![]() |
(12) |
式中: nd为电渗系数; iM为团聚体体积分数(膜内iM=1, 催化层内iM=0.2);cwd为溶解水浓度, mol/m3; DW-M为膜中水扩散系数; kp, M为水在膜中的渗透率, m2; Sw2为溶解水源项, mol/(m3·s); λ为PEM含水量; σm为PEM电导率, S/m; Mpem为PEM摩尔质量, kg/mol; ρmem为PEM密度, kg/m3。
有关数值模型的边界条件为:
阴阳极两侧流道入口为气体入口, 其入口质量流量定义为[10]
![]() |
(13) |
![]() |
(14) |
式中: Iref为参考电流密度, A/m2; ρga和ρgc分别代表阳极和阴极混合气密度, kg/m3; Ra和Rc分别为阳极和阴极入口混合气的化学计量比; A为燃料电池的活化面积, m2; cH2, in和cO2, in分别为流道入口氢气和氧气的摩尔浓度, mol/m3, 可以通过(15)~(16)式求得[10]
![]() |
(15) |
![]() |
(16) |
式中: pga和pgc分别为阳极和阴极入口处混合气压力, kPa; RH, a和RH, c分别为阳极与阴极入口湿度; psat为水蒸气的饱和蒸汽压, 是压力和温度的函数; T0为入口处温度, K。此外, 阴阳两极出口压力为常压条件。
上述数值模型的具体参数见表 2。
参数 | 数值 | 参数 | 数值 | |
扩散层孔隙率 | 0.6 | 催化层孔隙率 | 0.5 | |
扩散层渗透率/m2 | 6.2×10-12 | 催化层渗透率/m2 | 6.2×10-13 | |
扩散层接触角/(°) | 100 | 催化层接触角/(°) | 100 | |
扩散层电导率/(S·m-1) | 300 | 催化层电导率/(S·m-1) | 300 | |
H2参考浓度/(mol·m-3) | 40.88 | O2参考浓度/(mol·m-3) | 40.88 | |
阳极催化层比表面积/m-1 | 100 000 | 阴极催化层比表面积/m-1 | 100 000 | |
阳极化学计量 | 1.2 | 阴极化学计量 | 2.0 | |
阳极压力/Pa | 101 325 | 阴极压力/Pa | 101 325 | |
PEM等效摩尔质量/(kg·mol-1) | 1.1 | PEM密度/(kg·m-3) | 1 980 | |
阳极传递系数 | 0.5 | 阴极传递系数 | 0.7 | |
阳极铂载量/(mg·cm-2) | 0.1 | 阴极铂载量/(mg·cm-2) | 0.4 |
由于PEMFC常见工作温度范围为60~80℃, 且阴极进气加湿不宜超过100%。本文选取3种工作温度(60, 70, 80℃)和3种阴极湿度(40%, 70%, 100%)进行正交研究, 具体组合方案见表 3。
为了验证仿真模型的正确性, 将C3组合下的电池极化曲线的仿真与实验结果[22]进行对比。实验中, 使用的质子交换膜型号为Nafion112, 有效面积为1.0 cm×1.0 cm, 阳极铂载量为0.1 mg/cm2, 阴极铂载量为0.4 mg/cm2, 铂碳质量比为20%, 阴阳极化学计量比分别为2.0, 1.2, 工作温度为80℃, 阴阳极湿度均为100%, 工作压力为101.325 kPa。从图 2可以看出, 两者结果吻合良好, 它们之间的积分相对误差为4.07%。并且, 两者之间的最大误差出现在电压为0.45 V时, 此时, 两者之间的相对误差为6.74%。
![]() |
图 2 仿真模型验证 |
表 4给出了0.5 V电压、80℃工作温度、100%阴极湿度条件下, 电池的电流密度。从表中可以看出, 随着网格数的增加, 电流密度区域稳定。当网格数大于78 960, 网格数对计算结果影响较小; 且与网格数为157 920的电流密度值相比, 它们之间的相对误差小于0.002%。综合计算精度和计算效率, 在后续研究中, 选取网格数78 960。
网格数 | 20 160 | 39 480 | 78 960 | 131 600 | 157 920 |
电流密度/(mA·cm-2) | 942.91 | 940.45 | 934.12 | 934.13 | 934.13 |
图 3~4分别给出了C3组合下, 0.5 V电压时, PEMFC内部氧气与氢气的分布图。
![]() |
图 3 氧气摩尔分数分布 |
![]() |
图 4 氢气摩尔分数分布 |
从图 3可以看出, 在燃料电池内部多孔区域中, 氧气沿着流动方向逐渐被消耗, 其浓度分布呈阶梯降低趋势; 在阴极肋板下多孔区域的边界处氧气浓度分布较低, 最小值约为0.18%。从图 4中可以看出, 阳极侧氢气浓度分布较为均匀, 这保证了阳极侧催化层充分的氢氧化反应; 此外, 在100%阳极加湿的条件下, 电池阳极侧水相摩尔分数维持在50%附近, 避免了“阳极拖干”现象。
2.3 PEM电导率与含水量PEM电导率是影响电池欧姆损失的关键因素; 含水量可以反映出PEM内部湿润程度, 对PEM寿命具有重要影响。图 5~6分别表示PEM含水量和电导率与电压的关系曲线。为了更好分析9组电池的性能, 根据电压变化分为高电位(0.9~0.7 V电压)、中电位(0.7~0.5 V电压)和低电位(0.5~0.3 V电压)3个阶段。
![]() |
图 5 PEM含水量与电压关系 |
![]() |
图 6 PEM电导率与电压关系 |
从图 5可以看出, 随着电压的降低, 温度越高PEM含水量越低, 阴极湿度越高PEM含水量越高。在高电位阶段, 阴极湿度对PEM含水量影响最为明显。这是由于在高电位阶段, 电化学反应速率较慢, 电渗拖拽作用并不显著, 此时压差扩散是影响PEM含水量的主要因素。因此, 在高电位阶段, 阴极湿度对PEM含水量影响最显著。在中电位阶段, 随着电压的降低, PEM逐渐湿润, 阴极湿度对PEM含水量的影响逐渐削弱。在低电位阶段, 各组电池的PEM内部含水量逐步趋于稳定。
图 6表明, 温度越高, PEM电导率越大; 阴极湿度越大, PEM电导率越大。在高电位阶段, 与温度相比, 湿度对PEM电导率的影响更显著。而在低电位阶段, 电导率几乎不受湿度影响。这是含水量逐渐饱和导致的结果。
相比单一改变温度或湿度, 同时调节温度和阴极湿度可以使电池在各工作点有更小的欧姆内阻。在高电位阶段, 采用“阴极高加湿”方案, 可以使PEM含水量与电导率保持在较高的值; 在中电位阶段, 由于PEM逐渐湿润, 可以逐步降低阴极加湿度保证多孔域的透气性; 在低电位阶段, 由于PEM含水量已经达到饱和, 此时采用“阴极高温”可以使得PEM电导率更高, 并且同时采用“阴极低加湿”, 可有效缓解阴极多孔区域“水淹”现象。
为了更加直观地描述PEM含水量分布规律, 图 7给出了在0.5 V电压下, 9组PEM电池的含水量。A1电池(60℃、40%阴极湿度)和A2电池(70℃、40%阴极湿度)内部含水量最大差值为0.31和0.59,与之相比, A3电池(80℃、40%阴极湿度)内部含水量的最大差值为0.80。这容易导致PEM内部容易出现局部皱缩甚至脱水的现象, 增加电池的欧姆损耗, 并降低电池的耐久性。结合B3(80℃、70%阴极湿度)和C3(80℃、100%阴极湿度)电池含水量分布发现, 当工作温度较高时, 采用较高的阴极湿度, 可以有效避免PEM局部含水量偏低的现象。
![]() |
图 7 0.5 V电压下PEM含水量 |
在高电位阶段, 电池内部反应速率较慢, 多孔区域发生“水淹”的风险较低。并且, 电池反应产生的水会优先补充到PEM中, 当PEM含水量达到工作阈值后, 反应生成的水会进入扩散层, 增加扩散层液态水的饱和度。因此, 为分析电池内部多孔区域液态水分布对电池性能的影响, 图 8和图 9分别给出了在0.5 V电压下, 阴极扩散层的液态水饱和度分布以及对应的阴极催化层与扩散层交界面处的电流密度。
![]() |
图 8 阴极扩散层液态水饱和度 |
![]() |
图 9 阴极催化层与扩散层交界面电流密度 |
从图 8可以看出, A3电池扩散层流道下方与肋板(极板与多孔区域接触部分)下方的液态水饱和度分布最为均匀, 最大差值仅为0.06, 这有效保障了氧气在整个扩散层的均匀、稳定传质。其余8组电池内部水含量最大值均超过0.10, 在流道下方与肋板下方区域内液态水饱和度存在明显的断层, 扩散层内部流动不稳定, 部分液态水无法被及时送入流道, 导致肋板下方区域氧气的传质阻力增大。
对比图 8和图 9, 可以发现随着阴极湿度的增大, 电池肋板下方扩散层液态水饱和度逐渐提升, 导致这部分区域氧气传输受阻严重, 反应物得不到及时补充, 电化学反应变缓, 阻碍电池性能提升。在9组电池中, A2和A3电池具有更高电流密度。
总的来看, 阴极湿度越低, 扩散层液态水饱和度分布更加均衡。在中低电位阶段, 采用较低的阴极湿度可以有效减少多孔区域液态水含量, 降低“水淹”风险, 利于提升电池电流密度。另外, 温度的提升也会使液态水饱和度下降, 提升电池的输出性能。
2.5 功率密度功率密度可以直观反映燃料电池的输出性能, 图 10给出了9组电池的电压与功率密度变化曲线。从图中可以看出: 在高电位阶段, 9组电池电化学反应速率均较低, 其功率密度相差不大。在中电位阶段, A2电池(70℃工作温度、40%阴极湿度)的功率密度最大。这是由于经过高电位阶段, 电池内部逐渐湿润, PEM含水量与电导率迅速上升, 电池内部反应速率增大, 产水量相应增加。结合图 8可以看出, 高阴极湿度会使电池内部积水过多而堵塞多孔区域, 进而阻碍氧气的传输, 导致电池功率密度下降。此外, 较高的温度会制约阴极入口处氧气的比重, 导致催化层表面缺氧, 同样降低了电池的输出性能。在低电位阶段, 虽然阴极侧反应速率进一步提高, 产水量更多, 但由于电池内部温度逐步上升, 0.5 V时A2、A3、B2和B3电池内部液态水所占的比重更少, 氧气的传质过程不会受到阻碍。此外, 由图 5可知, 电压低于0.5 V时, PEM含水量已经达到饱和, 低电位阶段随着电压的降低, 电池内部“水淹”会加剧。因此, 中低湿度与中高温组合的这4组电池在低电位阶段展现出更好的性能。
![]() |
图 10 电压与功率密度变化图 |
从整个过程来看, 温度维持在60℃, 单一地改变电池阴极湿度无法使电池的功率密度有较大提升; 同样湿度维持在70%, 单一调节电池工作温度, 电池的功率密度也并不是在每个电压下都能达到9组电池中的最大值。因此, 单一地改变温度或湿度对提高电池输出性能具有一定的局限性。同时调节温度和阴极湿度, 可以在不同电压下, 获得更宽的电池功率密度调节范围。
3 结论本文构建了逆流形式的单直通道PEMFC三维模型, 分析了工作温度和阴极湿度对电池性能及内部水分布和水含量的协同影响, 结论如下:
1) 高中电位阶段, 相比温度的变化, 调节湿度对于PEM特性的影响更显著, 低电位阶段PEM特性对温度更敏感;
2) 从电池寿命角度考虑, 中低电位阶段, 在较高的工作温度下, 适当增加电池的阴极湿度, 可防止PEM靠近流道入口处局部干燥, 延长电池寿命。
3) 液态水多分布于肋板下方多孔区域, 这部分区域易发生“水淹”, 致使电池局部电流密度偏低;
4) 从电池输出性能的角度考虑, 相比单独改变燃料电池工作温度或阴极湿度, 协同改变两参数可获得更宽的电池功率密度的调节域。随着电压降低, 适当提高工作温度并降低阴极湿度, 可有效提升电池功率密度, 且可在一定程度上改善多孔区域内水分布特性。
[1] |
李建林, 李光辉, 郭丽军, 等. "十四五"规划下氢能应用技术现状综述及前景展望[J]. 电气应用, 2021, 40(6): 10-16.
LI Jianlin, LI Guanghui, GUO Lijun, et al. Overview and prospect of hydrogen energy application technology under the 14th five year plan[J]. Electrotechnical Application, 2021, 40(6): 10-16. (in Chinese) |
[2] |
马建, 刘晓东, 陈轶嵩, 等. 中国新能源汽车产业与技术发展现状及对策[J]. 中国公路学报, 2018, 31(8): 1-19.
MA Jian, LIU Xiaodong, CHEN Yisong, et al. Current status and countermeasures for China's new energy automobile industry and technology development[J]. China Journal of Highway and Transport, 2018, 31(8): 1-19. (in Chinese) |
[3] | DJILALI N. Computational modelling of polymer electrolyte membrane (PEM) fuel cells: challenges and opportunities[J]. Energy, 2006, 32(4): 269-280. |
[4] | ALI A S, EBRAHIM A, ELNAZ Z, et al. Three-dimensional simulation of different flow fields of proton exchange membrane fuel cell using a multi-phase coupled model with cooling channel[J]. Energy, 2021, 234: 121247. DOI:10.1016/j.energy.2021.121247 |
[5] | TURKMEN A C, CELIK C. The effect of different gas diffusion layer porosity on proton exchange membrane fuel cells[J]. Fuel, 2018, 222: 465-474. DOI:10.1016/j.fuel.2018.02.058 |
[6] | FARCAS A C. A Proposed control technique for water management in proton exchange membrane fuel cells[J]. Applied Mechanics and Materials, 2013, 2776: 145-152. |
[7] |
樊林浩, 张国宾, 焦魁. 高电流、低加湿工况下PEMFC的性能研究[J]. 工程热物理学报, 2019, 40(4): 870-876.
FAN Linhao, ZHANG Guobin, JIAO Kui. Investigations on PEMFC performance under high current density and low humidification[J]. Journal of Engineering Thermophysics, 2019, 40(4): 870-876. (in Chinese) |
[8] | XING L, LIU X T, ALAJE T, et al. A two-phase flow and non-isothermal agglomerate model for a proton exchange membrane(PEM) fuel cell[J]. Energy, 2014, 73: 618-634. DOI:10.1016/j.energy.2014.06.065 |
[9] |
陆佳斌, 申欣明, 陈明, 等. 阴极湿度与电流密度对PEMFC性能的协同影响[J]. 电源技术, 2021, 45(8): 1018-1022.
LU Jiabin, SHEN Xinming, CHEN Ming, et al. Synergistic effect of cathode humidity and current density on performance of PEMFC[J]. Chinese Journal of Power Sources, 2021, 45(8): 1018-1022. (in Chinese) |
[10] |
蒋杨, 焦魁. 质子交换膜燃料电池建模及水热传输特性分析[J]. 热科学与技术, 2019, 18(3): 200-205.
JIANG Yang, JIAO Kui. Modeling and analysis of water and heat transfer characteristics of proton exchange membrane fuel cell[J]. Journal of Thermal Science and Technology, 2019, 18(3): 200-205. (in Chinese) |
[11] | LU J B, WEI G H, ZHU F J, et al. Pressure effect on the PEMFC performance[J]. Fuel Cells, 2019, 19(3): 211-220. DOI:10.1002/fuce.201800135 |
[12] | SHEN J, ZENG L P, TU Z K, et al. Numerical investigation of temperature distribution of proton exchange membrane fuel cells at high current density[J]. Journal of Porous Media, 2019, 22(7): 813-829. DOI:10.1615/JPorMedia.2019028954 |
[13] | LIU Y F, GAO J H, PEI P C, et al. Effects of dynamic changes in inlet temperature on proton exchange membrane fuel cell[J]. Journal of Renewable and Sustainable Energy, 2019, 11(4): 044302. DOI:10.1063/1.5050300 |
[14] | YANG T F, SHEU B H, GHALAMBAZ M, et al. Effects of operating parameters and load mode on dynamic cell performance of proton exchange membrane fuel cell[J]. International Journal of Energy Research, 2020, 45(2): 2474-2487. |
[15] | WEI G H, LU J B, ZHANG Q L, et al. Analyze the effects of flow mode and humidity on PEMFC performance by equivalent membrane conductivity[J]. International Journal of Energy Research, 2019, 43(9): 4592-4605. DOI:10.1002/er.4592 |
[16] | RAJ A, SHAMIM T. Investigation of the effect of multidimensionality in PEM fuel cells[J]. Energy Conversion and Management, 2014, 86: 443-452. DOI:10.1016/j.enconman.2014.04.088 |
[17] |
王世学, 齐贺, 李桦. 一种新型进气方式对燃料电池水管理效果的影响[J]. 热科学与技术, 2013, 12(4): 354-359.
WANG Shixue, QI He, LI Hua. Influence on water management using new inlet gas humidification design for PEFC[J]. Journal of Thermal Science and Technology, 2013, 12(4): 354-359. (in Chinese) |
[18] |
高强, 张拴羊, 徐洪涛, 等. 肋片结构对质子交换膜燃料电池性能影响模拟研究[J]. 热能动力工程, 2020, 35(7): 215-222.
GAO Qiang, ZHANG Shuanyang, XU Hongtao, et al. Simulation study on the effect of fin structure on the performance of proton exchange membrane fuel cell[J]. Journal of Engineering for Thermal Energy and Power, 2020, 35(7): 215-222. (in Chinese) |
[19] | ZHANG G B, JIAO K. Three-dimensional multi-phase simulation of PEMFC at high current density utilizing eulerian-eulerian model and two-fluid model[J]. Energy Conversion and Management, 2018, 176: 409-421. |
[20] | HU G L, LI G N, ZHENG Y Q, et al. Optimization and parametric analysis of PEMFC based on an agglomerate model for catalyst layer[J]. Journal of the Energy Institute, 2014, 87(2): 163-174. DOI:10.1016/j.joei.2014.03.004 |
[21] | LIU H J. Research on liquid water distribution in PEMFC cathode porous media[J]. International Journal of Electrochemical Science, 2020, 15(7): 6717-6736. DOI:10.20964/2020.07.27 |
[22] | XING L, CAI Q, LIU X T, et al. Anode partial flooding modelling of proton exchange membrane fuel cells: optimization of electrode properties and channel geometries[J]. Chemical Engineering Science, 2016, 146: 88-103. |
2. School of Energy and Electrical Engineering, Chang'an University, Xi'an 710018, China