Talks and Poster Presentations (with Proceedings-Entry):
M. Feischl, M. Page, D. Praetorius:
"Convergence of adaptive FEM for elliptic obstacle problems";
Talk: GAMM Jahrestagung 2011,
Graz;
04-18-2011
- 04-21-2011; in: "PAMM: Proceedings in Applied Mathematics and Mechanics",
PAMM,
11
(2011),
ISSN: 1617-7061;
767
- 768.
English abstract:
We treat the convergence of adaptive lowest-order FEM for some elliptic obstacle problem with affine obstacle χ, i.e.
(\nabla u, \nabla (u − v)) ≤ (f, u − v)
for all v ∈ A := {v ∈ H^1_0 (Ω) : v ≥ χ a.e in Ω}. For steering the adaptive mesh-refinement, we use a residual error estimator
\eta_\ell := ( \rho_\ell^2 + osc_\ell^2 )^1/2
that consists of the estimator from [Braess et al., Numer. Math. ´07] and additionally controls the data oscillations. We extend recent ideas from [Cascon et al., Numer. Anal. ´08] for the unrestricted variational problem to overcome the lack of Galerkin orthogonality. The main result states that an appropriately weighted sum
\Delta_\ell := \epsilon_\ell + \gamma \mu_\ell^2
of energy error, edge residuals, and data oscillations satisfies a contraction property
\Delta_{\ell+1} ≤ \kappa \Delta_\ell for all \ell \in \N
within each step of the adaptive loop. Here, \mu_\ell denotes a second error estimator which is equivalent to \eta_\ell and
0 < γ, κ < 1. This result is superior to a prior result from [Braess
et al., Numer. Math. ´07] in two ways: First, it is unnecessary to control the decay of the
data oscillations explicitly. Second, our analysis avoids the use of discrete local efficiency so that the local mesh-refinement is fairly arbitrary. In addition, we discuss the generalization to the case of inhomogeneous Dirichlet data and non-affine obstacles χ ∈ H^2 (Ω) and obtain similar results.
"Official" electronic version of the publication (accessed through its Digital Object Identifier - DOI)
http://dx.doi.org/10.1002/pamm.201110373
Created from the Publication Database of the Vienna University of Technology.