流体力学 · 界面稳定性

Kelvin–Helmholtz 剪切层失稳

两层流体相互滑过时，界面上的波纹什么时候会长大？ 某一波长的波纹，只要剪切推它长大的速度快过重力和表面张力把界面拉平，又没被黏性阻尼抵掉，就会增长。重力的恢复作用来自上下两层的密度差。指数增长只发生在开头；波纹一卷成涡，运动变成非线性，线性增长率就不再适用。

在二维画面里看两层反向流动的流体、它们之间的界面、涡量和速度矢量，旁边是傅里叶谱和增长率图。五个预设工况分别突出密度、重力、表面张力、黏性和初始扰动谱的作用。

剪切时间

t\,\Delta U/L

$0.000$

模态增益

|\eta_k|/|\eta_k(0)|

$1.000$

线性预测

\gamma_{\mathrm{th}}

$0\,\mathrm{s}^{-1}$

测得拟合

\gamma_{\mathrm{fit}}

$0\,\mathrm{s}^{-1}$

最快已解析模态

n_*

$1$

当前模型阶段 $\mathrm{linear}$

稳定性实验

检验剪切层的稳定性

速度矢量在部分点上显示推动涡片标记点的那个速度场。线性理论包络进入非线性模型后，虚线理论边界会消失。速度跃变

\Delta U

0.56\,\mathrm{m\,s^{-1}}

密度比

\rho_u/\rho_\ell

0.92

黏性阻尼代理

\nu

6.0\times10^{-4}\,\mathrm{m^2\,s^{-1}}

剪切层厚度

\delta

0.034\,\mathrm{m}

重力

g

0.18\,\mathrm{m\,s^{-2}}

表面张力

\sigma

8.0\times10^{-4}\,\mathrm{N\,m^{-1}}

种子波长

\lambda_0

0.50\,\mathrm{m}

种子振幅

\eta_0

0.009\,\mathrm{m}

C 面板拿早期振幅的对数对时间做拟合。理论直线画到图中标出的模型切换处为止。

参考：Miles，JFM 3（1957），DOI 10.1017/S0022112057000567；Krasny，JCP 65（1986），DOI 10.1016/0021-9991(86)90210-X。

早期增长由 $\eta_k(t)=\eta_k(0)e^{\gamma_k t}$ 和重力—毛细 Kelvin–Helmholtz 色散关系给出。一旦 $|\eta_k|/\lambda_k\ge 0.075$ ，画面明确切换到正则化的 Birkhoff–Rott 涡片。这个后期模型不解析混合、黏性扩散和界面破碎。

物理教程

Kelvin–Helmholtz 失稳：从微扰增长到涡片卷起

背景两层不可压缩流体沿一条水平界面反向流动。界面上如果出现波数为 $k$ 的微小起伏，速度跃变会把它放大；稳定的密度分层、重力和表面张力则抵抗变形。Lab 先逐个预测每个 Fourier 模态的增长；小振幅假设一失效，就明确切换到另一个模型。

为什么值得理解云层、海洋内波、喷流边缘和天体等离子体里都能看到剪切卷起的波浪。样子相像，不代表背后是同一个模型；真正能检验的，是增长率和模型成立的范围。

先抓重点

核心问题: 速度差提供的放大，什么时候能压过重力、表面张力和黏性阻尼？
一句话直觉: 线性阶段里，每个波数各自独立演化。增长率 $\gamma_k>0$ 只说明微扰在指数增长；界面一卷起来，模态开始耦合，几何也变得非线性，线性预测就必须停下。

核心数学模型

两层重力—毛细色散关系

\omega=k\frac{\rho_u U_u+\rho_\ell U_\ell}{\rho_u+\rho_\ell}\pm\sqrt{\frac{(\rho_\ell-\rho_u)gk+\sigma k^3}{\rho_u+\rho_\ell}-\frac{\rho_u\rho_\ell}{(\rho_u+\rho_\ell)^2}(\Delta U)^2k^2}

根号里第一项来自稳定分层和表面张力的恢复作用，第二项是剪切的驱动。根号里的值一旦变负，频率就有了虚部，扰动开始指数增长。

有限厚度和黏性的教学修正

\gamma_k=\max\!\left[0,\sqrt{\frac{\rho_u\rho_\ell}{(\rho_u+\rho_\ell)^2}(\Delta U)^2k^2\operatorname{sech}(k\delta)-\frac{(\rho_\ell-\rho_u)gk+\sigma k^3}{\rho_u+\rho_\ell}}-2\nu k^2\right]

厚度衰减因子和黏性项是写得很直白的低阶教学近似，不是黏性 Orr–Sommerfeld 方程的本征值解。它们让你看到短波对厚度和耗散有多敏感。

线性阶段的振幅增长

\eta_k(t)=\eta_k(0)e^{\gamma_k t},\qquad \gamma_{\mathrm{fit}}=\frac{\mathrm{d}}{\mathrm{d}t}\ln|\eta_k|

增长率图拟合对数振幅的斜率，再和理论预测比较。进入非线性阶段后，理论直线停在模型切换的标记处。

正则化涡片阶段

\frac{\partial \overline{z}}{\partial t}=\frac{1}{2\pi i}\operatorname{PV}\!\int\frac{\gamma(\alpha')}{z(\alpha)-z(\alpha')}\,\mathrm{d}\alpha'

Birkhoff–Rott 速度把涡片自己诱导出的流动反馈给界面。程序用周期性的 vortex-blob 核和有限的核心半径把奇异点抹平；这个半径是数值处理，不代表真实的物理小尺度。

常见难点

一个增长率管不了整段演化

容易误解早期振幅按指数增长，那么后面的涡旋大小、混合和破碎，也能用同一个增长率算出来。

正确理解线性理论要求振幅远小于波长。振幅与波长之比一到 0.075，Lab 就停止线性预测，后期改标为正则化涡片。

有限厚度不是锐利界面

容易误解调大剪切层厚度只会让画面变模糊。

正确理解剪切层有厚度时，短波感受到的有效速度跃变会变弱。这里用 $\operatorname{sech}(k\delta)$ 展示这个趋势；要得到准确数值，得针对具体的速度剖面求解本征值问题。

动手实验

01
测量单模指数增长
选择“单模种子”，暂停，然后反复点“推进 30 步”，对比“增长率审计”图里的实线和虚线。
你会看到： 振幅还小的时候，对数振幅大致沿直线上升，拟合出的增长率接近理论值。模型切换以后，两条线不再被要求重合。
02
寻找稳定阈值
选择“稳定分层”，逐步提高速度跃变或降低重力，留意最快模态和当前阶段。
你会看到： 恢复作用占上风时，扰动只振荡并衰减。越过阈值后出现增长最快的模态，阶段标签变成“线性增长”。
03
让谱自己选出尺度
选择“噪声种子谱”，分别加大剪切层厚度和黏性阻尼，比较傅里叶谱。
你会看到： 短波受到的剪切驱动更强，但厚度、表面张力和黏性对它的压制也更强。最终胜出的尺度，来自这几种作用的较量。

概念灵感：FluidMotion。模拟和图形都是原创，模型的适用范围见上文。