温度和阴极湿度对质子交换膜燃料电池的影响

2024-01-03马菁马强王俊杰郭镇松孙亚松

西北工业大学学报 2023年6期

马菁,马强,王俊杰,郭镇松,孙亚松

(1.长安大学汽车学院,陕西西安 710018; 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几何参数

图1 PEMFC几何模型

1.2 PEMFC数值模型

为防止出现“阳极拖干”,阳极进气湿度假定为100%[10]。在保证各项传输与反应合理的前提下,对数值模型做出如下假设:

1) 工作压力为常压,混合气体遵守理想气体状态方程[18];

2) 在稳态条件下运行;

3) 由于电池内部流速较低,为层流流动[19];

4) 气体扩散层和催化剂层均为各向同性、均质、不可压缩的多孔材料[20];

5) PEM不允许气体通过[21];

6) 不考虑流道内液态水传输。

基于上述假设,建立PEMFC的三维数学模型。

连续方程[20-21]

·(ερu)=Sm

(1)

式中:ρ为混合气体密度,kg/m3;ε为多孔介质孔隙率;u为速度场,m/s;Sm为质量源项,kg/(m3·s)。

动量守恒方程[8,11]

·(ερuu)=-εp+·(εμu)+Su

(2)

式中:μ为混合物的动力黏度,kg/(m·s);Su为动量源项,kg·m/s。

能量守恒方程[21]

·(ερcpuT)=·(keffT)+SQ

(3)

式中:cp为气体比热容,J/(kg·K);T为温度,K;keff为有效导热系数,W/(m·K);SQ为能量源项,W/m3。

组分守恒方程[9,14]

(4)

电荷守恒方程[8]

(5)

(6)

电化学动力学守恒方程[19]

液态水传输方程[8]

(9)

溶解水传输方程[8]

(10)

(11)

σm=(0.513 9λ-0.326)·

(12)

有关数值模型的边界条件为:

阴阳极两侧流道入口为气体入口,其入口质量流量定义为[10]

(13)

(14)

(15)

(16)

上述数值模型的具体参数见表2。

表2 数值模型的参数

1.3 正交设计方案

由于PEMFC常见工作温度范围为60～80℃,且阴极进气加湿不宜超过100%。本文选取3种工作温度(60,70,80℃)和3种阴极湿度(40%,70%,100%)进行正交研究,具体组合方案见表3。

表3 正交设计方案

2 结果与讨论

2.1 网格无关性分析与有效性验证

为了验证仿真模型的正确性,将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。

表4 不同网格数下的电流密度

2.2 反应物浓度分布

图3～4分别给出了C3组合下,0.5 V电压时,PEMFC内部氧气与氢气的分布图。

图3 氧气摩尔分数分布

从图3可以看出,在燃料电池内部多孔区域中,氧气沿着流动方向逐渐被消耗,其浓度分布呈阶梯降低趋势;在阴极肋板下多孔区域的边界处氧气浓度分布较低,最小值约为0.18%。从图4中可以看出,阳极侧氢气浓度分布较为均匀,这保证了阳极侧催化层充分的氢氧化反应;此外,在100%阳极加湿的条件下,电池阳极侧水相摩尔分数维持在50%附近,避免了“阳极拖干”现象。

图4 氢气摩尔分数分布

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含水量与电压关系

从图5可以看出,随着电压的降低,温度越高PEM含水量越低,阴极湿度越高PEM含水量越高。在高电位阶段,阴极湿度对PEM含水量影响最为明显。这是由于在高电位阶段,电化学反应速率较慢,电渗拖拽作用并不显著,此时压差扩散是影响PEM含水量的主要因素。因此,在高电位阶段,阴极湿度对PEM含水量影响最显著。在中电位阶段,随着电压的降低,PEM逐渐湿润,阴极湿度对PEM含水量的影响逐渐削弱。在低电位阶段,各组电池的PEM内部含水量逐步趋于稳定。

图6表明,温度越高,PEM电导率越大;阴极湿度越大,PEM电导率越大。在高电位阶段,与温度相比,湿度对PEM电导率的影响更显著。而在低电位阶段,电导率几乎不受湿度影响。这是含水量逐渐饱和导致的结果。

图6 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含水量

2.4 阴极扩散层液态水饱和度

在高电位阶段,电池内部反应速率较慢,多孔区域发生“水淹”的风险较低。并且,电池反应产生的水会优先补充到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组电池中的最大值。因此,单一地改变温度或湿度对提高电池输出性能具有一定的局限性。同时调节温度和阴极湿度,可以在不同电压下,获得更宽的电池功率密度调节范围。