すすモデル
すすの理論基礎
概要
先生、すす（soot）はなぜ燃焼CFDで重要なんですか？
すすは不完全燃焼で生成される炭素微粒子（粒径10-100 nm）で、3つの理由で重要だ。(1) 排出ガス規制の対象（PM：Particulate Matter）、(2) 輻射伝熱への寄与が大きい（火炎の輻射はすすが支配的）、(3) 健康被害（発がん性）。ディーゼルエンジン、航空機エンジン、工業炉でのすす予測は必須の課題だ。
すすの生成メカニズムを教えてください。
すすの生成は4段階のプロセスで進行する。
1. 核生成（Nucleation）: PAH（多環芳香族炭化水素）の重合により最初のすす核が形成される。C2H2（アセチレン）がPAH成長の主要前駆体
2. 表面成長（Surface Growth）: HACA機構（H-Abstraction-C2H2-Addition）によりすす粒子表面に炭素が堆積
3. 凝集（Coagulation）: 粒子同士が衝突・合体して大きくなる
4. 酸化（Oxidation）: O2やOHによりすすが燃焼・消滅する
Moss-Brookeモデル
CFDで使われるすすモデルの支配方程式を教えてください。
Moss-Brooke 2変数モデルは、すす質量分率 $Y_s$ とすす数密度 $N$（粒子数/kg）の2つの輸送方程式を解く。
すす体積分率 $f_v$ は次式で求まる。
ここで $\rho_s \approx 1800$ kg/m$^3$ はすすの密度、$d_p$ は平均粒径だ。
すす生成の主要パラメータ
どういう条件ですすが多く出るんですか？
すす生成の主要因子をまとめよう。
| 因子 | すす増加方向 | 理由 |
|---|---|---|
| 当量比 | 過濃（$\phi > 1$） | 酸素不足で不完全燃焼 |
| 温度 | 1500-1800 K | 核生成の最適温度域 |
| 圧力 | 高圧 | 衝突頻度増加 |
| 燃料構造 | 芳香族 > 直鎖 | PAH前駆体生成しやすい |
| 滞留時間 | 長い | すす成長の時間確保 |
1500-1800 Kがすす生成の温度窓なんですね。
そうだ。これより低温では核生成速度が遅く、高温ではOH酸化が優勢になってすすが燃え尽きる。この「soot formation window」は$\phi$-T マップとして可視化される。ディーゼル燃焼の $\phi$-T マップ（Dec diagram）はすすとNOxの同時低減戦略の基礎になっている。
すすモデルは化学反応速度論と粒子力学の融合なんですね。
そうだ。気相のPAH化学とすす粒子のダイナミクスの両方を正確に記述する必要があるから、燃焼モデリングの中でも最も難易度が高い分野の一つだ。
すすはナノサイズの炭素粒子——直径1nm以下の物体の生成を方程式で記述する挑戦
すすの生成は直径1〜100nmの炭素粒子が、燃焼中わずか数ミリ秒で形成・成長・凝集するプロセスだ。これをCFDで表現するためには「核生成→表面成長→凝集→酸化」という4つのプロセスを数式化する必要がある。Fenimore-Jonesの2方程式モデルや、Frenklachのdetailed soot modelなど流派があるが、特に難しいのが粒子サイズ分布（PSD）の記述だ。粒径1nmと100nmでは挙動が全然違うのに、CFDの格子はミリメートルオーダー——つまり物体のサイズが計算格子の100万分の1という異次元の問題を扱っている。
すすの数値計算手法
数値手法の詳細
すすモデルの数値実装について教えてください。
CFDで使われるすすモデルは大きく3つに分類される。
| モデル | 精度 | 計算コスト | 特徴 |
|---|---|---|---|
| 経験的2変数モデル | 低-中 | 低 | Moss-Brooke、$Y_s$と$N$を輸送 |
| Method of Moments (MoM) | 中-高 | 中 | 粒径分布のモーメントを輸送 |
| Sectional法 | 高 | 高 | 粒径分布を離散セクションで解像 |
経験的2変数モデル
最もシンプルなモデルから教えてください。
Moss-Brooke（Fluent標準搭載）は$Y_s$（すす質量分率）と$N$（数密度）の2変数で粒子集団を記述する。核生成、表面成長、凝集、酸化の各プロセスにArrhenius型の速度式を使う。簡便だが粒径分布の情報は平均値のみになる。
Method of Moments (MOMIC)
Moment法とは何ですか？
粒径分布関数 $n(v,t)$（$v$は粒子体積）のモーメント $M_r = \int_0^\infty v^r n(v) dv$ の輸送方程式を解く手法だ。$M_0$ が数密度、$M_1$ が体積分率に対応する。Frenklach & Harrisが提唱したMOMIC（Method of Moments with Interpolative Closure）がFluent 2020以降で利用可能だ。
Sectional法
Sectional法の利点は？
粒径範囲を離散的なセクション（bin）に分割し、各セクションの数密度を個別に輸送する。粒径分布の形状が任意に表現できるため最も高精度だが、20-30セクションの追加スカラー輸送が必要で計算コストが高い。STAR-CCM+やCONVERGEで利用可能だ。
Fluentでの設定
Fluentでのすすモデル設定手順を教えてください。
1. Models > Species > Species Transport（燃焼モデル設定済み前提）
2. Models > Soot > Moss-Brooke（簡易）またはMOMIC（推奨）
3. PAH前駆体化学種の設定（C2H2, C6H6等）-- 反応機構に含まれている必要がある
4. 輻射モデルとの連成 -- すすの吸収係数を輻射モデルに渡す
重要な注意点として、すすモデルにはPAH前駆体（最低でもC2H2）を含む反応機構が必要だ。グローバル1段機構ではC2H2が含まれないため、すす計算はできない。DRM-19以上の機構を使おう。
すすと輻射の連成
すすと輻射はどう連成するんですか？
すす粒子は連続スペクトルの輻射を放射・吸収する。すすの吸収係数は次式で近似される。
ここで $C_0$ と $C_2$ は光学定数だ。ガス輻射（CO2, H2Oの帯状輻射）にすすの連続輻射が加わるため、すすが多い火炎では輻射損失が大幅に増加する。
すすモデルは燃焼モデル＋粒子モデル＋輻射モデルの三重連成なんですね。
そうだ。モデルの複雑さと計算コストのバランスが重要で、まずMoss-Brooke 2変数モデルで傾向を掴み、必要に応じてMOMICやSectional法に進むのが実務的だ。
すすモデルの数値手法——輸送方程式法とセクショナル法の選択
CFDにおけるすす解析には主に「2方程式モデル（Moss-Brookes：すす数密度N＋すす体積分率f）」と「セクショナル法（粒子径分布を離散ビンで追跡）」がある。2方程式モデルは計算コストが低くFluentに標準実装されているが、粒子径分布の詳細が得られない。セクショナル法は粒子の核形成・凝集・表面成長を粒径別に解像でき、光散乱・吸収特性の予測精度が高い。最近はMoM（Method of Moments）が両者の中間として注目されており、計算コストとモデリング精度のバランスが優れている。航空エンジンのNVPM（非揮発性粒子状物質）規制対応でセクショナル法の実用化が加速している。
すすの実務適用
実践ガイド
すす解析の実務手順を教えてください。
ディーゼルエンジンやガスタービンでのすす予測のフローだ。
解析フロー
1. 燃焼モデルの確立 -- まず温度場と化学種分布（特にC2H2, OH）を正確に計算する
2. すすモデルの追加 -- 燃焼モデルが安定した後にすすモデルをオンにする
3. 輻射モデルとの連成 -- すすの吸収係数を輻射に渡す（DO + WSGGM + すす吸収）
4. 後処理 -- すす体積分率 $f_v$, 粒径分布, 輻射熱流束を評価
ディーゼル噴霧火炎でのすす解析
ディーゼルエンジンのすす解析で特有の注意点はありますか？
ディーゼルのすす解析ではスプレーモデルとの連成が鍵だ。
| 要素 | 推奨設定 | 備考 |
|---|---|---|
| 噴霧モデル | KH-RT（Kelvin-Helmholtz / Rayleigh-Taylor） | 液滴分裂の標準モデル |
| 蒸発モデル | Frossling | マルチコンポーネント蒸発 |
| 燃焼モデル | SAGE (CONVERGE) or EDC (Fluent) | 詳細化学反応必須 |
| すすモデル | MOMIC or Sectional | 粒径分布が重要 |
| 反応機構 | n-heptane skeletal (44種/112反応) | ディーゼルサロゲート |
n-ヘプタンがディーゼル燃料の代理なんですね。
そうだ。実際のディーゼル燃料は数百成分の混合物だが、着火特性とすす生成傾向を再現するためにn-ヘプタンやn-ドデカンをサロゲート燃料として使う。ECN（Engine Combustion Network）のSpray Aベンチマークがディーゼル噴霧火炎の標準検証ケースだ。
$\phi$-T マップによるすす/NOx分析
$\phi$-T マップの使い方を教えてください。
CFD結果の各セルを当量比 $\phi$ と温度 $T$ の散布図にプロットする。すす生成領域（$\phi > 2$, $T = 1500-1800$ K）とNOx生成領域（$\phi \approx 1$, $T > 2200$ K）を重ね合わせると、soot-NOx tradeoffが可視化できる。
よくある失敗と対策
| 症状 | 原因 | 対策 |
|---|---|---|
| すすが全く生成されない | 反応機構にC2H2が含まれていない | DRM-19以上の機構に変更 |
| すすが過大 | 酸化モデルが弱い | OHによる酸化パスの確認 |
| すすの空間分布が合わない | 混合場の解像度不足 | 噴霧近傍メッシュの細分化 |
| 輻射損失が大きすぎる | すす吸収係数の過大評価 | すす光学定数の確認 |
すす解析は燃焼・スプレー・輻射の全てが絡むので、段階的に構築するのが鉄則ですね。
そのとおり。まず冷間スプレー → 着火・燃焼 → すす → 輻射の順に積み上げていく。一度に全部有効にすると収束しなくなることが多い。
ディーゼルエンジンの「NOx-PMトレードオフ」——すすモデルが規制対応の実務を支える理由
ディーゼルエンジンの排気規制対応でエンジニアが直面する「NOx-PMトレードオフ」は、燃焼温度を上げると窒素酸化物（NOx）が増えるが、温度を下げて酸素濃度を下げると今度はすす（PM）が増える、という根本的な矛盾だ。Euro 6・国内重量車規制への対応には、EGR（排気再循環）量・噴射タイミング・コモンレール圧力の3変数を同時に最適化する必要がある。すすモデル付きCFD（Fenimore-JonesモデルやTwo-Equation Sootモデル）は、この多変数最適化の試験回数を台上実験の1/10以下に削減できる。「まずCFDで方向性を絞り、最後だけエンジン台上で確認」という開発プロセスは、現在の自動車メーカーで標準的になっている。
すすのすすのソフトウェア比較
商用ツール比較
| ツール | 経験的モデル | MOMIC | Sectional | PAH詳細 | 最適用途 |
|---|---|---|---|---|---|
| Ansys Fluent | Moss-Brooke | 2020R1以降 | なし | CHEMKIN連携 | 汎用バーナー |
| STAR-CCM+ | 2変数 | 搭載 | 搭載 | DARS連携 | ガスタービン |
| CONVERGE | Hiroyasu | 搭載 | 搭載 | SAGE | ディーゼルエンジン |
| OpenFOAM | コミュニティ | コミュニティ | 研究実装 | Cantera連携 | 研究用途 |
ディーゼルエンジンのすす解析にはやはりCONVERGEですか？
CONVERGEはスプレー+燃焼+すすの全連成が最も成熟している。AMRが自動で噴霧・火炎・すす生成領域を追跡するため、手動メッシュ調整が不要な点が大きな利点だ。
Fluent固有の機能
FluentではMoss-Brookeモデルが長年の標準だったが、2020R1でMOMICが追加された。PAH前駆体としてC6H6（ベンゼン）やA4（ピレン）まで扱える。すす-輻射連成はDO/P-1モデルとの組み合わせで自動設定される。
STAR-CCM+固有の機能
STAR-CCM+はSectional法が標準搭載されている点が強みだ。20セクションで粒径1-500 nmの分布を解像でき、DPF（Diesel Particulate Filter）との連成解析にも対応する。
選定の指針
結局どれを選べばいいですか？
すすモデルの選択は、求める出力（総量 vs 粒径分布）で決まるんですね。
そうだ。排出規制値との比較だけなら2変数モデルで十分だが、粒子数規制（PN）や粒径分布評価が必要ならSectional法が必須だ。
すすCFDツール——Ansys Fluent vs OpenFOAM sootFoam
すすモデルのCFDツールはAnsys Fluent（Moss-Brookes/Enhanced Moss-Brookes）とOpenFOAMのsootFoam（セクショナル法）が代表的だ。Fluentは産業用バーナ・ガスタービン設計での実績が豊富で、詳細化学反応とすすモデルのカップリングが整備されている。OpenFOAMはsootFoamで粒子径分布追跡が可能で、航空機排気粒子（NVPM）の研究に多く使われる。Numecaは燃焼専用ツールFINE/Openで高精度の乱流燃焼-すす連成解析を提供している。どのツールも汎用の反応機構ファイル（CHEMKIN format）を取り込むことができ、燃料変更への対応が容易だ。
すすの先端研究
先端トピックと研究動向
すすモデルの最先端研究を教えてください。
すすモデリングは燃焼CFDの中で最も活発な研究分野の一つだ。3つの方向性がある。
PAH化学の詳細化
PAH（多環芳香族炭化水素）の化学がどう進展していますか？
すす核生成のメカニズムとして、従来のHACA（H-Abstraction-C2H2-Addition）に加えて、PAH同士のファンデルワールス力による結合（Physical Nucleation）やラジカル-ラジカル結合（CHRCR機構）が提案されている。
| 核生成機構 | 研究グループ | 特徴 |
|---|---|---|
| HACA (Chemical Growth) | Frenklach (Berkeley) | 古典的、高温で有効 |
| Physical Nucleation | Kraft (Cambridge) | ファンデルワールス結合 |
| CHRCR | Wang (Stanford) | ラジカル-ラジカル結合 |
| Resonance-stabilized | Mebel (FIU) | 共鳴安定化ラジカル経路 |
Hybrid Method of Moments
Moment法の発展はどうですか？
HMOM（Hybrid Method of Moments）はMueller & Pitsch（2012）が提案した手法で、MOMICの粒径分布仮定の限界を克服する。二峰性（bimodal）の粒径分布を扱えるため、核生成と表面成長が同時に進行する領域での精度が向上する。
LESとすすモデルの連成
LESでのすす予測は可能ですか？
LESのサブグリッドスケールですす生成速度が非線形に変動するため、フィルタード反応速度の取り扱いが課題だ。Conditional Moment Closure（CMC）やTransported PDFとSectional法を組み合わせたLES-CMC-Sectionalが最高精度の手法として研究されている。計算コストは極めて大きいが、Sandia flame Dや航空機エンジンモデル燃焼器での検証が進んでいる。
機械学習の適用
すす研究は分子レベルの化学と粒子物理学の最前線なんですね。
そうだ。すすの核生成メカニズムは60年以上議論されてきたがまだ決着していない。CFDモデルの精度向上には、基礎化学の進歩が不可欠だ。
すすモデルの最前線——PAH（多環芳香族炭化水素）成長の直接数値シミュレーション
すす生成の核心は「PAH（多環芳香族炭化水素）の成長・凝集・表面反応」だが、これを完全に再現する計算コストは膨大だ。2010年代からDNS（直接数値シミュレーション）でPAH形成の素反応が追跡され、アセチレン（C₂H₂）からベンゼン（C₆H₆）、さらにピレン（C₁₆H₁₀）を経る成長経路の反応速度定数が定量化された。現在のLES/RANSレベルのすすモデル（Moss-Brookes、Semi-empirical, HACA）はDNSから得られたパラメータで検証されている。2030年代の航空機エンジン設計に向け、非化石燃料（SAF）でのすす生成特性の精密予測がCAE開発の最前線課題だ。
すすのトラブル対応
トラブルシューティング
すす計算のトラブルを教えてください。
1. すすが全く生成されない
原因チェックリスト:
- 反応機構にC2H2が含まれているか確認。グローバル1段機構ではすす前駆体が存在しない
- すすモデルが正しく有効化されているか
- 火炎温度がすす生成窓（1500-1800 K）に入っているか
- 当量比が局所的に過濃（$\phi > 1$）になる領域があるか
完全予混合リーンバーン条件ではすすが出ないのは正常ですか？
正常だ。$\phi < 0.8$ の均質リーンバーン条件ではすす生成はほぼゼロになる。すすが問題になるのは過濃条件（ディーゼル噴霧、リッチバーン域）だ。
2. すすが実験の10-100倍過大
原因と対策:
- Moss-Brookeモデルの核生成定数が過大。デフォルト値はエチレン火炎に合わせてあり、メタン火炎では過大評価する場合がある
- OH酸化が弱い。OHによるすす酸化はFenimoreのモデルで記述されるが、OH質量分率が不正確だと酸化が不十分になる
- 乱流-すす相互作用のモデル化不足。RANSでは平均場からすす生成を計算するが、局所的な温度・濃度変動を考慮していない
3. すすと輻射の連成問題
| 症状 | 原因 | 対策 |
|---|---|---|
| 火炎温度が急降下 | すすの輻射が過大 | すす光学定数を確認、$f_v$が過大でないか |
| 輻射フィードバックで振動 | すす→輻射→温度→すすの非線形ループ | Under-Relaxation追加、段階的連成 |
| 壁面熱流束が合わない | すす輻射の壁面吸収モデル不足 | DO模型の角度分割を増やす |
4. 計算が遅い（Sectional法）
Sectional法で20-30セクション追加すると化学種数が大幅に増えてメモリと計算時間が増大する。
対策:
- セクション数を減らす（10セクションでも傾向は掴める）
- すすモデルのUpdate Frequencyを設定（毎反復ではなく5-10反復に1回）
- AMRですす生成領域のみ解像度を上げ、それ以外を粗くする（CONVERGE）
デバッグ手順
1. まず燃焼モデルのみで温度・C2H2・OH分布を確認
2. すすモデルを追加してすす体積分率$f_v$の分布を確認
3. 輻射モデルを追加してすす-輻射連成の影響を確認
4. 各段階で実験データと比較
すす計算は段階的構築と実験データとの照合が不可欠ですね。
そうだ。すすモデルのパラメータは燃料や条件に依存する部分が大きいから、対象条件に近い実験データでの検証なしに結果を信頼してはいけない。
「すす予測が実測と合わない」——反応機構と乱流モデルの相互作用
CFDすす計算と光学計測（LII法・消光法）の結果が1桁以上乖離するケースは珍しくない。原因の多くは①乱流-化学反応相互作用モデルの不正確さ、②すす前駆体（PAH）の不確かな反応速度定数、③放射冷却のすす形成への影響の三つだ。特に乱流フラクチュエーションがすすの局所形成速度を変調する「乱流すす相互作用（TSI）」の精確なモデル化が研究課題だ。実測との差異を「成分（核形成/凝集/酸化）」に分解するための感度解析とCFD-実験連携アプローチが、産業燃焼器設計での精度向上の標準手法になりつつある。
関連トピック
なった
詳しく
報告