固体氧化物燃料电池数值仿真模型的建立

2019-04-17汤水张茂贵贾明刘宇

中南大学学报（自然科学版） 2019年11期

汤水，张茂贵，贾明，刘宇

(中南大学冶金与环境学院，湖南长沙，410083)

近年来，我国资源与安全环境问题的矛盾日益凸显。在电动汽车大力发展的推动下，燃料电池作为一项清洁高效的尖端技术成为人们研究的焦点。燃料电池是继火电、水电和核电后的第四代发电方式[1]，发电不受卡诺循环限制，能量转化效率高达40%～60%，可直接将燃料与氧化剂中储存的化学能转化为电能。以氢氧燃料电池为例[2]，在电池阳极侧通入氢气，在电池阴极侧通入氧气，在阴极和阳极的电催化剂的帮助下，氧气在阴极反应区发生还原反应生成水，氢气在阳极反应区发生氧化反应生成氢离子和电子，通过电解质传导离子及电子沿外电路定向移动形成回路，从而形成电流带动负载工作。固体氧化物燃料电池SOFC(solid oxide fuel cell)[3]具有诸多优点，研发和应用大规模的SOFC发电装置系统能有效缓解目前人类面临的能源问题。然而，固体氧化物燃料电池工作温度高，结构紧凑、复杂，工作环境封闭且复杂，通过实际实验来获得其内部状态变化情况和具体位置参数往往极为困难且耗费大量财力和物力。随着计算机计算能力的提高，搭建固体氧化物燃料电池的数学模型、使用计算机数值分析求解来进行仿真实验成为可能。 CFD(computational fluid dynamics)方法是分析燃料电池性能的一种重要方法[4−6]，它使用数值方法将描述特定流场运动的方程组离散到网格节点上并求出其离散数值解，主要有有限差分法、有限元法、有限容积法等。WEI等[7]提出了基于流量通道设计对堆叠性能影响的数值研究、流路和堆叠布置的新设计，并分析了平面阳极支撑的固体氧化物燃料电池堆的热应力。MANGLIK 等[8]在具有不同形状的双极板互连流动通道、由平面阳极支撑的SOFC中同时模拟热和质量传递。黑盒模型基于人工智能，其实现主要依赖于实验数据[9]。SORRENTINO 等[10]综合白盒和黑盒模型，提出一种SOFC单电池的灰盒模型。它基于物理关系，根据经验建模，较灵活[11]。这些研究大多以温度、流速等为研究对象，未考虑到内部多物理场耦合情况。为此，本文基于COMSOL Multiphysics数值仿真平台，考虑板式固体氧化物燃料电池内部的气体流动、组分质量分数、电化学反应过程，建立固体氧化物燃料电池的三维模型。

1 固体氧化物燃料电池模型的建立

1.1 研究对象

SOFC 工作温度高(600～1 000 ℃)，可以快速达到热力学平衡，从而降低对催化剂的要求。其以氧化物离子导电陶瓷材料作为电解质，氧离子作为传输离子，工作原理如图1所示。

氧化剂(一般为氧气或富氧空气)在阴极反应区内发生还原反应生成O2−：

O2−通过固体电解质输入至阳极反应区，而电子则通过外电路从电池的阳极传输到阴极形成电流。燃料气在阳极反应区内发生氧化还原反应，见下式(以氢气为例)：

电池的总反应为

在标准状态下，SOFC单电池电压理论上能达到1.229 V，能量转化效率高达80%。电池的开路电压可用下式计算：

图1 固体氧化物燃料电池工作示意图Fig.1 Schematic diagram of solid oxide fuel cell operation

其中：

式中：PO2，PH2和PH2O分别为氧气、氢气和水蒸气分压；R为通用气体常数；F为法拉第常数；E0为电池标准电压；ΔG0为标准吉布斯自由能。

工作温度、燃料及氧化剂的种类对单电池的电压和能量转化效率都有很大影响[12]。此外，由于电池中的欧姆极化、电化学极化和浓差极化，电池的实际输出电压与理论值有很大差异。

1.2 模型建立

1.2.1 几何模型

三维模型可以模拟燃料电池内部行为，并且给出关于温度和燃料分布信息[13]。选取固体燃料电池单体作为研究对象，建立SOFC 单电池三维模型，如图2所示，具体参数如表1所示。

图2 SOFC单气体流道模型和SOFC单电池三维模型Fig.2 SOFC single gas channel model and 3D model of SOFC single cell

1.2.2 模型假设

为简化模型、提高模型的适用型，在模型建立的过程进行如下假设：

1)本模型反应过程是稳态的，且电池内部各种材料具有各向同性；

2)温度均匀且在运行过程中保持不变；

3)压力为标准大气压；

表1 SOFC模型参数Table 1 Parameters of SOFC model

4)反应气体均为理想气体；

5)电极内物质对流传递可以忽略。

1.2.3 SOFC单电池三维模型的建立

经分析，模型的反应过程可由Butler-Volmer公式、欧姆定律(描述电子电荷平衡、离子电荷平衡)、Navier-Stokes模型(描述气体通道中的流量分布)、Strokes-Brinkman 模型(描述多孔气体扩散电极中的流量)、Maxwell-Stefan 扩散和对流模型(描述气体通道和多孔电极中气相的质量平衡)，以及边界条件和初始条件构成封闭的方程组来描述。下面详述模型所使用的物理场接口及方程[14−16]。

1.3 物理场分析

1.3.1 二次电流分布

本模型在阴极、电解质和阳极中采用二次电流分布接口，通过二次电流分布接口与模拟质量传递的接口耦合，可描述模型中电流密度分布随浓度变化的变化；通过二次电流分布接口与欧姆定律和电荷守恒耦合，可描述模型中电极和电解质中的电流传导。

将阳极和阴极定义为多孔电极，给定电解质和电极材料电导率、多孔电极反应的温度、平衡电位、相关动力学参数和活性比表面积，设置电池的零电势点和电压初始值。在离子电荷守恒中，假设所有外部边界都是绝缘的，而在内部边界将电流和电势看作是连续变化的。假设阳极和阴极反应区中的电子和离子满足电荷平衡，电荷转移和过电位之间的关系可以用Butler-Volmer 公式来描述[17]。

在阳极，氢被氧化形成水，并假设第一电子转移为速率控制步骤，则满足

式中：i0,a为阳极交换电流密度，A/m2；cH2和cH2O分别为H2和H2O浓度，mol/m3；cH2,ref和cH2O,ref分别为H2和H2O 参考浓度，mol/m3；T为温度，K；η为过电位，V。

对阴极，满足

式中：i0,c为阴极交换电流密度，A/m2；xO2为氧气的摩尔分数，%；cO2,ref为O2参考浓度，mol/m3。

过电位η用下式计算：

式中：φele为电子电阻引起的欧姆电阻电位，V；φionic为离子电阻引起的欧姆电阻，Ω；Δφeq为平衡电势差，V。

以阳极的入口边界作为零电势点，阴极的入口边界作为电池电势点，Vcell通过下式计算：

式中：Vcell为SOFC工作电压；Vpol为极化电压，V；Δφeq,c和Δφeq,a分别为阴、阳极平衡电势差。

在此模型中，Δφeq,a为0 V，Δφeq,c为1 V。当压力一定时，电池电压随电流密度增加而降低，即电压损失会逐渐增大，因电流密度增加导致极化损失和欧姆损失增大；当电流密度较低时，电池反应受电极动力学控制，电压损失以活化极化为主；当电流密度较高时，电化学反应受气体扩散速率控制，电压快速下降。据式(4)，温度升高，电池的理想电压会降低，同时，高温有利于增大电解质的离子电导率，减小欧姆极化过电位[18]。

1.3.2 浓物质传递

在阳极和阴极各设1个浓物质传递接口，阳极因变量为H2和H2O的质量分数，阴极因变量为O2，H2O和N2的质量分数。浓物质传递用于研究气体和流体混合物，其中，物质浓度具有相同量级，且都不能被看作是溶剂。混合物的属性依赖于组成，需考虑所有物质的分子和离子之间的相互作用。此接口中包含多组分扩散模型，其中每种物质的扩散驱动力受组成、温度以及压力的影响，用于求解所有参与物质的质量分数，通过对流(粒子随溶液的流动而流动)、扩散(粒子在浓度梯度下的运动)以及电迁移(带电粒子在电场作用下的迁移)来传质。假设电极和气体通道的扩散机制为Maxwell-Stefan 模型，给定方程相关参数和初始值，加入多孔电极耦合，输入局部电流源，给定反应电子数和化学计量数，指定气体出口和入口。将气体通道和气体扩散电极壁的边界条件设为零质量通量，设组分和通量的连续性适用于气体扩散电极和气体通道所有界面。

1.3.3 自由和多孔介质流动

在阳极和阴极各设1个自由和多孔介质流动接口，阳极因变量为速度场ua、压力pa，阴极因变量为速度场uc、压力pc。本模型使用弱可压缩Navier-Stokes 方程控制气体通道中的气体，用Brinkman 方程描述多孔气体扩散电极中的流速，计算单相自由流动与多孔介质连接在一起时的速度场和压力场。

将电极和气体通道的自由和多孔介质流动的物理模型设为可压缩流(马赫数Ma<0.3)，给定流体的动力黏度和密度；设定多孔基体属性，给定孔隙率、渗透率；加入多孔电极耦合，输入局部电流源，给定反应电子数、化学计量数。指定流体出口和入口，设置压力，使入口处轻微过压以驱动气体流动，出口处选择法向流，设置边界条件为滑移。耦合以上几个接口，可描述SOFC单电池的内部电化学反应过程和传质过程。

2 结果讨论

2.1 模型验证

确定模型的物理场、内部参数和网格划分方式，分别对极化电压、入口氧气质量分数、工作温度、单电池长度进行参数化扫描和模拟计算。考察电池内部气体分布、电流密度分布、极化情况。将模拟所得结果与实验结果进行对比，以验证模型的有效性[19-20]。该三维模型剖分网格为36 288个，满足研究要求。

2.2 不同极化电压下电池内部状态

为研究不同极化值下(即电池电压不同时)电池内部气体分布情况以及电流密度分布情况，应用上述数学模型，在温度为800 ℃，压力为1×105Pa，单电池长度为0.01 m，入口燃料气体为40%(质量分数，下同)氢气与60%水蒸气，阴极气体为15%氧气与37%水蒸气与58%氮气，其他参数见表1，计算不同入口氧气质量分数下的阴阳极气体分布情况。

图3所示为不同极化值时电解质平面上电流密度及分布。从图3可见：当极化值为0.05 V即电池电压为0.95 V 时，电流集中分布在电池左端；随着极化值不断增大，电池电流密度也随之增大；电流开始集中分布于电池右端，且随极化值进一步增大，电池密度峰值范围进一步扩大，分布更加均匀，由右端向中部逐渐减少；当电流密度较低时，电池反应受电极动力学控制，电压损失以活化极化为主；当电流密度较高时，电化学反应受气体扩散速率控制，电压快速下降；因氧气耗尽，电流密度分布很不均匀，集中在右端阴极入口处。

图3 不同极化值下电解质平面上电流密度及分布Fig.3 Current density and distribution on electrolyte plane under different polarization values

在极化值为0.5 V即电池电压为0.5 V时，氧气摩尔分数在阴极和气体流道的分布情况如图4所示。从图4 可见：氧气由电池右端向电池左端扩散，扩散非均匀。氧气在阴极的消耗量表明阴极发生反应的程度。

当极化值为0.5 V 时，氢气摩尔分数在气体流道和阳极中的分布情况如图5所示。从图5 可见：与氧气的分布类似，氢气在阳极的分布也呈现沿输入口向出口方向减少的现象。

图4 氧气的摩尔分数分布Fig.4 Mole fraction distribution of oxygen

图5 氢气的摩尔分数分布Fig.5 Mole fraction distribution of hydrogen

2.3 入口氧气质量分数对电池的影响

为研究不同入口氧气质量分数下电池内部气体分布情况以及电流密度分布情况，在温度为800 ℃，压力为1×105Pa，极化值为0.5 V，单电池长度为0.01 m，入口燃料气体为40%氢气与60%水蒸气，其他参数按表1设置的条件下，计算不同入口氧气质量分数下的阴阳极气体分布情况，分析入口氧气质量分数对极化曲线的影响，仿真结果如图6和图7所示。

从图6和表2可看出：随着入口氧气质量分数增加，出口氧气最小质量分数也随之增加，且在本模型中明显呈现线性增加的趋势。

图6 氧气分布随入口氧气质量分数的变化Fig.6 Oxygen distribution varies with inlet oxygen mass fraction

图7 氢气分布随入口氧气质量分数的变化Fig.7 Hydrogen distribution varies with inlet oxygen mass fraction

表2 入口氧气质量分数对出口氧气最小质量分数的影响Table 2 Effect of inlet oxygen mass fraction on the minimum mass fraction of outlet oxygen

图7所示为不同入口氧气质量分数时氢气质量分数在阳极和气体流道中的分布。质量分数为40%的氢气从电池左端输入，从电池右端流出，如表3所示。

随入口氧气质量分数增加，出口氢气最小质量分数也随之减少，且在本模型中呈现线性减少的趋势，但变化不大。由表3可知：当氧气质量分数为11%～23%时，氢气消耗基本相同，即在模型中入口氢气的质量只要达到40%就可保证电化学反应过程不受氢气质量分数的限制。

表3 入口氧气质量分数对出口氢气最小质量分数的影响Table 3 Influence of inlet oxygen mass fraction on the minimum mass fraction of outlet hydrogen

不同入口氧气质量分数下极化曲线如图8所示。从图8可见：在高电流密度下，入口氧气的质量分数会显著影响电池性能；而在低电流密度下，氧气质量分数不是影响电池电压的主要因素。

图8 不同入口氧气质量分数下极化曲线Fig.8 Polarization curves at different inlet oxygen mass fractions

2.4 工作温度对电池的影响

应用上述数学模型，为研究SOFC在不同工作温度(800～1 000℃)下电池电压与功率随电流密度的变化趋势，在压力为1×105Pa，单电池长度为0.01 m，燃料气体为40%氢气与60%水蒸气，阴极气体为15%氧气、37%水蒸气和58%氮气，其他参数按表1设置的条件下进行计算，仿真结果分别如图9和图10所示。

图9 不同温度下平均电流密度与工作电压的关系Fig.9 Relationship between average current density and working voltage at different temperatures

图10 不同温度下平均电流密度与功率密度的关系Fig.10 Relationship between average current density and power density at different temperatures

分析图9可知：在同一温度下，随着电流密度增大，电池工作电压不断下降。这是因为随着电流密度增大，电池极化现象越来越严重；当电流密度较低(800 A/m2以下)时，电池工作电压受温度的影响较小，而当电流密度较高时，电池工作电压受温度的影响很明显；温度越高，随着电流密度增大，电池工作电压越小，且下降速度越快。

随着温度升高，电池的活化极化、扩散极化和欧姆损失都会减少，固体电解质的离子导电率增大，气体扩散速度也会加快，电池电压反而会增高。在高温(800～1 000 ℃)下，离子电导率和扩散系数变化难以得到准确值，因而，模型将电解质的离子电导率和气体扩散系数简化为常数，不考虑温度升高对于传质的影响。据能斯特公式，处在标准状态的电池可逆电动势随温度升高而降低，比率为0.27 mV/℃[7]。模型验证了温度升高、电池工作电压降低的现象。

由图10 可知：在模拟高温下(800～1 000 ℃)，电流密度与功率密度呈抛物线型变化，表明在低电流密度时，功率密度的变化可以忽略；在同一平均电流密度下，随温度升高，电池的功率密度下降；当电流密度超过2 000 A/m2时，随温度升高，功率密度的极大值点反而朝电流密度减小的方向移动。

2.5 单电池长度对电池性能的影响

为研究不同长度时单电池内部气体分布情况以及电流密度分布情况，在温度为800 ℃，压力为1×105Pa，极化值为0.5 V，燃料气体为40%氢气与60%水蒸气，阴极气体为15%氧气、37%水蒸气和58%氮气，其他参数按表1设置的条件下，计算氧气和氢气的质量分数分布情况，仿真结果分别如图11和图12所示。

质量分数为15%的氧气从电池右端输入，左端流出，沿输入方向氧气被不断消耗。具体消耗情况如表4所示。

当单电池长度较短如2 mm时，氧气的消耗量很小，利用率不高；当长度为10 mm 时，氧气消耗较充分。

由图12 可知：当电池长度改变时，氢气质量分数分布也会随之变化。质量分数为40%的氢气从电池右端输入，左端流出，沿输入方向氢气被不断消耗。具体消耗情况如表5所示。

因模型长度较短，氢气的利用率不高，电化学反应主要受氧气质量分数的影响。