すすモデル
理論と物理
概要
先生、すす（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という異次元の問題を扱っている。
各項の物理的意味
- 時間項 $\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で使われるすすモデルは大きく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（非揮発性粒子状物質）規制対応でセクショナル法の実用化が加速している。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
実践ガイド
すす解析の実務手順を教えてください。
ディーゼルエンジンやガスタービンでのすす予測のフローだ。
解析フロー
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で方向性を絞り、最後だけエンジン台上で確認」という開発プロセスは、現在の自動車メーカーで標準的になっている。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ 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）を取り込むことができ、燃料変更への対応が容易だ。
選定で最も重要な3つの問い
- 「何を解くか」：すすモデルに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端トピックと研究動向
すすモデルの最先端研究を教えてください。
すすモデリングは燃焼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-実験連携アプローチが、産業燃焼器設計での精度向上の標準手法になりつつある。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——すすモデルの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告