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

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

实时计算分层流中浮力 vs 惯性力比 Ri = gβΔT·L/V²。Brunt-Väisälä 振动数 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,显示「中性混合」,继续降低则变为「稳定分层」。设定 ΔT 为负值时,温度分层反转(上方暖),Ri 变负属于对流不稳定分类。

计算结果
理查德森数
稳定性
混合临界速度
Brunt-Väisälä 振动数
竖直分层与水平流

用彩色梯度表示竖直温度分层(下=暖色,上=冷色,当 Δ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-Väisälä 振动数:稳定分层内竖直位移流体粒子的浮力恢复固有振动频率(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 = 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 滑块能清晰展示这两个临界值的效应变化。RANS 湍流模型中 Ri 也是调节湍流粘性的重要参数。
🙋
Brunt-Väisälä 振动数 0.129 rad/s 是什么振动?水里看不出来振动啊…
🎓
好问题!N = √(gβΔT/L) 是「稳定分层中竖直位移的流体粒子在浮力作用下回振」的固有频率(角频率,rad/s)。不是要你手工晃动,而是内部重力波的最大频率。本工具默认 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 分)。
🙋
混合临界速度 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。

常见问题

理查德森数 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 为「稳定分层」。
Brunt-Väisälä 振动数 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 ≈ 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⁻³(水)算出 V_crit ≈ 0.314 m/s,要用大孔径低速扩散设计保层。超速会乱流混合,高温水变平均温,效率大幅下降。

常见误解与注意事项

最常见的误解是「Ri=0.25 时混合瞬间开始」。实际上 Ri 是线性稳定论的必要条件,边界附近(0.1~0.5)混合呈间断性,没有明确阈值。观测显示 Ri=0.2~0.3 范围内出现「湍流间歇性」(间断性湍流爆发),平均混合强度高度依赖 Ri。设计时应安全考虑 Ri < 0.5 为「可能混合」。

常见误解二:「Ri 用层平均值就能代表」。本工具用「体 Richardson 数」(层厚平均),实际流体在竖直方向上温度、速度梯度变化,局部「梯度 Richardson 数」Ri_g = N²/(dU/dz)² 更精准预测混合。CFD 在每个网格点算 Ri_g 判断层流→乱流,即使体 Ri > 1,薄剪切层的局部 Ri 若 < 0.25 也会混合。

误解三:「Ri > 1 时完全静止」。稳定分层内部仍有内部重力波(周期 2π/N)传播,远距离(大气数千 km、海洋数百 km)输送能量,最终在远处砕波产生混合。山岳波励起的内部波上升到平流层砕波,对行星级运动量输送重要。N 值是分层允许的最大波频,低于 N 的内部波均可存在。

使用指南

  1. 用密度梯度滑块(slDT)设置分层流的密度差。海洋典型值 0.5~2.0 kg/m³,湖泊 0.1~0.5 kg/m³。
  2. 输入竖直距离(slL),一般观测层厚 1~100 m,大气边界层 50~200 m。
  3. 设置风速或流速(slV),定量化剪切层的惯性作用。河川混合 0.1~1.5 m/s,沿岸流 0.2~0.8 m/s。
  4. 用盐度梯度系数(slBeta)调整浮力效应,同时计算 Ri、N、混合临界速度。

具体计算例

北太平洋沿岸案例:密度梯度 dρ/dz=0.8 kg/m³、层厚 L=30 m、水平流速 V=0.4 m/s、盐度梯度系数 β=0.78×10⁻³,输入得 Brunt-Väisälä 振动数 N≈0.0085 rad/s、理查德森数 Ri≈0.45。Ri<1/4 处剪切不稳定区,混合临界速度 V_mix≈0.51 m/s,超过此值层破坏、密度界面消失。

实务注意事项