4. 仿真工具与流程

功率电子液冷系统的仿真，本质上是把“冷却液带走热量的物理过程”翻译成计算机能求解的离散方程。这个翻译过程由一条固定的技术链路完成：几何清理 → 网格划分 → 求解器设置 → 迭代求解 → 后处理判读。工具只是这条链路上不同环节的载体，理解链路本身比记住软件名字更重要。

一、主流工具的分工与选型

液冷仿真涉及两类物理场：流体流动与传热（CFD），以及固体内部的导热（热仿真）。工程中常见的组合方式有三类。

• 通用CFD平台：ANSYS Fluent、Ansys CFX、STAR-CCM+、OpenFOAM。它们能同时求解流场与共轭传热（CHT），适合冷板流道、歧管分配、多支路并联等复杂流动问题。STAR-CCM+在表面包面与多面体网格上效率突出；Fluent的UDF与多相模型生态成熟；OpenFOAM开源、可脚本化批量计算，适合参数扫描。
• 电子散热专用工具：Ansys Icepak、Flotherm、6SigmaET。它们内置风扇曲线库、PCB等效导热、器件封装热阻模型，建模速度快，适合机柜级、板级风冷与液冷冷板耦合分析，但复杂流道细节能力弱于通用CFD。
• 系统级一维工具：Simcenter Amesim、GT-SUITE、Modelica类平台。用于整车或整机热管理回路，把冷板、泵、换热器、管路抽象为集中参数元件，计算快，用于控制策略与工况匹配，不解析局部流场。

选型判断可以归结为三个问题：是否需要看清流道内部的流速与温度分布？是否需要与电路、控制联合仿真？计算资源与周期是否允许？若三者答案偏向“局部细节”，选通用CFD；偏向“快速迭代与系统匹配”，选一维工具；介于两者之间，选电子散热专用工具。

二、几何清理：把CAD变成可网格化的“干净水体”

从CAD导出的模型通常带有大量对流动无影响的细节：倒角、螺纹、小孔、装配间隙、重复面。这些特征会迫使网格在这些位置加密，导致网格量暴涨而结果几乎不变。几何清理的目标是得到一个封闭、无干涉、无碎面的流体域与固体域。

典型操作包括：删除小于0.5 mm的倒角与圆角；封闭工艺孔与螺钉孔；把装配体中的微小间隙合并为实体接触；抽取冷板内部流道形成流体域；对进出口延长3～5倍水力直径，避免回流影响。清理后必须做水密性检查：流体域应是单一连通体，无自由边、无重复面。若存在泄漏面，求解时会报“反向体积”或直接发散。

三、网格划分：质量决定结果可信度

液冷仿真推荐使用多面体网格或六面体核心+边界层。多面体网格在相同精度下单元数约为四面体的1/3～1/2，且对复杂流道适应好。边界层是传热计算的关键：壁面附近速度梯度大，必须布置足够层数。

边界层设置的经验规则：第一层高度由目标y+反推。对壁面函数法，y+取30～300；对低雷诺数解析，y+≈1。以水为例，流速1 m/s、水力直径8 mm时，若目标y+≈1，第一层高度约0.01～0.02 mm，增长率1.2，层数8～12层，总边界层厚度约0.3 mm。网格质量指标需满足：正交质量>0.15，长宽比在流道内<100，无负体积。

网格无关性验证不可省略：将网格量按1.5倍递增，观察压降与最高温度变化。当压降变化小于2%、温度变化小于0.5 ℃时，可认为网格收敛。这一步是判断结果是否可信的直接依据。

四、求解器选择与设置要点

液冷问题通常采用压力基耦合求解器，稳态或瞬态取决于工况。稳态用于额定工况下的压降与温度分布；瞬态用于启停、负载阶跃、流量突变等过程。湍流模型方面，冷板流道内多为过渡区或湍流，SST k-ω 在近壁面与主流区都有较好表现，是液冷仿真的常用选择；若流速很低（Re<2300），可切换为层流模型，避免湍流模型过度耗散。

边界条件设置需与实物一致：入口给质量流量或体积流量，出口给压力出口（0 Pa表压），壁面给无滑移，固体域与流体域交界面自动生成共轭传热耦合面。材料属性中，冷却液的密度、比热、导热系数、粘度应随温度变化，尤其乙二醇水溶液在低温段粘度变化显著，用常数会带来明显偏差。

收敛判断不能只看残差。残差降到1e-4以下只是必要条件，还需监控出口温度、最高温度、压降三个物理量是否稳定。若残差下降但Tmax仍在漂移，说明尚未收敛。能量守恒检查是最后一道关：进出口焓差与总热耗的相对误差应小于1%。

五、后处理：从云图到工程结论

后处理不是“看图”，而是提取判据。液冷系统最关心的四个量是：

速度云图用于识别死区与回流区；温度云图用于定位热点；流线图用于判断是否存在短路流。若发现某支路流量明显偏低，应回到几何清理阶段检查该支路是否被误删或网格过密导致阻力偏大，而不是直接修改求解器参数。

把这条链路走通一次，后续的冷板改型、流量调整、工况切换都只是在这条链路上替换输入。工具会更新，链路不会变。