予混合火炎モデル
予混合火炎の理論基礎
概要
先生、予混合火炎って拡散火炎とどう違うんですか？
予混合火炎は、燃料と酸化剤が燃焼前に十分混合された状態で火炎が伝播する形態だ。ガソリンエンジン、ガスタービンのリーンバーン燃焼器、家庭用ガスコンロが該当する。火炎面は明確な境界を持ち、未燃混合気と既燃ガスを分ける。
火炎の伝播速度が重要なパラメータですね。
そうだ。層流燃焼速度 $S_L$ が予混合火炎の基本パラメータだ。メタン/空気（当量比1.0、常温常圧）で $S_L \approx 0.36$ m/s、水素/空気で $S_L \approx 2.1$ m/s だ。
進行変数 $c$ と支配方程式
予混合火炎のCFDではどんな変数を使うんですか？
予混合火炎では進行変数 $c$ を使って火炎面を追跡する。$c=0$ が未燃混合気、$c=1$ が既燃ガスだ。
ここで $\dot{\omega}_c$ は反応ソース項で、火炎面近傍でのみ非ゼロになる。
ソース項 $\dot{\omega}_c$ はどうモデル化するんですか？
これが予混合火炎モデリングの核心だ。主なアプローチは3つある。
主要な予混合燃焼モデル
| モデル | 原理 | 長所 | 短所 | ||
|---|---|---|---|---|---|
| G-equation (Level Set) | 火炎面を $G=0$ の等値面で追跡 | 幾何学的に明快 | 火炎面内部構造なし | ||
| TFC (Turbulent Flame Closure) | Zimont モデル。$S_T = A(u'/S_L)^n S_L$ | 実装容易 | $S_T$ の経験則依存 | ||
| FSD (Flame Surface Density) | 火炎面積密度 $\Sigma$ の輸送方程式 | 物理ベース | $\Sigma$ 方程式のモデル定数 | ||
| c-equation + 反応速度 | $\dot{\omega}_c = \rho_u S_L | \nabla c | $ | 直接的 | 数値的に火炎厚みが問題 |
乱流燃焼速度
乱流中での火炎速度はどうなりますか？
乱流燃焼速度 $S_T$ は乱流強度 $u'$ とともに増大する。Zimontの相関式が広く使われる。
ここで $\alpha$ は温度拡散係数、$l_t$ は乱流の積分スケール、$A$ はモデル定数（$A \approx 0.52$）だ。
乱流が強いほど火炎面がしわしわになって、見かけの燃焼速度が上がるということですね。
そのとおり。乱流が火炎面をwrinkle（しわ付け）して面積を増やすことで、単位断面積あたりの燃焼速度が増大する。これがDamkohler（1940）の古典的描像で、現代のCFDモデルもこの考え方が基礎にある。
予混合火炎は拡散火炎と違って「火炎面の追跡」が核心なんですね。
そうだ。混合分率ではなく進行変数で火炎を記述する点が根本的に異なる。
「炎の厚さ」は1mm以下——層流燃焼速度の測定がどれほど難しいか
予混合炎の理論で基礎になる「層流燃焼速度SL」は、実はとても測定が難しい量だ。アドバンスト化炎法、対向流法、球形伝播法など複数の実験手法があるが、手法によって同じガスでも10〜20%の値のばらつきが出る。理由は「伸長（strain）の影響をどう補正するか」で研究者ごとに考えが違うからだ。さらに炎の厚さ自体が0.1〜1mm程度しかなく、温度計を刺して測ると計測器が炎を乱す。この「計測そのものが対象を壊す」問題は量子力学の不確定性原理とは別だが、測定の難しさとしては似た構造を持つ。GRI-Mech 3.0の反応定数検証には、こうした慎重な測定データが使われている。
予混合火炎の数値計算手法
数値手法の詳細
予混合火炎をCFDで解くときの数値的な課題を教えてください。
最大の課題は火炎厚みの解像だ。層流予混合火炎の厚みは $\delta_L \approx \alpha/S_L$ で、メタン/空気で約0.5 mm、水素/空気で約0.2 mmだ。これをRANSの数mm単位のメッシュで直接解像することは不可能だ。
Thickened Flame Model（TFM）
それをどう解決するんですか？
LESで広く使われるのがThickened Flame Model（Colin et al., 2000）だ。火炎を人工的に厚くして、メッシュで解像可能にする。
拡散係数を $F$ 倍に増やし、反応速度を $1/F$ 倍にする。
これにより火炎厚さが $F\delta_L$ に増大するが、$S_L$ は変わらない。$F = 5-20$ が一般的だ。
でも、火炎を太くすると乱流との相互作用が変わりませんか？
鋭い指摘だ。厚くした火炎は小スケールの乱流wrinklingを解像できなくなる。そこで効率関数 $E$ を導入して補正する。
Charlette効率関数が代表的で、$E = E(\Delta/\delta_L, u'/S_L)$ の形で与えられる。
Fluentでの実装
Fluentで予混合火炎をどう設定しますか？
Fluentでは以下のモデルが使える。
1. Premixed Combustion (Zimont TFC model): c-equationベース。RANS向け
2. Partially Premixed Combustion: 予混合+非予混合のハイブリッド
3. FGM (Flamelet Generated Manifold): Progress Variable + 混合分率
Zimont TFCモデルの設定:
- Models > Species > Premixed Combustion
- Turbulent Flame Speed model: Zimont
- Laminar Flame Speed: 入力値 or 計算値（当量比依存）
- Flame Stretch Factor: デフォルト0.26
OpenFOAMでの実装
OpenFOAMではどうですか？
XiFoam が予混合燃焼用ソルバーだ。火炎しわ係数 $\Xi$（= $S_T/S_L$）の輸送方程式を解く。
| ソルバー | 対象 | モデル |
|---|---|---|
| XiFoam | 予混合圧縮性 | $\Xi$-equation |
| reactingFoam + PaSR | 予混合/部分予混合 | Species Transport |
| fireFoam | 火災 | EDM/拡散火炎 |
Thickened Flame ModelはOpenFOAMに標準搭載されていますか？
標準ディストリビューションにはないが、コミュニティ版（TFM4OpenFOAMなど）が利用可能だ。ガスタービンLESの研究では広く使われている。
予混合火炎の数値解法は「火炎を太くする」という大胆なアイデアが核心なんですね。
そうだ。TFMは物理的に洗練されたトリックで、LES予混合燃焼のデファクトスタンダードになっている。
「進行度変数c」の正体——予混合炎モデルで最も悩む変数の話
予混合炎モデルの数値実装で多くの人がつまずくのが、反応進行度変数cの定義だ。cは0（未燃）〜1（既燃）を表す変数だが、「どの化学種の質量分率でcを定義するか」はモデルや研究者によって異なる。CO2で定義する流派、温度を規格化して使う流派、複数成分の線形結合で定義する流派——それぞれ微妙に異なる結果になる。Fluentのデフォルトでは生成物の組み合わせが使われるが、燃料や当量比によっては別の定義のほうが精度が出ることもある。「cの定義を変えたら急に合った」という経験談はCAE学会でもたびたび語られる。
予混合火炎の実務適用
実践ガイド
予混合火炎解析の実務手順を教えてください。
ガスタービンのリーンバーン燃焼器やSIエンジンでの予混合火炎解析のフローだ。
解析フロー
1. 層流燃焼速度 $S_L$ の検証 -- Canteraで1D層流予混合火炎を計算し、当量比・圧力・温度依存性を確認
2. 燃焼レジームの確認 -- Peters diagramで条件がwrinkled flamelet / thin reaction zone のどちらに入るか判定
3. モデル選定 -- RANS: Zimont TFC or FGM、LES: TFM or FGM
4. メッシュ設計 -- 火炎帯に十分な解像度（TFM使用時は $F\delta_L$ に5点以上）
5. 着火方法 -- パッチ着火（$c=1$ の領域を設定）またはスパーク着火モデル
6. 後処理 -- 火炎面位置、$S_T$の評価、フラッシュバック/ブローオフのリスク評価
$S_L$ の設定
$S_L$ の設定が重要なんですよね？
そうだ。$S_L$ が10%狂うと $S_T$ も10%以上ずれる。温度・圧力・当量比の3変数に対する$S_L$の関数を正確に設定することが不可欠だ。
| 燃料 | $\phi=1.0$, 1atm, 300K での $S_L$ [cm/s] | 出典 |
|---|---|---|
| CH4/air | 36 | GRI-Mech 3.0 |
| C3H8/air | 39 | USC Mech II |
| H2/air | 210 | Li et al. |
| NH3/air | 7 | Otomo et al. |
アンモニアは7 cm/sですか。メタンの1/5ですね。
だからNH3単独では安定な火炎維持が難しく、H2を10-30%混合して$S_L$を改善する戦略が主流だ。
よくある失敗と対策
| 症状 | 原因 | 対策 |
|---|---|---|
| 火炎がフラッシュバック | $S_T$ が過大 | Zimont定数$A$を下げる、メッシュ細分化 |
| 火炎がブローオフ | $S_T$ が過小 or 混合不良 | $S_L$の設定確認、入口乱流強度の見直し |
| 火炎面が数値的に厚い | メッシュ解像度不足 | TFMの$F$を下げる（メッシュ細分化とセット） |
| 圧力振動（サーモアコースティック） | 火炎-音響連成 | 非反射境界条件の設定、LES推奨 |
サーモアコースティック不安定性もCFDで予測できるんですか？
LESを使えば圧力変動と発熱変動の相関（Rayleigh criterion）を直接評価できる。ガスタービンの燃焼振動問題はLES予混合燃焼の重要な適用分野だ。Fluentの圧縮性ソルバーやOpenFOAMのXiFoamで対応可能だ。
予混合火炎の解析は$S_L$の正確な設定が出発点ということですね。
そうだ。「$S_L$ を制する者が予混合燃焼CFDを制する」と言っても過言ではない。
「点火しない」トラブルの8割は初期条件——予混合炎実践のあるある
予混合炎のCFD実践で最も多い相談が「計算は収束するのにいつまでたっても着火しない」問題だ。Fluent・STAR-CCM+ともにソルバーは「反応なし流れ」として収束してしまい、ユーザーが気づかないことが多い。原因の大半は初期化（Patch）で高温ガス領域を炉の燃焼域に設定していないか、点火セルの温度が活性化エネルギー閾値を超えていないかのどちらかだ。現場の定石は「燃焼域に2000K以上の小さなパッチを当てる→数百ステップ進めて発熱を確認→本計算へ移行」という2段階着火手順。これを知っているかどうかだけで、初日の詰まり方が全然違う。
予混合火炎のソフトウェア比較
商用ツール比較
予混合火炎モデルのツール対応を教えてください。
| ツール | TFC/Zimont | TFM | FGM | G-equation | LES対応 |
|---|---|---|---|---|---|
| Ansys Fluent | 標準搭載 | UDF実装 | R2以降 | なし | あり |
| STAR-CCM+ | 搭載 | 搭載 | 標準推奨 | なし | あり |
| OpenFOAM | XiFoam($\Xi$-eq) | コミュニティ | コミュニティ | コミュニティ | あり |
| CONVERGE | 搭載 | 搭載 | 搭載 | G-equation搭載 | あり |
G-equationが使えるのはCONVERGEだけですか？
CONVERGEはSIエンジン向けにG-equation（Level Set法）を標準搭載している。火炎面をG=0の等値面として追跡し、ノック予測やサイクル変動の評価に使われる。他のツールではG-equationは研究レベルの実装だ。
ガスタービンLES
ガスタービンのLES予混合燃焼にはどれが適していますか？
SIエンジン
SIエンジンの火炎伝播解析は？
ライセンスとコスト
| ツール | 予混合モデルの追加コスト | LES計算の規模感 |
|---|---|---|
| Fluent | 基本パッケージに含む | 500万-2000万セル |
| STAR-CCM+ | Reacting Flowトークン | 1000万-5000万セル |
| CONVERGE | 基本に含む | AMRで実効500万-2000万セル |
| OpenFOAM | 無償 | メッシュ生成は手動 |
予混合燃焼はアプリケーションによって最適なツールが大きく異なるんですね。
そうだ。SIエンジンならCONVERGE、ガスタービンならSTAR-CCM+かOpenFOAM、というのが現状の棲み分けだ。
ANSYS vs STAR-CCM+ の予混合炎比較——どちらが正しい？という問いの難しさ
予混合炎モデルの商用ツール比較でよく話題になるのが「同じ燃焼器でFluentとSTAR-CCM+の結果が10%ずれる」問題だ。原因の大半はモデルのデフォルトパラメータ（乱流シュミット数、火炎伸長モデルの係数など）がベンダーごとに微妙に異なることにある。「どちらが正しい？」という問いは実は意味がなく、「どちらの検証データセットに近い条件か？」を確認するのが正解だ。現場のベテランは「同じ実験データに対してどのツールがフィットしているか確認してからデフォルト値を信用する」と口をそろえて言う。ツールの優劣より、自社の設計条件に合った検証の積み重ねが肝心なのだ。
予混合火炎の先端研究
先端トピックと研究動向
予混合火炎の最先端研究を教えてください。
3つの方向性が注目されている。(1) 水素予混合火炎のLES、(2) サーモアコースティック不安定性の予測、(3) DNSデータに基づくモデル改良だ。
水素予混合火炎
水素の予混合火炎で特有の問題は？
水素はルイス数 $Le \approx 0.3$ と極端に小さいため、差拡散効果で火炎面に強い局所的温度上昇（thermodiffusive instability）が生じる。これにより火炎面がcellular structureを形成し、通常の予混合燃焼モデルの仮定が破れる。
H2予混合火炎のLESでは以下の対応が必要だ。
- TFMの効率関数にルイス数効果を組み込む
- FGMテーブルを$Le \neq 1$で再計算
- 可能であればDNSとの比較で検証
サーモアコースティック不安定性
燃焼振動のCFD予測はどこまでできますか？
LESで圧力変動と発熱変動の時系列を取得し、Rayleigh indexを計算する。
$RI > 0$ の領域が不安定性を駆動する。この手法でガスタービンの燃焼振動周波数と振幅を定量的に予測できるようになってきた。CERFACS（仏）やTU Munich の研究グループが先導している。
DNS-informed モデル
DNSの成果がモデル改良に使われているんですか？
PeTA-FLOPSクラスのHPCにより、予混合火炎のDNS（$Ka$ = 100以上のthin reaction zone条件）が可能になった。DNSデータから以下の知見が得られている。
| DNS知見 | モデルへの影響 |
|---|---|
| 高Ka数で火炎面が断裂する | FSD方程式の消滅項修正 |
| 差拡散で局所$S_L$が変動 | TFM効率関数の修正 |
| 火炎-乱流の逆相関 | $\Xi$-equation の生成項修正 |
DNSがモデル改良のための「数値実験」として使われているんですね。
そうだ。実験では計測困難な火炎面の微細構造をDNSが明らかにし、RANSやLESのサブモデルを改良する。この「DNS → モデル改良 → 3D検証」のサイクルが現代の燃焼CFD研究の標準的な方法論だ。
水素の燃焼速度340m/sがエンジン設計者を悩ませる理由
ガソリンの層流燃焼速度が約0.4m/sなのに対し、水素は約3m/s、希薄条件でも燃え広がる。そして一部の研究では爆ごう近傍で340m/sに達することが報告されている。この異常な速さのせいで、水素を予混合エンジンに使うと「早期着火（プレイグニッション）」や「逆火（バックファイア）」が多発する。吸気ポートに火炎が逆流してくる現象だ。現在の水素エンジン開発では、これを防ぐために直噴＋高EGR戦略が主流になっているが、燃焼速度が速すぎる水素の扱いは今も先端的な研究課題であり続けている。
予混合火炎のトラブル対応
トラブルシューティング
予混合火炎計算のトラブルを教えてください。
1. 火炎が上流に伝播（フラッシュバック）
症状: 火炎が燃焼室からミキシングゾーンに逆流する。
原因と対策:
- 乱流燃焼速度 $S_T$ が過大に評価されている → Zimont定数$A$を下げる（0.52 → 0.4で試行）
- メッシュが粗すぎて数値拡散で火炎フロントが上流に広がる → 火炎面を含む領域を細分化
- 境界層内の低速領域を火炎が伝播している → 壁面近傍のメッシュ解像度を確認
2. 火炎が安定しない（ブローオフ）
症状: 火炎が燃焼室から流出して消滅する。
対策:
- $S_L$ の設定値が低すぎないか確認（特に高圧条件では$S_L$が減少する）
- パイロット火炎や保炎器の再循環域が十分解像されているか
- 乱流モデルが再循環域を正しく予測しているか（k-$\varepsilon$では過度に拡散する場合あり、SST推奨）
3. 進行変数 $c$ の数値問題
| 症状 | 原因 | 対策 |
|---|---|---|
| $c$ が [0,1] 外に | 数値スキームのオーバーシュート | Bounded schemes使用（QUICK避ける） |
| $c$ のフロントが厚い | メッシュ解像度不足 | TFMの$F$値見直し、メッシュ細分化 |
| $c$ がノイジー | 圧力-密度の結合不良 | 密度比補正、Under-Relaxation調整 |
4. 圧力振動（サーモアコースティック）
計算で非物理的な圧力振動が出る場合は？
5. 当量比変動がある場合
燃料濃度が均一でない場合はどうですか？
完全予混合ではなく部分予混合（Partially Premixed）条件だ。FluentのPartially Premixed Combustion（フレームレット+進行変数）やSTAR-CCM+のFGMで対応する。純粋な予混合モデルを適用すると、局所的に過濃・希薄な領域で非物理的な結果になる。
予混合火炎は「$S_L$の正確さ」と「火炎面の数値解像度」がトラブルの根本原因なんですね。
そのとおり。まずCanteraで$S_L$を検証し、次にメッシュ解像度を確認する。この2ステップを省略するとデバッグが泥沼化するから、必ず実行しよう。
予混合炎が「振動する」——サーモアコースティック不安定の怖い現実
予混合炎のトラブルシューティングで見落とされがちなのが「サーモアコースティック不安定」だ。燃焼室内で熱放出率が振動すると、圧力波と共鳴して振動が増幅される。低NOxバーナーで希薄予混合にすると安定燃焼の余裕が減り、この振動が起きやすくなる。実際、1990年代のガスタービン新型燃焼器でサーモアコースティック振動による燃焼器壁面の損傷事故が複数報告されている。CFDで予測するには非定常LESが必要で、定常RANSでは見えない問題だ。「定常計算で問題なし」→「実機でバーナー振動で部品破損」というパターンは、設計レビューで特に警戒すべき罠だ。
関連トピック
なった
詳しく
報告