质子交换膜燃料电池非等温多相准二维动态模型：面向实时控制的耦合传热传质机理研究

为解决质子交换膜燃料电池（PEMFC）动态运行时三维计算流体力学模型计算成本高、零维模型空间分辨率不足的问题，本研究开发了一种非等温多相准二维动态模型。该模型通过耦合热质传输、液相水生成与传输、膜水合作用等关键物理过程，实现了对极化曲线、气体分布、液态水积累和温度剖面的高精度预测。模型计算效率接近实时，为燃料电池控制系统设计、参数研究和状态监测提供了有效工具。

随着全球交通领域碳中和目标的推进，质子交换膜燃料电池（PEMFC）因其高能量密度和高效率成为备受关注的清洁能源技术。然而，燃料电池在动态运行过程中常面临性能衰减和寿命缩短的挑战，其中水热管理是关键难点。过高湿度会导致膜脱水，液态水过量积聚则引发 flooding（水淹），而温度分布不均又会加速材料降解。传统三维计算流体力学（CFD）模型虽能精确描述内部状态，但计算耗时严重，难以用于实时控制；零维或一维简化模型又缺乏空间分辨率，无法捕捉局部 flooding、膜干涸等关键现象。

为解决这一矛盾，Florian Altmann 等人开发了一种非等温多相准二维动态模型，旨在在计算效率与预测精度之间取得平衡。该研究发表于《BUILDING AND ENVIRONMENT》，文章通过耦合气体通道、气体扩散层（GDL）、催化层和质子交换膜（PEM）等多个物理域，建立了包含质量、物种、动量与能量守恒的偏微分方程组。模型创新性地引入线性化时间积分（LIT）与切比雪夫谱配点法，显著降低计算成本，并实现了对多组分扩散、相变传热、电化学反应以及液态水输运的联合模拟。

在方法层面，作者采用准二维建模思路，将三维结构简化为沿气流方向的“切片”连接，通过截面平均获取主要变量的一维分布；在膜水合模型中，考虑电渗拖曳与扩散的共同作用；在两相流处理中，引入广义菲克定律与毛细压力模型描述气液传输；热模型中则同时考虑导热、对流与反应热。模型参数通过粒子群优化算法与高分辨率CFD结果进行匹配，确保了物理合理性。

气体物种分布与温度预测

通过对比AVL Fire™ M的CFD模拟结果，本模型在多种出口压力条件下（1.3/1.5 bar、2.3/2.5 bar、3.3/3.5 bar）均能准确捕捉阴极与阳极通道内O₂、H₂O、N₂、H₂的摩尔分数变化趋势。尽管在低电流密度下因忽略反应物穿透导致一定偏差，但在主要工作区间（电流密度＞2000 A·m⁻²）预测误差控制在±25%以内。

液态水分布与膜水合状态

在阴极GDL内，模型成功再现了液态水饱和度（s_l）沿通道方向的积累特征。虽然局部数值与CFD结果存在差异，但其空间分布形态和随压力升高而增加的趋势与高精度仿真一致。膜内水含量（λ）的模拟则进一步显示，随着运行压力提高，膜水合程度显著改善，质子传导性增强，与Springer等人提出的水合模型吻合。

温度场与动态响应

非等温模拟表明，燃料电池温度分布受焦耳热、反应熵热和相变潜热共同影响。在阶跃电压实验中，模型捕捉到了温度超调现象，特别是在阴极催化层，因水吸附放热导致初始温度骤升。动态仿真进一步揭示，在高湿度条件下，电流密度响应存在明显延迟，时间常数可达百秒量级，主要受液态水传输过程控制。

研究结论指出，该准二维模型兼具计算效率与物理完备性，能够准确预测PEMFC在稳态与动态工况下的水热分布、物种浓度及极化行为。其意义在于为燃料电池实时控制、系统优化和寿命预测提供了可靠的数字孪生平台。此外，模型中采用的线性化与模型降阶方法，对能源电化学领域的多物理场建模具有普适参考价值。