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

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

歴史

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

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

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

p 0 = σ n = σ R 0 {\displaystyle p_{0}=\sigma \nabla \cdot {\bf {{n}={\frac {\sigma }{R_{0}}}}}}

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

R ~ = R 0 + ε e ω t + i k z {\displaystyle {\tilde {R}}=R_{0}+\varepsilon e^{\omega t+ikz}}

ここで、攪乱の振幅は εR0 であり、ω は不安定成長率、kz 方向の攪乱の波数である。節状の摂動の対応する波長2π/k である。速度の摂動の動径方向成分を ~ur、軸方向成分を ~uz、圧力の摂動を ~p で表す。これらの摂動場をナビエ・ストークス方程式に代入し、ε のオーダーの項のみを残すと

u ~ r t = 1 ρ p ~ r u ~ z t = 1 ρ p ~ z {\displaystyle {\begin{aligned}{\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{aligned}}}

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

u ~ r t + u ~ r r + u ~ z z = 0 {\displaystyle {\frac {\partial {\tilde {u}}_{r}}{\partial t}}+{\frac {{\tilde {u}}_{r}}{r}}+{\frac {\partial {\tilde {u}}_{z}}{\partial z}}=0}

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

u ~ r = R ( r ) e ω t + i k z , u ~ z = Z ( r ) e ω t + i k z , p ~ = P ( r ) e ω t + i k z {\displaystyle {\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つの式に代入することで、摂動場を支配する線形化された方程式は

ω R = 1 ρ d P d r ω Z = i k ρ P d R d r + R r + i k Z = 0 {\displaystyle {\begin{aligned}\omega R&=-{\frac {1}{\rho }}{\frac {dP}{dr}}\\\omega Z&=-{\frac {ik}{\rho }}P\\{\frac {dR}{dr}}+{\frac {R}{r}}+ikZ&=0\end{aligned}}}

となる。これらより、R(r) の微分方程式が以下のように得られる。

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

これは1次の修正ベッセル方程式に対応し、解はそれぞれ第一種 I1(kr)、第二種 K1(kr) の修正ベッセル関数で記述される。r → 0K1(kr) → ∞ なので、取りうる R(r) の形は

R ( r ) = C I 1 ( k r ) {\displaystyle R(r)=CI_{1}(kr)}

である。ここで、C は境界条件を適用することによって決定される定数である。

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

P ( r ) = ω ρ C k I 0 ( k r ) {\displaystyle P(r)=-{\frac {\omega \rho C}{k}}I_{0}(kr)}

ここで、修正ベッセル関数の性質 I '
0
 
(ξ) = I1(ξ)
を用いた。

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

R ~ t u ~ r {\displaystyle {\frac {\partial {\tilde {R}}}{\partial t}}\sim {\tilde {u}}_{r}}

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

C = ε ω I 1 ( k R 0 ) {\displaystyle C={\frac {\varepsilon \omega }{I_{1}(kR_{0})}}}

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

p 0 + p ~ = σ u {\displaystyle p_{0}+{\tilde {p}}=\sigma \nabla \cdot {\vec {u}}}

となる。

噴流表面の曲率半径R1R2 と書くと、σ∇⋅u = (1/R1 + 1/R2) で表される。ここで

1 R 1 = 1 R 0 + ε e ω t + i k z 1 R 0 ε R 0 2 e ω t + i k z {\displaystyle {\frac {1}{R_{1}}}={\frac {1}{R_{0}+\varepsilon e^{\omega t+ikz}}}\simeq {\frac {1}{R_{0}}}-{\frac {\varepsilon }{R_{0}^{2}}}e^{\omega t+ikz}}
1 R 2 = ε k 2 e ω t + i k z {\displaystyle {\frac {1}{R_{2}}}=\varepsilon k^{2}e^{\omega t+ikz}}

であり

p ~ = ε σ R 0 2 ( 1 k 2 R 0 2 ) e ω t + i k z {\displaystyle {\tilde {p}}=-{\frac {\varepsilon \sigma }{R_{0}^{2}}}\left(1-k^{2}R_{0}^{2}\right)e^{\omega t+ikz}}

が得られる。以上より、下記のような成長率 ω と波数 k の分散関係が得られる。

ω 2 = σ ρ R 0 3 k R 0 I 1 ( k R 0 ) I 0 ( k R 0 ) ( 1 k 2 R 0 2 ) {\displaystyle \omega ^{2}={\frac {\sigma }{\rho R_{0}^{3}}}kR_{0}{\frac {I_{1}(kR_{0})}{I_{0}(kR_{0})}}\left(1-k^{2}R_{0}^{2}\right)}

これにより kR0 < 1 のとき、つまり、円柱の円周より大きな波長の攪乱に対して不安定となることがわかる。下に分散関係のグラフ(横軸 kR0、縦軸 ω)を示す。ただし、縦軸は σ/ρR 3
0
 
で規格化してある。

身近な例

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

脚注

[脚注の使い方]
  1. ^ Plateau (1873).
  2. ^ Rayleigh (1878).
  3. ^ Rayleigh (1879).
  4. ^ Bohr (1909).
  5. ^ Geer & Strikwerda (1983).
  6. ^ Keller (1983).
  7. ^ Hurd et al. (2013).
  8. ^ "家族が涙!トイレ問題 大解決SP". ためしてガッテン. 2 September 2015. NHK総合. NHK. 2019年10月16日閲覧

参考文献

  • Plateau, J (1873), Experimental and theoretical statics of liquids subject to molecular forces only .
  • Rayleigh, L (1878). “On the instability of jets”. Proceedings of the London mathematical society s1-10 (1): 4–13. doi:10.1112/plms/s1-10.1.4. NAID 10004661794. 
  • Rayleigh, L (1879). “On the capillary phenomena of jets”. Proc. R. Soc. London 29: 71-97. Bibcode: 1879RSPS...29...71R. doi:10.1098/rspl.1879.0015. 
  • Bohr, N (1909). “Determination of the surface-tension of water by the method of jet vibration”. Phil. Trans. R. Soc. Lond. 209: 281-317. doi:10.1098/rspa.1909.0014. JSTOR 91039. 
  • Geer, JF; Strikwerda, JC (1983). “Vertical slender jets with surface tension”. Journal of Fluid Mechanics 135: 155-169. Bibcode: 1983JFM...135..155G. doi:10.1017/S0022112083003006. 
  • Keller, JB (1983). “Capillary waves on vertical jet”. Journal of Fluid Mechanics 135: 171-173. Bibcode: 1983JFM...135..171K. doi:10.1017/S0022112083003018. 
  • Hurd, R; Hacking, K; Haymore, B; Truscott, T (2013). “Urinal Dynamics”. 66th Annual Meeting of the American Physical Society Division of Fluid Dynamics 58 (18). http://meetings.aps.org/Meeting/DFD13/Event/202554. 

外部リンク

ウィキメディア・コモンズには、プラトー・レイリー不安定性に関連するカテゴリがあります。
  • Plateau–Rayleigh Instability – a 3D lattice kinetic Monte Carlo simulation
  • Savart–Plateau–Rayleigh instability of a water column – Adaptive numerical simulation
  • An MIT lecture on falling fluid jets, including the Plateau -Rayleigh instability (PDF)