ランキン・ユゴニオの式のソースを表示

'''ランキン・ユゴニオの式'''（ランキン・ユゴニオのしき、{{lang-en-short|Rankine–Hugoniot equation}}）、または'''ランキン・ユゴニオ関係式'''とは、垂直[[衝撃波]]の通過前後における[[物理量]]の関係を表す次の式のことである<ref name="tatsumi">{{cite|和書|page=224||author=巽友正|title=流体力学|edition=1|publisher=培風館|isbn=4-563-02421-X|year=1982}}</ref>：
: <math>\begin{align}
& \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{align}</math>
ここで
: ''&rho;''：流体の[[密度]]、[kg/m<sup>3</sup>]
: ''u'' ：流速、[m/s]
: ''p'' ：[[圧力]]、[ [[パスカル (単位)|Pa]] ]
: ''T'' ：[[温度]]、[ [[ケルビン|K]] ]
: ''a'' ：[[音速]]、[m/s]
: ''&gamma;''：[[比熱比]]
: 添字の1，2は衝撃波の上流、下流の意味
である。

これらの関係式は、衝撃波の前後の状態だけを、その内部構造に立ち入ることなく関係付けることができる点に特徴がある。

[[ウィリアム・ランキン]]が1870年に発表し、[[ピエール＝アンリ・ユゴニオ]]がそれを知らないまま1887年<!--巽の本には1889年という記述もあるが、本人死亡後？-->にランキンと同様の結果を報告した<ref name="nagata">{{cite|和書 |author=永田雅人 |title=高速流体力学 |publisher=森北出版 |year=2010 |isbn=978-4-627-67361-8 |page=5, 74}}</ref>。

== 物理的仮定 ==
ランキン・ユゴニオの式を導出するにあたっては、以下の仮定を置いている：
* 垂直衝撃波：平面衝撃波がその面に垂直な方向に伝播しており、流れは1次元的である。
* [[定常]]：波面に固定した座標系を用い、流れは時間変化しないものとする。
* 衝撃波の前後はいずれも[[一様]]な状態である。
* 流体は[[理想気体]]であり、かつ状態変化は[[断熱過程]]とする。

== 導出 ==
=== 磁場なしのランキン・ユゴニオの関係 ===
上述の物理的仮定のもとで、流体の状態は以下の[[連続の式]]、[[運動量保存則]]および[[エネルギー保存則]]によって記述される。
: <math>\begin{align}
  \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{align}</math>

ここで
: ''e'' ：比[[内部エネルギー]]もしくは比[[エンタルピー]]、[J/kg]
: <math>E = e+\frac{1}{2}u^2</math> ：総エネルギー、[J/kg]
である。さらに定常なので[[時間微分]]項は 0 になるなどの仮定を用いてこれらを積分すると、以下の式が得られる：
: <math>\begin{align}
& Q \equiv \rho u = \rm{constant} \\
& \rho u^2 + p = \rm{constant} \\
& pu + Q\left(\frac{u^2}{2}+e\right) = \mathrm{constant}
\end{align}</math>
簡単のため、衝撃波は平面として、<math>x</math>方向にのみ伝搬するものとする。衝撃波が通過する前の領域(衝撃波上流)と衝撃波が通過した後の領域(衝撃波下流)とで物理量は不連続になっており、上流側の密度、速度、単位質量あたりの内部エネルギー(specific internal energy)、圧力を<math>\rho_0, v_0,\varepsilon_0, P_0</math>とし、下流側の密度、速度、単位質量あたりの内部エネルギー、圧力を<math>\rho_1, v_1,\varepsilon_1,P_1 </math>とする。

質量(連続の式)、運動量、エネルギーの保存則から

<math>\begin{align}
\rho_0v_0 &= \rho_1v_1\\
\rho_0v_0^2+P_0 &= \rho_1v_1^2 + P_1\\
\left(\rho_0\varepsilon_0+\frac{1}{2}\rho_0v_0\right)v_0 + P_0v_0 &= \left(\rho_1\varepsilon_1+\frac{1}{2}\rho_1v_1\right)v_1+P_1v_1.
\end{align}</math>

第１式で第３式をわると

<math>\begin{align}

&\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{align} </math>


ここで、<math>h_i </math>は単位質量あたりのエンタルピーである。さらに第１式を

<math>\begin{align}

 \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{align} </math>

と変形する。ここで<math>V_i </math>は単位質量あたりの体積である。すると第２式から

<math>\begin{align}

&\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{align} </math>

より

<math>\begin{align}
 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{align} </math>

エンタルピーの表式に代入することで

<math>\begin{align}

h_1 -h_0 = \frac{1}{2}(P_1-P_0)(V_0+V_1)

\end{align} </math>

を得る。もしくは

<math>\begin{align}

\varepsilon_1 -\varepsilon_0 &= \frac{1}{2}(v_0^2-v_1^2)+V_0P_0 -V_1P_1 \\

&=\frac{1}{2}(P_1-P_0)(V_0+V_1)+V_0P_0 -V_1P_1\\

&=\frac{1}{2}(P_1+P_0)(V_0-V_1).

\end{align} </math>

これを'''ランキン・ユゴニオ関係(Rankin-Hugoniot relation)'''と呼ぶ。

1つの衝撃波による圧縮の限界を調べる。理想気体の場合、状態方程式<math>P = (\gamma-1)\rho\varepsilon </math>
を仮定すれば
 <math>
\begin{align}
     &\frac{P_1V_1}{\gamma-1} - \frac{P_0V_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{align}
 </math>

つまり、輻射優勢期 <math>\gamma=4/3</math>を考えれば元々の体積の1/7まで圧縮される。

=== 磁場ありのランキン・ユゴニオの関係<ref>{{Cite book|洋書 |title=Physics of the Interstellar and Intergalactic Medium |year=2011 |publisher=Princeton University Press |url=https://ui.adsabs.harvard.edu/abs/2011piim.book.....D |author=Draine, Bruce T.}}</ref> ===

磁場があるときの運動量保存則は次のようにかける：

<math>\begin{align}
    \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{align}</math>

これを<math>x</math>方向について書き下し、<math>\rho v_x=const</math>を用いると

<math>\begin{align}
    &\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_1v_1^2 + p_1+ \frac{B_{1y}^2+B_{1z}^2}{8\pi} = \rho_2v_2^2 + p_2+ \frac{B_{2y}^2+B_{2z}^2}{8\pi}.
\end{align}</math>

ここで粘性ストレステンソル<math>\sigma_{xx}</math>が衝撃波上流と下流で、波面に非常に近い領域でない限り、０に限りなく近いことを利用して項を落としている。また、重力項も無視している。

一方、磁場があるときのエネルギー保存則は

<math>
\begin{align}
    \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{align}</math>

同様に<math>x</math>方向について書き下し、

<math>\begin{align}
    \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_xB_yv_y}{4\pi} - \frac{B_xB_zv_z}{4\pi}- v_j\sigma_{jx} -\kappa \frac{dT}{dx}+\rho v_x \Phi_{\mathrm{grav}}
    \right)  = \Gamma -\Lambda.
\end{align}</math>

<math>\int(\Gamma-A)dx \approx0</math>、重力項を無視して

<math>
\begin{align}
    \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{align}</math>

を得る。加えて、磁束保存

<math>
\begin{align}
    \frac{\partial }{\partial t}\boldsymbol{B} = \nabla\times (\boldsymbol{v}\times \boldsymbol{B})
\end{align}</math>

について、<math>\partial_t \partial_y =\partial_z =0</math>より<math>y,z</math>成分について

<math>
\begin{align}
    0 &= - \partial_x(v_xB_y - v_yB_z)\\
    0 &= - \partial_x(v_zB_x - v_xB_z).
\end{align}</math>

さて、前節と同様の状況を考える。<math>v_{1x} = v_1, v_{2x}=v_2</math>、それ以外の速度は0とする。質量、運動量、エネルギー、磁束の保存を書き下せば

<math>
\begin{align}
\rho_1v_1 &= \rho_2v_2\\
\rho_1v_1^2 + p_1 + \frac{B_1^2}{8\pi} &= \rho_2v_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_1B_1 &= v_2B_2.
\end{align}</math>

理想気体の状態方程式<math>p = (\gamma-1)\varepsilon</math>を同様に仮定する。同様に計算を頑張ると

<math>
\begin{align}
    \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_1V_1}},\quad M_A \equiv \frac{v_1}{B_1}\sqrt{\frac{4\pi}{V_1}}.
\end{align}</math>

これを'''ランキン・ユゴニオ関係(Rankin-Hugoniot relation)'''と呼ぶ。
<math>v_1</math>の速度が磁気音速(magnetrosonic speed)<math>V_{\mathrm{ms}}</math>に比べ十分大きいとき(<math>v_1\gg V_{\mathrm{ms}}</math>:強衝撃波極限)、

<math>
\begin{align}
    \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_2k_B} \rightarrow \frac{2(\gamma-1)}{(\gamma+1)^2}\frac{\mu v_s^2}{k_B}.
\end{align}
</math>

== 参考文献 ==
{{reflist}}

== 関連項目 ==
* [[流体力学]]
* [[空気力学]]
* [[ソニックブーム]]
* [[熱力学]]
* [[CJ理論]]
* [[ZND理論]]
* [[燃焼工学]]

{{デフォルトソート:らんきんゆこにおのしき}}
[[Category:超音速]]
[[Category:流体力学の方程式]]
[[Category:人名を冠した数式]]