努塞尔数(Nusselt Number)——对流换热的增益倍数

努塞尔数(Nusselt number,简称 Nu)量化对流换热强度相对于纯导热的倍数关系,是传热学中最核心的无量纲参数。工程上真正关心的往往是”换热量到底有多少”,而 Nu 把答案归一成了一个数:Nu = 3 意味着对流把传热增强了 3 倍。换热器设计、芯片冷却方案评估、CFD 换热结果校验,最终几乎都落到这个数的比较上。

定义与表达式

$$
Nu = \frac{\text{对流换热强度}}{\text{纯导热强度}} = \frac{h L}{k}
$$

其中:

$h$:对流换热系数(W/(m²·K)),反映流体与固体表面的实际换热能力;
$L$:特征长度(如管道直径、平板长度,单位:m);
$k$:流体热导率(W/(m·K)),反映流体自身的导热能力。

物理意义:

Nu = 1:对流换热强度等于纯导热(仅存在于静止流体或对称热边界条件);

Nu > 1:对流增强传热(如流动流体冲刷壁面);

Nu ≫ 1:强湍流或高流速下显著强化换热。

物理意义:边界层与传热机制

边界层主导作用
- Nu 直接关联速度边界层与热边界层的相对发展:
  - 高 Nu:热边界层薄,温度梯度大(如湍流冲击冷却);
  - 低 Nu:热边界层厚,传热效率低(如层流自然对流)。
- 典型流体举例:

流体类型	Nu 范围	边界层特征
空气(强制对流)	10–100	热边界层薄,梯度显著
水(自然对流)	5–50	热边界层较厚
液态金属(高速流)	100–1000	极薄热边界层,高效传热

与普朗特数(Pr)的关联
- 对于强制对流,Nu 常通过 Re(雷诺数) 和 Pr(普朗特数) 关联:
  $$Nu = C \cdot Re^m \cdot Pr^n$$
  - $C, m, n$:经验系数(依赖几何与流动状态);
  - Pr 的影响:高 Pr 流体(如油)Nu 较低,低 Pr 流体(如液态金属)Nu 较高。

关键影响因素

流动状态
- 层流:$Nu \propto Re^{0.5}$(如管内层流 Nu ≈ 3.66);
- 湍流:$Nu \propto Re^{0.8}$(如光滑管湍流 $Nu \approx 0.023 Re^{0.8} Pr^{0.4}$)。
几何形状
- 平板:Nu 随长度增加而减小(入口效应);
- 圆管:Nu 随曲率与长径比变化(如弯曲管 Nu 比直管高 20–30%)。
热边界条件
- 恒壁温 vs 恒热流:后者通常给出更高 Nu(因温度梯度更陡)。

在 CFD 中的应用

壁面函数选择
- 高 Nu 场景需细网格(y⁺ ≈ 1)以解析热边界层;低 Nu 可用标准壁面函数(如 nutkWallFunction)。
数值稳定性约束
- 网格 Peclet 数需满足 Peᵢ = Reᵢ · Pr ≤ 2,否则需改用迎风格式避免振荡。
空化与相变模拟
- 在 interPhaseChangeFoam 中,Nu 影响气-液界面蒸发/冷凝速率。

工程应用场景

领域	典型问题	Nu 的作用
换热器设计	管壳式换热器优化	通过 Nu 评估翅片效率与压降平衡点
电子散热	芯片冷却液选型	高 Nu 流体(如水)优先用于高热流密度场景
能源工程	核反应堆液态金属冷却	低 Pr 流体 Nu ≈ 5–10,适合高效传热设计
化工反应器	反应物温度均匀性控制	通过 Nu 预测混合死区与热点风险

典型案例:

微通道散热器:通过微尺度涡流发生器将 Nu 提升至 200+,解决 5G 芯片散热瓶颈;

太阳能集热管:选择性涂层 + 真空层设计使 Nu 降低至 0.5,抑制热损失。

与其他无量纲数的关系

Nu 与 Bi(毕渥数):
- $Bi = hL/k_s$($k_s$ 为固体热导率),Bi 用于固体内部导热阻力分析,Nu 用于流体侧对流分析;
Nu 与 St(斯坦顿数):
- $St = Nu / (Re \cdot Pr) = h / (\rho c_p U)$,表征传热与流动动能的比率。

小结

努塞尔数 Nu 是连接流动特性与传热效率的核心桥梁:

量化传热增益:指导换热设备选型(如高 Nu 流体优先用于紧凑式换热器);
约束数值设计:确保 CFD 壁面解析精度(y⁺ 与 Nu 强关联);
优化热管理方案:从微电子冷却到反应堆安全均依赖 Nu 范围分析。

CFD 算完一个换热问题,先把截面 Nu 与经验关联式对一遍,是最基本的可信度检查。

本文整理自个人 CFD 学习笔记系列