Generalized Debye Sources-Based EFIE Solver on Subdivision Surfaces

The electric field integral equation is a well-known workhorse for obtaining fields scattered by a perfect electric conducting object. As a result, the nuances and challenges of solving this equation have been examined for a while. Two recent papers motivate the effort presented in this paper. Unlike traditional work that uses equivalent currents defined on surfaces, recent research proposes a technique that results in well-conditioned systems by employing generalized Debye sources (GDS) as unknowns. In a complementary effort, some of us developed a method that exploits the same representation for both the geometry (subdivision surface representations) and functions defined on the geometry, also known as isogeometric analysis (IGA). The challenge in generalizing GDS method to a discretized geometry is the complexity of the intermediate operators. However, thanks to our earlier work on subdivision surfaces, the additional smoothness of geometric representation permits discretizing these intermediate operations. In this paper, we employ both ideas to present a well-conditioned GDS-electric field integral equation. Here, the intermediate surface Laplacian is well discretized by using subdivision basis. Likewise, using subdivision basis to represent the sources results in an efficient and accurate IGA framework. Numerous results are presented to demonstrate the efficacy of the approach.

[1]  T. Hughes,et al.  Isogeometric analysis : CAD, finite elements, NURBS, exact geometry and mesh refinement , 2005 .

[2]  Zydrunas Gimbutas,et al.  The Decoupled Potential Integral Equation for Time‐Harmonic Electromagnetic Scattering , 2014, 1404.0749.

[3]  Charles L. Epstein,et al.  Debye Sources, Beltrami Fields, and a Complex Structure on Maxwell Fields , 2013, 1308.5425.

[4]  E. Michielssen,et al.  A Multiplicative Calderon Preconditioner for the Electric Field Integral Equation , 2008, IEEE Transactions on Antennas and Propagation.

[5]  Kathrin Abendroth,et al.  Computational Methods For Electromagnetics , 2016 .

[6]  Weng Cho Chew,et al.  Enhanced A-EFIE With Perturbation Method , 2010, IEEE Transactions on Antennas and Propagation.

[7]  W. Chew,et al.  Integral equation solution of Maxwell's equations from zero frequency to microwave frequencies , 2000 .

[8]  W. Chew Vector Potential Electromagnetics with Generalized Gauge for Inhomogeneous Media: Formulation (Invited Paper) , 2014 .

[9]  Karl F. Warnick,et al.  Numerical analysis for electromagnetic integral equations , 2008 .

[10]  G. Vecchi,et al.  A multiresolution method of moments for triangular meshes , 2005, IEEE Transactions on Antennas and Propagation.

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

[12]  P. Raviart,et al.  A mixed finite element method for 2-nd order elliptic problems , 1977 .

[13]  Charles T. Loop,et al.  Smooth Subdivision Surfaces Based on Triangles , 1987 .

[14]  Jos Stam,et al.  Evaluation of Loop Subdivision Surfaces , 2010 .

[15]  John A. Evans,et al.  Isogeometric analysis using T-splines , 2010 .

[16]  A high order locally corrected Nyström implementation of the Decoupled Potential Integral Equation , 2015, 2015 9th European Conference on Antennas and Propagation (EuCAP).

[17]  Andrew F. Peterson Mapped Vector Basis Functions for Electromagnetic Integral Equations , 2006, Mapped Vector Basis Functions for Electromagnetic Integral Equations.

[18]  Bert Jüttler,et al.  On numerical integration in isogeometric subdivision methods for PDEs on surfaces , 2016 .

[19]  A. Boag,et al.  Method of Generalized Debye Sources for the Analysis of Electromagnetic Scattering by Perfectly Conducting Bodies With Piecewise Smooth Boundaries , 2013, IEEE Transactions on Antennas and Propagation.

[20]  Giuseppe Vecchi,et al.  Loop-star decomposition of basis functions in the discretization of the EFIE , 1999 .

[21]  Yiying Tong,et al.  Subdivision based isogeometric analysis technique for electric field integral equations for simply connected structures , 2015, J. Comput. Phys..

[22]  L. J. Jiang,et al.  A new EFIE method based on Coulomb gauge for the low-frequency electromagnetic analysis , 2013 .

[23]  Weng Cho Chew,et al.  An integral equation method based on vector and scalar potential formulations , 2015, 2015 IEEE International Symposium on Antennas and Propagation & USNC/URSI National Radio Science Meeting.

[24]  Andrew F. Peterson,et al.  Higher-Order Techniques in Computational Electromagnetics , 2015 .

[25]  B. Kolundžija Electromagnetic modeling of composite metallic and dielectric structures , 1996, IEEE Antennas and Propagation Society International Symposium 1997. Digest.

[26]  B. Shanker,et al.  The Generalized Method of Moments for Electromagnetic Boundary Integral Equations , 2014, IEEE Transactions on Antennas and Propagation.

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

[28]  M. Taskinen,et al.  Current and charge Integral equation formulation , 2006, IEEE Transactions on Antennas and Propagation.

[29]  Jin-Fa Lee,et al.  A Discontinuous Galerkin Surface Integral Equation Method for Electromagnetic Wave Scattering From Nonpenetrable Targets , 2013, IEEE Transactions on Antennas and Propagation.

[30]  Charles L. Epstein,et al.  Debye Sources and the Numerical Solution of the Time Harmonic Maxwell Equations II , 2008, 0808.3369.

[31]  J. Volakis,et al.  Integral Equation Methods for Electromagnetics , 2012 .

[32]  B. Shanker,et al.  Generalized Method of Moments: A Novel Discretization Technique for Integral Equations , 2011, IEEE Transactions on Antennas and Propagation.