A posteriori error estimators for stabilized P1 nonconforming approximation of the Stokes problem

• Published in: Computer methods in applied mechanics and engineering Vol. 199; no. 45; pp. 2903 - 2912
• Main Authors: Lee, Hyung-Chun, Kim, Kwang-Yeon
• Format: Journal Article
• Language: English
• Published: Kidlington Elsevier B.V 15.11.2010; Elsevier
• Subjects: A posteriori error estimator; Approximation; Balancing; Computational techniques; Equilibrated residual method; Errors; Estimators; Exact sciences and technology; Finite-element and galerkin methods; Fluid dynamics; Fluid flow; Fundamental areas of phenomenology (including applications); General theory; Mathematical analysis; Mathematical methods in physics; Mathematical models; Physics; Solid mechanics; Stabilized P1 nonconforming finite element method; Static elasticity (thermoelasticity...); Stokes problem; Stresses; Structural and continuum mechanics; Equilibrated residual method; 65N15; A posteriori error estimator; Stokes problem; Stabilized P1 nonconforming finite element method; 65N50; 65N30; Neumann problem; Partial differential equations; Error estimation; method; Non conforming finite element; Error bound; Automatic mesh generation; Stabilized P1 nonconforming finite element; Galerkin-Petrov method; Incompressible flow; Normal stress; A posteriori estimation; Boundary-value problems; Solvability
• ISSN: 0045-7825; 1879-2138
• DOI: 10.1016/j.cma.2010.06.002

In this paper we derive and analyze some a posteriori error estimators for the stabilized P1 nonconforming approximation of the Stokes problem involving the strain tensor. This will be done by decomposing the numerical error in a proper way into conforming and nonconforming contributions. The error estimator for the nonconforming error is obtained in the standard way, and the implicit error estimator for the conforming error is derived by applying the equilibrated residual method. A crucial part of this work is construction of approximate normal stresses on interelement boundaries which will serve as equilibrated Neumann data for local Stokes problems. It turns out that such normal stresses can be simply computed by local weak residuals of the discrete system plus jumps of the velocity solution and that a stronger equilibration condition is satisfied to ensure solvability of local Stokes problems. We also derive a simple explicit error estimator based on the nonsymmetric tensor recovery of the normal stress error. Numerical results are provided to illustrate the performance of our error estimators.