化学反応速度論の基礎 — トラブルシューティングガイド

燃焼CFDの化学反応関連のトラブルは、大きく分けて発散問題、精度問題、性能問題の3カテゴリーに分かれる。

症状: 温度や化学種質量分率が非物理的な値（負の質量分率、10000 K超の温度）になり、ソルバーがクラッシュする。

対策チェックリスト:

• Stiff ODEソルバーを使っているか: FluentならStiff Chemistry Solver、OpenFOAMならodeソルバーを確認
• 時間刻みは適切か: 化学反応のタイムスケールに対してCFL条件を満たしているか。CFD側のΔtが化学反応の最小時間スケールを大きく超えるとoperator splitting誤差が増大する
• 初期条件に未燃混合気の組成が正しく設定されているか: $\sum Y_i = 1$ の制約を確認
• Species Boundingを有効にする: 質量分率を [0, 1] にクリップする（Fluentでデフォルト有効）

数値拡散と化学反応ソース項の競合が原因だ。特にマイナー化学種（ラジカル）で起きやすい。対策としてSpecies Boundingに加え、2次精度の差分スキーム（Fluent: Second Order Upwind、OpenFOAM: limitedLinear）を使うことが重要だ。

考えられる原因と対策:

• パッチ着火が不十分: 初期温度パッチの範囲と温度レベルが低い。着火カーネルとして2000-2500 Kの領域を数セル分設定する
• 数値拡散で火炎が消える: メッシュが粗すぎて火炎フロントが解像されていない。層流火炎厚 $\delta_L$ あたり最低5-10セル必要
• 反応機構に問題: グローバル1段機構では低温着火を再現できない。0Dで着火遅れを検証
• 乱流-化学反応相互作用: Finite Rate Onlyでは乱流の影響を無視している。EDCやPaSRモデルの使用を検討

原因と対策を表にまとめよう。

• ISATを有効にする: 既に有効なら許容誤差を緩める（$10^{-4} \to 10^{-3}$）
• 反応機構を縮約する: DRM-19やLu13など縮約版に切り替え
• 化学反応計算をGPUにオフロード: Fluent 2024R1以降でGPU Chemistry Solver対応
• AMR（Adaptive Mesh Refinement）: 火炎帯のみメッシュを細分化（CONVERGEの得意分野）

• "Temperature limited to xxx K in xx cells": 温度がユーザー設定の上下限を超えたセルがある。初期条件・境界条件を確認
• "Negative species in xx cells": 質量分率が負になった。メッシュ品質とスキームを確認
• "ISAT table full": ISATテーブルがメモリ上限に到達。テーブルサイズを増やすかmax-storage設定を調整

そうだ。3Dで問題が起きたときも、まず該当条件で0D計算を回して反応機構自体の問題かCFD側の問題かを切り分ける。これが最も効率的なデバッグ手順だ。