衝撃波管シミュレーター

衝撃波の前後では、質量、運動量、エネルギーの保存則が成り立ちます。これらをまとめたランキン-ユゴニオ関係式が衝撃波管問題の核心です。

$$ \frac{p_2}{p_1}= 1 + \frac{2\gamma}{\gamma+1}(M_s^2 - 1) $$ $$ \frac{\rho_2}{\rho_1}= \frac{(\gamma+1)M_s^2}{(\gamma-1)M_s^2 + 2}$$

ここで、$p$は圧力、$\rho$は密度、$\gamma$は比熱比、$M_s$は衝撃波マッハ数です。添字1, 2はそれぞれ衝撃波の前（低圧側）と後（高圧側）を表します。圧力比はマッハ数の2乗に比例して大きくなることがわかります。

衝撃波管の全体の解は、衝撃波・接触不連続面・希薄波が満たす条件から、領域間の圧力$p^ $と速度$u^ $を求める非線形方程式を解くことで得られます。

$$ \frac{p_4}{p_1}= \frac{p_2}{p_1}\left[ 1 - \frac{(\gamma_4-1)(\frac{a_1}{a_4})(\frac{p_2}{p_1}-1)}{\sqrt{2\gamma_1}\sqrt{2\gamma_1+(\gamma_1+1)(\frac{p_2}{p_1}-1)}}\right]^{-\frac{2\gamma_4}{\gamma_4-1}} $$

これは衝撃波関係式と希薄波の等エントロピー関係を、接触不連続面での圧力・速度連続条件で結びつけた式です。$a$は音速、添字1,4は初期状態の右側（低圧）と左側（高圧）の領域を表します。この式を数値的に解くことで、すべての波の強さと位置が決まります。

CFDソルバーの検証・ベンチマーク：衝撃波管問題（ソッド問題）は、圧縮性流体を扱う数値計算コードの精度を検証するための最も基本的で重要なテストケースです。新しいCFDソルバーを開発したら、まずこの厳密解と比較し、衝撃波のシャープネスや希薄波の解像度を確認します。

超音速・極超音速空力の研究：衝撃波管は、短時間ながら高エンタルピーの気流を発生させることができます。この特性を利用し、再突入カプセルや極超音速ミサイルの先端で発生する強い衝撃波とそれに伴う高温ガス流の研究に応用されています。

爆轟・爆発現象のモデル化：爆薬の爆轟やガス爆発では、非常に強い衝撃波（爆轟波）が発生します。衝撃波管の理論は、こうした爆轟波の伝播速度、背後に誘起される流れ、そして圧力分布を予測する物理モデルの基礎となっています。

航空宇宙エンジンのインテーク設計：超音速飛行時、エンジンインテーク内部では複雑な衝撃波系が形成されます。衝撃波と境界層の干渉は流れを不安定にします。衝撃波管の研究で得られる知見は、インテーク形状の最適化やショックコントロール技術の開発に役立てられています。

まず、「衝撃波は常に超音速で伝わる」と思いがちですが、これは衝撃波に対する相対速度の話です。シミュレーター上で観測される波の伝播速度は、初期状態の音速に対するマッハ数で決まります。例えば、左右の圧力差をほんの少しだけ大きくした「弱い衝撃波」では、マッハ数は1.05など1に近い値になり、ほぼ音速で伝わります。逆に「強い衝撃波」プリセットではマッハ数が3を超えることもあり、これは爆発に近い現象です。

次に、接触不連続面の重要性を見落としがちです。この面は圧力と速度が連続なので、グラフをぱっと見ただけでは見逃してしまいます。しかし、ここは密度と温度、さらにはエントロピーが不連続になる、流体の「履歴」が刻まれた面です。例えば、初期状態で左右のガスを変える（左：ヘリウム、右：空気）と、この面の前後で密度が大きく異なるため、後続の流れや混合現象に巨大な影響を与えます。実務のCFDでは、この面の数値的拡散をいかに抑えるかが計算精度の鍵になります。

最後に、比熱比γの設定は軽視できないという点です。デフォルトの空気（γ=1.4）から、二原子ガスや単原子ガス、高温ガスに変えると、結果は想像以上に変わります。例えば、単原子ガスのアルゴン（γ=1.67）では、同じ圧力比でも衝撃波後の温度上昇が空気の場合よりも顕著になります。これは、比熱比が気体の分子自由度を反映し、内部エネルギーへの配分を決めるためです。安易にデフォルト値を使い回さず、対象とする流体の物性を確認する習慣をつけましょう。