3. Properties of the residue estimator#
We note |
\({\eta }_{\text{EX}}(K)\) the exact error \({\parallel u-{u}_{h}\parallel }_{{H}^{1}(K)}\) on the element \(K\) (unknown at first) |
and |
\({\eta }_{\text{EX}}(\Omega )\) the exact global error \({\parallel u-{u}_{h}\parallel }_{{H}^{1}(O)}\) |
We then have the following properties ([bib1]):
regardless of element \(K\), the elementary error \(\eta (K)\) is increased by the exact local error (multiplied by a constant independent of triangulation),
Be \(\forall K\eta (K)\le {C}_{1}\times {\eta }_{\text{EX}}(K)\)
the exact global error is increased by the global estimated error \(\eta (\Omega )\) (multiplied by a constant independent of \(T\))
Be \({\eta }_{\text{EX}}(\Omega )\le {C}_{2}\times \eta (\Omega )\)
The constants \({C}_{1}\text{et}{C}_{2}\) depend in principle on the type of finite element and on the boundary conditions of the problem. Kelly and Gago [bib3] proposed in 2D a constant \({C}_{2}\) that only depends on the degree \(p\) of the interpolation polynomial used:
\(\begin{array}{cc}{C}_{2}={(\frac{1}{\text{24}{p}^{2}})}^{1/2}\text{soit}& {C}_{2}=\frac{1}{\mathrm{2p}\sqrt{6}}\text{pour les TRIA3 et QUAD4}(\text{degré 1})\\ & {C}_{2}=\frac{1}{\mathrm{4p}\sqrt{6}}\text{pour les TRIA6 et QUAD8}(\text{degré 2})\end{array}\)
For 3D, we do not have a constant evaluation. However, it can be said that the global estimated error overestimates the exact global error in all cases. This result is not necessarily true at the local level.