LESの基礎理論 -- トラブルシューティングガイド

LES特有の問題をいくつか見ていこう。RANSとは性質が異なる問題が多いんだ。

考えられる原因と対策:

• 数値散逸が大きすぎる: 1次上流差分を使っていないか確認。運動量方程式には中心差分系のスキームを使うこと
• メッシュが粗すぎる: LES_IQを確認して0.8未満なら格子を細分化
• 入口で乱流変動を与えていない: SEM等の変動生成法を適用すること
• SGS粘性が過大: Smagorinskyモデルの定数 $C_s$ が大きすぎないか確認

対策:

• CFL数の確認: 最大CFL数が1.0を超えていないか。局所的に高CFL領域がないかチェック
• 初期過渡の取り扱い: RANS解から初期化した直後は残差が大きく発散しやすい。最初の数百ステップは時間刻みを小さくするか、PIMPLE外側ループの反復回数を増やす
• メッシュ品質: 非直交性やスキューネスが大きいセルが存在すると不安定化する。特に非直交補正（non-orthogonal corrector）の回数を増やすこと
• 中心差分の不安定性: 純粋な中心差分で不安定な場合は、微量のblending（5-10%程度の上流差分成分）を加える

サンプリング時間が不十分な場合がほとんどだ。以下をチェックしてほしい。

• 流れ場が統計的に定常な状態（statistically stationary）に達してからサンプリングを開始しているか
• サンプリング時間が主要な渦の特徴的時間スケールの10倍以上あるか
• 低周波数の変動（流れの首振りなど）がある場合、さらに長いサンプリングが必要

壁面解像LESで壁面が合わない場合は、

• 壁面垂直方向の格子解像度: 第一セルの $y^+$ が1.0以下か確認
• スパン方向解像度: $\Delta z^+ < 40$ を満たしているか。壁面近傍の縦渦対を解像するために不可欠
• 流れ方向解像度: $\Delta x^+ < 100$ を満たしているか
• 壁面接線方向の等方性: 壁面近傍で過度にアスペクト比が大きいセルは避ける

WMLESを使っている場合は、マッチング位置の $y_{match}$ が適切かを確認しよう。

• DES/IDDESへの切り替え: 壁面近傍のコスト削減に有効
• 壁面モデルLES（WMLES）の適用: 壁面解像不要でコスト大幅削減
• 流れ方向に周期境界条件: 完全発達流れなら使える
• 対称条件の活用: 平均流が対称な場合に限り有効
• 並列計算の最適化: 領域分割のバランスを確認。OpenFOAMのdecomposeParでscotchやhierarchical法を試す

その通り。LESは「ゴミを入れてゴミを出す（garbage in, garbage out）」の原則がRANS以上に顕著だ。設定の一つ一つが結果に直結するので、丁寧なセットアップと品質評価を心がけてほしい。