The reliability of local error estimators for convection-diffusion equations

Abstract

We assess the reliability of a simple a posteriori error estimator for steady state convection-diffusion equations in cases where convection dominates. Our estimator is computed by solving a local Poisson problem with Neumann boundary conditions. It gives global upper and local lower bounds on the error measured in the semi-norm, except that the error may be over-estimated locally within boundary layers if these are not resolved by the mesh, that is, when the local mesh Péclet number is significantly greater than unity. We discuss the implications of this over-estimation in a practical context where the estimator is used as a local error indicator within a self-adaptive mesh refinement process.

This work was supported by EPSRC grants GR/K91262 and GR/L05617.