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

Resumo Técnico: Modelagem Cinética de Reforma Seca de CH4/CO2 em uma Descarga de Barreira Dielétrica

Definição do Problema
A conversão de metano (CH4) e dióxido de carbono (CO2) via reforma seca (DRM) assistida por plasma em descargas de barreira dielétrica (DBD) oferece uma via não térmica para a produção de hidrogênio e gás de síntese (syngas). No entanto, as ferramentas de modelagem cinética para este processo frequentemente dependem de plataformas especializadas e não modulares e utilizam frequentemente constantes de velocidade ilustrativas, de ordem de magnitude, em vez de dados rigorosamente validados. Além disso, há uma carência de frameworks abertos que acoplem explicitamente a química homogênea 0D com o transporte de dispersão axial 1D, o que dificulta a prototipagem rápida de parâmetros operacionais ou a validação de modelos contra benchmarks experimentais com precisão quantitativa.

Metodologia
Os autores desenvolveram um framework Python aberto e modular (0D/1D) para simular a conversão de CH4/CO2 em uma DBD.

• Mecanismo Cinético: Foi empregado um mecanismo reduzido compreendendo 23 espécies neutras e 38 reações. Este inclui dissociação por impacto eletrônico, reações bimoleculares neutro-neutro, associações de três corpos e decomposições unimoleculares.
• Cinética Eletrônica: Diferente de iterações anteriores que utilizavam ajustes analíticos, os coeficientes de taxa de impacto eletrônico (E1–E5) foram derivados de uma resolução genuína da equação de Boltzmann de dois termos (usando o solver BOLOS) alimentada com seções de choque reais de espalhamento elétron-neutro da base de dados LXCat/Morgan.
• Forçamento de Plasma: A DBD foi modelada como um trem de pulsos periódicos representando microdescargas filamentares. Uma aproximação de média temporal foi aplicada para substituir as taxas instantâneas resolvidas por pulso por valores médios por ciclo, validada para reproduzir os resultados resolvidos por pulso com erro <1%, reduzindo significativamente o custo computacional.
• Análise de Sensibilidade: Uma análise de sensibilidade global de Morris foi realizada sobre as 33 constantes de taxa neutro-neutro. Esta identificou seis reações específicas como os principais drivers da conversão de CO2, da razão H2/CO e da eficiência energética.
• Estratégia de Validação: As seis reações críticas identificadas (R2, R11–R14 e R33) foram substituídas por valores provenientes da literatura avaliada de cinética de combustão e atmosférica. As outras 27 constantes de taxa foram mantidas como parâmetros ilustrativos de ordem de magnitude.
• Estrutura do Modelo: O framework acopla um modelo Lagrangiano 0D (rastreando um pacote de gás) com um modelo de fluxo de pistão de dispersão axial 1D Euleriano para simular perfis espaciais.

Principais Resultados

• Validação do Mecanismo: Substituir suposições ilustrativas por valores da literatura corrigiu erros que variaram de 6 vezes a mais de 3000 vezes. Especificamente, a reação CO + O (+M) → CO2 (+M) (R33), que é proibida por spin e difícil de medir, foi ancorada a um valor de literatura de 200 K (o melhor disponível, embora o modelo opere a 400 K).
• Tendências de Conversão: Com o mecanismo validado, a conversão de CO2 respeita o teto de ~25% observado em 19 fontes de plasma de pressão atmosférica na literatura para um SEI (Input de Energia Específica) de até ~30.000 J/L. A conversão de CH4 sobe de forma mais íngreme, aproximando-se de uma conversão quase completa em altos SEI.
• Seletividade: A razão H2/CO estabiliza-se entre 1,1 e 2,1, refletindo um regime limitado por CH4, onde o CH4 se dissocia mais rápido que o CO2, consistente com observações experimentais.
• Comparação com Benchmark:
  • Em um SEI correspondente de 5,7 eV/molécula, o modelo inicialmente sobreprestou a conversão de CH4 em ~5,8 vezes e a de CO2 em ~2,7 vezes em comparação ao benchmark experimental de Ozkan et al. [55].
  • Aplicar um SEI "corrigido por potência" (contabilizando a transferência de potência total em vez de apenas o custo de energia por impacto eletrônico) reduziu essas lacunas para ~3,3 vezes (CH4) e ~1,6 vezes (CO2).
  • Crucialmente, em valores de SEI mais baixos e tipicamente experimentais (938–2.500 J/L), o modelo mostrou concordância próxima com o experimento, subprestando levemente a conversão de CH4 e correspondendo-se de perto com a de CO2. A divergência aumenta apenas em altos SEI, sugerindo que o modelo é robusto no regime de baixo a moderado SEI típico de DBDs laboratoriais.
• Fontes de Incerteza: A análise de sensibilidade confirmou que, mesmo após a validação, a R33 permanece como a principal fonte de incerteza para a conversão de CO2 e a razão H2/CO. A sobreprestação restante da conversão de CH4 é atribuída à química do pool de radicais de CH4 não validada (R1, R3–R10) e à suposição de razões de ramificação resolvidas por canal para a dissociação de CH4.

Significância e Alegações
O artigo apresenta um passo significativo em direção a uma ferramenta de código aberto e quantitativamente validada para DRM de plasma. Suas principais contribuições são:

1. Framework Aberto: Fornecer um código Python reprodutível que acopla química 0D com transporte 1D, acessível para prototipagem rápida e ensino.
2. Upgrade Metodológico: Mover de ajustes analíticos arbitrários para coeficientes de taxa genuínos resolvidos por Boltzmann baseados em seções de choque reais.
3. Validação Direcionada: Demonstrar que a substituição de apenas seis constantes de taxa chave por valores da literatura, guiada pela sensibilidade, melhora significativamente a fidelidade do modelo, removendo uma grande fonte de erro e aproximando as previsões da realidade experimental.
4. Limitações Realistas: Os autores alegam modestamente que, embora o modelo ainda não seja totalmente quantitativo em todas as faixas de SEI, ele fornece uma base "qualitativa-para-semi-quantitativa" para o regime de SEI baixo a moderado. Eles identificam explicitamente o roteiro restante para uma validação total: obter dados de baixa temperatura para R33, validar a química de radicais de CH4 e calibrar os parâmetros de forçamento de plasma contra medições elétricas específicas do reator.

O trabalho não alega ter resolvido o problema do DRM quantitativamente, mas sim fornece uma base computacional transparente e aberta com um roteiro específico e delimitado para alcançar um modelo totalmente validado.