燃料电池仿真

类别: 分析 | 集成版 2026-04-06

燃料电池的理论基础

🎓

PEMFC/SOFC的电化学-热-流体耦合。Butler-Volmer式。




支配方程




$$ V=E_{Nernst}-\eta_{act}-\eta_{ohm}-\eta_{conc} $$
$$ j=j_0\left[\exp(\frac{\alpha_a F\eta}{RT})-\exp(-\frac{\alpha_c F\eta}{RT})\right] $$



🧑‍🎓

听到这里，我终于明白了燃料电池仿真为什么重要了！


离散化方法

🧑‍🎓

这个方程在计算机上实际怎样求解呢？


🎓

用有限元法（FEM）进行空间离散化。组装单元刚度矩阵，构造整体刚度方程。


🎓

将微分方程转换为弱形式（变分形式），使用试函数和形状函数进行Galerkin方法定义。单元类型的选择（低阶单元 vs. 高阶单元、完全积分 vs. 降低积分）与解的精度和计算成本的折衷直接相关。




矩阵求解算法

🧑‍🎓

矩阵求解算法具体是什么意思？


🎓

通过直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解联立方程。大规模问题采用预处理迭代法效果显著。



求解法分类内存使用量适用规模
LU分解直接法O(n²)小~中等规模
Cholesky分解直接法（对称正定）O(n²)小~中等规模
PCG法迭代法O(n)大规模
GMRES法迭代法O(n·m)大规模·非对称
AMG预处理预处理O(n)超大规模
🧑‍🎓

也就是说，有限元法的部分如果做得不好，后面会吃亏。我记住了！


商用工具中的实现

🧑‍🎓

进行燃料电池仿真可以用什么软件呢？


工具名称开发者/现在主要文件格式
COMSOL MultiphysicsCOMSOL AB.mph
Ansys FluentANSYS Inc..cas, .dat, .msh, .jou
Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
ANSYS Mechanical (原ANSYS Structural)ANSYS Inc..cdb, .rst, .db, .ans, .mac

厂商谱系与产品整合历程

🧑‍🎓

各个软件的成长经历，有时是很有戏剧性的呢？



COMSOL Multiphysics

🧑‍🎓

请介绍一下「COMSOL Multiphysics」！


🎓

1986年在瑞典成立。作为MATLAB联动的FEMLAB起步，后来改名为COMSOL。多物理学方面强势。

现在的归属: COMSOL AB



Ansys Fluent

🧑‍🎓

接下来讲Ansys Fluent的内容吧。是什么样的呢？


🎓

由Fluent Inc.开发。2006年被ANSYS收购。基于非结构化网格的通用CFD求解器。

现在的归属: ANSYS Inc.




Simcenter STAR-CCM+

🧑‍🎓

接下来讲Simcenter STAR的内容吧。是什么样的呢？


🎓

由CD-adapco开发。2016年被西门子收购，纳入Simcenter品牌。多面体网格是其特色。

现在的归属: Siemens Digital Industries Software


🧑‍🎓

啊，是这样啊！ 在瑞典成立原来是这样的机制呢。


文件格式与互操作性

🧑‍🎓

在不同软件之间转换数据时有什么注意事项吗？


格式扩展名种类概述
STEP.stp/.step中立CAD符合ISO 10303的3D CAD数据交换格式。形状+PMI兼容。
IGES.igs/.iges中立CAD早期CAD数据交换规范。曲面数据兼容性有问题。正在向STEP迁移。
🎓

在不同求解器之间转换模型时，需要注意单元类型的对应关系、材料模型的兼容性、荷载·边界条件的表示差异。特别是高阶单元或特殊单元（内聚力单元、用户定义单元等）往往无法在求解器间直接转换。


🧑‍🎓

原来...格式看似简单，其实深度非常大呢。


实务注意事项

🧑‍🎓

教科书里没有的"现场的智慧"之类的东西有吗？


🎓

网格收敛性确认、边界条件的合理性验证、材料参数的敏感性分析非常重要。


🎓
  • 网格依赖性的验证: 至少用3个网格密度等级确认收敛性
  • 边界条件的合理性: 设置物理上有意义的约束条件
  • 结果的验证: 与理论解、实验数据、已知基准问题的比较


    • 🧑‍🎓

    燃料电池仿真的整体情况我已经掌握了！ 从明天开始我会在实务中认真应用。


    🎓

    嗯，很不错！ 实际动手操作才是最好的学习。有不懂的地方随时问我。


    Coffee Break 闲聊一下

    固体氧化物燃料电池（SOFC）为什么能像"热引擎"一样运作

    燃料电池仿真的理论魅力的象征是固体氧化物燃料电池（SOFC）的高温运作。与PEFC在80℃运作不同，SOFC在600~1,000℃下运作。这个高温看似"弱点"，其实是优点，温度越高，活化过电压越小（阿伦尼乌斯定律），不需要铂这样的昂贵催化剂。另外，废热温度高，与燃气轮机组合的联合循环（SOFC-GT）理论效率可达60~70%。仿真中热应力也很重要，温度急剧变化导致陶瓷电解质破裂是SOFC实用化的最大课题之一。"电化学×高温固体力学的耦合"——SOFC仿真正是多物理学的范例。

    燃料电池的数值计算方法

    数值方法的细节

    🧑‍🎓

    具体采用什么算法来求解燃料电池仿真呢？




    离散化的表述



    🎓

    用形状函数 $N_i$ 近似未知量:



    $$ u^h(\mathbf{x}) = \sum_{i=1}^{n} N_i(\mathbf{x}) \, u_i $$




    🎓

    用公式表示就是这样。


    $$ K_e = \int_{\Omega_e} B^T \, D \, B \, d\Omega \approx \sum_{g=1}^{n_g} w_g \, B^T(\xi_g) \, D \, B(\xi_g) \, |J(\xi_g)| $$

    基础方程离散形式


    🎓

    用公式表示就是这样。


    $$ V=E_{Nernst}-\eta_{act}-\eta_{ohm}-\eta_{conc} $$
    $$ j=j_0\left[\exp(\frac{\alpha_a F\eta}{RT})-\exp(-\frac{\alpha_c F\eta}{RT})\right] $$

    🧑‍🎓

    嗯，只有公式我还是有点不理解...那是表示什么呢？


    🎓

    连续体的控制方程离散化后，得到以下代数方程组:



    $$ [K]\{u\} = \{F\} $$


    🎓

    其中$[K]$是整体刚度矩阵（或等效的系统矩阵），$\{u\}$是未知节点变量向量，$\{F\}$是外力向量。


    🧑‍🎓

    啊，是这样！ 连续体的控制方程离散化原来是这种机制呢。


    单元技术

    🧑‍🎓

    我听说过"单元技术"这个说法，但可能没完全理解...。


    单元类型次数节点数(3D)精度计算成本
    四面体一次线性4低（剪切锁定）低
    四面体二次二次10高中
    六面体一次线性8中中
    六面体二次二次20非常高高
    棱柱线性/二次6/15中~高中

    积分方案

    🧑‍🎓

    积分方案具体是什么意思？


    🎓
    • 完全积分: 所有项都精确积分。刚度被过度估计的倾向（锁定）
    • 降低积分: 减少积分点数。提高计算效率，但有沙漏模式风险
    • 选择降低积分 (B-bar法): 体积项和偏差项分开积分。避免锁定

      • 🧑‍🎓

      听到这里，我终于明白了单元类型为什么重要！


      收敛性与稳定性

      🧑‍🎓

      不收敛了的时候，首先应该检查什么？


      🎓
      • h-细化: 将网格细分化（使单元大小h变小）来提高精度
      • p-细化: 提高单元的多项式次数来提高精度
      • hp-细化: 同时优化h和p

        • 🎓

        收敛速度: 二次单元的误差按 $O(h^2)$ 的阶数减少（光滑解的情况）


        🧑‍🎓

        原来...网格细分化看似简单，其实深度非常大呢。


        求解器设置建议

        🧑‍🎓

        具体采用什么算法来求解燃料电池仿真呢？


        参数推荐值备注
        迭代法收敛判定$10^{-6}$残差范数标准
        预处理手法ILU(0) or AMG根据问题规模
        最大迭代次数1000非收敛时要重新调整设置
        内存模式In-core尽可能使用

        单一方法

        将全部物理场作为1个联立方程系统同时求解。对强耦合稳定，但实现复杂，内存消耗大。

        分区法（分离迭代法）

        各物理场独立求解，在界面处进行数据交换。实现容易，可利用已有求解器。适于弱耦合。

        界面数据转移

        最近邻法（最简单但精度低），投影法（保守），RBF插值（对网格不匹配强）。保守性和精度的平衡很重要。

        子迭代

        各耦合步骤内进行充分的迭代，确保界面条件的一致性。残差准则应根据各物理场的典型值进行缩放。

        Aitken缓和

        自动调整耦合迭代的缓和系数。防止过缓和导致的发散，加速收敛的自适应方法。

        稳定性条件

        注意附加质量效应（流体-结构耦合时流体密度≈结构密度）。不稳定时应用罗宾型界面条件或IQN-ILS法。

        燃料电池的实务应用

        🎓

        解释燃料电池仿真的实务分析流程和注意事项。



        分析流程

        🧑‍🎓

        从最初开始请教我！ 应该从什么开始呢？


        🎓

        1. 前处理 (Pre-processing)

        • CAD数据的导入和形状简化
        • 材料特性的定义
        • 网格生成（单元类型·大小的决定）
        • 边界条件和荷载条件的设置

        🎓

        2. 求解 (Solving)

        • 求解器设置（解法、收敛准则、输出控制）
        • 作业提交和计算执行
        • 收敛监测

        🎓

        3. 后处理 (Post-processing)

        • 结果的可视化（位移、应力、其他物理量）
        • 结果的验证和合理性确认
        • 报告编制


        网格生成最佳实践

        🧑‍🎓

        怎样判断网格的好坏呢？



        单元品质指标

        🧑‍🎓

        请介绍一下「单元品质指标」！


        指标理想值容许范围影响
        宽高比1.0< 5.0精度下降
        雅可比比1.0> 0.3单元退化
        翘曲0°< 15°精度下降
        歪斜度0°< 45°收敛性恶化
        锥体比0< 0.5精度下降

        网格密度的决定

        🧑‍🎓

        网格密度的决定具体是什么意思？


        🎓
        • 应力集中部: 配置至少3层以上的单元
        • 应力梯度大的区域: 将单元大小设为周围的1/3~1/5
        • 荷载施加点附近: 局部细化
        • 远方区域: 采用粗网格来确保计算效率


          • 边界条件设置指南

          🧑‍🎓

          听说边界条件如果搞错，全部都会完蛋...。


          🎓
          • 过度约束要小心: 刚体移动的约束仅为6自由度
          • 对称条件的活用: 降低计算规模
          • 荷载的等效分配: 集中荷载 vs. 分布荷载的选择

            • 🧑‍🎓

            啊，是这样！ 过度约束要小心原来是这种机制呢。


            各商用工具的实现步骤

            🧑‍🎓

            有各种各样的软件吧？ 请分别介绍特点！


            工具名称开发者/现在主要文件格式
            COMSOL MultiphysicsCOMSOL AB.mph
            ANSYS FluentANSYS Inc..cas, .dat, .msh, .jou
            Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
            ANSYS Mechanical (原ANSYS Structural)ANSYS Inc..cdb, .rst, .db, .ans, .mac

            COMSOL Multiphysics

            🧑‍🎓

            请介绍一下「COMSOL Multiphysics」！


            🎓

            1986年在瑞典成立。作为MATLAB联动的FEMLAB起步，后来改名为COMSOL。多物理学方面强势。

            现在的归属: COMSOL AB



            ANSYS Fluent

            🧑‍🎓

            接下来讲ANSYS Fluent的内容吧。是什么样的呢？


            🎓

            由Fluent Inc.开发。2006年被ANSYS收购。基于非结构化网格的通用CFD求解器。

            现在的归属: ANSYS Inc.


            🧑‍🎓

            老师的解释很容易理解！ 工具名称的困惑解决了。


            常见失败与对策

            🧑‍🎓

            初学者容易犯什么失误呢？ 我想事先了解！


            症状原因对策
            计算不收敛网格品质不良、不适当的边界条件改善网格、重新检查约束条件
            应力异常大应力奇点、网格依赖避免奇点、局部网格细化
            位移非现实的材料常数错误、单位系统不一致确认输入数据
            计算时间过长不必要的细化、低效的求解方法网格优化、并行计算

            质量保证检查清单

            🧑‍🎓

            教科书里没有的"现场的智慧"之类的东西有吗？


            🎓
            • 用3个以上等级确认网格收敛性了吗
            • 确认力平衡（反力合计）了吗
            • 结果在物理的合理范围内吗
            • 与已知理论解或基准问题比较了吗


              • 🧑‍🎓

              燃料电池仿真的整体情况我已经掌握了！ 从明天开始我会在实务中认真应用。


              🎓

              嗯，很不错！ 实际动手操作才是最好的学习。有不懂的地方随时问我。


              Coffee Break 闲聊一下

              巴士·卡车的燃料电池系统——追随电力需求变化的仿真

              燃料电池仿真的实践应用正在急速扩大的是商用车（巴士·大型卡车）的燃料电池系统设计。与乘用车不同，商用车在坡道启动·负载变化·长距离连续运行这样的恶劣条件下运作。仿真发挥强大威力的地方是"对动态负载变化的追随性"的评估，急加速时燃料电池堆的应答跟不上会导致电压下降，对辅助电池的依赖增加。现场通常采用"燃料电池堆+锂离子电池的混合系统"的电力分配控制基于实际行驶数据的仿真优化方法。伊苏兹、日野、UD卡车等国内厂商正在推进FCV卡车的示范运行，这一领域的仿真需求在急速增长。

              燃料电池的软件比较

              🎓

              详述燃料电池仿真中支持的主要商用CAE工具的功能比较和各产品的历史背景。



              支持工具列表

              🧑‍🎓

              进行燃料电池仿真可以用什么软件呢？


              工具名称开发者/现在主要文件格式
              COMSOL MultiphysicsCOMSOL AB.mph
              ANSYS FluentANSYS Inc..cas, .dat, .msh, .jou
              Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
              ANSYS Mechanical (原ANSYS Structural)ANSYS Inc..cdb, .rst, .db, .ans, .mac

              COMSOL Multiphysics

              🧑‍🎓

              请介绍一下「COMSOL Multiphysics」！


              🎓

              1986年在瑞典成立。作为MATLAB联动的FEMLAB起步，后来改名为COMSOL。多物理学方面强势。

              现在的归属: COMSOL AB



              ANSYS Fluent

              🧑‍🎓

              接下来讲ANSYS Fluent的内容吧。是什么样的呢？


              🎓

              由Fluent Inc.开发。2006年被ANSYS收购。基于非结构化网格的通用CFD求解器。

              现在的归属: ANSYS Inc.




              Simcenter STAR-CCM+

              🧑‍🎓

              接下来讲Simcenter STAR的内容吧。是什么样的呢？


              🎓

              由CD-adapco开发。2016年被西门子收购，纳入Simcenter品牌。多面体网格是其特色。

              现在的归属: Siemens Digital Industries Software



              ANSYS Mechanical (原ANSYS Structural)

              🧑‍🎓

              请介绍一下「ANSYS Mechanical」！


              🎓

              1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL（ANSYS参数设计语言）。

              现在的归属: ANSYS Inc.


              🧑‍🎓

              原来…在瑞典成立看似简单，其实深度非常大呢。


              功能比较矩阵

              🧑‍🎓

              预算和时间都有限，成本效益最高的是哪个？


              功能COMSOLFluentStar-CCM+ANSYS Mechanical
              基本功能○○○○
              高级功能○○○△
              自动化/脚本○○○○
              并行计算○○○○
              GPU支持△△△○

              转换时的风险

              🧑‍🎓

              转换时的风险具体是什么意思？


              🎓
              • 单元类型的不兼容: 求解器固有单元无法用中立格式表现
              • 材料模型的差异: 同名但内部实现不同的情况有
              • 边界条件的重新定义: 多数情况需要手动重新设置
              • 结果数据的比较: 输出变量的定义（节点值 vs. 单元值、积分点值）有差异

                • 🧑‍🎓

                啊，是这样！ 不同工具之间的转移原来是这种机制呢。


                许可形式

                🧑‍🎓

                我听说过"许可形式"这个说法，但可能没完全理解...。


                工具许可特点
                商用FEA节点锁定/浮动价格高但有官方支持
                OpenFOAMGPL无偿但支持收费
                COMSOL节点锁定/浮动按模块购买
                Code_AsterGPLEDF开发的OSS求解器

                选择指南

                🧑‍🎓

                最终该选哪个呢，请教我判断标准！


                🎓

                在燃料电池仿真工具选择中，应考虑以下因素:


                🎓
                • 分析规模: 数万~数亿DOF的可扩展性
                • 物理模型: 所需构成关系·单元类型的兼容性
                • 工作流: CAD的联动、自动化的便利性
                • 成本: 初期投资 + 年度保守 + 教育成本
                • 支持: 技术支持的质量与响应


                  • 🧑‍🎓

                  燃料电池仿真的整体情况我已经掌握了！ 从明天开始我会在实务中认真应用。


                  🎓

                  嗯，很不错！ 实际动手操作才是最好的学习。有不懂的地方随时问我。


                  Coffee Break 闲聊一下

                  燃料电池仿真工具的"3个流派"——CFD·系统·专用工具

                  燃料电池仿真的工具大致分为3个流派。CFD系（ANSYS Fluent、Star-CCM+）流动·热的精度高，堆栈内的3D详细分析强势。系统系（MATLAB/Simulink、GT-Suite）电源管理·控制系统的整合强，车载系统的动态仿真强。专用工具（ANSYS燃料电池模块、COMSOL电化学模块）电化学物理模型充实，适合材料·电池设计的研究开发。实际开发现场，经常采用"材料开发用专用工具、详细设计用CFD、系统评估用系统仿真"的分工方式，交换数据进行协作。"自己的课题与哪个流派的工具相性最好"这个理解是选择的出发点。

                  燃料电池的先端研究

                  前沿课题与研究趋势

                  🧑‍🎓

                  燃料电池仿真分野今后会如何发展呢？


                  🎓

                  来看燃料电池仿真领域的最新研究动向和先进技术手法。



                  最新的数值方法

                  🧑‍🎓

                  接下来讲最新数值方法的内容吧。是什么样的呢？



                  🧑‍🎓

                  只有公式我还是有点不理解...那是表示什么呢？


                  🎓
                  • 等参数分析 (IGA): 直接使用NURBS基函数，实现CAD-CAE间的无缝联动
                  • 粒子法 (SPH, MPM): 无网格方法对大变形·破坏的追踪
                  • 相场法 (Phase-Field): 用界面的隐式表现实现复杂界面追踪
                  • 机器学习支援: 代理模型、物理信息神经网络 (PINN)


                    • 高性能计算 (HPC) 的应对


                    并行化手法概述支持求解器
                    MPI (领域分割)分布式内存型。大规模问题的标准全主要求解器
                    OpenMP共享内存型。节点内并行多数求解器
                    GPU (CUDA/OpenCL)GPGPU活用。特别在显式法中有效LS-DYNA, Fluent等
                    混合 MPI+OpenMP节点间+节点内并行大规模HPC环境

                    燃料电池的故障排除



                    常见错误与对策

                    🧑‍🎓

                    老师也在燃料电池仿真中经历过通宵调试吗？（笑）



                    1. 收敛失败

                    🧑‍🎓

                    收敛失败具体是什么意思？


                    🎓

                    症状: 求解器在指定迭代次数内不收敛，异常终止


                    🎓

                    可能的原因:

                    • 网格品质不足（过度扭曲的单元）
                    • 材料参数设置不当
                    • 不当的初始条件
                    • 非线性性过强（荷载步不足）

                    🎓

                    对策:

                    • 进行网格品质检查（宽高比、雅可比比）
                    • 确认材料参数的单位系统
                    • 将荷载分解为多个步骤（增加子步数）
                    • 放宽收敛判定基准（注意精度）

                    🧑‍🎓

                    也就是说，收敛失败的部分如果做得不好，后面会吃亏。我记住了！



                    2. 非物理的结果

                    🧑‍🎓

                    接下来讲非物理的结果的内容吧。是什么样的呢？


                    🎓

                    症状: 应力/位移/温度等出现物理上非现实的值


                    🎓

                    可能的原因:

                    • 边界条件误设置
                    • 单位系统混用（SI单位与工程单位混淆）
                    • 不适当的单元类型选择
                    • 应力奇点的存在

                    🎓

                    对策:

                    • 确认反力合计（力平衡）
                    • 确认单位系统的一致性
                    • 重新检讨单元类型的适切性
                    • 消除或进行子建模的奇点

                    🧑‍🎓

                    前辈说"收敛失败一定要做好"的意思我现在理解了。




                    3. 计算时间超过

                    🧑‍🎓

                    计算时间超过具体是什么意思？


                    🎓

                    症状: 计算花费想定时间的数倍


                    🎓

                    对策:

                    • 网格粗细分布的最优化
                    • 对称性的活用（1/2, 1/4模型）
                    • 求解器设置的最优化（迭代法、预处理的选择）
                    • 活用并行计算



                    4. 内存不足

                    🧑‍🎓

                    请介绍一下「内存不足」！


                    🎓

                    症状: Out of Memory 错误


                    🧑‍🎓

                    前辈说"收敛失败一定要做好"的意思我现在理解了。


                    🎓

                    对策:

                    • 使用超核求解法
                    • 削减网格规模
                    • 确认64位版