極超音速流れ

決定的に違うのは、衝撃波背後の温度が数千Kに達して気体の化学的性質が変わる点だ。空気が解離・電離を起こすから、理想気体の仮定が完全に崩れる。再突入カプセルやスクラムジェットの設計では避けて通れない領域だね。

そう。N₂やO₂が原子に解離し、さらに高温ではイオン化する。このため熱化学非平衡を考慮した支配方程式が必要になるんだ。

まず垂直衝撃波の基本、Rankine-Hugoniot関係式。理想気体の場合、衝撃波前後の圧力比は

理想気体で計算するとT₂/T₁は数百倍になるが、実際には解離反応がエネルギーを吸収するため温度上昇が抑えられる。これがまさに実在気体効果だ。Newtonの修正圧力係数も極超音速では重要で

5種（N₂, O₂, NO, N, O）または11種（イオン種と電子を含む）の化学種輸送方程式を追加する。各化学種 $Y_s$ の輸送方程式は

その通り。振動緩和時間はMillikan-Whiteの相関式で推定する。この非平衡領域が衝撃波のスタンドオフ距離や壁面加熱率に大きく影響するから、CFDでの適切なモデリングが不可欠なんだ。

淀み点の加熱率はFay-Riddellの式で概算できる。

そう。これが極超音速空力設計の基本思想だ。ブラントボディパラドックスとも呼ばれる。

基本的な有限体積法（FVM）のフレームワークは同じだが、フラックス計算スキームの選択が極めて重要になる。強い衝撃波を安定に捕獲するには、Roe法やAUSM+系のスキームにリミッター（Van Leer, Barth-Jespersen等）を組み合わせるのが標準的だ。

Advection Upstream Splitting Method。質量フラックスと圧力フラックスを分離して評価する手法で、極超音速域でのカーバンクル現象を低減できるんだ。特にAUSM+-upスキームは全速度域で安定性が高い。

その通り。衝撃波捕獲法（shock capturing）では、数値粘性で不連続面を数セル幅に滑らかにする。精度を上げるにはMUSCL再構築で2次精度にしつつ、TVD（Total Variation Diminishing）リミッターで振動を抑制する。

非常にシビア。M>10の強い弓状衝撃波では、Roe法単体だとカーバンクル不安定が発生しやすい。H-CUSP法やRotated Roeスキームなど、多次元効果を考慮した手法が有効だ。

定常解を求めるなら陰解法（LU-SGS, DPLR型point-implicit等）が効率的だ。CFL数を100以上に取れるから収束が速い。一方、非定常現象（カプセル振動等）を追うなら2次精度陰解法（BDF2）か陽解法の二段Runge-Kuttaを使う。ただし陽解法のCFL制限は

化学反応の時定数が流体の時定数より何桁も小さいため、ソース項が「stiff」になるんだ。これを解決するにはポイント陰解法で化学種ソース項のヤコビアンを陰的に扱うか、operator splitting法で流体と化学を分離して解く。NASAのDPLRコードやUS3Dコードではpoint-implicit法が標準的に使われている。

Data Parallel Line Relaxation。NASA Amesが開発した極超音速CFDコードで、Parkの熱化学モデルを内蔵している。再突入体解析の業界標準の一つだ。他にもLAURA（NASA Langley）、US3D（ミネソタ大）が有名だね。

Ansys Fluentも実在気体と有限速度化学反応に対応しているが、極超音速専用コードほど検証実績は多くない。STAR-CCM+もReacting Flowモデルがあるが、5種・11種空気モデルの実装はユーザー側でカスタマイズが必要な場合がある。

典型的なワークフローを示そう。

陰解法でもいきなり大きなCFL数で計算すると初期の強い非線形性で発散する。CFL=0.1から始めて、残差が下がるにつれて100や1000まで上げていく手法だ。Fluentでは「Courant Number」を自動調整するオプションがある。

いくつか重要なポイントがある。

事前にNewton法やMOC（特性曲線法）で衝撃波形状を推定し、格子の外側境界をその形状に合わせる。NASAのGRIDGEN2DやPointwiseの楕円型格子生成機能が使える。非構造格子の場合はAMR（Adaptive Mesh Refinement）で衝撃波を検出して自動細分化するアプローチもある。

極超音速では壁面近傍で解離した原子が壁面上で再結合するかどうかが加熱率に大きく影響する。完全触媒壁（FCW）では全ての原子が壁面で再結合してエネルギーを放出するから加熱率が最大になる。非触媒壁（NCW）では再結合しないから加熱率が低い。実際のTPS材料はこの中間で、触媒効率 $\gamma_{cat}$ をパラメータとして与える。

保守的な設計ではFCWを使う。実際のSiC表面やRCC（強化炭素複合材）の触媒効率はNASA TPで実測値が公開されている。CFDの検証では両方で計算して感度を確認するのが良い実務だ。

代表的なベンチマークを挙げよう。

その通り。極超音速CFDは物理モデルの不確かさが大きいから、V&V（Verification & Validation）が特に重要だ。

極超音速は化学反応・高温気体を扱うため、一般的なCFDソルバーでは不十分な場合がある。専用コードと汎用コードを比較してみよう。

DPLRやLAURAはNASA内部コードで一般公開されていない。ただしUS3Dは大学ライセンスがあり、hy2FoamはGitHubで公開されている。商用ソルバーで極超音速をやるなら、Fluentの実在気体モデル＋UDF（User Defined Function）でParkモデルを実装するケースが多い。

主なポイントを整理しよう。

標準機能では1温度モデルのみ。2温度モデルが必要ならUDFでエネルギー方程式を追加する必要がある。これはかなりの開発工数になるから、hy2Foamの方が手軽な場合もある。

ストラスクライド大学のVincent Casseau博士が開発したOpenFOAMベースのソルバーで、Parkの2温度モデルや多種空気モデルが実装済みだ。FIRE IIやRAM-Cの検証論文も出ている。研究用途なら十分実用的だし、ソースコードを読んで熱化学モデルの実装を学ぶ教材としても優秀だ。

ただしGUIはないから、全てテキストファイルベースの設定になる。また並列計算のスケーラビリティやメッシュサイズの実績は商用コードに及ばない点は理解しておく必要がある。

最大の課題の一つが境界層遷移だ。層流から乱流への遷移位置によって壁面加熱率が3-8倍変わるから、TPS重量に直結する。Mack第2モード不安定性が極超音速では支配的で、これは音響波が境界層内に捕捉される不安定性だ。

Mach 4以上の境界層で現れる高周波の不安定モードだ。壁面近傍の一般化変曲点に起因する第1モードとは異なり、第2モードは境界層の相対的超音速領域で音響波が反射・増幅されることで成長する。周波数は典型的に数百kHzで、境界層厚の約2倍の波長を持つ。

RANSでは遷移を直接予測できないから、eNメソッドやDNS/LESを使う。NASAのeNコード（LASTRAC）やPSE（Parabolized Stability Equations）で線形安定解析を行い、N因子が臨界値（通常N=5-10）に達する位置を遷移点とする。最近はDNSで非線形段階までフルに計算する研究も増えている。

そう。再突入体ではPICA（Phenolic Impregnated Carbon Ablator）やSiCなどのアブレーション材が使われる。材料の熱分解ガスが表面から噴出し（blowing effect）、境界層を押しのけて加熱率を低減するんだ。この連成解析にはCFDと材料応答コード（FIAT, TITAN等）のカップリングが必要になる。

その通り。ALE（Arbitrary Lagrangian-Eulerian）法やoverset grid法で表面後退に対応する。NASAのCHAR（Charring Ablation Response）コードがDPLRとカップリングされた実績がある。

活発に研究されている。主な応用分野は3つだ。

ただし外挿性（訓練データの範囲外への汎化）が課題だ。極超音速は実験データが限られているから、Physics-Informed Neural Network（PINN）で支配方程式の制約を組み込むアプローチが注目されている。

極超音速特有のトラブルパターンを整理しよう。

Roe法で鈍頭物体の計算をすると、淀み点付近の弓状衝撃波が非物理的に凸凹になる現象だ。格子の対称線上で特に発生しやすい。

強い膨張波や衝撃波干渉の近傍で保存量から原始変数への変換時に非物理的な値が出ることがある。特に化学反応込みの場合、質量分率が0未満や1超になるケースも。

これは極超音速CFDで最もよくある問題だ。壁面加熱率は温度勾配 $\partial T / \partial n |_{wall}$ に比例するから、壁面近傍のメッシュ解像度に強く依存する。

化学反応ソース項のstiffnessが原因であることが多い。

衝撃波スタンドオフ距離は気体モデルに敏感だ。理想気体モデルだとスタンドオフが大きすぎ、化学平衡モデルだと小さすぎることがある。非平衡モデルが正解に近い場合が多い。

その通り。数値手法のデバッグだけでなく、「この条件で本当に非平衡効果は重要なのか」「壁面触媒条件は妥当か」といった物理的な判断が求められる。それが極超音速CFDの難しさでもあり面白さでもある。