乱流モデルの曲率・回転補正 — トラブルシューティングガイド
この記事は統合版に移行しました
より充実した内容を curvature-correction.html でご覧いただけます。
より充実した内容を curvature-correction.html でご覧いただけます。
乱流モデルの曲率・回転補正 — トラブルシューティングガイド
トラブルシューティング
曲率補正を入れたら計算が不安定になったんですが、何が原因ですか？
よくあるケースを整理しよう。
よくある問題と対策
1. 収束性の悪化
残差が振動して収束しません。
原因: $f_{rot}$ が反復ごとに大きく変動している可能性がある。特に初期の流れ場が未発達な段階では $r^*$ や $\tilde{r}$ の計算が不安定になりやすい。
対策:
- 補正なしでまず収束させ、その後補正を有効にして再計算
- Fluentでは Under-Relaxation Factor を乱流量について 0.5-0.7 に下げる
- $f_{rot}$ の変化を緩和するために内部で時間平均を取る（一部ソルバーは自動で実施）
- 擬似時間ステップを小さくする
2. 旋回流の過小予測が改善しない
補正を入れたのに実験と合いません。
原因と対策:
- メッシュ不足: 周方向・半径方向の解像度を2倍にして再計算。数値散逸による旋回減衰はモデル補正では直せない
- 離散化スキーム: 一次風上差分は旋回を消す。Bounded Central Difference やQUICKに切り替え
- 乱流強度の初期条件: 入口の乱流強度が高すぎると旋回コアが早期に拡散する。実験値に合わせる
- 補正の限界: 強旋回流（$Sw > 2$）ではRANS+補正では限界がある。RSMかDDES/LESを検討
3. 非物理的な乱流粘性
渦粘性比（$\mu_t/\mu$）が局所的に10万を超えるセルがあります。
原因: $f_{rot} > 1$ の領域で乱流生成が過剰に増強されている。
対策:
- Production Limiter（$P_k \leq C_{lim} \cdot \rho \varepsilon$ や $P_k \leq 10 \beta^* \rho k \omega$）を有効にする
- $f_{rot}$ の上限を1.25から1.0に変更して試す
- 該当セルのメッシュ品質を確認（高スキューネスが原因の場合あり）
4. 定常計算での振動
定常解析なのにモニター値が周期的に振動します。
原因: 旋回流やカルマン渦列のような本質的に非定常な流れを定常解析で計算している。曲率補正により物理的な不安定性がより顕在化した。
対策:
- 非定常解析（URANS）に切り替え
- 時間平均量で比較
- 定常解析のままにする場合は残差の振動幅が十分小さいことを確認し、平均値を使用
問題の切り分けが大事なんですね。補正の問題なのかメッシュの問題なのか、まず補正なしで確認するのがポイントだ。
その通り。補正なしでベースラインを確立し、補正の効果を差分で評価する。二つの変数を同時に変えないのがデバッグの鉄則だ。
Coffee Break よもやま話
「曲率補正モデルで計算が発散した」——数値的不安定性の診断
曲率補正乱流モデル（SSTcc等）は標準SSTよりも数値的に不安定になることがある。補正係数r*が急変する領域（剥離・再付着点付近）で生成率が急増・急減し、k方程式が非物理的な解に飛ぶことがある。診断手順：①補正係数r*の空間分布を可視化して急変領域を特定、②補正係数の上限クリップ（多くの実装でr*∈[0, 1.25]に制限）を確認、③緩和係数を下げてより保守的な収束を促す。OpenFOAMではSSTccのソースコードを確認してリミッターの設定を調整することで安定性を改善できる。標準SSTで先に収束させてからSSTccに切り替えるアプローチも有効だ。
関連トピック
この記事の評価
ご回答ありがとうございます！
参考に
なった
なった
もっと
詳しく
詳しく
誤りを
報告
報告