理查森数模拟器 返回
环境流体力学模拟器

理查森数模拟器 — 分层流的稳定性

实时计算分层流中浮力与惯性之比 Ri = gβΔT·L/V²。可视化 Brunt-Vaisala 振荡频率 N、混合极限速度 V_crit 以及 Ri=0.25/1 边界的三区稳定性 (不稳定/中性/稳定)。

参数设置
温度差 ΔT
K
特征长度 L
m
流速 V
m/s
热膨胀系数 β
×10⁻³ /K

默认值 (ΔT=5 K、L=10 m、V=5 m/s、β=3.4×10⁻³ /K 为空气) 下,Ri ≈ 0.067,判定为「混合发生」,V_crit ≈ 2.58 m/s,N ≈ 0.129 rad/s。将 V 调低至约 2 m/s 时 Ri 超过 0.25 进入「中性混合」,继续调低至 1.6 m/s 以下进入「稳定分层」。ΔT 为负时层结反转 (上暖下冷),Ri 为非正时归入对流不稳定。

计算结果
理查森数
稳定性
混合极限速度
Brunt-Vaisala 频率
垂直分层与水平流动

用颜色渐变表示垂直温度分层 (下=暖色、上=冷色,当 ΔT>0)。水平方向的蓝色箭头表示流速 V。Ri 越大水平层越能维持 (稳定分层),Ri<0.25 时 Kelvin-Helmholtz 不稳定产生波状混合,界面起伏增大。中央波形即为界面的瞬时起伏。

Ri-稳定性曲线 (对数轴)

横轴:Ri (log10,0.01〜100) / 纵轴:稳定性区分 / 红带:Ri<0.25 不稳定 (混合发生) / 黄带:0.25≤Ri≤1 中性混合 / 绿带:Ri>1 稳定分层 / 红色竖线:Ri=0.25 与 Ri=1 边界 / 黄色标记:当前 Ri 工作点。拖动 V 滑块可见标记跨越各带。

理论与主要公式

理查森数:分层流中浮力 (由温差引起的密度差) 与惯性力之比的无量纲数。

$$\mathrm{Ri} = \frac{g\,\beta\,\Delta T\,L}{V^2}$$

Brunt-Vaisala 振荡频率:稳定分层内垂直位移的流体微团在浮力恢复下的固有角频率 (rad/s):

$$N = \sqrt{\frac{g\,\beta\,\Delta T}{L}}$$

混合极限速度:使 Ri = 0.25 的流速,超过即可能发生 Kelvin-Helmholtz 混合:

$$V_{\mathrm{crit}} = \sqrt{\frac{g\,\beta\,\Delta T\,L}{0.25}} = 2\sqrt{g\,\beta\,\Delta T\,L}$$

$g$ 是重力加速度 (9.81 m/s²),$\beta$ 是热膨胀系数 (理想气体取 $1/T$),$\Delta T$ 是层的上下温差,$L$ 是特征长度 (层厚),$V$ 是穿过层的流速。Ri < 0.25 为不稳定 (混合发生),0.25 ≤ Ri ≤ 1 为中性混合,Ri > 1 为稳定分层 (混合抑制)。

理查森数模拟器是什么

🙋
大气课上学过「夜晚风会停下来」,这真的能用理查森数解释吗?平时不会想到 Ri,它到底是衡量什么的?
🎓
好问题。理查森数 Ri = gβΔT·L/V² 是分层流中「浮力 vs 惯性」之比,夜间近地逆温层正是教科书例子。日落后地面被辐射冷却,地面附近变冷而高空仍较暖,形成逆温层 (上暖下冷的稳定分层)。当 ΔT 大、V 小时 Ri 超过 1,进入「稳定分层」,风 (惯性) 抗不过浮力,湍流混合停止——这就是「夜风停」的由来。本工具默认值 (ΔT=5 K、L=10 m、V=5 m/s、β=3.4×10⁻³ /K) 下 Ri ≈ 0.067,仍在「混合发生」区;将 V 调到 1 m/s 时 Ri ≈ 1.7,标签切换为「稳定分层」。
🙋
Ri=0.25 和 Ri=1 这两个边界为什么是这两个特定数值?0.25 看起来挺奇怪的。
🎓
本质性的问题。Ri=0.25 是 Miles 与 Howard 在 1961 年通过线性稳定性分析证明的必要条件:Ri 低于 0.25 时 Kelvin-Helmholtz 不稳定 (剪切层波动) 可以增长;反之 Ri ≥ 0.25 时线性理论保证稳定。Ri=1 是经验值,源自实验与观测:当 Ri 大于 1 时强稳定分层显著抑制湍流。在本工具拖动 V,会依次穿过 Ri=0.25 与 Ri=1,混合性质两阶段大幅改变。CFD 的 RANS 湍流封闭也用 Ri 修正稳定分层流的涡粘性。
🙋
Brunt-Vaisala 频率显示 0.129 rad/s,到底是什么在振荡?眼前的水并没有看到在抖动…
🎓
很有意思的点。N = √(gβΔT/L) 是稳定分层中一个被垂直位移的流体微团在浮力作用下回到平衡位置的固有角频率。宏观看不到什么在摇,N 给出的是内部重力波的最高频率上限。实际可观察的例子有:湖泊温跃层上的内部波 (湖面平静,但水面下数米的温度层在起伏)、大气山岳波 (山地下游气流中的重力波)。海洋中 N ≈ 0.001〜0.01 rad/s (周期 10〜100 分)、平流层 N ≈ 0.02 rad/s (周期 5 分)。本工具默认值下周期 T = 2π/N ≈ 48.7 秒。
🙋
混合极限速度 V_crit = 2.58 m/s,这个值怎么用?能当成什么设计参数?
🎓
非常实用的问题。V_crit 是 Ri=0.25 对应的流速,即「超过则发生混合,低于则维持分层」的临界。本工具默认值 (ΔT=5 K、L=10 m、β=3.4×10⁻³) 下 V_crit ≈ 2.58 m/s。比如数据中心冷却:冷气从机柜底部上升 10 m、上下温差 5 K 时,机柜内排热引起的空气流速若超过 2.58 m/s,精心建立的冷气分层就会被搅乱混合。同理在分层蓄热罐 (热水在上、冷水在下) 中,将循环泵流速控制在 V_crit 以下可保持分层、提升能效。试试将本工具调到 L=2 m、ΔT=20 K,看蓄热罐 V_crit。

常见问题

理查森数 Ri = gβΔT·L/V² 是分层流中浮力 (由温差引起的密度差) 与流动惯性力之比的无量纲数。g 是重力加速度,β 是热膨胀系数,ΔT 是垂直温差,L 是特征长度 (层厚),V 是流速。本工具默认值 (ΔT=5 K、L=10 m、V=5 m/s、β=3.4×10⁻³ /K、空气) 下 Ri ≈ 0.067,落入 Kelvin-Helmholtz 不稳定的混合发生区 (Ri < 0.25)。Ri 是判断大气边界层、海洋温跃层与建筑通风中能否发生湍流混合的基本参数。
Ri=0.25 是 Miles-Howard 线性稳定性定理的必要条件:当 Ri 低于 0.25 时,Kelvin-Helmholtz 不稳定 (剪切层波动→翻卷→湍流混合) 可能发生。Ri=1 是经验性的强稳定分层阈值,超过该值时湍流被大幅抑制,分层结构长时间维持。本工具将 Ri<0.25 标记为「混合发生」、0.25≤Ri≤1 标记为「中性混合」、Ri>1 标记为「稳定分层」,拖动 V 滑块可见标记在曲线上跨越各区。
Brunt-Vaisala 振荡频率 N = √(gβΔT/L) 是稳定分层内、被垂直微小位移的流体微团在浮力作用下回到平衡位置的固有角频率 (rad/s)。N 越大恢复力越强,给出内部重力波的最高频率上限。本工具默认值下 N ≈ 0.129 rad/s,周期 T = 2π/N ≈ 48.7 秒,是内部波的典型尺度。海洋中 N ≈ 0.001〜0.01 rad/s (周期 10〜100 分)、大气平流层 N ≈ 0.02 rad/s (周期 5 分) 为代表值。
Ri 在涉及分层流稳定性的广泛场合都会用到。具体例:(1) 大气边界层稳定度评价 (夜间辐射冷却→近地逆温层→Ri 增大→风停)、(2) 海洋温跃层混合 (Ri<0.25 时内部波破碎)、(3) 高层建筑温度成层自然通风、(4) 火灾安全 (烟气层稳定限)、(5) 室内空调冷气分层 (地板冷却效率)、(6) 河口盐楔 (淡水与海水界面混合)。提高本工具的 V 使 Ri 降至 0.25 以下,即可观察到混合发生的临界。

实际应用

夜间大气稳定度评价:地面 10 m 风速 V=2 m/s、地面〜50 m 高度温差 ΔT=3 K (上空更暖的逆温层) 时,空气 (β=3.4×10⁻³ /K) 下 Ri ≈ 1.25,本工具判定为「稳定分层」。这是大气污染物与花粉在近地累积、直到清晨起风才被带走的典型条件。环境标准用 Ri > 0.25 的持续时间来选定 PG 稳定度 D〜F 级。

海洋温跃层与内部波:夏季沿岸海域,表层 (20°C) 与深层 (10°C) 间形成温跃层 (厚度 L≈10 m、ΔT=10 K),海水 β≈0.21×10⁻³ /K 下 N ≈ 0.0144 rad/s,周期约 7.3 分。当潮流速度超过 V_crit = √(gβΔT·L/0.25) ≈ 0.575 m/s 时 Ri < 0.25,内部波破碎将深层冷水与营养盐带到表层,促发浮游植物增殖 (沿岸生态关键过程)。在本工具尝试 β=0.21、ΔT=10、L=10 即可再现。

高大空间 (中庭) 的分层通风:车站、机场、大型礼堂等高大中庭 (高度 L=20 m、ΔT=8 K) 常形成稳定温度分层。为不破坏地板冷气层,将气流速度 V 控制在 V_crit 以下是节能关键。本工具在 ΔT=8、L=20、β=3.4×10⁻³ 下得到 V_crit ≈ 4.62 m/s,将散流器出口流速控制在此值以下可保持分层,只冷却驻足区 (距地 1.8 m),显著节省空调能耗。

蓄热罐的分层维持:太阳能与区域供热蓄热罐 (高度 L=4 m、ΔT=30 K,热水在上) 中,进出口流速 V 控制在 V_crit 以下是热效率关键。在 ΔT=30、L=4、水 (β=0.21×10⁻³ /K) 下本工具给出 V_crit ≈ 0.314 m/s,需要低流速散流器 (大孔径、缓流) 来保持分层。一旦超出 V_crit,湍流混合将罐温均匀化、破坏高品位热——蓄热效率严重依赖分层保持。

常见误解与注意事项

最常见的误解是 「Ri 跨过 0.25 立即开始混合」。实际上 Ri=0.25 只是线性稳定性的必要条件,在阈值附近 (0.1〜0.5) 波的非线性增长产生间歇混合,并非锐变。Ri=0.2〜0.3 的现场观测呈「湍流间歇性」——湍流爆发,平均混合效率随 Ri 平滑变化。本工具为方便起见以 Ri=0.25 划分区域,但工程设计时宜将 Ri < 0.5 视为「可能混合」并留有裕度。

第二个误解是 「Ri 可用层内平均值代表局部梯度」。本工具使用层厚 L 内平均 ΔT 与 V 的「整体 Richardson 数」,实际中温度与速度沿垂直方向梯度变化,局部「梯度 Richardson 数」Ri_g = N²/(dU/dz)² 更准确地预测混合的发生。CFD 在每个网格点评估 Ri_g 来判定层流-湍流转捩。即便整体 Ri > 1,薄剪切层内局部 Ri < 0.25 也会触发局部混合。

最后一个误解是 「Ri > 1 即完全静止」。稳定分层内仍允许频率不超过 N 的内部重力波远距离 (大气数千公里、海洋数百公里) 传播,并在远处破碎并向外传输能量。被地形激发的山岳波抵达平流层破碎,对行星尺度的动量输运很重要。本工具的 N 是分层内允许的最高频率,许多更低频率的内部波同时存在。