理查德森数(Richardson number,简称 Ri)是流体力学中描述流动稳定性的关键无量纲参数,用于衡量浮力效应与剪切效应(惯性力)的相对重要性。大气边界层会不会翻腾混合、海洋分层流的湍流会不会被抑制、自然对流散热设计会不会失稳,这些问题第一眼该看的往往就是 Ri。它与瑞利数 Ra 名字相近,但物理含义不同,后者专注浮力驱动对流的强度,前者关心浮力与剪切的竞争。

基本定义与物理意义

Ri 的数学表达式为:

$$
Ri = \frac{\text{浮力效应}}{\text{剪切效应}} = \frac{g \beta \Delta T H}{U^2}
$$

其中:

  • $g$:重力加速度;
  • $\beta$:流体的体积膨胀系数;
  • $\Delta T$:特征温差(如流体层上下壁面温差);
  • $H$:特征长度(如流体层高度);
  • $U$:特征流速(如剪切流速)。

物理意义:

  • Ri > 0:浮力主导,抑制湍流,流动趋于稳定(如分层大气中热空气上升受阻);
  • Ri < 0:剪切力主导,促进湍流,流动失稳(如强风破坏温度分层);
  • Ri = 0:浮力效应消失,流动退化为纯剪切驱动(如理想槽道流)。

计算公式与相关参数

Ri 与雷诺数($Re$)、普朗特数($Pr$)、瑞利数($Ra$)的关系为:

$$
Ri = \frac{Ra}{Re^2 \cdot Pr}
$$

其中:

  • $Ra = g \beta \Delta T H^3 / (\nu \alpha)$:表征浮力驱动的对流强度;
  • $Re = U H / \nu$:表征惯性力与粘性力的比值;
  • $Pr = \nu / \alpha$:动量扩散与热扩散的比率。
参数 符号 物理含义
$Ra$ 瑞利数 浮力驱动的对流强度
$Re$ 雷诺数 惯性力与粘性力的比值
$Pr$ 普朗特数 动量扩散与热扩散的比率

物理意义深化

Ri 值对流动状态的影响:

  • Ri > 1:层流占优,热量传递以热传导为主(如平静大气层);
  • 0 < Ri < 1:过渡状态,可能发生间歇性湍流;
  • Ri < 0:强剪切导致湍流增强(如风速较大时的地表边界层)。

临界值意义:当 Ri = 0 时,系统退化为无重力影响的纯剪切流动(如理想管道流),此时能量方程中无量纲温度可能奇异,需谨慎处理。

应用场景与判据

  1. 气象学与海洋学:
    • 分析大气边界层稳定性(如逆温层中 Ri > 0 抑制污染物扩散);
    • 预测海洋分层流的混合效率(Ri > 0.25 时湍流被抑制)。
  2. 工程热力学:
    • 设计自然对流散热系统(如电子设备散热),需控制 Ri 以避免过热或流动失稳。
  3. 流体稳定性判据:
    • Ri 是判断开尔文-亥姆霍兹不稳定性(Kelvin-Helmholtz instability)的关键参数,当 Ri < 0.25 时易出现涡旋脱落。

变化因素与影响因素

主要变量:

  • 温差 $\Delta T$:温差增大 → Ra 增大 → Ri 增大,流动更稳定;
  • 流速 $U$:流速增大 → Re 增大 → Ri 减小,流动易失稳。

其他因素:

  • 流体性质(如 $Pr$ 高的油类 Ri 更易为正);
  • 几何尺度(如 $H$ 增大强化浮力效应)。

与相关概念的区别

  • Ri vs. Ra:
    • $Ra$ 仅表征浮力强度,而 Ri 综合了浮力与剪切力的竞争;
    • 例如:$Ra = 0$ 时 Ri 必为 0,但 Ri = 0 不要求 $Ra = 0$(如 $g=0$ 但 $\Delta T \neq 0$ 时仍可传导传热)。
  • Ri vs. 弗劳德数($Fr$):
    • $Fr = U / \sqrt{gH}$ 描述惯性力与重力的比值,而 $Ri = 1/Fr^2$,两者本质关联。

小结

理查德森数 Ri 是量化流动稳定性的核心参数,其值决定浮力与剪切力的平衡状态:

  • Ri > 0 → 稳定分层;
  • Ri < 0 → 湍流增强;
  • Ri = 0 → 纯剪切流动。

实际应用中需结合 $Ra$、$Re$、$Pr$ 等参数综合分析,尤其在气象预测、热工设备设计及环境流体力学中不可或缺。

工程备忘

在两相流模拟中,如果两相之间的气泡速度都很小、又采用了 PBM(群体平衡模型),就不能再用流速、雷诺数来衡量流动状态是否为湍流;PBM 与层流模型也是互斥的。(2025-07-16 备注记录,日后专文展开)

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