EDCモデル（Eddy Dissipation Concept） — トラブルシューティングガイド

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

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

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

対策:

• 乱流強度の入口境界条件を確認（$k$, $\varepsilon$ が低すぎないか）
• EDC定数 $C_\xi$ をやや大きくする（2.14 → 2.5程度で試行）
• EDC v2005モデルに切り替えて低Re数補正を有効化

• テーブルのhit rateが低い: retrieve の割合がmonitor出力で確認可能。80%以下なら組成空間が広すぎる。反応機構の縮約を検討
• テーブルがメモリ上限に達した: max-storage を増やすかISATをリセットして再構築
• 並列計算でのテーブル分散: 各プロセスが独立テーブルを持つため、プロセス間で重複が多い。Chemistry Agglomerationを有効化

• 段階的切替をしたか: cold flow → EDM → EDC の順に切り替える
• 化学反応ODE内部で発散: Stiff ODEソルバーのエラー許容値を緩める（$10^{-6}$ → $10^{-4}$）
• $\xi^*$ の計算でゼロ除算: $k$ がゼロの領域がないか確認。最小乱流エネルギーのリミットを設定

• "ISAT: no growth in xxx cells": ISATテーブルの成長が停止。テーブルリセットか許容誤差緩和
• "Species clipped in xxx cells": 質量分率が [0,1] 外に。Under-Relaxation低下で改善
• "ODE integration failed in cell xxx": 特定セルでODE積分が収束しない。そのセルの$T$, $p$, $Y_i$を確認し、初期条件に問題がないか調査

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