3. 热对流与热辐射：牛顿冷却定律、自然对流与强制对流、热辐射基本定律
在SiP多芯片封装中，芯片产生的热量除了通过固体传导（导热）路径传递外，还需要通过封装表面与周围流体（空气或冷却液）之间的对流换热，以及通过电磁波形式的热辐射散发到环境中。本章将系统阐述这两种基本换热模式的核心定律、工程分类及其在SiP热管理中的应用。
3.1 牛顿冷却定律：对流换热的基石
对流换热的计算基础是牛顿冷却定律（Newton's Law of Cooling），它描述了固体表面与周围流体之间的热交换速率。该定律的数学表达式为：
Q = h · A · (T_s - T_f)其中：
- Q —— 对流换热量（W）
- h —— 对流换热系数（W/(m²·K)），表征换热强度
- A —— 有效换热面积（m²）
- T_s —— 固体表面温度（K 或 °C）
- T_f —— 流体主流温度（K 或 °C）
对流换热系数h并非物性常数，而是受流体性质、流动状态、表面几何形状及温度差等多种因素影响的综合参数。下表给出了典型场景下h值的数量级范围：
| 换热场景 | h 典型范围 (W/(m²·K)) | SiP工程应用 |
|---|---|---|
| 空气自然对流 | 5 ~ 25 | 无风扇便携设备、低功耗IoT模块 |
| 空气强制对流 | 25 ~ 250 | 带散热风扇的通信基站、服务器SiP |
| 液体强制对流（水冷） | 500 ~ 15,000 | 高功率密度SiP、数据中心液冷模块 |
| 相变换热（沸腾/冷凝） | 2,500 ~ 100,000 | 均温板、热管与SiP集成散热方案 |
3.2 自然对流与强制对流
根据流体运动的驱动力来源，对流换热可分为自然对流和强制对流两大类。在SiP多芯片布局优化中，区分这两种模式对于选择散热策略和设定仿真边界条件至关重要。
3.2.1 自然对流
自然对流由流体内部密度差引起，密度差通常源于温度梯度。流体受热膨胀上升，冷流体下沉补充，形成自发的循环流动。其换热强度主要取决于格拉晓夫数（Gr）和普朗特数（Pr）。
- 特点：无运动部件、零噪音、高可靠性，但换热能力有限。
- SiP应用：智能手机、可穿戴设备、无风扇工业控制器等对噪音和可靠性要求高的场景。
- 仿真关键：需考虑封装朝向（水平/垂直）、周围空间尺寸对流动的约束效应。
3.2.2 强制对流
强制对流通过外部动力（风扇、泵、风机）迫使流体流经发热表面，流速远高于自然对流，从而显著提升换热系数。其强度由雷诺数（Re）和努塞尔数（Nu）关联描述。
- 特点：换热效率高、可控性强，但引入功耗、噪音和可靠性风险。
- SiP应用：高性能计算SiP、射频功率放大器模块、车载计算平台。
- 仿真关键：需准确设定入口流速/流量、湍流模型（如k-ε模型）及流道几何细节。
下表总结了两种对流模式在SiP热设计中的关键差异：
| 对比维度 | 自然对流 | 强制对流 |
|---|---|---|
| 驱动力 | 浮升力（密度差） | 外部机械动力（风扇/泵） |
| 典型h值 (W/(m²·K)) | 5 ~ 25 | 25 ~ 250 |
| 流速范围 (m/s) | 0.1 ~ 0.5 | 1.0 ~ 10.0 |
| 功耗与噪音 | 无 | 有（风扇功耗、风噪） |
| 可靠性影响 | 高（无运动部件） | 中（风扇寿命受限） |
| SiP典型场景 | 手机、穿戴、IoT传感器 | 服务器、基站、ADAS域控 |
3.2.3 自然对流与强制对流判定准则
工程中常用格拉晓夫数与雷诺数平方的比值（Gr/Re²）来判定对流主导模式：
- Gr/Re² ≪ 1：强制对流主导
- Gr/Re² ≈ 1：混合对流（需同时考虑）
- Gr/Re² ≫ 1：自然对流主导
3.3 热辐射基本定律
热辐射是物体因自身温度而向外发射电磁波的能量传递方式。与导热和对流不同，热辐射不需要介质，在真空中也能进行。在SiP封装中，当芯片温度较高（>80°C）或自然对流受限时，辐射散热的占比不可忽略。
3.3.1 黑体与辐射力
黑体是理想化的辐射体，能够吸收所有波长的入射辐射，其辐射力由斯特藩-玻尔兹曼定律（Stefan-Boltzmann Law）描述：
E_b = σ · T⁴其中：
- E_b —— 黑体辐射力（W/m²）
- σ —— 斯特藩-玻尔兹曼常数，5.67×10⁻⁸ W/(m²·K⁴)
- T —— 黑体温度（K）
3.3.2 实际物体辐射与发射率
实际物体的辐射力小于同温度下的黑体，引入发射率（ε）修正：
E = ε · σ · T⁴发射率ε取决于材料种类、表面粗糙度、温度及辐射波长，取值范围为0~1。常见SiP相关材料的发射率：
| 材料 | 表面状态 | 发射率 ε（典型值） |
|---|---|---|
| 硅芯片 | 抛光表面 | 0.15 ~ 0.25 |
| 硅芯片 | 未抛光/钝化层 | 0.60 ~ 0.80 |
| 铜散热盖 | 抛光 | 0.03 ~ 0.05 |
| 铜散热盖 | 氧化/发黑处理 | 0.70 ~ 0.85 |
| PCB FR4 | 绿色阻焊 | 0.85 ~ 0.95 |
| 铝散热器 | 阳极氧化 | 0.75 ~ 0.85 |
3.3.3 辐射换热网络与角系数
两个或多个表面之间的净辐射换热量计算涉及角系数（View Factor, F₁₂），表示表面1发射的辐射能到达表面2的比例。对于两个灰体表面之间的辐射换热：
Q_rad = (σ · (T₁⁴ - T₂⁴)) / ((1-ε₁)/(ε₁·A₁) + 1/(A₁·F₁₂) + (1-ε₂)/(ε₂·A₂))在SiP封装内部，芯片、基板、外壳之间构成复杂的辐射换热网络，需通过仿真软件（如FloEFD、Icepak）自动计算角系数矩阵。
3.4 对流与辐射的耦合效应
在SiP实际工作环境中，对流和辐射总是同时存在并相互影响。总散热量为两者之和：
Q_total = Q_conv + Q_rad = h · A · (T_s - T_f) + ε · σ · A · (T_s⁴ - T_surr⁴)其中T_surr为周围环境辐射温度。在低风速自然对流场景下，辐射散热量可能占总散热量的30%~50%，不可忽略。
3.5 本章小结
- 牛顿冷却定律是对流换热计算的基础，h值的准确选取是SiP热仿真的关键。
- 自然对流依靠浮升力驱动，适合低功耗、高可靠性场景；强制对流依靠外部动力，适合高功耗场景。
- 热辐射遵循斯特藩-玻尔兹曼定律，发射率ε是影响辐射换热的核心参数。
- 在SiP多芯片布局优化中，需综合考虑对流与辐射的耦合效应，尤其在自然对流和高温工况下，辐射散热占比显著。
- 通过提高表面发射率（如发黑处理）和优化流道布局，可有效增强封装整体的散热能力。