ノズル流れ — 先端技術と研究動向

大気圧が変化する上昇中にノズル膨張比を最適化する「高度補償ノズル」が注目されている。代表的な概念は3つだ。

• デュアルベルノズル: 内側ベルと外側ベルの2段構成。低高度では内側ベルの端で流れが剥離し、高高度では外側ベルまで完全に膨張する
• プラグノズル（エアロスパイクノズル）: 中心プラグの周りに膨張。外側は大気圧で自動的に調整される
• 伸展-偏向ノズル（E-Dノズル）: 壁面の一部を省略して高度補償効果を得る

Raptorは高膨張比のベル型ノズル（真空版）を使っている。デュアルベルは概念実証段階だが、DLR（ドイツ航空宇宙センター）やJAXAで実験的研究が進んでいる。CFDではノズル壁面の屈折点での流れ剥離-再付着の遷移を非定常で正確に予測する必要があり、DDES（Delayed DES）やWALE SGSモデルのLESが使われる。

そう。固体推進薬の燃焼ガスにはAl₂O₃粒子（直径1-10μm）が含まれ、比推力の低下やノズル侵食の原因になる。この二相流ノズルの計算にはEulerian-Lagrangianアプローチが使われる。

• 連続相（ガス）: Euler方程式
• 分散相（粒子）: Lagrangian追跡。粒子-ガス間の抗力・熱伝達のカップリング

FluentのDPM（Discrete Phase Model）やSTAR-CCM+のLagrangian Multiphaseで対応可能だ。粒子がスロートを通過するときの速度遅れ（velocity lag）が推力損失に寄与する。

推力偏向の方式は複数あり、CFDの扱いも異なる。

二次流体噴射は最も活発に研究されている。主流の超音速膨張流に横方向から亜音速or超音速ジェットを噴射し、衝撃波を誘起して横力を発生させる。この衝撃波-境界層干渉を正確に予測するにはDES以上の乱流モデルが望ましい。

アジョイント法による感度解析が強力だ。推力やIsp（比推力）を目的関数として、ノズル壁面形状パラメータ（ベジエ曲線の制御点等）の勾配を効率的に計算できる。

• Ansys Fluent: adjoint solver搭載。壁面形状のメッシュモーフィングと連携
• SU2: オープンソースCFD。アジョイント設計最適化に強い。Stanford大学で開発
• DAFoam: OpenFOAMベースの離散アジョイント最適化フレームワーク

そう。SU2はEuler/RANSのアジョイント法が洗練されていて、翼型やノズルの形状最適化で多数の論文がある。無料で使えるから大学の研究にも最適だ。