Analysis of a new stabilized finite element method for the reaction–convection–diffusion equations with a large reaction coefficient

Abstract In this paper, we propose and analyze a new stabilized finite element method using continuous piecewise linear (or bilinear) elements for solving 2D reaction–convection–diffusion equations. The equation under consideration involves a small diffusivity e and a large reaction coefficient σ , leading to high Peclet number and high Damkohler number. In addition to giving error estimates of the approximations in L 2 and H 1 norms, we explicitly establish the dependence of error bounds on the diffusivity, the L ∞ norm of convection field, the reaction coefficient and the mesh size. Our analysis shows that the proposed method is particularly suitable for problems with a small diffusivity and a large reaction coefficient, or more precisely, with a large mesh Peclet number and a large mesh Damkohler number. Several numerical examples exhibiting boundary or interior layers are given to illustrate the high accuracy and stability of the proposed method. The results obtained are also compared with those of existing stabilization methods.

[1]  Frédéric Valentin,et al.  Multiscale and Residual-Free Bubble Functions for Reaction-Advection-Diffusion Problems , 2005 .

[2]  L. R. Scott,et al.  The Mathematical Theory of Finite Element Methods , 1994 .

[3]  Howard C. Elman,et al.  Finite Elements and Fast Iterative Solvers: with Applications in Incompressible Fluid Dynamics , 2014 .

[4]  Isaac Harari,et al.  Stability of semidiscrete formulations for parabolic problems at small time steps , 2004 .

[5]  Arif Masud,et al.  Revisiting stabilized finite element methods for the advective–diffusive equation , 2006 .

[6]  Martin Stynes,et al.  Steady-state convection-diffusion problems , 2005, Acta Numerica.

[7]  Pavel B. Bochev,et al.  On stabilized finite element methods for the Stokes problem in the small time step limit , 2007 .

[8]  T. Hughes,et al.  The Galerkin/least-squares method for advective-diffusive equations , 1988 .

[9]  M. H. Doweidar,et al.  COMBINING ADJOINT STABILIZED METHODS FOR THE ADVECTION-DIFFUSION-REACTION PROBLEM , 2007 .

[10]  Philippe G. Ciarlet,et al.  The finite element method for elliptic problems , 2002, Classics in applied mathematics.

[11]  C. Farhat,et al.  Bubble Functions Prompt Unusual Stabilized Finite Element Methods , 1994 .

[12]  L. Franca,et al.  On an Improved Unusual Stabilized Finite Element Method for theAdvective-Reactive-Diffusive Equation , 1999 .

[13]  Lutz Tobiska,et al.  Numerical Methods for Singularly Perturbed Differential Equations , 1996 .

[14]  Volker John,et al.  On Finite Element Methods for 3D Time-Dependent Convection-Diffusion-Reaction Equations with Small Diffusion , 2008 .

[15]  L. Franca,et al.  Stabilized finite element methods. II: The incompressible Navier-Stokes equations , 1992 .

[16]  Alessandro Russo,et al.  CHOOSING BUBBLES FOR ADVECTION-DIFFUSION PROBLEMS , 1994 .

[17]  Vivette Girault,et al.  Finite Element Methods for Navier-Stokes Equations - Theory and Algorithms , 1986, Springer Series in Computational Mathematics.

[18]  T. Hughes,et al.  Stabilized finite element methods. I: Application to the advective-diffusive model , 1992 .

[19]  Po-Wen Hsieh,et al.  A Novel Least-Squares Finite Element Method Enriched with Residual-Free Bubbles for Solving Convection-Dominated Problems , 2010, SIAM J. Sci. Comput..

[20]  Alessandro Russo,et al.  Further considerations on residual-free bubbles for advective - diffusive equations , 1998 .

[21]  Po-Wen Hsieh,et al.  On efficient least-squares finite element methods for convection-dominated problems , 2009 .

[22]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[23]  T. Hughes Multiscale phenomena: Green's functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods , 1995 .

[24]  Claes Johnson Numerical solution of partial differential equations by the finite element method , 1988 .