超音速流れ — 先端トピックとスクラムジェットCFD

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

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

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

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

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

1. Mach波放射: ジェットせん断層の大規模渦が超音速で移動するときに発生

2. Screech tone: 衝撃波セル構造とフィードバックループによる離散周波数の騒音

3. Broadband shock-associated noise: 衝撃波セルと乱流の干渉

これらを正確に予測するには、衝撃波セル構造を2-3セルで解像できるメッシュ密度と、LESの空間フィルタ幅がジェットせん断層の渦構造を十分に解像するレベルが必要だ。

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

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

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