3. 强迫对流散热物理原理
强迫对流散热的核心，是风扇或风机驱动空气掠过PCB表面，把元件产生的热量从固体表面带走。要判断一块板在给定风量下能不能压住温升，不能只看“有没有风”，而要看空气与板面之间到底交换了多少热量。这个交换能力由对流换热系数 h 决定，而 h 又由边界层状态、层流或湍流、雷诺数与努塞尔数共同决定。理解这条链条，是后续建模和参数设置的基础。
一、牛顿冷却定律与对流换热系数
强迫对流换热的基本关系是牛顿冷却定律：
Q = h · A · (Ts − T∞)
Q ：对流换热量，W
h ：对流换热系数，W/(m²·K)
A ：有效换热面积，m²
Ts ：固体表面温度，K 或 ℃
T∞ ：来流空气温度，K 或 ℃其中 h 不是材料常数，而是随气流速度、流道几何、空气物性和表面状态变化的“工况量”。自然对流下，PCB表面的 h 通常只有 5～10 W/(m²·K)；强迫风冷时，h 可提升到 20～100 W/(m²·K)，局部高速区甚至更高。也就是说，同样一块板，从自然对流改为强迫对流，换热能力可能提高一个数量级，这正是风扇能显著压低温升的原因。
在仿真中，h 一般不由用户直接输入，而是由求解器根据流场和边界层计算得到。但理解 h 的物理含义，有助于判断结果是否合理：如果某区域风速接近零，h 会回落到自然对流水平；如果风速高且湍流强，h 会明显上升。
二、边界层：热量离开表面的必经通道
空气流过PCB表面时，紧贴表面的空气因黏性作用速度降为零，形成速度边界层；同时，表面与空气之间存在温差，形成热边界层。热量必须先穿过这层“贴壁空气”，才能被主流带走。边界层越薄，温度梯度越大，换热越强；边界层越厚，热阻越大，表面越容易积热。
边界层从板前缘开始，沿流动方向逐渐增厚。在PCB场景中，这意味着：
- 板前缘和上游元件附近边界层薄，换热好，温升相对低；
- 沿流动方向下游，边界层增厚，换热减弱，下游元件容易“吃热风”；
- 元件本体、散热器和板面凸起会扰动边界层，局部可能提前转捩为湍流，反而增强换热；
- 板间狭窄流道会限制边界层发展，使流速和压降重新分布。
因此，PCB布局不能只考虑电气连接，还要考虑气流方向：把发热量大、允许温升小的器件放在上游或风速较高区域，把耐热器件放在下游，是强迫风冷布局的基本原则。
三、层流与湍流：两种截然不同的换热状态
边界层内存在两种基本流态。层流时，空气分层有序流动，热量主要靠分子扩散穿过边界层，换热较弱；湍流时，流体出现随机脉动和涡旋，热量被剧烈掺混，边界层有效厚度减小，换热显著增强。湍流的 h 通常比同条件下层流高数倍。
在PCB风道中，常见流态判断如下：
- 低风速、板面平整、上游无扰动：多为层流，h 偏低；
- 风速较高、经过风扇出口、格栅、连接器或散热器后：容易转捩为湍流，h 升高；
- 板间狭窄通道内：若流速低、间隙小，可能保持层流；若风速高，则可能发展为湍流。
仿真中如果忽略湍流模型，只按层流计算，往往会低估高风速区的换热能力，导致预测温度偏高；反之，若在不该用湍流模型的地方强行使用，也可能高估换热。因此，流态判断是选择湍流模型和网格策略的前提。
四、雷诺数：判断流态的相似准则
雷诺数 Re 表示惯性力与黏性力之比，是判断层流与湍流的核心参数。对于PCB表面外流，常用局部雷诺数：
Re_x = (ρ · u · x) / μ = (u · x) / ν
ρ ：空气密度，kg/m³
u ：来流速度，m/s
x ：距前缘距离，m
μ ：动力黏度，Pa·s
ν ：运动黏度，m²/s以 25 ℃ 空气为例，ν ≈ 1.57×10⁻⁵ m²/s。若风速 u = 2 m/s，在距前缘 x = 0.1 m 处：
Re_x = (2 × 0.1) / (1.57×10⁻⁵) ≈ 1.27×10⁴平板外流通常认为临界雷诺数 Re_c ≈ 5×10⁵ 时发生转捩。上述算例中 Re 远低于临界值，说明该处仍为层流。若把风速提高到 10 m/s，同样位置 Re ≈ 6.4×10⁴，仍低于临界值；要到 x 接近 0.4 m 或风速更高时，才可能接近转捩。这说明普通PCB尺寸和常见风速下，板面外流往往以层流为主，但风扇出口、元件尾流和格栅后局部会提前湍流化。
对于板间通道内流，则用水力直径 D_h 计算：
Re_Dh = (ρ · u · D_h) / μ
D_h = 4A / P
A：通道截面积，m²
P：湿周，m矩形通道内流通常以 Re ≈ 2300 作为层流与湍流的分界。若板间距 5 mm、通道宽 100 mm，则 D_h ≈ 9.5 mm；当风速 3 m/s 时，Re_Dh ≈ 1.8×10³，仍属层流；风速升到 5 m/s 以上才可能进入湍流过渡区。
五、努塞尔数：把流场换算成换热系数
努塞尔数 Nu 表示对流换热与纯导热之比，是无量纲换热系数：
Nu = (h · L) / k
h：对流换热系数，W/(m²·K)
L：特征长度，m
k：空气导热系数，W/(m·K)反过来，已知 Nu 即可求 h：
h = (Nu · k) / L工程上常用经验关联式估算 Nu。例如平板层流外流：
Nu_x = 0.332 · Re_x^0.5 · Pr^(1/3)其中 Pr 为普朗特数，25 ℃ 空气约 0.71。仍以上述 u = 2 m/s、x = 0.1 m 为例，Re_x ≈ 1.27×10⁴：
Nu_x = 0.332 × (1.27×10⁴)^0.5 × 0.71^(1/3)
≈ 0.332 × 112.7 × 0.892
≈ 33.4取空气 k ≈ 0.026 W/(m·K)，则：
h = (33.4 × 0.026) / 0.1 ≈ 8.7 W/(m²·K)这个数值接近自然对流水平，说明低速下板面前缘换热并不强。若把风速提高到 10 m/s，Re_x ≈ 6.4×10⁴，Nu_x ≈ 75，h ≈ 19.5 W/(m²·K)，换热能力明显提升。若进入湍流，Nu 关联式指数变为约 0.8，h 会进一步大幅增加。
六、在PCB仿真中的实际含义
把上述关系落到仿真设置上，可以得到几条直接可用的判断：
- 风速是主导变量。 h 随 Re 的 0.5 次方（层流）或 0.8 次方（湍流）增长，提高风速比增大面积更有效，但风压和噪声代价也上升。
- 特征长度要选对。 板面外流用沿流向长度，板间通道用水力直径，元件表面用元件特征尺寸。选错会导致 Re 和 Nu 偏差数倍。
- 层流与湍流模型不能混用。 低风速板间通道优先层流；风扇出口、格栅后和元件尾流区应允许湍流模型生效。
- 下游热耦合不可忽略。 边界层增厚和上游热风使下游元件入口温度升高，仿真中应关注流线温度分布，而不是只看单点风速。
- 结果合理性可用 h 反查。 若仿真得到的表面 h 低于 5 W/(m²·K)，说明该处几乎无风；若高于 100 W/(m²·K)，应检查是否局部风速过高或网格过粗。