Sensitivity Analysis of Electromagnetic Scattering from Dielectric Targets with Polynomial Chaos Expansion and Method of Moments

In this paper,an adaptive polynomial chaos expansion method(PCE)based on the method of moments(MoM)is proposed to construct surrogate models for electromagnetic scattering and further sensitivity analysis.The MoM is applied to accurately solve the electric field integral equation(EFIE)of electromagnetic scattering from homogeneous dielectric targets.Within the bistatic radar cross section(RCS)as the research object,the adaptive PCE algorithm is devoted to selecting the appropriate order to construct the multivariate surrogate model.The corresponding sensitivity results are given by the further derivative operation,which is compared with those of the finite difference method(FDM).Several examples are provided to demonstrate the effectiveness of the proposed algorithm for sensitivity analysis of electromagnetic scattering from homogeneous dielectric targets.

EVALUATION OF ELECTROMAGNETIC SCATTERING OF TWO DIMENSIONAL CONDUCTING OBJECTS--TE CASE

A hybrid technique is developed for the evaluation of two dimensional electromagnetic scattering from electrically large conducting bodies with cracks on their surfaces (TE case). The edge based finite element method (FEM) is employed to compute the scattering from the cracks. Physical optics (PO) and physical theory of diffraction (PTD) are utilized to evaluate the scattering from the large bodies with the cracks filled with perfect conductors. These two methods are combined by an efficient coupling scheme. Some of numerical results are presented. It is shown that the hybrid technique has some advantages over other methods in regard to saving computer memory units and CPU time.

Electromagnetic scattering from two parallel 2D targets arbitrarily located in a Gaussian beam

In this paper based on the equivalence principle and the reciprocity theorem, the scattered field up to second-order by two parallel 2D targets arbitrarily located in a Gaussian beam is considered. The first-order solution can easily be obtained by calculating the scattered field from isolated targets when illuminated by a Gaussian beam. However, because of the difficulty in formulating the couple scattering field, it is almost impossible to find an analytical solution for the second-order scattered field if the shapes of 2D targets are not canonical geometries. In order to overcome this problem, in this paper, the second-order solution is derived by using the technique based on the reciprocity theorem and the equivalence principle. Meanwhile, the relation between the secondary scattered field from target #1 and target #2 is obtained. Specifically, the bi- and mono-static scattering of Gaussian beam by two parallel adjacent inhomogeneous plasma-coated conducting circular cylinders is calculated and the dependence of attenuation of the scattering width on the thickness of the coated layer, electron number density, collision frequency and radar frequency is discussed in detail.

The scattering fields for a spherical target irradiated by a plane electromagnetic wave in an arbitrary direction

The relation between corresponding trigonometric functions in two rotating coordinate systems is presented. The transformation formula for a vector in the two rotating spherical coordinate systems is obtained. The scattering fields for a spherical target irradiated by a plane electromagnetic wave in an arbitrary direction are derived. These fields in a particular case retrogress to those available in the literature. The obtained results have great potential in practical applications.

Polarization ratio characteristics of electromagnetic scattering fromsea ice inpolar areas

In the global climate system, the polar regions are sensitive indicators of climate change, in which sea ice plays an important role. Satellite remote sensing is a significant tool for monitoring sea ice. The use of synthetic aperture radar(SAR) images to distinguish sea ice from sea water is one of the current research hotspots in this topic. To distinguish sea ice from the open sea, the polarization ratio characteristics of sea ice and sea water are studied for L-band and C-band radars, based on an electromagnetic scattering model of sea ice derived from the integral equation method(IEM) and the radiative transfer(RT) model. Numerical experiments are carried out based on the model and the results are given as follows. For L-band, the polarization ratio for sea water depends only on the incident angle, while the polarization ratio for sea ice is related to the incident angle and the ice thickness. For C-band, the sea water polarization ratio is influenced by the incident angle and the root mean square(RMS) height of the sea surface. For C-band, for small to medium incident angles,the polarization ratio for bare sea ice is mainly determined by the incident angle and ice thickness. When the incident angle increases, the RMS height will also affect the polarization ratio for bare sea ice. If snow covers the sea ice, then the polarization ratio for sea ice decreases and is affected by the RMS height of snow surface, snow thickness, volume fraction and the radius of scatterers. The results show that the sea ice and the open sea can be distinguished by using either L-band or C-band radar according to their polarization ratio difference. However, the ability of L-band to make this differentiation is higher than that of C-band.

Fast and accurate analysis of electromagnetic scattering from targets situated in dielectric half space

Electromagnetic scattering from targets situated in half space is solved by applying fast inhomogeneous plane wave algorithm combined with a tabulation and interpolation method. The integral equation is set up based on derivation of dyadic Green＇s functions in this environment. The coupling is divided into nearby region and well-separated region by grouping. The Green＇s function can be divided into two parts： primary term and reflected term. In the well-separated region, the two terms are both expressed as Sommerfeld integral, which can be accelerated by deforming integral path and taking interpolation and extrapolation. For the nearby region, the direct Sommerfeld integral makes the filling of impedance matrix time-expensive. A tabulation and interpolation method is applied to speed up this process. This infinite integral is pre-computed in sampling region, and a two-dimensional table is then set up. The impedance elements can then be obtained by interpolation. Numerical results demonstrate the accuracy and efficiency of this algorithm.

Analytical solution for electromagnetic scattering from a sphere of uniaxial left-handed material

Based on the analytical solution of electromagnetic scattering by a uniaxial anisotropic sphere in the spectral domain, an analytical solution to the electromagnetic scattering by a uniaxial left-handed materials (LHMs) sphere is obtained in terms of spherical vector wave functions in a uniaxial anisotropic LHM medium. The expression of the analytical solution contains only some one-dimensional integral which can be calculated easily. Numerical results show that Mie series of plane wave scattering by an isotropic LHM sphere is a special case of the present method. Some numerical results of electromagnetic scattering of a uniaxial anisotropic sphere by a plane wave are given.

A Second-Order Method for the Electromagnetic Scattering from a Large Cavity

In this paper,we study the electromagnetic scattering from a two dimen- sional large rectangular open cavity embedded in an infinite ground plane,which is modelled by Helmholtz equations.By introducing nonlocal transparent boundary con- ditions,the problem in the open cavity is reduced to a bounded domain problem.A hypersingular integral operator and a weakly singular integral operator are involved in the TM and TE cases,respectively.A new second-order Toeplitz type approximation and a second-order finite difference scheme are proposed for approximating the hyper- singular integral operator on the aperture and the Helmholtz in the cavity,respectively. The existence and uniqueness of the numerical solution in the TE case are established for arbitrary wavenumbers.A fast algorithm for the second-order approximation is pro- posed for solving the cavity model with layered media.Numerical results show the second-order accuracy and efficiency of the fast algorithm.More important is that the algorithm is easy to implement as a preconditioner for cavity models with more general media.

Electromagnetic scattering from two-layered rough interfaces with a PEC object:vertical polarization

Electromagnetic （EM） scattering from a stack of two rough interfaces separating a homogeneous medium with a perfectly electric conducting （PEC） object has been calculated through the method of moments for vertical polarization. Theoretical formulations of EM scattering from multi-layered rough interfaces with a PEC object have been derived in detail and the total fields and their normal derivatives on the rough interfaces are solved. The two-layered model is a special case. In this work, a Gaussian rough surface was applied to simulate the rough interface. A cylinder was located above, between or below the two-layered rough interfaces. Through numerical simulations, the validity of this work is demonstrated by comparing it with existing scattering models, which are special cases that include a PEC object located above/below a single-layered rough interface and two-layered rough interfaces without an object. Subsequently, the influences of characteristic parameters, such as the relative permittivity of the medium, as well as the average height between the two rough surfaces, on the bistatic scattering coefficient are discussed.

Research on the Electromagnetic Scattering Near the Field

We calculate and analyze the scattering near the field from some simple and complex targets using the method of picture elements (PEL), based upon the method of high-frequency approximation. It introduces the critical distance of the near field and the far field which is related with the dimension of the target. The problem of the EMS near field from large size objects can be transformed to the problem of the far field by parting it to many very small size elements. By calculating the EMS near fields of some simple and complex targets based on the SCTE (scattering from complex targets and environments) system, the results show that there are much difference between the near field and the far field. And the characteristics of the near field are more complicated. This work has practical engineering value in the area of the electromagnetic compatibility (EMC), electromagnetic interference (EMI) prediction and electromagnetic scattering (EMS).

Numerical Method for Solving Electromagnetic Scattering Problem by Many Small Impedance Bodies

In this paper, we study electromagnetic (EM) wave scattering problem by many small impedance bodies. A numerical method for solving this problem is presented. The problem is solved under the physical assumptions ka??1, where a is the characteristic size of the bodies and k is the wave number. This problem is solved asymptotically and numerical experiments are provided to illustrate the idea of the method. Error estimate for the asymptotic solution is also discussed.

Numerical Method for Solving Electromagnetic Wave Scattering by One and Many Small Perfectly Conducting Bodies

In this paper, we investigate the problem of electromagnetic (EM) wave scattering by one and many small perfectly conducting bodies and present a numerical method for solving it. For the case of one body, the problem is solved for a body of arbitrary shape, using the corresponding boundary integral equation. For the case of many bodies, the problem is solved asymptotically under the physical assumptions a d a is the characteristic size of the bodies, d is the minimal distance between neighboring bodies, λ = 2π/k is the wave length and k is the wave number. Numerical results for the cases of one and many small bodies are presented. Error analysis for the numerical method is also provided.

Finite Volume Time Domain with the Green Function Method for Electromagnetic Scattering in Schwarzschild Spacetime

The finite volume time domain(FVTD) algorithm and Green function algorithm are extended to Schwarzschild spacetime for numerical simulation of electromagnetic scattering. The FVTD method in Schwarzschild spacetime is developed by filling the flat spacetime with an equivalent medium. The Green function in Schwarzschild spacetime is acquired by solving initial value problems. Both the FVTD code and the Green function code are validated by numerical results. Scattering in Schwarzschild spacetime is simulated with these methods.

Analysis of the Linear Sampling Method for Boundary Reconstruction of Chiral Obstacle in Electromagnetic Scattering Problems

In this paper, we consider the inverse scattering by chiral obstacle in electromagnetic fields, and prove that the linear sampling method is also effective to determine the support of a chiral obstacle from the noisy far field data.

Efficient Solution to Electromagnetic Scattering Problems of Bodies of Revolution by Compressive Sensing

Under the theory structure of compressive sensing （CS）, an underdetermined equation is deduced for describing the discrete solution of the electromagnetic integral equation of body of revolution （BOR）, which will result in a small-scale impedance matrix. In the new linear equation system, the small-scale impedance matrix can be regarded as the measurement matrix in CS, while the excited vector is the measurement of unknown currents. Instead of solving dense full rank matrix equations by the iterative method, with suitable sparse representation, for unknown currents on the surface of BOR, the entire current can be accurately obtained by reconstructed algorithms in CS for small-scale undetermined equations. Numerical results show that the proposed method can greatly improve the computgtional efficiency and can decrease memory consumed.

Electromagnetic backscattering from one-dimensional drifting fractal sea surface Ⅱ：Electromagnetic backscattering model

Sea surface current has a significant influence on electromagnetic（EM） backscattering signals and may constitute a dominant synthetic aperture radar（SAR） imaging mechanism. An effective EM backscattering model for a one-dimensional drifting fractal sea surface is presented in this paper. This model is used to simulate EM backscattering signals from the drifting sea surface. Numerical results show that ocean currents have a significant influence on EM backscattering signals from the sea surface. The normalized radar cross section（NRCS） discrepancies between the model for a coupled wavecurrent fractal sea surface and the model for an uncoupled fractal sea surface increase with the increase of incidence angle,as well as with increasing ocean currents. Ocean currents that are parallel to the direction of the wave can weaken the EM backscattering signal intensity, while the EM backscattering signal is intensified by ocean currents propagating oppositely to the wave direction. The model presented in this paper can be used to study the SAR imaging mechanism for a drifting sea surface.

Overview: Electromagnetic Scattering from Ocean Surface

Understanding the sea surface scattering process is very important in the development of models to detect the target above or under the surface. In this paper, both the analytical and the numerical methods applied in sea surface scattering are summarized. Some important problems concerned in this field are discussed. For numerical study, edge effect brings artificial nonrealistic scattering and therefore must be suppressed. Different edge treatment methods are compared in this paper. Scattering of breaking wave surface at very low grazing angle always needs more attentions than other scattering problems. Some numerical results show the existence of the special phenomena at very low grazing angle, for example, the ＂sea spikes＂ and the Doppler splitting.

Reducing the calculation workload of the Green function for electromagnetic scattering in a Schwarzschild gravitational field

When the finite difference time domain(FDTD) method is used to solve electromagnetic scattering problems in Schwarzschild space-time, the Green functions linking source/observer to every surface element on connection/output boundary must be calculated.When the scatterer is electrically extended, a huge amount of calculation is required due to a large number of surface elements on the connection/output boundary.In this paper, a method for reducing the calculation workload of Green function is proposed.The Taylor approximation is applied for the calculation of Green function.New transport equations are deduced.The numerical results verify the effectiveness of this method.

ELECTROMAGNETIC SCATTERING FROM A LOSSY DIELECTRIC CYLINDER WITH ARBITRARY CROSS SECTION

In this paper,the electromagnetic scattering of plane wave with TM and TEpolarizations from a lossy dielectric cylinder with arbitrary cross section is analyzed andcomputed by using the boundary element method.The numerical results of scattering pat-tern for a circular cylinder completely coincide with the analytical solution and that of thescattering properties in the near-or far-field region for elliptical and regular triangularcylinders are presented,which point out some valuable conclusions.

Scattering of Electromagnetic Waves by Drift Vortex in Plasma

In a quasi-two-dimensional model, the scattering of incident ordinary electromag- netic waves by a dipole-electrostatic drift vortex is studied with first-order Born approximation. The distribution of the scattering cross-section and total cross-section are evaluated analytically in different approximate conditions, and the physical interpretations are discussed. When the wavelength of incident wave is much longer than the vortex radius （kia〈〈1）, it is found that the angle at which the scattering cross-section reaches its maxim depends significantly on the approximation of the parameters of the vortex used. It is also found that the total scattering cross-section has an affinitive relation with the parameters of the plasma, while it is irrelevant to the frequency of the incident wave in a wide range of parameters of the vortex. In a totally different range of parameters when incident wave is in the radar-frequency range （then kia 〈〈 1, the wavelength of incident wave is much shorter than the vortex radius）, the numerical procedure is conducted with computer in order to obtain the distribution and the total expression of the scattering crosssection. Then it is found that the total scattering cross-section in the low frequency range is much larger than that in high frequency range, so the scattering is more effective in the low frequency range than in high frequency range.