Fluent流固耦合传热模拟：CHT分析与双向耦合的实战经验

共轭传热（Conjugate Heat Transfer，CHT）是指流体和固体之间的热量传递需要同时求解流体域的对流和固体域的导热。Fluent 天然支持 CHT——在固-流界面上自动创建 Wall/Wall-shadow 对，界面上的热边界条件设为 Coupled，求解器在每次迭代时让界面两侧的温度和热流密度自动匹配。

这个过程看似全自动，但中间有若干决策点直接决定了仿真结果的精度和可信度。

固-流界面——网格不连续是怎么回事

CHT 里最基础的配置是保证固体网格和流体网格在界面上”对上”。Fluent 的网格导入时，如果固体和流体在几何上是两个不同的 Body 但在 Workbench Meshing 里用了共享拓扑（Share Topology），界面上的网格天然是共节点的——节点一一对应，没有任何插值误差。这是 CHT 的理想网格状态。

但如果固体和流体是用不同网格工具导出的（例如固体用 ANSYS Mechanical 六面体、流体用 Fluent Meshing 多面体），界面两侧的网格拓扑不同、节点不对齐——Fluent 仍然能用非共形界面（Non-Conformal Interface）来做数据传递。但非共形界面是插值传递：一侧网格面上的温度场插值到另一侧——插值平滑了温度梯度，对于界面附近温度梯度极大的情况（例如电子芯片到 TIM 材料的界面，梯度可达 10°C/μm），插值带来的温度失真能达到 5-10°C，对于结温的评估是致命的。

结论很简单：但凡计算资源的壁厚允许，固体和流体在 CHT 界面必须用共节点网格。非共形界面只应该在几何极复杂、无法通过切割和共享拓扑对齐网格的不得已情况下使用，并且要在报告里明确标注非共形界面的位置和可能的插值误差。

界面热阻——Fluent 里没有”TIM 厚度”这个参数

工程中芯片到散热片之间往往涂有一层导热硅脂（TIM），厚度只有几十微米。直接在 Fluent 里建这个厚度的固体域，网格长宽比会极度恶化——一个 50 μm 厚的 TIM 层在 20 mm 宽的芯片下方，长宽比 400:1，四面体在做这种层时几乎必定出负体积。

物理上正确的处理方式不是建出 TIM 的实体模型，而是在固-流界面之间设置接触热阻（Thermal Contact Resistance）。Fluent 的界面边界条件里有 Contact Resistance 选项——输入一个热阻值（单位 K·m²/W），界面两侧就会有温差 ΔT = q” × R。

这里的技巧在于怎么获得一个合理的 R 值。厂商给的导热硅脂数据是体热导率 k（如 5 W/m·K），换算成接触热阻需要 R = t / k，t 是 TIM 的压实厚度。但这个”压实厚度”取决于装配压力、TIM 的流动性、表面粗糙度——比数据手册上的厚度厚 30-50% 是常态。如果直接按理论值设定，仿真温度会比实测低 5-8°C。更稳妥的做法是用叠层法（ASTM D5470）的实测等效热阻，或者从类似设计的实测芯片-散热片温差反算接触热阻，再代入仿真。

瞬态耦合的时间步长——流体和固体的”时差”

流固耦合传热如果涉及瞬态过程（如电机启动温升、电子设备热冲击测试），流体和固体的热响应时间尺度差异极大。流体侧的对流时间尺度是 L/U（特征长度/流速），可能只有毫秒级；固体侧的热扩散时间尺度是 L²/α（L²/热扩散率），可能长达数分钟。

如果按流体时间尺度来跑瞬态计算，那固体温度在这几千步里几乎没变化，浪费计算；如果按固体时间尺度来跑，那流体根本看不到瞬态流场的真实演化。

做瞬态 CHT 的一个实用策略是亚迭代时间推进：对每个物理时间步（按照固体热响应的需求设，例如 0.5 s 或 1 s），流体在这个步内的伪时间（inner iteration）里做充分多的亚迭代，让流场在每个物理时间步内达到准稳态。Fluent 的 Dual Time Stepping 正是处理这种情况的标准手段——物理时间步长为固体的热响应步长（1-10 s），伪时间步长（inner iteration 的时间推进步）由 CFL 数自动控制，在伪时间内把流场推进到收敛。

如果反过来——流体变化快、固体变化也快（如内燃机活塞-冷却油耦合）——那时间步长必须由流体决定（微秒到毫秒级），固体侧在每个时间步里跟着更新温度。这种情况下 CHT 的计算量极大——一个发动机工作循环（720°曲轴转角，约 0.04 s）在 0.1°步长下需要 7200 个时间步，湍流模型在每个时间步里需要 10-20 个亚迭代才能收敛。这是 CHT 计算中最烧资源的场景，也是必须要对”这个仿真值不值得用瞬态”做清醒判断的时刻——很多发动机热管理的分析，稳态 CHT 加上工程修正因子的估算精度已经足够用于设计决策。

从仿真到验证——热成像与热电偶的交叉比对

CHT 仿真最有说服力的验证方式是：热电偶（点测量）+ 热成像（面测量）+ 仿真（全场预测）三者交叉比对。 热电偶给出某几个关键位置的绝对温度基准（精度 ±0.5°C），热成像给出温度分布的整体形态（精度受表面发射率设置影响较大，±1-2°C），仿真给出任意位置的温度分布。

三者的差异模式能透露很多信息：如果仿真和热成像的温度分布形态一致但绝对偏差在所有点上都是 +5°C——说明可能是边界条件（如环境温度、对流换热系数的经验值）整体偏移，调整一个全局参数就能对齐。如果仿真和热成像的形态在某些局部偏离——比如固体角落仿真偏冷、热成像偏热——那说明局部的几何建模（如绝热边界条件假设错误、局部空气间隙未建模）或网格分辨率有问题。

流固耦合传热仿真的价值最终不在”一次跑完导出温度云图交差”，而在通过仿真-实验的闭环迭代，逐渐理解哪些传热路径是准确的、哪些还受限于简化假设——然后下一次建模仿真时，在之前薄弱的环节多投入精度。