衝撃波-境界層干渉 — トラブルシューティング

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

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

2. メッシュ品質不良: 衝撃波が通過するセルのアスペクト比が極端に大きいと数値振動が増幅する。特にプリズム層の頂部付近で衝撃波が境界層外縁と交差する場合、プリズム層からテトラ層への遷移部の品質が重要。

3. 初期条件の不適切さ: 一様流から始めると強い衝撃波が初期の非物理的な過渡状態を引き起こす。一段低いマッハ数の解を初期条件にするか、マッハ数をランプ関数で徐々に上げる方法が有効だよ。

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

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

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

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

2. 乱流モデルの選択: SA（Spalart-Allmaras）モデルは剥離を過小評価する傾向が強い。SST $k$-$\omega$ に切り替えるだけで改善することが多い

3. 2D解析の限界: SBLIの剥離は本質的に3次元的だから、2D RANSでは剥離長さが実験と合わないことがある。3Dメッシュでスパン方向周期条件を使うと改善する場合がある

4. 遷移の考慮: 層流境界層を仮定して解くべきところを完全乱流で解くと剥離が小さくなる。Fluent の$\gamma$-$Re_{\theta}$ 遷移モデルやSTAR-CCM+のGamma Transition Modelの使用を検討しよう

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