COMSOL锂电池热管理多物理场仿真——一个Pack级热失控预警项目的回顾

这个项目要回答的问题很直接：一个36串10并的21700圆柱电芯Pack，在1C快充条件下，最高温度会到多少？如果某一颗电芯发生内短路，热失控传播到相邻电芯需要多长时间？Pack的冷却流道设计能不能在热失控扩散前把温度压下来？

电化学-热-流体三场耦合，再加上热失控的副反应动力学——这已经不是”调一个参数跑一个结果”的节奏了。耦合的收敛性、副反应参数的标定、以及DOE设计空间的遍历——每一个环节都有一段不太顺畅的经历。

电化学-热耦合：P2D模型的网格与收敛

COMSOL的电池模块内置了Newman的P2D（pseudo-2D）模型，在电极厚度方向（x维）解析锂离子浓度和电势分布，在电极平面方向（r-θ维）做热传导。锂离子电池的P2D模型本身是一个刚性很强的DAE系统，加温度耦合后，锂浓度→电势→反应热→温度→扩散系数→锂浓度——这个闭环让收敛更加困难。

这个项目在网格划分上踩过一次坑：初期为了省计算量，电极厚度方向只画了8个单元。结果在2C倍率下，电极表面的电解液锂盐浓度出现了负值——不是物理的真实耗尽，而是网格太粗导致浓度梯度被数值弥散抹平了。

解法很直接但代价很大：电极厚度方向从8个单元加到20个（正极12+隔膜4+负极4），同时把时间步从固定0.1 s切换到自适应BDF（后向差分公式），允许求解器在浓度梯度陡峭时自动缩小步长。收敛了——但单个充放电循环的计算时间从8分钟跳到了45分钟。

做工程仿真有一个很现实的判断要下：你是要算得准，还是要迭代得快。这个项目选择了准——因为最终结论关系到冷却系统的选型成本，多等的几十分钟比起选错一套液冷方案的代价，不值一提。

热失控副反应动力学：参数从哪来？

COMSOL的热失控模型集成了四步副反应：SEI膜分解→负极-电解液反应→正极-电解液反应→电解液分解。每一步的动力学参数（指前因子A、活化能Ea、反应焓ΔH）需要从DSC/TGA实验数据中拟合。

这个项目没有自己跑DSC——用的是文献中NCM811/石墨体系的ARC（加速量热仪）数据。但问题来了：同一篇文献里的SEI分解活化能，不同拟合方法算出来的值差15 kJ/mol——这个差异在Arrhenius公式里对应的反应速率差超过两倍。

项目对四步副反应的动力学参数做了敏感性分析：把每个反应的Ea上下浮动10%，跑瞬态热失控传播模拟。结果很有指向性——前两步（SEI分解和负极反应）对传播速度的影响远远小于第三步（正极释氧反应）。正极反应的Ea变10%，热失控传播到相邻电芯的时间差出近40秒。

这意味着：正极副反应的动力学参数是你最需要标定准的。SEI和负极反应的参数用文献值±10%的误差，对工程判断的影响在可控范围内。

冷却流道优化的DOE设计

Pack的蛇形冷却流道有四个设计变量：流道宽度、流道间距、入口流速、入口温度。全因子实验需要跑81组仿真（3⁴全水平）——COMSOL多物理场单组就要跑两小时以上，81组是无法接受的。

项目用Box-Behnken设计（BBD）把仿真组数压到了27组——保留了二阶交互项的估计能力，但去掉了全因子中的极端组合点。27组跑完，拟合出响应面，最优参数组合落在：流道宽度3.0 mm、间距8.5 mm、流速0.3 m/s、入口温度25°C。

这个最优组合下的仿真结果：1C快充结束时Pack最高温度42.8°C，最大温差3.2°C——在锂电池热管理的安全窗口内（<45°C最高温，<5°C温差）。

但响应面的真正价值不只是找出最优——它把四个变量对温度和温差的敏感性排了序：入口流速 > 流道宽度 > 入口温度 > 流道间距。如果你只有一个变量可以调，调流速的回报最高。这种排序信息比”最优解”本身更实用——实际的Pack设计很少能完全按最优参数走，知道哪个变量是主力、哪些可以妥协，才是工程判断的核心。

这项目没做完的事

热失控传播模拟在27组参数下全部跑完了——但有一个假设没有验证：电芯之间的填充材料（绝缘泡棉）在高温下的热导率变化。仿真里用的是一个常数值0.05 W/(m·K)，但泡棉在100°C以上会热收缩，间隙增大，热传导路径从”泡棉传导”变成”空气间隙+辐射”——等效热导率可以掉到0.02以下。

如果真有热失控发生，这个常数值假设可能导致传播时间被低估。这是项目的遗留问题，后续如果要真的做热失控风险评估，需要先做泡棉高温热导率的实验标定。

引用来源：Newman & Tiedemann, JES 1993（P2D模型原始文献）；Feng et al., Energy Storage Materials 2018（锂电热失控副反应动力学综述）；COMSOL Battery Design Module文档。