Solving Electromagnetic Scattering from Conducting Objects with Novel Skeletonization Based Fast Direct Solver

In this paper, we present a novel skeletonization based fast direct solver for analyzing electromagnetic scattering from electrically conducting objects. The contributions of the present algorithm involve two aspects: reducing the cost for filling proxy matrices and enhancing the inverse process. To reduce the filling time of proxy matrices, the proxy surface is discretized by a series of points. The number of the points is considerably small compared with conventional method. Consequently, the dimensions of the proxy matrices for each group can be reduced remarkably in this way. On the other hand, a multiplicative factorization technique, recursive skeletonization factorization, is utilized in the inverse process. The projection matrices obtained by interpolative decomposition can be used directly with high efficiency. The remained block diagonal matrices after factorization can be easily solved by LU decomposition. The accuracy and efficiency of the proposed method will be validated by numerical results.

[1]  Jun Hu,et al.  A Butterfly-Based Direct Integral-Equation Solver Using Hierarchical LU Factorization for Analyzing Scattering From Electrically Large Conducting Objects , 2016, IEEE Transactions on Antennas and Propagation.

[2]  Per-Gunnar Martinsson,et al.  $\mathcal {O}(N)$ Nested Skeletonization Scheme for the Analysis of Multiscale Structures Using the Method of Moments , 2016, IEEE Journal on Multiscale and Multiphysics Computational Techniques.

[3]  Jin-Fa Lee,et al.  Multiscale electromagnetic computations using a hierarchical multilevel fast multipole algorithm , 2014 .

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

[5]  Per-Gunnar Martinsson,et al.  On the Compression of Low Rank Matrices , 2005, SIAM J. Sci. Comput..

[6]  Yaniv Brick,et al.  Fast Direct Solver for Essentially Convex Scatterers Using Multilevel Non-Uniform Grids , 2014, IEEE Transactions on Antennas and Propagation.

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

[8]  Jin-Fa Lee,et al.  A fast direct matrix solver for surface integral equation methods for electromagnetic wave scattering from non-penetrable targets , 2012 .

[9]  Lexing Ying,et al.  Hierarchical Interpolative Factorization for Elliptic Operators: Integral Equations , 2013, 1307.2666.

[10]  D. Wilton,et al.  Electromagnetic scattering by surfaces of arbitrary shape , 1980 .

[11]  Per-Gunnar Martinsson,et al.  O ( N ) Nested Skeletonization Scheme for the Analysis of Multiscale Structures Using the Method of Moments , 2017 .

[12]  Jun Hu,et al.  Fast Direct Solution of Integral Equations With Modified HODLR Structure for Analyzing Electromagnetic Scattering Problems , 2019, IEEE Transactions on Antennas and Propagation.

[13]  Jiming Song,et al.  Multilevel fast multipole algorithm for electromagnetic scattering by large complex objects , 1997 .