4. Fluent求解器基础与热管理控制方程

新能源汽车热管理仿真的本质,是在电池包、电机、电控、乘员舱等复杂几何中求解流场与温度场。Fluent 作为压力基/密度基通用 CFD 求解器,其底层依赖三大守恒方程:质量守恒、动量守恒与能量守恒。理解这三个方程在热管理场景中的物理含义,以及求解器类型、稳态/瞬态、松弛因子与收敛判据的搭配逻辑,是后续所有仿真设置的基础。

一、三大控制方程在热管理中的物理含义

质量守恒方程(连续性方程)描述的是:单位时间内流入控制体的质量减去流出的质量,等于控制体内质量的变化率。在电池液冷板流道中,它保证冷却液不会凭空产生或消失;若连续性残差长期不下降,往往意味着网格存在负体积、边界条件质量不平衡或进出口流量不匹配。

动量守恒方程(Navier-Stokes 方程)描述的是:控制体内动量的变化率等于流入流出的动量净通量、压力梯度力、黏性力与体积力之和。在热管理仿真中,动量方程决定了冷却液在流道内的速度分布、弯头处的二次流、以及风扇/泵驱动下的压降特性。对于液冷板,动量方程与能量方程的耦合决定了换热系数,进而影响电池模组的温度均匀性。

能量守恒方程描述的是:控制体内能量的变化率等于对流带入带出的能量、导热传入传出的能量、以及源项(如电池产热、电机铜损铁损)产生的能量。在电池包热管理中,电芯的产热率通常以体积热源形式加入能量方程;在电机冷却中,绕组铜损作为热源,冷却液通过对流带走热量。能量方程是热管理仿真的核心,其收敛质量直接决定温度场预测的可信度。

三方程并非独立求解,而是通过压力-速度耦合(如 SIMPLE、SIMPLEC、Coupled 算法)与温度-密度耦合(如 Boussinesq 假设或变物性)相互关联。热管理仿真中常见的浮升力驱动自然对流、温度依赖的黏度与密度,都会让三方程耦合更强,收敛难度上升。

二、压力基与密度基求解器的选择

Fluent 提供压力基(Pressure-Based)与密度基(Density-Based)两类求解器。压力基求解器以压力为基本变量,通过压力修正方程满足连续性,适用于不可压或弱可压流动,是新能源汽车热管理仿真的默认选择。液冷、风冷、乘员舱空调等场景,马赫数远低于 0.3,密度变化主要来自温度而非压力,压力基求解器配合变物性或 Boussinesq 假设即可准确捕捉浮升力效应。

密度基求解器以密度为基本变量,通过状态方程耦合压力与密度,适用于高速可压流动、激波捕捉与强压缩性场景。在热管理领域,密度基较少使用,但在某些高速电机风冷、压缩机内部流动或制冷剂两相流近似中可能涉及。若用户误用密度基求解低速液冷问题,会出现收敛慢、压力场振荡甚至发散。

选择建议:电池液冷、电机水冷、乘员舱风道、空调箱内部流动,一律优先压力基;只有明确存在高马赫数可压效应或需要耦合激波时,才考虑密度基。压力基内部还可选择分离式(Segregated)或耦合式(Coupled)压力-速度耦合:分离式内存占用低、适合大网格;耦合式收敛更快、适合强耦合浮升力或旋转机械场景。

热管理仿真求解器选择与三方程耦合关系 三大守恒方程 质量守恒 流入 − 流出 = 质量变化率 保证流道流量平衡 动量守恒 对流 + 压力 + 黏性 + 体积力 决定速度场与压降 能量守恒 对流 + 导热 + 源项 决定温度场与换热 压力-速度耦合 SIMPLE / SIMPLEC Coupled 温度-密度耦合 Boussinesq / 变物性 浮升力驱动自然对流 旋转机械 MRF 求解器选择决策 压力基求解器 Ma < 0.3,不可压/弱可压 密度基求解器 高速可压、激波 液冷板 / 电机水冷 乘员舱风道 / 空调箱 高速电机风冷 压缩机内部流动 热管理仿真默认:压力基 + Coupled 稳态 vs 瞬态设置 稳态 恒定工况、连续散热 液冷板额定流量压降 瞬态 启停、变工况、升温 电池快充温升过程 瞬态时间步:Δt = 0.5~2 s(液冷) 稳态:伪瞬态 + 自动时间步加速收敛 松弛因子:压力 0.3、动量 0.7、能量 0.9 收敛判据:连续性 1e-4、动量 1e-4、能量 1e-6;同时监控电池最高温与压降是否稳定

三、稳态与瞬态的适用场景

稳态求解假设流场与温度场不随时间变化,适用于连续恒定工况下的热平衡分析。例如:液冷板在额定流量 8 L/min、入口温度 25 ℃、电池产热 500 W 恒定条件下的压降与温度分布;电机在额定功率下持续运行的水冷套换热能力评估。稳态计算量小、收敛快,是热管理概念设计与参数扫描的首选。

瞬态求解则用于工况随时间变化的场景:电池快充过程中产热率随 SOC 变化、车辆起步/爬坡时电机功率阶跃、空调系统启停导致乘员舱温度波动、冬季低温加热策略下电池包升温过程。瞬态设置需要定义时间步长、总时长、内迭代次数与数据保存频率。液冷瞬态通常取 Δt = 0.5~2 s,内迭代 20~50 次;若时间步过大,会抹平温度响应细节,甚至导致能量方程发散。

实际操作中,常先用稳态求解获得合理初场,再切换瞬态继续计算,可显著缩短达到稳定周期的时间。Fluent 支持在稳态收敛后直接转为瞬态,保留已有流场与温度场作为初始条件。

四、松弛因子、收敛判据与残差监控

松弛因子控制每次迭代中变量更新的幅度。压力基分离式求解器默认压力松弛因子 0.3、动量 0.7、能量 0.9。若残差振荡或发散,可先将压力降至 0.2、动量降至 0.5,待残差稳定后再逐步恢复。耦合式求解器使用库朗数(Courant Number)控制,默认 200,强耦合浮升力场景可降至 50~100。

收敛判据通常以残差量级衡量:连续性 1e-4、动量 1e-4、能量 1e-6 是热管理仿真的常用阈值。但残差达标不等于结果正确,必须同时监控物理量:电池最高温度、模组温差、冷却液进出口温升、系统压降。若残差已降至 1e-5 但电池最高温度仍在缓慢上升,说明能量方程尚未真正平衡,应继续迭代或检查热源与边界条件。

残差监控建议:在求解过程中实时绘制残差曲线,观察是否单调下降;若出现周期性振荡,可能是流场非定常或松弛因子过大;若残差突然抬升,检查网格质量、时间步长或边界条件是否突变。对于热管理仿真,能量方程残差应比动量方程更严格,因为温度场直接决定电池寿命与热安全评估结论。

教学示例参数(合理假设,非实测):某方形电池模组液冷板,冷却液为 50% 乙二醇水溶液,入口流量 8 L/min,入口温度 25 ℃,电池总产热 500 W,环境温度 25 ℃。稳态计算采用压力基 Coupled 求解器,压力松弛 0.3、动量 0.7、能量 0.9,连续性残差收敛至 8e-5、能量残差 3e-7,电池最高温度约 38.5 ℃,模组温差约 4.2 ℃,液冷板压降约 12 kPa。该组数据仅用于说明设置与判读逻辑,实际结果需以真实网格与边界条件计算为准。

五、热管理仿真中的常见设置误区

误区一:所有工况都用稳态。电池快充、低温加热、空调启停等强瞬态过程必须用瞬态,否则会低估峰值温度。误区二:残差达标即停止。应同时检查进出口能量平衡与监测点温度是否稳定。误区三:松弛因子一味调小。过小会拖慢收敛,应先排查网格与边界条件,再适度调整。误区四:忽略浮升力。液冷板在低流量下自然对流影响不可忽略,应开启重力并采用变物性或 Boussinesq 假设。