A hybrid equation approach for the solution of electromagneticscattering problems involving two-dimensional inhomogeneous dielectriccylinders

A numerical procedure for the solution of electromagnetic scattering problems involving inhomogeneous dielectric cylinders of arbitrary cross section is discussed. The cases of illumination by both transverse magnetic (TM) and transverse electric (TE) plane waves are considered. The scattering problems are modeled via a hybrid integral-equation/partial-differential-equation approach. The method of moments is applied to obtain a system of simultaneous equations that can be solved for the unknown surface current densities and the interior electric field. The interior region partial differential equation and the exterior region surface integral equation are coupled in such a manner that many existing surface integral equation computer codes for treating problems involving scattering by homogeneous dielectric cylinders can be modified easily to generate the block of the matrix corresponding to the surface current interactions. The overall system matrix obtained using the method of moments is largely sparse. Numerical results are presented and compared with exact solutions for homogeneous and inhomogeneous circular cylinders     

Analysis of electromagnetic scattering from dielectrically coatedconducting cylinders using a multifilament current model

A method of moments solution is presented for the problem of transverse magnetic scattering from dielectrically coated conducting cylinders. The solution uses fictitious filamentary electric sources of yet unknown currents to simulate both the field scattered by the cylinder and the field inside the dielectric coating. The simulated fields obey the boundary conditions, namely, the continuity of the tangential components of the electric and magnetic fields across the air-dielectric interface and the vanishing of the tangential component of the electric field at the perfect conductor, at selected sets of points on these respective surfaces. The result is a matrix equation that is readily solved for the unknown current. The currents can be used to determine approximate values for the fields and field-related parameters of interest. The procedure is simple to implement and is general in that cylinders of smooth but otherwise arbitrary shape and coating of arbitrary complex permittivity can be handled. Illustrative examples are considered and compared with available data, demonstrating the efficiency of the solution     

Single integral equation for electromagnetic scattering bythree-dimensional homogeneous dielectric objects

A single integral equation formulation for electromagnetic scattering by three-dimensional (3-D) homogeneous dielectric objects is developed. In this formulation, a single effective electric current on the surface S of a dielectric object is used to generate the scattered fields in the interior region. The equivalent electric and magnetic currents for the exterior region are obtained by enforcing the continuity of the tangential fields across S. A single integral equation for the effective electric current is obtained by enforcing the vanishing of the total field due to the exterior equivalent currents inside S. The single integral equation is solved by the method of moments. Numerical results for a dielectric sphere obtained with this method are in good agreement with the exact results. Furthermore, the convergence speed of the iterative solution of the matrix equation in this formulation is significantly greater than that of the coupled integral equations formulation     

Surface integral formulation for calculating conductor anddielectric losses of various transmission structures

The power-loss method, along with a surface integral formulation, has been used to compute the attenuation constant in microstrip and coplanar structures. This method can be used for the analysis of both open and closed structures. Using the surface equivalence principle, the waveguide walls are replaced by equivalent electric surface currents and dielectric surfaces are replaced by equivalent electric and magnetic surface currents. Enforcing the appropriate boundary condition, and E-field integral equation (EFIE) is developed for these currents. Method of moments with pulse expansion and point matching testing procedure is used to transform the integral equation into a matrix one. The relationship between the propagation constant and frequency is found from the minimum eigenvalue of the moment matrix. The eigenvector pertaining to the minimum eigenvalue gives the unknown electric and magnetic surface currents     

Procedures for noninvasive electromagnetic property and dosimetrymeasurements

Five different formulations are presented that can be used to determine local values of the complex permittivity, electric field, polarization current density, and rate of energy deposition, within an inhomogeneous dielectric object of arbitrary shape. These formulations were obtained from the matrix equations for moment-method solutions of the electric field integral equation. Implementation of these formulations requires knowledge of the shape of the object and the incident electric field, and measurements of the scattered electric field at a number of points external to the object. The use of a new type of model, with antenna arrays for measurements, has caused a substantial improvement in matrix conditioning. Results of numerical simulations using 21-cell models with high dielectric contrast (bone/muscle) suggest that all five formulations have sufficient accuracy to warrant experimental testing using cylindrical scatterers with transverse magnetic polarization     

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.     

Combined field solution for TM scattering from multiple conducting and dielectric cylinders of arbitrary cross section

A simple moment solution is given for the problem of electromagnetic scattering from multiple conducting and dielectric cylinders of arbitrary cross section. The system of conducting and dielectric cylinders is excited by a plane-wave polarized transverse magnetic to the axis of the cylinders. The equivalence principle is used to obtain three coupled integral equations for the induced electric current on the conducting cylinders and the equivalent electric and magnetic currents on the surface of dielectric cylinders. The combined field integral equation (CFIE) formulation is used. Sample numerical results are presented. The agreement with available published data is excellent.     

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     

Fast convergent integral equation solution of strip gratings on dielectric substrate

A fast convergent integral equation solution to the scattering problem of a transverse electric/transverse magnetic (TE/TM) plane wave by a one-dimensional periodic array of thin metal strips on a dielectric substrate is described. The formulation of the integral equation is similar to that derived by Montgomery for a two-dimensional periodic array of thin conductors on a dielectric substrate. However, the basis functions which satisfy the appropriate edge conditions are incorporated here for the unknown current expansion on the strips. Following the standard Galerkin's procedure, one may readily determine the induced currents on the strips and thus the reflected and transmitted fields. Sample numerical results are given and good agreement with previously published data is obtained. It is found that the convergence rate of this method is improved by an order of magnitude. Also it is shown that the dielectric substrate has a strong effect on the scattering from the large spacing strip grating.     

On volume integral equations

The focus of this paper is on the volume integral representations to be used in constructing integral equations for composite volume media. The major thrust of the paper is to identify where derivatives of a discontinuous function arise in the derivation of the volume representation. Three different derivation methods are presented, resulting in identical representation independent of the derivation method. These representations agree with some in the existing literature and disagree with others. When an electric field formulation is considered, the source of disagreement manifests itself only when magnetic materials are present. Likewise, for the dual situation, the inconsistency appears for a magnetic field formulation of dielectric materials. This paper identifies the sources of error in the incorrect representations and its major contribution is the rigorously correct derivation of the representations to be used in volume integral equations. We also present numerical results for an integral equation derived from our representation. The numerical results employ only the E-field as the unknown and the singularity is handled in a manner analogous to a standard numerical treatment of the electric field integral equation.     

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.      

Numerical solution of integral equations for dielectric objects ofprismatic shapes

A numerical method to investigate scattering from dielectric geometries of prismatic shapes has been developed. The surface integral equations are formulated by Schelkunoff's equivalence principle in terms of equivalent surface electric and magnetic currents. To solve these integral equations for the unknown currents, the object's cross-section is mapped onto a circle. In the transformed space, Fourier type entire-domain basis functions are used in the cross section and triangular subdomain basis functions are selected along the generating curve to represent the currents. A moment method is then used to reduce the integral equations to a matrix equation to compute the current coefficients. It is found that the transformation of the object's surface to a circular shape improves the convergence of the current mode in the cross-section. However, the current modes are coupled on the surface and the matrix equation includes all the modes     

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     

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.     

RCS of two-dimensional structures consisting of both dielectricsand conductors of arbitrary cross section

The problem of determining the electromagnetic field scattered by two-dimensional structures consisting of both dielectric and conducting cylinders of arbitrary cross section is considered. The conductors may be in the form of strips and the dielectrics may be in the form of shells. The conductors may be partially or fully covered by dielectric layers, while the dielectrics may be partially covered by conductors. Only homogeneous dielectrics are studied. Both the transverse electric (TE) and the transverse magnetic (TM) cases are considered. The problem is formulated in terms of a set of coupled integral equations involving equivalent electric and magnetic surface currents radiating in unbounded media. The method of moments is used to solve the integral equations. Simple expansion and testing procedures are used. Numerical results for scattering cross sections are given for various structures     

Computation of cutoff wavenumbers of TE and TM modes in waveguidesof arbitrary cross sections using a surface integral formulation

A procedure is described for obtaining the cutoff wave numbers of transverse electric (TE) and transverse magnetic (TM) modes in waveguides of arbitrary cross section. A surface integral equation approach is used in which the E-field equation has been transformed into a matrix equation using the method of moments. An iterative technique is used to pick the eigenvalues of the solution matrix which corresponds to the waveguide cutoff wave numbers. The salient features of this technique are its speed, its simplicity, and the absence of any spurious modes when waveguides of arbitrary cross section are treated. The first four modes are tabulated for various waveguides, and the results are in very good agreement with published data     

应用互耦迭代法计算二维随机海面微波散射

海面对微波散射可用介质散射PMCHW方程描述,针对其离散后阻抗矩阵块的特点,引入电磁互耦迭代方法,结合多层快速多极子(MLFMA)求解,给出计算海面散射等效电流源和磁流源的算法.使用该方法求解分析得到不同起伏程度的海面对微波散射方向性的影响及规律,所得到的规律与文献报道实验测试结论相比较,对比结果验证方法的正确性.     

Green's function for arbitrarily shaped dielectric bodies of revolution

In this paper, a solution is developed to calculate the electric field at one point in space due to an electric dipole exciting an arbitrarily shaped dielectric body of revolution (BOR). Specifically, the electric field is determined from the solution of coupled surface integral equations (SIE) for the induced surface electric and magnetic currents on the dielectric body excited by an elementary electric current dipole source. Both the interior and exterior fields to the dielectric BOR may be accurately evaluated via this approach. For a highly lossy dielectric body, the numerical Green's function is also obtainable from an approximate integral equation (AIE) based on a surface boundary condition. If this equation is solved by the method of moments, significant numerical efficiency over SIE is realized. Numerical results obtained by both SIE and AIE approaches agree with the exact solution for the special case of a dielectric sphere. With this numerical Green's function, the complicated radiation and scattering problems in the presence of an arbitrarily shaped dielectric BOR are readily solvable by the method of moments.     

Transmission through dielectric-filled slots in a conductingcylindrical shell of arbitrary cross section: TE case

A simple moment solution is summarized for the problem of electromagnetic transmission through dielectric-filled slots in a conducting cylindrical shell of arbitrary cross section. The system is excited by a plane-wave polarized transverse electric (TE) to the axis of the shell. The equivalence principle is used to replace the shell and the dielectric by equivalent electric and magnetic surface currents radiating into an unbounded medium. Two different sets of coupled integral equations involving the surface currents are obtained by enforcing the boundary conditions on the tangential components of the total electric and magnetic fields. The method of moments is used to solve the integral equations. Pulses are used for both expansion and testing functions. Special attention is paid to circular and rectangular shells. Results for shell surface current, the internal field, and the aperture field are presented. For the case of air dielectric filling, the results computed using the electric field and/or the magnetic field formulation are in very good agreement with published data. In general, it is observed that the effect of filling a slot with a dielectric is not predictable from a simple theory     

On the determination of resonant modes of dielectric objects using surface Integral equations

Although surface integral equations have been extensively used for solving the scattering problem of arbitrarily shaped dielectric objects, when applied to the resonance problem, there are still some issues not fully addressed by the literature. In this paper, the method of moments with Rao-Wilton-Glisson basis functions is applied to the electric field integral equation (EFIE) for solving the resonance problem of dielectric objects. The resonant frequency is obtained by searching for the minimum of the reciprocal of the condition number of the impedance matrix in the complex frequency plane, and the modal field distribution is obtained through singular value decomposition (SVD). The determinant of the impedance matrix is not used since it is difficult to find its roots. For the exterior EFIE, the original basis functions are used as testing functions; for the interior EFIE, the basis functions rotated by 90/spl deg/ are used as testing functions. To obtain an accurate modal field solution, the impedance matrix needs to be reduced by half before SVD is applied to it. Numerical results are given and compared with those obtained by using the volume integral equation.