標準k-ωモデル（Wilcox）
標準k-ωモデル（Wilcox）の理論基礎
概要
先生！ 標準k-ωモデルとSST k-ωってどう違うんですか？
標準k-ωモデルはWilcox (1988, 2006)が開発したもので、SSTのベースとなったモデルだ。乱流の比散逸率（Specific Dissipation Rate）$\omega = \varepsilon/(C_\mu k)$ を輸送変数に使い、壁面近傍での処理がk-εモデルより自然だという利点がある。
壁面に強いってことですか？
そう。ωの壁面境界条件は解析的に定義でき（Dirichlet条件）、εのように壁面で特異性を持たない。壁関数なしで粘性底層まで直接解像できる。ただし自由流の境界条件に対する感度が高いという重大な弱点がある。
支配方程式
方程式を教えてください。
Wilcox (2006)の改訂版k-ω方程式:
k方程式:
ω方程式:
渦粘性: $\mu_t = \rho k / \omega$
モデル定数: $\alpha = 13/25$、$\beta = \beta_0 f_\beta$、$\beta^ = 9/100$、$\sigma = 1/2$、$\sigma^ = 3/5$、$\sigma_d = 1/8$（$\nabla k \cdot \nabla \omega > 0$ の場合のみ）。
$\sigma_d$ のクロス拡散項が条件付きなのが面白いですね。
Wilcox 2006版の重要な改良点だ。これにより自由流感度の問題がかなり改善された。1988年版ではこの項がなく、自由流の $\omega$ 値で解が大きく変わるという問題があった。
自由流感度問題
自由流感度って具体的にどういう問題ですか？
自由流（壁から離れた領域）での $\omega$ の境界値を変えると、壁面近傍の解まで変わってしまう。これは物理的に不合理だ。Menterがk-εとのブレンド（SST）を開発した主な動機がこの問題だった。
2006年版のクロス拡散項でかなり改善されたが、完全には解決していない。このため産業CFDではSST k-ωが標準k-ωより圧倒的に多く使われている。
Wilcoxが50年かけてアップデートし続けたモデル
標準k-ωモデルは1988年にDavid C. Wilcoxが提案して以来、1998年版、2006年版、2008年版と継続的に改訂されています。とくに自由流感受性（入口境界のω値に結果が敏感に依存する問題）が長年の弱点として知られており、Wilcox自身が「改善が必要だ」と認識して改訂を重ねてきました。研究者が自分のモデルの欠点を公言して改良し続けるというのは、工学の世界では珍しい誠実さです。使うときはどのバージョンを実装しているかを確認することが重要で、「標準k-ω」と一言で言っても年代によって式が異なります。
標準k-ωモデル（Wilcox）の数値計算手法
数値実装
標準k-ωの実装で特有の注意点はありますか？
ωの壁面境界条件が最大のポイントだ。
壁面境界条件
ωの壁面値はどう決めるんですか？
$y^+ \to 0$ の極限で:
これはDirichlet条件で、数値的にはいくつかの実装方法がある:
1. 壁面第1セル中心の値として設定: $\omega_P = \frac{6\nu}{\beta_1 (\Delta y_1)^2}$
2. 壁面Ghost Cellの値として設定: Wilcox推奨
3. 壁関数との切替え: $y^+ > 2.5$ では $\omega = u_\tau / (\sqrt{\beta^*} \kappa y)$
方法1だとメッシュが細かいほどωが大きくなりますよね。
そう。$y^+ = 1$ で $\omega \sim O(10^6)$、$y^+ = 0.1$ で $\omega \sim O(10^8)$ になる。これは行列の対角優位性を維持するのに注意が必要だ。
離散化スキーム
k-ωの離散化で推奨されるスキームは？
| 変数 | 推奨スキーム | 備考 |
|---|---|---|
| k | Second Order Upwind | TVD制限関数付き |
| ω | First/Second Order Upwind | ωは壁面で急激に変化するため一次精度でも許容 |
| 運動量 | Second Order Upwind以上 | |
| クロス拡散項 | Central Difference | 勾配の内積、条件分岐に注意 |
OpenFOAMでの設定
OpenFOAMでの使い方を教えてください。
OpenFOAMにはWilcox 2006版の kOmega が実装されている。
```
RAS
{
RASModel kOmega;
turbulence on;
printCoeffs on;
}
```
ただし実務では kOmegaSST を使うことが圧倒的に多い。kOmega を使うのは、特定のベンチマーク検証やSSTのブレンディングの影響を排除したい場合くらいだ。
収束性の改善
収束が悪いときのテクニックは？
k-ωとy+=1の相性——壁面解像の意外な恩恵
k-ωモデルは壁面近傍でω→∞という境界条件の性質を持ち、壁上でのωの解析解が存在するため、粘性底層を直接解像するy+=1程度の細かいメッシュと非常に相性が良いです。この特性のおかげで遷移流れや低Reynolds数域の境界層計算に強く、航空機翼の境界層解析に長く使われてきた経緯があります。逆に言えば「k-ωを選んだならy+は1以下を目指せ」という実務指針は、このモデルの物理的特性から自然に導かれます。メッシュが細かいほどk-ωの本領が発揮される、というのは実務エンジニアが知っておくべき重要な知識です。
標準k-ωモデル（Wilcox）の実務適用
実践ガイド
標準k-ωを使うべき場面はあるんですか？ SSTがあるなら不要では？
ほとんどの場合SSTで良いが、以下のケースでは標準k-ωが使われることがある。
使用が正当化されるケース
| ケース | 理由 |
|---|---|
| Wilcox論文との比較検証 | 同じモデルで再現する必要がある |
| 壁面主導の境界層流れ | 自由流の影響が小さい場合、SSTと同等の精度 |
| SSTのブレンディング影響の排除 | 研究用途でブレンディングの効果を切り分けたい |
| 旧来のバリデーション資産 | 過去の検証データがk-ωベースの場合 |
SST vs 標準k-ωの比較
具体的にどう違うんですか？
| 特性 | 標準k-ω | SST k-ω |
|---|---|---|
| 壁面近傍の精度 | 高 | 高（同等） |
| 自由流感度 | 中（2006版で改善） | 低（k-εブレンド） |
| 逆圧力勾配の剥離 | 過小予測傾向 | 改善（SSTリミッター） |
| 淀み点問題 | あり | Production Limiterで緩和 |
| DES/DDES化 | 可能だが稀 | 標準的 |
入口境界条件
入口のω値はどう設定しますか？
一般的な変換式:
または乱流粘性比から:
| パラメータ | 外部流れ | 内部流れ |
|---|---|---|
| 乱流強度 $I$ | 0.1-1% | 5-10% |
| 長さスケール $l_t$ | $0.01 L_{ref}$ | $0.07 D_h$ |
| 粘性比 $\mu_t/\mu$ | 1-10 | 10-100 |
自由流感度を考慮して、$\omega$ の入口値は感度解析で確認すべきだ。特に1988年版では結果が大きく変わりうる。
やはり実務ではSSTを使うのが安全そうですね。
そうだね。標準k-ωの存在意義はSSTの理解を深めるための基礎知識として重要ということだ。
宇宙機の熱防護解析でk-ωが使われる理由
再突入カプセルや超音速ミサイルの空力加熱解析では、k-ωモデルが標準的に使用されます。これは衝撃波後の高温境界層や壁面付近の熱流束予測でk-ωが比較的良好な精度を示すためです。JAXAやNASAの技術レポートでも超音速・極超音速の壁面熱流束評価にk-ωを用いた事例が多く見られます。「宇宙開発でも使われているモデルを普通の配管解析に使っている」というのは、k-ωの守備範囲の広さを示す面白いギャップです。
標準k-ωモデル（Wilcox）のソフトウェア比較
商用ツールでの対応
各ソルバーで標準k-ωはどう扱われていますか？
| ソルバー | 実装バージョン | 備考 |
|---|---|---|
| Ansys Fluent | Wilcox 1998版相当 | Low-Re補正、せん断流補正オプションあり |
| Ansys CFX | Wilcox版 | SSTを強く推奨。k-ωは後方互換性のため |
| STAR-CCM+ | Wilcox 2006版 | |
| OpenFOAM | Wilcox 2006版 | kOmega クラス |
Fluentでの注意点
Fluentの標準k-ωにはどんなオプションがありますか？
Fluent の k-omega (Standard) には:
- Low-Reynolds Corrections: 遷移領域での減衰関数
- Shear Flow Corrections: せん断流での拡散項修正
- Transitional Flows: Free shear layer での補正
これらはデフォルトでOFFなので必要に応じて有効にする。ただしFluentのマニュアルでもSSTの使用が推奨されている。
実装バージョンの違い
1988版と2006版で大きく違うんですか？
| 特性 | 1988版 | 2006版 |
|---|---|---|
| クロス拡散項 | なし | $\sigma_d/\omega (\nabla k \cdot \nabla \omega)$（条件付き） |
| 自由流感度 | 大きい | 改善 |
| $\beta$ | 定数 | $f_\beta$ で変動（渦伸張効果） |
| Realizability | なし | ストレスリミッター追加 |
2006版のクロス拡散項はSSTのクロス拡散項と形式は似ているが、ブレンディング関数を使わない点が異なる。
結局、標準k-ωを選ぶ合理的な理由は限られるということですね。
そうだ。学術的な比較研究、過去のバリデーションとの整合性確保、教育目的以外では、SSTを使うのが実務的に正しい選択だ。
「Wilcox 1988 vs 2006」——ソルバーのバージョン迷宮
同じ「標準k-ω」でも、ANSYS Fluent・STAR-CCM+・OpenFOAMでどのWilcox年代版を実装しているかが異なります。例えばOpenFOAMのkOmega辞書はWilcox 2006ベースですが、Fluent旧来バージョンは1988年式に近いものを使っていました。プロジェクト間でソルバーを変えると「同じ設定のはずなのに結果が違う」が起きやすく、チームで解析を引き継いだときにはソルバーのリリースノートを確認してモデル方程式の差異を把握することが欠かせません。
標準k-ωモデル（Wilcox）の先端研究
先端トピック
標準k-ωの研究的な発展はありますか？
Wilcox自身が2006年版以降も改良を続けている。
Wilcox 2006版の位置づけ
2006年版でどのくらい改善されたんですか？
Wilcoxは著書 "Turbulence Modeling for CFD" (2006, 3rd ed.) で大幅な改良を行った。主要な改善:
1. クロス拡散項: 自由流感度を大幅に低減
2. Stress Limiter: $\mu_t = \min\left(\frac{\rho k}{\omega}, \frac{\rho a_1 k}{\Omega F_2}\right)$ に類似したリミッター
3. 渦伸張補正: $\beta$ に渦度テンソルの不変量を含める
これらの改良によりSSTとの性能差は縮小したが、SSTの実績と普及度の前では巻き返しは難しい。
ω方程式の理論的意味
なぜ $\varepsilon$ じゃなくて $\omega$ を使うんですか？
歴史的には、Kolmogorov (1942) が最初に $\omega$（「乱れの周波数」）を提案した。$\omega = \varepsilon/(C_\mu k)$ は乱流の時間スケール $\tau = 1/\omega$ の逆数であり、壁面で有限値をとる（$\varepsilon$ は壁面で特異）。
物理的解釈:
- $k$: 乱流エネルギーの「大きさ」
- $\omega$: 乱流エネルギーの「散逸速度」（乱流の回転周波数）
- $\varepsilon$: 乱流エネルギーの「散逸量」
$\omega$ は壁面での振る舞いが $\varepsilon$ より好ましいが、自由流では $\varepsilon$ の方が安定。SSTはこの両者の良いところ取りだ。
圧縮性への拡張
超音速流れでk-ωは使えますか？
Wilcoxは圧縮性補正を提案している。高マッハ数での乱流圧縮性効果（dilatation dissipation）:
$M_t = \sqrt{2k}/a$ は乱流マッハ数、$a$ は音速。$M_t > M_{t0} \approx 0.25$ で圧縮性効果が発現する。超音速混合層の拡がり率抑制などに効果がある。
k-ωの歴史的・理論的な背景を知っておくとSSTの理解も深まりますね。
その通り。SSTの各項の意味を正しく理解するにはk-ωとk-εの両方の知識が必要だ。
k-ωの自由流感受性——悩ましい「無限遠の影響」
標準k-ωの最も知られた弱点は「自由流感受性」です。入口境界でのω（比消散率）の値をわずかに変えるだけで、離れた位置の境界層解への影響が出てしまう。これは物理的にはおかしい話で、「無限遠の静止した流体の乱流特性が壁面近傍に影響する」という非現実的な特性です。Wilcoxが1998年と2006年の改訂でこの問題を修正しようとしましたが、完全には解決していません。SSTモデルが生まれた動機の一つもこの問題を避けることで、k-ωの歴史的弱点がSSTという傑作を生んだとも言えます。
標準k-ωモデル（Wilcox）のトラブル対応
トラブルシューティング
標準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-ωユーザーの間でよく聞く話です。
関連トピック
なった
詳しく
報告