Fast, higher-order direct/iterative hybrid solver for scattering by Inhomogeneous media - with application to high-frequency and discontinuous refractivity problems

This paper presents a fast high-order method for the solution of two-dimensional problems of scattering by penetrable inhomogeneous media, with application to high-frequency configurations containing (possibly) discontinuous refractivities. The method relies on a combination of a differential volumetric formulation and a boundary integral formulation. Thus, in the proposed approach the entire computational domain is partitioned into large numbers of volumetric spectral approximation patches which are then grouped into patch subsets for local direct solution; the interactions with the exterior domain are handled by means of a boundary integral equation. The resulting algorithm can be quite effective: after a modestly-demanding precomputation stage (whose results for a given frequency can be repeatedly used for arbitrarily chosen incidence angles), the proposed algorithm can accurately evaluate scattering by configurations including large and complex objects and/or high refractivity contrasts, including possibly refractive-index discontinuities, in fast single-core runs.

[1]  Per-Gunnar Martinsson,et al.  A spectrally accurate direct solution technique for frequency-domain scattering problems with variable media , 2013, 1308.5998.

[2]  R. Kress,et al.  Inverse Acoustic and Electromagnetic Scattering Theory , 1992 .

[3]  G. Vainikko.  Fast Solvers of the Lippmann-Schwinger Equation , 2000 .

[4]  G. Schmidt,et al.  Numerical Solution of the Lippmann–Schwinger Equation by Approximate Approximations , 2004 .

[5]  Lexing Ying,et al.  Sparsify and Sweep: An Efficient Preconditioner for the Lippmann-Schwinger Equation , 2017, SIAM J. Sci. Comput..

[6]  Timothy A. Davis,et al.  Algorithm 832: UMFPACK V4.3---an unsymmetric-pattern multifrontal method , 2004, TOMS.

[7]  Ran Duan,et al.  High-order quadratures for the solution of scattering problems in two dimensions , 2009, J. Comput. Phys..

[8]  Thomas Hagstrom,et al.  A new auxiliary variable formulation of high-order local radiation boundary conditions: corner compatibility conditions and extensions to first-order systems , 2004 .

[9]  J. Waldvogel.  Fast Construction of the Fejér and Clenshaw–Curtis Quadrature Rules , 2006 .

[10]  Alfredo Bermúdez,et al.  An optimal perfectly matched layer with unbounded absorbing function for time-harmonic acoustic scattering problems , 2007, J. Comput. Phys..

[11]  A. Bayliss,et al.  On accuracy conditions for the numerical computation of waves , 1985 .

[12]  Akash Anand,et al.  An efficient high-order Nyström scheme for acoustic scattering by inhomogeneous penetrable media with discontinuous material interface , 2015, J. Comput. Phys..

[13]  Leslie Greengard,et al.  Fast, Adaptive, High-Order Accurate Discretization of the Lippmann-Schwinger Equation in Two Dimensions , 2015, SIAM J. Sci. Comput..

[14]  Eli Turkel,et al.  THE NUMERICAL SOLUTION OF THE HELMHOLTZ EQUATION FOR WAVE PROPAGATION PROBLEMS IN UNDERWATER ACOUSTICS , 1985 .

[15]  A. Majda,et al.  Absorbing boundary conditions for the numerical simulation of waves , 1977 .

[16]  Akash Anand,et al.  Improved convergence of fast integral equation solvers for acoustic scattering by inhomogeneous penetrable media with discontinuous material interface , 2018, J. Comput. Phys..

[17]  Oscar P. Bruno,et al.  A fast, higher-order solver for scattering by penetrable bodies in three dimensions , 2005 .

[18]  Oscar P. Bruno,et al.  A Chebyshev-based rectangular-polar integral solver for scattering by geometries described by non-overlapping patches , 2018, J. Comput. Phys..

[19]  Stefan A. Sauter,et al.  Is the Pollution Effect of the FEM Avoidable for the Helmholtz Equation Considering High Wave Numbers? , 1997, SIAM Rev..

[20]  B. Engquist,et al.  Sweeping preconditioner for the Helmholtz equation: Hierarchical matrix representation , 2010, 1007.4290.

[21]  Oscar P. Bruno,et al.  An efficient, preconditioned, high-order solver for scattering by two-dimensional inhomogeneous media , 2004 .

[22]  O. Bruno,et al.  A fast, high-order algorithm for the solution of surface scattering problems: basic implementation, tests, and applications , 2001 .

[23]  Clifford H. Thurber,et al.  Parameter estimation and inverse problems , 2005 .

[24]  Fredrik Andersson,et al.  A Fast, Bandlimited Solver for Scattering Problems in Inhomogeneous Media , 2005 .

[25]  Yu Chen,et al.  A high-order, fast algorithm for scattering calculation in two dimensions , 2004 .

[26]  Peter Kuchment,et al.  The Radon Transform and Medical Imaging , 2014, CBMS-NSF regional conference series in applied mathematics.

[27]  Lexing Ying,et al.  Sparsifying Preconditioner for the Lippmann-Schwinger Equation , 2014, Multiscale Model. Simul..

[28]  Oscar P. Bruno,et al.  Higher-Order Fourier Approximation in Scattering by Two-Dimensional, Inhomogeneous Media , 2004, SIAM J. Numer. Anal..

[29]  Hongkai Zhao,et al.  Fast Alternating BiDirectional Preconditioner for the 2D High-Frequency Lippmann-Schwinger Equation , 2016, SIAM J. Sci. Comput..

[30]  Marcus J. Grote,et al.  On nonreflecting boundary conditions , 1995 .

[31]  Martin J. Gander,et al.  Why it is Difficult to Solve Helmholtz Problems with Classical Iterative Methods , 2012 .

[32]  Peter Monk,et al.  An analysis of the coupling of finite-element and Nyström methods in acoustic scattering , 1994 .

[33]  Jianlin Xia,et al.  On 3D modeling of seismic wave propagation via a structured parallel multifrontal direct Helmholtz solver , 2011 .

[34]  A. Kirsch,et al.  Convergence analysis of a coupled finite element and spectral method in acoustic scattering , 1990 .