A priori and a posteriori analysis of non-conforming finite elements with face penalty for advection-diffusion equations

We analyse a non-conforming finite-element method to approximate advection-diffusion-reaction equations. The method is stabilized by penalizing the jumps of the solution and those of its advective derivative across mesh interfaces. The a priori error analysis leads to (quasi-)optimal estimates in the mesh size (sub-optimal by order ½ in the L2-norm and optimal in the broken graph norm for quasi-uniform meshes) keeping the Péclet number fixed. Then, we investigate a residual a posteriori error estimator for the method. The estimator is semi-robust in the sense that it yields lower and upper bounds of the error which differ by a factor equal at most to the square root of the Péclet number. Finally, to illustrate the theory we present numerical results including adaptively generated meshe

A posteriori energy-norm error estimates for advection-diffusion equations approximated by weighted interior penalty methods

22 pagesWe propose and analyze a posteriori energy-norm error estimates for weighted interior penalty discontinuous Galerkin approximations to advection-diffusion-reaction equations with heterogeneous and anisotropic diffusion. The weights, which play a key role in the analysis, depend on the diffusion tensor and are used to formulate the consistency terms in the discontinuous Galerkin method. The error upper bounds, in which all the constants are specified, consist of three terms: a residual estimator which depends only on the elementwise fluctuation of the discrete solution residual, a diffusive flux estimator where the weights used in the method enter explicitly, and a non-conforming estimator which is nonzero because of the use of discontinuous finite element spaces. The three estimators can be bounded locally by the approximation error. For moderate advection, it is shown that full robustness with respect to diffusion heterogeneities is achieved owing to the specific design of the weights in the discontinuous Galerkin method, while diffusion anisotropies remain purely local and impact the constants through the square root of the condition number of the diffusion tensor. For dominant advection, it is shown, in the spirit of previous work by Verfürth on continuous finite elements, that the constants are bounded by the square root of the local Péclet number