衝撃波-境界層干渉 — トラブルシューティング
より充実した内容を shock-boundary-layer.html でご覧いただけます。
収束しない・発散する
SBLI解析で計算が発散してしまうとき、まず何を疑えばいいですか？
圧縮性流れの発散原因のトップ3を教えよう。
1. CFL数が高すぎる: 密度ベースソルバーの陽的時間進行ではCFL < 1.0が安定条件。陰的ソルバーでもSBLIのような強い不連続では初期CFL=0.5から徐々に上げるのが安全だ。
2. メッシュ品質不良: 衝撃波が通過するセルのアスペクト比が極端に大きいと数値振動が増幅する。特にプリズム層の頂部付近で衝撃波が境界層外縁と交差する場合、プリズム層からテトラ層への遷移部の品質が重要。
3. 初期条件の不適切さ: 一様流から始めると強い衝撃波が初期の非物理的な過渡状態を引き起こす。一段低いマッハ数の解を初期条件にするか、マッハ数をランプ関数で徐々に上げる方法が有効だよ。
衝撃波がぼやける
衝撃波が何セルにも広がってしまって、シャープに捕らえられないんですが...
考えられる原因と対策を整理しよう。
| 症状 | 原因 | 対策 |
|---|---|---|
| 衝撃波が5セル以上に広がる | 1次精度スキーム使用 | 2次精度風上スキーム（Roe, AUSM+）に変更 |
| 衝撃波位置がメッシュ依存 | メッシュが粗すぎる | 衝撃波近傍のメッシュ密度を2倍に |
| 衝撃波前後に振動（Gibbs現象） | リミッターなし高次スキーム | TVDリミッター（van Leer, Superbee）を適用 |
| 斜め衝撃波がstaircase状 | 格子と衝撃波が非整列 | 格子線を衝撃波方向に整列させる |
リミッターの選び方にコツはありますか？
van Leerリミッターは安全だけどやや散逸が大きい。Superbeeリミッターはシャープだけど振動しやすい。Minmodは最も散逸が大きい。SBLIでは圧力フィールドにvan Leerを使い、もし衝撃波がまだぼやけるならSuperbeeに切り替えるというアプローチが定番だよ。
剥離バブルが再現されない
実験では明確な剥離があるのに、CFDでは剥離が出ないことがあるんですが...
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の使用を検討しよう
壁面の$y^+$を1以下にしても改善しない場合は、もうLES/DESに行くしかないですか？
RANSの範囲でもう一つ試せるのが、QCR（Quadratic Constitutive Relation）の有効化だ。標準のBoussinesq近似では表現できない二次の応力項を追加することで、コーナー流れや強い逆圧力勾配での予測精度が改善される。FluentではSST $k$-$\omega$ にQCRオプションを付けられるよ。
衝撃波-境界層干渉で「剥離バブル」が消えたり現れたりする謎
衝撃波-境界層干渉のCFD解析で初心者が戸惑う現象の一つが、メッシュを細かくするたびに剥離領域の大きさが変わることだ。粗いメッシュでは「剥離なし」と判定されたのに、精密メッシュにすると小さな剥離バブルが出現する。これはバグではなく、剥離のスケールが非常に小さく壁面第1層のY+が大きいと掴めないことが原因だ。衝撃波直後の逆圧力勾配は極めて急峻で、Y+≦1程度のメッシュが必要になる。実務では「まずY+マップを確認して境界層内に十分な解像度があるか確かめる」がトラブルシューティングの第一歩だ。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——衝撃波-境界層干渉の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告