A heterogeneous FMM for layered media Helmholtz equation. I: Two layers in R^2
From MaRDI portal
Publication:725487
DOI10.1016/J.JCP.2018.05.007zbMATH Open1392.78018arXiv1703.09136OpenAlexW2799865484MaRDI QIDQ725487FDOQ725487
Authors: Jingfang Huang, Dangxing Chen, Wei Cai, Min-Hyung Cho 
Publication date: 1 August 2018
Published in: Journal of Computational Physics (Search for Journal in Brave)
Abstract: In this paper, we will introduce a new heterogeneous fast multipole method (H-FMM) for 2-D Helmholtz equation in layered media. To illustrate the main algorithm ideas, we focus on the case of two and three layers in this work. The key compression step in the H-FMM is based on a fact that the multipole expansion for the sources of the free-space Green's function can be used also to compress the far field of the sources of the layered-media or domain Green's function, and a similar result exists for the translation operators for the multipole and local expansions. The mathematical error analysis is shown rigorously by an image representation of the Sommerfeld spectral form of the domain Green's function. As a result, in the H-FMM algorithm, both the "multipole-to-multipole" and "local-to-local" translation operators are the same as those in the free-space case, allowing easy adaptation of existing free-space FMM. All the spatially variant information of the domain Green's function are collected into the "multipole-to-local" translations and therefore the FMM becomes "heterogeneous". The compressed representation further reduces the cost of evaluating the domain Green's function when computing the local direct interactions. Preliminary numerical experiments are presented to demonstrate the efficiency and accuracy of the algorithm with much improved performance over some existing methods for inhomogeneous media. Furthermore, we also show that, due to the equivalence between the complex line image representation and Sommerfeld integral representation of layered media Green's function, the new algorithm can be generalized to multi-layered media with minor modification where details for compression formulas, translation operators, and bookkeeping strategies will be addressed in a subsequent paper.
Full work available at URL: https://arxiv.org/abs/1703.09136
Recommendations
hierarchical modelHelmholtz equationfast multipole methodimpedance boundary conditionlow-rank representationmultilayered media
Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Composite media; random media in optics and electromagnetic theory (78A48) Multipole methods applied to problems in optics and electromagnetic theory (78M16)
Cites Work
- A wideband fast multipole method for the two-dimensional complex Helmholtz equation
- DASHMM: dynamic adaptive system for hierarchical multipole methods
- Title not available (Why is that?)
- A sparse matrix arithmetic based on \({\mathfrak H}\)-matrices. I: Introduction to \({\mathfrak H}\)-matrices
- A sparse \({\mathcal H}\)-matrix arithmetic. II: Application to multi-dimensional problems
- An Algorithm for the Machine Calculation of Complex Fourier Series
- Integral equation methods in scattering theory
- Multi-Level Adaptive Solutions to Boundary-Value Problems
- Hybrid Gauss-Trapezoidal Quadrature Rules
- Rapid solution of integral equations of scattering theory in two dimensions
- High-Order Corrected Trapezoidal Quadrature Rules for Singular Functions
- Linear integral equations
- A wideband fast multipole method for the Helmholtz equation in three dimensions
- Quadrature by expansion: a new method for the evaluation of layer potentials
- Fast direct solvers for integral equations in complex three-dimensional domains
- A parallel fast algorithm for computing the Helmholtz integral operator in 3-D layered media
- A fast direct solver for structured linear systems by recursive skeletonization
- Preface
- Extending the fast multipole method to charges inside or outside a dielectric sphere
- On the efficient representation of the half-space impedance Green's function for the Helmholtz equation
- A fast volume integral equation solver for electromagnetic scattering from large inhomogeneous objects in planarly layered media
- Efficient and accurate computation of electric field dyadic Green's function in layered media
- Windowed Green function method for layered-media scattering
- Accurate and efficient Nyström volume integral equation method for electromagnetic scattering of 3-D metamaterials in layered media
- Computational software: simple FMM libraries for electrostatics, slow viscous flow, and frequency-domain wave propagation
- RECFMM: recursive parallelization of the adaptive fast multipole method for Coulomb and screened Coulomb interactions
Cited In (10)
- Spectrally-accurate numerical method for acoustic scattering from doubly-periodic 3D multilayered media
- Uniform far-field asymptotics of the two-layered Green function in two dimensions and application to wave scattering in a two-layered medium
- Exponential convergence for multipole and local expansions and their translations for sources in layered media: two-dimensional acoustic wave
- An \(O(N \log N)\) hierarchical random compression method for kernel matrices by sampling partial matrix entries
- A hybrid stochastic interpolation and compression method for kernel matrices
- Adapting free-space fast multipole method for layered media Green's function: algorithm and analysis
- A new numerical approach to the solution of the 2-D Helmholtz equation with optimal accuracy on irregular domains and Cartesian meshes
- On the construction of scattering matrices for irregular or elongated enclosures using Green's representation formula
- Fast multipole method for 3-D Laplace equation in layered media
- A parallel fast algorithm for computing the Helmholtz integral operator in 3-D layered media
Uses Software
This page was built for publication: A heterogeneous FMM for layered media Helmholtz equation. I: Two layers in \(\mathbb{R}^2\)
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q725487)
