---
title: "一条微小波纹怎样卷成剪切涡旋互动实验｜Sandbox Physics"
description: "两层流体相互滑过时，界面上的波纹什么时候会长大？某一波长的波纹，只要剪切推它长大的速度快过重力和表面张力把界面拉平，又没被黏性阻尼抵掉，就会增长。"
canonical_url: "https://sandboxphysics.com/zh/labs/kelvin-helmholtz-instability/"
language: "zh-CN"
page_type: "lab"
alternate_languages:
  en: "https://sandboxphysics.com/labs/kelvin-helmholtz-instability/"
last_modified: "2026-10-27"
---

流体力学 · 界面稳定性

# Kelvin–Helmholtz 剪切层失稳

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

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

> 交互仿真：在浏览器中打开本页，即可调整参数并运行模型。

早期增长由 $\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)$ 展示这个趋势；要得到准确数值，得针对具体的速度剖面求解本征值问题。

### 动手实验

1. 01

   #### 测量单模指数增长

   选择“单模种子”，暂停，然后反复点“推进 30 步”，对比“增长率审计”图里的实线和虚线。

   **你会看到：** 振幅还小的时候，对数振幅大致沿直线上升，拟合出的增长率接近理论值。模型切换以后，两条线不再被要求重合。

2. 02

   #### 寻找稳定阈值

   选择“稳定分层”，逐步提高速度跃变或降低重力，留意最快模态和当前阶段。

   **你会看到：** 恢复作用占上风时，扰动只振荡并衰减。越过阈值后出现增长最快的模态，阶段标签变成“线性增长”。

3. 03

   #### 让谱自己选出尺度

   选择“噪声种子谱”，分别加大剪切层厚度和黏性阻尼，比较傅里叶谱。

   **你会看到：** 短波受到的剪切驱动更强，但厚度、表面张力和黏性对它的压制也更强。最终胜出的尺度，来自这几种作用的较量。

物理结论

Kelvin–Helmholtz 失稳是各模态之间一场可以量化的较量：剪切让扰动增长，分层和表面张力把界面拉回，耗散筛掉短波；界面卷起之后，由另一个非线性模型接手。

概念灵感：[FluidMotion](http://www.bugman123.com/FluidMotion/index.html)。模拟和图形都是原创，模型的适用范围见上文。

场、流体与物质 · 入门路径

## 从电场到电子显微镜

先看懂电场、学会操控带电粒子，再算出电子的波长，用电子束观察晶体和表面。

[从这里开始 → **看见电场** 在偶极子周围移动探针，看电场是不是总指向电势地形的下坡方向。 观察 · 8 分钟](https://sandboxphysics.com/zh/labs/electric-field/)

1. [1 看见电场](https://sandboxphysics.com/zh/labs/electric-field/)
2. [2 电荷在场中转弯](https://sandboxphysics.com/zh/labs/charged-particle/)
3. [3 电子波长与加速电压](https://sandboxphysics.com/zh/labs/electron-wavelength/)
4. [4 电子衍射与倒空间](https://sandboxphysics.com/zh/labs/electron-diffraction/)
5. [5 TEM 明场、暗场与选区衍射](https://sandboxphysics.com/zh/labs/tem-mode-selector/)
6. [6 SEM 二次电子表面形貌](https://sandboxphysics.com/zh/labs/sem-topography/)

[查看完整路径 →](https://sandboxphysics.com/zh/topics/fields-and-matter/#starter-path)