Kinetic Modelling of CH4/CO2 Dry Reforming in a Dielectric Barrier Discharge for Hydrogen Production: An Open, Sensitivity- Guided Validation Study

Resumen Técnico: Modelado Cinético de la Reformación Seca de CH4/CO2 en una Descarga de Barrera Dieléctrica

Planteamiento del Problema
La conversión de metano (CH4) y dióxido de carbono (CO2) mediante la reformación seca (DRM) asistida por plasma en descargas de barrera dieléctrica (DBD) ofrece una vía no térmica para producir hidrógeno y gas de síntesis. Sin embargo, las herramientas de modelado cinético para este proceso suelen depender de plataformas especializadas y no modulares, y frecuentemente utilizan constantes de velocidad ilustrativas de orden de magnitud en lugar de datos rigurosamente validados. Además, existe una falta de marcos abiertos que acoplen explícitamente la química homogénea 0D con el transporte de dispersión axial 1D, lo que dificulta el prototipado rápido de parámetros operativos o la validación de modelos frente a puntos de referencia experimentales con precisión cuantitativa.

Metodología
Los autores desarrollaron un marco modular y abierto en Python (0D/1D) para simular la conversión de CH4/CO2 en una DBD.

• Mecanismo Cinético: Se empleó un mecanismo reducido que comprende 23 especies neutras y 38 reacciones. Este incluye la disociación por impacto electrónico, reacciones neutro-neutras bimoleculares, asociaciones de tres cuerpos y descomposiciones unimoleculares.
• Cinética Electrónica: A diferencia de iteraciones anteriores que utilizaban ajustes analíticos, los coeficientes de velocidad de impacto electrónico (E1–E5) se derivaron de la resolución de una ecuación de Boltzmann de dos términos genuina (usando el solver BOLOS) alimentada con secciones eficaces reales de choque electrón-neutro de la base de datos LXCat/Morgan.
• Forzamiento del Plasma: La DBD se modeló como un tren de pulsos periódicos que representa microdescargas filamentosas. Se aplicó una aproximación de promedio temporal para reemplazar las tasas instantáneas resueltas por pulsos con valores promediados por ciclo, validados para reproducir los resultados resueltos por pulso con un error <1%, reduciendo significamente el costo computacional.
• Análisis de Sensibilidad: Se realizó un análisis de sensibilidad global de Morris sobre las 33 constantes de velocidad neutro-neutras. Esto identificó seis reacciones específicas como los principales impulsores de la conversión de CO2, la relación H2/CO y la eficiencia energética.
• Estrategia de Validación: Las seis reacciones críticas identificadas (R2, R11–R14 y R33) fueron reemplazadas con valores de la literatura evaluada de cinética de combustión y atmosférica. El resto de las 27 constantes de velocidad se mantuvieron como parámetros ilustrativos de orden de magnitud.
• Estructura del Modelo: El marco acopla un modelo Lagrangiano 0D (que rastrea una parcela de gas) con un modelo de reactor de flujo pistón de dispersión axial 1D de tipo Euler para simular perfiles espaciales.

Resultados Clave

• Validación del Mecanismo: El reemplazo de conjeturas ilustrativas por valores de la literatura corrigió errores que variaban de 6 veces a más de 3000 veces. Específicamente, la reacción CO + O (+M) → CO2 (+M) (R33), que es prohibida por espín y difícil de medir, se ancló a un valor de la literatura de 200 K (el mejor disponible, aunque el modelo opera a 400 K).
• Tendencias de Conversión: Con el mecanismo validado, la conversión de CO2 respeta el techo de ~25% observado en 19 fuentes de plasma a presión atmosférica en la literatura para una Energía de Entrada Específica (SEI) de hasta ~30,000 J/L. La conversión de CH4 aumenta de forma más pronunciada, acercándose a una conversión casi completa a altas SEI.
• Selectividad: La relación H2/CO se estabiliza entre 1.1 y 2.1, lo que refleja un régimen limitado por CH4 donde el CH4 se disocia más rápido que el CO2, siendo consistente con las observaciones experimentales.
• Comparación de Referencia:
  • Con una SEI de 5.7 eV/molécula, el modelo inicialmente sobrepredijo la conversión de CH4 en ~5.8 veces y la de CO2 en ~2.7 veces en comparación con el punto de referencia experimental de Ozkan et al. [55].
  • La aplicación de una SEI "corregida por potencia" (que tiene en cuenta la transferencia de potencia total en lugar de solo el costo de energía por impacto electrónico) redujo estas brechas a ~3.3 veces (CH4) y ~1.6 veces (CO2).
  • Crucialmente, en valores de SEI más bajos y típicos de experimentos (938–2,500 J/L), el modelo mostró una concordancia cercana con el experimento, subprediciendo ligeramente la conversión de CH4 y coincidiendo estrechamente con la de CO2. La divergencia aumenta solo a altas SEI, lo que sugiere que el modelo es robusto en el régimen de baja a moderada SEI típico de las DBD de laboratorio.
• Fuentes de Incertidumbre: El análisis de sensibilidad confirmó que, incluso tras la validación, R33 sigue siendo la fuente dominante de incertidumbre para la conversión de CO2 y la relación H2/CO. La sobrepredicción restante de la conversión de CH4 se atribuye a la química del reservorio de radicales de CH4 no validada (R1, R3–R10) y al supuesto de relaciones de ramificación resueltas por canal para la disociación de CH4.

Significancia y Reivindicaciones
El artículo presenta un paso significativo hacia una herramienta de código abierto y cuantitativamente validada para la DRM por plasma. Sus contribuciones primarias son:

1. Marco Abierto: Proporcionar un código base en Python reproducible que acopla la química 0D con el transporte 1D, accesible para el prototipado rápido y la enseñanza.
2. Mejora Metodológica: Pasar de ajustes analíticos arbitrarios a coeficientes de velocidad genuinos resueltos por Boltzmann basados en secciones eficaces reales.
3. Validación Dirigida: Demostrar que el reemplazo de solo seis constantes de velocidad clave con valores de la literatura, guiado por la sensibilidad, mejora significativamente la fidelidad del modelo, eliminando una fuente importante de error y acercando las predicciones a la realidad experimental.
4. Limitaciones Realistas: Los autores afirman modestamente que, si bien el modelo aún no es totalmente cuantitativo en todos los rangos de SEI, proporciona una base "cualitativa a semi-cuantitativa" para el régimen de SEI de bajo a moderado. Identifican explícitamente la hoja de ruta restante para una validación completa: obtener datos de menor temperatura para R33, validar la química de los radicales de CH4 y calibrar los parámetros de forzamiento del plasma frente a mediciones eléctricas específicas del reactor.

El trabajo no pretende haber resuelto el problema de la DRM de manera cuantitativa, sino que proporciona una base computacional abierta y transparente con una hoja de ruta específica y delimitada para lograr un modelo totalmente validado.