4. 热设计仿真模型构建

热设计仿真模型不是把几何丢进求解器就能得到可信结果。真正决定仿真价值的，是能否把物理机制拆成可参数化的模块，并让每个模块的输入、边界和耦合关系都能被追溯、被扫描、被优化。本节围绕导热、对流、辐射、流固耦合四条主线，说明如何把工程对象转化为参数化模型，并给出可直接参考的配置示例。

一、先定义参照系与符号约定

在建模前必须固定坐标系与热流方向约定，否则后续参数扫描会出现符号混乱。本节统一采用：热流从高温指向低温为正；固体内部以傅里叶定律 q = -k∇T 描述，其中 k 为导热系数（W/(m·K)）；表面对流以牛顿冷却定律 q = h(T_s - T_f) 描述，h 为对流换热系数（W/(m²·K)），T_s 为表面温度，T_f 为流体参考温度；辐射以斯特藩-玻尔兹曼定律 q = εσ(T_s⁴ - T_sur⁴) 描述，ε 为发射率，σ = 5.67×10⁻⁸ W/(m²·K⁴)。所有图示与公式均遵循这一约定。

二、导热模块的参数化

导热参数化的核心是把材料属性、几何尺寸和热源强度变成可扫描变量。以一个 40 mm × 30 mm × 2 mm 的铝基板为例，其导热系数 k 可取 180～220 W/(m·K) 作为扫描区间，芯片热源功率 P 取 2～8 W，基板厚度 t 取 1.0～3.0 mm。稳态下的一维热阻可写为：

R_cond = t / (k · A)
A = 0.04 m × 0.03 m = 1.2e-3 m²
若 t = 2 mm, k = 200 W/(m·K)
R_cond = 0.002 / (200 × 1.2e-3) = 8.33e-3 K/W

这个 8.33 mK/W 的量级说明：对于小面积、高导热材料，导热热阻往往不是瓶颈，真正的温升来自扩散热阻和界面热阻。因此在参数化时，除了体导热，还应把界面材料（TIM）的等效热阻 R_TIM 作为独立参数，典型值 0.1～1.0 K·cm²/W，换算到 1.2e-3 m² 面积上约为 0.83～8.3 K/W，量级远大于体热阻。这一对比直接决定了优化重点应放在界面而非基板厚度。

三、对流边界参数化

对流边界的关键参数是 h 和参考温度 T_f。自然对流下 h 通常为 5～15 W/(m²·K)，强制风冷可达 25～100 W/(m²·K)，液冷可超过 1000 W/(m²·K)。参数化时不应只给一个常数，而应建立 h 与风速 v 的关系。对于平板强制对流，可用简化关联式：

h ≈ 10.45 - v + 10·sqrt(v)   (v 单位 m/s, h 单位 W/(m²·K))
v = 2 m/s  → h ≈ 10.45 - 2 + 14.14 ≈ 22.6 W/(m²·K)
v = 5 m/s  → h ≈ 10.45 - 5 + 22.36 ≈ 27.8 W/(m²·K)

可以看到，风速从 2 m/s 提升到 5 m/s，h 仅增加约 23%，而风扇功耗约按 v³ 增长。因此在多目标优化中，风速不能只作为提升换热的自由变量，必须同时纳入功耗与噪声目标。参数化设置时应把 v 作为设计变量，把 h(v) 作为派生参数，而不是直接扫描 h，否则优化结果无法映射回可执行的风扇选型。

四、辐射边界参数化

辐射在自然对流或真空环境中不可忽略。参数包括发射率 ε、环境温度 T_sur 和视角因子 F。对于两个表面间的净辐射，可用：

q_rad = εσF(T_s⁴ - T_sur⁴)
若 T_s = 350 K, T_sur = 300 K, ε = 0.85, F = 1
q_rad = 0.85 × 5.67e-8 × (350⁴ - 300⁴)
     ≈ 0.85 × 5.67e-8 × (1.50e10 - 8.10e9)
     ≈ 333 W/m²

333 W/m² 的辐射热流在自然对流 h = 10 W/(m²·K)、温差 50 K 时对应对流热流 500 W/m²，两者处于同一量级。因此当表面温度超过约 80 ℃ 且环境为自然对流时，辐射必须显式建模。参数化时，ε 可按表面处理取值：阳极氧化铝 0.80～0.88，抛光铝 0.05～0.10，喷涂黑漆 0.90～0.95。若忽略 ε 的差异，仿真可能高估或低估温升 10～20 K。

五、流固耦合的参数化接口

流固耦合（CHT）不是把流体和固体分别算完再拼起来，而是在交界面上同时满足热流连续和温度连续：

q_solid→interface = q_interface→fluid
T_solid,interface = T_fluid,interface

参数化设置时，交界面应单独定义为一个命名边界（如 interface_wall），并指定耦合类型为共轭传热。流体侧需要设置入口速度、入口温度、出口压力，固体侧需要设置热源和外部对流/辐射边界。一个常见的错误是把流体侧的对流系数人为指定为常数，这会破坏耦合的物理一致性。正确做法是让求解器从流场和温度场中计算局部 h，仅在需要降阶时才用等效 h 替代。

六、边界条件参数化配置示例

下表给出一个可直接参考的参数化边界配置，适用于上述铝基板加芯片的示例。所有数值均为教学示例，不代表实测结果。

七、从参数到优化目标的映射

完成上述参数化后，模型输出不应只是温度云图，而应提取可量化的目标：最高结温 T_j,max、平均表面温度 T_avg、热阻 R_total = (T_j,max - T_f)/P、压降 Δp 和风扇功耗 P_fan。这些量分别对应热性能、流动性能和能耗，构成多目标优化的目标向量。只有把边界条件参数化到这一层，后续的灵敏度分析和 Pareto 前沿才有意义。

上图把四个物理模块与边界条件层、优化目标层串联起来，强调参数化不是单点设置，而是从材料、几何、边界到目标的完整映射链。实际建模时，建议先冻结几何与材料，再逐层激活对流、辐射和耦合，每激活一层都检查能量守恒残差是否收敛到 1% 以内，这样才能保证后续多目标优化的结果可信。