Otimização estocástica de um mecanismo de reação de óxido de urânio usando medições de reator de fluxo de plasma

Jan 10, 2024

Scientific Reports volume 13, Número do artigo: 9293 (2023) Citar este artigo

Detalhes das métricas

Neste trabalho, uma abordagem acoplada do Algoritmo Genético de Monte Carlo (MCGA) é usada para otimizar um mecanismo de reação de óxido de urânio em fase gasosa baseado em medições de reator de fluxo de plasma (PFR). O PFR produz um plasma constante de Ar contendo espécies de U, O, H e N com regiões de alta temperatura (3.000–5.000 K) relevantes para observar a formação de UO por meio de espectroscopia de emissão óptica. Um tratamento cinético global é usado para modelar a evolução química no PFR e produzir sinais de emissão sintéticos para comparação direta com experimentos. O espaço de parâmetros de um mecanismo de reação de óxido de urânio é então explorado por meio de amostragem de Monte Carlo usando funções objetivas para quantificar a concordância modelo-experiência. Os resultados de Monte Carlo são subsequentemente refinados usando um algoritmo genético para obter um conjunto de caminhos de reação e coeficientes de velocidade corroborados experimentalmente. Dos 12 canais de reação direcionados para otimização, quatro canais são bem restritos em todas as execuções de otimização, enquanto outros três canais são restritos em casos selecionados. Os canais otimizados destacam a importância do radical OH na oxidação do urânio no PFR. Este estudo compreende um primeiro passo para a produção de um mecanismo de reação abrangente experimentalmente validado para a formação de espécies moleculares de urânio em fase gasosa.

A cinética de reação de óxidos metálicos na fase gasosa é de grande relevância para muitos campos de pesquisa, incluindo astrofísica, ciência da combustão, engenharia nuclear e química de materiais em ambientes extremos. Nos últimos anos, este último campo produziu numerosos trabalhos experimentais e computacionais sobre a química do vapor de óxido de urânio (\({{\mathrm{UO_x}}}\))1. Produtos de fase gasosa de óxidos refratários, como \({{\mathrm{UO_x}}}\), têm sido historicamente difíceis de produzir devido às altas temperaturas de vaporização dos óxidos originais. Mais recentemente, os sistemas de plasma térmico forneceram um caminho para a produção imediata de metais na fase gasosa e para o estudo de sua química em ambientes reativos. No entanto, os tempos de extinção rápidos, a presença de radicais de fundo e a formação de óxidos intermediários voláteis em tais sistemas dificultam o isolamento de canais de reação específicos para estudo. Problemas semelhantes surgem em outros sistemas reativos de alta temperatura, como combustíveis de combustão de metal. Como resultado, os mecanismos de oxidação de metais em fase gasosa são frequentemente baseados em dados experimentais esparsos e estimativas teóricas de primeira ordem, como para a formação de óxido de alumínio2,3,4. Da mesma forma, um mecanismo de reação \({{\mathrm{UO_x}}}\) foi construído usando uma metodologia comparável em nosso trabalho anterior5. Embora tais mecanismos produzam resultados qualitativamente razoáveis ​​que podem se alinhar com alguns observáveis ​​experimentais, a validação experimental detalhada é difícil de alcançar. Esta etapa de validação é crucial para garantir que o mecanismo cinético químico possa ser usado de maneira preditiva para informar modelos subsequentes. Aqui, exploramos um método de inferir coeficientes de taxa de óxido de urânio (\({{\mathrm{UO_x}}}\)) com base em medições experimentais de um sistema de plasma térmico.

Devido à natureza fortemente acoplada e não linear da cinética química em plasmas de urânio, a extração de informações sobre a taxa de reação requer a resolução de um problema de otimização. Neste problema, os parâmetros subjacentes do modelo (coeficientes de taxa) são determinados com base nas saídas observadas (ou seja, informações espectroscópicas). Resolver tal problema por métodos determinísticos baseados em gradientes é difícil devido ao espaço de parâmetros potencialmente complexo com numerosos mínimos locais. Nesse caso, deve-se utilizar um método de otimização capaz de explorar continuamente todo o espaço de parâmetros enquanto localiza o mínimo global. Um desses métodos usados ​​anteriormente para problemas de cinética química é o algoritmo genético de Monte Carlo (MCGA)6. Esta técnica é adequada para o problema atual devido à sua eficácia em evitar a convergência em mínimos locais e sua facilidade de implementação. Independentemente da metodologia, a solução de um problema de otimização requer avaliações repetidas do modelo associado, muitas vezes numerando de milhares a milhões de execuções. Embora tempos computacionais razoáveis ​​sejam alcançados ao resolver a cinética química em um sistema espacialmente uniforme, o problema rapidamente se torna inviável quando a química é acoplada ao transporte de fluidos complexos. Essa consideração se torna importante na escolha de um sistema experimental para informar o problema de otimização.