Song's Blog

外观

从 Reynolds 数到湍流

约 2596 字大约 9 分钟

地球流体力学

2026-05-23

流体力学里有一串经常一起出现、但初学时很容易各看各的概念：

  • Reynolds 数；
  • 黏性；
  • 扩散项；
  • 涡旋；
  • 湍流。

它们看上去分属不同层次，但其实可以被一条很清楚的逻辑线串起来。黏性决定流体内部会不会主动抹平速度差异，扩散项告诉我们这种“抹平”在方程里该怎么写，Reynolds 数刻画惯性和黏性谁更强，而当惯性远强于黏性时，流动就更容易卷出旋涡并进一步走向湍流。

这篇文章专门展开 Reynolds 数和湍流直觉。GFD 里的其他无量纲数只在 地球流体力学简介 中概览。

Reynolds 数：惯性效应/黏性效应

Reynolds 数通常写成

Re=ULν,\mathrm{Re}=\frac{UL}{\nu},

其中：

  • UU 是典型速度尺度；
  • LL 是典型长度尺度；
  • ν\nu 是运动黏性系数。

只背这个式子其实意义不大，关键是知道它为什么会长成这样。最直接的出发点是 Navier-Stokes 方程：

∂u∂t+(u⋅∇)u=−1ρ∇p+νΔu.\frac{\partial \bm{u}}{\partial t} +(\bm{u}\cdot\nabla)\bm{u} =-\frac{1}{\rho}\nabla p +\nu\Delta\bm{u}.

这里最值得比较的是两项：

  • (u⋅∇)u(\bm{u}\cdot\nabla)\bm{u}：惯性或平流项；
  • νΔu\nu\Delta\bm{u}：黏性扩散项。

若速度量级是 UU、长度尺度是 LL，那么

∇∼1L,Δ∼1L2.\nabla\sim\frac{1}{L}, \qquad \Delta\sim\frac{1}{L^2}.

于是两项的量级可以估成

(u⋅∇)u∼U2L,νΔu∼νUL2.(\bm{u}\cdot\nabla)\bm{u}\sim \frac{U^2}{L}, \qquad \nu\Delta\bm{u}\sim \nu\frac{U}{L^2}.

它们的比值就是

U2/LνU/L2=ULν=Re.\frac{U^2/L}{\nu U/L^2}=\frac{UL}{\nu}=\mathrm{Re}.

Reynolds 数可以理解成：

Reynolds 数衡量的是惯性效应和黏性效应的相对强弱。

Reynolds 数大和小，会把流动带向两种很不同的世界

当

Re≪1\mathrm{Re}\ll 1

时，黏性主导，流动往往比较平滑，扰动不容易被放大，旋涡也较难维持。蜂蜜流动、微流体器件、细菌游动这些问题，通常都更接近这个极限。若进一步考虑稳态、并忽略惯性项，那么 Navier-Stokes 方程可近似成

−∇p+νΔu=0,-\nabla p+\nu\Delta\bm{u}=0,

这类流动通常称为 Stokes flow 或 creeping flow。

而当

Re≫1\mathrm{Re}\gg 1

时，惯性项会更强。流体更容易把原有结构继续搬运、卷曲和拉伸，而不是立刻被黏性抹平。于是你会开始看到：

  • 小扰动被不断放大；
  • 旋涡更容易形成并维持；
  • 流动出现明显的多尺度结构；
  • 在合适的边界条件和扰动下，流动更容易进入湍流状态。

大气、海洋、飞机绕流和河流急流都属于这类问题。比如大气里若取

U∼10 m/s,L∼106 m,ν∼10−5 m2/s,U\sim 10\ \mathrm{m/s}, \qquad L\sim 10^6\ \mathrm{m}, \qquad \nu\sim 10^{-5}\ \mathrm{m^2/s},

则有

Re∼1012.\mathrm{Re}\sim 10^{12}.

这已经大到足以说明：对真实大气而言，分子黏性几乎不可能在大尺度上直接把流动维持得很“整齐”。

黏性到底意味着什么

黏性最直观的理解是：

流体内部抵抗相对滑动的能力。

也可以把它粗略理解成流体内部的“摩擦”。水的黏性小，蜂蜜的黏性大，所以两者在同样受力条件下会表现出完全不同的流动状态。

通常会见到两种黏性系数。

动力黏性系数

记作

μ,\mu,

单位是

Pa⋅s.\mathrm{Pa\cdot s}.

运动黏性系数

记作

ν,\nu,

并满足

ν=μρ,\nu=\frac{\mu}{\rho},

其中 ρ\rho 是流体密度，单位是

m2/s.\mathrm{m^2/s}.

在方程和尺度分析里，更常直接出现的是 ν\nu，因为它更自然地和扩散项联系在一起。

黏性项和 Laplace 算子

扩散方程最基本的形式是

ut=DΔu.u_t=D\Delta u.

对标量 uu 来说，这表示某个物理量在空间中的扩散；而在 Navier-Stokes 方程里，形式相似的 νΔu\nu\Delta\bm{u} 描述的则是动量扩散。两者作用在不同对象上，但背后的平滑机制是一样的。

从物理上看，它表达的是一件很简单的事：把局部过高的地方压低，把局部过低的地方抬高，让分布逐渐变平滑。

先看一维情形：

Δu=uxx.\Delta u=u_{xx}.

若某点附近是局部峰值，那么曲线向下凹，通常有

uxx<0.u_{xx}<0.

于是 ut<0u_t<0，也就是峰值会下降。

反过来，若某点附近是局部谷值，则通常有

uxx>0,u_{xx}>0,

于是 ut>0u_t>0，谷值会上升。

所以扩散在做的事情，本质上就是：

  • 高的地方往下掉；
  • 低的地方往上补；
  • 整个分布越来越均匀。

如果再用差分近似看一下，

uxx(x)≈u(x+h)−2u(x)+u(x−h)h2,u_{xx}(x)\approx\frac{u(x+h)-2u(x)+u(x-h)}{h^2},

就会发现它近似等于“邻居平均值减去当前位置”。这意味着 Laplace 算子本质上在推动函数向局部平均靠近。

回到速度方程里的

νΔu,\nu\Delta\bm{u},

它的物理意义就是：

动量会从速度较大的区域向速度较小的区域扩散。

这会带来两个直接后果。

第一，局部过强的速度梯度会被抹平，所以很细小的结构不容易长期保持得非常尖锐。

第二，动能最终会在小尺度上被耗散掉。这也是为什么黏性虽然看起来只是一个“平滑项”，却和耗散、混合、湍流尾端的结构密切相关。

涡旋和涡量：流动为什么会开始“卷”

当我们说流动里出现了旋涡，意思通常是：流体局部绕某个区域发生了明显旋转。台风、龙卷风、下水道漩涡、烟圈，都是很典型的例子。

描述这种局部旋转强度的量叫涡量。若二维速度场为

u=(u,v),\bm{u}=(u,v),

则相对涡量写成

ω=∂xv−∂yu.\omega=\partial_x v-\partial_y u.

若只是整体平移而不发生局部转动，那么涡量可以为零；而一旦速度场开始出现剪切、卷曲和绕转，涡量就会变得非零。

湍流不只是“乱”，而是多尺度涡旋的强相互作用

很多时候，人们会把湍流简单理解成“很乱的流动”。这不算错，但还不够准确。更接近本质的说法是：

湍流是大量不同尺度涡旋之间的强非线性相互作用。

它通常同时具有几个特征：

  • 非线性很强，平流项不能忽略；
  • 存在从大到小的一整串尺度；
  • 对初始扰动敏感，长期行为难以精确预测；
  • 小尺度最终承担主要耗散。

所以湍流不是“有没有旋涡”这么简单。一个单独的稳定旋涡并不等于湍流；真正的湍流意味着很多尺度的结构同时在生成、变形、相互作用和破碎。

为什么高 Reynolds 数更容易走向湍流

现在可以把前面的逻辑线完整接起来了。

黏性项试图平滑速度场，抹掉局部梯度；惯性项则试图继续搬运、拉伸和卷曲已有结构。若 Reynolds 数不大，黏性来得及把很多小扰动压下去，流动就更容易保持层流状态。

但当 Reynolds 数很大时，黏性太弱，来不及快速抹平这些结构。于是：

  • 小扰动不容易消失；
  • 旋涡更容易生成；
  • 已有旋涡会继续卷出更细的结构；
  • 流动也就更容易出现明显的多尺度混乱，甚至发生转捩并进入湍流。

这就是“高 Reynolds 数更容易湍流”背后的物理图像。它不是因为 Reynolds 数本身神秘地触发了某种开关，而是因为惯性和黏性的竞争关系彻底变了。至于流动是否真的转入湍流，还会受到边界条件、几何形状和初始扰动大小的影响。

能量级串：为什么湍流总会把事情越卷越细

湍流里最关键的过程之一，是能量在尺度之间传递。最常见的三维图像是：

大尺度→小尺度.\text{大尺度} \to \text{小尺度}.

大尺度运动先生成较大的涡旋，这些涡旋再被拉伸、卷曲、破碎，继续把能量传给更小的结构。这个过程常被称为 energy cascade，也就是能量级串。

直到尺度小到某个程度，黏性扩散终于变得足够强，能把这些细小结构耗散掉。也就是说，湍流并不是“完全没有黏性”，而是：

  • 大尺度主要由惯性主导；
  • 小尺度最终由黏性接管并耗散。

这也是为什么讨论湍流时，黏性虽然经常在大尺度分析里看起来很弱，却永远不能真的被忘掉。

二维和三维湍流并不完全一样

在三维湍流里，更典型的是能量向小尺度传递；而在理想二维、不可压、足够高 Reynolds 数的框架下，常讨论相反的一部分趋势，也就是所谓的 inverse cascade：

小尺度→大尺度.\text{小尺度} \to \text{大尺度}.

这也是为什么二维模型里更容易出现长期存在的大尺度旋涡结构。很多地球流体现象虽然不是真正严格二维，但由于强旋转、强分层和浅层结构的存在，常常会表现出某种“准二维”的特征，所以大尺度涡旋在大气海洋里格外重要。

总结

黏性通过 Laplace 型扩散项平滑速度场，Reynolds 数决定这种平滑能不能压住惯性产生的结构放大；当压不住时，旋涡会不断生成和破碎，流动就会走向湍流。

  • Reynolds 数衡量惯性和黏性的相对强弱；
  • 黏性本质上会抹平速度差异；
  • Laplace 算子描述的正是这种局部平均化趋势；
  • 涡量衡量局部旋转；
  • 湍流是多尺度涡旋的强非线性相互作用。

这几点基本构成了从 Navier-Stokes 方程走向湍流直觉的第一步。

版权归属：Guisong Wu

许可证：署名 4.0 国际 (CC-BY-4.0)