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.     

位积分方程组矩量法用于精确计算非均匀介质柱的电磁散射

位积分方程组的主要特点是以电磁位为未知函数,这些未知函数在具有不同电磁参数的介质分界面处是连续的,因而在矩量法的实现过程中能够非常方便地应用高阶插值基函数来展开未知函数,以便获得高精度的解。但是,经典的点匹配方案使该模型的数值稳定性较差。本文用位积分方程组矩量法模型计算任意截面非均匀介质柱的电磁散射,采用三角形离散方案和高阶插值基函数,在测试过程中应用新提出的测试方法,克服了原位方程组矩量法模型的数值不稳定性。对矩量法矩阵中自阻抗元素的奇异性处理方法也作了详细介绍。文中提供的数值结果表明,该方法是精确、稳定的。     

Wide-band electromagnetic scattering from a dielectric BOR buriedin a layered lossy dispersive medium

A method of moments (MoM) analysis is developed for electromagnetic scattering from a dielectric body of revolution (BOR) embedded in a layered medium (the half-space problem constituting a special case). The layered-medium parameters can be lossy and dispersive, of interest for simulating soil. To make such an analysis tractable for the wide-band (short-pulse) applications of interest here, we have employed the method of complex images to evaluate the Sommerfeld integrals characteristic of the dyadic layered-medium Green's function. Example wide-band scattering results are presented, wherein fundamental wave phenomenology is elucidated. Of particular interest, we consider wide-band scattering from a model plastic mine, buried in soil, with the soil covered by a layer of snow     

The use of the transfinite interpolation in the method of momentsapplied to electromagnetic scattering by dielectric cylinders

The method of moments (MoM) solution of electromagnetic scattering presents two major numerical difficulties: the number of unknowns and the computation time necessary to calculate the matrix elements. To circumvent these problems, a MoM using the transfinite interpolation and a reduced integration scheme is presented here. The so-called h and p versions of the new method are applied to the scattering of an electromagnetic wave by an infinite dielectric cylinder (TM case) in the Richmond formulation. The transfinite and classical methods are compared in terms of the convergence rates of the radar cross section and of the total electric field inside the dielectric. The results confirm the superiority of the new schemes as predicted by the theory     

Solution of combined-field integral equation using multilevel fastmultipole algorithm for scattering by homogeneous bodies

We present an accurate method of moments (MoM) solution of the combined field integral equation (CFIE) using the multilevel fast multipole algorithm (MLFMA) for scattering by large, three-dimensional (3-D), arbitrarily shaped, homogeneous objects. We first investigate several different MoM formulations of the CFIE and propose a new formulation, which is both accurate and free of interior resonances. We then employ the MLFMA to significantly reduce the memory requirement and computational complexity of the MoM solution. Numerical results are presented to demonstrate the accuracy and capability of the proposed method. The method can be extended in a straightforward manner to scatterers composed of different homogeneous dielectric and conducting objects     

Efficient computation for the general solution of a slot loaded by a hemispherical dielectric and/or backing cavity

The problem of a slot loaded by a hemispherical dielectric and/or a hemispherical backing cavity is solved using the method of moments (MoM). A novel recurrence formula is presented for fast calculations of the MoM admittance integrals associated with the spherical homogeneous solutions. The method enables the admittance integrals to be calculated without the need for any numerical integration. Measurements were carried out and good agreement between theory and experiment is obtained.     

Radiation and scattering from infinite periodic printed antennaswith inhomogeneous media

The hybrid method of moments (MoM)/Green's function method technique is applied to infinite periodic printed antenna arrays containing dielectric inhomogeneities. The solution uses an integral equation for an infinite periodic printed array on or over a homogeneous dielectric substrate, coupled with equivalent volume polarization currents for dielectric inhomogeneities on top of the homogeneous substrate. Volume pulse-basis functions were used to expand the volume polarization currents. A hybrid MoM/Green's function method solution was then obtained through the matrix form of the problem. The two-dimensional (2-D) solution of plane wave scattering from a grounded dielectric slab was used to validate the reaction impedance of the dielectric inhomogeneity. Several infinite periodic printed dipole arrays with dielectric supports and overlays were studied with this solution and good agreement was observed between the hybrid MoM/Green's function method and waveguide simulator experiments     

Electromagnetic scattering interference between two shallow objects buried under 2-D random rough surfaces

A rigorous electromagnetic model has been used to analyze the scattering from two dielectric shallow objects buried under the two-dimensional (2-D) random rough ground (3-D scattering problem) as a means of predicting false alarms. The method of moments (MoM) accelerated by the steepest descent fast multipole method (SDFMM) is used to compute the equivalent electric and magnetic surface currents on all scatterers (i.e., the rough ground and the two buried objects). The roughness parameters influence the scattering interference mechanism of the two objects, however, a large separation distance (e.g., several correlation lengths) showed stronger effect for small ground roughness.     

A Mesh-Adapted Closed-Form Regular Kernel for 3D Singular Integral Equations

The Green's functions employed in the method of moments (MoM) diverge when observation and source points coincide; this is at the origin of the difficulties in computing the MoM matrix entries, and in handling the near-field interactions in fast Fourier transform (FFT)-based fast methods and other sampling-based methods. In this paper, we show that this singularity can be avoided, and a modified regular Green's function can be used instead. This latter is obtained from the spectral representation of the usual Green's function via windowing of its spectrum; the width of the spectral window depends on the size of the mesh employed for discretizing the problem, so that the proposed regular Green's function is a mesh-adapted regular kernel. We address a general 3D problem; we relate the MoM reaction integrals to the 2D Fourier spectrum of the Green's function, that allows to discuss the necessary spectral bandwidth for the windowed Green's function. We employ a tapered window, and present a closed-form expression for the spatial Green's function. Numerical results are presented for 3D antenna and scattering problems discretized with Rao-Wilton-Glisson (RWG) functions, and for uniform and nonuniform meshing. They show that the proposed method yields accurate solutions also for the antenna input impedance. The meaning of the regularized Green's function is also discussed and put in perspective.     

On the solution of a class of large body problems with full orpartial circular symmetry by using the finite-difference time-domain(FDTD) method

This paper presents an efficient method to accurately solve large body scattering problems with partial circular symmetry. The method effectively reduces the computational domain from three to two dimensions by using the reciprocity theorem. It does so by dividing the problem into two parts: a larger 3-D region with circular symmetry, and a smaller 2-D region without circular symmetry. An finite-difference time-domain (FDTD) algorithm is used to analyze the circularly symmetric 3-D case, while a method of moments (MoM) code is employed for the nonsymmetric part of the structure. The results of these simulations are combined via the reciprocity theorem to yield the radiation pattern of the composite system. The advantage of this method is that it achieves significant savings in computer storage and run time in performing an equivalent 2-D as opposed to a full 3-D FDTD simulation. In addition to enhancing computational efficiency, the FDTD algorithm used in this paper also features one improvement over conventional FDTD methods: a conformal approach for improved accuracy in modeling curved dielectric and conductive surfaces. The accuracy of the method is validated via a comparison of simulated and measured results     

An accurate method for the calculation of singular integralsarising in time-domain integral equation analysis of electromagneticscattering

An accurate method for the evaluation of the Cauchy principal value integrals arising in time-domain electromagnetic wave scattering is presented. This is applied to a boundary integral equation (BIE) method employing quadratic curvilinear surface elements where such singularities do not vanish (as they generally do when simpler but less accurate discretizations are employed). The technique involves weakening the singularity in the original kernel to a degree where conventional integration methods may be employed and transforming the strong singularity to a line integral in a form which allows cancellation of its singular components. The effectiveness of the method is demonstrated by scattering calculations on a variety of targets and the computational cost of the approach is trivial     

金属介质混合目标电磁特性的快速分析

将自适应积分算法与基于体面混合积分方程的矩量法相结合快速分析任意结构金属/介质混合目标的电磁散射和辐射特性.通过将传统矩量法的阻抗矩阵分为两部分且采用不同的方法进行处理计算,提高了矩量法的计算速度并大幅度缩减了需要的计算机内存占用量.最后,分别用传统的矩量法与结合自适应积分快速算法的矩量法计算了三个典型例子,通过比较充分说明了文中方法的有效性.     

Some convergence considerations in space-domain moment-methodanalysis of a class of wide-band microstrip antennas

The method of moments (MoM) analysis of probe-fed rectangular microstrip patches requires the inclusion of a probe-to-patch attachment mode-expansion function when the substrate thickness d⩾0.02λ, where λ is the free-space wavelength. The results for the input impedance showed increased divergence with measurements when the attachment mode was omitted from the full-wave analysis. The attachment mode can be expressed as an infinite eigenfunction series that increases the fill time of the impedance matrix in an MoM analysis. In an earlier investigation, the infinite eigenfunction series was reduced to a residue series that required one or two terms compared to about 55 terms for the eigenfunction series. In this paper, the convergence properties of the eigenfunction and residue series are investigated in view of rigorous MoM analysis. The relative errors resulting from replacing the eigenfunction by the residue series for the attachment mode, are compared by numerically evaluating a class of two-dimensional (2-D) spatial integrals shown to be closely related to the elements of an MoM impedance matrix. Additionally, the computation times for the evaluation of these integrals for the two forms of the attachment mode-expansion function are also included. Based on the superior convergence properties of the residue series for the attachment mode-expansion function, it is mathematically justified that this form can readily be used for analytic reduction of the spatial, reaction integrals from four to 2-D forms. This feature allows further reduction of the fill time of the MoM impedance matrix, suggesting the possibility of developing an efficient space-domain MoM technique for modeling of wide-band microstrip antennas     

Efficient modelling of finite-size ungrounded PCBs with a full wave3D moment-method analysis

A method in which 2D sheet-impedance elements are used to represent the finite-size dielectric substrate of an ungrounded PCB within a full wave 3D moment-method analysis (MoM) is presented. The use of these 2D elements significantly reduces both the core-memory demands and the CPU requirements of the MoM by an order of magnitude, with little compromise in accuracy     

Scattering from a circular dielectric post embedded in a groundeddielectric sheet waveguide

A systematic method for obtaining the scattered electric field in a grounded dielectric sheet waveguide is presented. It is shown that the point-matching method can be used for an explicit calculation of the integral equation to estimate the scattering from a circular dielectric post embedded in the grounded sheet. The magnitude of the reflection coefficient as a function of dielectric constant is given     

Analytical technique to evaluate the asymptotic part of theimpedance matrix of Sommerfeld-type integrals

The propagation of electromagnetic waves in a grounded dielectric slab has numerous applications in printed antenna technology and in the analysis of microwave- and millimeter-wave integrated circuits. For the accurate analysis of microstrip dipoles and circuits based on the moment of method (MoM), a crucial step is the precise evaluation of the impedance matrix elements which contain the integration of Sommerfeld-type integrals. The integral transform method with the asymptotic extraction technique is formulated for calculating a Sommerfeld-type integral problem. This formulation allows the infinite double integral of the asymptotic part of the impedance matrix to be transformed into a finite one-dimensional (1-D) integral. This finite 1-D integral contains a spherical Legendre function and can be easily evaluated numerically after the singular part of the integral is performed analytically. It is shown that the proposed method dramatically reduces the computation time and improves the accuracy over the conventional method to evaluate the asymptotic part of impedance matrix     

Étude des structures planaires en très hautes fréquences par une nouvelle formulation de la méthode des moments dans le domaine spatial

An efficient analytical integration technique for computation of spatial method of moments (MoM) integrals in conjunction with numerical matched loads is presented. The current distribution on the device is solved by using the well-known Galerkin’s MoM procedure applied to mixed potential integral equation in the spatial domain. The scattering parameters are determined by considering infinite lines at each port where only the fundamental mode is assumed to propagate. The contribution of this work is the development of an integration technique for the computation of spatial domain integrals, that is fast and rigorous. This technique is based on a Taylor series expansion of the integrands involving only polynomial functions. The use of polynomial forms in the integrals leads an immediate analytical integration, and the computation time will be considerably reduced.     

Scattering from a double-strip grating: rigorous equivalent networkformulation

An exact solution for the problem of transverse electric (TE) or transverse magnetic (TM) plane-wave scattering from a periodic, planar double-strip grating at a dielectric interface is described. The metal-strip grating is assumed to be perfectly conductive and infinite in length, with two different strips within a unit-cell. The formulation is based on a multimode equivalent network representation, and uses a rigorous solution for the relevant integral equation that extends the novel solution developed previously for the single-strip grating. Expressions for the elements of the multimode coupling matrices are given, together with a comparison of results for power transmitted through the grating, obtained by using the networks developed with the present method and a simple point-matching solution. Results are presented to illustrate the differences between single and double-strip gratings     

Numerical scattering analysis for two-dimensional dense randommedia: characterization of effective permittivity

A new numerical method for determining effective permittivity of dense random media in two dimensions is presented. The core of the method is to compare the average scattered field of a random collection of scatterers confined within an imaginary boundary with the scattered field from a homogeneous dielectric of the same shape as the imaginary boundary. The two-dimensional (2-D) problem is aggressively studied to provide insight into the dependence of the method's convergence on particle size, boundary shape, and boundary dimension. A novel inverse scattering method is introduced based on the method of moments (MoM), which greatly reduces the computation time and increases the flexibility of the procedure to analyze a variety of geometries. Results from this 2-D method may be used directly to compare with theoretical methods for determining effective permittivity such as the Polder-Van Santen (1946) mixing formula or field techniques such as the quasi-crystalline approximation     

介质体电磁散射的偶极子模型法研究

提出用偶极子模型法来分析介质体的电磁散射.该方法以矩量法和Schaubert-Wilton-Glisson(SWG)基函数为基础,把介质体剖分成一定数量的四面体元.在介质体内,把含有公共面的体元对等效成电偶极子;在介质体表面,把边界面及其对应的体元等效成电偶极子.当等效偶板子单元离观察点大于临界距离时,用偶极子模型法计算阻抗矩阵元素.偶极子模型法简单易操作,不仅能大幅度降低阻抗矩阵的计算时间,还简化了边界条件的处理.数值结果表明了该方法的高效性及与原方法几乎相同的计算精度.