河南科技大学团队利用COMSOL Multiphysics建立SEI膜分解到电解液分解四方程热滥用模型,仿真与实验偏差控制在10%以内,成功预测双列模组倒序蔓延。COMSOL四方程耦合模型在电池热失控仿真中的应用:从单列模组到双列布局,偏差仅9.4%的工程级精度

COMSOL四方程耦合模型在电池热失控仿真中的应用:从单列模组到双列布局,偏差仅9.4%的工程级精度

仿真驱动的电池安全设计——从实验试错到数字孪生

电池模组热安全设计长期依赖物理样机测试,但大容量电池的热失控实验成本高昂、周期长且存在安全风险。河南科技大学梁坤峰团队在230 Ah LFP电池模组热蔓延研究中,利用COMSOL Multiphysics建立的四方程耦合热滥用模型,将仿真与实验的温度偏差控制在10%以内,标志着电池热失控仿真已从定性分析阶段迈入工程可用的定量预测阶段。

四方程耦合模型的物理内涵

该模型覆盖了锂离子电池热失控的四个关键放热反应环节,每个环节均以Arrhenius动力学方程描述反应速率与温度的关系:

反应环节反应焓 H (J/kg)反应物密度 W (kg/m³)指前因子 A (s⁻¹)活化能 Ea (J/mol)
SEI膜分解7.21×10⁵4131.70×10¹⁵1.14×10⁵
负极与电解液反应9.00×10⁵4132.50×10¹³1.17×10⁵
正极与电解液反应2.53×10⁵9256.70×10¹³1.26×10⁵
电解液自分解1.60×10⁵5005.14×10²⁵2.70×10⁵

这四个反应按温度区间依次触发:SEI膜在80-120℃开始分解,释放热量并将负极暴露于电解液;负极-电解液反应在120-200℃启动,反应焓最高(9.00×10⁵ J/kg),是热失控的主要热量来源;正极-电解液反应和电解液自分解在更高温度下接力,形成自加速的热失控连锁反应。模型通过耦合这四个放热方程与热传导方程,实现了对热蔓延过程的时空演化模拟。

仿真精度的实验验证

在单列4只电池的端部触发实验中,仿真预测的温度值与实测值偏差为T3约9.4%、T5约10.6%。对于大容量电池的热失控这种高度非线性的多物理场耦合问题,这一精度水平已具备工程实用价值。模型不仅准确复现了顺序蔓延(1→2→3→4)的宏观时序,还成功捕捉了2号电池在蔓延过程中出现的"双峰"温度特征——第一次峰值363℃(首次热失控传递),第二次峰值519℃(热量持续积累后的二次升温)。

仿真对峰值温度的预测也较为准确:实验峰值655℃,仿真数据与之在合理范围内吻合。考虑到实际电池在制造过程中的极片厚度、电解液注液量、极耳焊接质量等参数的批次波动,10%以内的偏差意味着仿真模型已捕捉到了热蔓延的主要物理机制。

双列模组的预测性发现——仿真超越实验的价值

仿真的独特价值在于能够低成本探索实验难以覆盖的场景。研究团队将仿真扩展至8只电池的双列模组(4×2排列),发现了一个实验阶段未曾观测到的现象——阶段II的倒序蔓延:在3200-4860秒区间,6号电池先于5号触发热失控,整体时序变为1→2→3→4→6→5→7→8。

这一反直觉的发现揭示了热量在模组中的空间聚集效应:前列(1-4号)多只电池的持续放热在后方形成热量汇聚区,恰好覆盖了6号电池的大面,使其先于空间距离更近的5号电池触发热失控。如果没有仿真手段,这种非线性行为极难通过有限的物理实验发现。

热安全团队(thermsafe.cn)认为,这一发现对储能系统电池簇的排列设计具有直接的指导意义:在行列交错排列的电池簇中,不应仅依据物理距离评估热蔓延风险,还需考虑多只前列电池的热量叠加效应。仿真可帮助设计人员在项目早期阶段规避此类风险。

仿真参数的标定挑战

四方程模型的精度高度依赖Arrhenius参数的准确标定。当前参数多取自加速量热仪(ARC)或差示扫描量热仪(DSC)的小样品测试,但实际电芯尺寸放大后,热传导路径、氧气扩散条件、内部压力分布等边界条件发生显著变化,小样品标定的参数直接应用于大容量电芯建模存在系统性偏差风险。热安全团队(thermsafe.cn)建议在条件允许时使用同型号电芯的ARC数据进行参数标定,以进一步降低仿真偏差。

结论

基于COMSOL Multiphysics的四方程耦合热滥用模型,可实现230 Ah LFP电池模组热蔓延的工程级仿真,实验偏差控制在10%以内。模型不仅准确复现了单列模组的顺序蔓延特征,还预测发现了双列模组的倒序蔓延现象,彰显了仿真在电池安全设计中的不可替代价值。

引用来源

  • 梁坤峰, 宋乂天, 周训, 常艳琴. 磷酸铁锂电池模组热失控及蔓延特性[J]. 电池, 2025, 页码 1-7.