猫咪的室友SOFC仿真做了几年，我的一个深切体会是：搞懂电化学机理不算难，真正让人头疼的是传质传热。就拿五层电堆模型来说，很多时候电池性能上不去、温度分布不均、局部热应力过大，根源都在气体输运和热量管理上。这篇内容就围绕我搭建SOFC五层电堆模型的过程展开，重点讲透里面的传质传热问题，包括模型怎么建、参数怎么定、结果怎么判读、坑怎么避。无论你是刚入门的燃料电池方向学生，还是正在做电堆设计的工程师，或者想用多物理场仿真工具分析电堆性能的研究人员，这篇都值得收藏。
1. 为什么盯着五层电堆模型不放：结构拆解与建模思路
1.1 五层结构到底是谁：每一层在传质传热里的角色
固体氧化物燃料电池的电堆，在微观上是一堆陶瓷和金属的叠合，但宏观上通常抽象成“五层”来处理。这五层依次是：阳极侧连接体、阳极、电解质、阴极、阴极侧连接体。有些文献会把密封件和流道也算进去，但核心传质传热分析基本围绕这五层展开。
每一层的职责完全不同，在传质传热中的角色也不一样。连接体（interconnect）在平板式SOFC里一般用铁素体不锈钢，它的主要任务是导电、导气、隔离相邻电池的燃气和空气，同时又承担着把电池产生的热量导出去的功能。它的热导率比较高，是电堆里重要的散热通道，但同时也是热应力最集中的地方，因为金属和陶瓷的热膨胀系数差异很大。
阳极为Ni-YSZ金属陶瓷，孔隙率一般30%-40%，燃料气体在这里扩散并发生氢氧化反应。阴极为LSM或LSCF材料，也是多孔结构，氧气在这里被还原成氧离子。电解质是致密的YSZ，只传导氧离子，不传导电子，电子必须走外电路做功。传质在电解质这里基本不发生，但传热是实打实的——电解质的热导率最低，是整个电堆里热阻最大的“瓶颈层”。
所以五层结构不只是“电化学活性区的夹心”，它实际上是一个完整的热输运和质量输运系统。每一层的厚度、孔隙率、热导率、电导率都在影响整体的传质传热行为。做模型时如果只把注意力放在阳极和阴极的反应动力学上，忽略连接体与流道的换热效应，结果必然失真。
1.2 建模选型的底层逻辑：为什么不做单电池而是做电堆模型
做SOFC仿真，有三类典型模型：单元级模型（单电池）、电堆级模型、系统级模型（与重整器、燃烧器、热交换器耦合）。很多最初接触SOFC仿真的人会问：单电池模型不是更简单吗？先把单元做准，再扩展到电堆不就行了？
这个思路在逻辑上没问题，但实际工程决策中，直接建模五层电堆模型的价值比单电池模型大得多，原因有三个。
第一，单电池模型通常假设边界条件均匀，入口流量、温度、浓度都是理想分布。但真实电堆里，气体通过流道分配后，沿着流道长度方向必然存在浓度梯度和温度梯度。电堆级模型才能把“气体分配不均”这种工程核心问题暴露出来。
第二，电堆里的热行为是三维的。单电池模型往往是二维甚至一维，把热流简化成沿着厚度方向传递，无法反映连接体肋片处的横向热传导及边缘散热。五层电堆模型可以在中间平面上捕捉“热点”的出现位置和扩展范围，这对密封设计和热应力评估都至关重要。
第三，电堆模型直接面向工程指标。电堆的输出电压、功率、燃料利用率、温度梯度上限，都是电堆级的概念。客户和设计部门关心的不是某一层电极的过电位，而是整堆的性能和寿命表现。五层电堆模型给出的结果，能直接对接电堆测试数据，起到标定和预测作用。
1.3 温度与浓度：两个贯穿全局的耦合变量
做SOFC仿真，最怕的就是“只盯着一个物理场”。电化学、传热、传质、流场、电场，这五个场在五层结构里高度耦合，其中温度和浓度是两个贯穿全局的耦合变量。
温度影响什么？影响电解质的离子电导率（Arrhenius关系）、电极的反应速率常数、气体的扩散系数、组分的热物性。温度高了，离子电导率指数上升，反应速率加快；但温度过高又加速材料烧结和退化。温度梯度大了，热应力变大，陶瓷电极可能开裂，密封件可能失效。
浓度影响什么？影响局部Nernst电压、反应速率、浓差极化大小。阳极侧氢气浓度太低，会导致燃料饥饿，局部电压骤降；阴极侧氧气浓度太低，浓差极化会急剧增大，尤其是在大电流工况下。
而且温度和浓度之间本身就是耦合的：反应热会改变温度，温度会改变扩散系数和电导率，进而改变浓度分布和产热功率。这是一个正反馈系统。比如在高电流密度区域，局部产热增大，温度升高，电导率上升，电流密度进一步集中，于是那个区域的温度更高。仿真中如果不控制这种正反馈，温度会“飞掉”，直接发散。
所以建模的首要任务，不是把方程写得多么复杂，而是把温度和浓度这两个变量的耦合路径搞清楚，抓住主要的反馈环节。后面各章展开的实操细节，本质上都是在处理这两个变量的行为。
需要模型API调用？ 免费领10W Token，多模型网关一键接入 Claude、DeepSeek 等主流模型。
2. 传质部分：气体怎么在五层结构里“挤”过去
2.1 多孔电极里的扩散机制：分子扩散与克努森扩散
SOFC的气体输运分为两部分：流道内的对流传质和电极多孔层内的扩散传质。流道内的传质相对简单，由流速和流道截面决定。真正的难点在多孔电极内部——气体分子在孔隙里一边碰撞一边挤向前，按孔径和压力的不同，会呈现出两种典型的扩散模式。
当孔径远大于气体分子平均自由程时，分子之间的碰撞占主导，这是分子扩散（Fick扩散），用二元扩散系数描述。但当孔径与分子平均自由程接近甚至更小的时候（典型多孔电极孔径在0.5-2微米），分子与孔壁的碰撞变得不可忽略，这种模式叫克努森扩散。
很多初次接触的人会问：什么时候该用哪种模型？答案是看克努森数和孔径。工程上为了让计算不过于复杂，常用Bosanquet公式把二者串联起来，得到有效扩散系数：
1/D_eff = 1/D_molecular + 1/D_knudsen
这里的关键是孔隙率ε和曲折因子τ的修正：D_eff = D * ε / τ。
孔隙率好理解，就是孔隙体积占比；曲折因子是描述孔道弯弯曲曲程度的量，一般为3-6。我曾经用默认值1算过一次，结果严重高估了扩散通量，导致在高电流密度下预测的浓差极化明显偏小，跟实验对不上。后来把曲折因子调到4.5，曲线才吻合。
具体到组分还要注意：氢气/水蒸气在阳极的二元扩散系数，和氧气/氮气在阴极的二元扩散系数相差很大，不能混用。而且温度升高时扩散系数按T^1.75左右增大，所以高温工况下即使反应消耗更快，扩散通常也不是瓶颈；反而是850℃以下或大电流密度时，浓差极化会快速上升。
2.2 流道设计对传质均匀性的影响
五层电堆模型里的流道，通常是在连接体上开槽形成的。常见的流道形式有平行流道、蛇形流道、网状流道。流道的几何参数——槽宽、槽深、肋宽、流道长度——直接决定气体分配的均匀性。
我建完几何模型，第一步都先做纯流场仿真，不开电化学反应，只看速度场和压力场分布。这一步能快速判断流道设计是否合理。如果流道入口和出口的压降太大，或者某个流道流速明显偏低，后面的电化学仿真结果基本没有意义，因为局部燃料供应已经出问题了。
举个直观的例子。对于一个平行流道结构，如果流道入口的分配腔设计得过窄，中间流道流速高，两侧流道流速低，导致阴极侧两侧氧气供应不足。反映到电化学上，就是两侧区域输出电压偏低、局部电流密度偏小、温度偏低。这个问题如果不做电堆级流场分析，永远发现不了。
流道设计里还有个容易忽略的点：肋片处（流道之间的实体部分）的气体输运。流道内气体直接接触电极表面对流传质很强，但肋片下方电极内部的气体只能靠侧面扩散过去，扩散路径长，常常出现“肋下缺氧”现象。肋片越宽，这个现象越严重。所以做三维模型时，如果网格在肋片下方不加密，这一块扩散阻力会被严重低估。
2.3 浓差极化的工程代价
浓差极化是传质受限的直接体现。阳极侧，氢气从流道扩散到反应位点，同时水蒸气从反应位点扩散回流道。这个双向往返过程如果受阻，反应位点处的局部氢气分压就会低于流道内分压，根据Nernst方程，局部电压就会下降。阴极侧类似，氧气被消耗后，氮气富集，氧气分压下降。
浓差极化的工程代价，在额定工况下可能不算大，但在两个场景下会变得致命。
第一个场景是高温差大电流。电流密度越高，单位面积消耗的燃料和氧气越多，浓度梯度越大，浓差极化呈非线性增长。我计算过，当电流密度从0.3 A/cm²提升到0.8 A/cm²时，阴极浓差极化的增量可以达到活化极化的两倍。也就是说，在大电流下瓶颈从“反应慢”变成了“气不够用”。
第二个场景是燃料利用率太高。当燃料利用率超过80%时，阳极尾部氢气浓度已经很低，局部会出现“燃料饥饿”。这时该区域为了维持总电流，只能拉大局部过电位，导致局部温度异常升高，甚至造成镍氧化。这也是为什么SOFC电堆很少长时间跑在85%以上燃料利用率的原因之一。
实际操作中，判断浓差极化是否过大的一个快速方法是看极化曲线（电压-电流密度曲线）在高电流密度区是否出现急剧掉电压的“膝部”。如果膝部出现得过早，优先检查电极孔隙率、曲折因子和流道形式，而不是去调电化学参数。
3. 传热部分：热量怎么在电堆里“走”起来
3.1 热源解析：三类极化的发热机制
SOFC电堆里的热量来源经常被简单概括为“反应放热”，但实际上热源至少有三类，每一类的成因和量级都不一样。
第一类是可逆热，对应电化学反应的热力学焓变与电功之差。氢气氧化反应的焓变约241.8 kJ/mol（低热值），如果按可逆电压1.1V计算，只有约一半焓变转化为电能，其余作为可逆热释放。这部分热量即使在零电流下也存在，对应的是电堆在开路状态下的保温需求。
第二类是活化极化热。电化学反应要跨越活化能垒，这部分过电位对应的能量损耗以热的形式释放。活化极化的数值在高电流密度下可能达到0.1-0.3V，乘以电流密度就是不小的热流密度。
第三类是欧姆热，也叫焦耳热。电流通过电解质、电极、连接体、接触界面时，离子和电子输运都要克服电阻，产生I²R热。电解质是主要的欧姆热来源，因为YSZ的离子电导率虽然随温度升高急剧增加，但在常规工作温度下仍然是电堆里电阻最大的部分。
把这些热源加起来，乘以对应的电流分布，才是真正的局部热源项。做模型时最忌讳的就是把热源均匀分摊到整个电堆体上。局部热源分布不均，正是温度场出现局部峰值的根源。
3.2 热传递路径：传导、对流、辐射怎么处理
热源产生后，热量沿着三条路径散走：固体的导热、流道内气体的对流传热、高温下的辐射传热。
固态导热是最主要的路径。连接体金属的热导率大约是20-30 W/(m·K)，Ni-YSZ阳极大约在10 W/(m·K)左右，YSZ电解质约2 W/(m·K)。热流在厚度方向上遇到电解质这个大热阻，会产生明显的温度梯度。这就是为什么电堆仿真中，电解质层内的温度变化往往是最大的。
对流换热发生在流道内。空气一般做过量系数3-5，流量大、换热能力强；燃料侧流量小，换热能力弱。因此在温度控制上，空气流量是主要调节手段。很多SOFC系统设计里把空气过量系数调得很大，不完全是用氧需要，更多是为了把电堆里的热带走。
辐射在800℃以上就不能忽略了。SOFC工作温度一般在700-900℃范围，此时固体表面之间的辐射换热量可能占总体换热的10%-30%。建模时至少要用表面-表面辐射模型，把相邻层之间的辐射换热算进去。忽略辐射的结果是电堆中部的温度明显偏高。
我个人的处理习惯是：初期快速验证模型时先用有效导热系数等效辐射效应，把辐射项“藏”进去，等模型骨架跑通了，再把辐射模型打开，对比两次的温度场差异。如果差异小于5%，说明辐射不敏感；如果差异大，说明这个工况下辐射必须显式建模。
3.3 温度场对电化学性能的反馈
很多人算完温度场就停了，觉得传热计算只是个“后处理”。实际上温度场对电化学性能的反馈是必须迭代的，因为电流密度分布反过来又决定热源分布。
YSZ电解质的离子电导率满足Arrhenius关系：σ = σ0 * exp(-Ea / RT)。这个关系意味着温度每升高100℃左右，离子电导率可以提升一个数量级。所以高温区域的电解质电阻小，电流密度自然更高，而更高的电流密度又带来更多的欧姆热和活化热，进一步推高该区域的温度。
这个正反馈如果不能收敛，就会在仿真里表现为温度飞升（thermal runaway）。实际电堆里也会有类似现象，只是材料热容、对流散热和系统控制会在一定程度上限制温度上升幅度。如果仿真里温度场不收敛，多半是热源项、边界条件或网格分辨率出了问题。
另一方面，温度太低的地方，电解质电导率下降，欧姆极化增大，电流密度被“推走”，导致该区域产热更少，进一步变冷。电堆里的热点和冷点就是这么形成的。
所以在五层电堆模型里，电化学模块和传热模块必须双向耦合。要么用全耦合求解器，要么做固定点迭代：先猜温度场，算电化学，得到热源，再算温度场，更新物性参数，再算电化学，如此反复直到收敛。我实测下来，固定点迭代如果阻尼设得好，速度和稳定性都不错，而且更容易定位每个物理场的收敛问题。
4. 实操过程：五层电堆模型的搭建与关键参数设置
4.1 几何建模与网格划分要点
五层电堆模型的几何建模原则就八个字：抓住主要矛盾，能简则简。
一个完整电堆可能有几十片重复单元，逐片建模意义不大。通常做法是取一个重复单元（包括两层连接体、一个阳极、一个电解质、一个阴极，以及流道），在垂直于平面的方向施加周期性边界条件，代表整堆行为。这个假设在电堆中部的电池上是合理的，边缘片则需要额外考虑边缘散热。
几何尺寸上的建议是：先参考典型文献，不要一上来就照抄企业图纸。比如单电池有效面积100 cm²左右，阳极厚度300-500微米，电解质10-20微米，阴极50-100微米，连接体1-3毫米，流道深度1-2毫米。电解质厚度对性能影响极大，模拟时如果厚度设错一个数量级，极化曲线会明显偏离实验。
网格划分是所有步骤里最影响效率和精度的一步。我的经验是：流道和电极的界面边界层必须加密，电解质层厚度方向至少5层网格，电极内部在靠近电解质一侧要细化，因为化学反应主要发生在三相界面附近。肋片下方的电极区域也要加密，否则扩散传质会被低估。全局网格数控制在几十万到几百万之间即可，配合相应算力，一次性算完。
网格划分前，先做一次“几何清理”：把圆角、倒角、螺纹孔等对传质传热影响极小的特征删掉。这些特征只会增加网格数、拖慢求解，却对结果没有任何实质性贡献。
4.2 边界条件与物性参数选取
边界条件设错了，后面一切白做。SOFC五层电堆模型的边界条件一般包括：
入口边界：氢气侧入口设定组分（典型为97% H₂ + 3% H₂O，或更高湿含量用于防止阳极积碳），温度（比如1073K），质量流量或流速。空气侧入口设定21% O₂ + 79% N₂，温度同样1073K，流量按空气过量系数3-5算。
出口边界：一般用压力出口，设定常压或略低于常压，避免回流。
壁面边界：流道壁面可设成无滑移，电堆外表面设对流换热系数，模拟电堆对环境的散热。这个散热系数的取值对边缘温度场影响很大，通常取5-20 W/(m²·K)之间。
电学边界：设定外加电压或电流密度，确定工作点。如果在恒流下跑，要指定总电流，让模型自动求出所需电压；如果在恒压跑，要指定电压，让模型求出电流分布。
物性参数方面，最容易被忽略的是“温度依赖性”。很多同学喜欢取常数，比如YSZ电导率直接给一个固定值。实际上电解质电导率在700-900℃范围内变化可达数倍，用常数会严重歪曲电流密度分布。同理，气体的比热容、热导率、粘度都随温度和组分变化，虽然在窄温度范围内变化不大，但电堆内部温差可能超过100℃、组分差异也大，常数假设会带来可观误差。
扩散系数的选值更是关键。前面提到Bosanquet公式，实际用的时候二元扩散系数要用Chapman-Enskog理论算，或者查表插件；克努森扩散系数则取决于孔径。孔径这个值不同文献差异很大，从亚微米到几微米都有，最好结合自己材料的压汞法孔径测试数据来定。
4.3 求解策略：稳态与瞬态的取舍
五层电堆模型的求解策略直接关系到迭代收敛和计算时间。我的习惯是分三步走。
第一步，先跑纯流场。关闭电化学反应和传热，只看速度场和压力场，检查流道内有没有回流、滞止区或明显分布不均。此步网格不需要太密，跑得快，主要用来验证几何和边界条件。
第二步，开等温电化学。固定温度（比如1073K），打开电化学反应和组分输运，求解浓度场和电流密度分布。此时没有热反馈，模型收敛较为容易，可以快速获得极化特性，并和文献里的极化曲线做对比，验证电化学参数是否标定正确。
第三步，全耦合。打开能量方程，让温度场与浓度场、电化学场双向耦合。这时有一个常用的技巧：用第二步的等温解作为初始值，可以大幅提高收敛速度。
求解器设置上，稳态求解时有必要给动量、能量、组分方程开启适当的阻尼因子。尤其是能量方程，温度场的惯性小，强非线性下容易振荡。我一般把能量方程的阻尼因子先调小，比如0.6左右，如果发现迭代曲线在振荡再进一步降低；等迭代稳定后再调回0.9-1.0跑最终结果。
瞬态仿真我只有在研究启动、停机、负载突变的温度响应时才用。稳态仿真的优点是可以快速扫描多个工作点，生成极化曲线和功率密度曲线，这在电堆性能工程评估中已经足够。
5. 仿真结果怎么看：典型温度场与浓度场特征
5.1 顺流/逆流布置下的温度分布差异
流道布置方式是SOFC电堆设计里一个关键自由度，最常见的对比是顺流与逆流。
顺流布置下，空气和燃料从同侧进入，沿着流道同向流动。顺着流动方向，反应物浓度都在下降，反应产物在积累，可逆电压逐渐降低；同时反应产热使温度沿着流道方向上升。热源最大值通常出现在入口附近，因为那里反应速率最高。顺流的温度场相对平滑，峰值温度出现在某个中间略靠前的位置，整体分布比较稳定。
逆流布置则不同，空气和燃料从相反方向进入。这样做的优点是换热效率更高，温度场更均匀，理论上可以获得更小的最大温度梯度。实际表现是：燃料入口处空气正在出口，局部温度可能偏低；中央区域容易形成较均匀的中高温区。逆流的峰值温度往往比顺流低一些，电堆寿命也通常更长。
但逆流也有它的麻烦。仿真中经常看到，逆流时高电流密度区域会被“压缩”到流道中段，局部热源密度大，处理不好容易出现中段温度峰值。顺流的情况则是尾部温度偏高、入口温度偏低，温差较大。
从工程经验来说，如果追求高温差下运行稳定，逆流更有优势；如果建模和实验对标方便，顺流更简单直观。我建议初学者先做顺流，跑通整个流程再尝试逆流，两种工况下温度场的对比本身就是一个很有价值的分析点。
5.2 电流密度与温度场的相互作用特征
仿真里最容易观察到的一个现象是：电流密度高的区域，温度一定偏高。这不是巧合，而是电流密度高说明那里反应速率快、产热多；产热多导致温度高；温度高又使电导率和反应速率进一步提升。这个正反馈是理解一切SOFC温度场特征的基础。
实际模型中，电流密度分布并不均匀。靠近流道入口的电极区域，反应物浓度最高，可逆电压最大，所以局部电流密度一般偏高。于是温度场也在这个区域出现一个峰值。如果再加上接触电阻的局部差异（比如肋片和电极之间的接触压力不均），还会出现更细小的局部热点，这种热点在宏观热像仪下都不一定能清晰分辨，但在仿真网格上可以看得清清楚楚。
观察电流密度和温度场耦合时，要格外注意一个量：局部燃料利用率。某些局部区域电流密度太高，燃料消耗速度超过扩散补给速度，局部燃料浓度就会很低。这时该区域的产热反而可能因为浓差极化增大而上升，形成“低浓度高热点”的异常现象。如果你在仿真结果里看到某个区域的温度峰值和燃料浓度谷值重合，就要警惕燃料分配是否合理。
5.3 如何判断模型是否可信：网格无关性与文献对比
任何一个仿真模型的终点不是跑出漂亮云图，而是确定结果可信。我判断SOFC五层电堆模型是否可信，有两个硬指标。
第一个是网格无关性验证。把网格数量从粗到细至少增加一倍，重新计算同一个工况，对比极化曲线和温度峰值。如果两者差异在1%-2%以内，说明网格已经足够；如果差异明显，必须继续加密。特别注意温度峰值对网格的敏感性——温度峰值最容易受网格分辨率影响，尤其是电解质层厚度方向网格不够时，峰值温度会被明显高估。
第二个是极化曲线对标。将模型预测的电压-电流密度曲线与公开实验数据对比，常温等温模型阶段目标误差应控制在5%以内；全耦合模型的极化曲线会有些许偏移，因为温度场分布改变了平均工作温度。对比时要特别注意实验测试的温度区间和模型设定的温度边界是否一致。
我见过太多人写完代码不验证，直接拿结果写报告。这不是仿真，是猜谜。网格无关性验证和对标实验数据这两步，无论多忙都应该做。即使文献里的实验条件与你模拟的工况不完全一致，也能提供量级上的参考，至少能暴露物性参数或边界条件设置的大问题。
6. 常见问题与排查技巧实录
6.1 迭代不收敛，先别急着调网格
遇到不收敛，新手第一反应是加密网格，结果往往更不收敛。正确的排查顺序应该是：先看残差曲线形态，判断是振荡型不收敛还是发散型不收敛。
振荡型不收敛（残差上下波动不下降）多半是阻尼因子设置不当，或者是能量方程与电化学方程之间的耦合过于紧密。对策是降低阻尼因子，或改用“等温-全耦合”的逐步求解路径。发散型不收敛（残差直线飞升）多半是初始值太离谱，比如初始温度场给的均匀值和入口边界温度相差太大，或者热源初始赋值过大。对策是先用等温解做初始值，让热源逐步加进去，不要一上来就全耦合。
此外，边界条件突变也是不收敛的常见原因。流道入口速度和组分在首个计算步内突然建立，会在入口附近造成数值振荡。这时可以通过逐步增加入口流速的“延续法”来解决，先算低流速工况，再以该解为初值加载高流速工况。
6.2 温度场出现不合理峰值，优先查热源和接触热阻
温度峰值位置异常，通常不是传热模块的问题，而是热源分布出了问题。一个典型的场景是：峰值出现在肋片边缘，温度梯度极大。排查时第一步检查局部电流密度分布，如果电流密度本身就在肋片边缘聚堆，那就是电学接触的建模有问题。
连接体和电极之间的接触电阻，在模型中常常被忽略，但实际电堆中接触电阻不能忽视。接触电阻大处局部欧姆热大，温度自然升高。可以把接触电阻等效成一层极薄的附加电阻层，在相应界面处设置接触热阻和接触电阻。很多时候温度场的“鬼峰”消失了，就是对接触电阻做了正确赋值。
另一个常见原因是对流换热系数设得太低。如果电堆外壁的对流换热系数设成极低值，热量散不出去，整个电堆的温度被整体推高，边缘处甚至可能出现不合理的隆起。此时把外壁对流换热系数调回合理范围（5-20 W/(m²·K)），温度场会立刻恢复合理形态。
6.3 参数敏感性分析：最值得多花时间的环节
建模的最终目标不是得到一组结果，而是搞清楚哪些参数对性能影响大。我强烈建议在模型跑通之后，做一轮参数敏感性分析，花的时间值回票价。
优先做这几个参数的敏感性：电解质电导率（或温度Arrhenius参数）、电极孔隙率和曲折因子、流道入口流速、接触电阻、外壁对流换热系数。方法很简单，把单个参数分别上调或下调10%-20%，记录输出电压和最大温度梯度的变化百分比。
从我的经验看，电解质电导率和入口温度的影响几乎总是最大的，这两个参数直接决定基本性能；孔隙率和曲折因子的影响在低电压高电流工况下会被放大；接触电阻的影响则取决于工艺水平，有时可以忽略，有时能占到总内阻的20%以上，完全取决于电堆封装质量。
敏感性分析结果也能指导实验：如果模型显示孔隙率是主要限制因素，那材料改性方向应该优先做孔结构优化，而不是一昧提升反应动力学。这种“仿真反哺实验”的思路，是SOFC多物理场仿真的真正价值所在。
6.4 一个容易被忽视的细节：水蒸气的物性
最后分享一个小细节。阳极侧的组分是氢气和高温水蒸气的混合物。水蒸气在高温下的比热容、粘度、扩散系数与氢气差异极大，而且高温水蒸气容易凝结的顾虑在计算时反而容易让人忽略它在气态下的传热效应。在计算混合气体的有效热导率时，如果忽略水蒸气的高比热容，阳极侧温度会被低估，进而影响阳极过电位的预测精度。
我踩过这个坑。有一版模型算出来的阳极极化偏小，跟实验对不上，查了三天发现是混合气体热物性模型里漏了水蒸气项。从那以后，所有混合气体物性我都用完整的多组分模型，不再做“简化为纯氢气”的假设。细节决定成败，这句话放在SOFC仿真里最贴切不过。
我个人在实际操作中的体会是：这个模型的价值不在于算得多么花哨，而在于帮你看清电堆内部那些实验测不到的分布——温度峰值在哪、燃料在哪不够用、热量从哪里散不掉。把这些看清楚了，设计改进的方向自然就出来了。后面如果你想往深了走，还可以在这个模型基础上叠加热应力分析、材料退化模型，甚至和系统级的启动策略联合仿真，那会是另一番天地了。