流体热仿真：共轭传热问题中流场与温度场的双向耦合策略

流体热仿真在这个动力电池包液冷项目里，最初给出的温度分布图让热管理工程师沉默了许久。项目背景是某纯电动车型的模组级液冷方案验证：冷却板内嵌于电池模组底部，冷却液为体积分数百分之五十的乙二醇水溶液，设计流量八升每分钟，目标是在一点五倍持续放电倍率下，电芯最大温差不超过五摄氏度。初始仿真给出的电芯最高温度约四十一摄氏度，与实验测试的四十七摄氏度相差六度——这个差距在热管理设计中意味着安全余量的严重误判，项目因此被要求重新核算全部仿真设置。

共轭传热的核心难点，在流体与固体交界面上的温度连续条件是否被正确处理。Fluent中默认开启能量方程后，流固界面会自动应用温度连续和热流连续条件，这一点不需要手动设置。但问题出在界面两侧的网格密度不匹配：流体侧第一层网格厚度约零点三毫米，固体侧网格尺寸约二毫米，跨越界面的温度梯度在固体侧只有一层网格来描述，导致温度分布在界面附近出现了不连续台阶。这个被忽略的数值问题，让冷却板到冷却液的传热阻力在仿真中被系统性低估。将固体侧近界面区域的网格尺寸收紧到零点五毫米，同时保证界面两侧网格节点的对齐（共节点网格），温度台阶消失，电芯最高温度的预测值上升到四十五点三摄氏度，与实验值的偏差缩小到约百分之四。

湍流普朗特数（Turbulent Prandtl Number）的设定，是影响流体热仿真精度的隐形因素。Fluent默认将湍流普朗特数设为零点八五，这个数值来源于对空气和水的经典湍流输运实验的归纳。但冷却液的物性与水有差异：百分之五十乙二醇水溶液在二十摄氏度时的普朗特数约二十五，明显高于水的七左右，对应的湍流普朗特数也偏离默认值。项目中引用了Churchill（2002）在International Journal of Heat and Mass Transfer上发表的变物性流体湍流普朗特数关联式，根据冷却液的当地温度和当地速度梯度动态计算湍流普朗特数，而非使用全局常数。这个改动对电芯表面换热系数的分布产生了可观测的影响：高流速区域的换热系数预测值上调了约百分之十二，对应的局部温度下降约一点八摄氏度。

求解顺序的设置，在流热耦合问题中存在一个容易被忽视的策略选择。Fluent提供了两种耦合方式：先后求解（Segregated）和耦合求解（Coupled）。先后求解先解流场至收敛，再固定流场解温度场，适合流速较高、温度场对流场影响可忽略的场景。但电池液冷系统中，冷却液温度沿流动方向可上升约八至十摄氏度，温度变化引起的黏度变化和密度变化会反过来影响流速分布——这个双向耦合效应在先后求解方式中被完全丢失。切换到耦合求解后，流场与温度场在每个迭代步中同步更新，同时开启浮升力修正（Boussinesq近似或完全可压密度模型），温度对流动的反向影响被纳入计算。这套改动让出口附近的流速分布与入口附近出现了可辨识的不对称性，最大流速偏低约百分之七，这个偏差的方向与实验测量的压力降偏大是一致的。

网格质量对传热系数预测的影响，比对阻力系数预测更为敏感。电池模组周围的流场含有多个尺度不一的涡结构：冷却板流道内的层流-过渡流区、模组棱角后的分离泡、出口扩张段的射流掺混区。这些区域的传热系数空间变化剧烈，需要足够的网格分辨率来捕捉。项目中在冷却板流道内布置了至少十五层边界层网格，增长率一点二，确保y+值在一点五以内以激活低雷诺数传热修正。模组周围的自由流体区域采用了局部加密，网格尺寸从主流区的四毫米收紧到模组表面的二毫米。这套网格策略让努塞尔数（Nu）沿流动方向的分布与Dittus-Boelter关联式的预测趋势一致，偏差在百分之八以内——这个比例在工程中通常可接受，但如果目标是电芯级的温度一致性优化，这个偏差仍需进一步压缩。

材料物性的温度依赖性，在这个项目中从第二个循环开始就不能再忽略。电芯在内阻发热的作用下，温度上升导致电解液电导率变化，反过来影响发热功率。这个电-热双向耦合在Fluent中无法直接实现，需要借助UDF（用户自定义函数）将发热功率定义为温度的函数。项目中采用了Bernardi等人提出的电化学产热模型，产热率由可逆熵热项和不可逆焦耳热项组成，其中熵热系数随温度的变化关系引用了实验拟合曲线。UDF每十迭代步更新一次产热率，同时限制温度变化率以避免收敛困难。引入这个耦合后，稳态解的温度分布出现了明显变化：先前均匀发热假设下温度最高的电芯中心区域，在实际耦合模型中因为局部温度高导致内阻增大、产热更多，形成了正反馈，最高温度比均匀发热假设高出约二点五摄氏度。这个差异说明，忽略产热率的温度依赖性，会系统性低估热失控风险。

边界条件的热环境设置，需要与实际运行场景逐一对应。项目中冷却板入口采用了恒温边界（三百摄氏度），这个设置在理论上简化了入口温度控制逻辑，但实际系统中冷却液入口温度沿流动方向是逐渐上升的。更准确的设置是质量流量入口加工质温度入口，让冷却液温度在求解过程中根据能量守恒自行演化。出口边界采用了压力出口，静压设为零表压，同时开启回流温度设定，防止出口出现非物理的温度回流。环境边界（电池包外壳）采用了对流换热边界，换热系数取自然对流的典型值约十瓦每平方米每开尔文，这个数值在电池包密闭环境中偏于保守，实际换热可能被高估。将这个换热系数下调到五，外壳温度上升了约三摄氏度，整个温度场的分布更为均匀——这个调整方向，与实验测量中包内温度梯度偏小的观察是一致的。

这个项目最终交付的仿真模型，在三个放电倍率（一C、一点五C、二C）下与实验数据的温差预测偏差均不超过三摄氏度，电芯温差（最大减最小）的预测误差在零点八摄氏度以内。回过头看，达到这个精度水平的关键，不在于使用了多高级的湍流模型或多细的网格，而在于逐一核查了流固热耦合链上的每个环节：界面网格对齐、湍流普朗特数设定、求解顺序选择、物性温度依赖性、边界条件物理合理性。流体热仿真的误差来源是链条式的，任何一个环节的疏忽，都会让整个温度预测偏离真实——理解这条链条，并在每个环节留下可追溯的设置记录，才是可交付的科研仿真工作。