4. 热仿真控制方程与边界条件
热管与均温板系统的热仿真，本质上是在求解一个带相变传热、接触界面和多种换热路径的稳态或瞬态温度场。无论使用 Icepak、FloTHERM、Ansys Mechanical 还是自编有限元程序，控制方程和边界条件的设置方式决定了仿真结果是否可信。本节从方程形式出发，说明导热、对流、辐射三类边界在软件中的落地方法，并给出接触热阻、热流密度、对流换热系数和环境温度的具体设置思路。
一、控制方程：能量守恒的统一形式
稳态、无内热源、各向同性材料的三维导热方程为：
∂/∂x(k·∂T/∂x) + ∂/∂y(k·∂T/∂y) + ∂/∂z(k·∂T/∂z) = 0当存在内热源 q'''（如芯片体热源、热管内部等效热源）时，右端改为 -q'''。瞬态问题还需加上 ρc∂T/∂t 项。热管与均温板内部存在工质相变与毛细回流，工程上通常不直接求解两相流，而是把热管等效为高导热率的固体域，或使用热阻网络模型替代。因此仿真中真正需要仔细设置的是边界条件，而不是把热管内部流动也纳入求解。
二、三类边界条件的物理含义与设置方法
第一类边界条件（Dirichlet）直接给定温度：T = T_w。常用于已知热源表面温度、恒温冷板或实验测得的参考点温度。设置时要注意，若把热源表面直接设为恒温，会掩盖芯片内部导热热阻，通常只在有实测数据时才这样处理。
第二类边界条件（Neumann）给定热流密度：-k·∂T/∂n = q''。芯片功耗除以有效发热面积即可得到 q''。例如 15 W 芯片、发热面积 10 mm × 10 mm，则 q'' = 15 / (1e-4) = 1.5e5 W/m²。设置时需明确热流方向：默认沿表面外法线流出为正，若软件采用流入为正，则需取负号。热流密度应施加在真实发热区域，而不是整个封装表面，否则会低估局部温升。
第三类边界条件（Robin）给定对流换热：-k·∂T/∂n = h(T_w - T_∞)。这里 h 是对流换热系数，T_∞ 是环境温度。自然对流下 h 通常为 5～15 W/(m²·K)，强制风冷为 25～100 W/(m²·K)，水冷可到 1000～10000 W/(m²·K)。设置时 h 和 T_∞ 必须成对出现，且参考温度要与实验或工况一致。
辐射边界通常以第四类方式叠加：q''_rad = εσ(T_w⁴ - T_sur⁴)。ε 为发射率，σ = 5.67e-8 W/(m²·K⁴)，T_sur 为环境辐射温度。自然对流散热器中辐射占比可达 20%～40%，不可忽略；强制风冷时辐射占比通常较小，但高温工况仍需保留。
三、接触热阻：最容易被低估的界面参数
均温板与芯片、均温板与散热器之间不可能理想接触，微观粗糙峰之间只有部分真实接触，其余空隙被空气或导热界面材料填充。接触热阻 R_c 定义为界面温差除以通过的热流密度：R_c = ΔT / q''。单位常用 K·cm²/W 或 K·m²/W。
典型取值参考：干接触（金属对金属）约 1～5 K·cm²/W；导热硅脂约 0.1～0.5 K·cm²/W；导热垫片约 0.5～2 K·cm²/W；焊接或烧结界面可低至 0.05 K·cm²/W 以下。设置时，若软件支持薄层单元，可把 R_c 折算为等效厚度和等效导热率；若只支持界面热阻，则直接输入 R_c 值。需要注意，R_c 与压力、表面粗糙度、温度和老化程度有关，仿真中应使用与装配工艺一致的数值，而不是随手取一个“看起来合理”的数。
四、对流换热系数与环境温度的设置方法
h 的获取有三条路径：一是查经验关联式，如平板自然对流 Nu = 0.54Ra^(1/4)，再由 h = Nu·k/L 计算；二是用 CFD 先算流场再提取局部 h；三是直接引用实验或供应商数据。对于风冷散热器，若已知风量和散热器热阻，也可反推等效 h。
环境温度 T_∞ 必须与 h 的参考温度一致。自然对流时 T_∞ 取远离散热器的环境空气温度；强制风冷时取进风口温度；水冷时取进水温度。若软件中同时存在多个温度参考（如辐射 T_sur 和对流 T_∞），要分别设置，不能混用。例如机箱内自然对流，T_∞ 可取 45 ℃，T_sur 可取机箱内壁温度 50 ℃，两者相差不大但物理含义不同。
五、一个可复现的设置示例
以 15 W 芯片、10 mm × 10 mm 发热面、均温板 40 mm × 40 mm、散热器 50 mm × 50 mm 为例。芯片表面热流密度 q'' = 15 / 1e-4 = 1.5e5 W/m²。芯片与均温板之间涂导热硅脂，R_c1 = 0.3 K·cm²/W；均温板与散热器之间用导热垫片，R_c2 = 1.0 K·cm²/W。散热器自然对流 h = 8 W/(m²·K)，T_∞ = 25 ℃，ε = 0.85，T_sur = 25 ℃。
若忽略辐射，散热器表面温度可由 q'' = h(T_s - T_∞) 估算：散热器有效散热面积按 0.05 m² 计，q'' = 15 / 0.05 = 300 W/m²，则 T_s ≈ 25 + 300/8 = 62.5 ℃。加入辐射后，等效 h_rad = εσ(T_s² + T_sur²)(T_s + T_sur)，在 60 ℃ 附近约为 6 W/(m²·K)，总等效 h 约 14 W/(m²·K)，T_s 降至约 46 ℃。这说明自然对流下辐射不可忽略。实际仿真中，这些数值会随温度变化迭代收敛，但设置逻辑不变：先确定热流密度，再逐段设置接触热阻，最后成对设置 h 与 T_∞、ε 与 T_sur。
完成设置后，检查三项：总热流是否等于芯片功耗；界面温差是否与 R_c × q'' 一致；散热器表面温度是否落在 h 和 T_∞ 决定的合理范围内。若三者对不上，优先检查单位制、面积和参考温度，而不是急于调整材料导热率。