4. 热仿真链路与误差来源

热仿真并不是“把模型丢进软件里算一下”这么简单。从 CAD 几何到最终温度云图，中间要经过几何简化、材料赋值、边界条件设定、网格划分、求解器配置与收敛判定六个关键环节。每一个环节都会向结果中注入误差，而这些误差并非彼此独立——它们会沿着仿真链路逐级传递、叠加甚至放大。理解这条链路，比记住某个软件按钮的位置重要得多。

一、仿真链路：误差是如何逐级传递的

可以把热仿真看作一条“信息加工流水线”：几何提供形状信息，材料提供热物性信息，边界条件提供与外界交换热量的信息，网格把连续方程离散化，求解器负责迭代求解，收敛判据决定何时停止。上游的任何近似，都会成为下游无法消除的偏差。

二、几何简化：最容易被低估的误差源

CAD 模型往往包含大量对热行为影响极小的细节：倒角、螺纹、小孔、装饰性凹槽。直接划分这些特征会导致网格量爆炸，因此工程上通常做简化。但简化必须“有依据地删”，而不是“方便地删”。

以某功率模块散热器为例：原始模型有 48 片鳍片，每片根部有 R0.5 的圆角。若保留圆角，网格量约 1200 万；去掉圆角后约 380 万。圆角对总散热量的影响通常小于 1%，因为圆角处增加的换热面积占比很小。但若把鳍片厚度从 1.2 mm 简化为 1.0 mm，热阻会明显上升——这属于“改变了主导热路径”，不能简化。

判断原则是：简化对象的热阻与系统总热阻之比小于 5% 时，可考虑简化；大于 10% 时必须保留或做等效处理。薄壁、窄缝、导热垫等往往承担关键热路径，即使尺寸小也不能随意删除。

三、材料参数：非线性与各向异性带来的偏差

材料参数误差常被简化为“查手册填个数”，但实际影响很大。导热系数 k 随温度变化，铝合金在 25°C 时约 237 W/(m·K)，到 150°C 时降到约 220 W/(m·K)，降幅约 7%。如果仿真中固定用 25°C 的值，高温区温度会被低估。

更隐蔽的是界面材料。导热硅脂的等效导热系数与厚度、压力、老化程度强相关，标称 3 W/(m·K) 的硅脂在 0.1 mm 厚度下实际可能只有 1.5～2.5 W/(m·K)。建议在仿真中把界面层单独建为薄层实体，而不是用“理想接触”代替。

各向异性材料（如 PCB、石墨片、碳纤维）必须输入方向相关的 k 值。PCB 面内导热系数可达 20～60 W/(m·K)，而厚度方向仅 0.3～0.8 W/(m·K)，差两个数量级。若按各向同性处理，热扩散路径会完全错误。

四、边界条件：换热系数与辐射的取值陷阱

自然对流换热系数 h 通常在 3～10 W/(m²·K)，强制风冷 20～100 W/(m²·K)，水冷可达 1000～10000 W/(m²·K)。取值偏差 20%，在自然对流主导的场景下可能导致结温偏差 5～10°C。

辐射在高温或低风速场景不可忽略。一个 80°C 的表面，辐射换热约 5～8 W/(m²·K) 等效，与自然对流同量级。若仿真中只设对流不设辐射，表面散热被低估，温度偏高。

边界条件还应与物理场景一致：风冷散热器用“对流系数”还是“风扇 P-Q 曲线 + 流场求解”，结果差异可能超过 15%。前者快但粗糙，后者慢但更接近真实。

五、网格与求解器：离散误差和迭代误差

网格质量直接决定离散误差。温度梯度大的区域（热源附近、界面层、边界层）必须加密。判断网格是否足够，不能只看“网格数量”，而要做网格无关性验证：将网格量增加 1.5～2 倍，若关键监测点温度变化小于 1°C 或 2%，可认为网格收敛。

求解器设置中，稳态问题常用 SIMPLE 或 Coupled 算法；瞬态问题的时间步应满足 CFL 条件，且每步内迭代残差下降至少两个数量级。松弛因子过大易发散，过小则收敛慢。收敛判据不能只看残差，还要监测关键点温度、热平衡误差等物理量是否稳定。

六、误差控制检查表（示例）