Reynolds応力モデル（RSM） — トラブルシューティングガイド
より充実した内容を reynolds-stress-model.html でご覧いただけます。
よくある問題と対策
RSMで計算がうまくいかないとき、何をチェックすればいいですか？
1. 計算が発散する
症状: 残差が発散、NaN/Infが発生
原因: RSMの初期値が不適切（ゼロ初期化した等）、Under-relaxationが高すぎる、メッシュ品質不良
対策:
- k-epsilon Realizableで収束解を得てからRSMに切替える（最重要）
- Under-relaxationを下げる: Reynolds Stresses = 0.3、Turbulent Dissipation Rate = 0.3
- Fluentでは Stabilization Method を "Implicit Body Force" にする
- メッシュのスキューネスを0.85以下に抑える
2. レイノルズ応力の実現可能性違反
実現可能性違反ってなんですか？
物理的に、レイノルズ応力テンソルの固有値は非負でなければならない（$k \geq 0$、各成分の分散は非負）。数値誤差でこの条件が破れると非物理的な解になる。
対策:
- Fluentでは Realizability Enforcement を有効にする
- SSGモデルはLRRモデルより実現可能性を保ちやすい
- ソース項の陰的線形化が適切か確認
3. 結果がk-epsilonと変わらない
RSMに切り替えたのに、k-epsilonとほとんど同じ結果になるんですが。
原因: 流れの異方性が弱い（単純な管路流、完全発達流れ等）。RSMの利点が出るのは旋回流や曲がり管のような異方性が強い流れだ。
対策:
- まず異方性指標 $A = 1 - \frac{27}{2}\det(a_{ij}^*)$ をポスト処理で確認。$A$ が0に近ければ等方的で、RSMの優位性がない
- 単純な流れなら k-epsilon 系に戻してコストを節約する
4. 壁面近傍で不安定
症状: 壁面近傍でレイノルズ応力が振動し、収束しない
対策:
- $\varepsilon$ベースのRSMから$\omega$ベースのBSL-RSMに変更
- 壁面処理をEnhanced Wall Treatment（$y^+ \approx 1$）に変更
- 壁面反射項のモデル（Gibson-Launder）が不適切な場合は無効にしてみる
RSMが必要かの判定方法
そもそもRSMが本当に必要かを事前に判断する方法はありますか？
以下のステップで判断できる。
1. k-epsilon Realizable で計算を実行
2. 渦粘性比 $\mu_t/\mu$ の分布を確認
3. ストレインレートとローテーションレートの比 $S/\Omega$ を計算。$S/\Omega$ が1から大きくずれる領域があればRSMの検討価値がある
4. 実験データがある場合、二次流れパターンがk-epsilonで再現できているか確認
闇雲にRSMを使うのではなく、渦粘性仮説の妥当性を先に確認するのが正しいアプローチですね。
壁面反射項が引き起こす謎の発散——RSMトラブルの深い沼
RSMのトラブルシューティングで最も厄介なのが圧力ひずみ相関の壁面反射項 $\Pi_{ij,w}$ による発散です。この項は壁に近いセルで大きな値を取り、粗いメッシュや複雑なジオメトリでは数値的に暴れることがあります。対処法は反射項をゼロにするLRRモデルの簡略版に一時切り替えることや、壁面反射係数を0.3から0.1に下げることですが、これは精度低下と引き換えです。「RSMが発散したら壁面反射項を疑え」は、CFDサポートエンジニアの間で語り継がれる診断の第一ステップです。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——Reynolds応力モデル（RSM）の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告