MÉTODO NUMÉRICO PARA SOLUÇÃO DE EDOs RÍGIDAS NA MODELAGEM DE REAÇÕES QUÍMICAS

Abstract

Problemas envolvendo reações químicas possuem uma característica importante denominada rigidez, indicando que as soluções dos sistemas de EDOs envolvidos variam em diferentes ordens de grandeza. Isso faz com que se escolham métodos numéricos adequados que tenham o intuito de obter soluções numéricas que sejam estáveis e convergentes a um custo computacional reduzido. Os métodos mais eficientes para tratar esse tipo de problema são os métodos implícitos, pois possuem uma região de estabilidade ilimitada no plano complexo que permite grandes variações no tamanho do passo. Neste trabalho propomos um método L-estável: o método de Rosenbrock de quarta ordem com quatro estágios. Simulam-se três problemas envolvendo reações químicas e os resultados apresentam boa concordância com a literatura.