A Nyström method for weakly singular integral operators on surfaces

  • Authors:
  • James Bremer;Zydrunas Gimbutas

  • Affiliations:
  • Department of Mathematics, University of California, Davis, CA 95616, USA;Courant Institute of Mathematical Sciences, New York University, New York, NY 10012, USA

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2012

Quantified Score

Hi-index 31.45

Visualization

Abstract

We describe a modified Nystrom method for the discretization of the weakly singular boundary integral operators which arise from the formulation of linear elliptic boundary value problems as integral equations. Standard Nystrom and collocation schemes proceed by representing functions via their values at a collection of quadrature nodes. Our method uses appropriately scaled function values in lieu of such representations. This results in a scheme which is mathematically equivalent to Galerkin discretization in that the resulting matrices are related to those obtained by Galerkin methods via conjugation with well-conditioned matrices, but which avoids the evaluation of double integrals. Moreover, we incorporate a new mechanism for approximating the singular integrals which arise from the discretization of weakly singular integral operators which is considerably more efficient than standard methods. We illustrate the performance of our method with numerical experiments.