Fully Reliable Localized Error Control in the FEM

If the first task in numerical analysis is the calculation of an approximate solution, the second is to provide a guaranteed error bound and is often of equal importance. The standard approaches in the a posteriori error analysis of finite element methods suppose that the exact solution has a certain regularity or the numerical scheme enjoys some saturation property. For coarse meshes those asymptotic arguments are difficult to recast into rigorous error bounds. The aim of this paper is to provide reliable computable error bounds which are efficient and complete in the sense that constants are estimated as well. The main argument is a localization via a partition of unity which leads to problems on small domains. Two fully reliable estimates are established: The sharper one solves an analytical interface problem with residuals following Babuska and Rheinboldt [SIAM J. Numer. Anal., 15 (1978), pp. 736--754]. The second estimate is a modification of the standard residual-based a posteriori estimate with explicit constants from local analytical eigenvalue problems. For some class of triangulations we show that the efficiency constant is smaller than 2.5. According to our numerical experience, the overestimation of our computable estimates proved to be reasonably small, with an overestimation by a factor between 2.5 and 4 only.

[1]  P. Clément Approximation by finite element functions using local regularization , 1975 .

[2]  Rüdiger Verfürth,et al.  A posteriori error estimation and adaptive mesh-refinement techniques , 1994 .

[3]  C. Carstensen,et al.  Constants in Clément-interpolation error and residual based a posteriori estimates in finite element methods , 2000 .

[4]  B. S. Garbow,et al.  Matrix Eigensystem Routines — EISPACK Guide , 1974, Lecture Notes in Computer Science.

[5]  L. R. Scott,et al.  Finite element interpolation of nonsmooth functions satisfying boundary conditions , 1990 .

[6]  W. Rheinboldt,et al.  Error Estimates for Adaptive Finite Element Computations , 1978 .

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

[8]  L. Hörmander Linear Partial Differential Operators , 1963 .

[9]  Claes Johnson,et al.  Introduction to Adaptive Methods for Differential Equations , 1995, Acta Numerica.

[10]  B. S. Garbow,et al.  Matrix Eigensystem Routines — EISPACK Guide , 1974, Lecture Notes in Computer Science.

[11]  R. Verfürth,et al.  Edge Residuals Dominate A Posteriori Error Estimates for Low Order Finite Element Methods , 1999 .

[12]  J. Lions,et al.  Non-homogeneous boundary value problems and applications , 1972 .

[13]  Ivo Babuška,et al.  Analysis of the efficiency of an a posteriori error estimator for linear triangular finite elements , 1992 .

[14]  Piotr Hajłasz,et al.  @ 1996 Kluwer Academic Publishers. Printed in the Netherlands. Sobolev Spaces on an Arbitrary Metric Space , 1994 .