A Higher Order Hybrid SIE/FEM/SEM Method for the Flexible Electromagnetic Simulation in Layered Medium

A novel hybrid method is developed for the flexible and accurate electromagnetic simulation of penetrable objects in a layered medium (LM). In this method, the original complex simulation domain is first divided into several subdomains, following the spirit of divide-and-conquer. Each subdomain is then meshed and solved independently, where nonconformal mesh is inevitable. The Riemann type transmission condition is utilized at the interfaces of each subdomain to correctly exchange information so that the solutions of all subdomains converge rapidly to the real solution of the original problem. More specifically, in our method, the surface integral equation (SIE) combined with the LM Green’s functions (LMGFs) is adopted for the boundary subdomain, while the finite-element method (FEM) and the spectral element method (SEM) are employed for all the other interior dielectric subdomains. The SIE with LMGFs truncates the simulation domain tightly within the object itself, which drastically decreases the number of unknowns. The interior subdomains are modeled by either FEM or SEM, depending on the geometry and material property of each subdomain. To further enhance the simulation capability, higher order approaches are adopted for all the subdomain solvers in this hybrid method. Several numerical examples are demonstrated, where a high convergence and accuracy of this method is observed. This paper will serve as an efficient and flexible simulation tool for the applications of geophysical exploration.

[1]  J. P. Webb Hierarchal vector basis functions of arbitrary order for triangular and tetrahedral finite elements , 1999 .

[2]  Jianming Jin,et al.  On the formulation of hybrid finite-element and boundary-integral methods for 3-D scattering , 1998 .

[3]  Wei Hong,et al.  A Fast Domain Decomposition Method for Solving Three-Dimensional Large-Scale Electromagnetic Problems , 2008, IEEE Transactions on Antennas and Propagation.

[4]  W. Chew Waves and Fields in Inhomogeneous Media , 1990 .

[5]  Qing Huo Liu,et al.  Spectral Methods and Domain Decomposition for Nanophotonic Applications , 2013, Proceedings of the IEEE.

[6]  Evert Slob,et al.  Fast evaluation of zero‐offset Green's function for layered media with application to ground‐penetrating radar , 2007 .

[7]  Weng Cho Chew,et al.  Fast inhomogeneous plane wave algorithm for scattering from objects above the multilayered medium , 2001, IEEE Trans. Geosci. Remote. Sens..

[8]  Antonio Iodice,et al.  Electromagnetic Wave Scattering From Layered Structures With an Arbitrary Number of Rough Interfaces , 2009, IEEE Transactions on Geoscience and Remote Sensing.

[9]  Jin-Fa Lee,et al.  A domain decomposition approach for non-conformal couplings between finite and boundary elements for unbounded electromagnetic problems in R3 , 2007, J. Comput. Phys..

[10]  W. Chew,et al.  An Augmented Electric Field Integral Equation for Layered Medium Green's Function , 2011, IEEE Transactions on Antennas and Propagation.

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

[12]  Carretera de Valencia,et al.  The finite element method in electromagnetics , 2000 .

[13]  Tapan K. Sarkar,et al.  Green’s Function Using Schelkunoff Integrals for Horizontal Electric Dipoles Over an Imperfect Ground Plane , 2016, IEEE Transactions on Antennas and Propagation.

[14]  Sébastien Lambot,et al.  Full-Wave Modeling of Near-Field Radar Data for Planar Layered Media Reconstruction , 2014, IEEE Transactions on Geoscience and Remote Sensing.

[15]  W. Chew,et al.  A New Green's Function Formulation for Modeling Homogeneous Objects in Layered Medium , 2012, IEEE Transactions on Antennas and Propagation.

[16]  Richard Liu,et al.  A fast integral equation solver for 3D induction well logging in formations with large conductivity contrasts , 2013 .

[17]  Qing Huo Liu,et al.  A Hybrid FEM/MoM Method for 3-D Electromagnetic Scattering in Layered Medium , 2016, IEEE Transactions on Antennas and Propagation.

[18]  A. Bossavit Whitney forms: a class of finite elements for three-dimensional computations in electromagnetism , 1988 .

[19]  Stephen D. Gedney,et al.  Convolution PML (CPML): An efficient FDTD implementation of the CFS–PML for arbitrary media , 2000 .

[20]  N. Kinayman,et al.  Efficient use of closed-form Green's functions for the analysis of planar geometries with vertical connections , 1997 .

[21]  Yun Lin,et al.  Spectral Integral Method and Spectral Element Method Domain Decomposition Method for Electromagnetic Field Analysis , 2011 .

[22]  Xin-Qing Sheng,et al.  Application of Asymptotic Waveform Evaluation to Hybrid FE-BI-MLFMA for Fast RCS Computation Over a Frequency Band , 2013, IEEE Transactions on Antennas and Propagation.

[23]  M. Djordjevic,et al.  Higher Order Hybrid FEM-MoM Technique for Analysis of Antennas and Scatterers , 2009, IEEE Transactions on Antennas and Propagation.

[24]  Jian-Ming Jin,et al.  A Nonconformal FEM-DDM With Tree–Cotree Splitting and Improved Transmission Condition for Modeling Subsurface Detection Problems , 2014, IEEE Transactions on Geoscience and Remote Sensing.

[25]  Roland Griesmaier,et al.  An Asymptotic Factorization Method for Inverse Electromagnetic Scattering in Layered Media , 2008, SIAM J. Appl. Math..

[26]  Lawrence Carin,et al.  Multilevel fast-multipole algorithm for scattering from conducting targets above or embedded in a lossy half space , 2000, IEEE Trans. Geosci. Remote. Sens..

[27]  Tage Røsten,et al.  A 2.5D finite-element-modeling difference method for marine CSEM modeling in stratified anisotropic media , 2008 .

[28]  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.

[29]  K. Michalski,et al.  Electromagnetic scattering and radiation by surfaces of arbitrary shape in layered media. I. Theory , 1990 .

[30]  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.

[31]  Mahta Moghaddam,et al.  Electromagnetic scattering from a buried cylinder in layered media with rough interfaces , 2006 .

[32]  B. M. Kolundzija,et al.  Maximally Orthogonalized Higher Order Bases Over Generalized Wires, Quadrilaterals, and Hexahedra , 2013, IEEE Transactions on Antennas and Propagation.

[33]  Z. Cendes,et al.  A FEM domain decomposition method for photonic and electromagnetic band gap structures , 2006, IEEE Transactions on Antennas and Propagation.

[34]  V.I. Okhmatovski,et al.  Large-scale broad-band parasitic extraction for fast layout verification of 3-D RF and mixed-signal on-chip structures , 2004, IEEE Transactions on Microwave Theory and Techniques.

[35]  P. Kosmas,et al.  FDTD simulation of TE and TM plane waves at nonzero incidence in arbitrary Layered media , 2005, IEEE Transactions on Antennas and Propagation.

[36]  Jian-Ming Jin,et al.  An Accurate and Efficient Finite Element-Boundary Integral Method With GPU Acceleration for 3-D Electromagnetic Analysis , 2014, IEEE Transactions on Antennas and Propagation.

[37]  Davis A. Kopriva,et al.  Computation of electromagnetic scattering with a non‐conforming discontinuous spectral element method , 2002 .

[38]  Qing Huo Liu,et al.  A Hybrid SIM-SEM Method for 3-D Electromagnetic Scattering Problems , 2009, IEEE Transactions on Antennas and Propagation.

[39]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[40]  Jesús González A discontinuous Galerkin finite element method for the Time-Domain solution of Maxwell equations , 2014 .

[41]  Q.H. Liu,et al.  A 3-D spectral-element method using mixed-order curl conforming vector basis functions for electromagnetic fields , 2006, IEEE Transactions on Microwave Theory and Techniques.