Multigrid-Augmented Deep Learning Preconditioners for the Helmholtz Equation using Compact Implicit Layers
arXiv:2306.17486 · doi:10.1137/23M1583302
Abstract
We present a deep learning-based iterative approach to solve the discrete heterogeneous Helmholtz equation for high wavenumbers. Combining classical iterative multigrid solvers and convolutional neural networks (CNNs) via preconditioning, we obtain a learned neural solver that is faster and scales better than a standard multigrid solver. Our approach offers three main contributions over previous neural methods of this kind. First, we construct a multilevel U-Net-like encoder-solver CNN with an implicit layer on the coarsest grid of the U-Net, where convolution kernels are inverted. This alleviates the field of view problem in CNNs and allows better scalability. Second, we improve upon the previous CNN preconditioner in terms of the number of parameters, computation time, and convergence rates. Third, we propose a multiscale training approach that enables the network to scale to problems of previously unseen dimensions while still maintaining a reasonable training procedure. Our encoder-solver architecture can be used to generalize over different slowness models of various difficulties and is efficient at solving for many right-hand sides per slowness model. We demonstrate the benefits of our novel architecture with numerical experiments on a variety of heterogeneous two-dimensional problems at high wavenumbers.
arXiv commentPublished in SIAM Journal on Scientific Computing Copper Mountain Special Section on Multigrid Methods 2023
References in corpus (11)
- DeepONet: Learning nonlinear operators for identifying differential equations based on the universal approximation theorem of operators
- DGM: A deep learning algorithm for solving partial differential equations
- Solving high-dimensional partial differential equations using deep learning
- Hidden Physics Models: Machine Learning of Nonlinear Partial Differential Equations
- Learning data driven discretizations for partial differential equations
- Deep learning meets nanophotonics: A generalized accurate predictor for near fields and far fields of arbitrary 3D nanostructures
- A physics-informed operator regression framework for extracting data-driven continuum models
- A parallel sweeping preconditioner for heterogeneous 3D Helmholtz equations
- Int-Deep: A Deep Learning Initialized Iterative Method for Nonlinear Problems
- Efficient Inversion of Multiple-Scattering Model for Optical Diffraction Tomography
- A Helmholtz equation solver using unsupervised learning: Application to transcranial ultrasound