ランキン・ユゴニオの式

ランキン・ユゴニオの式（ランキン・ユゴニオのしき、英: Rankine–Hugoniot equation）、またはランキン・ユゴニオ関係式とは、垂直衝撃波の通過前後における物理量の関係を表す次の式のことである： ρ 2 ρ 1 = u 1 u 2 = ( γ − 1 ) + ( γ + 1 ) p 2 / p 1 ( γ + 1 ) + ( γ − 1 ) p 2 / p 1 T 2 T 1 = a 2 2 a 1 2 = p 2 p 1 ⋅ ( γ + 1 ) + ( γ − 1 ) p 2 / p 1 ( γ − 1 ) + ( γ + 1 ) p 2 / p 1 {\displaystyle {\begin{aligned}&{\frac {\rho _{2}}{\rho _{1}}}={\frac {u_{1}}{u_{2}}}={\frac {(\gamma -1)+(\gamma +1)p_{2}/p_{1}}{(\gamma +1)+(\gamma -1)p_{2}/p_{1}}}\\&{\frac {T_{2}}{T_{1}}}={\frac {a_{2}^{2}}{a_{1}^{2}}}={\frac {p_{2}}{p_{1}}}\cdot {\frac {(\gamma +1)+(\gamma -1)p_{2}/p_{1}}{(\gamma -1)+(\gamma +1)p_{2}/p_{1}}}\end{aligned}}} ここで ρ：流体の密度、[kg/m3] u ：流速、[m/s] p ：圧力、[ Pa ] T ：温度、[ K ] a ：音速、[m/s] γ：比熱比 添字の1，2は衝撃波の上流、下流の意味 である。 これらの関係式は、衝撃波の前後の状態だけを、その内部構造に立ち入ることなく関係付けることができる点に特徴がある。 ウィリアム・ランキンが1870年に発表し、ピエール＝アンリ・ユゴニオがそれを知らないまま1887年にランキンと同様の結果を報告した。 ランキン・ユゴニオの式を導出するにあたっては、以下の仮定を置いている： 垂直衝撃波：平面衝撃波がその面に垂直な方向に伝播しており、流れは1次元的である。 定常：波面に固定した座標系を用い、流れは時間変化しないものとする。 衝撃波の前後はいずれも一様な状態である。 流体は理想気体であり、かつ状態変化は断熱過程とする。 上述の物理的仮定のもとで、流体の状態は以下の連続の式、運動量保存則およびエネルギー保存則によって記述される。 ∂ ρ ∂ t = − ∂ ( ρ u ) ∂ x ∂ ( ρ u ) ∂ t = − ∂ ∂ x ( ρ u 2 + p ) ∂ ( ρ E ) ∂ t = − ∂ ∂ x [ ρ u ( e + 1 2 u 2 + p ρ ) ] {\displaystyle {\begin{aligned}{\frac {\partial \rho }{\partial t}}&=-{\frac {\partial (\rho u)}{\partial x}}\\{\frac {\partial (\rho u)}{\partial t}}&=-{\frac {\partial }{\partial x}}(\rho u^{2}+p)\\{\frac {\partial (\rho E)}{\partial t}}&=-{\frac {\partial }{\partial x}}\left[\rho u\left(e+{\frac {1}{2}}u^{2}+{\frac {p}{\rho }}\right)\right]\end{aligned}}} ここで e ：比内部エネルギーもしくは比エンタルピー、[J/kg] E = e + 1 2 u 2 {\displaystyle E=e+{\frac {1}{2}}u^{2}} ：総エネルギー、[J/kg] である。さらに定常なので時間微分項は 0 になるなどの仮定を用いてこれらを積分すると、以下の式が得られる： Q ≡ ρ u = c o n s t a n t ρ u 2 + p = c o n s t a n t p u + Q ( u 2 2 + e ) = c o n s t a n t {\displaystyle {\begin{aligned}&Q\equiv \rho u={\rm {constant}}\\&\rho u^{2}+p={\rm {constant}}\\&pu+Q\left({\frac {u^{2}}{2}}+e\right)=\mathrm {constant} \end{aligned}}} 簡単のため、衝撃波は平面として、 x {\displaystyle x} 方向にのみ伝搬するものとする。衝撃波が通過する前の領域(衝撃波上流)と衝撃波が通過した後の領域(衝撃波下流)とで物理量は不連続になっており、上流側の密度、速度、単位質量あたりの内部エネルギー(specific internal energy)、圧力を ρ 0 , v 0 , ε 0 , P 0 {\displaystyle \rho _{0},v_{0},\varepsilon _{0},P_{0}} とし、下流側の密度、速度、単位質量あたりの内部エネルギー、圧力を ρ 1 , v 1 , ε 1 , P 1 {\displaystyle \rho _{1},v_{1},\varepsilon _{1},P_{1}} とする。 質量(連続の式)、運動量、エネルギーの保存則から ρ 0 v 0 = ρ 1 v 1 ρ 0 v 0 2 + P 0 = ρ 1 v 1 2 + P 1 ( ρ 0 ε 0 + 1 2 ρ 0 v 0 ) v 0 + P 0 v 0 = ( ρ 1 ε 1 + 1 2 ρ 1 v 1 ) v 1 + P 1 v 1 . {\displaystyle {\begin{aligned}\rho _{0}v_{0}&=\rho _{1}v_{1}\\\rho _{0}v_{0}^{2}+P_{0}&=\rho _{1}v_{1}^{2}+P_{1}\\\left(\rho _{0}\varepsilon _{0}+{\frac {1}{2}}\rho _{0}v_{0}\right)v_{0}+P_{0}v_{0}&=\left(\rho _{1}\varepsilon _{1}+{\frac {1}{2}}\rho _{1}v_{1}\right)v_{1}+P_{1}v_{1}.\end{aligned}}} 第１式で第３式をわると ε 0 + 1 2 v 0 2 + P 0 ρ 0 = ε 1 + 1 2 v 1 2 + P 1 ρ 1 ⇒ h 1 − h 0 = 1 2 ( v 0 2 − v 1 2 ) w h e r e h i ≡ ε i + P i ρ i {\displaystyle {\begin{aligned}&\varepsilon _{0}+{\frac {1}{2}}v_{0}^{2}+{\frac {P_{0}}{\rho _{0}}}=\varepsilon _{1}+{\frac {1}{2}}v_{1}^{2}+{\frac {P_{1}}{\rho _{1}}}\\\Rightarrow &h_{1}-h_{0}={\frac {1}{2}}(v_{0}^{2}-v_{1}^{2})\quad \mathrm {where} \quad h_{i}\equiv \varepsilon _{i}+{\frac {P_{i}}{\rho _{i}}}\end{aligned}}} ここで、 h i {\displaystyle h_{i}} は単位質量あたりのエンタルピーである。さらに第１式を ρ 1 ρ 0 = v 0 v 1 = V 0 V 1 w h e r e V i ≡ 1 ρ i {\displaystyle {\begin{aligned}{\frac {\rho _{1}}{\rho _{0}}}={\frac {v_{0}}{v_{1}}}={\frac {V_{0}}{V_{1}}}\quad \mathrm {where} \quad V_{i}\equiv {\frac {1}{\rho _{i}}}\end{aligned}}} と変形する。ここで V i {\displaystyle V_{i}} は単位質量あたりの体積である。すると第２式から ρ 0 v 0 2 + P 0 = ρ 0 2 ρ 1 v 0 2 + P 1 ⇒ v 0 2 = P 1 − P 0 ρ 0 − ρ 0 2 ρ 1 = P 1 − P 0 ρ 0 − ρ 0 2 ρ 1 = V 0 2 P 1 − P 0 V 0 − V 1 ⇒ v 1 2 = V 1 2 P 1 − P 0 V 0 − V 1 . {\displaystyle {\begin{aligned}&\rho _{0}v_{0}^{2}+P_{0}={\frac {\rho _{0}^{2}}{\rho _{1}}}v_{0}^{2}+P_{1}\\\Rightarrow &v_{0}^{2}={\frac {P_{1}-P_{0}}{\rho _{0}-{\frac {\rho _{0}^{2}}{\rho _{1}}}}}={\frac {P_{1}-P_{0}}{\rho _{0}-{\frac {\rho _{0}^{2}}{\rho _{1}}}}}=V_{0}^{2}{\frac {P_{1}-P_{0}}{V_{0}-V_{1}}}\\\Rightarrow &v_{1}^{2}=V_{1}^{2}{\frac {P_{1}-P_{0}}{V_{0}-V_{1}}}.\end{aligned}}} より v 0 − v 1 = ( 1 − V 1 V 0 ) v 0 = ( 1 − V 1 V 0 ) V 0 P 1 − P 0 V 0 − V 1 = ( P 1 − P 0 ) ( V 0 − V 1 ) v 0 2 − v 1 2 = ( P 1 − P 0 ) ( V 0 + V 1 ) {\displaystyle {\begin{aligned}v_{0}-v_{1}&=\left(1-{\frac {V_{1}}{V_{0}}}\right)v_{0}=\left(1-{\frac {V_{1}}{V_{0}}}\right)V_{0}{\sqrt {\frac {P_{1}-P_{0}}{V_{0}-V_{1}}}}={\sqrt {(P_{1}-P_{0})(V_{0}-V_{1})}}\\v_{0}^{2}-v_{1}^{2}&=(P_{1}-P_{0})(V_{0}+V_{1})\\\end{aligned}}} エンタルピーの表式に代入することで h 1 − h 0 = 1 2 ( P 1 − P 0 ) ( V 0 + V 1 ) {\displaystyle {\begin{aligned}h_{1}-h_{0}={\frac {1}{2}}(P_{1}-P_{0})(V_{0}+V_{1})\end{aligned}}} を得る。もしくは ε 1 − ε 0 = 1 2 ( v 0 2 − v 1 2 ) + V 0 P 0 − V 1 P 1 = 1 2 ( P 1 − P 0 ) ( V 0 + V 1 ) + V 0 P 0 − V 1 P 1 = 1 2 ( P 1 + P 0 ) ( V 0 − V 1 ) . {\displaystyle {\begin{aligned}\varepsilon _{1}-\varepsilon _{0}&={\frac {1}{2}}(v_{0}^{2}-v_{1}^{2})+V_{0}P_{0}-V_{1}P_{1}\\&={\frac {1}{2}}(P_{1}-P_{0})(V_{0}+V_{1})+V_{0}P_{0}-V_{1}P_{1}\\&={\frac {1}{2}}(P_{1}+P_{0})(V_{0}-V_{1}).\end{aligned}}} これをランキン・ユゴニオ関係(Rankin-Hugoniot relation)と呼ぶ。 1つの衝撃波による圧縮の限界を調べる。理想気体の場合、状態方程式 P = ( γ − 1 ) ρ ε {\displaystyle P=(\gamma -1)\rho \varepsilon } を仮定すれば P 1 V 1 γ − 1 − P 0 V 0 γ − 1 = 1 2 ( P 1 + P 0 ) ( V 0 − V 1 ) ⇒ V 1 V 0 = ( γ + 1 ) P 0 + ( γ − 1 ) P 1 ( γ − 1 ) P 0 + ( γ + 1 ) P 1 → P 1 / P 0 → ∞ γ − 1 γ + 1 . {\displaystyle {\begin{aligned}&{\frac {P_{1}V_{1}}{\gamma -1}}-{\frac {P_{0}V_{0}}{\gamma -1}}={\frac {1}{2}}(P_{1}+P_{0})(V_{0}-V_{1})\\\Rightarrow &{\frac {V_{1}}{V_{0}}}={\frac {(\gamma +1)P_{0}+(\gamma -1)P_{1}}{(\gamma -1)P_{0}+(\gamma +1)P_{1}}}\xrightarrow {P_{1}/P_{0}\rightarrow \infty } {\frac {\gamma -1}{\gamma +1}}.\end{aligned}}} つまり、輻射優勢期 γ = 4 / 3 {\displaystyle \gamma =4/3} を考えれば元々の体積の1/7まで圧縮される。 磁場があるときの運動量保存則は次のようにかける： ρ D v D t = − ∇ ( p + 1 8 π B 2 ) ⏟ pressure + 1 4 π ( B ⋅ ∇ ) B ⏟ EM − ρ ∇ Φ g r a v ⏟ gravitation + x ^ i ∂ j σ i j ⏟ v i s c o s i t y . {\displaystyle {\begin{aligned}\rho {\frac {D{\boldsymbol {v}}}{Dt}}=-\underbrace {\nabla \left(p+{\frac {1}{8\pi }}B^{2}\right)} _{\text{pressure}}+\underbrace {{\frac {1}{4\pi }}({\boldsymbol {B}}\cdot \nabla ){\boldsymbol {B}}} _{\text{EM}}-\underbrace {\rho \nabla \Phi _{\mathrm {grav} }} _{\text{gravitation}}+\underbrace {{\hat {x}}_{i}\partial _{j}\sigma _{ij}} _{viscosity}.\end{aligned}}} これを x {\displaystyle x} 方向について書き下し、 ρ v x = c o n s t {\displaystyle \rho v_{x}=const} を用いると ρ v x ∂ x v x = − ( p + 1 8 π B 2 ) + B x 4 π ∂ x B x − ρ ∂ x Φ g r a v + ∂ x σ x x ⇒ ∂ x ( ρ v x 2 + p + B y 2 + B z 2 8 π − σ x x ) = − ρ ∂ x Φ g r a v ⇒ ρ 1 v 1 2 + p 1 + B 1 y 2 + B 1 z 2 8 π = ρ 2 v 2 2 + p 2 + B 2 y 2 + B 2 z 2 8 π . {\displaystyle {\begin{aligned}&\rho v_{x}\partial _{x}v_{x}=-\left(p+{\frac {1}{8\pi }}B^{2}\right)+{\frac {B_{x}}{4\pi }}\partial _{x}B_{x}-\rho \partial _{x}\Phi _{\mathrm {grav} }+\partial _{x}\sigma _{xx}\\\Rightarrow &\partial _{x}\left(\rho v_{x}^{2}+p+{\frac {B_{y}^{2}+B_{z}^{2}}{8\pi }}-\sigma _{xx}\right)=-\rho \partial _{x}\Phi _{\mathrm {grav} }\\\Rightarrow &\rho _{1}v_{1}^{2}+p_{1}+{\frac {B_{1y}^{2}+B_{1z}^{2}}{8\pi }}=\rho _{2}v_{2}^{2}+p_{2}+{\frac {B_{2y}^{2}+B_{2z}^{2}}{8\pi }}.\end{aligned}}} ここで粘性ストレステンソル σ x x {\displaystyle \sigma _{xx}} が衝撃波上流と下流で、波面に非常に近い領域でない限り、０に限りなく近いことを利用して項を落としている。また、重力項も無視している。 一方、磁場があるときのエネルギー保存則は D D t ( 1 2 ρ v 2 + ε + B 2 8 π ) = − ( p + B 2 8 π ) ( ∇ ⋅ v ) + ( ∇ ⋅ B ) ( v ⋅ ) B 4 π + ∇ ⋅ ( v × σ ) − ρ v ⋅ ∇ Φ g r a v + Γ − Λ − ∇ ⋅ ( − κ ∇ T ) . {\displaystyle {\begin{aligned}{\frac {D}{Dt}}\left({\frac {1}{2}}\rho v^{2}+\varepsilon +{\frac {B^{2}}{8\pi }}\right)=-\left(p+{\frac {B^{2}}{8\pi }}\right)(\nabla \cdot {\boldsymbol {v}})&+{\frac {(\nabla \cdot {\boldsymbol {B}})({\boldsymbol {v}}\cdot ){\boldsymbol {B}}}{4\pi }}+\nabla \cdot ({\boldsymbol {v}}\times {\boldsymbol {\sigma }})\\&-\rho {\boldsymbol {v}}\cdot \nabla \Phi _{\mathrm {grav} }+\Gamma -\Lambda -\nabla \cdot (-\kappa \nabla T).\end{aligned}}} 同様に x {\displaystyle x} 方向について書き下し、 ∂ x ( 1 2 ρ v x v 2 + ε v x + p v x + B y 2 + B z 2 4 π v x − B x B y v y 4 π − B x B z v z 4 π − v j σ j x − κ d T d x + ρ v x Φ g r a v ) = Γ − Λ . {\displaystyle {\begin{aligned}\partial _{x}\left({\frac {1}{2}}\rho v_{x}v^{2}+\varepsilon v_{x}+pv_{x}+{\frac {B_{y}^{2}+B_{z}^{2}}{4\pi }}v_{x}-{\frac {B_{x}B_{y}v_{y}}{4\pi }}-{\frac {B_{x}B_{z}v_{z}}{4\pi }}-v_{j}\sigma _{jx}-\kappa {\frac {dT}{dx}}+\rho v_{x}\Phi _{\mathrm {grav} }\right)=\Gamma -\Lambda .\end{aligned}}} ∫ ( Γ − A ) d x ≈ 0 {\displaystyle \int (\Gamma -A)dx\approx 0} 、重力項を無視して ( 1 2 ρ 1 x v 1 2 + ε 1 + p 1 + B 1 y 2 + B 1 z 2 4 π ) v 1 x − B 1 x B 1 y v 1 y − B 1 x B 1 z v 1 z 4 π − κ 1 d T 1 d x = ( 1 2 ρ 2 x v 2 2 + ε 2 + p 2 + B 2 y 2 + B 2 z 2 4 π ) v 2 x − B 2 x B 2 y v 2 y − B 2 x B 2 z v 2 z 4 π − κ 2 d T 2 d x {\displaystyle {\begin{aligned}\left({\frac {1}{2}}\rho _{1x}v_{1}^{2}+\varepsilon _{1}\right.&\left.+p_{1}+{\frac {B_{1y}^{2}+B_{1z}^{2}}{4\pi }}\right)v_{1x}-{\frac {B_{1x}B_{1y}v_{1y}-B_{1x}B_{1z}v_{1z}}{4\pi }}-\kappa _{1}{\frac {dT_{1}}{dx}}\\&=\left({\frac {1}{2}}\rho _{2x}v_{2}^{2}+\varepsilon _{2}+p_{2}+{\frac {B_{2y}^{2}+B_{2z}^{2}}{4\pi }}\right)v_{2x}-{\frac {B_{2x}B_{2y}v_{2y}-B_{2x}B_{2z}v_{2z}}{4\pi }}-\kappa _{2}{\frac {dT_{2}}{dx}}\end{aligned}}} を得る。加えて、磁束保存 ∂ ∂ t B = ∇ × ( v × B ) {\displaystyle {\begin{aligned}{\frac {\partial }{\partial t}}{\boldsymbol {B}}=\nabla \times ({\boldsymbol {v}}\times {\boldsymbol {B}})\end{aligned}}} について、 ∂ t ∂ y = ∂ z = 0 {\displaystyle \partial _{t}\partial _{y}=\partial _{z}=0} より y , z {\displaystyle y,z} 成分について 0 = − ∂ x ( v x B y − v y B z ) 0 = − ∂ x ( v z B x − v x B z ) . {\displaystyle {\begin{aligned}0&=-\partial _{x}(v_{x}B_{y}-v_{y}B_{z})\\0&=-\partial _{x}(v_{z}B_{x}-v_{x}B_{z}).\end{aligned}}} さて、前節と同様の状況を考える。 v 1 x = v 1 , v 2 x = v 2 {\displaystyle v_{1x}=v_{1},v_{2x}=v_{2}} 、それ以外の速度は0とする。質量、運動量、エネルギー、磁束の保存を書き下せば ρ 1 v 1 = ρ 2 v 2 ρ 1 v 1 2 + p 1 + B 1 2 8 π = ρ 2 v 2 2 + p 2 + B 2 2 8 π ( 1 2 ρ 1 v 1 2 + ε 1 + p 1 + B 1 2 4 π ) v 1 = ( 1 2 ρ 2 v 2 2 + ε 2 + p 2 + B 2 2 4 π ) v 2 v 1 B 1 = v 2 B 2 . {\displaystyle {\begin{aligned}\rho _{1}v_{1}&=\rho _{2}v_{2}\\\rho _{1}v_{1}^{2}+p_{1}+{\frac {B_{1}^{2}}{8\pi }}&=\rho _{2}v_{2}^{2}+p_{2}+{\frac {B_{2}^{2}}{8\pi }}\\\left({\frac {1}{2}}\rho _{1}v_{1}^{2}+\varepsilon _{1}+p_{1}+{\frac {B_{1}^{2}}{4\pi }}\right)v_{1}&=\left({\frac {1}{2}}\rho _{2}v_{2}^{2}+\varepsilon _{2}+p_{2}+{\frac {B_{2}^{2}}{4\pi }}\right)v_{2}\\v_{1}B_{1}&=v_{2}B_{2}.\end{aligned}}} 理想気体の状態方程式 p = ( γ − 1 ) ε {\displaystyle p=(\gamma -1)\varepsilon } を同様に仮定する。同様に計算を頑張ると v 1 v 2 = 2 ( γ + 1 ) D + D 2 + 4 ( γ + 1 ) ( 2 − γ ) M A − 2 , w h e r e D ≡ γ − 1 + 2 M − 2 + γ M A − 2 , M ≡ v 1 γ p 1 V 1 , M A ≡ v 1 B 1 4 π V 1 . {\displaystyle {\begin{aligned}{\frac {v_{1}}{v_{2}}}={\frac {2(\gamma +1)}{D+{\sqrt {D^{2}+4(\gamma +1)(2-\gamma )M_{A}^{-2}}}}},\quad \mathrm {where} \quad D\equiv \gamma -1+2M^{-2}+\gamma M_{A}^{-2},\quad M\equiv {\frac {v_{1}}{\sqrt {\gamma p_{1}V_{1}}}},\quad M_{A}\equiv {\frac {v_{1}}{B_{1}}}{\sqrt {\frac {4\pi }{V_{1}}}}.\end{aligned}}} これをランキン・ユゴニオ関係(Rankin-Hugoniot relation)と呼ぶ。 v 1 {\displaystyle v_{1}} の速度が磁気音速(magnetrosonic speed) V m s {\displaystyle V_{\mathrm {ms} }} に比べ十分大きいとき( v 1 ≫ V m s {\displaystyle v_{1}\gg V_{\mathrm {ms} }} :強衝撃波極限)、 v 1 v 2 → V 1 V 2 → γ + 1 γ − 1 v 2 → γ − 1 γ + 1 v 1 T 2 = p 2 ρ 2 k B → 2 ( γ − 1 ) ( γ + 1 ) 2 μ v s 2 k B . {\displaystyle {\begin{aligned}{\frac {v_{1}}{v_{2}}}&\rightarrow {\frac {V_{1}}{V_{2}}}\rightarrow {\frac {\gamma +1}{\gamma -1}}\\v_{2}&\rightarrow {\frac {\gamma -1}{\gamma +1}}v_{1}\\T_{2}&={\frac {p_{2}}{\rho _{2}k_{B}}}\rightarrow {\frac {2(\gamma -1)}{(\gamma +1)^{2}}}{\frac {\mu v_{s}^{2}}{k_{B}}}.\end{aligned}}} ^ 巽友正『流体力学』（1版）培風館、1982年、224頁。ISBN 4-563-02421-X。  ^ 永田雅人『高速流体力学』森北出版、2010年、5, 74頁。ISBN 978-4-627-67361-8。  ^ Draine, Bruce T. (2011). Physics of the Interstellar and Intergalactic Medium. Princeton University Press. https://ui.adsabs.harvard.edu/abs/2011piim.book.....D  流体力学 空気力学 ソニックブーム 熱力学 CJ理論 ZND理論 燃焼工学