Abstract:In this paper, we present a multiscale framework for solving the Helmholtz equation in heterogeneous media without scale separation and in the high frequency regime where the wavenumber $k$ can be large. The main innovation is that our methods achieve a nearly exponential rate of convergence with respect to the computational degrees of freedom, using a coarse grid of mesh size $O(1/k)$ without suffering from the well-known pollution effect. The key idea is a non-overlapped domain decomposition and its associated coarse-fine scale decomposition of the solution space that adapts to the media property and wavenumber; this decomposition is inspired by the multiscale finite element method (MsFEM). We show that the coarse part is of \textit{low complexity} in the sense that it can be approximated with a nearly exponential rate of convergence via local basis functions, due to the compactness of a restriction operator that maps Helmholtz-harmonic functions to their interpolation residues on edges, while the fine part is \textit{local} such that it can be computed efficiently using the local information of the right hand side. The combination of the two parts yields the overall nearly exponential rate of convergence of our multiscale method. Our method draws many connections to multiscale methods in the literature, which we will comment in detail. We demonstrate the effectiveness of our methods theoretically and numerically; an exponential rate of convergence is consistently observed and confirmed. In addition, we observe the robustness of our methods regarding the high contrast in the media numerically. We specifically focus on 2D problems in our exposition since the geometry of non-overlapped domain decomposition is simplest to explain in such cases; generalizations to 3D will be outlined at the end.

A Neural Multigrid Solver for Helmholtz Equations with High Wavenumber and Heterogeneous Media

Multigrid-Augmented Deep Learning Preconditioners for the Helmholtz Equation using Compact Implicit Layers

A matrix-free parallel solution method for the three-dimensional heterogeneous Helmholtz equation

Multigrid-augmented deep learning preconditioners for the Helmholtz equation

Stand-alone Multigrid for Helmholtz Revisited: Towards Convergence Using Standard Components

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

Adaptive complex frequency with V-cycle GMRES for preconditioning 3D Helmholtz equation

An accurate and efficient multiscale finite-difference frequency-domain method for the scalar Helmholtz equation

Multigrid method applied to two kinds of Helmholtz equations

Deep Neural Helmholtz Operators for 3D Elastic Wave Propagation and Inversion

Exponentially convergent multiscale methods for high frequency heterogeneous Helmholtz equations

UGrid: An Efficient-And-Rigorous Neural Multigrid Solver for Linear PDEs

DNN-MG: A Hybrid Neural Network/Finite Element Method with Applications to 3D Simulations of the Navier-Stokes Equations

Robust Multiscale Methods for Helmholtz equations in high contrast heterogeneous media

A Robust and Scalable Multigrid Solver for 3-D Low-Frequency Electromagnetic Diffusion Problems

A deep learning algorithm to accelerate algebraic multigrid methods in finite element solvers of 3D elliptic PDEs

A scalable multidimensional fully implicit solver for Hall magnetohydrodynamics

Numerical Solution of the Helmholtz Equation with High Wavenumbers