2. 自然对流物理基础
自然对流散热不靠风扇，靠的是流体自身因温度差产生的密度差。理解这条“从热量到流动、再从流动到散热”的因果链，是后续所有散热设计计算的起点。本节把这条链拆成四个环节：浮升力从哪来、边界层如何长起来、流动何时从层流变湍流、以及如何用瑞利数和努塞尔数把这一切量化。
一、浮升力与密度差：自然对流的驱动力
当固体表面温度高于周围空气时，贴近表面的空气被加热、密度下降。以竖直热壁为例，设壁温为 Tw，远离壁面的环境温度为 T∞，则温差 ΔT = Tw − T∞。空气密度随温度近似满足理想气体关系 ρ ≈ p/(R T)，温度升高则密度降低。单位体积流体受到的净浮升力可写为：
F_b = (ρ∞ − ρ) · g ≈ ρ · g · β · ΔT
其中 β 为体积膨胀系数，理想气体 β ≈ 1/T_f（T_f 为膜温，单位 K）
g = 9.81 m/s²这个净力使热空气沿壁面向上加速，冷空气从下方补充，形成持续循环。需要注意两点：一是浮升力正比于 ΔT，温差越大驱动力越强；二是它同时正比于重力加速度 g，因此在微重力环境下自然对流几乎消失，这也是航天电子必须依赖导热和辐射的原因。此外，密度差不仅由温差产生，湿度、组分差异也会影响，但在电子散热场景中通常只考虑温度引起的密度差。
二、边界层发展：从壁面到主流
浮升力并非作用于整个空间，而是集中在紧贴壁面的一层薄流体中，这就是热边界层（同时存在速度边界层）。沿竖直热壁自下而上，边界层经历三个阶段：
- 起始段：壁面下缘边界层极薄，速度为零，热量主要靠导热传递，速度剖面尚未充分发展。
- 发展段：浮升力持续加速流体，边界层厚度 δ 随高度 x 增大，速度剖面出现峰值，卷吸周围冷空气。
- 充分发展/转捩段：边界层增厚到一定程度后，惯性力与黏性力之比增大，扰动被放大，流动开始向湍流过渡。
边界层厚度是判断散热能力的关键：δ 越薄，壁面附近的温度梯度越大，对流换热越强。因此竖直散热器通常不宜过高——超过一定高度后边界层过厚，上部区域的换热效率反而下降。工程上常把竖直平板高度控制在 150～300 mm 量级，或采用横向开槽打断边界层来“刷新”换热。
三、层流与湍流自然对流判据
自然对流同样存在层流与湍流两种流态，判据是瑞利数 Ra。对竖直平板，常用经验临界值：
| 流态 | 判据（竖直平板） | 换热特征 |
|---|---|---|
| 层流 | Ra < 10⁹ | Nu ∝ Ra1/4，换热随高度增长较慢 |
| 过渡 | 10⁹ ≤ Ra ≤ 10¹⁰ | 流动不稳定，换热波动 |
| 湍流 | Ra > 10¹⁰ | Nu ∝ Ra1/3，换热随高度增长更快 |
层流时边界层内热量主要靠导热穿过，湍流时涡流混合大幅增强换热。对电子散热而言，多数小型散热器处于层流自然对流范围（Ra 约 10⁶～10⁸），而大型机柜、大功率模块表面可能进入湍流区。判断流态的意义在于：层流区增大散热面积有效，湍流区则更依赖表面粗糙度和结构扰动。
四、瑞利数与努塞尔数：把自然对流算出来
瑞利数 Ra 是自然对流的“综合相似准数”，它把浮升力、黏性力和热扩散三者打包在一起：
Ra = Gr · Pr = (g · β · ΔT · L³ / ν²) · Pr
Gr：格拉晓夫数，表征浮升力与黏性力之比
Pr：普朗特数，空气约 0.71
ν：运动黏度，m²/s
L：特征长度，竖直平板取高度，水平板取面积/周长
β：体积膨胀系数，空气取 1/T_f努塞尔数 Nu 则把对流换热系数 h 无量纲化：Nu = h L / k，其中 k 为空气导热系数。工程上通过经验关联式由 Ra 求 Nu，再反求 h。以竖直平板层流为例：
Nu = 0.59 · Ra^(1/4) （10⁴ < Ra < 10⁹，层流）
Nu = 0.10 · Ra^(1/3) （Ra > 10⁹，湍流）
h = Nu · k / L
Q = h · A · ΔT下面这张图把本节四个环节串成一条完整的物理链：从温差产生密度差，到浮升力驱动边界层，再到流态判据与瑞利数、努塞尔数的计算路径，最后落到换热系数与散热量。
五、工程判断要点
把上述物理量落到实际散热设计上，可以形成几条直接可用的判断：
- 温差是根本：ΔT 每增大一倍，浮升力近似翻倍，Ra 随之增大，换热系数提升，但结温也会上升，需要在热阻与温升之间权衡。
- 特征长度要选对：竖直平板取高度，水平板取面积除以周长，水平圆柱取直径。选错 L 会让 Ra 偏差一个数量级，直接导致 h 计算错误。
- 物性取膜温：所有空气物性（ν、k、Pr、β）都应在膜温 Tf = (Tw + T∞)/2 下取值，否则误差可达 10% 以上。
- 层流区优先增大面积：层流时 Nu ∝ Ra1/4，换热系数对尺寸不敏感，增大散热面积是最有效的方向；湍流区则更应关注表面扰动和气流组织。