DES（Detached Eddy Simulation）—
理論と物理
概要
先生、RANSは全部モデル化して、LESは大スケールを直接解くって学んだんですが、DESはその中間ですか？
まさにそう！ DESは「壁近傍の境界層はRANS（計算コスト節約）、剥離した後流・渦はLES（高精度）」というハイブリッド戦略だ。RANSだけでは剥離・後流の非定常構造を正確に捉えられないが、全域LESは計算コストが膨大になる。DESはこのジレンマを実用的に解決する手法として、Spalart（1997年）が「bluff bodyの後流解析」を念頭に提案した。自動車・航空機・建築物の空力解析で広く使われている。
「剥離した後流」だと何が起きていて、なぜRANSだとダメなんですか？
例えば自動車のトランク後部では流れが剥離して大きな渦（Von Karman渦列に近い非定常構造）が形成される。RANSはこれを時間平均した「定常的な渦」としてしか表現できないから、非定常の圧力変動・空力騒音・揚力変動が正確に予測できない。DESのLES領域ではこの大スケール渦を直接解像する。
DESはどんな流れ現象に適していて、逆に苦手な問題はありますか？
DESが得意なのは「大きな分離渦・ウエイクが支配的な流れ」——ブラフボディ（円柱・角柱）後流、翼失速後流、自動車後部ウエイク、建物周囲風など。逆に苦手なのは「境界層が剥離せず薄いままで維持される」チャネル流・管路内流れや、剥離がほとんどない高アスペクト比翼。これらはRANSまたは全域LESの方が向いている。
DES97の定式化
DES97の数式を教えてください。
DES97はSpalart-Allmarasモデルをベースにしている。S-Aモデルの乱流長さスケール $d_w$（壁面距離）を次のように置き換えるだけだ：
$\Delta = \max(\Delta x, \Delta y, \Delta z)$ はメッシュ幅（最大辺長）。$C_{DES} \approx 0.65$（S-Aベース）。$d_w \ll C_{DES}\Delta$（壁近傍）では $\tilde{d} = d_w$ でRANS動作。$C_{DES}\Delta \ll d_w$（後流域）では $\tilde{d} = C_{DES}\Delta$ となりSmagorinskyLESと同等になる。これが境界層→後流へのスムーズな切り替えを実現する。
$C_{DES}$ ってキャリブレーションで決まるんですか？ デカルビークが変わると結果も変わりますか？
$C_{DES}$ は一様等方性乱流（HIT）のDNSデータとのキャリブレーションで決まる。S-Aベースは0.65、SST-k-ωベースは0.61が標準値。この値を大きくするとLESモードへの切り替えが遅くなり（境界層がより長く保護される）、小さくすると早くLESになる。実務では通常デフォルト値を使うが、非常に複雑な形状では感度テストが有用だ。
RANS/LES切り替え領域
切り替わる場所は自動で決まるんですか？ ユーザーが制御できますか？
基本的に自動だが、メッシュサイズ $\Delta$ が切り替えの主要パラメータだから間接的にユーザーが制御できる。壁近傍のメッシュを粗くすれば $C_{DES}\Delta$ が大きくなって切り替わりにくく（RANSを維持）、後流領域のメッシュを細かくすれば $C_{DES}\Delta$ が小さくなってLESモードになりやすい。ただしこの「メッシュで制御」という性質がDES97の弱点でもある——境界層内で意図せずLESモードに切り替わるMSD問題の原因だ（DDESで改善）。
$\Delta$ の定義は $\max(\Delta x, \Delta y, \Delta z)$ と言いましたが、他の定義もありますか？
いくつか提案されている。$\Delta_{max} = \max(\Delta x, \Delta y, \Delta z)$（DES97の標準）のほか、体積ベースの $\Delta_v = (\Delta x \Delta y \Delta z)^{1/3}$ も使われる。さらに近年はAnisotropic Measure $\Delta_\omega = 2\max(d_i \times d_j)$（渦量方向に基づく長さスケール）も提案されており、非等方性格子でのLES領域精度を改善する。OpenFOAMではcubeRootVolやmaxLengthが選択可能。
数値解法と実装
メッシュ要件
DES解析のメッシュ設計でよく使われる指針はありますか？
DESのメッシュ設計でよく使われる指針を示すよ。
| 領域 | 推奨y+ | 格子タイプ | 注意点 |
|---|---|---|---|
| RANS境界層 | y+ ≈ 1〜5 | プリズム層（壁面平行） | S-Aは y+≈1、SSTは壁面関数可 |
| RANS→LES遷移 | — | 均一サイズへ滑らかに移行 | 成長率 ≤ 1.2 |
| LES後流域 | — | 等方性格子（立方体に近い） | 非等方格子だとSGSモデルが劣化 |
SGSモデルとDESの関係
DESのLES領域はSmagorinskyモデルと同等と言っていましたが、SGS（Sub-Grid Scale）モデルはどうなっているんですか？
DES97のS-AベースではLES領域でのSGS粘性は：
$|S|$ は歪み速度テンソルのノルム。これはSmagorinskyモデル $\nu_{sgs} = (C_S \Delta)^2 |S|$（$C_S \approx 0.1$〜0.2）と同じ形式で、$C_{DES} \approx 0.65$ が $C_S$ に相当する。SST-DDESのLES領域では $k$-方程式の消散項の係数がSGSモデルの役割を果たす形になっており、「自然に」LESのSGSエネルギーを処理できる。
つまりDESを使えば別にSGSモデルを設定しなくていいということですか？
そう、それがDESの大きな利点のひとつだ。RANSモデルの中にSGSが自然に内包されている形になっている。動的Smagorinskyや壁適応局所渦粘性（WALE）モデルのような独立したSGSモデルを設定する必要はない。ただしこれは同時に、DESのLES精度が「埋め込まれたSmagorinskyレベル」に制限されることも意味する——より精密なSGSモデルが必要な場合はフルLESを使う必要がある。
実践ガイド
建築物の風工学解析
建築物の風荷重解析でDESを使う場合のポイントを教えてください。
建築外部環境（風工学）もDESの重要な応用分野だよ。
- 入口境界：大気境界層プロファイル（対数則）を入口に設定。乱流変動はSynthetic Turbulence Generator（STG）やSEM（Synthetic Eddy Method）で追加。
- 建物周囲メッシュ：建物高さに対して$\delta / H \approx 0.05$以下の等方性格子。
- 壁面処理：建物表面は壁面関数（y+ 30〜300）の方が現実的。プリズム層は地面（境界層底面）に設ける。
- 統計量：平均圧力係数 $C_p$ に加え、ピーク $C_p$ の確率分布を取得。風荷重の極値評価には十分な統計時間が必要（実時間で10分以上相当）。
実時間で10分相当って、計算ステップ数で言うとどのくらいになりますか？
スケールによって違うが、例として高さ100mのビル、風速10m/sで $\Delta t = 0.01$s とすると、10分 = 600s で 6万ステップが必要。これは計算機で数日〜数週間に相当する。実務では風洞試験との比較検証を先に行い、DESの計算スコープを絞り込んでから本計算に移る。建築規準での採用には検証ドキュメントが必須な国も多い。
実践チェックリスト
- 境界層 y+ をRANSモデル要件（S-A: y+≈1、SST: y+≈1または壁面関数）に合わせて設定
- 後流域メッシュを等方性（アスペクト比≤5）に設計
- RANS→LES遷移領域の格子成長率 ≤1.2
- LES移流スキームをLinear/LUSTに変更
- 非定常タイムスキームを2次精度以上に設定
- LES域のCourant数 ≤1 を確認（定期モニタリング）
- 初期場はRANS定常解から取得（warm up）
- 統計収集前に2L/U以上のウォームアップ実施
- 統計量の収束（平均CD・CL変動 <±1%）を確認してから本採用
ソフトウェア比較
DESを使えるツールと選択肢を教えてください。
主要ツールを比較してみよう。
| ツール | DES実装 | ベースモデル | 特徴 |
|---|---|---|---|
| Ansys Fluent | DES97, DDES, IDDES | S-A, SST k-ω | 豊富な乱流モデル選択肢、GUIセットアップ容易 |
| OpenFOAM | SpalartAllmarasDES, kOmegaSSTDES | S-A, SST k-ω | オープンソース、カスタマイズ自由、DDES/IDDES搭載 |
| Simcenter STAR-CCM+ | DES97, DDES | SST k-ω | 後処理・可視化が充実、DES切り替え領域の視覚化容易 |
| SU2 | SA-DES, SA-DDES | S-A | 航空宇宙最適化に特化、オープンソース |
| PowerFLOW（Exa） | VLES相当 | LBM基礎 | 自動車空力に特化、格子ボルツマン法ベース |
先端技術
DES系の最新発展と、将来どうなっていくかを教えてください。
DESファミリーの発展の流れをまとめると：DES97（1997）→DDES（2006）→IDDES（2008）→最新のSLA（Stress-Limited Approach, 2023）。IDDESではWMLESとDDESを統合し、境界層の薄い領域でも正確な乱流輸送が得られるようになった。ML（機械学習）との統合では、RANSとLESのブレンド係数をデータドリブンで決定する研究が2020年代から急増している。
機械学習との組み合わせで具体的にどんなことができるようになっていますか？
いくつかのアプローチがある。
- Physics-Informed Neural Network（PINN）：DESの時間平均場をNNで補間し、計算コストを1/10以下に削減する研究。
- 強化学習ベースのタイムステップ制御：Courant数と精度のトレードオフをRLで自動最適化。
- データ同化（Data Assimilation）：風洞・LiDAR計測データとDES結果を融合して高精度予測。
- 代替モデル（Surrogate Model）：DES計算の設計パラメータ感度をNNで学習し、設計探索を高速化。特に自動車・航空機設計で実用化が進む。
DESが生まれた背景
DESはPhilippe Spalartが1997年のAIAA論文「Comments on the Feasibility of LES for Wings, and on a Hybrid RANS/LES Approach」で提案した。当時の計算機ではRANSは実用的だったが、翼や車体の複雑な剥離流れの予測精度が問題になっていた。「LESを使いたいが計算コストが10倍以上かかる——なら壁近くだけRANSにすれば」というシンプルな発想から生まれた画期的な手法だ。
トラブルシューティング
DESでよくある問題はどんなものがありますか？
MSD（早期剥離）、グレーゾーン問題、LES領域での乱流変動不足、統計量の収束不足などが代表的なトラブルだ。
| 症状 | 主な原因 | 対策 |
|---|---|---|
| 早期境界層剥離 | MSD問題（DES97固有） | DDESまたはIDDESに変更、または境界層メッシュを粗くしてLES侵入を防ぐ |
| 後流に渦が現れない | グレーゾーン問題・LES域メッシュ粗すぎ | 後流メッシュを細化、LUS/LUSTスキームへ変更 |
| 揚力・抗力が収束しない | 統計収集時間不足 | 統計収集を10L/U以上に延長 |
| 入力変動がLES領域に入らない | 入口STG不使用 | Synthetic Turbulence Generatorを入口に追加 |
MSD問題、グレーゾーン問題、統計収束不足、境界条件設定など詳細解説