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

We analyze a nonlinear shock-capturing scheme for H-1-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.