Header

UZH-Logo

Maintenance Infos

Hybrid Galerkin boundary elements: theory and implementation


Graham, I; Hackbusch, W; Sauter, S (2000). Hybrid Galerkin boundary elements: theory and implementation. Numerische Mathematik, 86(1):139-172.

Abstract

In this paper we present a new quadrature method for computing Galerkin stiffness matrices arising from the discretisation of 3D boundary integral equations using continuous piecewise linear boundary elements. This rule takes as points some subset of the nodes of the mesh and can be used for computing non-singular Galerkin integrals corresponding to pairs of basis functions with non-intersecting supports. When this new rule is combined with standard methods for the singular Galerkin integrals we obtain a `hybrid' Galerkin method which has the same stability and asymptotic convergence properties as the true Galerkin method but a complexity more akin to that of a collocation or Nyström method. The method can be applied to a wide range of singular and weakly-singular first- and second-kind equations, including many for which the classical Nyström method is not even defined. The results apply to equations on piecewise-smooth Lipschitz boundaries, and to non-quasiuniform (but shape-regular) meshes. A by-product of the analysis is a stability theory for quadrature rules of precision 1 and 2 based on arbitrary points in the plane. Numerical experiments demonstrate that the new method realises the performance expected from the theory.

Abstract

In this paper we present a new quadrature method for computing Galerkin stiffness matrices arising from the discretisation of 3D boundary integral equations using continuous piecewise linear boundary elements. This rule takes as points some subset of the nodes of the mesh and can be used for computing non-singular Galerkin integrals corresponding to pairs of basis functions with non-intersecting supports. When this new rule is combined with standard methods for the singular Galerkin integrals we obtain a `hybrid' Galerkin method which has the same stability and asymptotic convergence properties as the true Galerkin method but a complexity more akin to that of a collocation or Nyström method. The method can be applied to a wide range of singular and weakly-singular first- and second-kind equations, including many for which the classical Nyström method is not even defined. The results apply to equations on piecewise-smooth Lipschitz boundaries, and to non-quasiuniform (but shape-regular) meshes. A by-product of the analysis is a stability theory for quadrature rules of precision 1 and 2 based on arbitrary points in the plane. Numerical experiments demonstrate that the new method realises the performance expected from the theory.

Statistics

Citations

8 citations in Web of Science®
7 citations in Scopus®
Google Scholar™

Altmetrics

Additional indexing

Item Type:Journal Article, refereed, original work
Communities & Collections:07 Faculty of Science > Institute of Mathematics
Dewey Decimal Classification:510 Mathematics
Language:English
Date:2000
Deposited On:29 Nov 2010 16:27
Last Modified:06 Dec 2017 20:54
Publisher:Springer
ISSN:0029-599X
Publisher DOI:https://doi.org/10.1007/PL00005400
Related URLs:http://www.zentralblatt-math.org/zbmath/search/?q=an%3A0966.65091

Download

Full text not available from this repository.
View at publisher