质子交换膜燃料电池仿真模型并行计算与优化

doi:10.16183/j.cnki.jsjtu.2023.636

摘要/Abstract

摘要：

为实现质子交换膜燃料电池与电堆的高效和高精度仿真,基于Fluent求解器开发一套三维耦合一维的混合维度仿真程序,深入分析Fluent内置的global summations宏PRF_GRSUM对燃料电池仿真程序并行求解效率的影响,并基于混合维度模型的数据耦合特性,提出数据压缩同步法和数据筛选标记法两种新的并行策略.另外,详细阐述了3种数据交换方法的原理和流程,并通过简化模型测试其对程序求解效率的影响,结果表明采用数据筛选标记法相比PRF_GRSUM宏方法,计算包含百节电池的燃料电池堆的加速比达到112.5倍.

关键词: 质子交换膜燃料电池, 并行计算, 混合维度模型, 数据筛选标记法

Abstract:

To achieve highly efficient and precise simulation of proton exchange membrane fuel cells (PEMFCs) and stacks, this paper develops a three-dimensional coupled one-dimensional hybrid dimension modeling strategy based on the Fluent solver and conducts an in-depth analysis of the impact of the Fluent built-in global summations macro PRF_GRSUM on the parallel solving efficiency of fuel cell simulation. It also proposes two new parallel strategies including the data compression synchronization method and the data filter marker method, based on the data coupling characteristics of the hybrid dimension model. Then, this paper describes the principles and processes of three data synchronization methods in detail, and investigates the impact of these methods on the solving efficiency of parallel computation of PEMFC models. The results show that the data filter marker method can achieve an acceleration ratio up to 112.5 times compared with PRF_GRSUM macro when simulating a PEMFC stack with a hundred cells.

Key words: proton exchange membrane fuel cell (PEMFC), parallel computation, hybrid dimension model, data filter marker method

中图分类号:

刘克, 郝亮, 耿珺, 陈洵, 陆维. 质子交换膜燃料电池仿真模型并行计算与优化[J]. 上海交通大学学报, 2025, 59(11): 1754-1762.

LIU Ke, HAO Liang, GENG Jun, CHEN Xun, LU Wei. Parallel Computation and Optimization of Proton Exchange Membrane Fuel Cell Models[J]. Journal of Shanghai Jiao Tong University, 2025, 59(11): 1754-1762.

图/表 13

图1

图2

图3

图4

表1

图5

图6

图7

图8

表2

图9

图10

图11

参考文献 25

[1]	BARBIR F. PEM 燃料电池:理论与实践[M]. 第2版. 北京: 机械工业出版社, 2016.
	BARBIR F. PEM Fuel cells: Theory and practice[M]. 2nd ed. Beijing: China Machine Press, 2016.
[2]	YUAN S, ZHAO C F, CAI X Y, et al. Bubble evolution and transport in PEM water electrolysis: Mechanism, impact, and management[J]. Progress in Energy & Combustion Science, 2023, 96: 101075.
[3]	YOU J B, ZHENG Z F, CHENG X J, et al. Insight into oxygen transport in solid and high-surface-area carbon supports of proton exchange membrane fuel cells[J]. ACS Applied Materials & Interfaces, 2023, 15(17): 21457-21466.
[4]	SHEN J, TU Z K, CHAN S H. Enhancement of mass transfer in a proton exchange membrane fuel cell with blockage in the flow channel[J]. Applied Thermal Engineering, 2019, 149: 1408-1418. doi: 10.1016/j.applthermaleng.2018.12.138 URL
[5]	郑宇. PEMFC传热传质研究[D]. 济南: 山东大学, 2021.
	ZHENG Yu. Research on heat and mass transfer of PEMFC[D]. Jinan: Shandong University, 2021.
[6]	胡经纬. 质子交换膜燃料电池的电化学和数值模拟研究[D]. 大连: 中国科学院大连化学物理研究所, 2006.
	HU Jingwei. Electrochemical and mathematical model studies on PEMFC[D]. Dalian: Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 2006.
[7]	LIU Z X, MAO Z Q, WANG C. A two dimensional partial flooding model for PEMFC[J]. Journal of Power Sources, 2006, 158(2): 1229-1239. doi: 10.1016/j.jpowsour.2005.10.060 URL
[8]	SHEN J, XU L, CHANG H W, et al. Partial flooding and its effect on the performance of a proton exchange membrane fuel cell[J]. Energy Conversion & Management, 2020, 207: 112537.
[9]	HAO L, CHENG P. Lattice Boltzmann simulations of anisotropic permeabilities in carbon paper gas diffusion layers[J]. Journal of Power Sources, 2009, 186(1): 104-114. doi: 10.1016/j.jpowsour.2008.09.086 URL
[10]	DU Z P, STEINDL C, JAKUBEK S. Efficient two-step parametrization of a control-oriented zero-dimensional polymer electrolyte membrane fuel cell model based on measured stack data[J]. Processes, 2021, 9(4): 713. doi: 10.3390/pr9040713 URL
[11]	ROWE A, LI X G. Mathematical modeling of proton exchange membrane fuel cells[J]. Journal of Power Sources, 2001, 102(1/2): 82-96. doi: 10.1016/S0378-7753(01)00798-4 URL
[12]	BERNARDI D M, VERBRUGGE M W. A mathematical model of the solid-polymer-electrolyte fuel cell[J]. Journal of the Electrochemical Society, 1992, 139(9): 2477-2491. doi: 10.1149/1.2221251
[13]	SPRINGER T E, ZAWODZINSKI T A, GOTTESFELD S. Polymer electrolyte fuel cell model[J]. Journal of the Electrochemical Society, 1991, 138(8): 2334-2342. doi: 10.1149/1.2085971
[14]	SINGH D, LU D M, DJILALI N. A two-dimensional analysis of mass transport in proton exchange membrane fuel cells[J]. International Journal of Engineering Science, 1999, 37(4): 431-452. doi: 10.1016/S0020-7225(98)00079-2 URL
[15]	BAPAT C J, THYNELL S T. Effect of anisotropic electrical resistivity of gas diffusion layers (GDLs) on current density and temperature distribution in a polymer electrolyte membrane (PEM) fuel cell[J]. Journal of Power Sources, 2008, 185(1): 428-432. doi: 10.1016/j.jpowsour.2008.06.016 URL
[16]	UM S, WANG C Y. Three-dimensional analysis of transport and electrochemical reactions in polymer electrolyte fuel cells[J]. Journal of Power Sources, 2004, 125(1): 40-51. doi: 10.1016/j.jpowsour.2003.07.007 URL
[17]	ROSLI M I, LIM B H, MAJLAN E H, et al. Performance analysis of PEMFC with single-channel and multi-channels on the impact of the geometrical model[J]. Energies, 2022, 15(21): 7960. doi: 10.3390/en15217960 URL
[18]	BAO C, BESSLER W G. A computationally efficient steady-state electrode-level and 1D+1D cell-level fuel cell model[J]. Journal of Power Sources, 2012, 210: 67-80. doi: 10.1016/j.jpowsour.2012.03.023 URL
[19]	SCHUMACHER J O, ELLER J, SARTORIS G, et al. 2+1D modelling of a polymer electrolyte fuel cell with glassy-carbon microstructures[J]. Mathematical & Computer Modelling of Dynamical Systems, 2012, 18(4): 355-377.
[20]	CORDINER S, MULONE V, ROMANELLI F. Thermal-fluid-dynamic simulation of a proton exchange membrane fuel cell using a hierarchical 3D-1D approach[J]. Journal of Fuel Cell Science & Technology, 2007, 4(3): 317-327.
[21]	XIE B, ZHANG G B, JIANG Y, et al. “3D+1D” modeling approach toward large-scale PEM fuel cell simulation and partitioned optimization study on flow field[J]. eTransportation, 2020, 6: 100090. doi: 10.1016/j.etran.2020.100090 URL
[22]	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. doi: 10.1016/j.apenergy.2017.06.048 URL
[23]	ANSYS. ANSYS Fluent customization manual[M]. Canonsburg, USA: ANSYS, 2023.
[24]	LIU K, STEFANI F, WEBER N, et al. Numerical and experimental investigation of electro-vortex flow in a cylindrical container[J]. Magnetohydrodynamics, 2020, 56(1): 27-42. doi: 10.22364/mhd URL
[25]	ANSYS. ANSYS Fluent user’s guide[M]. Canonsburg, USA: ANSYS, 2023.

方法	预处理	一维数组及个数	重复初始化	初始化时刻	初始化对象	各节点一致性	优点	缺点
方法1	无	1个:数组_0	需要	每次数据同步前	数组_0	一致	易实现	数据交换量最大
方法2	需编号、统计数据等	2个:数组_0、数组_1	需要	每次数据同步前	数组_1	数组_1一致	效率较方法1高	编程过程繁琐
方法3	需编号、标识数据	1个:数组_0	不需要	每次数据同步前	数组_0	不一致	效率最高	易出错,难度较大

需要交换数据的量	所用核数	消耗时长/s
需要交换数据的量	所用核数	方法1	方法2	方法3
56	2	1261.240	839.320	836.860
135	4	993.270	437.332	436.217
358	8	1017.454	256.242	256.474
710	16	1098.927	145.131	149.206
1271	32	1255.663	115.212	123.052
1823	48	2040.343	107.112	121.867
2491	64	2462.074	132.015	156.739
3938	96	3026.511	132.948	185.341