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.

Fast Fourier Solvers for the Tensor Product High-Order FEM for a Poisson Type Equation

A Novel Matrix-Free Finite Element Method for Time-Harmonic Maxwell’s Equations

Fast Poisson solvers for spectral methods

Fast Algorithms for Finite Element nonlinear Discrete Systems to Solve the PNP Equations

Non-equispaced Fourier Neural Solvers for PDEs

Fourier pseudospectral methods for the spatial variable-order fractional wave equations

A High Performance Approach for Solving the High Voltage Direct Current Ion Flow Field Problem by Tensor‐structured Finite Element Method

TetraFEM: Numerical Solution of Partial Differential Equations Using Tensor Train Finite Element Method

High order solution of Poisson problems with piecewise constant coefficients and interface jumps

A simple GPU implementation of spectral-element methods for solving 3D Poisson type equations on rectangular domains and its applications

FFT-Based High Order Central Difference Schemes for Poisson’s Equation with Staggered Boundaries

High-order Finite Element--Integral Equation Coupling on Embedded Meshes

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

A Fast and Accurate Solver for the Fractional Fokker-Planck Equation with Dirac-Delta Initial Conditions

Fast Algorithm for the Three-Dimensional Poisson Equation in Infinite Domains

Semi-analytical Exponential RKN Integrators for Efficiently Solving High-Dimensional Nonlinear Wave Equations Based on FFT Techniques.

Acceleration of Perturbation-Based Electric Field Integral Equations Using Fast Fourier Transform

Quasi-optimal complexity $hp$-FEM for Poisson on a rectangle

Tensor neural networks for high-dimensional Fokker-Planck equations

Fast finite difference solvers for singular solutions of the elliptic Monge-Ampère equation

Fast Sparse Eigen Solution In Finite Element Analysis