水柱为何会断裂成水滴？—— Plateau–Rayleigh 不稳定性与碎裂机制

厨房水龙头拧开一半。开头一掌是光滑液柱。大约 10 cm 处断成规则水滴。每天都看见，连续体为什么自己断开，并不是马上答得出来。

1849 年 Plateau 只留下表面积会减小的波长。1879 年 Rayleigh 用色散关系从中挑出长得最快的一支。发动机喷嘴和喷墨今天仍用的刻度是 韦伯数（惯性力/表面张力）和 奥内佐格数（粘性/惯性与表面张力）。有的射流留下粗液滴，有的变成雾。

拨动波数 $kR$ 滑块，液柱上种下的波长和增长率曲线上的点一起移动。

把 $kR$ 从 0.2 抬到 0.7，增长变快，液柱更早碎成滴；$kR$ 超过 1.1 时曲线贴在底部，液柱保持光滑。奥内佐格数从 0 提到 0.4，只是变慢，不稳定的 $kR$ 范围不动。

为什么连续体选择断裂？

表面张力是一种倾向于减小表面积的力。对于等体积的液体，球形的表面积最小。因此，在无重力环境下，一团水会形成球体。

那么，对于等体积的长液柱和将其按一定间隔切开形成的液滴链，哪种情况的表面积更小呢？与直觉相反，液滴链的表面积更小 —— 但前提是切割间隔长于液柱的周长。这就是液柱自发碎裂的原因：这是向低能量状态演化的自然结果。

Plateau 的几何直觉#

Plateau 计算了半径为 $R$ 的圆柱体在受到轴向波长为 $\lambda$ 的微弱变形时表面积的变化。

这里 $2\pi R/\lambda$ 是无量纲波数 $kR$（波数 $k = 2\pi/\lambda$ 与半径的乘积）。如果该值小于 1，则 $\Delta A < 0$，即随着扰动增大，表面积减小，这在能量上是有利的。反之，如果大于 1，则 $\Delta A > 0$，扰动会消散。

简单来说，该条件为：

只有长于圆柱周长的波长才会增长。 这一准则被称为 Plateau 极限。

Rayleigh 的线性稳定性#

Plateau 仅回答了“哪些波长可以增长”。30 年后，Rayleigh 通过在无粘性 Navier–Stokes 方程中引入微小扰动 $\eta = \epsilon e^{\sigma t} \cos(kx)$ 并寻找色散关系，找到了“哪个波长增长最快”。

其中 $\sigma$ 是增长率（单位 $s^{-1}$），$\sigma_{s}$ 是表面张力（N/m），$\rho$ 是液体密度，$I_{0}, I_{1}$ 是修正贝塞尔函数。右侧括号中的 $1-(kR)^{2}$ 正是 Plateau 极限。

最大值出现在 $kR \approx 0.697$ 处，对应的波长 $\lambda_{\max} \approx 9.02 R$。液柱更倾向于以其半径的约 9 倍周期 进行切割。这一周期几乎完美地解释了实际水龙头水滴的间距。

如果考虑粘性（由 Chandrasekhar 总结），$\sigma$ 会减小，但 $kR < 1$ 这一不稳定边界本身并不会改变。粘性仅会 减缓增长速度，而不会影响碎裂的必要条件。

韦伯与奥内佐格：碎裂机制图

对于水龙头这类低速流动，上述解释已足够，但在火箭燃烧室或柴油喷油器的高速射流中，周围气体的惯性也会参与其中。这里引入了两个无量纲数：

$We$ 是气体惯性与表面张力的比值，$Oh$ 是液体粘性与（惯性及表面张力）的比值。综合实验数据，射流随着速度增加会经历四个机制阶段。

当 $Oh$ 较小（$< 0.1$）时，机制仅由 $We$ 决定。当 $Oh$ 增大时，需要进行粘性修正，工程实践中常用 $Oh(We)^{0.5}$ 等复合数来划定边界。水龙头和喷墨通常处于 Rayleigh breakup 区域。柴油喷油器则追求 Atomization，将 $We$ 提高到 $10^{5}$ 以上。射流已经变成液滴之后，气体再把液滴撕开的计算，由 韦伯数与 Taylor 类比破碎接着做。这里的 $We$ 用来选机制，那篇笔记跟踪的是单个液滴的振动。

一个波长的体积决定液滴直径

最快波长是 $\lambda_{\max} \approx 9.02 R$。被掐下的一段体积按圆柱 $\pi R^{2}\lambda$ 算，就是 $9.02\pi R^{3}$。同体积的球满足 $\frac{4}{3}\pi a^{3} = 9.02\pi R^{3}$，于是 $a^{3}/R^{3} = 9.02 \times 3/4 = 6.765$。立方根约为 1.89。液滴直径大约是射流直径的 1.9 倍。表里只有 Rayleigh 一行写“液滴比射流粗”，数字上的理由就在这段体积。

不稳定开启的长度是周长，而且只有 $\lambda > 2\pi R$ 会增长，大约 $6.28 R$。$9.02/6.28 \approx 1.44$，最快周期只比刚开始失稳的波长长 44%。喷嘴若只给更短的皱褶，液柱就保持光滑。

增长率平方前面的因子是 $\sigma_{s}/(\rho R^{3})$。乘在 $kR$ 上的部分只依赖波数，所以峰值 $\sigma$ 按 $R^{-3/2}$ 下降。代码写下的情形是半径 1 mm，$\sigma \approx 380\,\mathrm{s}^{-1}$，时间约 18 ms。乘以 1 m/s，断裂位置是 1.8 cm。这不是水槽下看到的 10–15 cm。液柱更粗时，同一公式会把时间拉长。

使用 Python 寻找增长最快的波长#

我们可以通过数值搜索 Rayleigh 色散关系来寻找最优波数。使用 scipy 中的 i0 和 i1。

import numpy as np
from scipy.special import i0, i1
from scipy.optimize import minimize_scalar
 
def rayleigh_sigma(kR, sigma_s=0.072, rho=1000.0, R=1e-3):
    """无粘 Plateau-Rayleigh 增长率 (s^-1)。"""
    if kR <= 0 or kR >= 1:
        return 0.0
    coeff = sigma_s / (rho * R**3)
    return np.sqrt(coeff * kR * (1 - kR**2) * i1(kR) / i0(kR))
 
def fastest_mode(R, sigma_s=0.072, rho=1000.0):
    """增长最快的波数及对应波长。"""
    res = minimize_scalar(
        lambda x: -rayleigh_sigma(x, sigma_s, rho, R),
        bounds=(0.01, 0.99), method='bounded'
    )
    kR_star = res.x
    wavelength = 2 * np.pi * R / kR_star
    sigma = rayleigh_sigma(kR_star, sigma_s, rho, R)
    return kR_star, wavelength, sigma
 
if __name__ == "__main__":
    R = 1e-3          # 1 mm 半径液柱
    kR_star, lam, sig = fastest_mode(R)
    t_break = np.log(R / 1e-6) / sig   # 假设初始扰动 ε_0 = 1 μm
    print(f"kR* = {kR_star:.3f}")
    print(f"lambda_max = {lam*1e3:.2f} mm  (= {lam/R:.2f} R)")
    print(f"sigma_max = {sig:.1f} /s")
    print(f"breakup time ~ {t_break*1e3:.1f} ms")

运行结果为 kR* ≈ 0.697，λ_max ≈ 9.02 R，σ_max ≈ 380 /s，碎裂时间约 18 毫秒。乘以 $U = 1\,\text{m/s}$ 得到的下落距离是 1.8 cm。把水龙头下方 10~15 厘米当成这个算例的核对，长度会对不上。初始凹凸取 1 μm 时，振幅要长到 1000 倍，把 $\ln 1000 \approx 6.9$ 除以 380，18 ms 会再出现一次。

用 1 mm 的 18 ms 去对水龙头高度，位置会对错#

1.8 cm 和 10 cm 不是同一根液柱。$\sigma$ 按 $R^{-3/2}$ 变小，还没量半径就把 18 ms 写成设计值，碎裂长度会偏短。$\lambda_{\max} \approx 9 R$ 选的是波长，不是下落距离的通用尺子。

要留下粗液滴，就停在 Rayleigh，把波长固定在大约 9 倍半径附近。要雾，就把 $We$ 抬到 $10^{5}$ 以上，进入雾化。粘性把这条边界歪多少，单看 $Oh$ 不够。惯性相对粘性的刻度在 Reynolds 数如何改写流动尺度里。