EDCモデル（Eddy Dissipation Concept）

Eddy Dissipation Conceptの略で、MagnussenがEddy Dissipation Model（EDM）を発展させた乱流燃焼モデルだ。EDMは化学反応が無限に速い仮定だったが、EDCは有限速度の詳細化学反応機構を乱流場で扱えるよう拡張したものだ。

そうだ。EDMでは $\dot{\omega} = A\,\rho\,\frac{\varepsilon}{k}\min(Y_F, Y_O/s)$ のように反応速度が乱流混合で決まり、Arrhenius速度論を無視する。これはDamkohler数が大きい（反応が十分速い）場合には妥当だが、CO酸化やNOx生成のような有限速度反応では不正確だ。EDCはこの限界を克服する。

EDCでは乱流場の微細構造（fine structure）内で化学反応が進行すると考える。微細構造の体積分率 $\xi^$ と滞留時間 $\tau^$ は、乱流の $k$, $\varepsilon$ から以下のように求める。

鋭い指摘だ。$\xi^$ はKolmogorovスケールの体積分率に対応し、$\tau^$ はKolmogorov時間スケールのオーダーだ。物理的には「乱流の最小渦の中で化学反応が進行する」というイメージだ。

化学種 $i$ の平均反応速度は次のように書ける。

微細構造内で定容の0Dリアクターを $\tau^*$ だけ時間積分して求める。この0D積分にCVODE等のStiff ODEソルバーが使われる。つまりEDCの計算コストの大部分はこの0D化学反応積分にある。

そのとおり。乱流と化学反応の相互作用を物理的に明快なモデルで表現している点がEDCの強みだ。

EDCの数値コストは化学反応ODE積分が支配的だ。各CFDセルで毎反復（RANS）または毎タイムステップ（LES）に0Dリアクターを解く必要がある。

目安を示そう。

だからこそISATとの併用が必須なんだ。ISATを使えばGRI-Mech 3.0でも実用的な時間で回せる。FluentではEDC + Stiff Chemistry Solver + ISATの組み合わせがデフォルトの推奨設定だ。

1. Models > Species > Species Transport を有効化

OpenFOAMの reactingFoam ソルバーで combustionProperties にEDCを指定する。

OpenFOAMのEDC実装はv2005（Magnussen 2005改訂版）にも対応している。改訂版ではReynolds数依存の $\xi^*$ 補正が入り、低Re数領域での精度が改善されている。

$C_\xi$ を大きくすると微細構造体積が増え、反応速度が上がる。$C_\tau$ を大きくすると滞留時間が延び、やはり反応が進む。通常はデフォルト値で十分だが、火炎リフトオフ高さの調整で$C_\xi$を±20%程度変える研究報告はある。ただしこの調整はケース依存であり、万能な推奨値はない。

そうだ。ISAT未使用のEDC計算は研究用途でも非常に時間がかかる。ISATの精度設定と合わせてチューニングしよう。

EDCは以下のようなケースで威力を発揮する。

判断基準はDamkohler数だ。

LESではサブグリッドの乱流-化学反応相互作用が小さいから、EDCの微細構造仮定がやや過剰になる。PaSR（Partially Stirred Reactor）モデルのほうがLESに適している。OpenFOAMの reactingFoam ではPaSRが標準搭載だ。

そうだ。cold flow → EDM → EDC の3段ロケットで確実に収束させよう。

EDCは広く実装されているが、細かい仕様が異なる。

SAGE（Stochastic Adaptive Grid Engine）はEDCとは異なるアプローチだ。各セルでArrhenius反応速度を直接評価し、乱流-化学反応相互作用は「Well-Stirred Reactor限界」で処理する。EDCの微細構造概念はないが、実用的にはEDCと同程度の結果を与える場面が多い。

そのとおり。NOx規制対応やCO排出管理が求められる工業燃焼器の設計には、EDCは非常に信頼性の高い選択だ。

EDCは1981年のMagnussenの原論文以来40年以上の歴史があるが、近年も改良が続いている。

Magnuessen自身が2005年にReynolds数依存の修正版（EDC v2005）を提案した。低Re数領域（層流-乱流遷移域）でのオリジナルEDCの過大評価を補正するものだ。さらにParente（2016）はDNSデータに基づいて $C_\xi$ と $C_\tau$ の最適値を再評価し、MILD燃焼条件では標準値からの調整が必要なことを示した。

Moderate or Intense Low-oxygen Dilution燃焼の略で、高温排ガス再循環により酸素濃度が低い条件で燃焼する技術だ。温度が均一化されNOxが大幅に低減する。EDCは微細構造内のPSR仮定があるためMILD条件に比較的適しているが、標準定数ではReactedness（反応進行度）を過大評価する傾向がある。

EDCの最大のボトルネックである0Dリアクター積分をニューラルネットワークで置き換える研究が進んでいる。Blanchoらは、DNNをISATの代替として使い、5-10倍の高速化を達成している。ただし汎化性能と数値安定性が課題で、production-readyにはまだ至っていない。

原理的には使えるが、LESではフィルタ幅がKolmogorovスケールより大きいという仮定が成り立たない場合がある。そこでフィルタ幅依存のEDC（filtered EDC）が研究されている。Fureby（2009）はフィルタ幅に応じて$\xi^$と$\tau^$を修正するLES-EDCを提案した。

そうだ。理論がシンプルで拡張しやすい点がEDCの長所だ。MILD燃焼やアンモニア燃焼など新しい燃焼形態への適用拡大が今後も続くだろう。

EDCは強力だがトラブルも起きやすい。代表的な問題を整理しよう。

症状: EDMでは正常な火炎温度が出ていたのに、EDCに切り替えると温度が数百K低下する。

原因: EDCの微細構造体積 $\xi^*$ が小さすぎて反応速度が不足している。特に壁面近傍や乱流強度が低い領域で顕著。

対策:

そうだ。EDCは乱流モデリングと化学反応積分の二重の複雑さを持つから、問題の切り分けが重要だ。まず非反応流で流れ場を確立し、EDMで着火を確認してからEDCに移行する手順を厳守しよう。