目录1引言 (2)1.1燃料电池的概念及特点 (2)1.2固休氧化物燃料电池 (4)1.2.1固休氧化物燃料电池的结构类型及其特点 (4)1.2.2 SOFC工作原理 (5)2固体燃料电池多物理场模拟 (6)2.1控制方程 (6)2.1.1动量守恒方程 (6)2.1.2能量守恒方程 (6)2.1.3质量守恒方程 (6)2.1.4导电方程 (7)2.2物理模型 (7)2.3数学模型 (8)2.3.1气体输运控制方程 (8)2.3.2导电控制方程 (8)2.4边界条件 (9)3结果与讨论 (11)3.1电势分布 (12)3.2不同阳极厚度燃料电池的浓度分布 (12)3.2.1不同阳极厚度燃料电池的电势分布 (14)3.3阴极厚度对燃料电池性能影响 (15)3.4连接体宽度变化对浓度、电势分布的影响 (18)4 结论 (19)固体氧化物燃料电池仿真摘要燃料电池是将化学反应的化学能直接转变为电能的装置。

和传统的热机相比,燃料电池具有更高的电效率,并且燃料电池是一种环境友好的发电方式。

固体氧化物燃料电池(SOFC)属于高温燃料电池,除具有燃料电池的一般特点外,其高温排气也可以进一步加以利用。

本文建立了描述平板式SOFC的物理数学模型,使用多物理场耦合模拟软件Comsol对其进行模拟计算。

通过改变阳极和阴极厚度、连接体rib宽度等，研究其对固体氧化物燃料电池内燃料浓度、电势分布等的影响。

模拟结果显示：当燃料沿燃料通道方向流动未出现低燃料浓度区或产物浓度区时，电池电势在燃料流动方向上变化不大；阳极厚度的增加对反应物在垂直于燃料流动方向的分布几乎没有影响，随着阳极和阴极厚度及连接体宽度的增加，燃料电池的性能更好。

本模拟可以为燃料电池的设计提供参考。

关键词：固体氧化物燃料电池Comsol1引言随着全球工业化的加速及人们生活水平的不断提高,人类对能源的需求持续增长。

目前全球能源的大部分来自化石燃料的燃烧过程，全世界对化石燃料利用的持续增长导致了温室气体排放的增加,美国能源部预计,2015年全球的排放量要比1990年增加60%;燃料燃烧过程产生的氮氧化物,硫氧化物,未燃尽的碳氢化合物等是主要的大气污染物。

因此,解决能源需求的增长和由此造成的环境问题的关键就是改善能源结构问题,研究开发清洁能源技术。

而燃料电池技术正是符合这一需求的高效洁净能源。

1.1燃料电池的概念及特点燃料电池是把化学反应的化学能直接转化为电能的装置。

与传统的发电方式相比较,关键的区别是燃料电池的能量转化过程是直接的。

燃料电池需要清洁的燃料比如天然气作为反应气体,像煤和生物质这样的燃料可以通过部分氧化或气化转化成合成气后供燃料电池使用，与一般电池一样,燃料电池也是由阴极，阳极，电解质构成的。

燃料通入燃料电池的阳极后被氧化成阳离子并释放出电子,氧气或空气在燃料电池的阴极得到电子后形成阴离子,阴离子通过电解质向阳极扩散,并同阳极的阳离子发生电化学反应。

电子由电池的阳极通过外电路到达电池的阴极从而形成回路。

图1给出了氢一氧燃料电池的工作原理示意图图1 燃料电池的工作原理示意图虽然燃料电池和电池都是将化学能转化成电能的装置,但不同之处在于燃料电池是能量转化装置,而电池是能量存储装置。

对于电池,其化学能被存储在电池中,当电池内储存的化学反应物耗尽时,电池就不能再发出电能"而对于燃料电池,其所需的燃料和氧化剂只有当燃料电池发电时,才从外部供给燃料电池,燃料和氧化剂是独立于燃料电池本身的,只要不断向其供给燃料和氧化剂,就可以连续不断的发电。

燃料电池的种类很多,分类方法也有多种。

可以根据其温度分为低温(25100)中温(100500)和高温(5001100)燃料电池最常用的分类方法是根据电解质的性质,将燃料电池划分为五大类:碱性燃料电池(AFC),磷酸燃料电池(PAFC),质子交换膜燃料电池(PEMFC),熔融碳酸盐燃料电池(MCFC)和固体氧化物燃料电池(SOFC)表1给出了主要燃料电池极其特性表1 主要燃料电池极其特性固体氧化物燃料电池(SOFC)是八十年代发展起来的一种高温燃料电池,通常被视为固定发电厂的候选,是一种最有前途的高效，洁净发电技术之一。

SOFC的自身特点决定了它在许多领域有着广泛的应用前景,如：(1)固定式大型电站；(2)中小型分散型电站或备用电源；(3)移动式辅助电源，便携式电源；(4)军事，航天航空应用等。

1.2固休氧化物燃料电池1.2.1固休氧化物燃料电池的结构类型及其特点SOFC的结构包括几何构型和组装方法,它是电池制备工艺的基础和SOFC技术创新的主要方面。

合理的SOFC结构应具有四个基本特征,即性能可靠，便于放大，方便维修和价格低廉。

SOFC的构型不同,对其研究的内容和电池的制作方法也不同，虽然SOFC的研究过程中先后出现了多种构型的SOFC,但经过长期的发展和优化,目前采用的SOFC结构类型主要有管式，平板式和瓦楞式三种。

三种SOFC电池结构如图2，图3，图4所示,三种SOFC电池结构的优缺点如表2所示。

图2 管式SOFC电池组的结构图3 平板式SOFC电池组的结构图 4 瓦楞式SOFC电池组的结构1.2.2 SOFC工作原理SOFC与其他燃料电池的主要区别在于它采用固体氧化物电解质，其内输运的离子为氧离子。

为了减小离子在电解质的运动阻力。

一般在比较高的温度下工作(500～1000℃)。

在阴极，空气中的氧气结合外电路输运过来的电子，在电极催化作用下，形成氧离子：1+2e−→ −在电化学势的驱动下,氧离子通过固态电解质被传输至阳极。

在那里氧离子在阳极催化作用下与那里的燃料（比如H2，CO）发生电化学反应，生成相应的产物（比如水蒸气和CO2，同时伴随着电子的释放：H+ −→H+2e−+ −→+2e−产生的产物和未反应的燃料通过气道排出，电子则通过外电路回到电池的阴极，形成完整回路。

并伴随着向外输出电功。

在闭合电路下，只要燃料和空气持续不断地供给，那么电池就能连续将化学能转化为电能。

固体氧化物燃料电池的基本工作原理如图5所示：图5 工作原理示意图2固体燃料电池多物理场模拟2.1控制方程2.1.1动量守恒方程Navier-Stokes方程可以用来模拟燃料气道和空气气道中的动量传输过程。

Navier-Stokes方程和连续性方程相结合可表示为：(ρu∙∇)u=−∇p+∇∙[μ(∇u+(∇u)T)−23μ(∇∙u)I]∇∙(ρu)=0式中，ρ是密度，I是单位矩阵，μ是流体的粘滞系数。

2.1.2能量守恒方程温度是最重要的参数之一，因为SOFC很多性质与温度密切相关。

在流体区域热传导和热对流必须同时考虑，可以用方程表示为：∇∙N T(−k f∇T+f p Tu)=Q式中，f是流体的摩尔浓度，p是流体的摩尔热容，k f是流体的热导率。

2.1.3质量守恒方程对于SOFC 来说，电极中孔的直径典型的分布范围为0.2～1um 。

SOFC 多孔电极的质量运输模型最准确的是尘气模型。

尘气模型的摩尔形式如下：N iD iK+∑x j N i −x i N jD ijnj=1=−1(P∇x i +x i ∇p +x i ∇p kp D iKμ)式中，N i 是物质i 的总摩尔流量，x i 是物质i 的摩尔分数，R 是气体常数，T 是绝对温度，p 是气体总压强，k 是渗透系数，μ是粘滞系数，D iK eff是物质i 的有效Knudsen 扩散系数，D ij eff 是物质i 与物质j 的二元扩散系数。

2.1.4导电方程电子电流密度和离子电流密度由电荷守恒决定，结合欧姆定律，电荷守恒方程可以表示为：∇∙i el =∇∙(−σel eff∇φel ){0ASL−S currentAFL S currentFL∇∙i io =∇∙(−σio eff∇φio ){S currentAFL电解质−S current FL式中，i el 和i io 分别是电子电流密度和离子电流密度，∇φel （φio ）为局域的电子（离子）电势，S current 为电流源，σel eff （σio eff）是有效电子（离子）电导率。

2.2物理模型典型的阴极支撑的SOFC 由五层组成，如图6所示。

a.阳极电流收集层（ACCL ）。

为了减小气体输运的阻力，此层具有相对比较大的Ni 和YSZ 颗粒和高的孔隙率。

b.阳极功能层（AFL ），为了减小活化极化，此层具有相对比较小的Ni 和YSZ 颗粒和低的孔隙率以增大单位体积中的三相线长度。

c.电解质层d.阴极功能层（CFL ）,此层具有相对比较小的LSM 和YSZ 颗粒和低的孔隙率，增大了单位体积中的三相线长度，减小了阴极的活化损失。

e.阴极支撑层（CSL ），此层具有相对比较大的LSM 颗粒和高的孔隙率。

图6 典型的阴极支撑的SOFC 物理模型考虑到对称性，我们的计算区域仅仅包括重复单元的一半，如上图中处于中间框中的区域，气道及连接体rib 可以用气体浓度和电势边界条件的形式替代。

2.3数学模型2.3.1气体输运控制方程气体在SOFC 多孔电极中的输运过程非常复杂，包括三种截然不同的输运机制：Knudsen 扩散、分子扩散和黏滞流。

因次为了准确描述气体在SOFC 多孔电极中的输运过程就必须使用尘气模型，这是因为尘气模型不仅考虑了分子间的碰撞，同时也考虑了分子与孔壁的碰撞。

对于两组分的体系尘气模型经过推导可以给出每种组分流量的解析表达式。

物质i （i=1,2）的总的摩尔流量可以表达为：N i =−D 1 eff D iKeffD 1 eff +x 1D K eff+x D 1Keff ∇c i −c i[D 1K eff D KeffRTc tot (D 1 eff +x 1D K eff +x D 1K eff )+kμ]∇p =−D i ∇c i −c i k ′∇p=N i diffusion +N i convection2.3.2导电控制方程电子电流密度和离子电流密度由电荷守恒决定，结合欧姆定律，电荷守恒方程可以表示为：∇∙i el =∇∙(−σel eff ∇φel ){0ASL−S currentAFL S currentFL∇∙i io =∇∙(−σio eff∇φio ){S currentAFL0电解质−S currentFL2.4边界条件正如上面所描述的那样，此模型考虑了气体输运过程、离子导电过程、电子导电过程和Butler-Volmer 方程。

为了求解这些耦合的偏微分方程需要设置合理的边界条件。

表2列出了质量输运方程的边界设置。

由于在电极/电解质交界面发生电化学反应，所以在此边界上设置了物质的摩尔流量。

表3列出了电子导电方程和离子导电方程的边界设置。

表4列出了模型的基本参数。