On the evaluation of layer potentials close to their sources
Författare
Summary, in English
When solving elliptic boundary value problems using integral
equation methods one may need to evaluate potentials represented by
a convolution of discretized layer density sources against a kernel.
Standard quadrature accelerated with a fast hierarchical method for
potential field evaluation gives accurate results far away from the
sources. Close to the sources this is not so. Cancellation and
nearly singular kernels may cause serious degradation. This paper
presents a new scheme based on a mix of composite polynomial
quadrature, layer density interpolation, kernel approximation,
rational quadrature, high polynomial order corrected interpolation
and differentiation, temporary panel mergers and splits, and a
particular implementation of the GMRES solver. Criteria for which
mix is fastest and most accurate in various situations are also
supplied. The paper focuses on the solution of the Dirichlet problem
for Laplace's equation in the plane. In a series of examples we
demonstrate the efficiency of the new scheme for interior domains
and domains exterior to up to 2000 close-to-touching contours.
Densities are computed and potentials are evaluated, rapidly and
accurate to almost machine precision, at points that lie arbitrarily
close to the boundaries.
equation methods one may need to evaluate potentials represented by
a convolution of discretized layer density sources against a kernel.
Standard quadrature accelerated with a fast hierarchical method for
potential field evaluation gives accurate results far away from the
sources. Close to the sources this is not so. Cancellation and
nearly singular kernels may cause serious degradation. This paper
presents a new scheme based on a mix of composite polynomial
quadrature, layer density interpolation, kernel approximation,
rational quadrature, high polynomial order corrected interpolation
and differentiation, temporary panel mergers and splits, and a
particular implementation of the GMRES solver. Criteria for which
mix is fastest and most accurate in various situations are also
supplied. The paper focuses on the solution of the Dirichlet problem
for Laplace's equation in the plane. In a series of examples we
demonstrate the efficiency of the new scheme for interior domains
and domains exterior to up to 2000 close-to-touching contours.
Densities are computed and potentials are evaluated, rapidly and
accurate to almost machine precision, at points that lie arbitrarily
close to the boundaries.
Avdelning/ar
- Matematik LTH
- Harmonic Analysis and Applications
Publiceringsår
2008
Språk
Engelska
Sidor
2899-2921
Publikation/Tidskrift/Serie
Journal of Computational Physics
Volym
227
Issue
5
Fulltext
Länkar
Dokumenttyp
Artikel i tidskrift
Förlag
Elsevier
Ämne
- Mathematics
Status
Published
Forskningsgrupp
- Harmonic Analysis and Applications
- Harmonic Analysis and Applications
ISBN/ISSN/Övrigt
- ISSN: 0021-9991