Numerical Modelling of High-Temperature Methane Pyrolysis Under Plasma-Inspired Thermal Gradients

技術要約：プラズマ由来の熱勾配下における高温メタン熱分解の数値モデリング

問題提起
本論文は、遠隔地の天然ガス資源（主にメタン）を輸送可能な液体炭化水素へと変換するという課題に取り組んでいる。マイクロ波プラザ技術は、従来のエネルギー集約的な間接経路（蒸気メタン改質とその後のフィッシャー・トロプシュ合成など）に代わる直接的な選択肢を提供するものであるが、重大な運用上の障壁に直面している。具体的には、高温プラズマ領域における固体すす粒子の同時生成が、リアクターの閉塞や目標とする液体原料への選択性の低下を引き起こす。さらに、既存の数値シミュレーションは、プラズマリアクターに見られる複雑な熱勾配を過度に簡略化し、均一な熱ゾーンや基本的な一次元近似を仮定する傾向がある。これらの簡略化は、分子の真の化学的軌跡を捉えることができず、生成物の収率やすす生成率の不正確な予測を招く。核心となる問題は、高温度の定常状態プラズマコアと、リアクター出口における急速な熱クエンチ（消熱）ゾーンを厳密に切り離しつつ、多次元流体力学結合に対して計算効率を維持できる、厳密で予測可能なフレームワークが欠如していることである。

手法
本研究では、熱分解キネティック・フレームワークを用い、電子キネティクスを明示的に排除した中性ラジカルネットワークに基づき、大気圧リアクター環境におけるメタン変換をモデリングしている。

• 反応機構の拡張: 基盤となるメカニズムは、メタン熱分解に関するディーン（Dean）メカニズムに基づいており、マイクロ波プラズマアークで典型的な300～4000 Kの温度範囲をカバーするように拡張されている。極めて重要な追加要素として、アルゴン・バスガスとの衝突によって駆動される、電子励起状態の一重項メチレン（$^1CH_2$）から三重項メチレン（$^3CH_2$）への衝突誘起系間交差（ISC）が含まれている。この量子力学的遷移は、特に急速な熱クエンチ時におけるラジカルプールの発達に大きな影響を与える。
• 熱プロファイルのデカップリング: モデルはリアクターを以下の2つの異なるゾーンに分離している。
  1. 活性プラズマコア: 電磁エネルギー吸収によって、高度に励起されたラジカル種が準定常状態となる2500 Kから4000 Kの間で維持される、滞留時間約15.54 ms（表3では155.4 sと記載されているが、原文のフォーマットエラーと思われる）の連続加熱ゾーン。
  2. 冷却ジェット領域: ガスが約76 ms（表3に記載された76.9 ms）かけて2500 Kから300 Kへと急速に冷却される、拡散限定のクエンチゾーン。このプロファイルは、レイリー散乱およびOESデータによって検証された、1次元径方向熱伝導方程式を解くことで導出されている。
• 電子キネティクスの無視: 0次元ボルツマン解析により、1 barおよび2500 Kにおいて、熱解離速度が電子衝突解離速度を7桁上回ることが確認された。したがって、本モデルは純粋な熱分解プロセスとして扱い、中性ラジカル化学によって制御されるものとする。
• メカニズムの簡約化: 多次元計算流体力学（CFD）への統合を可能にするため、KINALCソフトウェアを用いて2段階の簡約手順を実行した。これには、冗長な反応の骨格的除去と、反応速度における構造的多重共線性に対処するための主成分分析（PCA）が含まれる。

主な結果

• 転換率と収率: モデルは98.7%の純メタン転換率を予測している。理想的な単相（ガスのみ）の仮定下では、主要生成物は水素とアセチレン（収率88.0%）であり、少量のエチレン（7.0%）、プロペン（1.1%）、ブタジエン（1.0%）が含まれる。
• 逐次的なキネティック・フリーズ（速度論的凍結）: 本研究は、クエンチ中の生成物形成を支配する明確な熱的閾値を特定している。
  • アセチレンは3900 K以下で安定化し、凍結する。
  • エチレンは1520 K以下で著しく蓄積する。
  • プロペンおよびブタジエンは、それぞれ1060 Kおよび870 Kで出現する。
  • キネティック・ロック: すべての分子変換は500 K以下で停止し、最終的な生成物分布を固定する。
• すす沈着の影響: 固相炭素析出を含む多段階シミュレーションによれば、気相モデルはアセチレン収率を著しく過大評価することが明らかになった。2500 Kで固相炭素の析出を許容すると、炭素質量の99%がすすとして脱落する。この炭素が枯渇した水素リッチな環境では、その後のクエンチ化学により、最終生成物が不飽和種（アセチレン/エチレン）から安定なエタンおよび未反応メタンへとシフトする。これは単相モデルが予測する熱力学的極小値とは対照的な、多相モデル特有の結果である。
• 計算効率: メカニズムの簡約化により、詳細なスキーム（数百の種/ステップ）を、26種の成分と102個の不可逆反応からなる軽量なメカニズムへと圧縮することに成功した。抄録では61倍の加速が引用されているが、セクション3.5では実行プロファイルにおける14.2倍の高速化係数が明記されている。この簡約化されたモデルは、フルメカニズムの種濃度を0.1%以内の誤差で再現している。

意義および主張
本論文は、反応ガスの混合物が経験する複雑な温度－時間プロファイルのより厳密な推定値を確立することにより、数値プラズマリアクターの予測能力を向上させると主張している。高温度コアと急速なクエンチゾーンを正確にデカップリングすることで、詳細な化学キネティクス・モデリングのための物理的に現実的なフレームワークを提供している。

主な意義は、**キネティック・フリーズ（速度論的凍結）**が特定の熱的閾値によって支配される逐次的なプロセスであり、それが最終的な生成物構成を決定するということを示した点にある。さらに、本研究は、**不均一相変化（すす生成）**が、単相モデルによって予測される熱力学的極小値から生成物分布をシフトさせ、気相反応経路を根本的に変えることを強調している。

最後に、本研究は計算効率の高いフレームワークを提示している。フルキネティックスキームに対して検証されたこの簡約メカニズムは、詳細な化学キネティクスを伴う流体コードにおける膨大な計算コストという課題を克服し、将来的なエンジニアリング応用における高精度な3次元フルスケール工業用プラズマリアクターのシミュレーションを可能にするものである。