Volume 6, pp. 133-152, 1997.

Multigrid method for $H({\rm div})$ in three dimensions

R. Hiptmair


We are concerned with the design and analysis of a multigrid algorithm for ${\bf H}(\mbox{div};\Omega)$–elliptic linear variational problems. The discretization is based on ${\bf H}(\mbox{div};\Omega)$–conforming Raviart–Thomas elements. A thorough examination of the relevant bilinear form reveals that a separate treatment of vector fields in the kernel of the divergence operator and its complement is paramount. We exploit the representation of discrete solenoidal vector fields as $\mbox{curl}$s of finite element functions in so-called Nédélec spaces. It turns out that a combined nodal multilevel decomposition of both the Raviart–Thomas and Nédélec finite element spaces provides the foundation for a viable multigrid method. Its Gauß–Seidel smoother involves an extra stage where solenoidal error components are tackled. By means of elaborate duality techniques we can show the asymptotic optimality in the case of uniform refinement. Numerical experiments confirm that the typical multigrid efficiency is actually achieved for model problems.

Full Text (PDF) [214 KB], BibTeX

Key words

multigrid, Raviart–Thomas finite elements, Nédélec's finite elements, multilevel, mixed finite elements.

AMS subject classifications

65N55, 65N30.

ETNA articles which cite this article

Vol. 11 (2000), pp. 1-24 Andrea Toselli: Neumann-Neumann methods for vector field problems
Vol. 15 (2003), pp. 132-151 B. Chang and B. Lee: A multigrid algorithm for solving the multi-group, anisotropic scattering Boltzmann equation using first-order system least-squares methodology

< Back