可再生能源

质子交换膜燃料电池热管理系统建模及故障仿真

PDF下载

雍加望 ¹ , 赵倩倩 ² , 冯能莲 ²

可再生能源 | 2024,42(3): 308-316

收起

可再生能源 | 2024, 42(3): 308-316

质子交换膜燃料电池热管理系统建模及故障仿真

全屏

雍加望¹, 赵倩倩², 冯能莲²

作者信息

¹ 北京工业大学交通工程系北京 100124

² 北京工业大学汽车工程系北京 100124

雍加望(1988-),男,博士,讲师,研究方向为自动驾驶汽车、汽车动力学与稳定性控制。E-mail: yongjw@bjut.edu.cn。

Modeling and malfunction simulation of thermal management system for proton exchange membrane fuel cells

Jiawang Yong¹, Qianqian Zhao², Nenglian Feng²

Affiliations

¹ Department of Traffic Engineering Beijing University of Technology Beijing 100124 China

² Department of Vehicle Engineering Beijing University of Technology Beijing 100124 China

出版时间: 2024-03-20

文章导航

摘要

收起

文章以质子交换膜燃料电池(Proton Exchange Membrane Fuel Cell, PEMFC)系统为研究对象,分析系统中各组成部件的工作机理,在MATLAB/Simulink 软件中搭建了燃料电池电堆电压、阴极、质子交换膜、阳极以及温度的数学机理模型,同时使用 Simulink/Simscape 物理建模平台搭建热管理系统的物理模型,将两种模型集成为完整的PEMFC系统仿真模型。在所搭建的热管理系统模型中注入典型故障:散热器风扇故障和冷却液流量不足故障,对各故障模式下的燃料电池性能进行了分析。仿真结果表明:PEMFC 系统模型结果与试验结果基本一致,验证了模型的合理性;通过对热管理系统故障进行仿真,可以清楚地了解故障发生的机理,为其故障诊断提供了参考依据。

关键词

质子交换膜燃料电池 / 热管理系统 / Simulink/Simscape / 建模 / 故障仿真

Abstract

收起

This paper analyzes the working mechanism of each components of the proton exchange membrane fuel cell (PEMFC) system, and the mathematical mechanism models of stack, anode, cathode, proton exchange membrane, and temperature of the PEMFC are established using MATLAB/Simulink software, furthermore, the physical model of thermal management system is built in the Simulink/Simscape environment, and the mathematical and physical models are integrated into a complete PEMFC system simulation model. Typical malfunctions, including the radiator fan failure and insufficient coolant flow failure, are injected into the PEMFC system simulation model to analyze the influence of malfunctions on the performance of PEMFC. The simulation results are basically consistent with the experimental results, which indicates that the proposed model is reasonable and accurate. Moreover, the generation mechanism of malfunctions is figured out through malfunction simulation of the thermal management system, which provides a reference for malfunction diagnosis.

Key words

proton exchange membrane fuel cell / thermal management system / simulink/ simscape / modeling / malfunction simulation

引用本文

雍加望, 赵倩倩, 冯能莲. 质子交换膜燃料电池热管理系统建模及故障仿真. 可再生能源, 2024 , 42 (3) : 308 -316 .

Jiawang Yong, Qianqian Zhao, Nenglian Feng. Modeling and malfunction simulation of thermal management system for proton exchange membrane fuel cells[J]. Renewable Energy Resources, 2024 , 42 (3) : 308 -316 .

正文

收起

0 引言

收起

能源危机和环境污染问题使得潜力巨大的燃料电池技术得到快速发展。在所有类型的燃料电池中, 质子交换膜燃料电池 (Proton Exchange Membrane Fuel Cell, PEMFC) 具有广阔的应用前景。在催化剂的作用下, PEMFC 中的氢气与空气中的氧气发生电化学反应只生成水, PEMFC 具有能量转换效率高、低温快速启动、环境友好、噪声小以及与车辆中其他子系统兼容集成性能好等优势 ^{[ 1 , 2 ]} 。由于这些潜在的优势,使得 PEMFC 被开发应用在新能源汽车领域, 并得到了许多政府和汽车企业的认可。

针对燃料电池系统模型的搭建, 国内外研究人员进行了广泛而深入的研究。Omran A ^{[ 3 ]} 建立了 PEMFC 数学模型(包括电堆模型、升压转换器模型、外部负载模型和性能参数计算模型),并通过仿真与实验手段验证了外部负载变化工况下的系统整体效率。Le A D ^{[ 4 ]} 建立了考虑流体流动、传热、电化学和电流密度分布等因素的通用 PEMFC 模型, 并研究了液态水对蛇形通道 PEMFC 的影响。房鑫 ^{[ 5 ]} 将机理模型和辨识模型结合,建立了 PEMFC 混合动态模型。郭建忠 ^{[ 6 ]} 采用 MATLAB/ Simulink 中的 S-Function 模块搭建了 PEMFC 的电堆模型以及气路和水路模型。卫超强 ^{[ 7 ]} 和石磊 ^{[ 8 ]} 通过 COMSOL 软件搭建了 PEMFC 单体电池模型, 分别研究了工作温度和流道结构对电堆内水分布和电池输出性能的影响。

基于 Simscape 的物理建模方法具有模型结构简单、计算准确性高的特点。Simscape 建模是根据实际系统的物理结构搭建模型, 不需要按照传统的 Simulink 数学计算模块进行建模,并且通过数学信号-物理信号转换模块, 可将构建的 Simscape 系统模型与普通 Simulink 模型直接集成 ^{[ 9 ]} 。现有的燃料电池系统研究主要是基于数学原理进行建模仿真, 很少有学者基于 Simscape 物理建模方法对其进行研究, 并且对故障工况下燃料电池性能影响的研究较少。因此, 本文在文献 [ 10 ]中模型的基础上, 采用 Simulink 软件搭建了燃料电池温度模型, 并基于 Simscape 物理建模平台建立热管理系统模型, 之后对热管理系统进行了故障仿真,可为 PEMFC 故障诊断提供依据,有助于延长 PEMFC 使用寿命并提高安全性。

1 PEMFC 系统

收起

燃料电池系统由燃料电池堆、氢气供给系统、空气供给系统、增湿系统、热管理系统以及控制系统组成, 图 1 为其组成示意图。在实际应用中, 燃料电池堆是根据负载的功率需求由若干片单片电池串联组成；氢气供给系统采用氢气瓶供给燃料, 经过减压阀输入燃料电池堆中进行电化学反应, 反应剩余的氢气通过氢气循环泵再次送入电堆中；空气供给系统由空气过滤器及压缩机组成, 可向燃料电池系统输送氧气；增湿系统可控制电堆内的含水量(若含水量过多则会造成电极水淹和淹膜,含水量太少则会造成膜过干,两者都会使燃料电池性能下降), 该子系统的主要设备是增湿器; 热管理系统的主要任务是保持电堆的工作温度处于最佳区间内,其主要包括散热器、冷却液泵和冷却液箱; 燃料电池系统在工作时会产生大量的热,需要冷却液循环及时带走多余的热量,否则电池过热会导致膜过干,使燃料电池系统无法正常工作, 当环境温度较低时, 为了防止反应生成物 (水)的冻结,需要加热冷却液进而对电堆加热；控制系统中的控制策略用来控制各种阀门、泵、散热器风扇等,进而调节气、水、热的供给,保证燃料电池的输出性能。

2 热管理系统建模及验证

收起

2.1 PEMFC 电堆模型

PEMFC 电堆模型采用文献[ 10 ]中的模型, 其中包括电堆电压模型、阴极模型、质子交换膜水模型和阳极模型,其模型框图见图 2 。

2.2 温度模型

燃料电池电堆内的能量平衡方程可表示为

(1)${q}_{\text{theo }}= {q}_{\text{elec }}+ {q}_{\text{sens }}+ {q}_{\text{latent }}+ {q}_{\text{cool }}$

式中: ${q}_{\text{theo }}$ 为理论上燃料电池电化学反应产生的能量, kW； ${q}_{\text{elec }}$ 为燃料电池产生的电能, kW； ${q}_{\text{sens }}$ 为燃料电池堆中流体的温度转化为标准温度时需要的能量,即显热, kW； ${q}_{\text{latent }}$ 为燃料电池堆中水相变产生的能量,即潜热, $\mathrm{{kW}};{q}_{\mathrm{{cool}}}$ 为冷却液带走的电堆中产生的热量, kW。

${q}_{\text{theo }}$ 的计算式为

(2)${q}_{\text{theo }}= \frac{{W}_{{\mathrm{H}}_{2},\text{ reacted }}}{{M}_{{\mathrm{H}}_{2}}}\Delta {H}_{\mathrm{{rxn}}}$

式中: ${W}_{{\mathrm{H}}_{2},\text{ reacted }}$ 为电堆中反应消耗氢气的质量流量, $\mathrm{{kg}}/\mathrm{s};{M}_{{\mathrm{H}}_{2}}$ 为氢气的摩尔质量, $\mathrm{{kg}}/\mathrm{{mol}};\Delta {H}_{\mathrm{{rxn}}}$ 为电化学反应中氢气反应时的焓变,标准条件下,该值为 285.8 kJ/mol。

${q}_{\text{elec }}$ 的计算式为

(3)${q}_{\text{elec }}= \frac{n{V}_{\mathrm{{fc}}}I}{1000}$

式中: $n$ 为燃料电池堆中单电池的数量; ${V}_{\mathrm{{fc}}}$ 为燃料电池堆的输出电压, $\mathrm{V};I$ 为燃料电池堆的负载电流, A。

燃料电池堆中阳极流体的显热 ${q}_{\mathrm{{sens}},\mathrm{{an}}}$ 和阴极流体的显热 ${q}_{\mathrm{{sens}},\mathrm{{ca}}}$ 的计算式分别为

(4)${q}_{\text{sens, an }}= \left({{W}_{{\mathrm{H}}_{2},\text{ an, out }}{C}_{\mathrm{p},{\mathrm{H}}_{2},\mathrm{\;g}}+ {W}_{\mathrm{w},\mathrm{g},\text{ an, out }}{C}_{\mathrm{p},{\mathrm{H}}_{2}\mathrm{O},\mathrm{g}}+ {W}_{\mathrm{w},1,\text{ an, out }}.}\right.\\\left.{C}_{\mathrm{p},{\mathrm{H}}_{2}\mathrm{O},1}\right)\left({{T}_{\mathrm{{an}},\text{ out }}- {T}_{\text{atm }}}\right)- \left({{W}_{{\mathrm{H}}_{2},\mathrm{{an}},\mathrm{{in}}}{C}_{\mathrm{p},{\mathrm{H}}_{2},\mathrm{\;g}}+ {W}_{\mathrm{w},\mathrm{g},\mathrm{{an}},\mathrm{{in}}}.}\right.\\{C}_{\mathrm{p},{\mathrm{H}}_{2}\mathrm{O},\mathrm{g}}+ {W}_{\mathrm{w},1,\mathrm{{an}},\text{ in }}{C}_{\mathrm{p},{\mathrm{H}}_{2}\mathrm{O},1})\left({{T}_{\mathrm{{an}},\text{ in }}- {T}_{\text{atm }}}\right)$

(5)${q}_{\text{sens, ca }}= \left({{W}_{{\mathrm{O}}_{2},\text{ca, out }}{C}_{\mathrm{p},{\mathrm{O}}_{2},\mathrm{\;g}}+ {W}_{\mathrm{w},\mathrm{g},\text{ ca, out }}{C}_{\mathrm{p},{\mathrm{H}}_{2}\mathrm{O},\mathrm{g}}+ {W}_{\mathrm{w},1,\text{ ca, out }}.}\right.\\{C}_{\mathrm{p},{\mathrm{H}}_{2}\mathrm{O},1}+ {W}_{\mathrm{w},{\mathrm{N}}_{2},\text{ ca, out }}{C}_{\mathrm{p},{\mathrm{N}}_{2},\mathrm{\;g}})\left({{T}_{\mathrm{{ca}},\text{ out }}- {T}_{\text{atm }}}\right)- \left({{W}_{{\mathrm{O}}_{2},\text{ ca, in }}\cdot }\right.\\{C}_{\mathrm{p},{\mathrm{O}}_{2},\mathrm{\;g}}+ {W}_{\mathrm{w},\mathrm{g},\mathrm{{ca}},\text{ in }}{C}_{\mathrm{p},{\mathrm{H}}_{2}\mathrm{O},\mathrm{g}}+ {W}_{\mathrm{w},1,\mathrm{{ca}},\text{ in }}{C}_{\mathrm{p},{\mathrm{H}}_{2}\mathrm{O},1}+ \\\left.{{W}_{{\mathrm{N}}_{2},\text{ in }}{C}_{\mathrm{p},{\mathrm{N}}_{2},\mathrm{\;g}}}\right)\left({{T}_{\mathrm{{ca}},\text{ in }}- {T}_{\text{atm }}}\right)$

式中: ${W}_{{\mathrm{H}}_{2},\text{ an, out }}$ 为阳极出口氢气质量流量, $\mathrm{{kg}}/\mathrm{s}$ ; ${W}_{\mathrm{w},\mathrm{g},\mathrm{{an}},\mathrm{{out}}}$ 为阳极出口水蒸气质量流量, $\mathrm{{kg}}/\mathrm{s}$ ; ${W}_{\mathrm{w},1,\mathrm{{an}},\mathrm{{out}}}$ 为阳极出口液体水质量流量, $\mathrm{{kg}}/\mathrm{s}$ ; ${W}_{{\mathrm{H}}_{2},\mathrm{{an}},\mathrm{{in}}}$ 为阳极进口氢气质量流量, $\mathrm{{kg}}/\mathrm{s};{W}_{\mathrm{w},\mathrm{g},\mathrm{{an}},\mathrm{{in}}}$ 为阳极进口水蒸气质量流量, $\mathrm{{kg}}/\mathrm{s};{W}_{\mathrm{w},1,\mathrm{{an}},\mathrm{{in}}}$ 为阳极进口液体水质量流量, $\mathrm{{kg}}/\mathrm{s};{C}_{\mathrm{p},{\mathrm{H}}_{2},\mathrm{\;g}}$ 为氢气比热容, $\mathrm{{kJ}}/\left({\mathrm{{kg}}\cdot \mathrm{K}}\right);{C}_{\mathrm{p},{\mathrm{H}}_{2}\mathrm{O},\mathrm{g}}$ 为水蒸气比热容, $\mathrm{{kJ}}/\left({\mathrm{{kg}}\cdot \mathrm{K}}\right)$ ; ${C}_{\mathrm{p},{\mathrm{H}}_{2}\mathrm{O},1}$ 为液态水比热容, $\mathrm{{kJ}}/\left({\mathrm{{kg}}\cdot \mathrm{K}}\right);{T}_{\mathrm{{an}},\text{ out }}$ 为阳极出口温度, $\mathrm{K};{T}_{\mathrm{{an}},\mathrm{{in}}}$ 为阳极进口温度, $\mathrm{K};{T}_{\mathrm{{atm}}}$ 为环境温度, $\mathrm{K};{W}_{{\mathrm{O}}_{2},\mathrm{{ca}},\text{ out }}$ 为阴极出口氧气质量流量, $\mathrm{{kg}}/\mathrm{s}$ ; ${C}_{\mathrm{p},{\mathrm{O}}_{2},\mathrm{\;g}}$ 为氧气比热容, $\mathrm{{kJ}}/\left({\mathrm{{kg}}\cdot \mathrm{K}}\right);{W}_{\mathrm{w},\mathrm{g},\mathrm{{ca}},\mathrm{{out}}}$ 为阴极出口水蒸气质量流量, $\mathrm{{kg}}/\mathrm{s};{W}_{\mathrm{w},1,\mathrm{{ca}},\text{ out }}$ 为阴极出口液体水质量流量, $\mathrm{{kg}}/\mathrm{s};{\mathrm{W}}_{\mathrm{W},{\mathrm{N}}_{2},\mathrm{{ca}},\mathrm{{out}}}$ 为阴极出口氮气质量流量, $\mathrm{{kg}}/\mathrm{s};{\mathrm{C}}_{\mathrm{p},{\mathrm{N}}_{2},\mathrm{g}}$ 为氮气比热容, $\mathrm{{kJ}}/\left({\mathrm{{kg}}\cdot \mathrm{K}}\right)$ ; ${W}_{{\mathrm{O}}_{2},\mathrm{{ca}},\mathrm{{in}}}$ 为阴极进口氧气质量流量, $\mathrm{{kg}}/\mathrm{s};{W}_{\mathrm{w},\mathrm{g},\mathrm{{ca}},\mathrm{{in}}}$ 为阴极进口水蒸气质量流量, $\mathrm{{kg}}/\mathrm{s};{W}_{\mathrm{w},1,\mathrm{{ca}},\mathrm{{in}}}$ 为阴极进口液体水质量流量, $\mathrm{{kg}}/\mathrm{s};{\mathrm{W}}_{{\mathrm{N}}_{2},\text{ in }}$ 为阴极进口氮气质量流量, $\mathrm{{kg}}/\mathrm{s}$ ； ${T}_{\mathrm{{ca}},\text{ out }}$ 为阴极出口温度, $\mathrm{K}$ ； ${T}_{\mathrm{{ca}},\text{ in }}$ 为阴极进口温度, $\mathrm{K}$ 。

${q}_{\text{sens }}$ 的计算式为

(6)${q}_{\text{sens }}= {q}_{\text{sens, an }}+ {q}_{\text{sens, ca }}$

阳极中水蒸发产生的潜热量 ${q}_{\text{latent, an }}$ 的计算式为

(7)${q}_{\text{latent },\mathrm{{an}}}= \frac{\left({{W}_{\mathrm{w},\mathrm{g},\mathrm{{an}},\mathrm{{out}}}- {W}_{\mathrm{w},\mathrm{g},\mathrm{{an}},\mathrm{{in}}}+ {W}_{\text{membrane }}}\right){H}_{\text{vaporization }}}{{M}_{\text{vaporization }}}$

式中: ${W}_{\text{membrane }}$ 为通过质子交换膜的水流量, $\mathrm{{kg}}/\mathrm{s}$ ; ${H}_{\text{vaporization }}$ 为水的摩尔焓, $\mathrm{{kJ}}/\mathrm{{mol}};{M}_{\text{vaporization }}$ 为水的摩尔质量, kg/mol。

阴极中的潜热量的计算要比阳极复杂, 因为阴极中存在水的生成和水的相变, 如果阴极中水蒸气达到了饱和状态,开始有液态水的存在,即液态水质量大于 0 ,则阴极中水蒸发产生的潜热量 ${q}_{\text{latent },\text{ ca }}$ 的计算式为

(8)${q}_{\text{latent },\mathrm{{ca}}}= \frac{\left({{W}_{\mathrm{w},\mathrm{g},\mathrm{{ca}},\mathrm{{out}}}- {W}_{\text{membrane }\mathrm{e}}+ {W}_{\mathrm{w},\mathrm{g},\mathrm{{ca}},\mathrm{{in}}}}\right){H}_{\text{vaporization },\mathrm{c}1}}{{M}_{\text{vaporization }}}\text{(}$

否则, 阴极中反应生成的水和流入阴极的液态水均会蒸发为水蒸气, 阴极中水蒸发产生的潜热量 ${q}_{\mathrm{{latent}},\mathrm{{ca}}}$ 的计算式为

(9)${q}_{\text{latent },\mathrm{{ca}}}= \left\lbrack {{W}_{\mathrm{w},1,\mathrm{{ca}},\text{ in }}{H}_{\text{vaporization },\mathrm{{cl}}}+ \left({{W}_{\mathrm{w},\mathrm{g},\mathrm{{ca}},\text{ out }}- {W}_{\text{membrane }}- }\right.}\right.\\\left.{{W}_{\mathrm{w},\mathrm{g},\mathrm{{ca}},\mathrm{{in}}}- {W}_{\mathrm{w},1,\mathrm{{ca}},\mathrm{{in}}}){H}_{\text{vaporization },\mathrm{c}2}}\right\rbrack /{M}_{\text{vaporization }}$

水的摩尔焓 ${H}_{\text{vaporization }}$ 的计算式为

(10)${H}_{\text{vaporization }}= {45070}- {41.9T}+ {3.44}\times {10}^{-3}{T}^{2}+ \\{2.54}\times {10}^{-6}{T}^{3}- {8.98}\times {10}^{-{10}}{T}^{4}$

式(8)~(10)中:下标 c1 和 c2 代表水的不同状态, 也对应着不同温度,即计算 ${H}_{\text{vaporization, cl }}$ 时, $T$ 取环境温度 ${T}_{\mathrm{{atm}}}$ ; 计算 ${H}_{\text{vaporization, c }2}$ 时, $T$ 取电堆温度 ${T}_{\mathrm{{fc}}}$ 。

${q}_{\text{latent }}$ 的计算式为

(11)${q}_{\text{latent }}= {q}_{\text{latent, an }}+ {q}_{\text{latent, ca }}$

${q}_{\text{cool }}$ 的计算式为

(12)${q}_{\mathrm{{cool}}}= {W}_{\mathrm{w},\text{ in }}{C}_{\mathrm{p},{\mathrm{H}}_{2}\mathrm{O},1}\left({{T}_{\mathrm{{fc}}}- {T}_{\mathrm{w},\text{ in }}}\right)$

(13)${m}_{\mathrm{{st}}}{C}_{\mathrm{p},\mathrm{{st}}}\frac{d{T}_{\mathrm{{fc}}}}{dt}= {q}_{\text{theo }}- {q}_{\text{elec }}- {q}_{\text{sens }}- {q}_{\text{latent }}- {q}_{\text{cool }}$

式中: ${W}_{\mathrm{w},\mathrm{{in}}}$ 为冷却液流量, $\mathrm{{kg}}/\mathrm{s};{T}_{\mathrm{w},\mathrm{{in}}}$ 为冷却液进堆温度, $\mathrm{K};{m}_{\mathrm{{st}}}$ 为燃料电池电堆的质量, $\mathrm{{kg}};{C}_{\mathrm{p},\mathrm{{st}}}$ 为燃料电池电堆的平均比热容,取 ${35}\mathrm{\;{kJ}}/{\left(\mathrm{{kg}}\cdot \mathrm{K}\right)}^{\left\lbrack {11}\right\rbrack }$ 。

至此,本文所搭建的温度模型如图 3 所示。

2.3 热管理系统模型及参数设置

采用 MATLAB/Simulink/Simscape 搭建燃料电池堆的热管理系统模型( 图 4 ),构建其散热器、冷却水泵、冷却水箱等关键模块( 图 5 )。

图 5(a)为管道模块, 可以表示流体通过管道壁面的传热情况,其中, $\mathrm{A}$ 为入口, $\mathrm{B}$ 为出口, $\mathrm{H}$ 为管壁换热接口。在热管理系统模型中, 使用该模块模拟散热器管道。图 5(b)为控制质量流量源模块, 可以表示热流体系统中理想的机械能量来源, 其中,流体从 $\mathrm{A}$ 端口流向 $\mathrm{B}$ 端口为正方向; $\mathrm{M}$ 为控制信号, 控制冷却液通过的质量流量大小。该模块不考虑进出口的压差, 并且假设其与流体之间没有热交换。在本模型中, 使用该模块模拟冷却水泵。图 5(c)为平移机械转换器模块,模拟了流体和机械平移系统之间的界面, 可以用作线性执行器的构件, 该转换器装有可变容积的液体, 温度的变化取决于液体的热容, 将流体动态可压缩性设置为打开模式后,压力也会根据液体体积的动态可压缩性来变化。该模块可以用来模拟热管理系统中的冷却水箱,本模型设定水箱横截面积为 0.1 ${\mathrm{m}}^{2}$ ,初始冷却液液位为 ${0.07}\mathrm{\;m}$ 。

燃料电池的发电效率为 40%~60%,其余的能量是以热能的形式散出去的。由式(3)计算得到所搭建的燃料电池堆模型可以发出 ${36.8}\mathrm{\;{kW}}$ 的电能,假设燃料电池的发电效率为 ${50}\%$ ,则燃料电池可以同时产生大约 ${36.8}\mathrm{\;{kW}}$ 的废热。用于散热的冷却液流量的计算式为

(14)${W}_{\mathrm{w},\mathrm{{in}}}= \frac{{q}_{\mathrm{{cool}}}}{{C}_{\mathrm{p},\mathrm{w}}{\Delta T}}$

式中: ${\Delta T}$ 为冷却液的进出堆温差, ${}^{\circ }\mathrm{C}$ 。

当 ${\Delta T}$ 为 ${10}^{\circ }\mathrm{C}$ 时, ${W}_{\mathrm{w},\mathrm{{in}}}$ 约为 ${0.88}\mathrm{\;{kg}}/\mathrm{s}$ ; 当 ${\Delta T}$ 为 ${5}^{\circ }\mathrm{C}$ 时, ${W}_{\mathrm{w},\mathrm{{in}}}$ 约为 ${1.75}\mathrm{\;{kg}}/\mathrm{s}$ 。

热管理系统模型中选用翅片式散热器,散热器管道数目 $n$ 与散热管的横切面积和冷却液流量有关 ^{[ 12 ]} ,其计算式为

(15)$ n =\frac{{1000}{W}_{\mathrm{w},\text{ in }}}{{\rho }_{\mathrm{w}}{vS}}$

式中: ${\rho }_{\mathrm{w}}$ 为冷却液的密度,选用水作为冷却液,则该值为 ${1000}\mathrm{\;{kg}}/{\mathrm{m}}^{3};v$ 为散热器管道中冷却液的流速,参照文献[ 13 ],本文取 $1\mathrm{m}/\mathrm{s};S$ 为散热器管道的横截面, 模型中散热器管道横截面的尺寸为 ${1.5}\mathrm{\;{mm}}\times {25}{\mathrm{\;{mm}}}_{\circ }$

根据式 (15), 计算得出该热管理系统的散热器管道数目应大于 47 ,本模型选取 50 。

热侧流体的总流通面积 ${A}_{\mathrm{w}}$ 为

(16)${A}_{\mathrm{w}}= S \times n ={1875}{\mathrm{\;{mm}}}^{2}$

根据上文确定的散热器尺寸, 求得当量直径 ${d}_{\mathrm{e}}$ 为

(17)${d}_{\mathrm{e}}= 4 \times \frac{A}{{l}_{\mathrm{p}}}= 4 \times \frac{S}{4 \times \left({{W}_{\mathrm{{rad}}}+ {H}_{\mathrm{{rad}}}}\right)} ={2.8}\mathrm{\;{mm}}$

考虑空气侧的换热情况,空气的体积流量 ${V}_{\text{air }}$ 的计算式为

(18)${V}_{\mathrm{{air}}}= \frac{{q}_{\mathrm{{air}}}}{{\rho }_{\mathrm{{air}}}{C}_{\mathrm{p},\mathrm{{air}}}\Delta {T}_{\mathrm{{air}}}}$

式中: ${q}_{\mathrm{{air}}}$ 为散热器散发到空气中的热量, $\mathrm{{kW}}$ ,由于要保持热平衡,该值与 ${q}_{\text{cool }}$ 相等; ${\rho }_{\text{air }}$ 为空气密度,取 ${1.2}\mathrm{\;{kg}}/{\mathrm{m}}^{3};\Delta {T}_{\text{air }}$ 为空气与散热器换热前后的温差,℃。

空气与散热器换热前的初始温度为 ${25}^{\circ }\mathrm{C}$ ,与散热器换热后的温度取散热器空气入口温度和出口温度的平均值 $\left({{52.5}^{\circ }\mathrm{C}}\right)$ ,则 $\Delta {T}_{\text{air }}$ 的计算结果为 27.5 °C。

冷却空气与散热器换热的面积 ${A}_{\text{rad, air }}$ 的计算式为

(19)${A}_{\text{rad, air }}= \frac{{V}_{\text{air }}}{{v}_{\text{air }}}$

式中: ${v}_{\mathrm{{air}}}$ 为空气流速,取 ${2.5}\mathrm{\;m}/\mathrm{s}$ 。

由式 (18),(19) 可计算出 ${A}_{\text{rad, air }}$ 为 ${0.44}{\mathrm{\;m}}^{2}$ 。

2.4 热管理系统控制

PEMFC 热管理系统可使燃料电池维持在适宜的工作温度, 从而使燃料电池具有较好的输出性能。本文采用传统的 PID 控制方法对 PEMFC 热管理系统进行控制, PID 控制被广泛应用于工业过程控制,具有算法简单、鲁棒性好和可靠性高等优点。

PEMFC 热管理系统的具体控制方法:在热管理系统模型中加入温度传感器模块, 该模块可以测量电堆温度, 并将其与设定的电堆目标温度进行比较, PID 根据得出的温度差值进行运算调节处理,得到经过冷却水泵的冷却液流量,从而通过改变冷却液流量来调节电堆的工作温度。因此, PID 控制器的输入量为电堆的目标工作温度 ${T}_{\text{set }}$ 与测量的工作温度 ${T}_{\text{meas }}$ 的差值 $e$ ,输出量为冷却水泵的冷却液流量 ${W}_{\mathrm{w},\mathrm{{in}}}$ 。PID 控制框图如图 6 所示。

将控制模型在 Simulink 下进行仿真, 可得到燃料电池电堆温度在该 PID 控制下的响应曲线, 如图 7 所示。

由图 7 可以看出:采用 PID 控制燃料电池热管理系统后,当系统的负载电流为 ${50}\mathrm{\;A}$ 时,电堆可以在 ${0.9}\mathrm{\;s}$ 左右达到设定的目标温度 $\left({{80}^{\circ }\mathrm{C}}\right)$ ,之后上升到 ${82}^{\circ }\mathrm{C}$ 并在约 ${100}\mathrm{\;s}$ 时逐渐稳定到 ${80}^{\circ }\mathrm{C}$ , 在达到稳定状态后,振荡几乎消除；当负载电流由 ${50}\mathrm{\;A}$ 增大到 ${80}\mathrm{\;A}$ 时以及由 ${80}\mathrm{\;A}$ 增大到 ${100}\mathrm{\;A}$ 时, 在 PID 控制下, 电堆温度突增后瞬间回到目标值, 表明 PID 控制方法的控制效果很好; 当负载电流从 ${100}\mathrm{\;A}$ 下降到 ${60}\mathrm{\;A}$ 时,电堆温度出现突降后同样稳定上升到设定的目标值。综上可知, PID 控制方法可以将燃料电池电堆温度维持在目标值 $\pm 2{}^{\circ }\mathrm{C}$ ,同时能对负载波动快速响应,使电堆温度迅速达到稳定。

2.5 模型验证

为了验证所搭建模型的合理性与准确性, 将燃料电池堆模型和热管理系统模型的模拟数据与实验数据进行对比验证。针对燃料电池堆模型的验证, 采用文献[ 14 ]中设置的运行参数, 即环境温度为 ${25}^{\circ }\mathrm{C}$ ,电堆工作温度为 ${72}^{\circ }\mathrm{C}$ ,氢气供给压力为 ${3.039}\times {10}^{5}\mathrm{\;{Pa}}$ ,空气供给压力为 ${3.039}\times {10}^{5}\mathrm{\;{Pa}}$ ,过氧比为 2。将燃料电池堆模型得出的极化曲线与文献[ 14 ]中的实验结果进行对比验证, 结果如图 8 所示。

由图 8 可以看出: 本文所搭建模型得出的极化曲线与实验得出的极化曲线基本一致, 当电流密度为 $0 \sim {0.05}\mathrm{\;A}/{\mathrm{{cm}}}^{2}$ 时,模型输出的电压略低于实验值,其中最大误差为 ${0.023}\mathrm{\;V}$ ,相对误差约为 2.1%；当电流密度为 ${0.05}\sim {0.8}\mathrm{\;A}/{\mathrm{{cm}}}^{2}$ 时,模型输出的电压略高于实验值, 其中最大误差约为 0.036 V,相对误差约为 ${4.6}\%$ ; 当电流密度为 ${0.8}\sim {1.0}$ $\mathrm{A}/{\mathrm{{cm}}}^{2}$ 时,模型输出的电压略低于实验值,其中最大误差约为 ${0.014}\mathrm{\;V}$ ,相对误差约为 ${2.5}\%$ 。综上所述, 电堆电压的仿真结果与实验数据的相对误差在 5%以内,保证了本文所建模型的合理性和准确性。

针对热管理系统模型的验证, 首先将模型中冷却液的流量设定为与文献[ 15 ]中一样, 为一恒定值,且负载电流设定为从 ${80}\mathrm{\;A}$ 逐渐上升到 200 $\mathrm{A}$ ,步长为 ${40}\mathrm{\;A}$ ,对比模型与实验的冷却水进、出堆温差结果见图 9 。由图 9 可知:随着负载电流的增加, 冷却水进、出堆温差也随之增加, 这是由于电堆电流的增加会导致电池发电时产热量的增加,从而影响冷却水进、出电堆的温度；冷却水进、出堆温差的模型输出结果与实验结果的变化趋势一致。当负载电流为 ${80}\mathrm{\;A}$ 时,最大误差约为 ${0.212}{}^{\circ }\mathrm{C}$ ,相对误差约为 ${3.8}\%$ ；当负载电流为 120 A 时,最大误差约为 ${0.544}{}^{\circ }\mathrm{C}$ ,相对误差约为 ${8.7}\%$ ; 当负载电流为 ${160}\mathrm{\;A}$ 时,最大误差约为 ${0.240}{}^{\circ }\mathrm{C}$ ,相对误差约为 ${3.0}\%$ ；当负载电流为 200 A 时,最大误差约为 ${0.273}{}^{\circ }\mathrm{C}$ ,相对误差约为 2.7%。综上所述,冷却水进、出堆温差的仿真结果与实验数据的相对误差大部分在 5% 以内, 说明模型可以有效反映燃料电池的温度特性, 能够较好地模拟电堆的产热情况。

3 热管理系统故障仿真

收起

燃料电池系统中的热管理系统主要是通过控制冷却液进入电堆流量大小的方式使燃料电池电堆保持在最佳工作温度,进而使其输出性能最佳。热管理系统是燃料电池系统中重要的子系统,在实际情况中, 系统中各设备由于使用时间过长而出现老化或损坏(直接失效),由此可能出现冷却液出堆温度变高或冷却液压力降低等现象,直接导致燃料电池输出性能的降低。在稳态工况下, 本文对热管理系统中可能出现的典型故障进行仿真和分析。

3.1 散热器风扇故障

图 10 为散热器风扇故障模式下燃料电池的性能曲线。在热管理系统中, 通过控制散热器风扇使冷却液降温从而冷却电堆,因此,本文通过改变散热器风扇效率参数来模拟散热器风扇的健康状态。当效率为 ${80}\%$ 时,散热器风扇为正常工作状态；当效率降为 30% 时,散热器风扇已老化使效率降低; 当效率为 0 时, 散热器风扇出现了严重故障,已完全失效。由图 10 (a) 可知:随着散热器风扇效率的降低, 冷却液出堆温度升高, 当散热器风扇效率为 ${80}\%$ 时,经过约 ${40}\mathrm{\;s}$ ,冷却液出堆温度趋于稳定值 ${346}\mathrm{\;K}$ ; 当效率降为 ${30}\%$ 时,经过约50 s,冷却液出堆温度趋于稳定值 ${353}\mathrm{\;K}$ ;当散热器风扇完全故障后 (效率为 0 ),冷却液出堆温度升高到 ${372}\mathrm{\;K}$ ,远高于 PEMFC 的最佳工作温度 353 ${\mathrm{K}}_{\circ }$ 从图 10 (b)可以看出:当散热器风扇发生故障, 热管理系统效果变差, 膜中含水量也随之下降; 当散热器风扇完全故障后 (效率为 0 ), 膜中含水量的变化不再是一条稳定曲线,而是从 ${15}\mathrm{\;s}$ 左右开始,从 14%急剧下降到 7%左右,膜从初始的湿润状态下降到干燥状态,出现了膜干现象。膜中含水量的多少直接会影响到膜的传导率,即影响燃料电池内质子穿过交换膜的能力,因此,燃料电池的输出电压也出现了下降,在大约 ${50}\mathrm{\;s}$ 时,电压从 ${352}\mathrm{\;V}$ 左右下降到了 ${345}\mathrm{\;V}$ 左右[ 图 10 (c)]。

3.2 冷却液流量不足

图 11 为冷却液流量不足模式下燃料电池的性能曲线。

在热管理系统中, 造成冷却液出现流量不足的原因有:管道破裂、管道堵塞、水箱破裂、水箱缺水、水泵出现故障等。通过改变冷却液流量对以上故障进行分析, 冷却液流量也间接表示了故障严重程度,选取充足的冷却液流量为 ${1.5}\mathrm{\;{kg}}/\mathrm{s}$ ,表示各设备为正常状态；冷却液流量为 ${1.0}\mathrm{\;{kg}}/\mathrm{s}$ ,表示系统出现了轻度故障；冷却液流量下降到 ${0.5}\mathrm{\;{kg}}/\mathrm{s}$ , 表示系统出现了重度故障。由图 11 可知:随着冷却液流量的降低,冷却液压力从 ${1.29}\times {10}^{5}\mathrm{\;{Pa}}$ 下降到了 ${1.13}\times {10}^{5}\mathrm{\;{Pa}}$ ; 冷却液流量不足也会造成冷却液出堆温度过高,当冷却液流量为 ${0.5}\mathrm{\;{kg}}/\mathrm{s}$ 时,冷却液出堆温度一直在上升,约 ${80}\mathrm{\;s}$ 时,温度趋于平稳,此时温度已上升到 ${370}\mathrm{\;K}$ ,超过了燃料电池的最佳工作温度,从而影响燃料电池的输出性能, 导致输出电压的下降。

4 结论

收起

本文以 PEMFC 系统为研究对象, 将数学机理模型和物理模型集成构建了完整的 PEMFC 系统仿真模型, 并基于该模型对 PEMFC 热管理系统进行了故障仿真研究。

① 将基于 MATLAB/Simulink 软件搭建的 PEMFC 电堆数学机理模型和基于 Simscape 平台搭建的热管理系统物理模型集成了 PEMFC 系统仿真模型, 经与实验结果进行比较, 验证了模型的合理性。

②本文搭建的热管理系统模型可以进行故障仿真。当热管理系统的散热器风扇发生故障后,随着故障程度的加深, 冷却液出堆温度逐渐升高, 远高于最佳工作温度,同时,燃料电池电堆的交换膜的含水量快速下降,出现了膜干现象,导致燃料电池的输出电压急剧下降。随着热管理系统中冷却液流量的不足,冷却液压力会有显著的下降,同样会造成冷却液出堆温度过高而影响燃料电池的输出性能。

③从仿真结果可知,可以通过监测冷却液压力、冷却液出堆温度、电堆电压等参数及时了解热管理系统的工作状态,及时对其进行故障诊断。

基金

收起

北京市自然科学基金面上项目(3222003)
国家自然科学基金青年项目(52002009)
汽车安全与节能国家重点实验室开放基金课题(KF2010)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

Song

, Wang

, Ding

Y H

, et al. Assembly techniques for proton exchange membrane fuel cell stack: A literature review[J]. Renewable and Sustainable Energy Reviews, 2022, 153: 111777.

[2]

, Li

, Gao

, et al. Design and simulation of proton exchange membrane fuel cell system[J]. Energy Reports, 2021, 7(6): 522-530.

[3]

Omran

, Lucchesi

, David

, et al. Mathematical model of a proton exchange membrane (PEM) fuel cell[J]. International Journal of Thermofluids, 2021, 11: 100110.

[4]

A D

, Zhou

. A general model of proton exchange membrane fuel cell[J]. Journal of Power Sources, 2008, 182(1): 197-222.

[5]

房鑫, 张洪伟, 隋宗强, 等. 车用质子交换膜燃料电池系统建模[J]. 青岛大学学报(工程技术版), 2018, 33(1): 6-10.

[6]

郭健忠, 王家欣, 孙凯华, 等. 基于 S-Function 的 PEM 燃料电池建模仿真分析[J]. 中国科技论文, 2020, 15(12): 1364-1370.

[7]

卫超强, 武志斐, 蒋栋. 温度对质子交换膜燃料电池输出性能的影响[J]. 可再生能源, 2021, 39(12): 1588-1593.

[8]

石磊, 郑明刚, 孔祥利. 8通道复合蛇形流道质子交换膜燃料电池综合优化研究[J]. 可再生能源, 2020, 38(4): 434-439.

[9]

曹新宇. 基于Simscape 的超临界CO₂换热器建模仿真研究[J]. 装备制造技术, 2021(6): 57-59.

[10]

Pukrushpan

J T

. Modeling and control of fuel cell systems and fuel processors[D]. Michigan: The University of Michigan, 2003.

[11]

X C

, Zhou

, Sobiesiak

. Water and thermal management for Ballard PEM fuel cell stack[J]. Journal of Power Sources, 2005, 147(1-2): 184-195.

[12]

赵兴强. 水冷型质子交换膜燃料电池热管理系统研究[D]. 成都: 西南交通大学, 2015.

[13]

蒋诚. 车用燃料电池热管理系统设计开发与性能验证[D]. 吉林: 吉林大学, 2020.

[14]

展茂胜. 质子交换膜燃料电池热管理系统的优化与控制[D]. 济南: 山东大学, 2020.

[15]

史青. 水冷型 PEMFC 热管理系统建模与控制研究[D]. 成都: 西南交通大学, 2017.

2024年第42卷第3期

PDF下载

253

126

引用本文

BibTeX

文章信息

接收时间：2022-04-18
首发时间：2025-07-22
出版时间：2024-03-20

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2022-04-18

基金

北京市自然科学基金面上项目(3222003)

国家自然科学基金青年项目(52002009)

汽车安全与节能国家重点实验室开放基金课题(KF2010)

作者信息

¹ 北京工业大学交通工程系北京 100124

² 北京工业大学汽车工程系北京 100124

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/kzsny/CN/1154429038199296928

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT