基于单层电堆形式的质子交换膜燃料电池仿真模拟研究及优化

doi:10.3969/j.issn.2097-0706.2022.08.010

摘要/Abstract

摘要：

质子交换膜燃料电池（PEMFC）由于其较高的发电效率和常温工作稳定性而被广泛关注。传统PEMFC由于电堆内部发电负载不均匀,容易发生“水淹”等问题。大多数研究基于在单电池层面上通过改变流场结构的方法优化PEMFC,但PEMFC工作时往往是组装电堆的形式。为此,设计带有竖直流道作为反应物进出口的单电池电堆结构,并提出一种新型的串联补偿流场,采用3个进气口,并具有首尾互相串联的3个流道,流道之间可以互相补偿反应物。使用ANSYS软件基于单层电堆层面进行仿真模拟研究,保证阳极是相同蛇形流场的情况下,阴极分别采用传统单流道蛇形流场和新型串联补偿流场,分析单层PEMFC电堆中多孔介质层（即扩散层（GDL）和催化层（CL））的电流密度、氧气质量分数、液态水饱和度、速度矢量和流道压力的分布规律。仿真结果表明：传统单流道蛇形流场单层PEMFC电堆中,整体管道压降较高,电流密度和氧气质量分数较高的区域集中在入口段,越靠近出口段分布越少;出口段液态水饱和度较高而速度矢量没有明显升高,液态水容易囤积在出口段;新型串联补偿流场,电流密度分布相对均匀,氧气质量分数有多个较高的区域,分布更加广泛,电化学反应更加充分,整体压降较低,在氧气质量分数较高的区域,速度矢量明显较高,促使氧气在多个流道间快速补充和扩散,整体液态水饱和度相对较低,出口段同样速度矢量较高,排水效果增强。

关键词: 质子交换膜燃料电池, 流场, 电堆, 结构, 仿真, 氢能

Abstract:

Proton exchange membrane fuel cells（PEMFCs） attract extensive attention due to their high power-generation efficiency and working stability at normal temperatures. Traditional PEMFCs are vulnerable to flooding because of the uneven power generation loads among stacks. Most researches on PEMFC optimization are made by improving the flow field of a single cell, but PEMFCs are usually assembled in the form of stacks. Thus, a single-cell stack structure with a vertical channel as the inlet and outlet of reactants is designed, and a new series compensation flow field is proposed which has three air inlets and three flow channels connected in series. Reactants can compensate each other in the flow channel. A simulation study is conducted on the singe-layer PEMFC stack by ANSYS. On the premises that the anode adopts the same serpentine flow field and the cathode adopts the traditional single-channel serpentine flow field or the new series compensation flow field, the distributions of current density, oxygen mass fraction, liquid water saturation, velocity vector and flow channel pressure in the porous medium layers（GDL and CL） of the singe-layer PEMFC stack are analyzed. The simulation results show that the single-layer PEMFC stack in traditional single-channel serpentine flow field is of a higher overall channel pressure, and the closer it is to the inlet, the higher the current density and oxygen mass fraction are, and vice versa. The liquid water saturation around the outlet is relatively high and the velocity vector shows no significant increase as liquid water tends to hoard around the outlet. However, the current density distribution in the new series compensation flow field is relatively even, and there are multiple and extensive-distributed areas with high oxygen mass fractions. The electrochemical reaction is more sufficient and the overall pressure drop is lower. The velocity vector is significantly higher in the areas with higher oxygen mass fraction, which promotes the rapid compensation and diffusion of oxygen between flow channels. The overall liquid water saturation is relatively low. The high velocity vector around the outlet facilitates the drainage effect.

Key words: proton exchange membrane fuel cell, flow field, stack, structure, simulation, hydrogen energy

中图分类号:

TK01

胡翀, 赵袁, 陈代芬. 基于单层电堆形式的质子交换膜燃料电池仿真模拟研究及优化[J]. 综合智慧能源, 2022, 44(8): 91-96.

HU Chong, ZHAO Yuan, RAZA Ali, CHEN Daifen. Simulation and optimization for the PEMFC based on single-cell stack structure[J]. Integrated Intelligent Energy, 2022, 44(8): 91-96.

图/表 13

图1

图2

表1

表2

表3

控制方程

方程	公式及描述
连续性方程	$∇ ⋅ ρ u → = S m$ $∇ ⋅ 1 ε 1 - S ρ u → = S m$ 式中： $ρ$ 和 $u →$ 为气体的密度与速度; $S$ 为液态水的饱和度; $ε$ 为多孔介质层的孔隙率; $S m$ 为质量源项
动量方程	$∇ ⋅ ρ u → u → = - ∇ p + ∇ ⋅ μ ∇ u → + S u ∇ ⋅ 1 ε 2 (1 - s) 2 ρ u → u → = - ∇ p + ∇ ⋅ 1 ε 1 - s μ ∇ u → + S u$ 式中： $μ$ 为气体动力黏度; $p$ 为气体压力; $S$ 为液态水的饱和度; $ε$ 为多孔介质层的孔隙率; $S u$ 为动量源项
组分方程	$∇ ⋅ 1 ε 1 - s c k u → = ∇ ⋅ D k e f ∇ c k + S k$ 式中： $c k$ 为对应组分 $k$ 的浓度; $D k e f$ 为对应组分的有效扩散系数; $S k$ 为对应组分的组分源项
液态水方程	$p c = σ c o s θ c ε K J 1 - s θ c < 90 ∘ σ c o s θ c ε K J s θ c > 90 ∘$ $J x = 1.417 x - 2.12 x 2 + 1.263 x 3$ 式中： $S 为$ 液态水的在多孔介质中的饱和度,其定义为液态水体积占多孔介质中传质通道体积的百分比; $θ c$ 为液滴与壁面的接触角; $σ 为$ 表面张力
溶解态方程	$∇ ⋅ i → m n d F = ∇ ⋅ D w i ∇ λ + S λ + S g d + S l d$ 式中： $i → m$ 为离子电流密度; $λ$ 为计算区域的溶解态水含量; $n d$ 为电渗透系数; $D w i$ 为溶解态水的扩散系数; $S λ$ 为电化学反应产生水所带来的溶解态水的源项; $S g d$ 为溶解态水域气态水转化所产生的源项; $S l d$ 为溶解态水与液态水相互转化所产生的源相
能量守恒方程	$∇ ⋅ ε ρ c p u - T = ∇ ⋅ k e f f ∇ T + S Q$ 式中： $k e f f$ 为有效热导率; $T$ 为温度; $S Q$ 为热量源项
电化学方程	$∇ ⋅ σ s o l ∇ φ s o l + R s o l = 0 ∇ ⋅ σ m e m ∇ φ m e m + R m e m = 0$ 式中： $φ s o l$ 为固相电势即电子产生的电势; $φ m e m$ 为膜相里势即离子电势; $σ$ 为对应导体中对电子或离子的电导率; $R$ 为对应相中电子或离子的源项,即实际产生的电流密度大小

表3

图3

图4

图5

图6

图7

图8

图9

图10

参考文献 13

[1]	郭天民. 准对称固体氧化物燃料电池电极材料的研究与表征[D]. 徐州: 中国矿业大学, 2019.
[2]	胡小夫, 汪洋, 田立, 等. 中高温SOFC/MGT联合发电技术研究进展[J]. 华电技术, 2019, 41(8): 1-5.
	HU Xiaofu, WANG Yang, TIAN Li, et al. Progress in intermediate and high temperature SOFC/MGT combined power generation technology[J]. Huadian Technology, 2019, 41(8): 1-5.
[3]	郭雅婷, 邓甜音, 刘艳莹, 等. 碱性电解水制氢隔膜和阳极材料性能研究[J]. 综合智慧能源, 2022, 44(5): 64-68. doi: 10.3969/j.issn.2097-0706.2022.05.007
	GUO Yating, DENG Tianyin, LIU Yanying, et al. Research on the performance of membranes and anode materials in alkaline water electrolysis[J]. Integrated Intelligent Energy, 2022, 44(5): 64-68. doi: 10.3969/j.issn.2097-0706.2022.05.007
[4]	刘旺玉, 何芋钢, 罗远强, 等. 仿猪笼草结构的质子交换膜燃料电池流道设计[J]. 电源技术, 2022, 46(3): 280-283.
	LIU Wangyu, HE Yugang, LUO Yuanqiang, et al. Design of PEMFC flow channel imitating Nepenthes alata structure[J]. Chinese Journal of Power Sources, 2022, 46(3): 280-283
[5]	于建平, 魏慧利, 许思传. PEMFC蛇形流道几何结构的多目标优化[J]. 佳木斯大学学报(自然科学版), 2022, 40(1): 63-68.
	YU Jianping, WEI Huili, XU Sichuan. Multi-objective optimization of PEMFC serpentine flow channel geometry structure[J]. Journal of Jiamusi University (Natural Science Edition), 2022, 40(1): 63-68.
[6]	陈佳浩, 苏丹丹, 梁玉娇, 等. 基于多物理场的PEMFC流场结构优化[J]. 电池工业, 2021, 25(6): 285-290.
	CHEN Jiahao, SU Dandan, LIANG Yujiao, et al. Optimization of flow field structure of PEMFC based on multiple physical fields[J]. Chinese Battery Industry, 2021, 25(6): 285-290.
[7]	梁凤丽, 闻冉冉, 毛军逵, 等. PEMFC流道结构研究现状及发展趋势[J]. 南京航空航天大学学报, 2021, 53(4): 477-503.
	LIANG Fengli, WEN Ranran, MAO Junkui, et al. Status and development trend of PEMFC flow channel structure[J]. Journal of Nanjing University of Aeronautics and Astronautics, 2021, 53(4): 477-503.
[8]	KURNIA J C, SASMITO A P, SHAMIM T. Performance evaluation of a PEM fuel cell stack with variable inlet flows under simulated driving cycle conditions[J]. Applied Energy, 2017, 206: 751-764. doi: 10.1016/j.apenergy.2017.08.224
[9]	SURESH P V, JAYANTI S, DESHPANDE A P, et al. An improved serpentine flow field with enhanced cross-flow for fuel cell applications[J]. International Journal of Hydrogen Energy, 2011, 36(10) : 6067-6072. doi: 10.1016/j.ijhydene.2011.01.147
[10]	董琪琪. 质子交换膜燃料电池流道内促排水方法研究[D]. 西安: 西北工业大学, 2019.
[11]	王季康, 李华, 彭宇飞, 等. PEMFC建模及性能分析控制[J]. 电子测量技术, 2022, 45(8):27-34.
	WANG Jikang, LI Hua, PENG Yufei, et al. PEMFC modeling and performance analysis control[J]. Electronic Measurement Technology, 2022, 45(8):27-34.
[12]	张立栋, 陈怡冰, 龚明, 等. 质子交换膜电解水制氢影响因素的过程模拟[J]. 综合智慧能源, 2022, 44(5): 88-94. doi: 10.3969/j.issn.2097-0706.2022.05.010
	ZHANG Lidong, CHEN Yibing, GONG Ming, et al. Process simulation of factors affecting proton exchange membrane water electrolysis for hydrogen production[J]. Integrated Intelligent Energy, 2022, 44(5): 88-94. doi: 10.3969/j.issn.2097-0706.2022.05.010
[13]	邹雨廷. 质子交换膜燃料电池动热质电耦合仿真研究与优化设计[D]. 镇江: 江苏科技大学, 2021.

部件	值	单位
流道高度	0.800	mm
流道宽度	0.800	mm
肋宽	0.800	mm
BP厚度	1.600	mm
GDL厚度	0.190	mm
CL厚度	0.015	mm
PEM厚度	0.050	mm
活性面积	510.400	mm²

项目		值	单位
操作温度		60	℃
操作压力		101 325	Pa
BP	材料密度	1 990	kg/m³
	比热容	710	J/(kg·K)
	电导率	92 600	S/m
	热导率	120	W/(m·K)
GDL	材料密度	321.5	kg/m³
	电导率	280	S/m
	液水接触角	110	(°)
	孔隙率	0.6
	渗透率	1.76×10^-11	m^-2
CL	孔隙率	0.4
	渗透率	1.76×10^-11	m^-2
	比表面积	1.127×10⁷	m^-1
PEM	热导率	0.16	W/(m·K)
	干膜密度	1 980	kg/m³
	干膜当量质量	1 100	kg/kmol
电化学反应相关参数	孔隙堵塞饱和度指数	2
	阳极浓度指数	1
	阴极浓度指数	1
	开路电压	1.05	V
	阳极参考浓度	0.881 4	mol/m³
	阴极参考浓度	0.881 4	mol/m³
	阳极电荷转移系数	1
	阳极参考电流密度	7.17	A/m²
	阴极电荷转移系数	1
	阴极参考电流密度	7.17×10^-5	A/m²

[1]	高圆, 李智, 李甲鸿, 高九涛, 李成新, 李长久. 金属支撑固体氧化物燃料电池技术进展[J]. 综合智慧能源, 2022, 44(8): 1-24.
[2]	曹加锋, 李欣然, 邵善德, 冀月霞, 邵宗平. 质子陶瓷燃料电池稳定性研究综述[J]. 综合智慧能源, 2022, 44(8): 58-67.
[3]	许阳森, 张磊, 毕磊. 中温质子导体固体氧化物燃料电池的发展与挑战[J]. 综合智慧能源, 2022, 44(8): 68-74.
[4]	颜畅, 黄晟, 屈尹鹏. 面向碳中和的海上风电制氢技术研究综述[J]. 综合智慧能源, 2022, 44(5): 30-40.
[5]	李梓丘, 乔颖, 鲁宗相. 计及效率与寿命的海上风电-多堆氢能系统运行优化[J]. 综合智慧能源, 2022, 44(5): 69-77.
[6]	张立栋, 陈怡冰, 龚明, 赵桦粮, 王欣, 黄宏艳. 质子交换膜电解水制氢影响因素的过程模拟[J]. 综合智慧能源, 2022, 44(5): 88-94.
[7]	陈逸珲, 林令淇, 田鑫, 张栋梁, 吴军, 刘子晨. 三级式风电AVC协调控制策略[J]. 综合智慧能源, 2022, 44(4): 20-27.
[8]	陈泽泓, 廖锷, 吴磊, 曹淇, 陈国强, 杜广瀚, 刘桂秀, 李根. 太阳能光热电站熔融盐管道支吊架应力分析[J]. 综合智慧能源, 2022, 44(4): 85-91.
[9]	李旭炯, 孙林花, 杨郭明. 基于改进粒子群算法的光伏系统附加向心属性最大功率跟踪研究[J]. 综合智慧能源, 2022, 44(3): 70-76.
[10]	姚哲豪, 郑莆燕, 袁言周. 基于两级策略的双源分布式供能系统运行优化研究[J]. 综合智慧能源, 2022, 44(1): 56-62.
[11]	蒋文坤, 韩颖慧, 薛智文, 朱勇奇, 徐艳梅. 多能互补能源系统中储能原理及其应用[J]. 综合智慧能源, 2022, 44(1): 63-71.
[12]	汪宗恒, 熊鸿韬, 尚磊. 光伏电站半波与全波有功注入阻尼技术研究[J]. 华电技术, 2021, 43(9): 31-36.
[13]	陈洋, 程乐峰, 邹涛. 一种基于蚁群算法的串联电池组路径规划策略[J]. 华电技术, 2021, 43(8): 27-32.
[14]	位帅洁, 李帅辉, 赵志鹏, 李丹. 花状二硫化锡的储钠性能研究[J]. 华电技术, 2021, 43(7): 37-41.
[15]	喻小宝, 郑丹丹, 杨康, 孔杰, 章天浩. “双碳”目标下能源电力行业的机遇与挑战[J]. 华电技术, 2021, 43(6): 21-32.