レイリー・プレセット方程式

流体力学においてレイリー・プレセット方程式（レイリー・プレセットほうていしき、英語: Rayleigh–Plesset equation）とは、無限遠点まで満たされた液体内における球形の気泡の動力学を記述する常微分方程式である 。 この方程式は一般的には次のように書かれる: P B ( t ) − P ∞ ( t ) ρ L = R d 2 R d t 2 + 3 2 ( d R d t ) 2 + 4 ν L R d R d t + 2 S ρ L R {\displaystyle {\frac {P_{B}(t)-P_{\infty }(t)}{\rho _{L}}}=R{\frac {d^{2}R}{dt^{2}}}+{\frac {3}{2}}\left({\frac {dR}{dt}}\right)^{2}+{\frac {4\nu _{L}}{R}}{\frac {dR}{dt}}+{\frac {2S}{\rho _{L}R}}} ここで、 P B ( t ) {\displaystyle P_{B}(t)} は気泡内部の圧力(均一とする) P ∞ ( t ) {\displaystyle P_{\infty }(t)} は気泡外部の無限遠点における圧力 ρ L {\displaystyle \rho _{L}} は気泡外部の液体の密度(一定とする) R ( t ) {\displaystyle R(t)} は気泡の半径 ν L {\displaystyle \nu _{L}} は気泡外部の液体の動粘度(一定とする) S {\displaystyle S} は気泡外部の液体と気泡内部の気体との間の表面張力 である。 P B ( t ) {\displaystyle P_{B}(t)} が既知で P ∞ ( t ) {\displaystyle P_{\infty }(t)} が与えられているとすると、レイリー・プレセット方程式は時間変動する気泡の半径 R ( t ) {\displaystyle R(t)} について解くことができる。 レイリー・プレセット方程式は球対称の仮定の元でナビエ–ストークス方程式から求められる。 この方程式は1917年に表面張力と粘度を無視することで、ジョン・ウィリアム・ストラット (第3代レイリー男爵)によって初めて求められ、1949年にミルトン・スピノザ・プレセットによって初めて移動するキャビテーション気泡に対して応用された。 レイリー・プレセット方程式は気泡半径を動的パラメータとして完全に第一原理で求められる。 時間 t {\displaystyle t} に依存する半径 R ( t ) {\displaystyle R(t)} を持つ球面の気泡を考えよう。 温度 T B ( t ) {\displaystyle T_{B}(t)} および圧力 P B ( t ) {\displaystyle P_{B}(t)} が均一で均質に分布した水蒸気または気体を含む気泡を仮定する。 気泡外部は一定の密度 ρ L {\displaystyle \rho _{\mathrm {L} }} 、動粘度 μ L {\displaystyle \mu _{\mathrm {L} }} で無限の領域に満たされた液体があるとする。 気泡から無限遠点における温度と圧力をそれぞれ T ∞ {\displaystyle T_{\infty }} 、 P ∞ ( t ) {\displaystyle P_{\infty }(t)} とする。 気泡の中心からの半径方向距離 r {\displaystyle r} において、圧力 P ( r , t ) {\displaystyle P(r,t)} 、温度 T ( r , t ) {\displaystyle T(r,t)} 、および外向きの速度 u ( r , t ) {\displaystyle u(r,t)} が液体の変動するパラメータとなる。 ただし、それらの液体のパラメータは気泡の外部 r ≥ R ( t ) {\displaystyle r\geq R(t)} においてのみ定義されることに注意する。 質量保存の法則により、半径方向の外向きの速度 u ( r , t ) {\displaystyle u(r,t)} は気泡中心からの距離の2乗に反比例するという逆2乗の法則が要求される。 すなわち、 F ( t ) {\displaystyle F(t)} を何らかの時間の関数として、 u ( r , t ) = F ( t ) r 2 {\displaystyle u(r,t)={\frac {F(t)}{r^{2}}}} と表わせる。 気泡表面での質量輸送が無い場合には、表面の速度は u ( R , t ) = d R d t = F ( t ) R 2 {\displaystyle u(R,t)={\frac {\mathrm {d} R}{\mathrm {d} t}}={\frac {F(t)}{R^{2}}}} とならなくてはならないので、 F ( t ) {\displaystyle F(t)} は F ( t ) = R 2 d R d t {\displaystyle F(t)=R^{2}{\frac {\mathrm {d} R}{\mathrm {d} t}}} と求められる。 質量輸送が起こる場合には、気泡内部の質量増加率は d m V d t = ρ V d V d t = ρ V d ( 4 π R 3 / 3 ) d t = 4 π ρ V R 2 d R d t {\displaystyle {\frac {\mathrm {d} m_{V}}{\mathrm {d} t}}=\rho _{V}{\frac {\mathrm {d} V}{\mathrm {d} t}}=\rho _{V}{\frac {\mathrm {d} (4\pi R^{3}/3)}{\mathrm {d} t}}=4\pi \rho _{V}R^{2}{\frac {\mathrm {d} R}{\mathrm {d} t}}} 表される。ここで V {\displaystyle V} は気泡の体積である。 もし、 u L {\displaystyle u_{L}} が r = R {\displaystyle r=R} における液体の気泡に対する相対速度とすると、気泡に入り込む質量は d m L d t = ρ L A u L = ρ L ( 4 π R 2 ) u L {\displaystyle {\frac {\mathrm {d} m_{L}}{\mathrm {d} t}}=\rho _{L}Au_{L}=\rho _{L}(4\pi R^{2})u_{L}} で与えられる。ここで、 A {\displaystyle A} は気泡の表面積である。 すると質量保存 d m v / d t = d m L / d t {\displaystyle \mathrm {d} m_{v}/\mathrm {d} t=\mathrm {d} m_{L}/\mathrm {d} t} から、 u L = ( ρ V / ρ L ) d R / d t {\displaystyle u_{L}=(\rho _{V}/\rho _{L})\mathrm {d} R/\mathrm {d} t} であるから、 u ( R , t ) = d R d t − u L = d R d t − ρ V ρ L d R d t = ( 1 − ρ V ρ L ) d R d t {\displaystyle u(R,t)={\frac {\mathrm {d} R}{\mathrm {d} t}}-u_{L}={\frac {\mathrm {d} R}{\mathrm {d} t}}-{\frac {\rho _{V}}{\rho _{L}}}{\frac {\mathrm {d} R}{\mathrm {d} t}}=\left(1-{\frac {\rho _{V}}{\rho _{L}}}\right){\frac {\mathrm {d} R}{\mathrm {d} t}}} となるので、 F ( t ) {\displaystyle F(t)} は F ( t ) = ( 1 − ρ V ρ L ) R 2 d R d t {\displaystyle F(t)=\left(1-{\frac {\rho _{V}}{\rho _{L}}}\right)R^{2}{\frac {\mathrm {d} R}{\mathrm {d} t}}} と求められる。 多くの場合、液体の密度はその蒸気の密度に比べて大変大きい( ρ L ≫ ρ V {\displaystyle \rho _{L}\gg \rho _{V}} )ので、 F ( t ) {\displaystyle F(t)} は前者の質量輸送の無い場合 F ( t ) = R 2 d R / d t {\displaystyle F(t)=R^{2}\mathrm {d} R/\mathrm {d} t} に近似され、 u ( r , t ) = F ( t ) r 2 = R 2 r 2 d R d t {\displaystyle u(r,t)={\frac {F(t)}{r^{2}}}={\frac {R^{2}}{r^{2}}}{\frac {\mathrm {d} R}{\mathrm {d} t}}} となる。 液体がニュートン流体であると仮定すると、極座標における半径方向の運動に対する非圧縮性ナビエ-ストークス方程式は ρ L ( ∂ u ∂ t + u ∂ u ∂ r ) = − ∂ P ∂ r + μ L [ 1 r 2 ∂ ∂ r ( r 2 ∂ u ∂ r ) − 2 u r 2 ] {\displaystyle \rho _{L}\left({\frac {\partial u}{\partial t}}+u{\frac {\partial u}{\partial r}}\right)=-{\frac {\partial P}{\partial r}}+\mu _{L}\left[{\frac {1}{r^{2}}}{\frac {\partial }{\partial r}}\left(r^{2}{\frac {\partial u}{\partial r}}\right)-{\frac {2u}{r^{2}}}\right]} と与えられる。 動粘度 ν L = μ L ρ L {\displaystyle \nu _{L}={\frac {\mu _{L}}{\rho _{L}}}} を代入し変形すると − 1 ρ L ∂ P ∂ r = ∂ u ∂ t + u ∂ u ∂ r − ν L [ 1 r 2 ∂ ∂ r ( r 2 ∂ u ∂ r ) − 2 u r 2 ] {\displaystyle -{\frac {1}{\rho _{L}}}{\frac {\partial P}{\partial r}}={\frac {\partial u}{\partial t}}+u{\frac {\partial u}{\partial r}}-\nu _{L}\left[{\frac {1}{r^{2}}}{\frac {\partial }{\partial r}}\left(r^{2}{\frac {\partial u}{\partial r}}\right)-{\frac {2u}{r^{2}}}\right]} ここで質量保存から求めた u ( r , t ) {\displaystyle u(r,t)} を代入し、 − 1 ρ L ∂ P ∂ r = 2 R r 2 ( d R d t ) 2 + R 2 r 2 d 2 R d t 2 − 2 R 4 r 5 ( d R d t ) 2 = 1 r 2 ( 2 R ( d R d t ) 2 + R 2 d 2 R d t 2 ) − 2 R 4 r 5 ( d R d t ) 2 {\displaystyle -{\frac {1}{\rho _{L}}}{\frac {\partial P}{\partial r}}={\frac {2R}{r^{2}}}\left({\frac {\mathrm {d} R}{\mathrm {d} t}}\right)^{2}+{\frac {R^{2}}{r^{2}}}{\frac {\mathrm {d} ^{2}R}{\mathrm {d} t^{2}}}-{\frac {2R^{4}}{r^{5}}}\left({\frac {\mathrm {d} R}{\mathrm {d} t}}\right)^{2}={\frac {1}{r^{2}}}\left(2R\left({\frac {\mathrm {d} R}{\mathrm {d} t}}\right)^{2}+R^{2}{\frac {\mathrm {d} ^{2}R}{\mathrm {d} t^{2}}}\right)-{\frac {2R^{4}}{r^{5}}}\left({\frac {\mathrm {d} R}{\mathrm {d} t}}\right)^{2}} を得る。 代入の過程で粘性項が消えたことに注意する。 変数分離して、気泡の境界 r = R {\displaystyle r=R} から r → ∞ {\displaystyle r\rightarrow \infty } まで積分 − 1 ρ L ∫ P ( R ) P ∞ d P = ∫ R ∞ [ 1 r 2 ( 2 R ( d R d t ) 2 + R 2 d 2 R d t 2 ) − 2 R 4 r 5 ( d R d t ) 2 ] d r {\displaystyle -{\frac {1}{\rho _{L}}}\int _{P(R)}^{P_{\infty }}\mathrm {d} P=\int _{R}^{\infty }\left[{\frac {1}{r^{2}}}\left(2R\left({\frac {\mathrm {d} R}{\mathrm {d} t}}\right)^{2}+R^{2}{\frac {\mathrm {d} ^{2}R}{\mathrm {d} t^{2}}}\right)-{\frac {2R^{4}}{r^{5}}}\left({\frac {\mathrm {d} R}{\mathrm {d} t}}\right)^{2}\right]\mathrm {d} r} を求めると P ( R ) − P ∞ ρ L = [ − 1 r ( 2 R ( d R d t ) 2 + R 2 d 2 R d t 2 ) + R 4 2 r 4 ( d R d t ) 2 ] R ∞ = R d 2 R d t 2 + 3 2 ( d R d t ) 2 {\displaystyle {\frac {P(R)-P_{\infty }}{\rho _{L}}}=\left[-{\frac {1}{r}}\left(2R\left({\frac {\mathrm {d} R}{\mathrm {d} t}}\right)^{2}+R^{2}{\frac {\mathrm {d} ^{2}R}{\mathrm {d} t^{2}}}\right)+{\frac {R^{4}}{2r^{4}}}\left({\frac {\mathrm {d} R}{\mathrm {d} t}}\right)^{2}\right]_{R}^{\infty }=R{\frac {\mathrm {d} ^{2}R}{\mathrm {d} t^{2}}}+{\frac {3}{2}}\left({\frac {\mathrm {d} R}{\mathrm {d} t}}\right)^{2}} が得られる。 液体中における気泡の中心から半径方向外向きの垂直応力を σ r r {\displaystyle \sigma _{rr}} とする。 球面座標系では、一定の密度および一定の粘度を持つ流体に対して、その垂直応力は σ r r = − P + 2 μ L ∂ u ∂ r {\displaystyle \sigma _{rr}=-P+2\mu _{L}{\frac {\partial u}{\partial r}}} と求められる。 ゆえに気泡表面の微小部分では、単位面積当たりの境界の薄い層にかかる力は σ r r ( R ) + P B − 2 S R = − P ( R ) + 2 μ L ∂ u ∂ r | r = R + P B − 2 S R = − P ( R ) + 2 μ L ∂ ∂ r ( R 2 r 2 d R d t ) r = R + P B − 2 S R = − P ( R ) − 4 μ L R d R d t + P B − 2 S R {\displaystyle {\begin{aligned}\sigma _{rr}(R)+P_{B}-{\frac {2S}{R}}&=-P(R)+2\mu _{L}\left.{\frac {\partial u}{\partial r}}\right|_{r=R}+P_{B}-{\frac {2S}{R}}\\&=-P(R)+2\mu _{L}{\frac {\partial }{\partial r}}\left({\frac {R^{2}}{r^{2}}}{\frac {\mathrm {d} R}{\mathrm {d} t}}\right)_{r=R}+P_{B}-{\frac {2S}{R}}\\&=-P(R)-{\frac {4\mu _{L}}{R}}{\frac {\mathrm {d} R}{\mathrm {d} t}}+P_{B}-{\frac {2S}{R}}\\\end{aligned}}} となる。ここで、 S {\displaystyle S} は表面張力である。 この境界で質量輸送が無いならば、面積当たりのこの力はゼロになる必要があるので、 P ( R ) = P B − 4 μ L R d R d t − 2 S R {\displaystyle P(R)=P_{B}-{\frac {4\mu _{L}}{R}}{\frac {\mathrm {d} R}{\mathrm {d} t}}-{\frac {2S}{R}}} ゆえに運動量保存の結果から P ( R ) − P ∞ ρ L = P B − P ∞ ρ L − 4 μ L ρ L R d R d t − 2 S ρ L R = R d 2 R d t 2 + 3 2 ( d R d t ) 2 {\displaystyle {\frac {P(R)-P_{\infty }}{\rho _{L}}}={\frac {P_{B}-P_{\infty }}{\rho _{L}}}-{\frac {4\mu _{L}}{\rho _{L}R}}{\frac {\mathrm {d} R}{\mathrm {d} t}}-{\frac {2S}{\rho _{L}R}}=R{\frac {\mathrm {d} ^{2}R}{\mathrm {d} t^{2}}}+{\frac {3}{2}}\left({\frac {\mathrm {d} R}{\mathrm {d} t}}\right)^{2}} ν L = μ L / ρ L {\displaystyle \nu _{L}=\mu _{L}/\rho _{L}} で置き換えることで、レイリー・プレセット方程式を得る P B ( t ) − P ∞ ( t ) ρ L = R d 2 R d t 2 + 3 2 ( d R d t ) 2 + 4 ν L R d R d t + 2 S ρ L R {\displaystyle {\frac {P_{B}(t)-P_{\infty }(t)}{\rho _{L}}}=R{\frac {\mathrm {d} ^{2}R}{\mathrm {d} t^{2}}}+{\frac {3}{2}}\left({\frac {\mathrm {d} R}{\mathrm {d} t}}\right)^{2}+{\frac {4\nu _{L}}{R}}{\frac {\mathrm {d} R}{\mathrm {d} t}}+{\frac {2S}{\rho _{L}R}}} 時間についてドット記法を用いると、レイリー・プレセット方程式はより簡潔に書ける P B ( t ) − P ∞ ( t ) ρ L = R R ¨ + 3 2 ( R ˙ ) 2 + 4 ν L R ˙ R + 2 S ρ L R {\displaystyle {\frac {P_{B}(t)-P_{\infty }(t)}{\rho _{L}}}=R{\ddot {R}}+{\frac {3}{2}}\left({\dot {R}}\right)^{2}+{\frac {4\nu _{L}{\dot {R}}}{R}}+{\frac {2S}{\rho _{L}R}}} 最近、空または気体で満たされた気泡のレイリー・プレセット方程式に対する解析的閉形式解（英語版）が見つかり、さらにN次元の場合まで一般化された。 毛細管現象による表面張力が存在する場合も研究されている。 また、表面張力と粘性が無視できる特殊な場合に対して、高次の解析的近似も知られている。 静的な場合、レイリー・プレセット方程式は単純化されてヤング・ラプラス方程式（英語版）になる: P B − P ∞ = 2 S R {\displaystyle P_{B}-P_{\infty }={\frac {2S}{R}}} 気泡の半径と圧力における微小な周期的変化のみが考慮されるときは、レイリー・プレセット方程式から気泡振動（英語版）の固有振動数が得られる。 [脚注の使い方] ^ Rayleigh, Lord (1917). “On the pressure developed in a liquid during the collapse of a spherical cavity”. Phil. Mag. 34: 94–98.  ^ Plesset, M.S. (1949). “The dynamics of cavitation bubbles”. ASME J. Appl. Mech. 16: 228–231.  ^ a b Leighton, T. G. (17 April 2007). Derivation of the Rayleigh–Plesset equation in terms of volume. Southampton, UK: Institute of Sound and Vibration Research. http://eprints.soton.ac.uk/45698/.  ^ a b Lin, Hao; Brian D. Storey; Andrew J. Szeri (2002). “Inertially driven inhomogeneities in violently collapsing bubbles: the validity of the Rayleigh–Plesset equation”. Journal of Fluid Mechanics 452. doi:10.1017/S0022112001006693. ISSN 0022-1120. http://digitalcommons.olin.edu/cgi/viewcontent.cgi?article=1002&context=facpub_2002&sei-redir=1&referer=http%3A%2F%2Fscholar.google.com%2Fscholar%3Fq%3DRayleigh–Plesset%2Bequation%26hl%3Den%26as_sdt%3D0%26as_vis%3D1%26oi%3Dscholart%26sa%3DX%26ei%3Du_2rT8q8GsqLgwe8zsnfAQ%26ved%3D0CBgQgQMwAA#search=%22Rayleigh–Plesset%20equation%22.  ^ a b c d e Brennen, Christopher E. (1995). Cavitation and Bubble Dynamics. Oxford University Press. ISBN 0-19-509409-3  ^ Kudryashov, Nikolay A.; Sinelshchikov, Dnitry I. (18 September 2014). Analytical solutions of the Rayleigh equation for empty and gas-filled bubble. arXiv:1409.6699v1. Bibcode: 2014JPhA...47N5202K. doi:10.1088/1751-8113/47/40/405202.  ^ a b Kudryashov, Nikolay A.; Sinelshchikov, Dnitry I. (31 December 2014). Analytical solutions for problems of bubble dynamics. arXiv:1608.00811. Bibcode: 2016arXiv160800811K. doi:10.1016/j.physleta.2014.12.049.  ^ Mancas, Stefan C.; Rosu, Haret C. (7 August 2015). Cavitation of spherical bubbles: closed-form, parametric, and numerical solutions. arXiv:1508.01157. Bibcode: 2016PhFl...28b2009M. doi:10.1063/1.4942237.  ^ Obreschkow, D.; Bruderer M.; Farhat, M. (5 June 2012). “Analytical approximations for the collapse of an empty spherical bubble”. Physical Review E 85. arXiv:1205.4202. Bibcode: 2012PhRvE..85f6303O. doi:10.1103/PhysRevE.85.066303.