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

技术摘要：受等离子体启发的热梯度下高温甲烷热解的数值模拟

问题陈述
本文旨在解决将偏远天然气储量（主要为甲烷）转化为可运输液体烃类的挑战。虽然微波等离子体技术为传统的、高能耗的间接路径（如蒸汽甲烷重整结合费托合成）提供了一种直接的选择，但它面临着显著的运行障碍。具体而言，在高热等离子体区内固体烟炱颗粒的同步形成会导致反应器堵塞并降低目标液体原料的选择性。此外，现有的数值模拟往往过度简化了等离子体反应器中复杂的温度梯度，假设存在均匀的热区或基础的一维近似。这些简化处理无法捕捉分子的真实化学轨迹，导致对产物收率和烟眜形成速率的预测不准确。核心问题在于缺乏一个严谨的、具有预测能力的框架，能够准确地将高温稳态等离子体核心与反应器出口处的快速热淬火区进行解耦，同时保持足够的计算效率以进行多维流体动力学耦合。

方法论
本研究采用热解动力学框架来模拟大气压反应器环境下的甲烷转化，通过明确排除显式电子动力学，转而采用中性自由基网络。

• 动力学模型扩展： 基础机制基于用于甲烷热解的 Dean 机制，并将其扩展至微波等离子体弧典型的 300–4000 K 温度范围。一个关键的补充是纳入了由氩气稀释气体驱动的电子激发单线态亚甲基 ($^1CH_2$) 向三线态亚甲基 ($^3CH_2$) 的碰撞诱导系间窜越 (ISC) 过程。这种量子力学跃迁显著影响了自由基池的发展，特别是在快速热淬火期间。
• 热剖面解耦： 模型将反应器分为两个截然不同的区域：
  1. 活性等离子体核心： 一个连续加热区，停留时间约为 15.54 ms（表 3 中标注为 155.4 s，可能为原文格式错误），维持在 2500 K 至 4000 K 之间，其中电磁能吸收驱动系统达到高度激发的自由基物种的准稳态。
  2. 冷却喷射区： 一个扩散限制的淬火区，气体在此处经历从 2500 K 到 300 K 的快速冷却，历时约 76 ms（表 3 中明确为 76.9 ms）。该剖面通过求解一维径向热传导方程得出，并经过瑞利散射和 OES 数据的验证。
• 电子动力学的忽略： 零维玻尔兹曼分析证实，在 1 bar 和 2500 K 条件下，热解离速率比电子碰撞解离速率高出七个数量级。因此，模型将系统视为受中性自由基化学控制的纯热解过程。
• 机制简化： 为了能够集成到多维计算流体力学 (CFD) 中，使用 KINALC 软件执行了二步简化程序。这包括通过骨架消除冗余反应，随后利用主成分分析 (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% 以内。

意义与主张
本文声称通过建立对反应气体经历的复杂温度-时间剖面的更严谨估计，提升了数值等离子体反应器模拟的预测能力。通过准确解耦高温核心与快速淬火区，作者为详细的化学动力学建模提供了一个符合物理实际的框架。

其主要意义在于证明了动力学冻结是一个受特定热阈值支配的序列过程，它决定了最终的产物构型。此外，研究强调了异相变化（烟眜形成）从根本上改变了气相反应路径，使产物分布偏离了单相模型所预测的热力学极小值。

最后，这项工作提供了一个具有计算效率的框架，用于未来的工程应用。经过验证且符合全机制的简化机制，使得在高保真三维工业级全尺寸等离子体反应器中进行流体动力学模拟成为可能，克服了以往在流体动力学代码中使用如此详细的化学动力学时所面临的极高计算成本。