基于压力速度耦合算法的PEMFC电堆空气分布

2021-01-25罗马吉

江苏大学学报（自然科学版） 2021年1期

赵岩，罗马吉，陈奔

(1. 武汉理工大学现代汽车零部件技术湖北省重点实验室，湖北武汉 430070； 2. 武汉理工大学汽车零部件技术湖北省协同创新中心，湖北武汉 430070)

质子交换膜燃料电池(proton exchange membrane fuel cell，PEMFC)具有高能效、无污染、启动速度快和工作温度低等优点，在汽车、固定和便携式发电等领域有广泛的应用[1].单个质子交换膜燃料电池由于各种损失只能产生小于1 V的电压.损失包括活化损失、内部渗透电流损失、欧姆损失和浓度损失[2].为了有效满足车辆实际应用中的电压和功率需求，可将多个电池单元串联起来装配成质子交换膜燃料电池堆.

燃料电池的低可靠性、耐久性差和高成本严重影响燃料电池汽车的广泛商业应用[3-4].单电池性能的衰减会导致整个电堆性能的降低，这就是著名的“短板效应”[5].电堆歧管流体的不均匀分布会导致单电池在反应物饥饿情况下工作或产生水淹，这将会降低燃料电池堆的性能，加速催化剂和膜的衰减[6].

国内外对于燃料电池堆内的流体分布有诸多研究.詹志刚等[7]通过简化的电池堆模型，研究了等截面和变截面进气总管在U形和Z形2种进气方式下的电堆内部压力分布及各单电池流量分布，结果表明，具有一定锥度结构的变截面进气总管在U形进气方式时表现较好.YANG Z. R.等[8]发现增加歧管的横截面积能够改善燃料电池堆流体分布的均匀性.崔欢[9]分析了质子交换膜燃料电池电堆阴极气体分配规律，提出了气体分配不均匀度估算公式.HUANG F. X.等[10]研究了商业尺寸的质子交换膜燃料电池电堆的空气分布，结果表明：当电堆单体燃料电池数量增加时，流体分布的均匀性会降低，电堆的性能会下降.

燃料电池中的流体分布依赖于电堆进口歧管的流速和出口歧管的压力，且燃料电池流场的结构复杂，所以需对进出口歧管的特性和电堆的流体分配进行研究.笔者基于压力速度耦合算法，提出1个简化的模型，对U形和Z形结构质子交换膜燃料电池电堆空气分配进行研究，比较电堆歧管宽度和燃料电池单体数量对空气分配的影响.

1 数学模型

1.1 变量数和自由度分析

燃料电池堆中2种结构的流动模式如图1所示.

图1 数学模型中变量的指定

进口歧管体积流率和压力方程分别为

Qin(i)=Qin(i-1)+qx(i),

(1)

pin(i)=pin(i-1)+Δpin(i),

(2)

式中：Qin为进口歧管体积流率；qx为燃料电池流道体积流率；pin为进口歧管压力；Δpin为进口歧管内的压降；i为单体燃料电池序号，i=1,2，…，N，N为单体燃料电池数量.

对于Z形流场，出口歧管体积流率和压力方程分别为

Qout(i)=Qout(i+1)+qx(i),

(3)

pout(i)=pout(i-1)+Δpout(i),

(4)

式中：Qout为出口歧管体积流率；pout为出口歧管压力；Δpout为出口歧管内的压降.

对于U形流场，出口歧管体积流率和压力方程分别为

Qout(i)=Qout(i-1)+qx(i),

(5)

pout(i-1)=pout(i)+Δpout(i).

(6)

总入口体积流率Qin(N)等于U形结构中的总出口体积流率Qout(N)和Z形结构中的总出口体积流率Qout(1)，并且已给出.未知数共有5N个，共有(2N+2)个Qin、Qout、qx相关联的方程，只需要再知道N个qx的方程及(2N-2)个pin、pout和流速方程，模型就可以计算.

1.2 模型假设与边界条件

质子交换膜燃料电池电堆有着复杂的几何结构，并且建立1个完整的模型需要很多的计算资源，因此提出1个简化的模型十分必要.通过设置1个等效阻力系数将各单体燃料电池流道简化成1个直流道，流道与歧管的横截面形状都是正方形，流道的横截面积为1 mm2.等效流动阻力系数可以参考文献[11].

为了简化计算，假设如下: ① 流体为连续不可压，流动为层流； ② 不考虑流道中气体的电化学反应消耗； ③ 燃料电池流道内的流体流动过程用达西定律来描述； ④ 流体流动过程中温度恒定.

电堆的空气进口体积流率为

(7)

式中：I为电流，I=300 A;γ为过量系数，γ=2.5；R为理想气体常数；T为运行温度，T=338.15 K;n为粒子带电数，n=4 个；F为法拉第常数；φ为空气中氧的体积分数，φ=0.21；p为工作压力，p=0.202 MPa；RH为相对湿度，RH=50%；psat为水的饱和蒸汽压.

1.3 压力和流速关系

燃料电池流道内的压降Δpx用达西定律来描述，即

(8)

(9)

式中：μ为气体黏度，μ=1.9×10-5Pa·s-1；vx为燃料电池流道内气体流速;1/K为等效阻力系数，1/K=4×107m-2；Lx为直流道长度，Lx=0.1 m；Ax为流道横截面积.

电堆歧管间的压降Δp可由伯努利方程推导出来，即

(10)

(11)

式中：ρ为气体密度，ρ=1.24 kg·m-3；vy为电堆歧管内气体流速；g为重力加速度；h为高度；f为摩擦因子；Ly为相邻流道间歧管的距离，Ly=3 mm;Dh为流道的特征长度；Kc为局部阻力系数，Kc=0.5；Q为电堆歧管体积流率；Ay为歧管横截面积.

通常，重力势能项比其他项至少小1个数量级，在这里忽略不计.电堆歧管间的压降为

(12)

此时方程的未知数和方程的个数相等，自由度为0，模型已经可以计算.

1.4 压力速度耦合算法

参考文献[12]提出的求解算法，假设总气体流量平均分配到每个燃料电池内.ref为参考电池序号，对于Z形流场，ref=1；对于U形流场，ref=N.燃料电池体积流率、流道内压降及歧管间的体积流率计算方法如下：

(13)

(14)

进口歧管体积流率为

(15)

Z形流场的出口歧管体积流率为

(16)

U形流场的出口歧管体积流率为

(17)

对于Z形流场和U形流场来说，根据式(10)歧管的流速和压降的关系以及出口压力设置为0的已知条件，只需要计算出压降就可以知道pin和pout，计算方法如下：

pin(ref)=pout(ref)+Δpx(ref)，

(18)

(19)

(20)

将式(18)-(20)代入到式(2)、(4)、(6)中，可以求得pin和pout.

由于初始假设歧管流体分布均匀，因此定义参数β对qx进行修正.

(21)

当β值较大时，这意味着此时燃料电池的体积流率与歧管间的压降不匹配，因此β称为流动分布因子，为了使迭代计算更加准确，需要对初始值进行调整.调整因子为

(22)

将调整因子代入到式(14)中得新的参考燃料电池流道压降为

(23)

(24)

(25)

(26)

此时计算得到的新燃料电池体积流率为

(27)

残差为

(28)

将残差值的收敛标准设置为1×10-4,如果ERR≥1×10-4,进入下一次迭代，如果ERR<1×10-4，则模型收敛，输出qx(i)，计算结束.使用Matlab 2016b对整个算法进行编程求解.

2 结果与讨论

2.1 U形和Z形结构空气体积流率分布

燃料电池空气体积流率变化曲线如图2所示，U形结构和Z形结构的燃料电池空气分布趋势有明显的差别，在单体燃料电池数量为10个、歧管的宽度为6 mm时，Z形结构的电堆空气体积流率分布更均匀，在实际工作中为了方便安装，大多采用U形结构.

图2 燃料电池空气体积流率变化曲线

2种结构的歧管进、出口压力变化曲线如图3所示，造成2种结构空气体积流率分布趋势不同的主要原因是出口流动方向不同导致歧管出口压力分布不同.

图3 电堆歧管进、出口压力变化曲线

2.2 歧管宽度对空气体积流率的影响

歧管宽度对U形和Z形结构歧管空气体积流率的影响分别如图4、5所示.歧管宽度增大会明显改善空气体积流率分布的均匀性，这是因为空气通过歧管的压降减小，迫使燃料电池的空气分布变得均匀.当歧管宽度由6 mm增加到7 mm时(当歧管横截面的宽度改变时，歧管横截面的长度也做同样改变)，空气体积流率分布均匀性得到了较大改善，歧管宽度从7 mm增加到8 mm时，空气体积流率分布的均匀性虽有提升但变化不大.从性能角度分析，歧管的宽度不能过小，但是当空气体积流率分布差异不大时，再增加歧管宽度，效果会降低，且会增加电堆成本.

图4 歧管宽度对U形结构空气体积流率的影响

图5 歧管宽度对Z形结构空气体积流率的影响

2.3 单体燃料电池数量对空气体积流率的影响

歧管宽度为20 mm，单体燃料电池数量不同时，U形和Z形结构电堆空气体积流率变化曲线分别如图6、7所示.单体燃料电池数量从30个增加到50个，2种结构的空气分布均匀性都开始降低，但U形结构降低得更多，这会严重降低电堆性能，导致输出的平均电压变低，会对燃料电池造成较大损害.

图6 燃料电池单体数量不同时U形结构电堆空气体积流率变化曲线

图7 单体燃料电池数量不同时Z形结构电堆空气体积流率变化曲线

歧管宽度为20 mm、单体燃料电池数量为30个时，电堆空气体积流率变化曲线如图8所示.歧管宽度为30 mm、单体燃料电池数量为50个时，电堆空气体积流率变化曲线如图9所示.从图8、9可以看出：当单体燃料电池数量为30个、歧管宽度为20 mm时，空气分配较均匀;电池数量为50个时，歧管宽度需要达到30 mm，空气分配才较均匀.