Fast integral equation methods for the modified Helmholtz equation

  • Authors:

  • Mary Catherine A. Kropinski;Bryan D. Quaife

  • Affiliations:

  • Department of Mathematics, Simon Fraser University, Burnaby, British Columbia, Canada V5A 1S6;Department of Mathematics, Simon Fraser University, Burnaby, British Columbia, Canada V5A 1S6

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

Quantified Score

Hi-index 31.45

Visualization

Abstract

We present integral equation methods for the solution to the two-dimensional, modified Helmholtz equation, u(x)-@a^2@Du(x)=0, in bounded or unbounded multiply-connected domains. We consider both Dirichlet and Neumann problems. We derive well-conditioned Fredholm integral equations of the second kind, which are discretized using high-order, hybrid Gauss-trapezoid rules. Our fast multipole-based iterative solution procedure requires only O(N) operations, where N is the number of nodes in the discretization of the boundary. We demonstrate the performance of our methods on several numerical examples, and we show that they have both the ability to handle highly complex geometry and the potential to solve large-scale problems.