Multilevel diagonal scaling preconditioners for boundary element equations on locally refined meshes

M. Ainsworth, W. McLean

Research output: Contribution to journalArticle

13 Citations (Scopus)

Abstract

We study a multilevel preconditioner for the Galerkin boundary element matrix arising from a symmetric positive-definite bilinear form. The associated energy norm is assumed to be equivalent to a Sobolev norm of positive, possibly fractional, order m on a bounded (open or closed) surface of dimension d, with $0<2mle d$. We consider piecewise linear approximation on triangular elements. Successive levels of the mesh are created by selectively subdividing elements within local refinement zones. Hanging nodes may be created and the global mesh ratio can grow exponentially with the number of levels. The coarse-grid correction consists of an exact solve, and the correction on each finer grid amounts to a simple diagonal scaling involving only those degrees of freedom whose associated nodal basis functions overlap the refinement zone. Under appropriate assumptions on the choice of refinement zones, the condition number of the preconditioned system is shown to be bounded by a constant independent of the number of degrees of freedom, the number of levels and the global mesh ratio. In addition to applying to Galerkin discretisation of hypersingular boundary integral equations, the theory covers finite element methods for positive-definite, self-adjoint elliptic problems with Dirichlet boundary conditions.
Original languageEnglish
Pages (from-to)387-413
Number of pages26
JournalNumerische Mathematik
Volume93
Issue number3
DOIs
Publication statusPublished - Jan 2003

Fingerprint

Boundary integral equations
Preconditioner
Boundary Elements
Boundary conditions
Mesh
Scaling
Finite element method
Galerkin
Positive definite
Refinement
Degree of freedom
Multilevel Preconditioners
Grid
Norm
Hypersingular Integral Equation
Piecewise Linear Approximation
Adjoint Problem
Local Refinement
Triangular Element
Bilinear form

Keywords

  • mathematics
  • diagonal scaling
  • boundary element equations
  • numerical analysis

Cite this

@article{26b5fe9c65aa4e52aae1bcf0eeaad6c1,
title = "Multilevel diagonal scaling preconditioners for boundary element equations on locally refined meshes",
abstract = "We study a multilevel preconditioner for the Galerkin boundary element matrix arising from a symmetric positive-definite bilinear form. The associated energy norm is assumed to be equivalent to a Sobolev norm of positive, possibly fractional, order m on a bounded (open or closed) surface of dimension d, with $0<2mle d$. We consider piecewise linear approximation on triangular elements. Successive levels of the mesh are created by selectively subdividing elements within local refinement zones. Hanging nodes may be created and the global mesh ratio can grow exponentially with the number of levels. The coarse-grid correction consists of an exact solve, and the correction on each finer grid amounts to a simple diagonal scaling involving only those degrees of freedom whose associated nodal basis functions overlap the refinement zone. Under appropriate assumptions on the choice of refinement zones, the condition number of the preconditioned system is shown to be bounded by a constant independent of the number of degrees of freedom, the number of levels and the global mesh ratio. In addition to applying to Galerkin discretisation of hypersingular boundary integral equations, the theory covers finite element methods for positive-definite, self-adjoint elliptic problems with Dirichlet boundary conditions.",
keywords = "mathematics, diagonal scaling, boundary element equations, numerical analysis",
author = "M. Ainsworth and W. McLean",
year = "2003",
month = "1",
doi = "10.1007/s002110100391",
language = "English",
volume = "93",
pages = "387--413",
journal = "Numerische Mathematik",
issn = "0029-599X",
number = "3",

}

Multilevel diagonal scaling preconditioners for boundary element equations on locally refined meshes. / Ainsworth, M.; McLean, W.

In: Numerische Mathematik, Vol. 93, No. 3, 01.2003, p. 387-413.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Multilevel diagonal scaling preconditioners for boundary element equations on locally refined meshes

AU - Ainsworth, M.

AU - McLean, W.

PY - 2003/1

Y1 - 2003/1

N2 - We study a multilevel preconditioner for the Galerkin boundary element matrix arising from a symmetric positive-definite bilinear form. The associated energy norm is assumed to be equivalent to a Sobolev norm of positive, possibly fractional, order m on a bounded (open or closed) surface of dimension d, with $0<2mle d$. We consider piecewise linear approximation on triangular elements. Successive levels of the mesh are created by selectively subdividing elements within local refinement zones. Hanging nodes may be created and the global mesh ratio can grow exponentially with the number of levels. The coarse-grid correction consists of an exact solve, and the correction on each finer grid amounts to a simple diagonal scaling involving only those degrees of freedom whose associated nodal basis functions overlap the refinement zone. Under appropriate assumptions on the choice of refinement zones, the condition number of the preconditioned system is shown to be bounded by a constant independent of the number of degrees of freedom, the number of levels and the global mesh ratio. In addition to applying to Galerkin discretisation of hypersingular boundary integral equations, the theory covers finite element methods for positive-definite, self-adjoint elliptic problems with Dirichlet boundary conditions.

AB - We study a multilevel preconditioner for the Galerkin boundary element matrix arising from a symmetric positive-definite bilinear form. The associated energy norm is assumed to be equivalent to a Sobolev norm of positive, possibly fractional, order m on a bounded (open or closed) surface of dimension d, with $0<2mle d$. We consider piecewise linear approximation on triangular elements. Successive levels of the mesh are created by selectively subdividing elements within local refinement zones. Hanging nodes may be created and the global mesh ratio can grow exponentially with the number of levels. The coarse-grid correction consists of an exact solve, and the correction on each finer grid amounts to a simple diagonal scaling involving only those degrees of freedom whose associated nodal basis functions overlap the refinement zone. Under appropriate assumptions on the choice of refinement zones, the condition number of the preconditioned system is shown to be bounded by a constant independent of the number of degrees of freedom, the number of levels and the global mesh ratio. In addition to applying to Galerkin discretisation of hypersingular boundary integral equations, the theory covers finite element methods for positive-definite, self-adjoint elliptic problems with Dirichlet boundary conditions.

KW - mathematics

KW - diagonal scaling

KW - boundary element equations

KW - numerical analysis

UR - http://dx.doi.org/10.1007/s002110100391

U2 - 10.1007/s002110100391

DO - 10.1007/s002110100391

M3 - Article

VL - 93

SP - 387

EP - 413

JO - Numerische Mathematik

JF - Numerische Mathematik

SN - 0029-599X

IS - 3

ER -