Premium
Weakly symmetric stress equilibration and a posteriori error estimation for linear elasticity
Author(s) -
Bertrand Fleurianne,
Kober Bernhard,
Moldenhauer Marcel,
Starke Gerhard
Publication year - 2021
Publication title -
numerical methods for partial differential equations
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.901
H-Index - 61
eISSN - 1098-2426
pISSN - 0749-159X
DOI - 10.1002/num.22741
Subject(s) - mathematics , estimator , linear elasticity , compressibility , finite element method , mathematical analysis , upper and lower bounds , elasticity (physics) , a priori and a posteriori , constant (computer programming) , mechanics , physics , thermodynamics , philosophy , statistics , epistemology , computer science , programming language
This paper proposes and analyzes a posteriori error estimator based on stress equilibration for linear elasticity with emphasis on the behavior for (nearly) incompressible materials. It is based on an H (div)‐conforming, weakly symmetric stress reconstruction from the displacement‐pressure approximation computed with a stable finite element pair. Our focus is on the Taylor‐Hood combination of continuous finite element spaces of polynomial degrees k + 1 and k for the displacement and the pressure, respectively. This weak symmetry allows us to prove that the resulting error estimator constitutes a guaranteed upper bound for the error with a constant that depends only on local constants associated with the patches and thus on the shape regularity of the triangulation. It does not involve global constants like those from Korn's in equality which may become very large depending on the location and type of the boundary conditions. Local efficiency, also uniformly in the incompressible limit, is deduced from the upper bound by the residual error estimator. Numerical results for the popular Cook's membrane test problem confirm the theoretical predictions.