Computable a posteriori error estimates in the finite element method based on its local conservativity: improvements using local minimization

¹ Univ. Grenoble and CNRS, Laboratoire Jean Kuntzmann, 51 rue des Mathématiques, 38400 Saint Martin d'Hères & INRIA Grenoble-Rhône-Alpes, Inovallée, 655 avenue de l'Europe, Montbonnot 38334 Saint Ismier Cedex, France;
² Department of Mathematics, Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University in Prague, Trojanova 13, 12000 Prague, Czech Republic;
³ Departamento de Matemática, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Intendente Güiraldes 2160, Ciudad Universitaria, C1428EGA, Argentina;
⁴ UPMC Univ. Paris 06, UMR 7598, Laboratoire Jacques-Louis Lions, 75005, Paris, France & CNRS, UMR 7598, Laboratoire Jacques-Louis Lions, 75005, Paris, France;

Corresponding authors: ibrahim.cheddadi@imag.fr fucik@fjfi.cvut.cz mprieto@dm.uba.ar vohralik@ann.jussieu.fr

Abstract

We investigate in this paper improvements of the a posteriori error estimates in the finite element method discretization of the Poisson equation, introduced in [M. Vohralík, A posteriori error estimation in the conforming finite element method based on its local conservativity and using local minimization, C. R. Math. Acad. Sci. Paris 346 (2008), 687–690] and [M. Vohralík, Guaranteed and fully robust a posteriori error estimates for conforming discretizations of diffusion problems with discontinuous coefficients, submitted]. The estimates presented in these references are guaranteed in the sense that they feature no undetermined constants and fully computable but numerical experiments show that the effectivity index, i.e., the ratio of the estimated and actual error, does not approach the optimal value of one but rather a slightly bigger value. We identify in this paper the reason for this and introduce a possible remedy, which consists in performing a local minimization of the values of the estimators over patches of simplicial submesh elements. We then present a set of numerical experiments showing the improvements achieved and compare our estimators, both theoretically and numerically, with the classical residual ones.

Résumé

Nous étudions dans cet article des améliorations des estimations a posteriori pour la discrétisation de l'équation de Poisson par la méthode des éléments finis, introduites dans [M. Vohralík, A posteriori error estimation in the conforming finite element method based on its local conservativity and using local minimization, C. R. Math. Acad. Sci. Paris 346 (2008), 687–690] et [M. Vohralík, Guaranteed and fully robust a posteriori error estimates for conforming discretizations of diffusion problems with discontinuous coefficients, soumis]. Les estimations présentées dans ces références sont garanties au sens où elles ne comportent pas de constantes indéterminées et sont entièrement calculables, mais des essais numériques montrent que l'indice d'efficacité, donné par le rapport entre l'erreur estimée et l'erreur exacte, n'approche pas la valeur optimale de un mais une valeur légèrement plus élevée. Dans cet article, nous en identifions la raison et nous introduisons une correction basée sur une minimisation locale des valeurs des estimateurs sur des ensembles de triangles ou de tétraèdres d'un sous-maillage. Nous présentons ensuite une série d'essais numériques et comparons, tant du point de vue théorique que numérique, nos estimations avec les estimations classiques par résidu.