Integral Equation Based Domain Decomposition Method for Solving Electromagnetic Wave Scattering From Non-Penetrable Objects

The integral equation (IE) method is commonly utilized to model time-harmonic electromagnetic (EM) problems. One of the greatest challenges in its applications arises in the solution of the resulting ill-conditioned matrix equation. We introduce a new domain decomposition method (DDM) for the IE solution of EM wave scattering from non-penetrable objects. The proposed method is a non-overlapping/non-conformal DDM and it provides a computationally efficient and effective preconditioner for the IE matrix equations. Moreover, the proposed approach is very suitable for dealing with multi-scale electromagnetic problems since each sub-domain has its own characteristics length and will be meshed independently. Furthermore, for each sub-domain, we are free to choose the most effective IE sub-domain solver based on its local geometrical features and electromagnetic characteristics. Additionally, the multilevel fast multi-pole algorithm (MLFMA) is utilized to accelerate the computations of couplings between sub-domains. Numerical results demonstrate that the proposed method yields rapid convergence in the outer Krylov iterative solution process. Finally, simulations of several large-scale examples testify to the effectiveness and robustness of the proposed IE based DDM.

[1]  Jin-Fa Lee,et al.  Preconditioned Electric Field Integral Equation Using Calderon Identities and Dual Loop/Star Basis Functions , 2009, IEEE Transactions on Antennas and Propagation.

[2]  Yijun Lu,et al.  A domain decomposition finite-difference method for parallel numerical implementation of time-dependent Maxwell's equations , 1997 .

[3]  Jin-Fa Lee,et al.  Non-conformal domain decomposition method with second-order transmission conditions for time-harmonic electromagnetics , 2010, J. Comput. Phys..

[4]  R.J. Adams,et al.  Physical and analytical properties of a stabilized electric field integral equation , 2004, IEEE Transactions on Antennas and Propagation.

[5]  Jun Zhang,et al.  Sparse inverse preconditioning of multilevel fast multipole algorithm for hybrid Integral equations in electromagnetics , 2004 .

[6]  W. Wendland,et al.  Boundary Element Tearing and Interconnecting Domain Decomposition Methods , 2006 .

[7]  Jian-Ming Jin,et al.  A New Dual-Primal Domain Decomposition Approach for Finite Element Simulation of 3-D Large-Scale Electromagnetic Problems , 2007, IEEE Transactions on Antennas and Propagation.

[8]  B. Stupfel,et al.  A domain decomposition method for the vector wave equation , 2000 .

[9]  Eric de Sturler,et al.  Recycling Krylov Subspaces for Sequences of Linear Systems , 2006, SIAM J. Sci. Comput..

[10]  Jian-Ming Jin,et al.  An efficient algorithm for analyzing large-scale microstrip structures using adaptive integral method combined with discrete complex-image method , 2000 .

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

[12]  Jin-Fa Lee,et al.  Fast Computation of Angular Responses of Large-Scale Three-Dimensional Electromagnetic Wave Scattering , 2010, IEEE Transactions on Antennas and Propagation.

[13]  Yassine Boubendir,et al.  A FETI-like domain decomposition method for coupling finite elements and boundary elements in large-size problems of acoustic scattering , 2007 .

[14]  Vineet Rawat,et al.  One way domain decomposition method with second order transmission conditions for solving electromagnetic wave problems , 2010, J. Comput. Phys..

[15]  Weng Cho Chew,et al.  The use of Huygens' equivalence principle for solving 3-D volume integral equation of scattering , 1993 .

[16]  A. Buffa,et al.  A multiplicative Calderón preconditioner for the electric field integral equation , 2008, 2008 IEEE Antennas and Propagation Society International Symposium.

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

[18]  W. Chew,et al.  Wave-Field Interaction With Complex Structures Using Equivalence Principle Algorithm , 2007, IEEE Transactions on Antennas and Propagation.

[19]  E. Sturler,et al.  Truncation Strategies for Optimal Krylov Subspace Methods , 1999 .

[20]  Weng Cho Chew,et al.  Multiscale Simulation of Complex Structures Using Equivalence Principle Algorithm With High-Order Field Point Sampling Scheme , 2008, IEEE Transactions on Antennas and Propagation.

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

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

[23]  Weng Cho Chew,et al.  A multilevel fast multipole algorithm for electrically small composite structures , 2004 .

[24]  Jin-Fa Lee,et al.  A single-level dual rank IE-QR algorithm to model large microstrip antenna arrays , 2004 .

[25]  M. Maischak A multilevel additive Schwarz method for a hypersingular integral equation on an open curve with graded meshes , 2009 .

[26]  M. Jensen,et al.  A recursive Green's function method for boundary integral analysis of inhomogeneous domains , 1998 .

[27]  Patrick Amestoy,et al.  Task Scheduling in an Asynchronous Distributed Memory Multifrontal Solver , 2005, SIAM J. Matrix Anal. Appl..

[28]  Olaf Steinbach,et al.  Domain decomposition methods via boundary integral equations , 2000 .

[29]  Jin-Fa Lee,et al.  A Domain Decomposition Method With Nonconformal Meshes for Finite Periodic and Semi-Periodic Structures , 2007, IEEE Transactions on Antennas and Propagation.

[30]  Snorre H. Christiansen,et al.  A Preconditioner for the Electric Field Integral Equation Based on Calderon Formulas , 2002, SIAM J. Numer. Anal..

[31]  Olaf Steinbach,et al.  Inexact Fast Multipole Boundary Element Tearing and Interconnecting Methods , 2007 .

[32]  Jin-Fa Lee,et al.  A fast IE-FFT algorithm for solving PEC scattering problems , 2005 .