A Fast Direct Solver for High Frequency Scattering from a Large Cavity in Two Dimensions

We present a fast direct solver for the simulation of electromagnetic scattering from an arbitrarily shaped, large, empty cavity embedded in an infinite perfectly conducting half-space. The governing Maxwell equations are reformulated as a well-conditioned second kind integral equation, and the resulting linear system is solved in nearly linear time using a hierarchical matrix factorization technique. We illustrate the performance of the scheme with several numerical examples for complex cavity shapes over a wide range of frequencies.

[1]  Jianming Jin.  The Finite Element Method , 2010 .

[2]  Shivkumar Chandrasekaran,et al.  A Fast ULV Decomposition Solver for Hierarchically Semiseparable Representations , 2006, SIAM J. Matrix Anal. Appl..

[3]  J. Tukey,et al.  An algorithm for the machine calculation of complex Fourier series , 1965 .

[4]  V. Rokhlin,et al.  A fast direct solver for boundary integral equations in two dimensions , 2003 .

[5]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[6]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[7]  S. Chandler-Wilde,et al.  A boundary integral equation formulation for the Helmholtz equation in a locally perturbed half‐plane , 2005 .

[8]  Shivkumar Chandrasekaran,et al.  A Fast Solver for HSS Representations via Sparse Matrices , 2006, SIAM J. Matrix Anal. Appl..

[9]  Mario Bebendorf,et al.  F ¨ Ur Mathematik in Den Naturwissenschaften Leipzig Hierarchical Lu Decomposition Based Preconditioners for Bem Hierarchical Lu Decomposition Based Preconditioners for Bem , 2022 .

[10]  Leslie Greengard,et al.  Quadrature by expansion: A new method for the evaluation of layer potentials , 2012, J. Comput. Phys..

[11]  Piet Hut,et al.  A hierarchical O(N log N) force-calculation algorithm , 1986, Nature.

[12]  Johan Helsing,et al.  Corner singularities for elliptic problems: Integral equations, graded meshes, quadrature, and compressed inverse preconditioning , 2008, J. Comput. Phys..

[13]  Leslie Greengard,et al.  A Fast Direct Solver for Structured Linear Systems by Recursive Skeletonization , 2012, SIAM J. Sci. Comput..

[14]  H. Ammari,et al.  An integral equation method for the electromagnetic scattering from cavities , 2000 .

[15]  Weiwei Sun,et al.  A Fast Algorithm for the Electromagnetic Scattering from a Large Cavity , 2005, SIAM J. Sci. Comput..

[16]  M. Saunders,et al.  Solution of Sparse Indefinite Systems of Linear Equations , 1975 .

[17]  J. Bremer.  On the Nyström discretization of integral equations on planar curves with corners , 2012 .

[18]  P. Werner,et al.  The Helmholtz equation in disturbed half‐spaces , 1987 .

[19]  W. Hackbusch,et al.  On the fast matrix multiplication in the boundary element method by panel clustering , 1989 .

[20]  Per-Gunnar Martinsson,et al.  A Fast Direct Solver for a Class of Elliptic Partial Differential Equations , 2009, J. Sci. Comput..

[21]  Ralph E. Kleinman,et al.  Electromagnetic scattering by indented screens , 1994 .

[22]  John F. Shaeffer,et al.  Radar Cross Section , 2004 .

[23]  Calvin H. Wilcox,et al.  Integral Equation Method in Scattering Theory (David Colton and Rainer Kress) , 1984 .

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

[25]  Junshan Lin,et al.  Mode Matching for the Electromagnetic Scattering From Three-Dimensional Large Cavities , 2012, IEEE Transactions on Antennas and Propagation.

[26]  Sergej Rjasanow,et al.  Adaptive Cross Approximation of Dense Matrices , 2000 .

[27]  W. Arnoldi.  The principle of minimized iterations in the solution of the matrix eigenvalue problem , 1951 .

[28]  R. Freund,et al.  QMR: a quasi-minimal residual method for non-Hermitian linear systems , 1991 .

[29]  Leslie Greengard,et al.  A fast direct solver for scattering from periodic structures with multiple material interfaces in two dimensions , 2014, J. Comput. Phys..

[30]  Bradley K. Alpert,et al.  Hybrid Gauss-Trapezoidal Quadrature Rules , 1999, SIAM J. Sci. Comput..

[31]  Eric Darve,et al.  An O(NlogN)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal O (N \log N)$$\end{document} Fast Direct Solver fo , 2013, Journal of Scientific Computing.

[32]  James Bremer,et al.  A Nonlinear Optimization Procedure for Generalized Gaussian Quadratures , 2010, SIAM J. Sci. Comput..

[33]  Gang Bao,et al.  Optimal Shape Design of a Cavity for Radar Cross Section Reduction , 2014, SIAM J. Control. Optim..

[34]  Jian-Ming Jin,et al.  A special higher order finite-element method for scattering by deep cavities , 2000 .

[35]  Jun Lai,et al.  Radar Cross Section Reduction of a Cavity in the Ground Plane , 2014 .

[36]  W. Hackbusch,et al.  Introduction to Hierarchical Matrices with Applications , 2003 .

[37]  Per-Gunnar Martinsson,et al.  Fast direct solvers for integral equations in complex three-dimensional domains , 2009, Acta Numerica.

[38]  G..  A Theory for Multiresolution Signal Decomposition : The Wavelet Representation , 2004 .

[39]  R. Coifman,et al.  The fast multipole method for the wave equation: a pedestrian prescription , 1993, IEEE Antennas and Propagation Magazine.

[40]  M. Vouvakis,et al.  The adaptive cross approximation algorithm for accelerated method of moments computations of EMC problems , 2005, IEEE Transactions on Electromagnetic Compatibility.

[41]  James Bremer,et al.  An adaptive fast direct solver for boundary integral equations in two dimensions , 2009 .

[42]  Stéphane Mallat,et al.  A Theory for Multiresolution Signal Decomposition: The Wavelet Representation , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

[43]  Eric Darve,et al.  A fast block low-rank dense solver with applications to finite-element matrices , 2014, J. Comput. Phys..

[44]  Per-Gunnar Martinsson,et al.  Randomized algorithms for the low-rank approximation of matrices , 2007, Proceedings of the National Academy of Sciences.

[45]  Leslie Greengard,et al.  Stable and accurate integral equation methods for scattering problems with multiple material interfaces in two dimensions , 2012, J. Comput. Phys..

[46]  Eric Darve,et al.  An $$\mathcal O (N \log N)$$O(NlogN)  Fast Direct Solver for Partial Hierarchically Semi-Separable Matrices , 2013 .

[47]  H. Ammari,et al.  Analysis of the electromagnetic scattering from a cavity , 2002 .

[48]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[49]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[50]  Gang Bao,et al.  Stability of the Scattering from a Large Electromagnetic Cavity in Two Dimensions , 2012, SIAM J. Math. Anal..

[51]  W. Hackbusch,et al.  An introduction to hierarchical matrices , 2001 .

[52]  Jian-Ming Jin,et al.  The Finite Element Method in Electromagnetics , 1993 .

[53]  Jianming Jin,et al.  A finite-element-boundary integral formulation for scattering by three-dimensional cavity-backed apertures , 1991 .

[54]  Roland W. Freund,et al.  A Transpose-Free Quasi-Minimal Residual Algorithm for Non-Hermitian Linear Systems , 1993, SIAM J. Sci. Comput..

[55]  L. Greengard,et al.  Regular Article: A Fast Adaptive Multipole Algorithm in Three Dimensions , 1999 .

[56]  L. Greengard,et al.  A new version of the Fast Multipole Method for the Laplace equation in three dimensions , 1997, Acta Numerica.

[57]  H. Ammari,et al.  A cavity problem for Maxwell's equations , 2002 .

[58]  G. Chandler,et al.  Galerkin's method for boundary integral equations on polygonal domains , 1984, The Journal of the Australian Mathematical Society. Series B. Applied Mathematics.

[59]  Gang Bao,et al.  An improved mode-matching method for large cavities , 2005 .