Source-linked AI summary
From multi-layered problems to multiple two-layered problems: a novel frequency-time hybrid multiple-scattering integral equation solver
Shuai Pan, Tao Yin, Lu Zhang
TL;DR
Accurate and efficient simulation of wave scattering in general multilayered media remains challenging because interface interactions and trapped-wave effects complicate direct broadband solvers. The paper decomposes the problem into N−1 two-layered subproblems, solves them with Fourier-transform and PML-BIE techniques, and assembles the resulting fields through multiple scattering. The analysis and numerical experiments show finite-time equivalence, high accuracy, and fast convergence of the resulting FTH-MS solver.
Problem
Direct frequency-domain treatment of general multilayered media remains challenging because possible trapped waves make frequency-robust broadband integral-equation solvers difficult.
Method
The method uses finite propagation speed to decompose an N-layered problem into N−1 two-layered subproblems, applying Fourier transforms, PML-BIE solvers, and multiple-scattering assembly.
Results
For fixed terminal time T, a finite number of multiple-scattering iterations provides a highly accurate approximation over the interval, with numerical tests verifying fast convergence and high accuracy.
Takeaways & Limitations
The FTH-MS formulation avoids constructing a global discretized system for the entire multilayered structure while providing an equivalent multiple-scattering representation supported by finite propagation speed.
Takeaways & Limitations
The presented formulation is two-dimensional, and low-frequency numerical evaluation remains challenging for broadband incident fields because of logarithmic frequency-domain singularities.
Abstract
from arXiv · showhide
This paper proposes a novel frequency-time hybrid multiple scattering (FTH-MS) integral equation solver for time-dependent wave equation problems in general multi-layered media, with remarkable scalability with respect to the number of layers $N$. In light of the finite speed of wave propagation, the new methodology provides an innovative multiple-scattering idea of re-modeling the original $N$-layered problem into a sequence of $N-1$ two-layered sub-problems, for which the main advantages lie in that (i) each sub-problem enjoys much simpler wave scattering properties compared with the complicated problem in a multi-layered medium, (ii) it enables to develop high-accuracy solver utilizing Fourier transform and frequency-domain boundary integral equation (BIE) method; and (iii) numerical evaluation of the sub-problems in each multiple scattering step can be parallelized. Both multiplicative- and additive-type strategies are developed and equivalence results, which indicate that the $M$-th order multiple scattering sums can provide equivalent representations of the solutions up to a certain time $T(M)$, are rigorously derived. Owing to the existed result of exponential convergence of the perfectly-matched-layer (PML) truncation for two-layered problem, all the sub-problems is numerically resolved by means of the FTH method based on the Fourier transform and the PML-BIE method whose numerical evaluation is addressed utilizing the Chebyshev-based rectangular-polar solver with high accuracy. Numerical examples are presented to validate the efficiency and accuracy of the proposed method.
1 Introduction
Accurate wave-scattering simulation in general multilayered media is difficult because repeated interface interactions create complex interference, while existing frequency- and time-domain approaches face truncation, conditioning, or integration challenges. The paper addresses this by decomposing the problem into two-layered subproblems and solving them with a frequency–time hybrid multiple-scattering method.
- Multilayered wave scattering is challenging because repeated reflections, transmissions, phase accumulation, interference, and resonance complicate accurate simulation.
- Frequency-domain volumetric methods require truncating unbounded domains, but rigorous transparent-boundary and PML treatments remain difficult for general multilayered interfaces.
- Time-domain boundary integral methods involve complicated time-dependent integration regions, while convolution quadrature inherits accuracy limits from inverse Z-transform contour integration.
- The proposed method uses finite propagation speed to represent an N-layered problem through N−1 two-layered scattering subproblems and reconstructs the global field from interface interactions.
- Fourier transformation, PML-BIE solves, and high-order inverse Fourier reconstruction yield the FTH-MS solver, with additive formulations suited to parallel computation.
2 Preliminaries
The paper formulates time-domain transmission problems in general N-layered media and develops domain-of-influence foundations for decomposing them into simpler subproblems. Finite propagation speed supports finite-order multiple-scattering approximations over finite time intervals, while Fourier transformation enables frequency-domain solution of each component.
- 2.1 Wave equation problem in a multi-layered medium: The model considers a two-dimensional general N-layered medium with N−1 locally perturbed interfaces, positive layer wave speeds, and separated adjacent interfaces.
- 2.1 Wave equation problem in a multi-layered medium: Incident fields are causal point sources or plane waves that reach the interfaces at time t0, after which reflected and transmitted fields satisfy wave equations and interface transmission conditions.
- 2.1 Wave equation problem in a multi-layered medium: The transmission problems are uniquely solvable for Lipschitz interfaces, providing the well-posedness basis for the subsequent scattering constructions.
- 2.3 FTH-MS solver for problem with multiple bounded obstacles: Direct broadband frequency-domain treatment of general multilayered media remains challenging because trapped waves may occur in intermediate layers.
- 2.2 Domain-of-influence property: The domain-of-influence property states that causal waves remain zero at a point until the incident wave and interlayer propagation paths can reach it.
- 2.3 FTH-MS solver for problem with multiple bounded obstacles: For well-separated bounded obstacles, the M-term multiple-scattering approximation equals the exact solution through time c−1Mδmin.
- 2.3 FTH-MS solver for problem with multiple bounded obstacles: Fourier transformation converts each time-domain subproblem into frequency-domain Helmholtz problems that can be solved with boundary integral methods and reconstructed from positive-frequency solutions.
3 Multiple-scattering strategy
The paper reformulates an N-layered wave problem as multiple-scattering combinations of two-layered subproblems. Multiplicative and additive constructions yield equivalent solution representations over finite time intervals, with additive steps enabling parallel evaluation.
- General strategy: The method replaces the general N-layered problem with a multiple-scattering series of N−1 two-layered medium subproblems.Each subproblem corresponds to one interface between adjacent media.
- Multiplicative-type strategy: For N = 3, multiplicative scattering models wave propagation as ping-pong interactions between the two interfaces.The resulting time-domain solution is formed by summing wave-equation solutions generated through repeated scattering.
- Equivalence results: The domain-of-influence property establishes that the M-th order multiplicative representation equals the original solution through a finite time T(M).This equivalence follows from finite propagation speed and recursively vanishing boundary data over the relevant time interval.
- Additive-type strategy: The additive strategy rewrites the multiplicative process so that the two two-layered subproblems can be solved in parallel at each scattering step.Their solutions are then summed in each layer to form the next-order approximation.
- Equivalence results: The additive construction likewise reproduces the three-layered solution on [0,T(M)] and extends to general N-layered media.The general result preserves the same finite-time equivalence structure across all layers.
- Generalization to N ≥3: The general N-layer formulation defines interface-wise boundary data and two-layered solutions recursively, then verifies the transmission conditions for the resulting series.The verification applies to both field and normal-derivative transmission conditions.
4 FTH-MS integral equation solver
The FTH-MS solver decomposes an N-layered wave problem into two-layered sub-problems, solves them in the frequency domain, and reconstructs time-domain solutions through Fourier transforms. Its implementation combines multiple-scattering iteration, FFT-based processing, and PML-BIE discretization, with low-frequency and long-time regimes identified as limitations.
- Solver construction: The method remodels a general N-layered problem into a union of two-layered sub-problems and extends the FTH solver to each sub-problem.This remodeling underlies the proposed FTH-MS integral equation solver.
- Solver construction: Each two-layered sub-problem is transformed into a frequency-domain scattering problem by applying Fourier transforms to the time-domain boundary data.The transformed data are defined inductively for scattering order m, with the first order using the original boundary data.
- Time reconstruction: The time-domain solution is recovered from the frequency-domain sub-problem solutions and agrees with the original solution on [0, T(M)] after M multiple-scattering orders.The numerical procedure requires frequency truncation, direct and inverse Fourier transforms, and a suitable BIE solver.
- Algorithm: Algorithm 1 iterates over multiple-scattering orders and interfaces, evaluates transformed boundary data, and solves each two-layered problem with a boundary integral equation solver.The algorithm initializes sub-problem fields, loops over m = 1, · · · , M and j = 1, · · · , N −1, then updates the solution.
- Fourier transform algorithm: For specified Gaussian pulses with ω0 − W > 0, the frequency interval [ω0 − W, ω0 + W] can be processed with FFT because the spectral amplitude is negligible outside it.The pulse is assumed to have exponentially decaying frequency content, allowing the interval to be selected so omitted solutions fall below machine precision.
- Limitations: Low-frequency problems remain challenging because PML-BIE accuracy is not uniformly high, while long-time windowing enlarges the required frequency range; both issues are left for future work.The paper notes that low-frequency inverse transforms can use modified Filon-Clenshaw-Curtis quadrature, but frequency-domain layered-media evaluation remains difficult.
- PML-BIE discretization: The PML-BIE method truncates the infinite interface and uses complex coordinate stretching, vanishing exterior boundary data, and PML-transformed layer potentials.Exponential convergence of PML truncation has been established for the two-layered medium problem.
5 Numerical experiments
Numerical experiments evaluate FTH and FTH-MS across two-, three-, and five-layered problems in two dimensions, plus a three-dimensional three-layered case. The results report fast convergence, high accuracy, finite-iteration approximation, and a low-frequency error-growth effect.
- Experimental setup: The experiments measure maximum time-domain error ε∞ over [0,T] at specified observation points, using exact or finely discretized reference solutions.All computations used a C++ implementation on a laptop with an AMD Ryzen 7 8845H processor and 32 GB of RAM.
- Two-layered problems: The two-layered 2D tests show fast convergence and high accuracy of the FTH and PML-BIE methods as the number of Fourier frequencies increases.Plane-wave and point-source cases use total fields, scattered-field traces, and maximum errors for T = 10.
- Three-layered problems: The three-layered 2D test shows fast frequency convergence and high accuracy, while finite multiple-scattering iterations approximate the solution accurately over a fixed terminal time.The reported equivalence result states that only finitely many iterations are required for highly accurate solutions for all t ∈ [0,T].
- Five-layered problems: The five-layered 2D test reports high accuracy for the FTH-MS solver and presents total fields at multiple times.The test uses a point source and examines the scattered-field trace and numerical error for different iteration counts.
- Low-frequency case: In the low-frequency three-layered 2D case, increasing numerical error behaves approximately as O(e^t) because the deformed Fourier path introduces exponential growth in the inverse transform.This behavior is attributed to the largest distance of the new integration path from zero.
- Three-dimensional problem: The three-dimensional three-layered example solves 51 equi-spaced frequencies over ω ∈ [1,19] with PML-BIE subproblems and uses M = 5 multiple-scattering iterations.The numerical solutions are shown at several times; further acceleration is left for future work.
6 Conclusion and further extensions
The paper develops FTH-MS by decomposing multilayered wave problems into two-layered sub-problems, avoiding a global discretized system. Further work targets stable long-time simulation, three-dimensional acceleration, and more complicated wave settings.
- FTH-MS decomposes complicated multilayered wave problems into sequences of two-layered problems, avoiding construction of a global discretized system.The approach combines multiple-scattering series with Fourier transform and PML-BIE methods.
- Further extensions: Stable long-time simulation for general incidences containing low-frequency information remains an extension of the proposed solver.
- Further extensions: Efficient acceleration of the solver for three-dimensional numerical evaluation remains an open extension, including fast multipole methods for PML-BIE.
- Further extensions: The method is also intended for more complicated elastic and electromagnetic waves and layered media with deep cavities.