A novel grid-robust higher order vector basis function for the method of moments

A set of novel, grid-robust, higher order vector basis functions is proposed for the method-of-moments (MoM) solution of integral equations for three-dimensional (3-D) electromagnetic (EM) problems. These basis functions are defined over curvilinear triangular patches and represent the unknown electric current density within each patch using the Lagrange interpolation polynomials. The highlight of these basis functions is that the Lagrange interpolation points are chosen to be the same as the nodes of the well-developed Gaussian quadratures. As a result, the evaluation of the integrals in the MoM is greatly simplified. Additionally, the surface of an object to be analyzed can be easily meshed because the new basis functions do not require the side of a triangular patch to be entirely shared by another triangular patch, which is a very stringent requirement for traditional vector basis functions. The proposed basis functions are implemented with point matching for the MoM solution of the electric-field integral equation, the magnetic-field integral equation, and the combined-field integral equation. Numerical examples are presented to demonstrate the higher order convergence and the grid robustness for defective meshes using the new basis functions.

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

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

[3]  Jan S. Hesthaven,et al.  From Electrostatics to Almost Optimal Nodal Sets for Polynomial Interpolation in a Simplex , 1998 .

[4]  Stephen M. Wandzura,et al.  Electric Current Basis Functions for Curved Surfaces , 1992 .

[5]  M. G. Duffy,et al.  Quadrature Over a Pyramid or Cube of Integrands with a Singularity at a Vertex , 1982 .

[6]  J. Nédélec Mixed finite elements in ℝ3 , 1980 .

[7]  Jiming Song,et al.  Multilevel fast‐multipole algorithm for solving combined field integral equations of electromagnetic scattering , 1995 .

[8]  Elise de Doncker,et al.  New euler-maclaurin expansions and their application to quadrature over the s-dimensional simplex , 1979 .

[9]  D. A. Dunavant High degree efficient symmetrical Gaussian quadrature rules for the triangle , 1985 .

[10]  James N. Lyness,et al.  Moderate degree symmetric quadrature rules for the triangle j inst maths , 1975 .

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

[12]  J.M. Song,et al.  Higher-order MoM implementation to solve integral equations , 1999, IEEE Antennas and Propagation Society International Symposium. 1999 Digest. Held in conjunction with: USNC/URSI National Radio Science Meeting (Cat. No.99CH37010).

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

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