超音速流れ

最も根本的な違いは情報伝播の方向だ。亜音速では圧力擾乱が全方向に伝わるから、下流の障害物の情報が上流にも伝わる。しかし超音速では流れが音速を超えているため、擾乱は下流にしか伝播しない。この性質が衝撃波の形成、マッハ円錐の存在、そして支配方程式の型の変化（楕円型から双曲型へ）を生み出す。

マッハ数 $M = U/a$ が1を超えると、擾乱はマッハ角 $\mu$ の円錐内にしか伝わらない。

直接関係するよ。くさび型物体の先端に形成される斜め衝撃波の角度 $\beta$ は、くさびの半角 $\theta$ とマッハ数 $M$ の関係式（$\theta$-$\beta$-$M$ 関係）で決まる。

超音速流れが凸角で曲がるとき、連続的な膨張波（Prandtl-Meyer expansion fan）が形成される。膨張波を通過すると流れは加速し、マッハ数が上昇する。転向角 $\Delta\theta$ とマッハ数の関係はPrandtl-Meyer関数 $\nu(M)$ で記述される。

数値的にテーブルを作るか、Newton法で逆解きする。CFDでは当然ソルバーが自動的に膨張波を解像してくれるけど、結果の検証には不可欠な式だよ。超音速ノズル設計のMOC（Method of Characteristics、特性線法）でも中心的な役割を果たす。

ダイヤモンド翼型は超音速線形理論の代表例だ。前縁と後縁にそれぞれ衝撃波と膨張波が形成され、上下面の圧力差から揚力が、圧力の流れ方向成分から造波抗力が生じる。超音速線形理論（Ackeret理論）による圧力係数は、

主な違いは3つある。第一に衝撃波捕捉のための風上系スキームが必要。第二に化学反応（高温ガス効果）の考慮が高マッハ数で必要になる。第三に境界条件の設定が異なる。超音速入口では全変数を規定し、超音速出口では全変数を外挿する。

Shock fitting（衝撃波適合法）は衝撃波の位置を未知数として扱い、衝撃波を格子の境界に正確に配置する方法だ。衝撃波が鋭く解像されるメリットがある。一方、shock capturing（衝撃波捕捉法）は衝撃波を数値拡散の範囲で自動的に捕捉する。

実務ではshock capturingが圧倒的に主流だ。理由は以下の通り。

MOCは超音速ノズルの輪郭設計では今でも標準的な手法だよ。定常2次元超音速流れでは特性線に沿ってリーマン不変量が保存されるから、特性線を追跡していけば流れ場を構築できる。

具体的には、$C^+$ 特性線（マッハ線）に沿って、

MOCは無粘性・等エントロピーの仮定のもとで厳密解を与えるから、CFDの検証に使える。またノズル設計では逆問題（所望のマッハ数分布→壁面形状）を直接解けるのがCFDにはない強みだ。

$M > 5$ 程度（極超音速域）では衝撃波後方の温度が数千Kに達し、以下の効果が無視できなくなる。

これらを扱うには、理想気体の状態方程式 $p = \rho R T$ を化学平衡モデルまたは化学非平衡モデルに置き換える必要がある。FluentのSpecies Transportモデル、OpenFOAMのhy2Foamソルバー、NASA CEAデータベースを用いたテーブル参照法などが使われるよ。

典型的なワークフローを示すよ。

Step 1: 計算領域設計

マッハ数、静圧、静温度、そして流れ方向のベクトル成分を指定する。この境界条件はRiemann不変量に基づいて、局所マッハ数に応じて流入と流出を自動判定してくれるんだ。超音速外部流れでは最も使い勝手がいい境界条件だよ。

インテークでは以下が重要だ。

出口断面の質量流量加重平均全圧を入口全圧で割る。

以下が超音速外部流れの標準設定だ。

衝撃波が弱い場合や格子品質が良ければCFL=50以上も可能だ。ただしSBLIがある場合や格子品質が悪い場合はCFL=5-10程度に留めた方が安全だよ。残差モニターを見ながら段階的に上げるのが定石だ。

用途で使い分ける。

Krassilnikovらが開発したもので、GitHub上で公開されている。低マッハ数域ではPISO的に振る舞い、高マッハ数域ではKT的に振る舞うスイッチング機構を持つ。遷音速ノズルのように局所的にマッハ数が大きく変わる問題に向いているよ。

15度半角のくさびに $M = 3.0$ の流れを当てる問題が最も基本的だ。斜め衝撃波理論から $\beta = 32.2°$、衝撃波後方のマッハ数 $M_2 = 2.25$、圧力比 $p_2/p_1 = 2.82$ と解析的に求まる。

CFD結果との比較ポイントは以下の通り。

その通り。次のステップとしてはダイヤモンド翼型の揚力・抗力計算、双円錐物体（double cone）の衝撃波干渉、超音速ノズル内の流れなどが定番の検証問題だよ。

超音速燃焼の最大の難しさは、燃焼と流れの相互作用が非常に強いことだ。燃焼による熱解放が衝撃波構造を変え、衝撃波が燃焼を促進する。このフィードバックループを正確に捉えるには、圧縮性流体力学と化学反応速度論の両方が高精度で必要になる。

主要な物理モデルは以下の通りだ。

水素/空気なら9化学種・19反応（Mueller et al.）が定番。炭化水素燃料だと詳細メカニズムは数千反応になるから、Skeletal mechanismやReduced mechanismに縮約して使う。FluentではEDC+有限速度化学の組み合わせが使いやすい。

ジェット騒音のCFD予測は活発な研究分野だ。LESで近傍場の音源を解像し、FWH（Ffowcs Williams-Hawkings）面積分で遠方場の音圧を計算する手法が標準的だよ。

超音速ジェット特有の騒音メカニズムとして、

単一円形ノズルの超音速ジェット騒音LESで典型的に50-200Mセル、シミュレーション時間は100-500個の非定常周期分が必要。HPC上で数十万コア時間のオーダーだ。NASAやJAXAが精力的に取り組んでいる分野だよ。

ソニックブーム低減のための形状最適化は超音速旅客機開発の核心技術だ。アジョイント法（adjoint method）を使えば、設計変数の数に依存しない効率でコスト関数の勾配を計算できる。

FluentのAdjoint Solverは超音速にも対応していて、抗力最小化やソニックブームシグネチャの最適化に使える。SU2もアジョイントベースの形状最適化機能を標準搭載しているよ。NASAのFun3Dコードは超音速ソニックブーム最適化の研究で実績がある。

超音速特有の発散原因をまとめよう。

超音速出口に静圧を固定すると、出口から非物理的な情報が上流に伝わり、数値的に不安定になる。超音速出口は全変数が外挿されるべきだから、FluentではSupersonic/Initial Gauge Pressureの「Supersonic」設定を使う。出口にごく一部でも亜音速域がある場合（例えば壁面境界層内）は、Pressure Outletに小さなbackflow圧力を設定するのが安全だ。

まず確認すべきは以下だ。

それは数値拡散による数値的なエントロピー生成だ。風上スキームは本質的に人工散逸を含むので、衝撃波がない領域でもわずかなエントロピー生成が起こる。対策としては、

等エントロピー過程（膨張波、滑らかな加速）での数値的全圧損失が0.1%以下であれば十分な精度だ。衝撃波を含む場合は理論的なRankine-Hugoniot全圧損失と比較して、数値的な追加損失が数%以下であれば実用的な精度と言えるよ。格子収束性の確認は必須だ。