圧縮性乱流モデリング — トラブルシューティングガイド

実務でよく遭遇する問題と対策を整理しよう。

症状: 超音速混合層の拡散角が実験値の2〜3倍になる

原因: 圧縮性補正が無効。標準k-epsilonやk-omegaは非圧縮で校正されているため、圧縮性効果（成長率抑制）を予測できない。

対策:

• Sarkar or Zeman圧縮性補正を有効化
• Fluentの場合: Viscous Model → Options → Compressibility Effects にチェック
• 対流マッハ数 $M_c = (U_1 - U_2)/(a_1 + a_2)$ を確認。$M_c > 0.3$ なら補正必須

症状: 垂直衝撃波の直後で乱流運動エネルギー $k$ が数倍にスパイクする

原因: 衝撃波による圧力勾配を乱流生成項 $P_k$ が過大評価。Production limiterが未設定。

対策:

• SST k-omega の Production Limiter を有効化: $P_k = \min(P_k, C_{lim} \cdot \bar{\rho} \beta^* k \omega)$、$C_{lim} = 10$
• Fluentでは「Production Limiter」または「Kato-Launder modification」を使用
• メッシュを衝撃波近傍で細かくして、衝撃波の数値的な厚みを減らす

症状: 計算開始から数十ステップで残差が発散、NaN発生

原因と対策:

• 初期場が不適切: 全領域をM=0の静止場で初期化すると、超音速入口との不整合で衝撃波が初期に発生し発散する。→ 等エントロピー関係で適切な初期場を設定
• CFL数が大きすぎる: 初期は CFL = 0.5〜1.0 から始めて徐々に上げる
• 境界条件の不整合: 超音速出口にpressure outletを使っている場合、背圧が流れ場と矛盾する可能性。→ 超音速出口ではextrapolation（外挿）境界条件を使用

症状: 境界層が物理的根拠なく剥離する。特に付着流の領域で発生。

原因: DESのRANS→LES切り替えが境界層内で発生し、RANS的に解くべき領域がLES格子解像度不足で解かれてしまう。

対策:

• DDES（Delayed DES）を使用。$f_d$ シールド関数で境界層をRANS領域として保護
• IDDES（Improved DDES）に切り替え。WMLES機能を含み、より安定
• 格子設計を見直し: 壁面平行方向の格子アスペクト比が大きすぎないか確認

多くの場合はそうだけど、メッシュの遷移領域（RANS→LES）の設計は依然として重要だ。Log-layer mismatch（対数層の不一致）の問題は完全には解決されていないから、壁面摩擦係数や速度プロファイルを必ず検証すること。

そうだね。最低限、乱流マッハ数の確認、圧縮性補正の有効化、Production limiterの設定、境界条件の整合性チェック、この4点は毎回確認する習慣をつけておくといいよ。