標準k-ωモデル（Wilcox） — トラブルシューティングガイド
より充実した内容を k-omega-standard.html でご覧いただけます。
トラブルシューティング
標準k-ωで特有の問題はありますか？
SSTと共通するものも多いが、標準k-ω特有の問題がある。
1. 自由流感度（1988版）
入口のω値を変えたら結果が大きく変わりました。
原因: 1988版のk-ωは自由流の $\omega$ 境界値に敏感。$\omega_{inlet}$ を10倍に変えると壁面Cf が20%以上変わることがある。
対策:
- 2006版に切り替える（クロス拡散項で改善）
- SSTに切り替える（根本的解決）
- どうしても1988版を使う場合、$\omega$ の感度解析を実施し結果の不確かさを評価
2. 遷移のない全乱流仮定
低Reynolds数の流れで層流領域が計算されません。
原因: k-ωは全乱流を仮定している。Low-Re修正を使っても自然遷移の予測は不十分。
対策:
- γ-Reθ遷移モデル（Transition SST）を使用
- 層流領域を手動で指定（乱流粘性をゼロに固定）
- 層流-乱流の領域を分けてメッシュを作成
3. 後流域での過剰拡散
円柱後流のカルマン渦が早く消えてしまいます。
原因: 定常RANSの限界。渦粘性が後流域で過大になり渦を拡散させる。k-ωもSSTも同じ問題を持つ。
対策:
4. ωの初期値設定
ωの初期値をどう設定すればいいかわかりません。
以下の手順で推定する。
1. 乱流強度 $I$ と長さスケール $l_t$ を仮定
2. $k = \frac{3}{2}(U I)^2$
3. $\omega = \frac{k^{0.5}}{C_\mu^{0.25} l_t}$ または $\omega = \frac{\rho k}{\mu \cdot (\mu_t/\mu)}$
粘性比 $\mu_t/\mu$ は外部流れで1-10、内部流れで10-100が目安。ωが小さすぎる（$\mu_t$ が大きすぎる）と初期の発散、大きすぎる（$\mu_t$ が小さすぎる）と乱流効果が弱くなる。
標準k-ωのトラブルは自由流感度に起因するものが多いんですね。
そう。SSTが広く使われる最大の理由がこの自由流感度の解消だ。
壁上でのω→∞問題——初期値設定の悩み
k-ωモデルの実装で初学者が困るのが「壁面でのωが理論上無限大になる」という問題です。実際の数値計算では無限大は扱えないため、壁から最初のセルでのωに対して解析解に基づく近似値を設定します。この値が大きすぎると数値的に不安定になり、小さすぎると境界層が正しく解けません。OpenFOAMのomegaWallFunctionが実は非常に繊細な実装で、バージョンによって振る舞いが変わった歴史があります。「壁面ωの初期化に丸一日悩んだ」という経験談は、k-ωユーザーの間でよく聞く話です。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——標準k-ωモデル（Wilcox）の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告