衝撃波-境界層干渉

圧縮性流れの中で衝撃波が壁面の境界層と交差するときに起こる強い相互作用だよ。衝撃波による急激な圧力上昇が境界層に伝わると、境界層内の低速流体は逆圧力勾配に耐えられずに剥離してしまう。この剥離バブルと衝撃波が互いに影響し合って、非常に複雑な流れ場を作るんだ。

代表的なのはインテーク内部の衝撃波/境界層干渉、翼面上の遷音速衝撃波による剥離、そしてスクラムジェットエンジン内の衝撃波反射だね。いずれも空力性能の急激な劣化、非定常荷重、あるいは機体振動に直結するから、航空宇宙分野では設計上の最重要課題の一つなんだ。

基本は圧縮性Navier-Stokes方程式だけど、SBLIで特に重要なのは入射衝撃波前後の関係式と、境界層の剥離判定だよ。まず斜め衝撃波の圧力比はこう表せる。

その通り。自由干渉理論（Free Interaction Theory）によれば、剥離が始まる圧力上昇は上流の境界層状態で決まる。層流境界層の場合、剥離に必要な圧力上昇は非常に小さいけど、乱流境界層は壁面近くの運動量が大きいから剥離しにくい。干渉の強さを示す指標として、干渉圧力パラメータがよく使われるよ。

まさにそこがSBLIの難しさだ。衝撃波が境界層の遷移を誘起する場合もあるし、衝撃波の上流で既に乱流に遷移している場合もある。衝撃波誘起遷移（Shock-Induced Transition）はスクラムジェットのインテーク設計で特に重要で、遷移位置の予測精度がエンジン性能の推算に直結する。

Edneyの分類が有名で、衝撃波と境界層の交差形態によってType I からType VIまで6種類に分けられる。代表的なものを紹介しよう。

非干渉時の数倍から十数倍の局所熱流束が発生することがある。スペースシャトルの前縁部で問題になった事例が有名だよ。設計段階でこれを見逃すと熱防御系が耐えられないから、CFDでの事前予測が不可欠なんだ。

これがSBLIのCFDにおける最大の課題と言ってもいい。RANS乱流モデルは衝撃波による急激な逆圧力勾配のもとでの境界層挙動を正確に予測するのが難しいんだ。特に剥離バブルのサイズと再付着位置の精度が問題になる。

大きく変わるよ。代表的なモデルの傾向をまとめると以下の通りだ。

Wall-Resolved LESが最も高精度だけど、壁面近傍のメッシュ要件が $\Delta x^+ \sim 50$, $\Delta y^+ < 1$, $\Delta z^+ \sim 20$ と非常に厳しい。SBLI問題のように高Re数の壁面流れではセル数が膨大になる。そこで実務的にはDES（Detached Eddy Simulation）やWMLES（Wall-Modeled LES）が有効だよ。

DESの考え方は、壁面近傍はRANS、剥離領域やせん断層はLESで解くハイブリッド手法だ。Ansys FluentではDES、DDES（Delayed DES）、SBES（Stress-Blended Eddy Simulation）が選択できる。STAR-CCM+ではIDDES（Improved DDES）が標準的に使われているよ。

DDESは遮蔽関数でRANS領域とLES領域を分けるけど、切り替え領域で「Grey Area」と呼ばれる曖昧な領域が生じやすい。SBESはこのGrey Area問題を改善した手法で、RANS応力テンソルからLES応力テンソルへの遷移がよりスムーズになる。Ansys Fluent 2020R1以降で利用可能だよ。

衝撃波の捕捉には風上系スキーム（upwind scheme）が基本だ。代表的な選択肢を示すよ。

LES/DESでは数値拡散が乱流構造を過度に減衰させるため、衝撃波がない領域では低散逸の中心差分系スキームを使い、衝撃波付近のみ風上系に切り替えるハイブリッドスキームが理想的だ。OpenFOAMのlocalBlendedやFlunetのHybrid Bounded Central Differencingがこのアプローチに該当する。

SBLIではメッシュの質が結果に極端に効くから、以下の3領域を重点的に細分化する必要がある。

1. 衝撃波足元近傍: 衝撃波が壁面と交差する領域。ここは剥離バブルの起点だから、流れ方向に衝撃波厚みの数倍以内のセルサイズが必要。具体的には衝撃波角度の予測位置前後に $\Delta x / \delta \lesssim 0.05$（$\delta$ は境界層厚さ）のメッシュを配置する。

衝撃波の方向が予測可能なら構造格子が圧倒的に有利だ。衝撃波に対して格子線を平行・垂直に揃えられるため、数値拡散を最小化できる。非構造格子の場合はAMR（Adaptive Mesh Refinement）を使って衝撃波近傍を自動的に細分化するアプローチが実用的だよ。Fluent、STAR-CCM+ともにAMRをサポートしている。

超音速流れでは上流から全ての情報が流入するから、入口には全温度 $T_0$、全圧 $p_0$、流れ方向を指定する。もしくは超音速入口ではリーマン不変量に基づいて全変数を固定する。出口は超音速なら外挿、遷音速なら静圧指定だ。

壁面は断熱壁（$\partial T/\partial n = 0$）か等温壁（$T_w = \text{const}$）を選ぶ。衝撃波/境界層干渉では壁面温度条件が剥離バブルサイズに影響するので注意が必要だよ。冷却壁（$T_w / T_{aw} < 1$）は境界層内の密度を上げて壁面近くの運動量を増すため、剥離が抑制される傾向がある。

以下がSBLI研究で標準的に使われるベンチマーク問題だよ。

その通り。Settlesの24度圧縮ランプ問題は格子収束性の確認に最適だ。3水準以上のメッシュ密度で壁面圧力分布とプラトー圧力を比較し、Richardson外挿で収束次数を確認するのが定石だよ。

FluentでSBLIを解くときの推奨設定はこうだ。

密度ベースソルバーは保存変数を直接解くため、衝撃波前後の保存則が厳密に満たされる。圧力ベースソルバーでもCoupled Schemeで圧縮性流れは解けるけど、強い衝撃波がある場合は密度ベースの方がロバストだ。ただしFluent 2023R1以降では圧力ベースのCoupled Schemeもかなり改善されていて、遷音速域では十分使えるよ。

OpenFOAMではrhoCentralFoamが圧縮性流れの標準ソルバーだ。中心差分ベースのKurganov-Tadmorスキームを使っていて、衝撃波捕捉能力と低散逸性のバランスが良い。

設定のポイントはこうだよ。

sonicFoamはPISO法ベースの圧力-速度連成で、rhoCentralFoamはRiemann解法ベースだ。衝撃波の解像度はrhoCentralFoamの方が高いが、亜音速域の精度はsonicFoamの方が良い場合もある。遷音速域ではどちらも使えるけど、強い衝撃波がある場合はrhoCentralFoamを推奨する。

STAR-CCM+では以下の設定が推奨される。

本当だよ。SBLIの非定常性は二つの周波数帯に分けられる。高周波成分は上流境界層内の乱流構造に起因するもので、$St = fL/U_{\infty} \sim O(1)$ 程度。低周波成分は衝撃波全体が前後に動く現象で、$St \sim O(0.01-0.05)$ と非常にゆっくりした振動だ。

これは長年の議論の的で、二つの仮説がある。一つ目は上流メカニズム説で、上流境界層内の大規模構造（superstructures）が衝撃波位置を変動させるという考え。二つ目は下流メカニズム説で、剥離バブル内の循環流が不安定性を持ち、それが衝撃波にフィードバックするという考えだ。

最近のDNS研究（Priebe & Martin, 2012; Pasquariello et al., 2017）では、両方のメカニズムが関与しているという見方が主流になってきている。低周波振動の振幅は壁面圧力変動の20-30%に達することもあるから、構造の疲労設計には大きなインパクトがある。

現時点ではまだ低Reynolds数（$Re_{\delta} \sim O(10^4)$）に限られるけど、メカニズムの理解には極めて有用だ。代表的な成果を紹介しよう。

2D翼型の遷音速SBLIなら、スパン方向周期境界で50-100Mセル程度。3D形状の場合は数億セルになるから、HPCクラスタが必須だ。計算時間は数万〜数十万コア時間のオーダーだよ。

活発に研究されているよ。主な応用は以下の三つだ。

Fluent 2024R1以降ではUDF（User Defined Function）を通じて機械学習モデルをソルバーに組み込むフレームワークが整備されてきている。ただし、学習データの範囲外への外挿には慎重でなければいけないよ。

圧縮性流れの発散原因のトップ3を教えよう。

1. CFL数が高すぎる: 密度ベースソルバーの陽的時間進行ではCFL < 1.0が安定条件。陰的ソルバーでもSBLIのような強い不連続では初期CFL=0.5から徐々に上げるのが安全だ。

考えられる原因と対策を整理しよう。

van Leerリミッターは安全だけどやや散逸が大きい。Superb​eeリミッターはシャープだけど振動しやすい。Minmodは最も散逸が大きい。SBLIでは圧力フィールドにvan Leerを使い、もし衝撃波がまだぼやけるならSuperbeeに切り替えるというアプローチが定番だよ。

RANS解析でよくある問題だね。チェックポイントは以下の通りだ。

1. 壁面メッシュは十分か: $y^+ > 5$ だと壁関数が働いて壁面せん断応力を過大評価し、剥離が抑制される。$y^+ < 1$ を厳守すること

RANSの範囲でもう一つ試せるのが、QCR（Quadratic Constitutive Relation）の有効化だ。標準のBoussinesq近似では表現できない二次の応力項を追加することで、コーナー流れや強い逆圧力勾配での予測精度が改善される。FluentではSST $k$-$\omega$ にQCRオプションを付けられるよ。