ランキン・ユゴニオの式

**ランキン・ユゴニオの式**（ランキン・ユゴニオのしき、）、または**ランキン・ユゴニオ関係式**とは、垂直衝撃波の通過前後における物理量の関係を表す次の式のことである：
: \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}
ここで
: *&rho;*：流体の密度、[kg/m3]
: *u* ：流速、[m/s]
: *p* ：圧力、[ Pa ]
: *T* ：温度、[ K ]
: *a* ：音速、[m/s]
: *&gamma;*：比熱比
: 添字の1，2は衝撃波の上流、下流の意味
である。

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

ウィリアム・ランキンが1870年に発表し、ピエール＝アンリ・ユゴニオがそれを知らないまま1887年にランキンと同様の結果を報告した。

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

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

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

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

\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}

第１式で第３式をわると

\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} 

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

\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} 

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

\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} 

より

\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} 

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

\begin{align}

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

\end{align} 

を得る。もしくは

\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} 

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

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

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

### 磁場ありのランキン・ユゴニオの関係

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

\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}

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

\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}

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

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


\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}

同様にx方向について書き下し、

\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}

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


\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}

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


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

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


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

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


\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}

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


\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}

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


\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}


## 参考文献

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