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

Resumo Técnico: Modelagem Numérica de Pirólise de Metano em Alta Temperatura sob Gradientes Térmicos Inspirados por Plasma

Definição do Problema
O artigo aborda o desafio de converter reservas remotas de gás natural (principalmente metano) em hidrocarbonetos líquidos transportáveis. Embora a tecnologia de plasma de micro-ondas ofereça uma alternativa direta às rotas indiretas tradicionais e energeticamente intensivas (como a reforma a vapor do metano seguida pela síntese de Fischer-Tropsch), ela enfrenta obstáculos operacionais significativos. Especificamente, a formação simultânea de partículas de fuligem sólida em zonas de plasma de alta temperatura leva ao entupimento do reator e à redução da seletividade para os alvos de matéria-prima líquida. Além disso, as simulações numéricas existentes frequentemente simplificam excessivamente os complexos gradientes térmicos encontrados em reatores de plasma, assumindo zonas térmicas uniformes ou aproximações unidimensionais básicas. Essas simplificações falham em capturar as verdadeiras trajetórias químicas das moléculas, levando a previsões imprecisas de rendimentos de produtos e taxas de formação de fuligem. O problema central é a falta de um arcabouço preditivo rigoroso que consiga desacoplar com precisão o núcleo de plasma de estado estacionário de alta temperatura da zona de resfriamento rápido (quenching) na saída do reator, mantendo-se computacionalmente eficiente o suficiente para o acoplamento com a dinâmica de fluidos multidimensional.

Metodologia
O estudo emprega um arcabouço cinético de pirólise térmica para modelar a conversão de metano em um ambiente de reator de pressão atmosférica, excluindo explicitamente a cinética de elétrons em favor de redes de radicais neutros sob gradientes térmicos agressivos.

• Expansão do Modelo Cinético: O mecanismo fundamental baseia-se no mecanismo de Dean para a pirólise do metano, expandido para cobrir o espectro de temperatura de 300–4000 K típico de arcos de plasma de micro-ondas. Uma adição crítica é a inclusão do cruzamento entre sistemas (intersystem crossing - ISC) induzido por colisão do metileno singlete eletronicamente excitado ($^1CH_2$) para o metileno triplete ($^3CH_2$), impulsionado por colisões com o gás de banho de argônio. Esta transição de mecânica quântica influencia significamente o desenvolvimento do reservatório de radicais, particularmente durante o resfriamento térmico rápido.
• Desacoplamento do Perfil Térmico: O modelo separa o reator em duas zonas distintas:
  1. Núcleo de Plasma Ativo: Uma zona de aquecimento contínuo com um tempo de residência de aproximadamente 15,54 ms (notado como 155,4 s na Tabela 3, provavelmente um erro de formatação no texto de origem) mantida entre 2500 K e 4000 K, onde a absorção de energia eletromagnética conduz o sistema a um estado de quase-estacionário de espécies de radicais altamente excitadas.
  2. Região de Jato de Resfriamento: Uma zona de quench limitada por difusão onde o gás sofre um resfriamento rápido de 2500 K para 300 K ao longo de aproximadamente 76 ms (especificamente 76,9 ms, conforme listado na Tabela 3). Este perfil é derivado pela resolução de uma equação de condução de calor radial 1D, validada contra dados experimentais de espalhamento Rayleigh e OES.
• Negligenciabilidade da Cinética de Elétrons: Uma análise de Boltzmann zero-dimensional confirma que, a 1 bar e 2500 K, as taxas de dissociação térmica excedem as taxas de dissociação por impacto eletrônico em sete ordens de magnitude. Consequentemente, o modelo trata o sistema como um processo de pirólise puramente térmica governado pela química de radicais neutros.
• Redução de Mecanismo: Para permitir a integração em Dinâmica de Fluidos Computacional (CFD) multidimensional, um procedimento de redução em duas etapas foi executado usando o software KINALC. Isso envolveu a eliminação esquelética de reações redundantes seguida por Análise de Componentes Principais (PCA) para abordar a multicolinearidade estrutural nas taxas de reação.

Principais Resultados

• Conversão e Rendimento: O modelo prevê uma conversão líquida de metano de 98,7%. Sob a suposição idealizada de fase única (apenas gás), os produtos primários são hidrogênio e acetileno (rendimento de 88,0%), com frações menores de etileno (7,0%), propeno (1,1%) e butadieno (1,0%).
• Congelamento Cinético Sequencial: O estudo identifica limiares térmicos distintos que governam a formação de produtos durante o quench:
  • O acetileno estabiliza e congela abaixo de 3900 K.
  • O etileno acumula-se significativamente abaixo de 1520 K.
  • O propeno e o butadieno emergem em 1060 K e 870 K, respectivamente.
  • O Bloqueio Cinético: Todas as transformações moleculares cessam abaixo de 500 K, fixando a distribuição final dos produtos.
• Impacto da Deposição de Fuligem: Uma simulação comparativa de múltiplas etapas incorporando a precipitação de carbono sólido revela que modelos de fase gasosa superestimam significativamente os rendimentos de acetileno. Quando o carbono sólido é permitido precipitar a 2500 K, 99% da massa de carbono decanta como fuligem. Neste ambiente rico em hidrogênio e empobrecido em carbono, a química subsequente de quenching altera os produtos finais de espécies insaturadas (acetileno/etileno) para etano estável e metano não reagido, um resultado específico do modelo multifásico que contrasta com a linha de base de fase gasosa.
• Eficiência Computacional: A redução do mecanismo conseguiu comprimir o esquema detalhado (centenas de espécies/etapas) em um mecanismo enxuto de 26 espécies e 102 reações irreversíveis. Embora o resumo cite uma aceleração de 61 vezes, a Seção 3.5 reporta explicitamente um fator de aceleração de 14,2 vezes para o perfil de execução. O modelo reduzido reproduz as concentrações de espécies dentro de 0,1% do mecanismo completo.

Significância e Alegações
O artigo afirma avançar a capacidade preditiva de simulações de reatores de plasma numéricos ao estabelecer uma estimativa mais rigorosa dos complexos perfis de temperatura-tempo experimentados por misturas gasosas reagentes. Ao desacoplar com precisão o núcleo de alta temperatura da zona de quench rápido, o autor fornece um arcabouço fisicamente realista para a modelagem de cinética química detalhada.

A principal significância reside na demonstração de que o congelamento cinético é um processo sequencial governado por limiares térmicos específicos, o que dita a configuração final dos produtos. Além disso, o estudo destaca que as mudanças de fase heterogêneas (formação de fuligem) alteram fundamentalmente as vias de reação em fase gasosa, deslocando a distribuição de produtos para longe dos mínimos termodinâmicos previstos por modelos de fase única.

Finalmente, o trabalho oferece um arcabouço computacionalmente eficiente para futuras aplicações de engenharia. O mecanismo reduzido, validado contra o esquema cinético completo, permite simulações de alta fidelidade e três dimensões de reatores de plasma industriais em escala total, superando os custos computacionais proibitivos anteriormente associados a tais cinéticas químicas detalhadas em códigos de dinâmica de fluidos.