考虑扩散层表面润湿性的流道内两相流对燃料电池性能影响研究

doi:10.11949/0438-1157.20250452

摘要/Abstract

摘要：

质子交换膜燃料电池流道内的两相流是影响其性能的关键因素，但考虑气液界面演化等详细两相流特性的多物理场耦合研究仍然不足。通过构建耦合流体体积模型的燃料电池模型，定量解析了气体扩散层（gas diffusion layer，GDL）表面润湿性对燃料电池传质及性能的影响。结果表明，忽略流道内详细两相流会高估燃料电池性能及物理量分布均匀性，而随着GDL表面疏水性的增强，GDL表面的液态水覆盖面积下降，在前期有利于增强氧气输运及燃料电池性能，但也会导致排水变慢，在后期导致燃料电池性能恢复较慢。在本研究中，GDL表面接触角为160°时燃料电池性能最佳，而接触角为110°时性能最低。本研究为理解燃料电池流道内两相流影响机制、设计高性能燃料电池提供了理论及模型支撑。

关键词: 燃料电池, 流道, 气液两相流, 性能, 均匀性, 计算流体力学

Abstract:

The two-phase flow within flow channels of proton exchange membrane fuel cells is a critical factor affecting the performance, but the multi-physics coupling research considering detailed two-phase flow characteristics such as the evolution of the gas-liquid interface is still insufficient. By constructing a fuel cell model coupled with the volume of fluid model, this study quantitatively analyzes the impact of the surface wettability of the gas diffusion layer (GDL) on mass transfer and performance in fuel cells. The results show that neglecting the detailed two-phase flow in channels can overestimate both the fuel cell performance and the uniformity of physical quantity distribution. As the hydrophobicity of the GDL surface increases, the coverage area of liquid water on the GDL surface decreases, which initially enhances oxygen transport and fuel cell performance. However, this also leads to slower water drainage, causing slower performance recovery in the later stage. In this study, the fuel cell performance is best when the contact angle of the GDL surface is 160°, while the performance is lowest when the contact angle is 110°. The study provides theoretical and model support for understanding the mechanisms underlying the influence of two-phase flow within flow channels and for the design of high-performance fuel cells.

Key words: fuel cell, flow channel, gas-liquid flow, performance, uniformity, computational fluid dynamics

中图分类号:

TM 911

张晓卿, 马骁, 帅石金. 考虑扩散层表面润湿性的流道内两相流对燃料电池性能影响研究[J]. 化工学报, 2025, 76(11): 5574-5583.

Xiaoqing ZHANG, Xiao MA, Shijin SHUAI. Study on influence of two-phase flow in channels considering surface wettability of gas diffusion layers on fuel cell performance[J]. CIESC Journal, 2025, 76(11): 5574-5583.

图/表 14

图1 计算域及网格划分ACH—阳极流道;CCH—阴极流道;AGDL—阳极气体扩散层;CGDL—阴极气体扩散层;AMPL—阳极微孔层; CMPL—阴极微孔层;ACL—阳极催化层;CCL—阴极催化层;MEM—膜

Fig.1 Computational domain and mesh

表1 主要几何参数及工作条件

Table 1 Main physical parameters and operating conditions

参数	数值
流道长度/mm	50
流道高度/mm	0.5
厚度(GDL、CL、MEM)/mm	0.2、0.01、0.0254
接触角(GDL、CL)/(°)	100~160、100
渗透率(GDL、CL)/m²	1.0×10^-12、2.0×10^-20
空气流速/(m/s)	1
工作压力/atm	1.5
工作温度/K	353.15

表2 流体体积模型

Table 2 Volume of fluid model

名称	方程
质量守恒	$∂ ρ ∂ t + ∇ ⋅ ρ u = 0$
动量守恒	$∂ ∂ t ρ u + ∇ ⋅ ρ u ⋅ u = - ∇ P + μ ∇ ⋅ ∇ u + ∇ u T + ρ g + F s$
相守恒	$∂ f 1 ∂ t + ∇ ⋅ f 1 u = 0$
表面张力源项	$F s = - σ 2 ρ κ ∇ f 1 ρ 1 + ρ 2$
表面曲率	$κ = ∇ n w c o s θ + t w s i n θ$

表2 流体体积模型

Table 2 Volume of fluid model

名称	方程
质量守恒	$∂ ρ ∂ t + ∇ ⋅ ρ u = 0$
动量守恒	$∂ ∂ t ρ u + ∇ ⋅ ρ u ⋅ u = - ∇ P + μ ∇ ⋅ ∇ u + ∇ u T + ρ g + F s$
相守恒	$∂ f 1 ∂ t + ∇ ⋅ f 1 u = 0$
表面张力源项	$F s = - σ 2 ρ κ ∇ f 1 ρ 1 + ρ 2$
表面曲率	$κ = ∇ n w c o s θ + t w s i n θ$

表3 守恒方程

Table 3 Conservation equations

名称	方程
质量守恒	$∂ ∂ t ε 1 - s ρ g + ∇ ⋅ ρ g u g = S m$
动量守恒	$∂ ∂ t ρ g u g ε 1 - s + ∇ ⋅ ρ g u g u g ε 2 1 - s 2 = - ∇ P g + μ g ∇ ⋅ ∇ u g ε 1 - s + ∇ u g T ε 1 - s - 23 μ g ∇ ⋅ ∇ u g ε 1 - s + S u$
气相组分守恒	$∂ ∂ t ε 1 - s ρ g Y i + ∇ ⋅ ρ g u g Y i = ∇ ⋅ ρ g D e f f i ∇ Y i + S i$
液态水守恒	$∂ ∂ t ρ l ε s = ∇ ⋅ ρ l K k l μ l ∇ P l + S l$
膜态水守恒	$ρ M E M E W ∂ ∂ t ω λ + ∇ ⋅ n d F I i o n = ρ M E M E W ∇ ⋅ D d e f f ∇ λ + S m w$
电子守恒	$0 = ∇ ⋅ k e e f f ∇ φ e + S e$
离子守恒	$0 = ∇ ⋅ k i o n e f f ∇ φ i o n + S i o n$

表3 守恒方程

Table 3 Conservation equations

名称	方程
质量守恒	$∂ ∂ t ε 1 - s ρ g + ∇ ⋅ ρ g u g = S m$
动量守恒	$∂ ∂ t ρ g u g ε 1 - s + ∇ ⋅ ρ g u g u g ε 2 1 - s 2 = - ∇ P g + μ g ∇ ⋅ ∇ u g ε 1 - s + ∇ u g T ε 1 - s - 23 μ g ∇ ⋅ ∇ u g ε 1 - s + S u$
气相组分守恒	$∂ ∂ t ε 1 - s ρ g Y i + ∇ ⋅ ρ g u g Y i = ∇ ⋅ ρ g D e f f i ∇ Y i + S i$
液态水守恒	$∂ ∂ t ρ l ε s = ∇ ⋅ ρ l K k l μ l ∇ P l + S l$
膜态水守恒	$ρ M E M E W ∂ ∂ t ω λ + ∇ ⋅ n d F I i o n = ρ M E M E W ∇ ⋅ D d e f f ∇ λ + S m w$
电子守恒	$0 = ∇ ⋅ k e e f f ∇ φ e + S e$
离子守恒	$0 = ∇ ⋅ k i o n e f f ∇ φ i o n + S i o n$

表4 控制方程源项

Table 4 Source terms

源项	单位
$S m = - S v - l G D L s - S v - l + S d - v M H 2 O - J a 2 F M H 2 A C L - S v - l + S d - v M H 2 O + J c 2 F M H 2 O - J c 4 F M O 2 C C L$	kg/(m³·s)
$S u = - μ g K k g u g G D L s, C L s$	kg/(m²·s)
$S H 2 = - J a 2 F M H 2 A C L$	kg/(m³·s)
$S O 2 = - J c 4 F M O 2 C C L$	kg/(m³·s)
$S H 2 O = - S v - l G D L s - S v - l + S d - v M H 2 O A C L - S v - l + S d - v M H 2 O + J c 2 F M H 2 O C C L$	kg/(m³·s)
$S l = S v - l G D L s, C L s$	kg/(m³·s)
$S m w = S v - d + ρ l K M E M μ l M H 2 O P l C C L - P l A C L δ M E M δ C L A C L S v - d - ρ l K M E M μ l M H 2 O P l C C L - P l A C L δ M E M δ C L C C L$	kg/(m³·s)
$S v - l = γ c o n d ε 1 - s C - H 2 O C s a t M H 2 O C H 2 O > C s a t γ e v a p ε s C - H 2 O C s a t M H 2 O C H 2 O < C s a t$	kg/(m³·s)
$S d - v = 1.3 ρ M E M E W λ - λ e q$	mol/(m³·s)
$S e l e = - J a A C L J c C C L$	A/m³
$S i o n = J a A C L - J c C C L$	A/m³

表4 控制方程源项

Table 4 Source terms

源项	单位
$S m = - S v - l G D L s - S v - l + S d - v M H 2 O - J a 2 F M H 2 A C L - S v - l + S d - v M H 2 O + J c 2 F M H 2 O - J c 4 F M O 2 C C L$	kg/(m³·s)
$S u = - μ g K k g u g G D L s, C L s$	kg/(m²·s)
$S H 2 = - J a 2 F M H 2 A C L$	kg/(m³·s)
$S O 2 = - J c 4 F M O 2 C C L$	kg/(m³·s)
$S H 2 O = - S v - l G D L s - S v - l + S d - v M H 2 O A C L - S v - l + S d - v M H 2 O + J c 2 F M H 2 O C C L$	kg/(m³·s)
$S l = S v - l G D L s, C L s$	kg/(m³·s)
$S m w = S v - d + ρ l K M E M μ l M H 2 O P l C C L - P l A C L δ M E M δ C L A C L S v - d - ρ l K M E M μ l M H 2 O P l C C L - P l A C L δ M E M δ C L C C L$	kg/(m³·s)
$S v - l = γ c o n d ε 1 - s C - H 2 O C s a t M H 2 O C H 2 O > C s a t γ e v a p ε s C - H 2 O C s a t M H 2 O C H 2 O < C s a t$	kg/(m³·s)
$S d - v = 1.3 ρ M E M E W λ - λ e q$	mol/(m³·s)
$S e l e = - J a A C L J c C C L$	A/m³
$S i o n = J a A C L - J c C C L$	A/m³

图2 模型验证

Fig.2 Model validation

图3 流道内两相流输运特性对比

Fig.3 Comparison of two-phase flow in the flow channel

图4 流道内排水过程对比

Fig.4 Comparison of drainage processes within flow channels

图5 氧气输运特性对比

Fig.5 Comparison of oxygen transport characteristics

图6 考虑流道内详细两相流时CL内氧气输运特性的对比

Fig.6 Comparison of oxygen transport characteristics in the CL

图7 GDL内氧气浓度演化对比

Fig.7 Comparison of oxygen concentration evolution in the GDL

图8 电流密度及反应速率对比

Fig.8 Comparison of current density and reaction rate

图9 考虑流道内详细两相流时电流密度及反应均匀性的对比

Fig.9 Comparison of current density and reaction uniformity

图10 CL内反应速率演化对比

Fig.10 Comparison of reaction rate evolution in the CL

参考文献 31

[1]	Zhang X Q, Ma X, Zhang Z H, et al. Review and analysis of thermal management for proton exchange membrane fuel cell hybrid power system[J]. Renewable Energy, 2025, 244: 122716.
[2]	Jiao K, Xuan J, Du Q, et al. Designing the next generation of proton-exchange membrane fuel cells[J]. Nature, 2021, 595(7867): 361-369.
[3]	Liu H, Chen J, Hissel D, et al. Prognostics methods and degradation indexes of proton exchange membrane fuel cells: a review[J]. Renewable and Sustainable Energy Reviews, 2020, 123: 109721.
[4]	Liu L N, Guo L Y, Zhang R Y, et al. Numerically investigating two-phase reactive transport in multiple gas channels of proton exchange membrane fuel cells[J]. Applied Energy, 2021, 302: 117625.
[5]	Wang Y, Ruiz Diaz D F, Chen K S, et al. Materials, technological status, and fundamentals of PEM fuel cells—a review[J]. Materials Today, 2020, 32: 178-203.
[6]	Xiong K N, Wu W, Wang S F, et al. Modeling, design, materials and fabrication of bipolar plates for proton exchange membrane fuel cell: a review[J]. Applied Energy, 2021, 301: 117443.
[7]	Xu S N, Liao P Y, Yang D J, et al. Liquid water transport in gas flow channels of PEMFCs: a review on numerical simulations and visualization experiments[J]. International Journal of Hydrogen Energy, 2023, 48(27): 10118-10143.
[8]	Zhou Y L, Chang H, Qi T Y. Gas-liquid two-phase flow in serpentine microchannel with different wall wettability[J]. Chinese Journal of Chemical Engineering, 2017, 25(7): 874-881.
[9]	Zhang X Q, Ma X, Shuai S J. Impact of detailed liquid water transport in channel on mass transfer and performance of proton exchange membrane fuel cell[J]. International Communications in Heat and Mass Transfer, 2025, 162: 108515.
[10]	Bazylak A. Liquid water visualization in PEM fuel cells: a review[J]. International Journal of Hydrogen Energy, 2009, 34(9): 3845-3857.
[11]	Ferreira R B, Falcão D S, Oliveira V B, et al. Numerical simulations of two-phase flow in proton exchange membrane fuel cells using the volume of fluid method—a review[J]. Journal of Power Sources, 2015, 277: 329-342.
[12]	Li M J, Zhang E R, Zhang M F, et al. Ex-situ experimental study on extraction of droplet dynamic parameters based on droplet shape in PEMFC[J]. International Journal of Green Energy, 2025, 22(5): 858-865.
[13]	Xu Y F, Peng L F, Yi P Y, et al. Numerical investigation of liquid water dynamics in wave-like gas channels of PEMFCs[J]. International Journal of Energy Research, 2019, 43(3): 1191-1202.
[14]	Qiu D, Xu Z, Shao H, et al. Analytical modelling of water droplet behavior at the gas channel corner for proton exchange membrane fuel cells[J]. Journal of Electrochemical Energy Conversion and Storage, 2025, 22(1): 011001.
[15]	Chen J X, Bao Z M, Xu Y F, et al. Investigation of liquid retention behavior in the flow field plate of large-size proton exchange membrane fuel cells: effects of sub-distribution zone[J]. Applied Energy, 2024, 358: 122651.
[16]	Zhang X Q, Yang J P, Ma X, et al. Numerical investigation of water dynamics in a novel wettability gradient anode flow channel for proton exchange membrane fuel cells[J]. International Journal of Energy Research, 2020, 44(13): 10282-10294.
[17]	Ding Y L, Bi X T, Wilkinson D P. 3D simulations of the impact of two-phase flow on PEM fuel cell performance[J]. Chemical Engineering Science, 2013, 100: 445-455.
[18]	Zhang G B, Wu L Z, Qin Z K, et al. A comprehensive three-dimensional model coupling channel multi-phase flow and electrochemical reactions in proton exchange membrane fuel cell[J]. Advances in Applied Energy, 2021, 2: 100033.
[19]	Zhang X Q, Ma X, Qin Y Z, et al. Modeling research of the impact of liquid water in channel on gas and water transport in proton exchange membrane fuel cell[J]. International Journal of Heat and Mass Transfer, 2024, 235: 126127.
[20]	Penga Ž, Bergbreiter C, Barbir F, et al. Numerical and experimental analysis of liquid water distribution in PEM fuel cells[J]. Energy Conversion and Management, 2019, 189: 167-183.
[21]	Ding Y J, Xu L F, Zheng W B, et al. Characterizing the two-phase flow effect in gas channel of proton exchange membrane fuel cell with dimensionless number[J]. International Journal of Hydrogen Energy, 2023, 48(13): 5250-5265.
[22]	Chen H, Guo H, Ye F, et al. Improving two-phase mass transportation under non-Darcy flow effect in orientated-type flow channels of proton exchange membrane fuel cells[J]. International Journal of Hydrogen Energy, 2021, 46(41): 21600-21618.
[23]	Ding Y, Bi X T, Wilkinson D P. Numerical investigation of the impact of two-phase flow maldistribution on PEM fuel cell performance[J]. International Journal of Hydrogen Energy, 2014, 39(1): 469-480.
[24]	Ferreira R B, Falcão D S, Oliveira V B, et al. 1D+3D two-phase flow numerical model of a proton exchange membrane fuel cell[J]. Applied Energy, 2017, 203: 474-495.
[25]	Le A D, Zhou B. A generalized numerical model for liquid water in a proton exchange membrane fuel cell with interdigitated design[J]. Journal of Power Sources, 2009, 193(2): 665-683.
[26]	Zhang X Q, Ma X, Yang J P, et al. Effect of liquid water in flow channel on proton exchange membrane fuel cell: focusing on flow pattern[J]. Energy Conversion and Management, 2022, 258: 115528.
[27]	Liu H C, Tan J, Cheng L S, et al. Enhanced water removal performance of a slope turn in the serpentine flow channel for proton exchange membrane fuel cells[J]. Energy Conversion and Management, 2018, 176: 227-235.
[28]	Zhang G B, Qu Z G, Tao W Q, et al. Porous flow field for next-generation proton exchange membrane fuel cells: materials, characterization, design, and challenges[J]. Chemical Reviews, 2023, 123(3): 989-1039.
[29]	Ma X, Zhang X Q, Yang J P, et al. Impact of gas diffusion layer spatial variation properties on water management and performance of PEM fuel cells[J]. Energy Conversion and Management, 2021, 227: 113579.
[30]	Jiao K, Li X G. Water transport in polymer electrolyte membrane fuel cells[J]. Progress in Energy and Combustion Science, 2011, 37(3): 221-291.
[31]	Theodorakakos A, Ous T, Gavaises M, et al. Dynamics of water droplets detached from porous surfaces of relevance to PEM fuel cells[J]. Journal of Colloid and Interface Science, 2006, 300(2): 673-687.