Convergence rates for an adaptive dual weighted residual finite element algorithm

K. S. Moon*, E. Von Schwerin, A. Szepessy, Raul Tempone

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

7 Scopus citations

Abstract

Basic convergence rates are established for an adaptive algorithm based on the dual weighted residual error representation, error = Σ elements error density x mesh size2+d, applied to isoparametric d-linear quadrilateral finite element approximation of functionals of multi scale solutions to second order elliptic partial differential equations in bounded domains of ℝd. In contrast to the usual aim to derive an a posteriori error estimate, this work derives, as the mesh size tends to zero, a uniformly convergent error expansion for the error density, with computable leading order term. It is shown that the optimal adaptive isotropic mesh uses a number of elements proportional to the d/2 power of the Ld/d+2 quasi-norm of the error density; the same error for approximation with a uniform mesh requires a number of elements proportional to the d/2 power of the larger L 1 norm of the same error density. A point is that this measure recognizes different convergence rates for multi scale problems, although the convergence order may be the same. The main result is a proof that the adaptive algorithm based on successive subdivisions of elements reduces the maximal error indicator with a factor or stops with the error asymptotically bounded by the tolerance using the optimal number of elements, up to a problem independent factor. An important step is to prove uniform convergence of the expansion for the error density, which is based on localized averages of second order difference quotients of the primal and dual finite element solutions. The averages are used since the difference quotients themselves do not converge pointwise for adapted meshes. The proof uses weak convergence techniques with a symmetrizer for the second order difference quotients and a splitting of the error into a dominating contribution, from elements with no hanging nodes or edges on the initial mesh, and a remaining asymptotically negligible part. Numerical experiments for an elasticity problem with a crack and different variants of the averages show that the algorithm is useful in practice also for relatively large tolerances, much larger than the small tolerances needed to theoretically guarantee that the algorithm works well.

Original languageEnglish (US)
Pages (from-to)367-407
Number of pages41
JournalBIT Numerical Mathematics
Volume46
Issue number2
DOIs
StatePublished - Jun 1 2006

Keywords

  • A posteriori error estimate
  • Adaptive methods
  • Computational complexity
  • Finite elements
  • Mesh refinement algorithm

ASJC Scopus subject areas

  • Software
  • Computer Networks and Communications
  • Computational Mathematics
  • Applied Mathematics

Cite this