Abstract:This article presents novel numerical algorithms based on pseudodifferential operators for fast, direct, solution of the Helmholtz equation in 1D, 2D, and 3D inhomogeneous unbounded media. The proposed approach relies on an Operator Fourier Transform (OFT) representation of pseudodifferential operators ({\Psi}DO) which frame the problem of computing the inverse Helmholtz operator with a spatially-dependent wave speed in terms of two sequential applications of an inverse square root {\Psi}DO. The OFT representation of the action of the inverse square root {\Psi}DO, in turn, can be effected as a superposition of solutions of a pseudo-temporal initial-boundary-value problem for a paraxial equation. The OFT framework offers several advantages over traditional direct and iterative approaches for the solution of the Helmholtz equation. The operator integral transform is amenable to standard quadrature methods and the required pseudo-temporal paraxial equation solutions can be obtained using any suitable numerical method. A specialized quadrature is derived to evaluate the OFT efficiently and an alternating direction implicit method, used in conjunction with standard finite differences, is used to solve the requisite component paraxial equation problems. Numerical studies in 1D, 2D, and 3D are presented to confirm the expected OFT-based Helmholtz solver convergence rate. In addition, the efficiency and versatility of our proposed approach is demonstrated by tackling nontrivial wave propagation problems, including two-dimensional plane wave scattering from a geometrically complex inhomogeneity, three-dimensional scattering from turbulent channel flow and plane wave transmission through a spherically-symmetric gradient-index Luneburg lens. All computations, 3D problems which involve solving the Helmholtz equation with more than one billion complex unknowns, are performed in a single workstation.

Partial balayage for the Helmholtz equation

Multi-phase $k$-quadrature domains and applications to acoustic waves and magnetic fields

High-Order Nyström/Convolution-Quadrature Solution of Time-Domain Scattering from Closed and Open Lipschitz Boundaries with Dirichlet and Neumann Boundary Conditions

Method of virtual sources using on-surface radiation conditions for the Helmholtz equation

On the theory of balayage on locally compact spaces and its applications

Integral representation of balayage on locally compact spaces and its application

A multi-domain spectral-Galerkin method for the Neumann problem on quadrilaterals

High-Order Quadratures for the Solution of Scattering Problems in Two Dimensions

A High-Order Well-Balanced Discontinuous Galerkin Method for Hyperbolic Balance Laws Based on the Gauss-Lobatto Quadrature Rules

An iterative method for Helmholtz boundary value problems arising in wave propagation

An Additive Overlapping Domain Decomposition Method for the Helmholtz Equation

On the theory of pseudo-balayage on a locally compact space and its applications

Balayage of Radon measures of infinite energy on locally compact spaces

Convolution quadrature for the wave equation with a nonlinear impedance boundary condition

A novel direct Helmholtz solver in inhomogeneous media based on the operator Fourier transform functional calculus

Spectral solver for Cauchy problems in polar coordinates using discrete Hankel transforms

Hybridizable Discontinuous Galerkin Methods for Helmholtz Equation with High Wave Number. Part I: Linear case

Non-homogeneous initial boundary value problems and small data scattering of the Hirota equations posed on the half line

The Calderón problem for the Schrödinger equation in transversally anisotropic geometries with partial data

Approximating partial differential equations without boundary conditions

Extended Smoothed Boundary Method for Solving Partial Differential Equations with General Boundary Conditions on Complex Boundaries