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

Resumen Técnico: Modelado Numérico de la Pirólisis de Metano a Alta Temperatura bajo Gradientes Térmicos de Inspiración de Plasma

Planteamiento del Problema
El artículo aborda el desafío de convertir las reservas remotas de gas natural (principalmente metano) en hidrocarburos líquidos transportables. Si bien la tecnología de plasma de microondas ofrece una alternativa directa a las rutas indirectas tradicionales que consumen mucha energía (como la reformación de metano con vapor seguida de la síntesis de Fischer-Tropsch), esta enfrenta obstáculos operativos significativos. Específicamente, la formación simultánea de partículas de hollín sólido en las zonas de plasma de alta temperatura provoca la obstrucción del reactor y reduce la selectividad hacia los precursores líquidos objetivo. Además, las simulaciones numéricas existentes suelen simplificar excesivamente los complejos gradientes térmicos presentes en los reactores de plasma, asumiendo zonas térmicas uniformes o aproximaciones unidimensionales básicas. Estas simplificaciones no logran capturar las verdaderas trayectorias químicas de las moléculas, lo que conduce a predicciones inexactas de los rendimientos de los productos y de las tasas de formación de hollín. El problema central es la falta de un marco predictivo riguroso que desacople con precisión el núcleo de plasma de estado estacionario a alta temperatura de la zona de enfriamiento rápido (quench) en la salida del reactor, manteniendo al mismo tiempo la eficiencia computacional necesaria para el acoplamiento con la dinámica de fluidos multidimensional.

Metodología
El estudio emplea un marco cinético de pirólisis térmica para modelar la conversión de metano en un entorno de reactor a presión atmosférica, excluyendo explícitamente la cinética electrónica en favor de redes de radicales neutros bajo gradientes térmicos agresivos.

• Expansión del Modelo Cinético: El mecanismo fundacional se basa en el mecanismo de Dean para la pirólisis de metano, expandido para cubrir el espectro de temperatura de 300–4000 K típico de los arcos de plasma de microondas. Una adición crítica es la inclusión del cruce de sistemas intersistema (ISC) inducido por colisiones de la metileno singlete electrónicamente excitado ($^1CH_2$) al metileno triplete ($^3CH_2$), impulsado por colisiones con el gas de baño de argón. Esta transición de mecánica cuántica influye significativamente en el desarrollo del reservorio de radicales, particularmente durante el enfriamiento térmico rápido.
• Desacoplamiento del Perfil Térmico: El modelo separa el reactor en dos zonas distintas:
  1. Núcleo de Plasma Activo: Una zona de calentamiento continuo con un tiempo de residencia de aproximadamente 15.54 ms (notado como 155.4 s en la Tabla 3, probablemente un error de formato en el texto fuente) mantenida entre 2500 K y 4000 K, donde la absorción de energía electromagnética impulsa el sistema hacia un estado cuasi-estacionario de especies radicales altamente excitadas.
  2. Región de Jet de Enfriamiento: Una zona de enfriamiento rápido (quench) limitada por difusión donde el gas experimenta un enfriamiento rápido de 2500 K a 300 K durante aproximadamente 76 ms (específicamente 76.9 ms según la Tabla 3). Este perfil se deriva de la resolución de una ecuación de conducción de calor radial 1D, validada contra datos de dispersión de Rayleigh y OES.
• Inelegibilidad de la Cinética Electrónica: Un análisis de Boltzmann de cero dimensiones confirma que, a 1 bar y 2500 K, las tasas de disociación térmica exceden las tasas de disociación por impacto electrónico en siete órdenes de magnitud. En consecuencia, el modelo trata el sistema como un proceso de pirólisis puramente térmica gobernado por la química de radicales neutros.
• Reducción del Mecanismo: Para permitir la integración en la dinámica de fluidos computacional (CFD) multidimensional, se ejecutó un procedimiento de reducción de dos pasos utilizando el software KINALC. Esto implicó la eliminación esquelética de reacciones redundantes seguida de un Análisis de Componentes Principales (PCA) para abordar la multicolinealidad estructural en las tasas de reacción.

Resultados Clave

• Conversión y Rendimiento: El modelo predice una conversión neta de metano del 98.7%. Bajo la suposición idealizada de una sola fase (solo gas), los productos primarios son hidrógeno y acetileno (rendimiento del 88.0%), con fracciones menores de etileno (7.0%), propeno (1.1%) y butadieno (1.0%).
• Congelamiento Cinético Secuencial: El estudio identifica distintos umbrales térmicos que gobiernan la formación de productos durante el enfriamiento:
  • El acetileno se estabiliza y congela por debajo de 3900 K.
  • El etileno se acumula significativamente por debajo de 1520 K.
  • El propeno y el butadieno emergen a 1060 K y 870 K, respectivamente.
  • El bloqueo cinético: Todas las transformaciones moleculares cesan por debajo de 500 K, fijando la distribución final de los productos.
• Impacto de la Deposición de Hollín: Una simulación multi-etapa comparativa que incorpora la precipitación de carbono sólido revela que los modelos de fase gaseosa sobrepredicen significativamente los rendimientos de acetileno. Cuando se permite que el carbono sólido precipite a 2500 K, el 99% de la masa de carbono cae en forma de hollín. En este entorno enriquecido en hidrógeno y empobrecido en carbono, la química de enfriamiento posterior desplaza los productos finales de especies insaturadas (acetileno/etileno) hacia etano estable y metano no reaccionado, un resultado específico del modelo multifásico que contrasta con la línea base de fase gaseosa.
• Eficiencia Computacional: La reducción del mecanismo logró comprimir el esquema detallado (cientos de especies/pasos) en un mecanismo ligero de 26 especies y 102 reacciones irreversibles. Mientras que el resumen cita una aceleración de 61 veces, la Sección 3.5 reporta explícitamente un factor de aceleración de 14.2 veces para el perfil de ejecución. Este modelo reducido reproduce las concentraciones de especies dentro de un 0.1% del mecanismo completo.

Significancia y Reivindicaciones
El artículo afirma avanzar la capacidad predictiva de las simulaciones numéricas de reactores de plasma mediante el establecimiento de una estimación más rigurosa de los complejos perfiles de temperatura-tiempo experimentados por las mezclas de gases reactivos. Al desacoplar con precisión el núcleo de alta temperatura de la zona de enfriamiento rápido, el autor proporciona un marco físicamente realista para el modelado de la cinética química detallada.

La significancia principal radica en la demostración de que el congelamiento cinético es un proceso secuencial gobernado por umbrales térmicos específicos, lo que dicta la configuración final de los productos. Además, el estudio destaca que los cambios de fase heterogéneos (formación de hollín) alteran fundamentalmente las rutas de reacción en fase gaseosa, desplazando las distribuciones de los productos lejos de los mínimos termodinámicos predichos por los modelos de una sola fase.

Finalmente, el trabajo ofrece un marco computacionalmente eficiente para futuras aplicaciones de ingeniería. El mecanismo reducido, validado contra el esquema cinético completo, permite simulaciones de alta fidelidad en tres dimensiones de reactores de plasma industriales a escala completa, superando los costos computacionales prohibitivos asociados anteriormente con tales cinéticas químicas detalladas en códigos de dinámica de fluidos.