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

试着将厨房的水龙头拧开一半。起初，水流是一根平滑的柱体。但从大约 10 厘米处开始，它就会断裂成规则的水滴。虽然这是每天都能看到的场景，但要回答为什么连续体户自发断裂并不容易。本文将追随 1849 年 Joseph Plateau 和 1879 年 Lord Rayleigh 给出的经典解释，并探讨在当今发动机喷嘴和喷墨设计中仍在使用的 韦伯数 (Weber number，惯性力/表面张力比) 和 奥内佐格数 (Ohnesorge number，粘性力/惯性与表面张力之积的平方根比) 如何扩展了这些结论。阅读完本文，你将能用一张机制图（Regime Map）回答“为什么有的射流会断裂成大水滴，而有的射流会像雾一样粉碎”。

为什么连续体选择断裂？#

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

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

Plateau 的几何直觉#

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

\Delta A \propto \left[1 - \left(\frac{2\pi R}{\lambda}\right)^{2}\right]

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

简单来说，该条件为：

\lambda > 2\pi R \ \Rightarrow \ \text{不稳定}, \qquad \lambda < 2\pi R \ \Rightarrow \ \text{稳定}

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

Rayleigh 的线性稳定性#

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

\sigma^{2} = \frac{\sigma_{s}}{\rho R^{3}} \, kR \, \bigl(1 - (kR)^{2}\bigr) \, \frac{I_{1}(kR)}{I_{0}(kR)}

其中 $\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 = \frac{\rho_{g} U^{2} D}{\sigma_{s}}, \qquad Oh = \frac{\mu_{l}}{\sqrt{\rho_{l} \sigma_{s} D}}

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

机制	主导力	液滴尺寸 vs 射流直径	碎裂位置
Rayleigh breakup	表面张力	液滴 > 射流直径	远离喷嘴
First wind-induced	表面张力 + 气体惯性	液滴 ≈ 射流直径	距喷嘴较近
Second wind-induced	气体惯性占优	液滴 < 射流直径	紧邻喷嘴
Atomization	纯惯性破坏	液滴 ≪ 射流直径	喷嘴出口处立即发生

当 $Oh$ 较小（ $< 0.1$ ）时，机制仅由 $We$ 决定。当 $Oh$ 增大时，需要进行粘性修正，工程实践中常用 $Oh(We)^{0.5}$ 等复合数来划定边界。水龙头和喷墨通常处于 Rayleigh breakup 区域。柴油喷油器则追求 Atomization，将 $We$ 提高到 $10^{5}$ 以上。

直观观察增长率曲线#

在下方的模拟中亲自操作一下。通过滑块改变波数 $kR$ ，为射流植入该波长的扰动，增长率曲线上的点将随之移动。

파수 kR =0.70Ohnesorge =0.02

kR은 파장 λ를 반지름으로 나눈 값(λ = 2πR/kR). kR < 1이면 Plateau 조건을 만족해 불안정, 가장 빠른 성장은 kR ≈ 0.697에서 발생한다. Ohnesorge를 올리면 점성이 커져 성장이 느려진다.

将 $kR$ 从 0.2 提高到 0.7，增长率会剧增，液柱会更快地碎裂成液滴。当 $kR$ 超过 1.1 时，曲线贴近零点，扰动消散，液柱保持平滑。将奥内佐格数从 0 提高到 0.4，虽然增长速度变慢，但不稳定 $kR$ 范围保持不变 —— 你可以亲眼确认粘性“无法抹除边界”。

使用 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 毫秒。这与日常观察到的水龙头下方 10~15 厘米处产生水滴的现象基本一致。乘以速度 $U = 1\,\text{m/s}$ 即可得出落下距离。

核心三要点#

液柱仅在长于周长的波长（ $\lambda > 2\pi R$ ）下增长，且增长最快的波长固定为 $\lambda_{\max} \approx 9 R$ 。
粘性会降低增长率，但不会改变不稳定边界 —— 奥内佐格数调节的是“碎裂速度”而非“是否碎裂”。
工程喷雾设计在 $We$ 和 $Oh$ 平面上导航四个机制 —— Rayleigh / first·second wind-induced / atomization，通过调节 $We$ 来控制目标液滴尺寸。