ラージ・エディ・シミュレーション (LES) — CAE用語解説
LESとは — 空間フィルタリングの考え方
いい質問だ。RANSは乱流の「全スケール」を時間平均してモデル化する。つまり渦の情報はすべてモデルに押し込む。一方LESは、大きな渦（エネルギーを多く持つスケール）はNavier-Stokes方程式で直接計算して、メッシュで解像できない小さな渦だけをモデル化する。だからLESは「大きい渦は直接、小さい渦だけモデル」というハイブリッドな考え方だね。
「大きい渦」と「小さい渦」をどうやって分けるんですか？
それが空間フィルタリングだ。Navier-Stokes方程式にフィルタ関数 $G$ を畳み込み積分（コンボリューション）して、格子幅 $\Delta$ より大きい成分だけを取り出す。数式で書くとこうなる：
フィルタ後のNS方程式にはSGS応力テンソル $\tau_{ij}^{\mathrm{sgs}} = \overline{u_i u_j} - \bar{u}_i \bar{u}_j$ が現れる。これが未知量になるから、何らかのモデルで閉じてやる必要がある。実際にはメッシュ自体がフィルタの役割を果たすことが多くて、「暗黙的フィルタリング（implicit filtering）」と呼ばれるよ。
SGSモデル（Smagorinsky・WALE・動的モデル）
SGS応力をモデル化する方法にはどんな種類があるんですか？
代表的なものを3つ紹介しよう。
1. Smagorinskyモデル（1963年）
最も古典的なSGSモデルだ。SGS渦粘性を次のように定義する：
$C_s$ はSmagorinsky定数で、等方性乱流だと約0.17、チャネル流れだと0.1くらいに下げることが多い。シンプルで安定だけど、壁面近くで渦粘性が過大になる欠点がある。壁面ではvan Driest減衰関数 $f_d = 1 - \exp(-y^+/A^+)$ をかけて補正するのが一般的だ。
壁の近くで渦粘性が大きすぎるってことは、乱流が実際より強く拡散されちゃうってことですか？
そのとおり。壁面近くでは乱流の渦は小さくて弱いのに、Smagorinskyモデルは $|\bar{S}|$ が大きい（速度勾配が急だから）という理由で渦粘性を大きくしてしまう。これを改善したのが次のモデルだ。
2. WALEモデル（Wall-Adapting Local Eddy-viscosity, Nicoud & Ducros 1999）
速度勾配テンソルのトレースフリー対称部分を使う：
ここで $\mathcal{S}_{ij}^d$ は速度勾配テンソルの二乗のトレースフリー対称部分。$C_w \approx 0.325$ が標準値だ。WALEの利点は、壁面で $\nu_{\mathrm{sgs}} \to 0$ が自然に達成されること、つまりvan Driest減衰関数が不要なんだ。層流領域でもゼロに落ちるから使いやすい。
3つ目の動的モデルというのは？
3. 動的Smagorinskyモデル（Germano et al. 1991, Lilly 1992）
から $C_s^2$ を局所的に求める。流れ場に応じてモデル定数が変わるから、剥離流れでも後流でも壁面近くでも適切な値になる。ただし計算が不安定になりやすいので、空間方向に平均を取ったり、$C_s^2$ が負にならないようクリッピングしたりする工夫が必要だ。例えばOpenFOAMだと dynamicKEqn や dynamicLagrangian として実装されている。
LESのメッシュ要件と計算コスト
LESをやるとき、メッシュってどのくらい細かくすればいいんですか？ RANSとは全然違いますか？
全然違う。LESでは渦を直接解像しないといけないから、壁面近くのメッシュ要件がかなり厳しい。目安としてはこんな感じだ：
- 壁面垂直方向：$y^+ \approx 1$（第一層セル）
- 流れ方向：$\Delta x^+ \approx 50 \sim 100$
- スパン方向：$\Delta z^+ \approx 15 \sim 40$
例えば自動車の外部空力で $\text{Re}_L = 10^6$ くらいだと、RANS なら数百万セルで済むところ、壁面解像LES（WRLES）だと数億セルが必要になる。しかも3次元の非定常計算で、時間ステップも $\text{CFL} \le 1$ で回さないといけないから、計算コストはRANSの100〜1000倍にもなる。
え、1000倍!? それでも使う価値があるんですか？
あるよ。RANSだと「時間平均された」流れ場しかわからないけど、LESなら渦のダイナミクスが見える。例えば自動車の後方乱流域での非定常的な圧力変動、建物まわりの風荷重のピーク値、燃焼器内の火炎とスワール流の相互作用——こういう「非定常性が本質的な現象」はRANSでは原理的に捉えられない。航空機エンジンの燃焼CFDではLESがほぼ標準になっているし、風工学でも急速に普及してるよ。
LES vs RANS vs DNS — 使い分けの判断基準
DNS（直接数値シミュレーション）もありますよね。LESとDNSはどう違うんですか？
DNSは乱流の全スケール——最大渦からKolmogorovスケール $\eta = (\nu^3/\varepsilon)^{1/4}$ まで——をすべてメッシュで解像する。モデルは一切使わない。だから最も正確だけど、メッシュ数は $N \propto \text{Re}^{9/4}$ でスケールするから、工業的な高Reynolds数流れには事実上適用不可能だ。
整理するとこうなる：
| 手法 | モデル化範囲 | メッシュ数スケーリング | 実用的なRe上限 |
|---|---|---|---|
| DNS | なし（全スケール解像） | $\sim \text{Re}^{9/4}$ | $\sim 10^4$ |
| LES | SGSスケールのみ | $\sim \text{Re}^{13/7}$（壁面解像） | $\sim 10^6$ |
| RANS | 全乱流スケール | Re依存小 | 制限なし |
じゃあ、実務ではどうやって使い分けるんですか？ いつLESを選ぶべきなんだろう。
判断基準をシンプルにまとめるとこうだ：
- RANSで十分：定常的な流れ、添付流れの抗力・揚力予測、パイプ内の圧損計算など。コストが圧倒的に安い。
- LESが必要：大規模剥離、後流の非定常変動、混合・燃焼、空力騒音（エアロアコースティクス）、風工学での突風応答など。
- DNSが必要：乱流の基礎研究、モデルの検証データ生成、遷移現象の詳細解明。学術目的がメイン。
実務では「RANSで試して精度が不足 → LESに切り替える」というステップが現実的だ。最初からLESで始めるとコストが読めなくて苦労する。
DES/DDES — LESとRANSのハイブリッド手法
壁面のメッシュがそんなに大変なら、壁の近くだけRANSで、離れたところはLESで計算する方法ってないんですか？
まさにその発想で生まれたのがDES（Detached Eddy Simulation）だ。Spalart（1997年）が提案した手法で、壁面近くではRANSモデル（SA, SST等）として振る舞い、壁から離れた領域ではLES的に渦を解像する。切り替えはメッシュサイズと壁面距離で自動的に判定される。
ただし初期のDESには問題があった。「グレーゾーン問題」といって、RANS領域とLES領域の境界で、メッシュが中途半端に細かいとRANSの渦粘性が不十分に下がり、LESとしても渦を解像できない「どっちつかず」の領域ができてしまう。
それを解決する方法はあるんですか？
改良版としてDDES（Delayed DES）とIDDES（Improved DDES）がある。DDESはSpalart et al.（2006年）が提案したもので、境界層内部ではRANSモードを強制するシールド関数 $f_d$ を導入した：
これにより、壁面添付流の境界層がうっかりLESモードに切り替わる問題（MSD: Modelled Stress Depletion）を防げる。IDDESはさらに壁面モデルLES（WMLES）の機能も統合していて、壁面近くのメッシュが十分細かければWMLESとして動く。
OpenFOAMでは SpalartAllmarasDDES や kOmegaSSTDDES として簡単に使える。FLUENTでもDDESはサポート済みだ。実務では「とりあえずRANS → 精度不足ならDDES → それでも足りなければ壁面解像LES」というステップアップが王道だよ。
実務でのLES運用Tips
実務でLESを使うとき、注意すべきポイントを教えてもらえますか？
よく聞かれるので、重要なものをまとめるよ。
- 初期条件と助走時間：LESは非定常計算だから、統計的に定常になるまでの「助走時間（wash-out time）」が必要。流れが領域を数回通過する時間（フロースルータイム）の5〜10倍は回してから統計を取り始める。
- 入口境界条件：乱流の入口条件が非常に重要。一様流を入れると境界層が発達するまで長い助走区間が必要。Synthetic Eddy Method（SEM）やRecycling/Rescalingで人工的な乱流変動を与えると効率的だ。
- 数値スキーム：中心差分系が基本。風上差分は数値拡散が大きくて渦を潰してしまう。ただし中心差分は不安定になりやすいので、少量のブレンド（5〜10%の風上成分）を入れるか、低拡散スキーム（例：LUST, Gamma）を使う。
- 時間刻み：CFL数を0.5〜1.0以下に保つ。陰解法でもCFL > 5だと時間精度が落ちてLESの意味がなくなる。
- 統計処理：平均値・RMS値・Power Spectral Density（PSD）を計算して、結果が時間窓に依存しないことを確認する。統計が不十分だと「たまたまのスナップショット」を結論にしてしまう危険がある。
すごくよくわかりました。結局、LESは「コストは高いけど、渦のダイナミクスが見えるから非定常現象には不可欠」ってことですね。
そのまとめは完璧だ。もう一つ付け加えると、GPUコンピューティングの進歩でLESのコストは年々下がっている。NVIDIAのAmgXやPETScのGPU対応、Lattice Boltzmann法ベースのGPUソルバー（PowerFLOW, XFlow, ProLB等）も出てきて、自動車OEMの空力開発では既にLES系手法が標準になりつつある。5年後にはRANSよりLESが主流になっている分野も増えるだろう。
CAE用語の正確な理解は、チーム内のコミュニケーションの基盤です。 — Project NovaSolverは実務者の学習支援も視野に入れています。
LESの実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
お問い合わせ（準備中）関連トピック
なった
詳しく
報告