Supersonic jet impingement on concave surfaces

技術要約：凹面への超音速噴流の衝突

問題提起
非理想的な膨張状態にある超音速噴流が表面に衝突する際のエアロアコースティック共鳴（スクリーチ）は、よく知られた現象であるが、壁面の曲率がこの共鳴に与える影響については未だ十分に理解されていない。自由噴流および平面壁への衝突におけるフィードバックループ機構（下流へ対流する不安定波と、上流へ伝播する音響波または誘導波からなるもの）は確立されているが、凹面がフィードバック経路、トーン振幅、およびモード選択をどのように変化させるかという具体的な役割は不明である。凹面壁は音響ミラーとして機能し、上流波を収束させる可能性がある一方、有効なループ長やソース強度を変化させる可能性もある。本研究では、ガウス型凹状の窪みの広がり（spread）が、流体の非定常性、表面荷重、およびスクリーチ特性に与える影響を調査し、グローバルなループ長と上流波の分散の役割を区別する。

手法
本研究では、理想的な膨張数 $M_j = 1.56$、レイノルズ数 $6 \times 10^4$ において、圧縮性ラージエディシミュレーション（LES）を用いて6つの衝突構成を調査する。幾何学的形状には、ノズル・壁間距離 $L/D = 2.08$ および $2.58$の2種類の平面板と、固定された深さを持ちつつ広がりパラメータ（$\sigma \in {0.4, 0.8, 1.6, 4.0}$）を変化させた4種類のガウス型凹面が含まれる。

数値計算には、構造化マルチブロックメッシュ（約1億2700万セル）上で3次元非定常圧縮性ナビエ・ストークス方程式を解くための、高次かつ低分散な差分法ソルバーを用いる。これらのシミュレーションは、せん断層の発達、衝撃波セル構造、およびマッハ・ディスクのダイナミクスを含む完全なフィードバックループを捉える。

主な解析手法は以下の通りである：

• 固有直交分解 (POD): 主要なコヒーレント構造および振動モードを特定するため。
• Powell-Tam ソース・トランスファー・バジェット: 曲率による増幅を、マッハ・ディスクのソース強度（$\eta_g$）と上流への伝達／受容効率（$\eta_u \eta_r$）の寄与に分割するため。
• 渦シートモデリング: 誘導ジェットモード（GJM）の分散関係を予測し、LESデータと比較するため。
• 波数・周波数解析: 内部（噴流内に閉じ込められた）成分と外部（周囲環境）の上流伝播波を区別し、それらの分散特性を決定するため。

主要な結果

1. 曲率依存的なトーン増幅:
   窪みの広がり（$\sigma$）が狭くなる（すなわち、凹面壁がジェット軸の周囲により局在化する）につれて、一次スクリーチ・トーンの振幅は著しく増加する。最も狭い窪み（$\sigma=0.4$）は、$L/D=2.6$ の平面壁参照値に対して最大23 dBの一次トーン増幅を示す。この増幅は、以下の2つの相補的なメカニズムに起因する：

   • ソース増幅: 閉じ込められた再循環領域における衝撃波と渦の相互作用の強化による、マッハ・ディスクの振幅増大。
   • 上流伝達効率: 上流波のノズルへの帰還の強化。Powell-Tam バジェットによれば、最も狭い窪みにおいて、全増幅に対する寄与はソースが約62%、上流伝達が約38%であった。
   • 収束メカニズム: 上流へ伝播する波は、外部の周囲波ではなく、ジェット噴流内に閉じ込められた成分によって支配されていることが判明した。凹面がジェットの柱の下に限定されている場合、周囲のチャネルと比較して、カラム・チャネル内の振幅は最大15.8 dBまで大幅に増大する。これは、凹面壁が主に噴流を通じて戻ってくる波に作用していることを示唆している。

2. モード選択と対称性:
   窪みの広がりに基づく振動モードの明確な分岐が特定された：

   • 軸対称モード (Mode A): 最も狭い窪み（$\sigma=0.4$）および $L/D=2.08$ の平面壁で見られる。これらは、マッハ・ディスクの強い周期的な軸方向の脈動と、トロイダルなジェット柱の振動を示す。
   • ヘリカルモード (Mode C): より広い窪み（$\sigma \ge 0.8$）および $L/D=2.58$ の平面壁で見られる。これらは、弱く広帯域なマッハ・ディスクの運動と、回転するヘリカルなジェット柱の振動を示す。
   • 軸対称モードからヘリカルモードへの遷移は、窪みがジェット柱の半径を超えて広がることと一致している。

3. 表面荷重特性:
   スクリーチ・モードの対称性は、非定常な表面荷重の性質を決定する：

   • 軸対称スクリーチ: トナールな法線（軸方向）方向の力を生成するが、横方向の力や曲げモーメントは無視できる。
   • ヘリカルスクリーチ: 広帯域な法線方向の力を生成するが、トナールなエネルギーを回転する面内曲げモーメントへと集中させる。このモーメントはジェット軸の周りを楕円軌道を描いて歳差運動し、軸方向荷重とは異なる特定の疲労リスクをもたらす。最大曲げモーメントは、必ずしも最大の曲率においてではなく、中間的な窪みの広がりにおいて発生する。

4. 周波数選択メカニズム:
   研究により、スクリーチ周波数を支配する2つの異なるメカニズムが明らかになった：

   • ヘリカルモード: ヘリカルなケース（凹面および平面の両方）の周波数は、渦シートモデルによって予測される $H_1$ 誘導ジェットモード（GJM）枝の低周波限界と密接に一致する。これは、等価な理想膨張噴流のせん断層プロファイルによって決定され、壁の幾何学的形状やループ長にはほとんど影響されない。
   • 軸対称モード: 軸対称のケースの周波数は、Powell の古典的なループ長基準（$N/f = \ell/u_c + \ell/a_\infty$）に従い、渦シートモデルのニュートラルな $A_n$ 枝とは一致しない。これは、これらの構成における軸対称スクリーチについては、ループ長が主要な選択基準であり、誘導モードの制約は適用されないか、あるいは満たされていないことを示唆している。

意義と主張
本論文は、壁の曲率がスクリーチ振幅および表面荷重の有効な制御メカニズムを提供することを実証しており、最も顕著な増幅は、凹状の窪みが半径方向にジェット柱に閉じ込められている場合に発生すると主張している。

極めて重要な点として、本研究は、衝突噴流におけるヘリカル・スクリーチと軸対称スクリーチの背後にある物理的メカニズムを区別している。ヘリカル・トーンは誘導ジェットモード（特に $H_1$ 枝）の分散特性によって支配されており、壁の形状の変化に対して頑健である一方、軸対称トーンはグローバルなフィードバックループ長によって支配されている。また、本研究の結果は、たとえ凹面壁が理論的に音響ミラーとして機能する可能性があったとしても、上流へのフィードバック経路は主に（周囲ではなく）ジェット噴流内を伝わっていることを明らかにしている。

これらの知見は、衝突噴流におけるフィードバックループの理解を精緻化し、ソース強度（マッハ・ディスクのダイナミクス）と上流伝達効率の効果を分離するとともに、選択された振動モードに基づいて、疲労に関連する荷重パターン（軸方向荷重か回転曲げ荷重か）を予測するための基礎を提供するものである。