努谢尔特数(Nusselt Number)详解

在热对流研究中,为了定量描述流体运动对热传递强度的影响,工程师和科学家通常不直接使用具有单位的物理量,而是使用无量纲数(Dimensionless Numbers)。在众多无量纲数中,努谢尔特数(Nusselt Number, $Nu$) 是衡量对流换热能力最核心的参数。它直接反映了对流换热与纯导热之间的相对强度。
努谢尔特数定义为流体边界层内的对流换热系数与流体导热系数的比值,并乘以特征长度。其数学表达式为:

$$Nu = \frac{h L}{k}$$

其中各参数的物理意义如下:

  • $h$ (Convective Heat Transfer Coefficient):对流换热系数,单位为 $\text{W}/(\text{m}^2 \cdot \text{K})$。它代表了流体与固体表面之间热量传递的效率。
  • $L$ (Characteristic Length):特征长度,单位为 $\text{m}$。特征长度的选择取决于几何形状。例如,在管内流动中,通常取管径 $D$;在平板流动中,通常取平板的长度 $x$。
  • $k$ (Thermal Conductivity):流体的导热系数,单位为 $\text{W}/(\text{m} \cdot \text{K})$。它描述了流体本身通过分子热运动传递热量的能力。

物理意义详解

从公式的结构来看,$Nu$ 实际上是一个比值,它比较了“对流换热”与“纯导热”的效率。

  1. 当 $Nu = 1$ 时:
    这意味着 $h = k/L$。在这种情况下,流体层的热传递完全通过导热完成,流体的运动(对流)对热传递没有起到任何增强作用。这通常发生在流体处于完全静止的状态时。

  2. 当 $Nu > 1$ 时:
    这意味着对流换热的效果优于纯导热。$Nu$ 的数值越大,说明流体的运动(如湍流、扰动等)对热传递的增强作用越显著。在工业换热器设计中,我们的目标通常是通过增加流速或改变几何结构来提高 $Nu$ 值,从而获得更高的换热效率。

通过引入 $Nu$,我们可以将复杂的对流换热问题转化为一个寻找经验关联式的过程,从而避开直接求解复杂的纳维-斯托克斯方程(Navier-Stokes equations)和能量方程。

努谢尔特数与其他无量纲数的关系

在实际的对流换热分析中,$Nu$ 很少孤立存在。它通常是流体流动状态(由雷诺数 $Re$ 表征)和流体物性(由普朗特数 $Pr$ 表征)的函数。这种关系可以用通用的经验关联式表示为:

$$Nu = C \cdot Re^m \cdot Pr^n$$

其中 $C, m, n$ 是由几何形状、流动状态(层流或湍流)以及边界条件决定的常数。

1. 雷诺数 ($Re$) 的作用

雷诺数 $Re = \frac{\rho v L}{\mu}$ 描述了惯性力与粘性力的比值。它决定了流动的性质(层流或湍流)。当 $Re$ 增大时,流体运动更加剧烈,边界层变薄,从而显著提升 $Nu$ 值。

2. 普朗特数 ($Pr$) 的作用

普朗特数 $Pr = \frac{\nu}{\alpha} = \frac{\mu c_p}{k}$ 描述了动量扩散率(粘性)与热扩散率(导热)的比值。它反映了流体物性的特征。$Pr$ 决定了速度边界层与温度边界层的相对厚度。

因此,通过 $Re$ 和 $Pr$ 的组合,我们可以建立起一套完整的对流换热预测体系。

常见应用场景与关联式示例

在工程实践中,针对不同的几何构型,研究人员已经总结出了大量的 $Nu$ 关联式。

平板边界层流动

对于在平板上流过的流体,在层流状态下,努谢尔特数可以通过以下形式近似:
$$Nu_x \approx 0.332 \cdot Re_x^{1/2} \cdot Pr^{1/3}$$

圆管内流动

对于管内流动,情况会根据雷诺数的变化而不同。在完全发展的湍流状态下,常用的关联式为:
$$Nu_D = 0.023 \cdot Re_D^{0.8} \cdot Pr^{1/3} \quad (\text{Dittus-Boelter 公式})$$
该公式广泛应用于化工和热力学设备的换热计算中。

实例计算

问题描述:
假设有一根直径 $D = 0.05,\text{m}$ 的圆管,内部流经水。已知水的物性参数如下:

  • 导热系数 $k = 0.6,\text{W}/(\text{m} \cdot \text{K})$
  • 动力粘度 $\mu = 0.001,\text{Pa} \cdot \text{s}$
  • 密度 $\rho = 1000,\text{kg}/\text{m}^3$
  • 运动粘度 $\nu = 1 \times 10^{-6},\text{m}^2/\text{s}$
  • 比热容 $c_p = 4180,\text{J}/(\text{kg} \cdot \text{K})$
  • 流速 $v = 2,\text{m}/\text{s}$

请计算该管内流动的对流换热系数 $h$。

计算步骤:

  1. 计算雷诺数 ($Re$):
    $$Re = \frac{v D}{\nu} = \frac{2 \times 0.05}{1 \times 10^{-6}} = 100,000$$
    由于 $Re > 4000$,判定为湍流流动。

  2. 计算普朗特数 ($Pr$):
    $$Pr = \frac{\nu}{\alpha} = \frac{\nu \cdot \rho \cdot c_p}{k} = \frac{1 \times 10^{-6} \times 1000 \times 4180}{0.6} \approx 6.97$$

  3. 利用 Dittus-Boelter 公式计算 $Nu$:
    $$Nu = 0.023 \cdot (100,000)^{0.8} \cdot (6.97)^{1/3}$$
    $$Nu \approx 0.023 \cdot 10,000 \cdot 1.91 \approx 439.3$$

  4. 求解对流换热系数 $h$:
    根据 $Nu = \frac{h D}{k}$,得:
    $$h = \frac{Nu \cdot k}{D} = \frac{439.3 \times 0.6}{0.05} = 5271.6,\text{W}/(\text{m}^2 \cdot \text{K})$$

结论:
通过计算得出,该工况下的对流换热系数约为 $5271.6,\text{W}/(\text{m}^2 \cdot \text{K})$。

总结

努谢尔特数是连接微观流体物性与宏观换热性能的桥梁。通过理解 $Nu$ 的定义及其与 $Re$、$Pr$ 的内在联系,工程师可以有效地预测复杂流动环境下的热传递强度。在设计冷却系统、热交换器或分析工业反应器时,准确选取合适的 $Nu$ 关联式是确保热管理系统可靠性的关键步骤。