2. 液冷电子系统热源建模

液冷系统的散热能力再强，如果热源模型本身失真，后续的流道设计、流量分配和结温预测都会失去意义。热源建模的任务，是把芯片、器件和封装层面的发热行为，转换成CFD可用的体积热源、面热源或等效热阻边界。本节围绕四个相互衔接的问题展开：功耗如何映射到几何体、热流密度如何在空间上分布、器件内部热量如何经热阻网络传导、以及接触界面如何用参数化方式表达。

一、芯片功耗映射：从数据手册到CFD体热源

芯片功耗通常以总功耗P_total给出，但CFD需要的是单位体积的发热率q'''（W/m³）。映射的第一步是确定发热体积。对于裸die，发热体积取有源区体积；对于带封装盖的器件，若只关心封装外壳到冷板的热路，可将功耗施加在die对应的等效块体上，而不是整个封装体。映射公式为：

q''' = P_total / V_active
V_active = L_die × W_die × t_active

例如某AI加速芯片，数据手册给出典型功耗700 W，die尺寸26 mm × 32 mm，有源层等效厚度0.15 mm，则V_active≈1.248×10^-7 m³，q'''≈5.61×10⁹ W/m³。这个量级说明芯片内部热流密度极高，必须配合高导热界面和高效冷板才能控制结温。

实际建模中还要区分稳态与瞬态。稳态热源可直接用平均功耗；瞬态工况下，若功耗随时间变化，应把P(t)以表格或表达式形式施加到体热源上，并注意时间步长要能分辨功耗变化的主要频率成分。对于多芯片模块，每个die单独建体热源，避免用总功耗除以总体积造成局部热流密度被平均化。

二、热流密度分布：均匀假设的边界与修正

把die当作均匀体热源是最常见的简化，但真实芯片的发热并不均匀。CPU的核、GPU的SM阵列、HBM的I/O区，热流密度可能相差数倍。若只做封装级冷却设计，均匀化通常可接受；若目标是预测局部热点或评估微通道冷板，就必须引入非均匀热流密度分布。

工程上常用两种修正方式。第一种是分区法：把die划分为若干功能区块，每块给定不同的q'''，例如计算核区取1.5倍平均值，I/O区取0.6倍平均值。第二种是热点叠加法：在均匀背景热源上叠加一个或多个高斯分布热点，热点峰值和半宽由芯片热仿真或红外热像数据标定。分区法适合结构规则、功能分区明确的芯片；热点叠加法适合热点位置随机或需要参数化扫描的场景。

无论采用哪种方式，都要保证总功耗守恒：各分区功耗之和等于P_total。这是检查热源映射是否正确的第一道关卡。

三、器件热阻网络：从结到冷板的串联路径

芯片热量从结区到冷板，依次经过die、TIM1、封装盖、TIM2、冷板壁面，最后进入冷却液。每一段都可以用热阻表示，形成串联热阻网络。热阻定义为单位热流下的温升：R = ΔT / P，单位K/W。对于一维导热，R = t / (kA)，其中t为厚度，k为导热系数，A为传热面积。

以典型封装为例：die厚度0.15 mm、k=120 W/(m·K)、面积832 mm²，R_die≈0.0015 K/W；TIM1厚度0.05 mm、k=3 W/(m·K)，R_TIM1≈0.020 K/W；封装盖厚度1.0 mm、k=390 W/(m·K)，R_lid≈0.0031 K/W；TIM2厚度0.1 mm、k=5 W/(m·K)，R_TIM2≈0.024 K/W。可见TIM1和TIM2合计约0.044 K/W，占总热阻的绝大部分。这意味着界面材料的选择和装配压力控制，往往比提高冷板导热系数更有效。

在CFD中，热阻网络有两种用法。一是作为边界条件：把总热阻折算成等效换热系数施加在封装表面；二是作为验证工具：用CFD算出的冷板表面温度和热阻网络反推结温，与实测或芯片热仿真对比。后者更适合需要分离各层贡献的场合。

四、接触热阻与界面材料参数化

两个固体表面即使看起来平整，微观上也只有少量凸峰接触，其余空隙被空气填充。热量通过接触点时收缩扩散，形成接触热阻R_c。接触热阻不是材料常数，而是压力、表面粗糙度、平面度、温度和界面材料共同作用的结果。在CFD中，接触热阻通常以薄层热阻或等效导热系数的方式参数化。

对于导热硅脂、相变材料、导热垫等界面材料，常用参数包括：厚度t_TIM、本体导热系数k_TIM、接触热阻R_c1和R_c2（两侧各一个）。总界面热阻为：

R_interface = R_c1 + t_TIM / (k_TIM × A) + R_c2

参数化时，建议把R_c1和R_c2合并为有效接触热阻R_c,eff，并通过压力-热阻曲线或供应商提供的压缩曲线来标定。例如某导热垫在30 psi下厚度从1.0 mm压缩到0.82 mm，k_TIM=3.5 W/(m·K)，面积832 mm²，则本体热阻约0.028 K/W；若R_c,eff=0.015 K/W，总界面热阻约0.043 K/W。压力不足时，厚度和接触热阻同时增大，总热阻可能翻倍。

在CFD中实现接触热阻有三种常见方式：一是把TIM建成薄实体层，赋予k_TIM和压缩后厚度；二是在两个固体界面之间插入薄壳热阻，直接给定R_c,eff；三是用等效换热系数h_c=1/(R_c,eff·A)施加在界面上。第一种最直观但网格量最大，第二种适合参数扫描，第三种适合快速估算。选择时需权衡精度、网格成本和参数不确定性。

最后需要强调，热源建模不是一次性工作。芯片功耗随工作负载变化，界面材料随时间和温度老化，装配压力在批量生产中也有波动。因此，参数化建模的价值在于：当这些输入发生变化时，可以快速重新计算结温和热流分布，而不是重新搭建整个模型。把功耗、热流密度分布、热阻网络和接触热阻都写成可调参数，是液冷电子系统CFD设计走向工程化的关键一步。