プラトー・レイリー不安定性

プラトー・レイリー不安定性（プラトー・レイリーふあんていせい、テンプレート:Lang-en-short）とは真空中や空気中に円柱状の流体が流れる際に、表面張力の効果により噴流の周長に対応した特定の波長の攪乱が成長する現象である。この効果によって、水道の蛇口から出る水が下にいくほど小さな粒状の液滴に分裂する。ジョゼフ・プラトーによって実験観測された後、1878年に初めてレイリーによって理論的に研究された。この理論は様々な液体微粒化技術に応用され、例えばインクジェットプリンターの技術にも多大な影響を与えている。

歴史

1800年代にジョゼフ・プラトーが実験的に円形噴流が小さな液滴に分裂する現象を発見したテンプレート:Sfnp。プラトーは下向きに流れる水が分裂する波長（液滴の大きさ）と円柱直径の関係を調べた。その後、1878年と翌1879年にレイリーは理論的にこの関係を導いたテンプレート:Sfnpテンプレート:Sfnp。レイリーは重力と粘性の効果を無視したモデルで計算をしたが、1909年にボーアは粘性の効果を少し考慮したモデルを発表したテンプレート:Sfnp。ギーアとストリークヴェルダテンプレート:Sfnp、またはケラーテンプレート:Sfnpは1983年に重力の効果を取り入れた解析を発表した。

円形噴流における不安定性

レイリーによって導かれた、表面張力による円形噴流の不安定性を以下に示す。ここでは、半径 テンプレート:Math、密度 テンプレート:Mvar、表面張力係数 テンプレート:Mvar の無限に長い円柱を流れる非粘性流体を考え、重力の影響は無視する。圧力 テンプレート:Math は円柱内で一定であり、境界における表面張力による法線応力のバランスによって

$p_0=\sigma \mathrm{\nabla }\cdot 𝐧\mathbf{=}\frac{\sigma }{R_0}$

と計算できる。ここで、界面において微小な節状の摂動の発達を考える。これにより、支配方程式の線形化ができる。攪乱を加えた柱状表面は以下の形で書ける。

$\tilde{R}=R_0+\epsilon e^{\omega t+ikz}$

ここで、攪乱の振幅は テンプレート:Math であり、テンプレート:Mvar は不安定成長率、テンプレート:Mvar は テンプレート:Mvar 方向の攪乱の波数である。節状の摂動の対応する波長は テンプレート:Math である。速度の摂動の動径方向成分を テンプレート:Mvar、軸方向成分を テンプレート:Mvar、圧力の摂動を テンプレート:Mvar で表す。これらの摂動場をナビエ・ストークス方程式に代入し、テンプレート:Mvar のオーダーの項のみを残すと

$\begin{array}{cc}\frac{\partial {\tilde{u}}_r}{\partial t}& =-\frac{1}{\rho }\frac{\partial \tilde{p}}{\partial r}\\ \frac{\partial {\tilde{u}}_z}{\partial t}& =-\frac{1}{\rho }\frac{\partial \tilde{p}}{\partial z}\end{array}$

となる。また、線形化された連続の方程式は

$\frac{\partial {\tilde{u}}_r}{\partial t}+\frac{{\tilde{u}}_r}{r}+\frac{\partial {\tilde{u}}_z}{\partial z}=0$

となる。ここで、速度や圧力の攪乱は表面攪乱の式と同じ形をとるとすると、速度と圧力の攪乱は

${\tilde{u}}_r=R(r)e^{\omega t+ikz},{\tilde{u}}_z=Z(r)e^{\omega t+ikz},\tilde{p}=P(r)e^{\omega t+ikz}$

と書ける。これを上の3つの式に代入することで、摂動場を支配する線形化された方程式は

$\begin{array}{cc}\omega R& =-\frac{1}{\rho }\frac{dP}{dr}\\ \omega Z& =-\frac{ik}{\rho }P\\ \frac{dR}{dr}+\frac{R}{r}+ikZ& =0\end{array}$

となる。これらより、テンプレート:Math の微分方程式が以下のように得られる。

$r^2\frac{d^{2}R}{d{r}^2}+r\frac{dR}{dr}-(1+{\left(kr\right)}^2)R=0$

これは1次の修正ベッセル方程式に対応し、解はそれぞれ第一種 テンプレート:Math、第二種 テンプレート:Math の修正ベッセル関数で記述される。テンプレート:Math で テンプレート:Math なので、取りうる テンプレート:Math の形は

$R(r)=C{I}_1(kr)$

である。ここで、テンプレート:Mvar は境界条件を適用することによって決定される定数である。

圧力は、以下のようになる。

$P(r)=-\frac{\omega \rho C}{k}{I}_0(kr)$

ここで、修正ベッセル関数の性質 テンプレート:Math を用いた。

境界条件を適用する。自由表面における運動論的境界条件は

$\frac{\partial \tilde{R}}{\partial t}\sim {\tilde{u}}_r$

であり、この条件を用いると

$C=\frac{\epsilon \omega }{I_1(k{R}_0)}$

が得られる。次に、自由表面における法線応力のつり合いを考えると

$p_0+\tilde{p}=\sigma \mathrm{\nabla }\cdot \overrightarrow{u}$

となる。

噴流表面の曲率半径を テンプレート:Math、テンプレート:Math と書くと、テンプレート:Math で表される。ここで

$\frac{1}{R_1}=\frac{1}{R_0+\epsilon e^{\omega t+ikz}}\simeq \frac{1}{R_0}-\frac{\epsilon }{R_0^2}{e}^{\omega t+ikz}$

$\frac{1}{R_2}=\epsilon k^2{e}^{\omega t+ikz}$

であり

$\tilde{p}=-\frac{\epsilon \sigma }{R_0^2}(1-k^2{R}_0^2)e^{\omega t+ikz}$

が得られる。以上より、下記のような成長率 テンプレート:Mvar と波数 テンプレート:Mvar の分散関係が得られる。

${\omega }^2=\frac{\sigma }{\rho R_0^3}k{R}_0\frac{I_1(k{R}_0)}{I_0(k{R}_0)}(1-k^2{R}_0^2)$

これにより テンプレート:Math のとき、つまり、円柱の円周より大きな波長の攪乱に対して不安定となることがわかる。下に分散関係のグラフ（横軸 テンプレート:Math、縦軸 テンプレート:Mvar）を示す。ただし、縦軸は テンプレート:Mathで規格化してある。

身近な例

プラトー・レイリー不安定性は液体が分裂して液滴となる現象であり、わかりやすい身近な例として尿を対象として紹介されているテンプレート:Sfnp。2015年9月2日、NHK「ためしてガッテン」の「科学で尿ハネをコントロール特集」で、アメリカのランディ・ハード研究員とテッド・トラスコット准教授の尿ハネに関する研究が紹介された。そこでは テンプレート:Val 程度で便器に尿をぶつければ、プラトー・レイリー不安定性により液滴化する前に便器に達することができ飛沫が抑えられると紹介された。

参考文献

• テンプレート:Citation.
• テンプレート:Cite journal
• テンプレート:Cite journal
• テンプレート:Cite journal
• テンプレート:Cite journal
• テンプレート:Cite journal
• テンプレート:Cite journal

外部リンク

• Plateau–Rayleigh Instability – a 3D lattice kinetic Monte Carlo simulation
• Savart–Plateau–Rayleigh instability of a water column – Adaptive numerical simulation
• テンプレート:PDFlink