Spalart-Allmarasモデル — トラブルシューティングガイド
より充実した内容を spalart-allmaras.html でご覧いただけます。
よくある問題と対策
SAモデルで計算がうまくいかないとき、何を確認すればいいですか？
1. nuTilda が負値になる
症状: 計算中に $\tilde{\nu}$ が負値に発散し、NaN/Infが発生
原因: 初期条件が不適切（$\tilde{\nu}$ の初期値がゼロまたは非常に小さい）、メッシュ品質が悪い、逆流が入口境界で発生
対策:
- 初期値を $\tilde{\nu} = 3\nu$ 〜 $5\nu$ に設定
- SA-negバリアントを使用（Fluentでは自動対応）
- 入口の $\tilde{\nu}$ 値を上げる（$10\nu$ 程度に増やして安定化後に戻す）
- メッシュ品質指標（スキューネス < 0.85）を確認
2. 壁面距離の計算エラー
壁面距離がおかしいとどうなりますか？
症状: 散逸項 $c_{w1}f_w(\tilde{\nu}/d)^2$ が異常値を取り、渦粘性が過大または過小
原因: 複数の壁面が近接する領域（隙間、薄い構造）で壁面距離が正しく計算されない
対策:
- FluentではSolve > Initialize > Wall Distance を再計算
- OpenFOAMでは
wallDistの method をmeshWave（デフォルト）からPoissonに変更してみる - 隙間が非常に狭い場合は形状を簡略化
3. 自由せん断層の過大予測
噴流のスプレッド角がおかしいんですけど。
症状: ジェットの拡がり率が実験より20〜30%大きい
原因: SAモデルは自由せん断流に対してキャリブレーションされていない。壁面距離 $d$ が大きい領域で散逸が不足し、渦粘性が過大になる
対策:
- 噴流・混合層の解析にはSST k-omegaまたはk-epsilon Realizableを使う
- DES/DDESに切り替えて自由せん断層をLESで解く
4. 遷移位置が実験と合わない
症状: 翼面の層流-乱流遷移位置が実験と異なる
原因: 標準SAモデルは完全乱流を仮定。遷移を予測する機能がない
対策:
- 遷移を考慮したいならSA + $\gamma$-$Re_\theta$ 遷移モデルを使う
- Fluentでは Transition SST が利用可能
- trip条件で遷移位置を手動指定する方法もある（NASAケースではよく使われる）
検証用ベンチマーク
SAモデルの結果を検証するのに適したベンチマークは？
NASA Turbulence Modeling Resource (https://turbmodels.larc.nasa.gov/) が最も信頼できる。以下のケースが推奨だ。
| ケース | 検証項目 | 実験データ |
|---|---|---|
| Flat Plate | $C_f$ 分布 | Wieghardt (1951) |
| NACA 0012 | $C_p$, $C_L$-$\alpha$ | Gregory-O'Reilly (1970) |
| RAE 2822 | 遷音速 $C_p$ | Cook et al. (1979) |
| DPW (CRM) | 巡航条件の $C_D$ | NASA CRM実験 |
まずFlat Plateで $C_f$ が理論値（$C_f = 0.058 Re_x^{-0.2}$）と一致するか確認して、それからNACA翼に進むのが王道ですね。
ν̃が負になる問題——OpenFOAMユーザーが最初に踏む地雷
SAモデルの作業変数 ν̃（修正渦粘性）は理論上は非負ですが、数値計算中に負値になることがあります。これが起きると渦粘性が負になり計算が発散します。OpenFOAMのSpalartAllmarasソルバーでは、負値を対処するSA-neg変形が導入されましたが、古いバージョンではこの対処が不十分で初心者が頻繁に詰まりました。「OpenFOAMのSAでnu_tがおかしい」という質問はCFDフォーラムの定番トピックで、最初のセルのy+を正しく設定してν̃の初期値を適切にすることが解決への近道です。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——Spalart-Allmarasモデルの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告