Stability and convergence of the Method of Fundamental Solutions for Helmholtz problems on analytic domains
arXiv:0708.3533 · doi:10.1016/j.jcp.2008.04.008
Abstract
The Method of Fundamental Solutions (MFS) is a popular tool to solve Laplace and Helmholtz boundary value problems. Its main drawback is that it often leads to ill-conditioned systems of equations. In this paper we investigate for the interior Helmholtz problem on analytic domains how the singularities (charge points) of the MFS basis functions have to be chosen such that approximate solutions can be represented by the MFS basis in a numerically stable way. For Helmholtz problems on the unit disc we give a full analysis which includes the high frequency (short wavelength) limit. For more difficult and nonconvex domains such as crescents we demonstrate how the right choice of charge points is connected to how far into the complex plane the solution of the boundary value problem can be analytically continued, which in turn depends on both domain shape and boundary data. Using this we develop a recipe for locating charge points which allows us to reach error norms of typically 10^{-11} on a wide variety of analytic domains. At high frequencies of order only 3 points per wavelength are needed, which compares very favorably to boundary integral methods.
arXiv comment25 pages, 14 figures, submitted to J. Comp. Phys
References in corpus (1)
Cited by in corpus (26)
- The AAA algorithm for rational approximation
- Trefftz discontinuous Galerkin methods on unstructured meshes for the wave equation
- A Survey of Trefftz Methods for the Helmholtz Equation
- A new integral representation for quasiperiodic fields and its application to two-dimensional band structure calculations
- Efficient numerical solution of acoustic scattering from doubly-periodic arrays of axisymmetric objects
- New Laplace and Helmholtz solvers
- Robust fast direct integral equation solver for quasi-periodic scattering problems with a large number of layers
- Flexibly imposing periodicity in kernel independent FMM: A Multipole-To-Local operator approach
- A Multiphase Shape Optimization Problem for Eigenvalues: Qualitative Study and Numerical Results
- General-purpose kernel regularization of boundary integral equations via density interpolation
- A robust solver for elliptic PDEs in 3D complex geometries
- A hybrid approach to solve the high-frequency Helmholtz equation with source singularity in smooth heterogeneous media
- Stable approximation of Helmholtz solutions in the disk by evanescent plane waves
- Comparable upper and lower bounds for boundary values of Neumann eigenfunctions and tight inclusion of eigenvalues
- A Time-Domain Method of Auxiliary Sources for Efficient Analysis of Transient Electromagnetic Scattering by Moderately Conductive Cylinders
- Evaluation of Abramowitz functions in the right half of the complex plane
- A unified integral equation scheme for doubly-periodic Laplace and Stokes boundary value problems in two dimensions
- Minimization of the eigenvalues of the dirichlet-laplacian with a diameter constraint
- A fast direct solver for two dimensional quasi-periodic multilayered media scattering problems
- On polynomial interpolation in the monomial basis
- A well conditioned Method of Fundamental Solutions
- Numerically stable computation of embedding formulae for scattering by polygons
- Preemptive Termination of Suggestions during Sequential Kriging Optimization of a Brain Activity Reconstruction Simulation
- Method of fundamental solutions for Neumann problems of the modified Helmholtz equation in disk domains
- Evaluation of layer potentials close to the boundary for Laplace and Helmholtz problems on analytic planar domains
- On the use of rational-function fitting methods for the solution of 2D Laplace boundary-value problems