在STAR-CCM+中进行热流固耦合(CHT)计算时,工程师常常面临计算资源与精度的权衡。对于包含复杂流动和长时间瞬态传热的问题(如电池包热管理),完全耦合求解流动与能量方程可能消耗大量时间。
为此,**解耦计算**(或称“冻结流动”方法)提供了一种高效的替代方案:在流动达到准稳态后,将流动求解器冻结,仅求解能量方程,从而大幅加速计算。本文以1P52S电池包为对象,对比了耦合与解耦两种计算方法在7200 s和12000 s时刻的温度结果,验证了解耦计算的精度。
耦合计算:同时求解流动控制方程(N-S方程)和能量方程。流动与温度场双向耦合,适用于流动剧烈变化或自然对流占主导的场景,但计算成本高。
解耦计算:在流动基本稳定后,通过“冻结”流动求解器(如将流动松弛因子设为0或直接勾选“冻结流”选项),使速度场不再更新,仅由能量方程求解温度分布。该方法假设后续传热对流动影响可忽略,适用于流动变化缓慢的长时间瞬态问题。

本次验证基于1P52S电池包(即1并52串的电池模组),分别采用耦合与解耦方式计算了7200 s和12000 s时的温度分布。关键温度指标包括:温度最大值、温度最小值、出口水温。结果汇总如下表(单位:℃):

1. 绝对偏差极小:在所有对比项中,解耦与耦合计算的偏差均小于 0.003 ℃。例如7200 s温度最大值的偏差仅0.0023 ℃,出口水温偏差绝对值不超过0.0013 ℃。这表明两种方法在工程关注的温度尺度上几乎无差别。
2. 相对误差可忽略:以温度最大值为例,最大相对误差仅为 0.006%(0.0023/37.87),远低于工程通常允许的1%误差范围。
3. 偏差随时间无明显放大:从7200 s到12000 s,偏差并未累积增大,说明在长时间计算中解耦方法依然保持稳定。
4. 原因分析:电池包内部流动(如冷却液循环)在初始阶段后可能已达到充分发展状态,速度场变化微弱,此时冻结流动对热边界层及换热系数的影响可忽略。因此,仅求解能量方程即可准确捕捉温度演化。
基于1P52S电池包的对比验证表明,在STAR-CCM+中采用解耦(冻结流动)方法进行热流固耦合计算,能够在保证亚毫摄氏度级精度的前提下,显著缩短计算周期。该方法尤其适用于流动基本稳定的长时间瞬态仿真,如电池包恒功率放电、长时间热失控蔓延等场景。建议工程师在应用前先进行短期耦合计算确认流动稳定,然后转为解耦模式以提升效率。
