SOME NUMERICAL STUDIES OF OSCILLATING CHEMICAL REACTIONS USING DISCONTINUOUS FINITE ELEMENTS

Abstract

La modelación numérica de la interacción entre reacciones químicas oscilantes y procesos difusivos requiere métodos numéricos sofisticados. La rigidez de estos modelos es eludida usando un método de Runge-Kutta diagonal para incrementar la estabilidad, y el tiempo de paso es controlado con un esquema incrustado de orden inferior. Dado que los patrones espaciales involucran grandes cambios de concentración, conviene realizar la discretización espacial con el método de Galerkin discontinuo (dG). Adicionalmente, mejoramos el método por medio de un precondicionamiento de multimalla. La eficiencia es evaluada con una reacción de un solo paso con solución exacta y se muestran algunos resultados con un modelo típico de reacción oscilante, el Brusselador.

Keywords

• difusión
• método de Galerkin discontinuo
• reacción química oscilante