Transient electromagnetic scattering from dielectric objects using the electric field Integral equation with Laguerre polynomials as temporal basis functions

In this paper, we propose a time-domain electric field integral equation (TD-EFIE) formulation for analyzing the transient electromagnetic response from three-dimensional (3-D) dielectric bodies. The solution method in this paper is based on the Galerkin's method that involves separate spatial and temporal testing procedures. Triangular patch basis functions are used for spatial expansion and testing functions for arbitrarily shaped 3-D dielectric structures. The time-domain unknown coefficients of the equivalent electric and magnetic currents are approximated using a set of orthonormal basis function that is derived from the Laguerre functions. These basis functions are also used as the temporal testing functions. Use of the Laguerre polynomials as expansion functions for the transient portion of response enables one not only to handle the time derivative terms in the integral equation in an analytic fashion but also completely separates the space and the time variables. Thus, the time variable along with the Courant condition can be eliminated in a Galerkin formulation using this procedure. We also propose an alternative formulation using a different expansion of the magnetic current. The total computational cost for this new method is similar to that of an implicit marching-on in time (MOT)-EFIE scheme, even though at each step this procedure requires more computations. Numerical results involving equivalent currents and far fields computed by the two proposed methods are presented and compared.     

导体目标瞬态电磁散射的时空矩量法研究

时域积分方程的矩量法是求解瞬态电磁散射的方法之一。研究了基于加权Lagurre函数和RWG基函数分别作为时间、空间基函数的时空矩量法,给出了时域磁场积分方程时空矩量法的全部计算公式,编制了相应串行和并行计算程序。计算结果表明：该方法具有很好的时域稳定特性,为宽带电磁散射分析提供了可能,同时也指出了其应用的局限性,为改进其方法提供了参考。     

Novel monopolar MFIE MoM-discretization for the scattering analysis of small objects

We present a novel method of moments (MoM)-magnetic field integral equation (MFIE) discretization that performs closely to the MoM-EFIE in the electromagnetic analysis of moderately small objects. This new MoM-MFIE discretization makes use of a new set of basis functions that we name monopolar Rao-Wilton-Glisson (RWG) and are derived from the RWG basis functions. We show for a wide variety of small objects -curved and sharp-edged-that the new monopolar MoM-MFIE formulation outperforms the conventional MoM-MFIE with RWG basis functions.     

Calculation of CFIE impedance matrix elements with RWG and n/spl times/RWG functions

The method of moments (MoM) solution of combined field integral equation (CFIE) for electromagnetic scattering problems requires calculation of singular double surface integrals. When Galerkin's method with triangular vector basis functions, Rao-Wilton-Glisson functions, and the CFIE are applied to solve electromagnetic scattering by a dielectric object, both RWG and n/spl times/RWG functions (n is normal unit vector) should be considered as testing functions. Robust and accurate methods based on the singularity extraction technique are presented to evaluate the impedance matrix elements of the CFIE with these basis and test functions. In computing the impedance matrix elements, including the gradient of the Green's function, we can avoid the logarithmic singularity on the outer testing integral by modifying the integrand. In the developed method, all singularities are extracted and calculated in closed form and numerical integration is applied only for regular functions. In addition, we present compact iterative formulas for computing the extracted terms in closed form. By these formulas, we can extract any number of terms from the singular kernels of CFIE formulations with RWG and n/spl times/RWG functions.     

Solution of time domain electric field Integral equation using the Laguerre polynomials

In this paper, we propose a numerical method to obtain a solution for the time domain electric field integral equation (TD-EFIE) for arbitrary shaped conducting structures. This method does not utilize the customary marching-on in time (MOT) solution method often used to solve a hyperbolic partial differential equation. Instead we solve the wave equation by expressing the transient behaviors in terms of Laguerre polynomials. By using these causal orthonormal basis functions for the temporal variation, the time derivatives can be handled analytically. In order to solve the wave equation, we introduce two separate testing procedures, a spatial and temporal testing. By introducing first the Galerkin temporal testing procedure, the MOT procedure is replaced by a recursive relation between the different orders of the weighted Laguerre polynomials. The other novelty of this approach is that through the use of the entire domain Laguerre polynomials for the expansion of the temporal variation of the current, the spatial and the temporal variables can be separated and the temporal variables can be integrated out. For convenience, we use the Hertz vector as the unknown variable instead of the electric current density. To verify our method, we compare the results of a TD-EFIE and inverse Fourier transform of a frequency domain EFIE.     

Integral equation solution for analyzing scattering fromone-dimensional periodic coated strips

A set of integral equations based on the surface/surface formulation are developed for analyzing electromagnetic scattering by one-dimensional periodic structures. To compare the accuracy, efficiency, and robustness of the formulation, the electric field integral equation (EFIE), magnetic field integral equation (MFIE), and combined field integral equation (CFIE) are developed for analyzing the same structure for different excitations. Due to the periodicity of the structure, the integral equations are formulated in the spectral domain using the Fourier transform of the integrodifferential operators. The generalized-biconjugate-gradient-fast Fourier transform method with subdomain basis functions is used to solve the matrix equation     

The use of curl-conforming basis functions for the magnetic-field integral equation

Divergence-conforming Rao-Wilton-Glisson (RWG) functions are commonly used in integral-equation formulations to model the surface current distributions on planar triangulations. In this paper, a novel implementation of the magnetic-field integral equation (MFIE) employing the curl-conforming n~/spl times/RWG basis and testing functions is introduced for improved current modelling. Implementation details are outlined in the contexts of the method of moments, the fast multipole method, and the multilevel fast multipole algorithm. Based on the examples of electromagnetic modelling of conducting scatterers, it is demonstrated that significant improvement in the accuracy of the MFIE can be obtained by using the curl-conforming n~/spl times/RWG functions.     

Electromagnetic radiation from structures consisting of combinedbody of revolution and arbitrary surfaces

A formulation is developed to treat radiation from structures consisting of a body of revolution (BOR) in the presence of multiple arbitrarily shaped three-dimensional objects. An electric field integral equation is set up on the surface of the combined structure. The resulting integro-differential equation is solved using the method of moments. On the BOR, harmonic entire domain expansion functions are used for the circumferential dependence, while overlapping subdomain functions are used to model the axial curvature. The arbitrarily shaped portions of the structure are modeled using triangular surface patches. The resulting system matrix has a partial block diagonal nature, which provides a more economical solution for structures that have some rational symmetry. Numerical results are presented and compared to measurements of a unique cavity-backed patch fed antenna     

利用时域电场积分方程分析天线辐射问题

研究了基于矩量法和RWG三角基函数的隐式电场积分方程的时域算法,引入了一种激励源的时域设置方法.利用直接在时域加激励源的方法来分析天线辐射问题,所得时域数据经傅立叶变换可得到很宽频带的频域数据,与在频域逐个频点求解相比大大节省了计算时间.通过对几种典型天线的分析计算,验证了方法的正确性和算法的稳定性.     

基于RWG基函数的伽略金法中奇异性积分 的精确快速计算

利用电磁场积分方程的伽略金法求解理想导体电磁散射问题时需要计算奇异性的二重面积分(即4维积分).伽略金法的基函数和检验函数广泛采用RWG(Rao-Wilton-Glisson)矢量基函数.传统上采用奇异值提取技术和Duffy坐标变换法处理该奇异性积分,本文提出了一种更为精确和高效的计算方法,该新方法通过参数坐标变换、相对坐标变换、积分区域分解和广义Duffy坐标变换相结合的技术消除了被积函数的奇异性并降低了原4维奇异性积分的数值积分维数.通过计算实例证明该方法的精确性和高收敛特性.     

Marching-on-in-Degree Based Time-Domain Magnetic Field Integral Equation Method for Bodies of Revolution

A marching-on-in-degree (MOD) based time-domain magnetic field integral equation method for bodies of revolution (BOR) is proposed and applied to obtain the induced currents on perfectly electric conducting BOR. Before this work, the time-domain integral equation method for BOR based on a marching-on-in-time procedure cannot really reduce the computational cost, since the number of unknowns cannot really be reduced. But it is the unknown reduction that serves as the key point of cost saving in BOR-problems. The method implemented in this letter can really utilize the symmetric property of BOR by applying two sets of entire domain basis functions. One is a set of scaled Laguerre polynomials inherited from common MOD method and used as temporal basis functions. The other is a Fourier series which comes from frequency domain method for solving BOR-problems. The validity, efficiency, and stability of the method are verified by several numerical examples.     

The transient response of finite open circular cylinders

An eigenmode expansion formulation of the singularity expansion method based on the electric field integral equation is developed for the transient response of conducting finite open cylyinders. The eigenvalues and eigenfunctions of the impedance operator are calculated by the Galerkin method using entire domain expansion functions. The transient surface current density and backscattered far field in response to an incident electromagnetic pulse are calculated for cylinders of various aspect ratios.     

Loop star basis functions and a robust preconditioner for EFIE scattering problems

An electric field integral equation (EFIE) formulation using the loop-star basis functions has been developed for modeling plane wave scattering from perfect conducting objects. A stability analysis at the DC limit shows that the use of the Rao-Wilton-Glisson (RWG) basis functions results in a singular matrix operator. However, the use of the loop-star basis functions results in a well-conditioned matrix. Moreover, a preconditioner constructed from a two-step process, based on near interactions and an incomplete factorization with a heuristic drop strategy, has been proposed in conjunction with the conjugate gradient method to solve the resulting matrix equation. The approach is shown to be effective for resolving both the low frequency instability and the bad conditioning of the EFIE method. The computational complexity of the proposed approach is shown to be O(N/sup 2/).     

An efficient algorithm for electromagnetic scattering from rough surfaces using a single integral equation and multilevel sparse-matrix canonical-grid method

An efficient algorithm for wave scattering from two-dimensional lossy rough surfaces is proposed. It entails the use of a single magnetic field integral equation (SMFIE) in conjunction with a multilevel sparse-matrix canonical-grid (MSMCG) method. The Rao-Wilton-Glisson (RWG) triangular discretization is adopted to better model the rough surface than the pulse basis functions used in the well-established SMCG method. Using the SMFIE formulation, only one unknown per interior edge of the triangular mesh approximating the rough surface is required, and the iterative solution to the moment equation converges more rapidly than that of the conventional coupled equations for dielectric rough surfaces. The MSMCG method extends the applicability of the SMCG method to rougher surfaces. Parallel implementation of the proposed method enables us to model dielectric surfaces up to a few thousand square wavelengths. Simulation results are presented as bistatic scattering coefficients for Gaussian randomly rough surfaces.     

Propagation and coupling properties of integrated opticalwaveguides-an integral equation formulation

The propagation and coupling properties of integrated optical waveguides are analyzed by means of the electric field integral equation approach. The kernel of the integral equation is the Green's function of a two-layered medium. The Galerkin's method is then employed to solve the integral equation numerically. The set of basis and test functions consists of entire domain plane wave functions. Fast convergence and superior accuracy are the advantages of the chosen set of basis and test functions. The method is used to compute the propagation and coupling properties of several structures. Very good agreement is observed with previously published results. Field distributions of several coupled mode structures, such as the symmetrical and asymmetrical coupler are also investigated and presented. Finally, the same method is used to produce the field distribution of waveguides having more complex cross section like the trapezoidal waveguide     

空域MPIE/MoM全波分析三维微带电路

矩量法在空域分析三维微带电路的关键是闭式格林函数的求解。本文首先介绍谱域格林函数一种新的表达式,使得源点和场点在同一层时,离散复镜像法可提取出与它们无关的闭式,从而避免了插值;不在同一层时,可提取与场点无关的闭式,此时只须对源点进行一维插值,因而提高了计算效率。然后,利用闭式格林函数和RWG基函数,基于混合位积分方程的矩量法就可以在空域精确、有效地分析三维任意形状微带电路。给出了几个典型实例,表明本文方法的有效性。     

Integral equation analysis of dielectric and conducting bodies ofrevolution in the presence of arbitrary surfaces

A formulation is developed to treat radiation from structures consisting of conducting and/or dielectric bodies of revolution (BOR) in the presence of multiple arbitrary shaped three-dimensional objects. A set of integral equations is developed on the surfaces of the combined structure and the resulting integro-differential equations are solved using the method of moments. On the BOR, harmonic entire domain expansion functions are used for the circumferential dependence, while overlapping sub-domain functions are used to model axial curvature. The arbitrary shaped portions of the structure are modeled using triangular surface patch basis functions. The resulting matrix has a partial block diagonal nature which provides a more economical solution for structures which have some rotational symmetry. The accuracy of the BOR and arbitrary surface formulation is verified using the self-consistency method and measured data     

The three dimensional weak form of the conjugate gradient FFTmethod for solving scattering problems

The problem of electromagnetic scattering by a three-dimensional dielectric object can be formulated in terms of a hypersingular integral equation, in which a grad-div operator acts on a vector potential. The vector potential is a spatial convolution of the free space Green's function and the contrast source over the domain of interest. A weak form of the integral equation for the relevant unknown quantity is obtained by testing it with appropriate testing functions. The vector potential is then expanded in a sequence of the appropriate expansion functions and the grad-div operator is integrated analytically over the scattering object domain only. A weak form of the singular Green's function has been used by introducing its spherical mean. As a result, the spatial convolution can be carried out numerically using a trapezoidal integration rule. This method shows excellent numerical performance     

Improved Stable Scheme for the Time Domain Integral Equation Method

A novel technique for improving the stability of time domain electric field integral equation (TD-EFIE) method is presented. The method involving the introduction of a new class of temporal basis functions employs a modified form of TD-EFIE to generate stable transient responses from arbitrarily shaped conducting objects, which possesses the advantage of simplicity in implementation. Numerical examples are given to demonstrate the accuracy and the efficiency of the proposed method     

A stable solution of time domain electric field Integral equation for thin-wire antennas using the Laguerre polynomials

In this paper, a numerical method to obtain an unconditionally stable solution of the time domain electric field integral equation for arbitrary conducting thin wires is presented. The time-domain electric field integral equation (TD-EFIE) technique has been employed to analyze electromagnetic scattering and radiation problems from thin wire structures. However, the most popular method to solve the TD-EFIE is typically the marching-on in time (MOT) method, which sometimes may suffer from its late-time instability. Instead, we solve the time-domain integral equation by expressing the transient behaviors in terms of weighted Laguerre polynomials. By using these orthonormal basis functions for the temporal variation, the time derivatives can be handled analytically and stable results can be obtained even for late-time. Furthermore, the excitation source in most scattering and radiation analysis of electromagnetic systems is typically done using a Gaussian shaped pulse. In this paper, both a Gaussian pulse and other waveshapes like a rectangular pulse or a ramp like function have been used as excitations for the scattering and radiation of thin-wire antennas with and without junctions. The time-domain results are compared with the inverse discrete Fourier transform (IDFT) of a frequency domain analysis.