An overfilled cavity problem for Maxwell's equations

This paper is concerned with the mathematical analysis of the scattering of a time‐harmonic electromagnetic plane wave by an open and overfilled cavity that is embedded in a perfect electrically conducting infinite ground plane, where the electromagnetic wave propagation is governed by the Maxwell equations. Above the flat ground surface and the open aperture of the cavity, the space is assumed to be filled with a homogeneous medium with a constant permittivity and permeability, whereas the interior of the cavity is filled with some inhomogeneous medium with a variable permittivity and permeability. The scattering problem is modeled as a boundary value problem over a bounded domain, with transparent boundary condition proposed on the hemisphere enclosing the inhomogeneity represented by the cavity. The existence and uniqueness of the weak solution for the model problem are established by using a variational approach. The perfectly matched layer (PML) method is investigated to truncate the unbounded electromagnetic cavity scattering problem. It is shown that the truncated PML problem attains a unique solution. An explicit error estimate is given between the solution of the original scattering problem and that of the truncated PML problem. The error estimate implies that the PML solution converges exponentially to the original cavity scattering problem by increasing either the PML medium parameter or the PML layer thickness. The convergence result is expected to be useful for determining the PML medium parameter in the computational electromagnetic scattering problem. Copyright © 2012 John Wiley & Sons, Ltd.     

Low-Cost Numerical Method for Solving a Coefficient Inverse Problem for the Wave Equation in Three-Dimensional Space

For the acoustic-sensing problem of determining the characteristics of a local inhomogeneity scattering a wave field in three-dimensional space, a numerical algorithm is proposed and justified that is efficient in terms of computational resources and CPU time. The algorithm is based on the fast Fourier transform, which is used under certain a priori assumptions on the character of the inhomogeneity and the observation domain of the scattered field. Typical numerical results obtained by solving this inverse problem with simulated data on a personal computer are presented, which demonstrate the capabilities of the algorithm.     

Preconditioning techniques for the iterative solution of scattering problems

We consider a time-harmonic electromagnetic scattering problem for an inhomogeneous medium. Some symmetry hypotheses on the refractive index of the medium and on the electromagnetic fields allow to reduce this problem to a two-dimensional scattering problem. This boundary value problem is defined on an unbounded domain, so its numerical solution cannot be obtained by a straightforward application of usual methods, such as for example finite difference methods, and finite element methods. A possible way to overcome this difficulty is given by an equivalent integral formulation of this problem, where the scattered field can be computed from the solution of a Fredholm integral equation of second kind. The numerical approximation of this problem usually produces large dense linear systems. We consider usual iterative methods for the solution of such linear systems, and we study some preconditioning techniques to improve the efficiency of these methods. We show some numerical results obtained with two well known Krylov subspace methods, i.e., Bi-CGSTAB and GMRES.     

Recent developments in the qualitative approach to inverse electromagnetic scattering theory

We consider the inverse scattering problem of determining both the shape and some of the physical properties of the scattering object from a knowledge of the (measured) electric and magnetic fields due to the scattering of an incident time-harmonic electromagnetic wave at fixed frequency. We shall discuss the linear sampling method for solving the inverse scattering problem which does not require any a priori knowledge of the geometry and the physical properties of the scatterer. Included in our discussion is the case of partially coated objects and inhomogeneous background. We give references for numerical examples for each problem discussed in this paper.     

Method of Integral Equations for the Three-Dimensional Problem of Wave Reflection from an Irregular Surface

The problem on the reflection of the field of a plane H-polarized three-dimensional electromagnetic wave from a perfectly conducting interface between media which contains a local perfectly conducting inhomogeneity is considered. To construct a numerical algorithm, the boundary value problem for the system of Maxwell equations in an infinite domain with irregular boundary is reduced to a system of singular integral equations, which is solved by the approximation–collocation method. The elements of the resulting complex matrix are calculated by a specially developed algorithm. The solution of the system of singular integral equations is used to obtain an integral representation for the reflected electromagnetic field and computational formulas for the directional diagram of the reflected electromagnetic field in the far region.     

Surface wave scattering by a rectangular impedance cylinder located on a reactive plane

The electromagnetic scattering of the surface wave by a rectangular impedance cylinder located on an infinite reactive plane is considered for the case that the impedances of the horizontal and vertical sides of the cylinder can have different values. Firstly, the diffraction problem is reduced into a modified Wiener–Hopf equation of the third kind and then solved approximately. The solution contains branch‐cut integrals and two infinite sets of constants satisfying two infinite systems of linear algebraic equations. The approximate analytical or numerical evaluations of corresponding integrals and numerical solution of the linear algebraic equation systems are obtained for various values of parameters such as the surface reactance of the plane, the vertical and horizontal wall impedances, the width and the height of the cylinder. Copyright © 2004 John Wiley & Sons, Ltd.     

On the Two-dimensional Inverse Scattering Problems in Electromagnetics

In recent years, inverse scattering problems have received much attention because of their important applications. Given the incident and scattered waves, an inverse scattering problem in general is to determine the properties of the scatterer. In radar or sonar a known incident wave and observed scattered wave are used to detect the properties and the presence of aircraft or submarine objects; in MRI scanning, tomography X-rays and ultrasound, scattered waves are used to determine the presence or properties of tumors by detecting density variations, to name a few. In this article, we are concerned with the two-dimensional electromagnetic inverse scattering problem. An iterative algorithm for the transversal electric waves will be given based on a singular domain integral equation formulation. Basic features of a scattering object such as shape, location and index of refraction will be recovered from measurements of the field scattered by the object (when illuminated by electromagnetic waves with the magnetic vector polarized along the cylinder axis). Some numerical experiments are included to illustrate the efficiency of the algorithm.     

Efficient Spectral Methods for Transmission Eigenvalues and Estimation of the Index of Refraction

An important step in estimating the index of refraction of electromagnetic scattering problems is to compute the associated transmission eigenvalue problem. We develop in this paper efficient and accurate spectral methods for computing the transmission eigenvalues associated to the electromagnetic scattering problems. We present ample numerical results to show that our methods are very effective for computing transmission eigenvalues (particularly for computing the smallest eigenvalue), and together with the linear sampling method, provide an efficient way to estimate the index of refraction of a non-absorbing inhomogeneous medium.     

关于电磁散射问题的一类边界元方法

Problems for electromagnetic scattering are of significant importance in many areas of technology. In this paper we discuss the scattering problem of electromagnetic wave incident by using boundary element method associated with splines. The problem is modelled by a boundary value problem for the Helmholtz eouation     

Numerical studies of time-independent and time-dependent scattering by several elliptical cylinders

A numerical solution to the problem of time-dependent scattering by an array of elliptical cylinders with parallel axes is presented. The solution is an exact one, based on the separation-of-variables technique in the elliptical coordinate system, the addition theorem for Mathieu functions, and numerical integration. Time-independent solutions are described by a system of linear equations of infinite order which are truncated for numerical computations. Time-dependent solutions are obtained by numerical integration involving a large number of these solutions. First results of a software package generating these solutions are presented: wave propagation around three impenetrable elliptical scatterers. As far as we know, this method described has never been used for time-dependent multiple scattering.     

On the far field patterns for electromagnetic scattering by a chiral obstacle in a chiral environment

A time-harmonic plane electromagnetic wave is scattered by a chiral body in a chiral environment. The body is either a perfect conductor, or a dielectric, or a scatterer with an impedance surface. Using the Huygens's principle, we construct in closed forms both the left-circularly polarized and right-circularly polarized electric far field patterns for such chiral media. We prove reciprocity relations and general scattering theorems for chiral materials which are a generalization of those obtained by Twersky for achiral electromagnetic scattering. In the special case when the directions of incidence and observation are the same we prove the associated forward scattering theorems.     

Numerical solution of electromagnetic scattering from a large partly covered cavity

The paper focuses on the numerical study of electromagnetic scattering from two-dimensional (2D) large partly covered cavities, which is described by the Helmholtz equation with a nonlocal boundary condition on the aperture. The classical five-point finite difference method is applied for the discretization of the Helmholtz equation and a linear approximation is used for the nonlocal boundary condition. We prove the existence and uniqueness of the numerical solution when the medium in the cavity is y-direction layered or the number of the mesh points on the aperture is large enough. The fast algorithm proposed in Bao and Sun (2005) [2] for open cavity models is extended to solving the partly covered cavity problem with (vertically) layered media. A preconditioned Krylov subspace method is proposed to solve the partly covered cavity problem with a general medium, in which a layered medium model is used as a preconditioner of the general model. Numerical results for several types of partly covered cavities with different wave numbers are reported and compared with those by ILU-type preconditioning algorithms. Our numerical experiments show that the proposed preconditioning algorithm is more efficient for partly covered cavity problems, particularly with large wave numbers.     

Incremental unknowns preconditioning for solving the Helmholtz equation

An efficient preconditioner is developed for solving the Helmholtz problem in both high and low frequency (wavenumber) regimes. The preconditioner is based on hierarchical unknowns on nested grids, known as incremental unknowns (IU). The motivation for the IU preconditioner is provided by an eigenvalue analysis of a simplified Helmholtz problem. The performance of our preconditioner is tested on the iterative solution of two‐dimensional electromagnetic scattering problems. When compared with other well‐known methods, our technique is shown to often provide a better numerical efficacy and, most importantly, to be more robust. Moreover, for the best performance, the number of IU levels used in the preconditioner should be designed for the coarsest grid to have roughly two points per linear wavelength. This result is consistent with the conventional sampling criteria for wave phenomena in contrast with existing IU applications for solving the Laplace/Poisson problem, where the coarsest grid comprises just one interior point. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Scattering of electromagnetic plane waves by a spheroid uniformly moving in free space

We investigate the scattering process, generated by a plane electromagnetic field that is incident upon a moving perfectly conducting spheroid. An accurate treatment of the electromagnetic waves interaction with scatterers in uniform motion is based on the special relativity principle. In the object's frame the incident wave is assumed to have a wavelength which is much larger than the characteristic dimension of the scatterer and thus the low‐frequency approximation method is applicable to the scattering problem. For the near electromagnetic field we obtain the zeroth‐order low‐frequency coefficients, while in the far field we calculate the leading terms for the scattering amplitude and scattering cross‐section. Finally, using the inverse Lorentz transform, we obtain the same approximations in the observer's frame. Copyright © 2006 John Wiley & Sons, Ltd.     

Combined far‐ield operators in electromagnetic inverse scattering theory

We consider the inverse scattering problem of determining the shape of a perfect conductor D from a knowledge of the scattered electromagnetic wave generated by a time‐harmonic plane wave incident upon D. By using polarization effects we establish the validity of the linear sampling method for solving this problem that is valid for all positive values of the wave number. We also show that it suffices to consider incident directions and observation angles that are restricted to a limited aperture. Copyright © 2003 John Wiley & Sons, Ltd.     

周期手性结构电磁散射问题的积分方程方法

In this paper, we consider the electromagnetic scattering by a periodic chiral structure. The media is homogeneous and the structure is periodic in one direction and invariant in another direction. The electromagnetic fields inside the chiral medium are governed by Maxwell equations together with the Drude-BornFedorov equations. We simplify the problem to a two-dimensional scattering problem and discuss the existence and the uniqueness of solutions by an integral equation approach. We show that for all but possibly a discrete set of wave numbers, the integral equation has a unique solution.     

Analysis of Direct and Inverse Cavity Scattering Problems

Consider the scattering of a time-harmonic electromagnetic plane wave by an arbitrarily shaped and filled cavity embedded in a perfect electrically conducting infinite ground plane.A method of symmetric coupling of finite element and boundary integral equations is presented for the solutions of electromagnetic scattering in both transverse electric and magnetic polarization cases.Given the incident field,the direct problem is to determine the field distribution from the known shape of the cavity; while the inverse problem is to determine the shape of the cavity from the measurement of the field on an artificial boundary enclosing the cavity.In this paper,both the direct and inverse scattering problems are discussed based on a symmetric coupling method.Variational formulations for the direct scattering problem are presented,existence and uniqueness of weak solutions are studied,and the domain derivatives of the field with respect to the cavity shape are derived.Uniqueness and local stability results are established in terms of the inverse problem.     

Numerical solution of integral equations of first and second kind in scalar and vector problems of diffraction on a cube

We consider the problem of numerical simulation of the scattering of acoustic and electromagnetic waves on a cube whose edge ha s length up to 8 wave lengths of the incident wave. We describe a scheme using a representation of the boundary integral equation in the form of an operator convolution equation on the symmetry group of the cube. We compare the results of numerical solution of integral equations of first and second kind for scalar and vector problems of diffraction of a plane wave on a cube. Translated fromProblemy Matematicheskoi Fiziki, 1998, pp. 36–45.     

A stable well-conditioned integral equation for electromagnetism scattering

Finding a formulation for electromagnetic scattering of surfaces which is both well-posed and produces a well-conditioned linear system is still a challenging problem. We here propose one such formulation valid in the high-frequency regime. The mathematical analysis is provided and numerical results on rather complex geometries show the performance of the method.     

Diagonal forms of the translation operators in the fast multipole algorithm for scattering problems

The integral equations of acoustic and electromagnetic scattering generate large dense systems of linear equations. These systems are efficiently solved with iterative methods where the matrix-vector multiplication is computed using a special fast method, such as the fast Fourier transform or the fast multipole method (FMM). In this paper, the so called diagonal forms of the translation operators for the fast multipole method are derived starting from integral representations of certain special functions. Error analysis of the FMM is given, considering both the truncation error of potential expansions and the errors from the use of numerical integration in the diagonal translation theorem. The implications of the error bounds on the FMM algorithm are discussed.This work has been financially supported by the Jenny and Antti Wihuri Foundation and by the Cultural Foundation of Finland.