Abstract:The recent results presented in <a class="link-https" data-arxiv-id="2202.05608" href="https://arxiv.org/abs/2202.05608">arXiv:2202.05608</a> have led to significant developments in achieving stable approximations of Helmholtz solutions by plane wave superposition. The study shows that the numerical instability and ill-conditioning inherent in plane wave-based Trefftz methods can be effectively overcome with regularization techniques, provided there exist accurate approximations in the form of expansions with bounded coefficients. Whenever the target solution contains high Fourier modes, propagative plane waves fail to yield stable approximations due to the exponential growth of the expansion coefficients. Conversely, evanescent plane waves, whose modal content covers high Fourier regimes, are able to provide both accurate and stable results. The developed numerical approach, which involves constructing evanescent plane wave approximation sets by sampling the parametric domain according to a probability density function, results in substantial improvements when compared to conventional propagative plane wave schemes. The following work extends this research to the three-dimensional setting, confirming the achieved results and introducing new ones. By generalizing the 3D Jacobi$-$Anger identity to complex-valued directions, we show that any Helmholtz solution in a ball can be represented as a continuous superposition of evanescent plane waves. This representation extends the classical Herglotz one and provides a relevant stability result that cannot be achieved with the use of propagative waves alone. The proposed numerical recipes have been tailored for the 3D setting and extended with new sampling strategies involving extremal systems of points. These methods are tested by numerical experiments, showing the desired accuracy and bounded-coefficient stability, in line with the two-dimensional case.

A variant of the plane wave least squares method for the time-harmonic Maxwell's equations

Plane Wave Discontinuous Galerkin methods for the Helmholtz equation and Maxwell equations in Anisotropic Media

A geometric optics ansatz-based plane wave method for two dimensional Helmholtz equations with variable wave numbers

A plane wave method based on approximate wave directions for two dimensional Helmholtz equations with large wave numbers

A discontinuous plane wave neural network method for Helmholtz equation and time-harmonic Maxwell's equations

A Novel Least Squares Method for Helmholtz Equations with Large Wave Numbers

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

Stable approximation of Helmholtz solutions in the ball using evanescent plane waves

Stable approximation of Helmholtz solutions in the disk by evanescent plane waves

A virtual element method for the solution of 2D time-harmonic elastic wave equations via scalar potentials

Evanescent Plane Wave Approximation of Helmholtz Solutions in Spherical Domains

A Discontinuous Least Squares Finite Element Method for the Helmholtz Equation

New discretization schemes for time-harmonic Maxwell equations by weak Galerkin finite element methods

An exact in time Fourier pseudospectral method with multiple conservation laws for three-dimensional Maxwell’s equations

An Investigation Of Wave Propagation With High Wave Numbers Via The Regularized Lbiem

A Novel Hybrid Implicit Explicit - Pseudospectral Time Domain Method for TMz Waves

Babich-Like Ansatz for Three-Dimensional Point-Source Maxwell's Equations in an Inhomogeneous Medium at High Frequencies.

Boussinesq-type Modeling in Surf Zone Using Mesh-Less Least-Square-based Finite Difference Method

A Modified weak Galerkin finite element method for the Maxwell equations on polyhedral meshes

Solving a Fully Nonlinear Highly Dispersive Boussinesq Model with Mesh‐less Least Square‐based Finite Difference Method

Numerical investigation of a 3D hybrid high-order method for the indefinite time-harmonic Maxwell problem