电池热管理系统的多物理场耦合仿真

一个36颗18650电池组成的模组，客户要求仿真验证液冷方案在3C放电条件下的温升是否超过安全阈值（电池表面温度≤55°C）。问题描述本身不复杂——电池发热、冷却液带走热量、温度场分布、压力损失，四个物理量串在一起。但真正上手之后发现，串在一起容易，算在一起几乎不可能——因为电化学场的时间尺度（毫秒级）、热传导场（秒级）和流场（毫秒到秒）差了四到六个数量级。

时间尺度错配是根本问题

电化学反应在电极/电解液界面上的电荷转移过程发生在微秒到毫秒量级。电池单体的热响应时间常数在分钟级（一颗18650的集总热容约60 J/K，3C下的发热功率约2.5W，温升特征时间约500-1000秒）。冷却液在流道里的对流换热时间约几秒。

如果直接做瞬态的全耦合计算：电化学场需要微秒级步长来捕捉界面动力学，热场需要秒级步长来覆盖整个放电循环，流场在两者之间。微秒级的电化学步长去跑3000秒（~50分钟）的3C放电，需要30亿个时间步——在现有计算资源下不现实。

解决方案不是把步长统一到一个值，而是分离求解——把三个场拆成两组：电化学+热做第一组（时间尺度匹配，都是相对慢的过程），流体做第二组（准稳态近似），两组之间通过间歇的数据交换来耦合。

电化学-热耦合：发热源项的分布化

第一步是把电池的发热源项从电化学模型里抽出来。用Newman的P2D（伪二维）电化学模型把沿电极厚度方向的电流密度、过电位和欧姆压降算出来，然后对电极厚度积分得到每个电池单元的发热功率密度（单位W/m³）。

P2D模型的核心参数来自电池厂提供的电极配方数据：正极NMC811（厚度70μm，孔隙率0.32），负极石墨（厚度85μm，孔隙率0.35），电解质1M LiPF₆/EC:DMC。用COMSOL的锂电池模块跑P2D，在3C倍率下预测的单体平均发热功率约2.6W——可逆热（熵变）约0.4W，不可逆热（欧姆+极化）约2.2W。不可逆热集中在负极-隔膜界面附近（锂离子浓度梯度最大的区域），因此在热传导模型里不能用一个均匀体热源，需要在电极叠层方向给热源做分层分布。

把P2D输出的热源空间分布映射到三维热传导模型的电池域上，热传导方程用Fourier定律，电池的导热系数取各向异性（轴向约30 W/m·K，径向约1.0 W/m·K——圆柱电池的卷绕结构导致轴向导热远好于径向）。这个映射之后，电化学场就退出了瞬态迭代——后续的瞬态热传导只需要秒级的步长，脱离了毫秒级的电化学约束。

流体-热耦合：准稳态假设的边界

冷却液（50%乙二醇水溶液，入口温度25°C，流速0.5 m/s）在冷板流道里流动，对流换热系数和流动压力损失需要CFD来算。但如果把瞬态CFD耦合上去，计算量又炸了——CFD稳态求解一个工况需要几分钟到十几分钟，乘以热传导的几千个时间步，不可行。

这里用了一个准稳态假设：在每一个热传导时间步内（如10秒），冷却液温升不超过1°C，流动状态可以认为在10秒内是准稳态的。所以每10个热传导时间步（即100秒）更新一次CFD求解，把更新后的对流换热系数和冷却液温度作为边界条件回传给热传导模型。

100秒的耦合间隔选得有物理依据：50%乙二醇溶液以0.5 m/s流过整个冷板流道（长约0.3m）所需时间约0.6秒，10秒内流道内的冷却液已经刷新了16轮，热交换充分。100秒对应温升约0.5-1°C——在这个温升区间内对流换热系数的变化（主要来自黏度随温度变化）不到2%，准稳态假设在这个间隔内是成立的。

结果：温场预测与流道压损

3C放电结束时刻（3000秒），电池模组中心位置的最高表面温度达到了53.4°C——低于55°C的安全阈值，但留的余量不大（1.6°C）。温度分布呈现中心高、边缘低的特征（中心电池被两侧电池的辐射热”保温”），芯-壳温差约8°C（中心电池表面温度比最外侧电池高8°C）。

冷板流道的压力损失约1.8 kPa，在合理的泵功率范围内（对应泵功约0.5W）。但流道的温度均匀性不够理想：入口端冷却液温度25°C，出口端约29°C，沿流道方向形成了4°C的冷却液温差，导致了模组入口侧电池平均温度约45°C、出口侧约49°C的不均匀分布。

优化方向有两个：增加流道并联分支（降低压损+提升温度均匀性），或者在电池之间增加导热硅胶垫片（压低芯-壳温差）。两个优化的代价分别是泵功率增加和热阻增加，需要重新跑一轮耦合计算来定量判断。分离式耦合方法的优势在这里体现出来了——改流道设计只需要重新跑CFD（几小时），不用重新算电化学P2D（几小时到一天），两个场的分离让设计迭代的成本降低了近一个数量级。