层流与湍流自然对流的判别标准

自然对流（Natural Convection），又称自由对流，是指流体在温度梯度引起的密度差作用下，在重力场中产生的流动现象。与强制对流不同，自然对流不需要外部机械动力（如风扇或泵），其驱动力完全源于浮力。

自然对流（Natural Convection），又称自由对流，是指流体在温度梯度引起的密度差作用下，在重力场中产生的流动现象。与强制对流不同，自然对流不需要外部机械动力（如风扇或泵），其驱动力完全源于浮力。

在典型的自然对流过程中，靠近热表面的流体被加热，密度降低并向上运动；而较冷的流体则由于密度较大而向下流动，从而形成循环流动。在这种流动过程中，流体状态可分为**层流（Laminar Flow）和湍流（Turbulent Flow）**两种基本形式。判别流体处于哪种状态对于准确计算对流换热系数（$h$）至关重要，因为湍流状态下的热量传输效率远高于层流状态。

核心判别参数：瑞利数 (Rayleigh Number)

在自然对流中，判定流态的核心无量纲参数是瑞利数（Rayleigh Number, $Ra$）。瑞利数本质上是浮力与粘性力及热扩散力之比。

1. 瑞利数的数学表达式

对于一个特征长度为 $L$ 的物体，瑞利数的定义公式如下：
$$Ra_L = \frac{g \beta (T_s - T_\infty) L^3}{\nu \alpha}$$

其中各参数的物理含义为：

  • $g$：重力加速度 ($\text{m/s}^2$)
  • $\beta$：流体的体膨胀系数 ($\text{1/K}$)，对于理想气体 $\beta = 1/T_{avg}$
  • $T_s - T_\infty$：壁面温度与环境流体温度之差 ($\Delta T$)
  • $L$：特征长度 ($\text{m}$)
  • $\nu$：流体的运动粘度 ($\text{m}^2/\text{s}$)
  • $\alpha$：流体的热扩散率 ($\text{m}^2/\text{s}$)

2. 物理意义

瑞利数可以看作是格拉晓夫数（Grashof number, $Gr$）与普朗特数（Prandtl number, $Pr$）的乘积：$Ra = Gr \cdot Pr$。

  • $Gr$ 代表浮力与粘性力的比值。
  • $Pr$ 代表动量扩散率与热扩散率的比值。
    当 $Ra$ 较小时，粘性力占据主导，流体呈层状流动；当 $Ra$ 增大到一定程度，浮力引起的扰动无法被粘性力抑制，流动便转变为湍流。

不同几何形状的判别标准

流态的转变并非发生在同一个 $Ra$ 值上，而是取决于物体的几何形状和边界条件。

1. 垂直平板

对于垂直平板，流体沿壁面上升或下降。其临界瑞利数 $Ra_{cr}$ 通常被认为在 $10^9$ 左右。

  • 层流区：$Ra_L < 10^9$。流动平稳，速度分布呈抛物线状。
  • 过渡区：$10^9 < Ra_L < 10^{10}$。流动开始出现不稳定性。
  • 湍流区：$Ra_L > 10^{10}$。流动剧烈，出现大量涡流，换热能力显著增强。

2. 水平圆柱或球体

对于水平放置的圆柱，由于流动在圆周方向上的不对称性，临界值通常低于垂直平板。通常在 $Ra_L \approx 10^7 \sim 10^8$ 范围内开始发生转变。

3. 封闭腔体

在封闭空间（如双层玻璃窗）内，自然对流受限于壁面约束。在这种情况下，临界瑞利数通常较低，且流动模式会经历从纯导热 $\rightarrow$ 稳定环流 $\rightarrow$ 湍流的演变。

层流与湍流的特性能对比

为了在工程计算中正确选择努塞尔数（Nusselt number, $Nu$）的相关式，必须明确两种流态的差异：

特性 层流自然对流 湍流自然对流
流动形态 流体呈层状平行流动，无宏观混合 存在随机的涡流，剧烈的径向混合
速度分布 速度梯度平缓，边界层较厚 速度分布较平坦，边界层极薄
换热效率 较低，主要依赖分子扩散 较高，由涡流输运主导
$Nu$ 相关式 $Nu \propto Ra^{1/4}$ $Nu \propto Ra^{1/3}$
稳定性 对微小扰动不敏感 对扰动极其敏感

实际计算示例

场景描述：
一块高度 $L = 0.5\text{m}$ 的垂直铝板，表面温度 $T_s = 80^\circ\text{C}$，置于温度 $T_\infty = 20^\circ\text{C}$ 的静止空气中。请判定该自然对流过程是层流还是湍流。

已知参数（空气在平均温度 $50^\circ\text{C}$ 时）：

  • $\nu \approx 1.8 \times 10^{-5} \text{m}^2/\text{s}$
  • $\alpha \approx 2.6 \times 10^{-5} \text{m}^2/\text{s}$
  • $\beta \approx 1/323 \text{K}^{-1} \approx 3.1 \times 10^{-3} \text{K}^{-1}$
  • $g = 9.81 \text{m/s}^2$
  • $\Delta T = 80 - 20 = 60\text{K}$

计算步骤：

  1. 代入瑞利数公式：
    $$Ra_L = \frac{9.81 \times 3.1 \times 10^{-3} \times 60 \times (0.5)^3}{1.8 \times 10^{-5} \times 2.6 \times 10^{-5}}$$
  2. 计算数值：
    $$Ra_L \approx \frac{0.228}{4.68 \times 10^{-10}} \approx 4.87 \times 10^8$$
  3. 判定：
    由于 $Ra_L \approx 4.87 \times 10^8 < 10^9$，该流动状态判定为层流。

结论：
在这种情况下，工程师应选用针对层流设计的 $Nu$ 相关式来计算换热系数，若错误地使用了湍流公式，将会导致严重高估散热能力，从而引发设备过热风险。