水流はなぜ滴に分裂するのか — プラトー・レイリー不安定性と分裂レジーム

表面張力が水柱を液滴へと切り分ける原理と、噴霧設計のための無次元数マップ

キッチンの水道の蛇口を半分だけ開けてみてください。最初は滑らかな水柱が見えます。しかし、10cmほど落ちた地点から、規則的な液滴へと分裂し始めます。毎日目にする光景ですが、なぜ連続体が自ら途切れるのかという問いに答えるのは簡単ではありません。この記事では、1849年のジョゼフ・プラトー（Joseph Plateau）と1879年のレイリー卿（Lord Rayleigh）による古典的な説明を辿り、エンジン噴射ノズルやインクジェットの設計で今日でも使われている ウェーバー数（慣性力／表面張力比）と オーネゾルゲ数（粘性力／（慣性・表面張力）比）が、どのようにその結果を拡張しているのかを見ていきます。最後まで読めば、「なぜあるジェットは太い滴に分裂し、あるジェットは霧のように砕けるのか」を、一つのレジームマップで説明できるようになります。

なぜ連続体は分裂する道を選ぶのか#

表面張力は表面積を減らそうとする力です。同じ体積の液体であれば、球形であるときが最も表面積が小さくなります。そのため、無重力状態の水の塊は球形になります。

それでは、同じ体積の長い水柱と、その水を一定の間隔で切り分けて作った液滴の鎖では、どちらの表面積が小さいでしょうか。直感に反して、液滴の鎖の方が小さくなります。ただし、切り分ける間隔が水柱の円周よりも長い場合に限ります。水柱が自ら分裂するのはこれが理由です。エネルギーの低い状態へと転がり落ちようとする、自然な結果なのです。

プラトーの幾何学的直感#

プラトーは、半径 $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 limit）です。

レイリーの線形安定性#

プラトーは「どの波長が成長できるか」という問いにのみ答えました。「どの波長が最も速く成長するか」については、その30年後、レイリーが非粘性ナビエ・ストークス方程式に微小な摂動 $\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}$ は、まさにプラトー限界そのものです。

最大値は $kR \approx 0.697$ で得られ、対応する波長は $\lambda_{\max} \approx 9.02 R$ です。水柱は、自身の半径の約 9倍の周期 で切り分けられることを好みます。この周期が、実際の水道の蛇口から出る液滴の間隔をほぼ正確に説明しています。

粘性を考慮に入れると、チャンドラセカール（Chandrasekhar）が整理した通り、 $\sigma$ は減少しますが、 $kR < 1$ という不安定性の境界自体は変わりません。粘性は 成長を遅らせる だけであり、分裂の必要条件には影響を与えません。

ウェーバーとオーネゾルゲ：分裂レジームのマップ#

水道のようにゆっくり流れる場合には上記の説明で十分ですが、ロケット燃焼器やディーゼルインジェクターの高圧噴射では、周囲の気体の慣性も関わってきます。ここで、2つの無次元数が登場します。

We = \frac{\rho_{g} U^{2} D}{\sigma_{s}}, \qquad Oh = \frac{\mu_{l}}{\sqrt{\rho_{l} \sigma_{s} D}}

$We$ は気体の慣性と表面張力の比、 $Oh$ は液体の粘性と（慣性・表面張力）の比です。実験データを統合すると、ジェットの速度が上がるにつれて、4段階のレジームを通過します。

レジーム	主な力	液滴サイズ 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 で最も速い波長を探す#

ベッセル関数を用いたレイリーの分散関係をそのままコードに落とし込み、最適な波数を数値的に探索します。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):
    """非粘性プラトー・レイリー成長率 (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 ms です。蛇口から 10～15 cm 下で液滴ができる日常の観察結果とほぼ一致します。速度 $U = 1\,\text{m/s}$ を掛ければ、落下距離がすぐに計算できます。

要約#

水柱は、円周よりも長い波長（ $\lambda > 2\pi R$ ）だけが成長し、最も速い波長は $\lambda_{\max} \approx 9 R$ で固定されます。
粘性は成長率を減少させますが、不安定性の境界は変えません。つまり、オーネゾルゲ数は「分裂の速度」を調節するだけであり、「分裂するかどうか」を変えることはできません。
実務的な噴霧設計では、 $We$ と $Oh$ の平面上の4つのレジーム（Rayleigh / first・second wind-induced / atomization）を使い分け、狙いの液滴サイズに合わせて $We$ を調整します。