SUPG-YZβ computation of chemically reactive convection-dominated nonlinear models

Abstract

In this computational study, we deal with the stabilized finite element solutions of convection-dominated models having nonlinear reaction mechanisms. The existence of advection terms in these models causes the numerical solutions obtained by standard discretization methods to exhibit nonphysical oscillations. The Galerkin finite element method is stabilized using the Streamline-Upwind/Petrov–Galerkin formulation to avoid such spurious oscillations. The stabilized formulation is also complemented with the YZβ shock-capturing technique to resolve steep gradients and discontinuities accurately. The nonlinear equation systems arising from the spatial discretizations are solved with the Newton–Raphson method supplemented with the ILU-preconditioned GMRES search technique. The proposed methods are tested on a comprehensive set of nonlinear reactive models. Numerical experiments show that the presented methods and techniques eliminate spurious oscillations significantly and resolve strong gradients accurately. All computations are performed in the FEniCS environment.