Civil Engineering Reference
In-Depth Information
There are five popular approaches to building a posteriori estimators.
1. Residual estimators.
We bound the error on an element T in terms of the size of the residual R T
and the jumps R e on the edges e ∂T . These estimators are due to Babuska and
Rheinboldt [1978a].
2. Estimators based on local Neumann problems
On every triangle T we solve a local variational problem which is a discrete
analog of
z = R T
in T,
( 8 . 5 )
∂z
∂n = R e
on e ∂T.
We choose the approximating space to contain polynomials whose degrees are
higher than those in the underlying finite element space. These estimators are
obtained using the energy norm
z 1 ,T , and are due to Bank and Weiser [1985].
See also the comment before Theorem 9.5.
3. Estimators based on a local Dirichlet problem.
For every element T , we solve a variational problem on the set ω T :
z = f
in ω T ,
( 8 . 6 )
z = u h
on ∂ω T .
Again, we expand the approximating space to include polynomials of higher degree
than in the actual finite element space. Following Babuska and Rheinboldt [1978b],
the norm of the difference
z u h 1 T
provides an estimator.
4. Estimators based on averaging.
We construct a continuous approximation σ h of
u h by a two-step process.
At every node of the triangulation, let σ h be a weighted average of the gradients
u h on the neighboring triangles, where the weight is proportional to the areas
of the triangles. We then extend σ h to the whole element by linear interpolation.
Then following Zienkiewicz and Zhu [1987], we use the difference between
u h
and σ h as an estimator. An analysis without restrictive assumptions was done by
Rodriguez [1994] and by Carstensen and Bartels [2002].
5. Hierarchical estimators.
In principle the difference from a finite element approximation on an ex-
panded space is estimated. The difference can be estimated by using a strengthened
Cauchy inequality; see Deuflhard, Leinen, and Yserentant [1989]. The procedure
fits into Carl Runge's old and general concept (Runge's rule) . The error of a nu-
merical result is estimated by comparing it with the result of a more accurate
formula.
Search WWH ::




Custom Search