1、第 12 卷 第 9 期2023 年 9 月Vol.12 No.9Sept.2023储能科学与技术Energy Storage Science and Technology固体氧化物燃料电池气流和热场的宏观尺度数值模拟研究进展栾凯夫1,2,3，蔡长焜1,2,3，谢满意1,2,3，张纯1,2,3，郑坤灿4，安胜利1,2,3（1内蒙古科技大学材料与冶金学院；2内蒙古先进陶瓷材料与器件重点实验室；3轻稀土资源绿色提取与高效利用教育部重点实验室；4内蒙古科技大学能源与环境学院，内蒙古 包头 014010）摘要：宏观数值仿真方法通过耦合多个物理场可对固体氧化物燃料电池(solid oxide fuel
2、cell，SOFC)进行仿真，在研究SOFC内部机理、外部输出性能等方面存在优势，可为电池优化设计提供依据。SOFC内部气体流动的均匀性直接影响电池效率，热场的分布会影响电池的发电性能和长期稳定性。本文总结了SOFC内部流道、外部歧管对流场的影响、优化传统流道以及设计新型流道中的宏观数值仿真模拟，以及宏观尺度数值仿真在热量传输和热稳定性方面的研究进展；概括了宏观数值模拟在分析燃料效率，结合多尺度模型，以及设计SOFC组件及新结构方面的研究进展。对应用宏观数值仿真方法研究SOFC进行总结和展望，认为目前需要统一SOFC结构设计的评价标准，以便量化对比。关键词：固体氧化物燃料电池；数值模拟；流场设
3、计；热模拟；性能优化doi：10.19799/ki.2095-4239.2023.0303 中图分类号：TM 911.4 文献标志码：A 文章编号：2095-4239（2023）09-2985-18Research progress of macroscale numerical simulation of fluid and thermal fields of solid oxide fuel cellsLUAN Kaifu1,2,3,CAI Changkun1,2,3,XIE Manyi1,2,3,ZHANG Chun1,2,3,ZHENG Kuncan4,AN Shengli1,2,3(
4、1School of Materials and Metallurgy,Inner Mongolia University of Science and Technology;2Inner Mongolia Key Laboratory of Advanced Ceramic Materials and Devices;3Key Laboratory of Green Extraction&Efficient Utilization of Rare-Earth Resources,Ministry of Education;4Energy and Environment School,Inne
5、r Mongolia University of Science and Technology,Baotou 014010,Inner Mongolia,China)Abstract:The macroscopic-scale numerical simulation method can simulate solid oxide fuel cells(SOFCs)by coupling multiple physical fields.It has advantages in studying the internal mechanism and external output perfor
6、mance of SOFCs,providing a basis for the optimization design of cells.The gas flow uniformity inside SOFCs directly affects the batterys efficiency,and the thermal field distribution will affect the power generation performance and long-term stability of the battery.This paper summarizes the researc
7、h progress of the influence of the internal flow channel and external manifold of SOFCs on the flow field,optimization of the traditional flow channel,and design of a new flow channel,and macroscopic-scale numerical 储能测试与评价收稿日期：2023-05-04；修改稿日期：2023-06-23。基金项目：国家自然科学基金项目（51974167）。第一作者：栾凯夫（1997），男，硕
8、士研究生，研究方向为固体氧化物燃料电池，E-mail：；通讯作者：安胜利，教授，研究方向为固体氧化物燃料电池，E-mail：shengli_。引用本文：栾凯夫,蔡长焜,谢满意,等.固体氧化物燃料电池气流和热场的宏观尺度数值模拟研究进展J.储能科学与技术,2023,12(9):2985-3002.Citation：LUAN Kaifu,CAI Changkun,XIE Manyi,et al.Research progress of macroscale numerical simulation of fluid and thermal fields of solid oxide fuel ce
9、llsJ.Energy Storage Science and Technology,2023,12(9):2985-3002.2023 年第 12 卷储能科学与技术simulation in heat transport and thermal stability.Furthermore,the research progress of macroscopic-scale numerical simulation in analyzing fuel efficiency,coupled multiscale model,and designing SOFC components and ne
10、w structures are reviewed.The application of macroscopic numerical simulation methods in studying SOFCs is summarized and prospected.It is considered necessary to unify the evaluation criteria of SOFC structure design for quantitative comparison.Keywords:solid oxide fuel cell;numerical simulation;fl
11、ow field design;thermal simulation;performance optimization随着化石燃料大量燃烧造成的环境污染、温室效应和能源短缺等问题的不断加剧，开发清洁可再生能源是目前的研究热点。固体氧化物燃料电池(solid oxide fuel cell，SOFC)作为一种发电器件可将燃料中的化学能直接转化为电能1-3。它具有高效、零污染、低成本和广泛的燃料来源等优点4-6，这使SOFC在动力辅助装置、可再生能源系统、分布式发电等领域的应用前景十分广阔。SOFC的工作温度一般在8001000，高温环境有利于提高SOFC的发电效率及燃料兼容性，但SOFC内部气体
12、流速、气体浓度、电流及温度的分布难以通过实验手段测量7-10。此外，SOFC涉及电化学反应、热量传递、质量动量传递、应力等多种因素，是一个内部工作过程复杂、多物理场耦合的非线性系统。考虑到实验手段研究内部物理场分布的复杂性11-13，为方便进行SOFC结构设计和系统控制，采用数值仿真的手段预测内部物理场分布可以节约大量实验成本14-15。SOFC模型建立和物理现象的研究可以基于不同层次的长度尺度，如微观尺度(10 m)16。微观尺度模型用于研究原子级或分子级的热力学和动力学反应机理，一般采用密度泛函理论、分子动力学等研究方法17。中尺度模型研究尺度介于微观与宏观之间，如研究多孔电极孔隙内部的反
13、应和传输过程18。对于宏观尺度模型来说，它与全局流场相匹配。宏观模型的经验参数可以基于微观或中尺度模型的结果，耦合多物理场19。目前，大多采用计算流体力学(computational fluid dynamics，CFD)软件 COMSOL 及 ANSYS Fluent 对宏观尺度 SOFC 模型的数值进行仿真分析。可实现研究SOFC的内部气体、热量、电流密度等分布，分析优化内部流道和外部歧管设计，探究不同操作参数对其性能产生的影响，以及全局系统热管理等20-21。此外，随着数值模拟与计算机计算技术的发展，更为复杂的SOFC数值模型的数值仿真结果与实验结果更加相符。这对于深入研究SOFC的工作
14、机理至关重要22。因此，利用数值模拟准确揭示SOFC内部热分布和不同流场对SOFC性能的影响、辅助加速研究进程以提高研究效率和效果已成为目前研究的热点。本文综述了近年来宏观尺度数值模拟在SOFC模型中的研究，包括数值仿真对SOFC中气流分布的模拟及其在SOFC流场设计中的应用，热模拟在SOFC中的应用及研究，以及宏观数值仿真分析SOFC其他方面的应用研究。本文总结了宏观尺度数值仿真研究SOFC领域的最新研究进展，并对未来通过数值仿真SOFC的应用设计提出展望。1 SOFC宏观尺度数值模拟SOFC 宏观数值模拟一般包括二维和三维模型。二维模型能够模拟电池内部的传质传热现象并简化计算量，但无法反映
15、准确的扩散情况23。而三维模型能够模拟SOFC内部的传质传热以及电化学反应，还可以对SOFC运行时涉及的物理和化学现象进行详细描述，但会有更大的计算量24。采用三维数值仿真模拟SOFC的步骤一般为建立几何模型，网格剖分，模型假设，建立数值模型，设定边界条件，求解计算以及输出结果。几何模型是根据实际计算需要而建立的，使用软件对建立好的几何模型进行网格剖分，其目的是将几何体分成有限的单元，便于对单元进行一阶或二阶计算。网格的数量、大小与质量将直接影响计算成本与精度，是模拟仿真中十分重要的一步。所以，一般会对网格进行调试，以达到成本与精度的平衡25。此外，还需要进行网格无关性检验来验证网格的有效性。
16、图1列举了几种不同的SOFC几何模型。2986第 9 期栾凯夫等：固体氧化物燃料电池气流和热场的宏观尺度数值模拟研究进展模型假设需根据具体情况确定，一般包括设置为稳态或瞬态，视气体流速和压力情况而设定为层流或湍流，气体混合物均视为理想气体，热边界条件假设，是否忽略SOFC内的热辐射等。对于边界条件的设定，通过求解软件以及添加用户自定义函数(user define function，UDF)，可以准确地设定边界条件。目前最为常用的CFD商用软件COMSOL和ANSYS Fluent分别采用有限元法(finite element method，FEM)和有限体积法(finite volume me
17、thod,FVM)在宏观尺度上对SOFC模型的边界图1SOFC几何模型(a)乙醇燃料SOFC二维模型29；(b)管式SOFC二维模型30；(c)交错式单通道模型31；(d)管式SOFC单元三维模型32；(e)平板式SOFC单元三维模型33；(f)双面阴极SOFC三维模型34；(g)平板式SOFC堆三维模型35；(h)管式SOFC堆三维模型36Fig.1Geometric model of SOFC:(a)2D model of ethanol SOFC29;(b)2D model of tubular SOFC30;(c)Staggered single channel model31;(d)
18、3D model of tubular SOFC cell32;(e)3D model of flat panel SOFC cell33;(f)3D model of SOFC with double-sided cathodes34;(g)3D model of panel SOFC stack35;(h)3D model of tubular SOFC stack3629872023 年第 12 卷储能科学与技术条件进行求解16。同时，SOFC数值模型的建立将直接影响计算量以及计算的准确性26。目前，研究SOFC宏观数值模型有两个方向：在保证准确性的情况下，简化大型SOFC堆的模型；随着
19、计算机算力的发展，建立更为复杂的精确模型。如Tseronis等人27建立了一个3D的非等温平板SOFC模型，首次在COMSOL多物理模型框架中描述多孔电极中的质量输运。该模型中引入了Dusty-gas模型(DGM)，相比基于Fick定律的Stefane-Maxwell扩散(SMG)建立的质量传输模型，在高电流密度的情况下更准确地预测了SOFC多孔电极内部的组分分布情况。ANSYS Fluent 在 2023R1 版本的 SOFC 模块中引入了Knudsen扩散系数的设定，以提高计算的准确性。多物理场耦合的数值模型多用于模拟SOFC单电池或涉及热的SOFC堆，在简化模型方面，He等人28通过统计
20、各项参数对电池性能的综合影响引入Monte Carlo模拟，将统计的电化学参数对电池的影响耦合到计算模型中，极大地缩小了计算量。在多物理场耦合的情况下，可将计算复杂的3D模型简化为多个 2D 模型。He 等人23验证了平板式SOFC的2D模型与3D非等温模型的空间平滑性和渐进简约性，简化后的模型在各方面与3D模型非常吻合，还可应用于各种流道形状细长且平行的SOFC，但计算成本降低了3个数量级以上。多物理场的耦合使得SOFC的数值模型较为复杂，若只研究电堆单元之间的气流均匀性，如电堆外部歧管的设计，可以将单元内部的物理量简化为边界条件进行计算，将极大地减小计算成本。2 气流模拟由于SOFC的工作
21、温度较高，内部的组分均为气相，因此SOFC的气流模拟均为单相流模拟。以离子导体电解质的SOFC为例：阴极侧通空气，以O2为氧化剂；阳极侧通燃料，并产生水蒸气、二氧化碳和一氧化碳气体。气体在多孔电极内部的扩散情况直接影响电化学反应的效率，进而影响SOFC的输出性能。另外，外部歧管可将气体输入到电堆各层，歧管的流场设计决定了SOFC堆内分配到各单元的气体是否均匀37。因此，采用CFD理论模拟SOFC内部多孔电极及气流通道内的组分分布情况，以及外部歧管分配气体的均匀情况，是目前对SOFC进行宏观数值模拟的研究重点38。2.1SOFC多孔电极内部扩散SOFC的电极均为多孔结构，其中气体的传质以扩散为主
22、，并受到在电极与电解质交界的三相界面(three phase boundary，TPB)处氧化还原反应的影响39。在电流密度较低的情况下，基于菲克模型(Fick model,FM)建立的 Stefan-Maxwell 模型(SMM)可以描述电极中的气体扩散：(-wiDij，effxj+(xj-wj)ppu-DTiTT)+(u )wi=Qi(1)式中，w为质量分数；x为摩尔分数；DTi为热扩散系数；Qi为化学反应引起的源项；Dij,eff为多孔介质中的有效扩散系数。Dij，eff=Dij(2)式中，Dij为Stefan-Maxwell二元扩散系数；为孔隙率；为弯曲度。而在电流密度较高的浓差极化阶
23、段，分子的平均自由程较小，多孔电极的孔径不变，使得孔径与分子扩散的自由程差距增大。此时，Knudsen扩散可以通过融合动力学的Knudsen系数Dk描述扩散过程：Dk=ur6(3)式中，u为气体分子的速度；r为平均孔径半径。若孔隙近似为圆形，分子i的扩散系数为Dik=97.0r TMi(4)式中，Dik为分子i的Knudsen扩散系数；T为温度；Mi为分子量。为了计算平均孔径半径(r)，使用多孔固体的表面积和孔隙率：r=2SAs(5)式中，SA为多孔电极的表面积比；s为固体颗粒的体密度。通过计算有效Knudsen扩散系数Dik,eff，可以解释多孔介质内的分子运动：Dik，eff=Dik(6)
24、当分子扩散和Knudsen扩散同时作用，多孔介质中的有效扩散系数Dij,eff可以计算为Dij，eff=(DijDij，kDij+Dij，k)(7)这种引入了Knudsen扩散的SMM模型，就是所2988第 9 期栾凯夫等：固体氧化物燃料电池气流和热场的宏观尺度数值模拟研究进展谓的Dusty-gas模型(DGM)。根据数学模型，可以计算出SOFC多孔电极的参数对内部扩散的影响40。如图 2 所示，目前 SOFC 阳极-电解质-阴极(positive electrode-electrolyte-negative electrode，PEN)的组装方式分为两种：增加阳极功能层或是阳极扩散层厚度的组
25、装方式，称为阳极支撑型SOFC(anode-supported SOFC，ASC)；而两极均有较厚功能层或扩散层的组装方式，称为双极支撑SOFC(anode-cathode-supported SOFC，ACSC)。由于通常研究的 SOFC 多为氧离子导体电解质，阳极侧不仅要通入燃料气体，还会生成气相产物，内部的分子活动和反应相比阴极侧更复杂，所以大多SOFC的组装方式会选择阳极支撑，以确保燃料和生成气可以充分扩散41。为探究不同的支撑形式对SOFC性能的影响，Su等人42分别建立了阳极支撑和双极支撑的 SOFC 平板电池的模型，比较了不同条件下两种支撑架构的性能。结果表明，虽然电池的性能差异
26、会受到阴极电导率、孔隙率、弯曲系数和输出电势的影响，但由于ACSC在阴极侧的电势差更小如图3(a)、(b)，使其性能始终优于ASC，平均电流密度较ASC高出17%，如图3(c)。2.2内部流场气体组分在电池内部的分布情况直接影响电化学反应的效率，进而影响SOFC的输出特性。通过CFD数值模拟可以分析各组分在电池内部的分布情况。多孔电极内的动量守恒方程可以通过Darcy-Brinkmam方程描述：(u)=Qi(8)u u=-p+(u+u)-23u-u(9)式中，为气体密度；u为速度矢量；Qi为质量源项；为孔隙率；为渗透率；为黏度。由于SOFC内部通入的气体速度较慢，一般的模拟研究将通道气体流动假
27、设为层流，气体假设为不可压缩流体，其动量守恒的控制方程为u u=-p+(u+u)-23u(10)因为式(10)是当时的Darcy-Brinkmam方程，因此可以描述气体由通道区域进入多孔电极区域的运动。质量源项Qi由电化学反应的消耗或生成影响。电池内部的气体会通过浓度影响电化学反应的效率，也会传输热量。因此SOFC的内部流道对内部流场的影响是对其宏观模拟的研究重点。由于加工和组装简单，截面几何形状为矩形，平行排列的图2电极内的电荷和气体输送路径示意图42(a)ASC设计中的气体输送路径；(b)ASC设计中的电流导电路径；(c)ACSC中的气体输送路径；(d)ACSC中的电流导电路径Fig.2S
28、chematic diagram of the charge and gas transport paths within the electrodes42:(a)Gas transport paths in typical support designs;(b)Electric current conducting paths in typical support designs;(c)Gas transport paths in ACSC;(d)Electric current conducting paths in ACSC29892023 年第 12 卷储能科学与技术通道是目前平板型S
29、OFC的主流如图4(a)。Hanieh等人43比较了不同截面形状的通道如图 4(b)对SOFC性能的影响，结果表明，在截面底边相等的情况下，矩形流道的电流密度最高，温差更小。平行或类平行通道的平板式SOFC两极通道的排列方式可以分为共流、逆流和交叉流(图5)。Zhang等人44通过模拟比较了阳极支撑SOFC的3种排列方式，表明共流和逆流配置的电流密度明显大于交叉流配置，逆流配置的通道间气体均匀性更好(图6)。由于使用实验手段表征电池内部气体分布情况十分困难，并且燃料在SOFC内部分布的均匀性会影响电池工作的效率，因此，Lin等人45设定了一个有别于均匀系数Fu的非均匀系数CV，用以量化通道中心
30、速度之间的差异：Fu=umini=1nui(11)CV=STD1ni=1nui(12)式中，n是通道数；ui是i通道中心的气体流速；STD是ui的标准差。Lin等人45模拟了通入60%H2+40%CO混合燃料的平板式SOFC，结果表明CV=1.1010比CV=0.1570的H2和CO燃料利用率分别低2.8%和4.5%，平均电流密度降低5.6%。因此，改善电池内部的气体均匀性是提高电池效率的关键。为改善SOFC内部气体的均匀性，电极两侧气流通道的设计至关重要。Duhn等人46通过建立CFD模型设计了在单电池内部气体分配到平行通道的气体分流器(图7)，将流动均匀系数(Fu)提升到0.978，并使燃
31、料利用率提升了72.9%。对SOFC的内部流场模拟的主要目的是分析气流在SOFC内部的分布情况：单通道模型主要优化通道方向的浓度变化；多通道的单电池模型主要优化各通道间气流分配的均匀性。图3阴极内电势分布与输出电流密度42(a)ASC阴极内电势分布；(b)ACSC阴极内电势分布；(c)输出电压对平均电流密度的影响Fig.3Distribution of electric potential in cathode and the output average current density42:(a)Distribution of electric potential in cathode of
32、 the ASC;(b)Distribution of electric potential in cathode of the ACSC;(c)The effect of the output voltage on the average current density图4平板式SOFC气流通道43(a)阳极支撑结构；(b)不同截面形状的气流通道Fig.4Flow channels of planar-type SOFC43:(a)Schematics of an anode-supported;(b)Schematic of different flow channel cross sec
33、tion shape2990第 9 期栾凯夫等：固体氧化物燃料电池气流和热场的宏观尺度数值模拟研究进展图7优化后的SOFC气体分流器46Fig.7Optimization of a new flow design for solid oxide cells46图5平板式SOFC40(a)顺流和逆流；(b)交叉流Fig.5Planar-type SOFC44:(a)Co-and counter-flow;(b)Cross-flow configuration(not in scale)图6氢气摩尔分数预测44(a)共流；(b)逆流Fig.6Mole fraction of H2 predicte
34、d for44:(a)Co-flow;(b)Counter-flow29912023 年第 12 卷储能科学与技术2.3外部歧管的设计与优化外部歧管是SOFC堆的重要组成部分之一，它可将燃料气与氧气分别通入各单电池中。外部歧管的设计会影响SOFC的发电性能。SOFC堆由许多单电池单元堆叠在一起组成，燃料气和空气由入口进入，需要通过外部歧管运输到各个单元。电池堆内通道的分布方式和几何结构也是影响气流均匀性的关键，根据入口和出口的方向不同，电池堆的设计可分为U形流道(进出口流量方向相反)和Z形(进出口流量方向相同)流道两种(图8)。Zou等人47开发了30层平板型SOFC电池堆的两种结构模型，比较
35、了两种结构模型的流量分布。结果表明，U形结构每层的流量供应从下到上逐渐减少，气流主要在下层积累；而Z形结构每层的流量供应从下到上逐渐增加，流体主要分布在上层。在单层流场中，U形结构的气体分布更均匀。Su等人48探究了单元数量与通道分布对均匀性的影响，通过比较传统矩形肋结构、离散圆柱肋结构、离散交错肋结构的模拟结果，发现10层平板型SOFC堆受到肋结构的影响最小，并且矩形肋结构是燃料通道的最佳选择。外部歧管的设计直接影响单元间的气体分配，通过模拟可以便利地研究新设计的外部歧管对SOFC堆的影响。Rashid等人49建立了由80个单元组成的平板式SOFC堆的3D模型图9(a)，模拟了未优化歧管与优
36、化歧管(加装了分流器)的气体分布情况，通过比较常用歧管的不同尺寸参数对均匀性的影响、优化与未优化的性能表现、优化后分流器尺寸参数的影响，发现优化后的歧管使得平管SOFC堆性能提升了8%。Zhao等人50通过CFD模拟研究了40个单元的SOFC堆的流量分布和压力变化，经过对外部歧管的设计优化和参数比较后，总结出气流的均匀性强烈地依赖于流场的几何形状，包括管与流形之间的连接位置、流形的尺寸和管的数量。如图9(b)所示，通过优化流形的几何结构，可以降低用于描述气流分布均匀性的流速比。Dong等人51设计并优化一种树形模块化歧管，从气流分布均匀性和总压损失两方面比较了两种不同的通道尺寸，研究了通道长度
37、和空气流量对均匀性和压降的影响。结果表明树形模块化歧管设计在空气流量较大时，单元内的空气流动均匀度均可达到0.99左右。综上所述，在对SOFC的气流模拟研究中，主图8电池堆空气侧的设计形式47(a)U形流道示意图与氧气摩尔分数分布；(b)Z形流道示意图与氧气摩尔分数分布Fig.8Air-side47:(a)U-type configuration of air flow path on stack level;(b)Z-type configuration of air flow path on stack level2992第 9 期栾凯夫等：固体氧化物燃料电池气流和热场的宏观尺度数值模拟研
38、究进展要关注的是气体流动均匀性。单电池关注内部流场的均匀性，包括电极内的扩散和通道内的流动均匀性。单电池的模拟需要耦合传热和能量守恒方程，才能准确模拟内部扩散。单电池内部的均匀性主要受到两方面的影响，一是分配到各通道的气体是否均匀，二是通道内延通道方向的气体浓度变化。电堆则主要关注单元间的均匀性，通常认为气体通入量最小的单元决定了电堆的整体效率。受此影响，SOFC的设计不能一味地扩大活性面积或单电池数量，应存在一个发电功率和发电效率的最优区间。此外，由于操作参数、燃料、材料参数、电池活性面积等因素的不同，导致设计出的流场虽然在某些方面有突出的性能，但没有衡量SOFC设计优劣的标准，难以进行量化
39、对比。图9歧管优化前后歧管内流速(a)80层平板式SOFC堆49；(b)外部歧管的形状优化50Fig.9Manifold optimization before and after manifold flow rate:(a)80-layer flat SOFC stack49;(b)Shape optimization of external manifold5029932023 年第 12 卷储能科学与技术3 热模拟SOFC内部的热量运输包括多孔电极内固相表面与气相间的传热、固体和多孔结构中的热传递以及内部气流的对流换热52。SOFC的热源包括在电解质与电极的三相界面处的电化学反应热，一氧
40、化碳在电极内部重整反应产生的反应热、甲烷在多孔阳极中的蒸汽重整反应中吸收的热量、电流流动产生的焦耳热。燃料在SOFC各处的浓度影响局部的反应速率，导致温度分布不均匀53。SOFC内部的热分布直接影响长期稳定性，但难以通过实验手段测量SOFC内部的热分布。因此通过宏观数值模拟的方式对SOFC内的热分布进行预测，对于评估和优化电池的整体性能以及避免热机械失效至关重要。3.1热分布SOFC内的反应热在TPB附近产生，并通过以下方式传输：多孔电极固相中的传导，通过孔内的对流从固相传递到气相，气体通过微孔平流到流动通道。由于气相与固相之间的导热系数有很大的差异，采用局部温度非平衡(LTNE)方法来预测多
41、孔电极内固相和气相之间的传导。电池内固相的热量分布为(-ks Ts)=Qs(13)式中，ks为固体的导热系数；Ts为固相温度，Qs为热源。多孔电极内气相的热量分布为(-kg Tg)=Qg-gcp，gu Tg(14)式中，kg为气体的导热系数；cp,g为气相热容；Tg为气相温度；Qg为固相传输到气相的热及气相的产热。当电极内部不考虑辐射时，气固两相的传热系数hs,g,por为hs，g，por=2kgdp(15)式中，dp为多孔电极的颗粒直径。气固两相之间的传热取决于温差和颗粒表面积：Qg=SAhs，g，por(Tg-Ts)(16)SOFC的热源来自电化学反应和焦耳热，Baek等人54建立了多尺度
42、多物理场耦合的数值模型，此模型耦合了电极微结构模型、电荷传输和电化学反应的连续电极模型、电极和气体通道中气体扩散和流动的传质模型、反应热的传导和对流传输(包括焦耳热)的传热模型，并模拟了中温SOFC的热分布。结果表明，共流配置的中温 SOFC，热分布更为均匀。Lee等人55建立耦合局部产热的传热特性的模型，探究产热机理及其与内部热条件的关系，并考虑了金属连接体和电流收集材料的传热，详细阐明了对传热的影响：阳极支撑型SOFC的热量主要来自电解质和电极界面附近的活性电极层，即TPB中的氧化反应；沿电池厚度方向的主要传热途径为阳极内的气体对流，阴极内的主要传热途径为固态传导；互联件的高热阻会增加内部
43、温度的增量及其梯度；改变互联件的热阻会影响进入气体和重复单元之间的热传递。在对SOFC的热模拟研究中，一般只考虑对流和热传递，而热辐射影响的模拟研究较少。Zeng等人56建立了辐射传热的三维综合模型，重点分析了辐射传热过程及其影响。采用辐射度法对空气和燃料通道中的热辐射进行了分析，表明发射率对模拟结果的影响很小。对于不同的流动布置，共流辐射对SOFC性能的影响大于逆流的影响。因此，忽略辐射的较小影响有利于计算的简化。Lin等人57模拟了逆流的平板型SOFC堆内的热量分布(图10)，气流分布控制着电池整体的电流密度分布，热分布主要控制着SOFC堆靠近顶部和底部单元的电流密度。Kim等人58对一种
44、功率为1 kW，并由30个单元组成的SOFC堆进行了数值模拟。当空气和燃料传输到堆顶部时，气流对热分布的影响更加显著，单元内的温差较小，图1020个单元组成的SOFC堆栈内的温度分布57Fig.10Temperature distributions within a 20-cell SOFC stack572994第 9 期栾凯夫等：固体氧化物燃料电池气流和热场的宏观尺度数值模拟研究进展温度升高得较慢，在进气口附近的单元和连接体上的温度梯度较大。有别于稳态条件下的研究，SOFC的加热过程需要建立瞬态模型来模拟。Zheng等人59通过建立瞬态模型，模拟了单通道情况下互联体的不同肋与通道宽度比对热
45、分布的影响，并研究了比例的连接体在SOFC升温过程中对热场的影响。结果表明，由于连接体理论导热量最大，其体积对升温速率和升温时间起主导作用。垂直于通道方向的最大温度梯度则随着通道宽度和肋宽的增加而增大图11(a)。而由于连接体体积的减小，平行于通道方向图11(b)和垂直于电解质层方向图11(c)的最大温度梯度都随着肋宽的增加而减小。Al-Masri等60模拟了用于辅助动力单元的小型SOFC堆的瞬态加热过程和温度场(图12)，加热过程中的最高温度和温度梯度出现在入口区域附近，而在加热结束时，由于热气体流量的增加，出口区域也会出现更高的温度梯度。3.2热稳定性热稳定性一直是对SOFC研究的主要方向
46、。为探究导致热应力的因素，Guo等人33利用多物理耦合模型预测了SOFC的热分布，并对SOFC堆进行了瞬态的力学性能分析，分析得出导致热应力的主要因素是温度变化导致的热应力和热膨胀系数不匹配导致的蠕变。Xu等人61对阳极支撑平板型SOFC建立多物理耦合模型，分析了共流和逆流配置的通道内应力的分布(图13)，并分析得出，与逆流工况相比，共流工况的热应力略低，且电池的固定约束加剧了热应力行为。与平板式 SOFC 相比，管式 SOFC 由于结构优势，工作时产生的热应力较小，因此也是SOFC热模拟研究的重点。Gari等人32通过宏观数值模拟，研究了在电极内加入功能梯度材料(functionally g
47、raded materials)对管状 SOFC 热应力 分 布 的 影 响。结 果 表 明，材 料 分 级 指 数(material grading index)为2时，热膨胀系数不匹配导致的应力降低最明显。Hao 等人30在管式SOFC的燃料通道中放置管式分离和法兰分离器来控制局部热的产生(图14)，从而影响燃料浓度分布和多孔阳极中相关的热电化学反应。两种分离器均能有效降低槽内温度梯度，但法兰分离器的温度分布更均匀。长度为25 mm的法兰分离器可使电池在仅损失17%功率密度的情况下，最高温度梯度由50/cm降至18.7/cm。影响SOFC堆热稳定性的因素更复杂，且热应力在更大尺寸的 SOF
48、C 堆中的作用效果更加明显。Fang等人9的模拟结果发现，在活性面积为144 mm2的平板型电池堆中，非等温条件的应变比等温条件大40%。Pianko-Oprych等人62则是对微管式SOFC堆进行3D建模(图15)，模拟了不同微观排布下的热应力。结果表明，采用外部气流冷却设计时，应变最小，并由于均匀的温度分布，轴向应力和总应力最小，防止了SOFC堆叠组件中过大的热应力。Promsen等人63设计了一种带有水冷的SOFC堆，并为其建立模型进行模拟，将水冷堆与传统气冷堆进行对比，发现水冷堆显著改善了堆内的温度分布，使电化学性能得到改善。热量从多个方面对SOFC的性能造成影响。目前，宏观数值模拟的
49、稳态模型可用于分析平板、管式SOFC单电池与堆在内部的热分布、热应力分图11截面上单通道SOFC的热分布59(a)x=0；(b)y=0；(c)z=90Fig.11Heat distribution of single channel SOFC on cross section59:(a)x=0;(b)y=0;(c)z=9029952023 年第 12 卷储能科学与技术布；瞬态模型可用于分析SOFC的升温过程。部分研究的结果在模拟中结合了微观尺度模型，可以使得计算更准确。SOFC是一个发电并产热的热电联供装置，利用产出的热量，会使得电池的能量效率进一步提升，但鲜有研究讨论SOFC产生热量的导出和
50、应用。设计SOFC适用的热量管理装置，会成为提高燃料利用率的发展方向。4 其他因素对SOFC性能的影响除了对流场和热稳定性的数值模拟，有许多研究指向了其他可以影响SOFC性能的因素。通过宏观数值模拟，SOFC几何设计和各项参数设置在数值模型的适用范围内可以自由变动，这使得SOFC的结构设计十分便利。图12加热过程的温度和温度梯度分布60(a)2010s；(b)4180s；(c)6200s；(s)8015s60Fig.12Temperature field and temperature gradient in plane P1 at t=2010s,4180s,6200s,8015s60299