Sharp algebraic and total a posteriori error bounds for $h$ and $p$ finite elements via a multilevel approach. Recovering mass balance in any situation - Inria - Institut national de recherche en sciences et technologies du numérique Accéder directement au contenu
Article Dans Une Revue Computer Methods in Applied Mechanics and Engineering Année : 2020

Sharp algebraic and total a posteriori error bounds for $h$ and $p$ finite elements via a multilevel approach. Recovering mass balance in any situation

Résumé

We present novel H(div) and H1 liftings of given piecewise polynomials over a hierarchy of simplicial meshes, based on a global solve on the coarsest mesh and on local solves on patches of mesh elements around vertices on subsequent mesh levels. This in particular allows to lift a given algebraic residual. In connection with approaches lifting the total residual, we show how to obtain guaranteed, fully computable, and constant-free upper and lower a posteriori bounds on the algebraic, total, and discretization errors; here we consider the model Poisson equation discretized by the conforming finite element method of arbitrary order and including an arbitrary iterative solver. We next formulate safe stopping criteria ensuring that the algebraic error does not dominate the total error. We also prove efficiency, i.e., equivalence of our upper total and algebraic estimates with the total and algebraic errors, respectively, up to a generic constant; this constant is polynomial-degree-independent for the total error. Numerical experiments illustrate sharp control of all error components and accurate prediction of their spatial distribution in several test problems, including cases where some classical estimators fail. The H(div)-liftings at the same time allow to recover mass balance for any problem, any numerical discretization, and any situation such as inexact solution of (nonlinear) algebraic systems or algorithm failure, which we believe is of independent interest. We demonstrate this mass balance recovery in a simulation of immiscible incompressible two-phase flow in porous media.
Fichier principal
Vignette du fichier
MG_FE.pdf (4.01 Mo) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)
Loading...

Dates et versions

hal-01662944 , version 1 (13-12-2017)
hal-01662944 , version 2 (07-02-2018)
hal-01662944 , version 3 (24-12-2019)
hal-01662944 , version 4 (27-04-2020)
hal-01662944 , version 5 (01-09-2020)

Identifiants

Citer

Jan Papež, Ulrich Rüde, Martin Vohralík, Barbara I Wohlmuth. Sharp algebraic and total a posteriori error bounds for $h$ and $p$ finite elements via a multilevel approach. Recovering mass balance in any situation. Computer Methods in Applied Mechanics and Engineering, 2020, 371, pp.113243. ⟨10.1016/j.cma.2020.113243⟩. ⟨hal-01662944v5⟩
1248 Consultations
777 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More