Stabilized Galerkin approximation of convection-diffusion-reaction equations: discrete maximum principle and convergence

We analyze a nonlinear shock-capturing scheme for H<sup>1</sup>-conforming, piecewise-affine finite element approximations of linear elliptic problems. The meshes are assumed to satisfy two standard conditions: a local quasi-uniformity property and the Xu-Zikatanov condition ensuring that the stiffness matrix associated with the Poisson equation is an M-matrix. A discrete maximum principle is rigorously established in any space dimension for convection-diffusion-reaction problems. We prove that the shock-capturing finite element solution converges to that without shock-capturing if the cell Peclet numbers are sufficiently small. Moreover, in the diffusion-dominated regime, the difference between the two finite element solutions super-converges with respect to the actual approximation error. Numerical experiments on test problems with stiff layers confirm the sharpness of the a priori error estimates


Published in:
Mathematics of Computation, 74, 252, 1637-52
Year:
2005
Keywords:
Other identifiers:
Laboratories:




 Record created 2007-04-24, last modified 2018-03-17


Rate this document:

Rate this document:
1
2
3
 
(Not yet reviewed)