A Grid-Robust Higher-Order Multilevel Fast Multipole Algorithm for Analysis of 3-D Scatterers

Recently, a set of novel, grid-robust, higher-order vector basis functions were proposed for the method-of-moments (MoM) solution of integral equations of scattering. The evaluation of integrals in the MoM is greatly simplified due to the unique properties associated with these basis functions. Moreover, these basis functions do not require the edge of a given patch to be completely shared by another patch; thus, the resultant MoM is applicable even for defective meshes. In this article, these new basis functions are employed to solve integral equations for three-dimensional (3-D) mixed dielectric/conducting scatterers. The multilevel fast multipole algorithm (MLFMA) is incorporated to speed up the solution of the resultant matrix system, thereby leading to a grid-robust, higher-order MLFMA solution having an O(N log N) computational complexity, where N denotes the total number of unknowns. Numerical examples are presented to demonstrate the accuracy of the proposed method.

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

[2]  R. Harrington,et al.  A surface formulation for characteristic modes of material bodies , 1977 .

[3]  Roberto D. Graglia,et al.  Higher order interpolatory vector bases for computational electromagnetics," Special Issue on "Advanced Numerical Techniques in Electromagnetics , 1997 .

[4]  Andrew J. Poggio,et al.  CHAPTER 4 – Integral Equation Solutions of Three-dimensional Scattering Problems , 1973 .

[5]  Jian-Ming Jin,et al.  A novel grid-robust higher order vector basis function for the method of moments , 2000 .

[6]  Jian-Ming Jin,et al.  Fast and Efficient Algorithms in Computational Electromagnetics , 2001 .

[7]  Te-kao Wu,et al.  Scattering from arbitrarily‐shaped lossy dielectric bodies of revolution , 1977 .

[8]  Stephen D. Gedney.  High‐order method‐of‐moments solution of the scattering by three‐dimensional PEC bodies using quadrature‐based point matching , 2001 .

[9]  Stephen D. Gedney,et al.  High-Order NystrÖm Solution of the Volume - EFIE for TE-Wave Scattering , 2001 .

[10]  John J. Ottusch,et al.  Numerical Solution of the Helmholtz Equation in 2D and 3D Using a High-Order Nyström Discretization , 1998 .

[11]  Jian-Ming Jin,et al.  A novel implementation of multilevel fast multipole algorithm for higher order Galerkin's method , 2000 .

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

[13]  Jian-Ming Jin,et al.  A higher order multilevel fast multipole algorithm for scattering from mixed conducting/dielectric bodies , 2003 .

[14]  A. C. Woo,et al.  Benchmark radar targets for the validation of computational electromagnetics programs , 1993 .