FFT-accelerated computation of the Dirichlet-to-Neumann map for inhomogeneous exterior Helmholtz problems using the method of fundamental solutions

TL;DR

This paper introduces an efficient FFT-accelerated method using the method of fundamental solutions to compute the Dirichlet-to-Neumann map for inhomogeneous exterior Helmholtz problems, enabling accurate solutions on unbounded domains.

Contribution

It combines MFS and FFT to efficiently compute the DtN map, reducing the exterior problem to a bounded domain for finite element analysis.

Findings

High-accuracy computation of the DtN map using MFS and FFT.

Finite element solutions incorporating the DtN map are highly accurate.

The method efficiently handles inhomogeneous terms in exterior Helmholtz problems.

Abstract

An efficient numerical method is proposed for computing the Dirichlet-to-Neumann (DtN) map associated with the exterior Dirichlet problem for the two-dimensional Helmholtz equation with an inhomogeneous term. The exterior solution is approximated by the method of fundamental solutions (MFS). When the source and collocation points are equally spaced on concentric circles, the coefficient matrices arising in the discretization become circulant, which enables efficient evaluation of the discrete DtN map by the fast Fourier transform (FFT). By incorporating the boundary condition defined by the DtN map into a finite element formulation, the original exterior problem posed on an unbounded domain is reduced to an equivalent problem on a bounded computational domain, which can then be solved by the finite element method (FEM). Numerical examples show that the proposed MFS-based approach computes the DtN map with high accuracy, and that the FEM incorporating the boundary condition defined by the DtN map yields accurate solutions for exterior Helmholtz problems with an inhomogeneous term posed on unbounded domains.