阻抗劈散射的MM-PO混合算法研究

本文系统地研究了各向异性阻抗劈绕射的矩量法-物理光学(MM-PO)混合算法.首先研究了任意各向异性阻抗面的物理光学模型,推导出表面物理光学等效电磁流计算式.其次,提出了一种有效的含Hankel函数的弱振荡被积函数无穷积分处理方法.最后,将作者已公开发表的修正绕射电流基函数用于各向异性阻抗劈散射场研究,数值结果和已知的一致性绕射理论结果高度吻合.     

Moment method analysis of the electric shielding factor of aconducting TM shield at ELF

This paper presents an integral equation and method of moments (MM) solution to the problem of TM transmission by a metallic conducting shield at extremely low frequencies (ELF). To obtain an accurate and efficient solution, the equivalent volume polarization currents representing the shield are expanded in terms of physical basis functions, corresponding to two planewaves propagating normal to the surface of the shield. ELF approximations are used to obtain closed form expressions for certain crucial elements in the MM matrix equation where extreme accuracy between the relative magnitudes of self and mutual impedance terms are required. Numerical data will illustrate that despite the fact that the equivalent polarization currents are being computed very accurately, the method is not capable of computing the extreme near zone electric fields of these currents with sufficient accuracy to compute the electric shielding factor. An alternate method, based upon the use of the volume equivalence theorem to directly compute the total electric field in the shield is devised, and is found to accurately compute the electric shielding factor     

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.     

Electromagnetic Scattering From Arbitrarily Shaped Dielectric Bodies Using Paired Pulse Vector Basis Functions and Method of Moments

A pair of orthogonal pulse vector basis functions is demonstrated for the calculation of electromagnetic scattering from arbitrarily-shaped material bodies. The basis functions are intended for use with triangular surface patch modeling applied to a method of moments (MoM) solution. For modeling the behavior of dielectric materials, several authors have used the same set of basis functions to represent equivalent electric and magnetic surface currents. This practice can result in zero-valued or very small diagonal terms in the moment matrix and an unstable numerical solution. To provide a more stable solution, we have developed orthogonally placed, pulse basis vectors: one for the electric surface current and one for the magnetic surface current. This combination ensures strongly diagonal moment matrices. The basis functions are suitable for electric field integral equation (EFIE), magnetic field integral equation (HFIE), and combined field formulations. In this work, we describe the implementations for EFIE and HFIE formulations and show example results for canonical figures.      

Hybrid solutions for large-impedance coated bodies of revolution

Electromagnetic scattering solutions are developed for coated perfectly conducting bodies of revolution (BOR) that satisfy the impedance boundary condition. The integral equation arising from the impedance (Leontovich) boundary condition is solved by use of the method of moments (MM) technique along with an Ansatz for the surface currents that is derived from physical optics (PO) and the Fock theory that is modified for imperfectly conducting surfaces. The MM solution is expressed in terms of two integral (Galerkin) operators. The form of the Galerkin expansion used results in a symmetric MM system matrix. The hybrid solution is specialized for BOR's although the approach is applicable to a broader class of scatterers as well. The results are compared with the Mie solution for penetrable spherical scatterers, which satisfy the impedance boundary condition, and with recently published MM solutions for nonspherical scatterers.     

Electromagnetic Transmission through a Filled Slit in a Conducting Plane of Finite Thickness, TE Case

A solution is developed for computing the transmission characteristics of a slit in a conducting screen of finite thickness placed between two different media. The slit may be filled with Iossy material while the two regions on either side of the screen are assumed Iossless. A magnetic line source excitation is used (TE case) which is parallel to the axis of the slit. The equivalence principle is invoked to replace the two slit faces by equivalent magnetic current sheets on perfect electric conductors. Two coupled integral equations containing the magnetic currents as unknowns are then obtained and solved for by the method of moments. Pulses are used for the expansion and testing functions. Quantities computed are equivalent magnetic currents, the transmission coefficient, the gain pattern, and the normalized far field pattern.     

TM scattering by an impedance sheet extension of a paraboliccylinder

An integral equation and method of moments (MM) solution are presented for the two-dimensional (2-D) problem of transverse magnetic (TM) scattering by an impedance-sheet extension of a perfectly conducting parabolic cylinder. An integral equation is formulated for a dielectric cylinder of general cross section in the presence of a perfectly conducting parabolic cylinder. It is then shown that the solution for a general dielectric cylinder considerably simplifies for the special case of TM scattering by a thin multilayered dielectric strip that can be represented as an impedance sheet. The solution is termed an MM/Green's function solution, where the unknowns in the integral equation are the electric surface currents flowing in the impedance sheet; the presence of the parabolic cylinder is accounted for by including its Green's function in the kernel of the integral equation. The MM solution is briefly reviewed, and expressions for the elements in the matrix equation and the scattered fields are given. Sample numerical results are provided     

Electromagnetic scattering from objects composed of multiplehomogeneous regions using a region-by-region solution

The paper presents an efficient procedure to calculate the electromagnetic field scattered by an inhomogeneous object consisting of N+1 linear isotropic homogeneous regions. The procedure is based on surface integral equation (SIE) formulations and the method of moments. The method of moments (MM) is used to reduce the integral equations for each homogeneous dielectric region into individual matrices. These matrices are each solved for the equivalent electric current in terms of the equivalent magnetic current. A simple algebraic procedure is used to combine these solutions and to solve for the magnetic current on the outer dielectric surfaces of the scatterer. With the magnetic current determined, the electric current on the outer surface of the scatterer is calculated. Because the matrix corresponding to each dielectric region is solved separately, the authors call this procedure the region-by-region method. The procedure is simple and efficient. It requires less computer storage and less execution time than the conventional MM approach, in which all the unknown currents are solved for simultaneously. To illustrate the use of the procedure, the bistatic and monostatic radar cross sections (RCS) of several objects are computed. The computed results are verified by comparison with results obtained numerically using the conventional numerical procedure as well as via the series solution for circular cylindrical structures. The possibility of nonunique solutions has also been investigated     

Evaluation of reflected fields at caustic regions using a set of GO equivalent line currents

A set of equivalent electric and magnetic line currents is derived which supplements the geometrical optics (GO) solution in the far zone whenever one of the surface principal radii becomes very large. These hypothetical currents lie along the specular line of the surface and are shown to produce the same result as the stationary phase contribution of the physical optics integral. An example of a systematic application of such equivalent currents for the computation of the scattered field from a complex structure is also demonstrated.     

TM and TE scattering by a dielectric/ferrite cylinder in the presence of a half-plane

An integral equation solution to the problem of transverse magnetic (TM) or transverse electric (TE) scattering by an isotropic dielectric/ferrite material cylinder in the presence of a perfectly conducting half-plane is presented. The technique is termed a method of moments (MM)/Green's function solution since the method of moments is used to determine the electric and magnetic polarization currents representing the material cylinder, while the presence of the half-plane is accounted for by including the half-plane Green's function in the kernel of the integral equations. Numerical results are presented for the echo width, material cylinder interior fields, and the surface impedance of a material slab on the surface of a half-plane.     

Radiation from wire antennas attached to bodies of revolution: The junction problem

A method of moments (MM) formulation for multiple wires attached to biodes of revolution (BOR) is presented. This formulation builds on previous MM analyses for BOR's and wires. A special junction basis function is introduced to represent the physical currents in the BOR-wire attachment (junction) region. In the formulation, currents are represented by 1) mode-dependent and piecewise-continuous expansion functions on the BOR, 2) piecewise-continuous functions on the wires, and 3) a continuity-preserving basis function for the junction region. The block diagonal character of the MM/BOR matrix formulation and the computational economy of the MM analysis for wires are retained. The analysis is demonstrated for surface-attached radiators such as monopoles and loops as well as surface-attached passive elements.     

Integral equation formulations for imperfectly conducting scatterers

Integral equation formulations are presented for characterizing the electromagnetic (EM) scattering interaction for nonmetallic surfaced bodies. Three different boundary conditions are considered for the surfaces: namely, the impedance (Leontovich), the resistive sheet, and its dual, the magnetically conducting sheet boundary. The integral equation formulations presented for a general geometry are specialized for bodies of revolution and solved with the method of moments (MM). The current expansion functions, which are chosen, result in a symmetric system of equations. This system is expressed in terms of two Galerkin matrix operators that have special properties. The solutions of the integral equation for the impedance boundary at internal resonances of the associated perfectly conducting scatterer are examined. The results are compared with the Mie solution for impedance-coated spheres and with the MM solutions of the electric, magnetic, and combined field formulations for impedance-coated bodies.     

Adaptive multiscale moment method (AMMM) for analysis of scatteringfrom three-dimensional perfectly conducting structures

The adaptive multiscale moment method (AMMM) is presented for the analysis of scattering from three-dimensional (3D) perfectly conducting bodies. This algorithm employs the conventional moment method (MM) using the subsectional triangular patch basis functions and a special matrix transformation, which is derived from solving the Fredholm equation of the first kind by the multiscale technique. This methodology is more suitable for problems where the matrix solution time is much greater than the matrix fill time. The widely used triangular patch vector basis functions developed by Rao et al., (1982), are used for expansion and testing functions in the conventional MM. The objective here is to compress the unknowns in existing MM codes, which solves for the currents crossing the edges of the triangular patch basis functions. By use of a matrix transformation, the currents, source terms, and impedance matrix can be arranged in the form of different scales. From one scale to another scale, the initial guess for the solution can be predicted according to the properties of the multiscale technique. AMMM can reduce automatically the size of the linear equations so as to improve the efficiency of the conventional MM. The basic difference between this methodology and the other wavelet-based techniques that have been presented so far is that we apply the compression not to the impedance matrix but to the solution itself directly in an iterative fashion even though it is an unknown. Two numerical results are presented, which demonstrate that the AMMM is a useful method for analysis of electromagnetic scattering from arbitrary shaped 3D perfectly conducting bodies     

A sheet impedance approximation for electrically thick multilayered shields

This paper describes an extension of the sheet impedance concept to treat inhomogeneous or multilayered shields that may be thick in terms of material shield wavelengths. For shields with magnetic materials, a simple relation between the equivalent electric and magnetic currents representing the shield is obtained. This allows the magnetic current to be treated as a dependent unknown and the electric current to be found as the solution of a single surface integral equation shown to be a perturbation of that for a perfect electric conducting (PEC) surface. By using the proper interior equivalent problem, it is shown that the method produces accurate and stable results for shielding by a rectangular box.     

Dynamic analysis of a microstrip line over a perforated groundplane

A full wave analysis is presented to compute the characteristic impedance and propagation constant of a microstrip line over a perforated ground plane. The perforations in the ground plane are modeled by equivalent magnetic currents. The method of moments is applied to solve the coupled integral equations for the unknown electric current on the microstrip line and the unknown magnetic currents in the apertures. The fields are formulated using the space domain Sommerfeld type Green's functions. The matrix pencil technique is used to obtain the amplitude and the propagation constant of the fundamental modes for both current and the voltage on the microstrip line. Typical numerical results are given     

A sheet impedance approximation for electrically thick materialshields

This paper describes a simple extension of the sheet impedance concept to treat electromagnetic (EM) shields that may be thick in terms of material shield wavelengths. For magnetic shields, a simple relation between the equivalent electric and magnetic currents representing the shield is obtained, and the electric current is found as the solution of a single surface integral equation that is shown to be a simple perturbation of that for a perfect electric conducting (PEC) surface. Finally, it is shown that the computation of the small interior fields of good shields requires the use of the proper PEC interior equivalent problem     

Equivalent currents for an aperture in an impedance surface

Expressions of equivalent magnetic and electric currents for an aperture in an impedance surface is derived in a self-consistent manner. Each equivalent current consists of the combination of the tangential electric and magnetic fields in the aperture, and is placed in front of the original aperture where the surface is extended to close the aperture. The result is particularly useful for problems involving apertures in an impedance ground plane     

利用高阶矢量基函数求解时域磁场积分方程

本文利用一种新的高阶矢量基函数求解了三维时域磁场积分方程,该基函数定义在一个曲边三角形贴片上并用拉格朗日插值多项式来表示每一个贴片内的未知电流密度.该基函数的实质就是将拉格朗日插值多项式的插值点选为高斯积分结点,极大地简化和加快了时域积分方程矩量法的繁琐的时间和空间积分运算;另外,该基函数不要求网格为规范网格,给复杂目标的网格剖分带来很大方便.在空间上利用点匹配方法求解了时域磁场积分方程,数值计算结果表明了该方法求解时域积分方程的精确性和高效性.     

Higher order impedance boundary conditions applied to scattering bycoated bodies of revolution

In this paper higher order impedance boundary conditions will be employed in the solution of scattering by coated conducting bodies of revolution. The higher order impedance solution reduces the total number of unknowns relative to the exact solution, and produces a system matrix which is less dense than that of the exact solution. The construction of the solution involves two distinct steps. In the first step the body of revolution is replaced by an equivalent set of electric and magnetic currents on its exterior surface which generate the true fields outside the body. An integral equation relating these currents through the free space Green's function is derived. Step two employs the higher order impedance boundary condition to relate the electric and magnetic currents on the surface of the body. This replaces the rigorous solution of the interior problem. The higher order impedance boundary conditions are derived by obtaining an exact impedance boundary condition in the spectral domain for the coated ground plane, approximating the impedances as ratios of polynomials in the transform variables, and employing the Fourier transform. The resulting spatial domain differential equations are solved in conjunction with the integral equation using the method of moments. Several examples of bistatic and monostatic radar cross section for coated bodies of revolution are used to illustrate the accuracy of the higher order impedance boundary condition solution relative to the standard impedance boundary condition solution and the exact solution. The effects of coating thickness, loss, and curvature on the accuracy of the solution are discussed     

Electromagnetic scattering from an arbitrarily shaped three-dimensional homogeneous chiral body

The method of moments technique for analyzing electromagnetic scattering from an arbitrarily shaped three-dimensional homogeneous chiral body is presented based on the combined field integral equations. The body is assumed to be illuminated by a plane wave. The surface equivalence principle is used to replace the body by equivalent electric and magnetic surface currents. These currents radiating in unbounded free space produce the correct scattered field outside. The negatives of these currents produce the correct total internal field, when radiating in an unbounded chiral medium. By enforcing the continuity of the tangential components of the total electric and magnetic fields on the surface of the body, a set of coupled integral equations is obtained for the equivalent surface currents. The surface of the body is modeled using triangular patches. The triangular rooftop vector expansion functions are used for both equivalent surface currents. The coefficients of these expansion functions are obtained using the method of moments. The mixed potential formulation for a chiral medium is developed and used to obtain explicit expressions for the electric and magnetic fields produced by surface currents. Numerical results for bistatic radar cross sections are presented for three chiral scatterers - a sphere, a finite circular cylinder, and a cube.