通过可视化理解 Reynolds 数与流动转捩

什么是 Reynolds 数#

理解流体流动时最重要的无量纲数之一是 Reynolds 数(Re)。它表示惯性力与粘性力之比。

$Re = \frac{\rho U L}{\mu} = \frac{U L}{\nu}$

其中:

• $\rho$: 流体密度 $[\text{kg/m}^3]$
• $U$: 特征速度 $[\text{m/s}]$
• $L$: 特征长度 $[\text{m}]$
• $\mu$: 动力粘度 $[\text{Pa·s}]$
• $\nu = \mu/\rho$: 运动粘度 $[\text{m}^2/\text{s}]$

流动状态划分#

流动转捩的物理机制#

从层流到湍流的转捩始于开尔文-亥姆霍兹(Kelvin-Helmholtz)不稳定性。在存在速度梯度的边界上,小扰动放大形成涡(vortex),进而引发能量级串。

速度剖面的变化由 Navier-Stokes 方程描述:

$\frac{\partial \mathbf{u}}{\partial t} + (\mathbf{u} \cdot \nabla)\mathbf{u} = -\frac{1}{\rho}\nabla p + \nu \nabla^2 \mathbf{u}$

左边的非线性项 $(\mathbf{u} \cdot \nabla)\mathbf{u}$ 承担惯性,右边的 $\nu \nabla^2 \mathbf{u}$ 承担粘性衰减。Reynolds 数表示这两项的相对大小。

速度场可视化:Reynolds 数下的流动变化#

通过下方模拟直接查看 Reynolds 数变化对速度场的影响:

观察要点:

• 低 Re (Re ≈ 1~50): 矢量平滑而规则 — 粘性立即抑制扰动
• 中 Re (Re ≈ 100~500): 障碍物尾流出现非对称性 — 惯性力与粘性力相当
• 高 Re (Re > 1000): 矢量方向变得不规则,涡结构发育

圆柱绕流流线:卡门涡街#

圆柱(cylinder)绕流中 Reynolds 数增加时,会出现周期性涡脱落,即卡门涡街(Kármán vortex street)。

脱落频率用 Strouhal 数 $St$ 归一化:

$St = \frac{f D}{U_\infty}$

其中 $f$ 为涡脱落频率,$D$ 为圆柱直径,$U_\infty$ 为来流速度。圆柱情况下 $St \approx 0.2$($100 < Re < 10^5$)几乎保持恒定。

通过下方模拟查看圆柱尾流的流线模式:

观察要点:

• 圆柱前缘**驻点(stagnation point)**处流线分离
• 尾流区域对称性破缺,出现周期性涡脱落
• 提高 speed(增加流速)会加剧尾流不稳定性

数值方法中的 Reynolds 数:网格需求#

湍流模拟(DNS, Direct Numerical Simulation)需要解析最小尺度,即 Kolmogorov 微尺度:

$\eta = \left(\frac{\nu^3}{\varepsilon}\right)^{1/4}$

总网格数随 Reynolds 数的标度关系:

$N \sim Re^{9/4}$

即 $Re$ 增加 10 倍,网格数需要增加约 178 倍。这是实际湍流仿真依赖 RANS、LES 等湍流模型的原因。

总结#

• Reynolds 数是惯性力/粘性力之比,决定流动特性的核心无量纲数
• 超过临界 $Re$ 后,流动经历层流 → 转捩 → 湍流
• 圆柱尾流出现卡门涡街,用 Strouhal 数特征化
• DNS 湍流仿真的网格成本按 $Re^{9/4}$ 标度,高 Re 流动必须用湍流模型

下周将介绍如何用**有限体积法(FVM)**离散此 Navier-Stokes 方程,并比较迎风差分与中心差分格式对精度的影响。