热流体仿真

热流体仿真，是每一个要同时对付”流动”和”传热”这两个相互纠缠的物理场的人必须掌握的耦合本领。我刚做液冷板那会儿，先单独算流场、再单独算温度，自以为两件事都搞清楚了，结果把流场速度拿来当已知边界喂给热分析，算出的温度分布和实测差了快一倍。后来才明白，热流体仿真真正的难点在于流动和传热是同时发生的——流速决定换热，温度又改物性反作用于流动，分开算等于砍掉了反馈链。把这条链路理顺，共轭传热才真正能给你可靠的温场。

一、背景：为什么流动与传热必须一起算热流体仿真

热流体仿真（共轭传热，CHT）在同一个求解里同时解流体运动的 N-S 方程和能量方程，并跨越流固界面传递温度与热流。它决定了液冷板的换热效率、电子机箱的温升分布、换热器的总传热系数。可以说，热流体仿真这件事，直接决定了你的冷却方案在打样前能不能被证伪。很多新手沿用”先流后热”的拆分思路，或只在流体域设温度忽略固体导热，等到和实测温差巨大，才发现固体内的导热路径和界面换热被整个漏掉，只能返工。更现实的是，热流体对边界（换热系数、入口温度、物性随温度变化）极其敏感，评审最看重你的流固界面和耦合设置是否自洽。

二、核心原理：流固耦合的能量传递

共轭传热的核心是流体域和固体域共享界面上的温度连续与热流连续：流体侧通过对流把热量带到壁面，壁面把热流传给固体，固体内部靠导热把热散开。在 CFD 求解里，能量方程同时作用于流体和固体区域（固体区速度为0），界面处自动满足 T 连续和 k∇T 连续。温度还会反馈：液体粘度随温度升高而下降，流速随之改变，进一步影响换热——这是单向”流影响热、热影响流”的双向耦合。核心认知是：只要存在固体参与导热（如芯片封装、冷板基板），就必须用 CHT 而非纯流体传热；否则固体当绝热处理，热路被截断，温度必然失真。物性（粘度、导热、密度）随温度变化时要用温度相关模型，否则高温区物性偏差放大误差。

三、关键技术要点：界面、物性与求解顺序

流固界面处理是地基：几何上固体和流体网格在界面要对齐（共节点或用 interface 耦合），否则热流传不过去或产生虚假跳跃。我早年算一个冷板，流体和固体网格在界面没对齐、还用了默认零热阻，结果固体侧温度均匀得不正常，实际该有的梯度全丢了。物性方面，液体粘度随温度变（如水温升粘度降约 2%/°C 量级），不设为温度相关会系统性偏差。求解上，先算流场收敛再开能量方程是常见稳妥顺序（或可同时解），但强耦合问题需同时解以保证反馈。湍流模型与近壁处理同前，换热对边界层解析极敏感，y+ 控制不能松。还有，辐射在高温或真空场景要和流动一起解，不可忽略。

四、实操流程：从几何到温场热流体仿真

第一步，建立包含流体域和固体域的完整几何，确认界面贴合，赋予各自材料参数（含温度相关物性）；第二步，划分网格，在壁面加边界层、在梯度区加密，界面处网格对齐；第三步，设物理模型（湍流、能量方程），定边界：流体入口速度/温度、出口压力、固体外表面对流/绝热、热源用体功率；第四步，求解（先流场后能量或同时），监控出口温度、压降和关键界面热流至收敛；第五步，后处理看流场、温度云图和换热系数分布。我后来养成一个习惯：跑前先画一张”热路图”——从热源经固体到流体到环境，标出每段的导热/对流阻力，确认模型覆盖了整条链，绝不漏掉固体导热那一段。界面网格对齐我也单独检查，这是 CHT 最容易翻车的点。

五、常见问题

Q1：温度比实测低很多？检查是否漏算固体导热、界面热阻或热源功率。 Q2：流固界面温度跳变？网格未对齐或interface设置错，热流不连续所致。 Q3：先流后热可以吗？弱耦合可，但强耦合或物性随温变需同时解。 Q4：粘度随温度变要不要管？高温或大温升必须设温度相关物性，否则偏差放大。 Q5：换热系数分布不均？看边界层 y+ 和流速分布，近壁解析不足会抹平。 Q6：压降和温度互相影响？这是双向耦合，物性随温变时务必开启反馈。 Q7：辐射何时开？高温、真空或温差大必须和流动一起解辐射。

讲一个具体的坑：有次仿真一个功率模块液冷，流场和温度都”收敛”了，可基板最高温比实测低 10°C。排查发现是我在固体域外表面随手设了绝热，而实际模块外壳通过对流和环境换热，这部分散热被我砍掉，导致温度虚高误判——不对，是虚低（少了散热路径温度应更高，但实测更高说明我漏了某处热阻）……复盘最终定位：冷板和壳体之间的导热硅脂界面热阻被忽略，热从芯片到 coolant 的主路瓶颈在这层，漏掉它温度自然算低。重新加上界面热阻后吻合。从那以后我对任何”固体-固体贴合”和”固体-流体”界面都先问热阻，绝不默认完美。另一个常被问的问题：为什么不能先单独流再单独热。因为真实换热系数沿壁面分布由流场决定，你把均匀 h 喂给热分析，等于用平均量代替分布量，遇到强非均匀流（如分离、冲击射流）误差极大，CHT 才是正解。还有，热源如果是瞬态脉冲，必须瞬态 CHT，稳态会漏掉峰值温度这个最危险的量。

再说一个细节：验证策略。热流体比纯流或纯热都更需要基准比对，因为两个场都可能错。我交付必附：网格无关性、界面对齐确认、与至少一个基准（经典算例或实验）的流速和温度双对比。如果你手上的热流体仿真总对不上实测，多半是界面或物性反馈的基础没打牢，而不是程序的问题。需要可复现的热流体仿真模板和共轭传热设置建议，可以联系我们拿到现成方案，少走两年弯路。

Q8：共轭传热的网格在界面怎么处理？流体和固体网格在界面要对齐或正确建立 interface 映射，否则热流不连续、温度跳变；我习惯在界面两侧都加密，并检查法向热流是否连续，这是 CHT 最容易翻车也最易被忽略的点。 Q9：辐射能和流动一起解吗？高温、真空或温差大时必须把辐射和流动耦合同解，单独后处理加辐射会漏掉温度场对流动的反馈；辐射视角的计算成本不低，但跳过它的热流体结论在大温跨下不可信。

六、复盘总结

回过头看，热流体仿真并不神秘，难的从来不是把两个场加起来，而是把流固界面、物性反馈和完整热路这条链路想清楚。我现在的习惯是每接一个新冷却结构，先画热路图确认耦合完整再 setup，基本一次出可信温场。如果你手上的热流体总对不上实验，多半是前面界面或物性的基础没打牢，与其盲目加密网格，不如回到耦合本身重新审视。需要可复现的参数文件和流程脚本，可以联系我们拿到现成模板，少走两年弯路。热流体最易漏的是固体导热那一段，只要固体参与散热，就必须用共轭传热而非纯流体传热，否则热路被截断，温度必然失真，这是和实测对不上的头号元凶。