混合層 — トラブルシューティングガイド

原因と対策:

• 数値拡散が大きい: 1次風上スキームや粗いメッシュ。中心差分（2次以上）に変更し、メッシュを細かくする
• 初期擾乱がない: tanh プロファイルだけでは擾乱が成長するまで非常に長い時間がかかる。最不安定モードの擾乱を振幅 $0.01 \Delta U$ 程度で加える
• 計算領域が狭い: 最不安定モードの波長が計算領域に収まっていない

確認ポイント:

• 計算領域長さ: 基本モードの波長の4倍以上が必要。ペアリングにはサブハーモニック波長の空間が要る
• サブハーモニック擾乱: 初期条件に基本モードの半分の波数の擾乱を加えると、ペアリングが促進される
• 計算時間: ペアリングには時間がかかる。渦度場のスナップショットを時系列で確認

対策:

• 標準 $k$-$\varepsilon$ の場合、$C_{\varepsilon 2}$ を $1.83$ から $1.92$ に変更するフィッティングが報告されている（ただし汎用性は失われる）
• SST $k$-$\omega$ のほうが平面混合層の予測精度が良い場合がある
• Reynolds Stress Model (RSM) は異方性を捉えるので、混合層のReynolds応力の予測精度が向上する

原因: スパン方向の初期擾乱が不十分か、スパン方向の計算領域が狭い。

対策:

• スパン方向に3Dモードの擾乱（oblique mode）を加える
• $L_z$ を最不安定3D波長の2倍以上に設定
• スパン方向の解像度を確認（$\Delta z < \lambda_{3D,min} / 10$）

時間的混合層では水平面（$x$-$z$平面）での空間平均が取れる。ただし以下に注意が必要だ。

• ペアリング後の混合層は非均一: ペアリングにより混合層厚さが場所によって異なる。特定の瞬間ではなく、多数の時間ステップの平均も組み合わせる
• 混合層中心の追跡: 混合層は $y$ 方向にドリフトすることがある（速度比 $r \neq 0$ の場合）。$\bar{u} = U_c$ の $y$ 座標を基準にプロファイルを揃える

K-H不安定性の理論、ペアリングのダイナミクス、圧縮性効果、密度比の影響と、流体力学の重要な概念が凝縮された問題だ。CFD初学者の次のステップとして最適な題材だと思う。