拡散火炎と混合分率
理論と物理
概要
先生、拡散火炎ってどういう燃焼形態ですか？
拡散火炎（非予混合火炎）は、燃料と酸化剤が別々に供給され、混合と同時に燃焼が進行する形態だ。ろうそくの炎、ディーゼルエンジン、ガスタービンの燃焼器がこれに該当する。火炎は燃料と酸化剤が化学量論比で出会う面（ストイキオメトリック面）に形成される。
予混合火炎とは本質的に違うんですね。
そうだ。予混合火炎では燃焼速度が支配パラメータだが、拡散火炎では混合速度が律速になる。この性質を利用して、混合分率 $Z$ という変数で火炎構造を記述するのがBurke-Schumann理論の核心だ。
混合分率
混合分率 $Z$ の定義を教えてください。
混合分率は「流体要素中にどれだけ燃料流由来の質量が含まれるか」を表す保存スカラーだ。Bilgerの定義が広く使われる。
ここで $Y_F$ は燃料質量分率、$Y_O$ は酸化剤質量分率、添字1は燃料流、添字2は酸化剤流の値を示す。$s$ は質量ベースの化学量論比だ。
$Z=0$ が純酸化剤、$Z=1$ が純燃料ということですか？
そのとおり。ストイキオメトリック混合分率 $Z_{st}$ は $Z_{st} = Y_{O,2}/(s\,Y_{F,1} + Y_{O,2})$ で、メタン/空気の場合 $Z_{st} \approx 0.055$ だ。火炎は $Z = Z_{st}$ の等値面に位置する。
Burke-Schumann解
Burke-Schumann解とは何ですか？
化学反応が無限に速い（$Da \to \infty$）と仮定すると、火炎面で燃料と酸化剤が瞬時に反応し、温度と化学種は混合分率のみの関数になる。
これがBurke-Schumann解で、温度は $Z_{st}$ でピークを取る三角形のプロファイルになる。この解はフレームレットモデルの出発点でもある。
混合分率の輸送方程式
混合分率 $Z$ はソース項なし（化学反応に依存しない）の輸送方程式に従う。
ソース項がないから、化学反応の詳細を知らなくても $Z$ の場を解けるわけですね。
そこが混合分率アプローチの大きな利点だ。化学反応の複雑さを $Z$ 空間のフレームレット方程式に押し込み、CFDでは $Z$ とその分散 $\widetilde{Z''^2}$ だけを輸送すればよい。
拡散火炎の理論は混合分率の概念が全ての基盤なんですね。
そうだ。混合分率は非予混合燃焼のCFDにおいて最も重要な変数と言っても過言ではない。
ろうそくの炎は「拡散炎」の教科書だった——Burke-Schumann解の限界
1928年にBurkeとSchumannが提唱した拡散炎の解析解は、燃料と酸化剤の混合が瞬時に完了するという「無限大反応速度」仮定に基づいている。ろうそくの炎はこのモデルで驚くほどよく再現できる。ところが、実際のガスタービン燃焼器で同じ仮定を使うと、NOx生成量が実測の10倍以上になることがある。混合が遅ければ「混合が律速」のはずなのに、高温域では反応速度も有限であることが露呈する。現場のエンジニアが「Burke-Schumann解は概念理解には最高だが、設計には使えない」と言う理由がここにある。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 速度 $u$ | m/s | 入口条件で体積流量から換算する際、断面積の単位に注意 |
| 圧力 $p$ | Pa | ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 |
| 密度 $\rho$ | kg/m³ | 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C |
| 粘性係数 $\mu$ | Pa·s | 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 |
| レイノルズ数 $Re$ | 無次元 | $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 |
| CFL数 | 無次元 | $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結 |
数値解法と実装
数値手法の詳細
混合分率ベースの燃焼モデルをCFDでどう実装するんですか？
非予混合燃焼モデルでは、$Z$ と $\widetilde{Z''^2}$（混合分率分散）の輸送方程式をCFDで解き、化学反応の情報はルックアップテーブルから取得する。
輸送方程式
乱流場でのFavre平均混合分率 $\widetilde{Z}$ と分散 $\widetilde{Z''^2}$ の方程式は次のようになる。
分散方程式の最後の項 $\widetilde{\chi}$ は何ですか？
スカラー散逸率（scalar dissipation rate）だ。$\chi = 2D|\nabla Z|^2$ で定義され、混合分率の微細構造がどのくらい速く消散するかを表す。乱流条件下では $\widetilde{\chi} = C_\chi \frac{\varepsilon}{k}\widetilde{Z''^2}$ のようにモデル化される（$C_\chi \approx 2.0$）。
PDF（確率密度関数）
PDFの役割を教えてください。
乱流場では $Z$ がセル内で揺らいでいるから、平均値だけでは火炎構造を正しく表せない。$Z$ の確率密度関数 $P(Z)$ を仮定して、平均温度や平均化学種を積分で求める。
PDFの形状には $\beta$ 関数分布が広く使われる。$\widetilde{Z}$ と $\widetilde{Z''^2}$ の2つのパラメータから $\beta$ 分布の形状が一意に決まるんだ。
ルックアップテーブルの構築
テーブルはどうやって作るんですか？
事前にフレームレット方程式やケミカルイクイリブリアムを解いて、$Z$ と $\chi_{st}$（ストイキオメトリック散逸率）をパラメータとした温度・化学種のテーブルを作る。これをPDFで積分したテーブルが、CFDのランタイムで参照される。
| テーブル変数 | 次元数 | 用途 |
|---|---|---|
| $\widetilde{Z}$, $\widetilde{Z''^2}$ | 2D | 平衡化学・薄い火炎 |
| $\widetilde{Z}$, $\widetilde{Z''^2}$, $\widetilde{\chi_{st}}$ | 3D | 定常フレームレット |
| $\widetilde{Z}$, $\widetilde{Z''^2}$, $\widetilde{C}$ | 3D | FGM (progress variable追加) |
Fluentではどう設定しますか？
Fluentでは Models > Species > Non-Premixed Combustion を選択し、CHEMKIN形式で反応機構をインポートしてPDFテーブルを自動生成させる。テーブル解像度（$Z$ 方向の分割数）は最低64点、できれば128点推奨だ。
混合分率 + PDFテーブルの手法は、化学反応を事前に解いてしまう点が巧みですね。
そうだ。3D CFDのランタイムでは化学反応を解く必要がないから、詳細化学反応でも計算コストがほぼ変わらないのが最大の強みだ。
スカラー散逸率χをどう解くか——拡散炎モデルの実装で最も議論になる変数
拡散炎のフレームレットモデルを実装するとき、最も実装者が頭を悩ませるのがスカラー散逸率χ（カイ）の計算だ。χは炎面での混合の激しさ（乱流による燃料-酸化剤の引き伸ばし率）を表し、χが高いほど炎が消えやすくなる。しかし乱流CFDでχを直接計算する式は複数あり、Favre平均・モデル定数の選択によって結果が変わる。特に再循環域（バーナー後流）でχが0に近い領域と、せん断層で高い領域が共存する場合、空間分布の計算精度がNOx・すす予測に直結する。「χをどう計算するか」がベンダー間の秘伝のタレになっていることも少なくない。
風上差分（Upwind）
1次風上: 数値拡散が大きいが安定。2次風上: 精度向上するが振動のリスク。高レイノルズ数流れでは必須。
中心差分（Central Differencing）
2次精度だが、Pe数 > 2で数値振動が発生。低レイノルズ数の拡散支配流れに適する。
TVDスキーム（MUSCL、QUICK等）
リミッタ関数により数値振動を抑制しつつ高精度を維持。衝撃波や急勾配の捕捉に有効。
有限体積法 vs 有限要素法
FVM: 保存則を自然に満足。CFDの主流。FEM: 複雑形状・マルチフィジックスに有利。SPH等のメッシュフリー法も発展中。
CFL条件（クーラン数）
陽解法: CFL ≤ 1が安定条件。陰解法: CFL > 1でも安定だが、精度と反復回数に影響。LES: CFL ≈ 1を推奨。物理的意味: 1タイムステップで情報が1セル以上進まないこと。
残差モニタリング
連続の式・運動量・エネルギーの各残差が3〜4桁低下で収束と判断。質量保存の残差は特に重要。
緩和係数
圧力: 0.2〜0.3、速度: 0.5〜0.7が一般的な初期値。発散する場合は緩和係数を下げる。収束後は上げて加速。
非定常計算の内部反復
各タイムステップ内で定常解に収束するまで反復。内部反復数: 5〜20回が目安。残差がタイムステップ間で変動する場合は時間刻みを見直す。
SIMPLE法のたとえ
SIMPLE法は「交互に調整する」手法。まず速度を仮に求め（予測ステップ）、その速度で質量保存が満たされるよう圧力を補正し（補正ステップ）、補正された圧力で速度を修正する——このキャッチボールを繰り返して正解に近づく。2人で棚を水平にする作業に似ている：片方が高さを合わせ、もう片方がバランスを取り、これを交互に繰り返す。
風上差分のたとえ
風上差分は「川の流れに立って上流の情報を重視する」手法。川の中にいる人が下流を見ても水の出所は分からない——上流の情報が下流を決めるという物理を反映した離散化手法。精度は1次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
実践ガイド
拡散火炎のCFD解析の実務手順を教えてください。
ガスタービン燃焼器やボイラーなどの非予混合燃焼解析の典型的なフローだ。
解析フロー
1. 燃料組成の定義 -- 天然ガス、灯油、水素など。多成分燃料の場合は代理混合物を定義
2. 反応機構の選定とPDFテーブル生成 -- Fluentの場合、Non-Premixed CombustionモデルでCHEMKIN形式をインポート
3. メッシュ生成 -- 混合領域（燃料噴射近傍）に十分な解像度を確保
4. 乱流モデル選定 -- k-$\varepsilon$系はスワール流で精度不足。SST k-$\omega$ または RSM を推奨
5. 境界条件設定 -- 燃料入口($Z=1$)と空気入口($Z=0$)。流量は全体当量比から逆算
6. 輻射モデル -- DO（Discrete Ordinates）またはP-1。WSGGMで吸収係数を設定
メッシュ設計の要点
混合分率モデルでメッシュの注意点は何ですか？
ミキシング精度がモデルの成否を決める。特に以下の領域を細分化する。
| 領域 | 推奨セルサイズ | 理由 |
|---|---|---|
| 燃料噴射孔出口 | 噴射孔径の1/10以下 | ジェット初期混合の解像 |
| 混合層（$0.01 < Z < 0.1$） | 1-2 mm | 火炎帯の解像 |
| スワール渦のコア | 2-3 mm | 再循環域の捕捉 |
| 壁面近傍 | y+ < 1 (LES) / 30-300 (RANS) | 壁面熱伝達 |
燃料噴射孔が1 mmだと、0.1 mmセルが必要になるんですね。
そうだ。だからガスタービン燃焼器の本格的なLESでは数千万セルが普通だ。RANSなら数百万セルで済むが、スワール誘起渦崩壊（VBB）の予測精度は限定的になる。
よくある失敗と対策
| 症状 | 原因 | 対策 |
|---|---|---|
| 火炎温度が断熱値を超える | PDFテーブルのZ方向分割が粗い | テーブル解像度を128点以上に |
| 排出口でZ>0の未燃燃料が残る | 混合不良（メッシュ粗い） | 噴射近傍メッシュを細分化 |
| NOxが実験の2倍以上 | thermal NOxの温度依存性が急峻 | 温度場精度を先に改善（輻射モデル） |
| 火炎長さが合わない | 乱流Sc数が不適切 | $Sc_t = 0.7$（デフォルト0.85を変更） |
乱流Schmidt数って変えていいんですか？
$Sc_t$ はモデル定数であり、燃焼流では0.7前後が実験と合う場合が多い。Fluentのデフォルト0.85はやや拡散が少ない。自由噴流火炎では0.7に下げると火炎長さの一致が改善されることが多い。
混合分率モデルは設定がシンプルだけど、ミキシングの解像度が成否を分けるんですね。
そうだ。「混合がすべてを支配する（Mixing controls everything）」は非予混合燃焼の鉄則だ。
ガスタービン燃焼器の「薄め燃焼」——拡散炎からの歴史的な転換
1960〜80年代のガスタービン燃焼器は、安定性重視の拡散炎設計が主流だった。燃料と空気を混ぜながら燃やせば確かに安定するが、高温域でNOxが激増する。排ガス規制の強化を受けて、1990年代以降は「予混合希薄燃焼（Lean Premixed, LPM）」が主流となった。GEのDLN（Dry Low NOx）バーナーや三菱重工のDLECバーナーはその代表例で、設計にはプリミックスドフレームレットモデルが欠かせない。ただし希薄予混合は吹き消えやすく、部分負荷での安定性確保が今もエンジニアを悩ませる実践的な課題だ。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
商用ツール比較
拡散火炎モデルのツール対応状況を教えてください。
非予混合燃焼（混合分率アプローチ）は主要CFDツールのほぼ全てで対応している。
| ツール | モデル名称 | PDFテーブル生成 | 多燃料流対応 | LES対応 |
|---|---|---|---|---|
| Ansys Fluent | Non-Premixed Combustion | 自動（CHEMKIN入力） | 2燃料流まで | あり |
| STAR-CCM+ | Flamelet/FGM | 自動（DARS連携） | 複数流対応 | あり |
| OpenFOAM | flameletFoam等 | 手動（Cantera等で事前生成） | カスタム | あり |
| CONVERGE | Flamelet model | SAGE連携 | 対応 | あり |
Ansys Fluent
Fluentの設定手順を教えてください。
1. Species > Non-Premixed Combustion を選択
2. 燃料組成（Rich Flammability Limit側）をCHEMKIN形式で定義
3. PDFテーブルを自動計算（Number of Grid Points: 128推奨）
4. 断熱/非断熱の選択（非断熱の場合、Enthalpy offsetテーブルが追加）
5. 境界条件で各Inletの$Z$値を設定（燃料入口: Z=1、空気入口: Z=0）
STAR-CCM+
STAR-CCM+ではどうですか？
STAR-CCM+ではFlamelets/FGMモデルとして実装されている。DARSライブラリと連携してフレームレットテーブルを自動生成する。Progress Variableを追加できるFGM拡張が強力で、局所消炎の再現が可能だ。
OpenFOAM
OpenFOAMでは標準搭載ですか？
標準ディストリビューションには基本的なフレームレットモデルがあるが、実務レベルではコミュニティ版の flameletFoam（TU Darmstadt開発）が広く使われている。PDFテーブルはCanteraで事前生成し、OpenFOAMのlookup tableとして読み込む。
選定の指針
どのツールがどういう場面に適していますか？
用途によって最適なツールが変わるんですね。
そうだ。非予混合燃焼モデル自体はどのツールでも同じ理論基盤だから、メッシュ生成の容易さ、LES対応、スプレーモデルとの連携など、周辺機能で差がつく。
Fluentの「Non-Premixed Combustion」——拡散炎モデルの実装がブラックボックスになる理由
ANSYSのFluentで拡散炎を計算する際に使う「Non-Premixed Combustion」モデルは、内部でPDF（確率密度関数）テーブルとフレームレット計算が走っている。ところがこのテーブル生成プロセスは、デフォルト設定のまま走らせると散逸率の範囲や格子点数がユーザーの意識外で決まってしまう。特に高圧条件（航空エンジン・高圧ボイラ）では実燃焼条件がテーブル範囲外になることがある。STAR-CCM+は同等の機能で設定の可視性が高い傾向があるが、その分パラメータが多くて初学者が混乱する。どちらも「デフォルトの中身を理解してから使う」という点は変わらない。
選定で最も重要な3つの問い
- 「何を解くか」：拡散火炎と混合分率に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端トピックと研究動向
拡散火炎の研究はどこまで進んでいますか？
拡散火炎のCFDモデリングは成熟した分野だが、まだ活発な研究が続いている。特に(1) 局所消炎と再着火、(2) 差圧噴霧火炎、(3) 水素拡散火炎の3分野が注目されている。
局所消炎と再着火
消炎現象は混合分率モデルでは扱えないんですか？
標準的なEquilibrium PDFモデルでは消炎を再現できない。スカラー散逸率が臨界値 $\chi_q$ を超えると火炎が局所的に消える現象を扱うには、Unsteady Flamelet Model（非定常フレームレット）やFGM（Flamelet Generated Manifold）にprogress variableを導入する必要がある。
消炎条件は次のDamkohler数で整理される。
$Da < Da_{\text{crit}} \approx 1$ で消炎が起きる。
Conditional Moment Closure (CMC)
CMCとは何ですか？
CMCは混合分率を条件変数として、化学種の条件付き平均を輸送する手法だ。Klimenkoと Bilgerが独立に提唱した。PDFテーブルの仮定に依存しないため、局所消炎やリフトオフ火炎の予測精度が高い。
| 手法 | 局所消炎 | 計算コスト | 実装難度 |
|---|---|---|---|
| Equilibrium PDF | 不可 | 低い | 容易 |
| Steady Flamelet | 限定的 | 低い | 容易 |
| FGM | 可能 | 中程度 | 中程度 |
| Unsteady Flamelet | 可能 | 中程度 | 中程度 |
| CMC | 良好 | 高い | 高い |
| Transported PDF | 最良 | 最高 | 高い |
水素拡散火炎
水素の拡散火炎で特別な考慮事項はありますか？
水素は質量拡散係数が他のガスの4-5倍大きい。このため等ルイス数仮定（$Le=1$）が成立せず、差拡散（differential diffusion）効果が顕著になる。通常の混合分率モデルは $Le=1$ を暗に仮定しているから、水素火炎では修正が必要だ。
具体的にはどう対処するんですか？
水素社会に向けて、水素拡散火炎のモデリングはますます重要になりそうですね。
そのとおり。水素のルイス数効果は火炎安定性やNOx排出にも影響するから、従来のメタン用モデルをそのまま転用するのは危険だ。
レーザー計測でCFDを壊す——LIF・PIVが暴いた拡散炎の「意外な真実」
プレーナーレーザー誘起蛍光（PLIF）やPIVが普及した2000年代、拡散炎のCFD検証は劇的に変わった。それまで「だいたい合ってる」と思われていたメタン-空気拡散炎のOHラジカル分布が、実測と比べると形状・ピーク値ともに20〜30%ずれることが明らかになった。原因はほとんどの場合「火炎伸長（strain rate）によるクエンチングの過小評価」だった。これが非定常フレームレット（Unsteady Flamelet）や消炎再着火モデルへの研究加速に直結した。実験計測がCFDの欠点を暴いて次世代モデルを生む——この循環が燃焼CFDの先端技術を動かしている。
トラブルシューティング
トラブルシューティング
混合分率モデルでよくあるトラブルを教えてください。
非予混合燃焼モデル特有のトラブルを整理しよう。
1. PDFテーブル生成エラー
症状: Fluent/STAR-CCM+でPDFテーブルの計算が失敗する。
原因と対策:
- 反応機構に問題: CHEMKIN形式の文法エラー。CanteraでYAML変換して検証
- 収束しないフレームレット: 極端に高い散逸率で解が発散。$\chi_{st}$ の上限を下げる
- 熱力学データの不整合: NASA7係数のフォーマットエラー。温度範囲を確認
2. 非物理的な温度
温度が断熱火炎温度を超えることがあるんですが…
原因: PDFテーブルの $Z$ 方向の離散化が粗い場合、$\beta$ PDF積分で補間誤差が生じる。特に $Z_{st}$ 近傍が重要だ。
対策:
- テーブル解像度を上げる（最低64点、推奨128点）
- $Z_{st}$ 近傍に格子点を集中させる（Fluentでは自動的に行われる）
- 非断熱モデルではEnthalpy Defect方向も十分な分割数を確保（最低20点）
3. 火炎形状が実験と合わない
| 症状 | 考えられる原因 | 対策 |
|---|---|---|
| 火炎が長すぎる | 乱流混合が弱い | $Sc_t$を下げる(0.7)、メッシュ細分化 |
| 火炎が短すぎる | 数値拡散による過度な混合 | 2次精度スキーム使用 |
| 火炎がリフトオフしない | Equilibriumモデルでは消炎不可 | FGMまたはPartially Premixedモデルに変更 |
| 火炎の半径方向広がりが大きい | 乱流モデルの拡散過大 | k-$\varepsilon$からSSTやRSMに変更 |
4. 多燃料流問題
燃料が複数箇所から入る場合はどうですか？
FluentのNon-Premixed Combustionモデルは基本的に1燃料流+1酸化剤流の2流体系だ。パイロット燃料と主燃料が異なる組成の場合、Secondary Streamを定義するか、Partially Premixedモデルに切り替える必要がある。STAR-CCM+のFGMモデルは複数燃料流にネイティブ対応している。
混合分率モデルのトラブルは、テーブル品質とミキシング解像度に集約されることが多いんですね。
そうだ。化学反応をテーブルに任せている以上、テーブルの品質とCFDでの$Z$場の精度がすべてだ。0D/1Dでテーブルの挙動を確認してから3Dに進む手順を守ろう。
「混合分率が1.0を超えた！」——拡散炎CFDの定番エラーの原因と対策
拡散炎のシミュレーションで初学者がよく遭遇するのが「混合分率Zが0〜1の範囲を外れる」問題だ。Zは定義上0（純酸化剤）〜1（純燃料）の範囲に収まるはずなのに、数値計算では境界条件の不整合や急激な勾配があると1.01や-0.02などの非物理的な値が出てしまう。こうなるとフレームレットテーブルの参照が範囲外になり、温度・密度が急変して計算が発散する。対処法は境界条件の再確認に加え、ScalarやZのクリッピング処理を有効にすること。Fluentでは「Bounded scalar」オプション、OpenFOAMでは`min/max`制約の設定がある。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——拡散火炎と混合分率の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告