The efficient solution of electromagnetic scattering for inhomogeneous media

We consider an electromagnetic scattering problem for inhomogeneous media. In particular, we focus on the numerical computation of the electromagnetic scattered wave generated by the interaction of an electromagnetic plane wave and an inhomogeneity in the corresponding propagation medium. This problem is studied in the VV polarization case, where some special symmetry requirements for the incident wave and for the inhomogeneity are assumed. This problem is reformulated as a Fredholm integral equation of second kind, which is discretized by a linear system having a special form. This allows to compute efficiently an approximate solution of the scattering problem by using iterative techniques for linear systems. Some numerical examples are reported.     

A Galerkin approximation for integro-differential equations in electromagnetic scattering from a chiral medium

We consider the scattering of time-harmonic electromagnetic waves from a chiral medium. It is known for the Drude–Born–Fedorov model that the forward scattering problem can be described by an integro-differential equation. In this paper we study a Galerkin finite element approximation for this integro-differential equation. Our Galerkin scheme, which relies on a suitable periodization of the integral equation, enables the use of the fast Fourier transform and a simple numerical implementation. We establish a quasi-optimal convergence analysis for the Galerkin method. Explicit formulas for the discrete scheme are also provided.     

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

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.     

A preconditioned GMRES for complex dense linear systems from electromagnetic wave scattering problems

We consider the numerical solution of linear systems arising from the discretization of the electric field integral equation (EFIE). For some geometries the associated matrix can be poorly conditioned making the use of a preconditioner mandatory to obtain convergence. The electromagnetic scattering problem is here solved by means of a preconditioned GMRES in the context of the multilevel fast multipole method (MLFMM). The novelty of this work is the construction of an approximate hierarchically semiseparable (HSS) representation of the near-field matrix, the part of the matrix capturing interactions among nearby groups in the MLFMM, as preconditioner for the GMRES iterations. As experience shows, the efficiency of an ILU preconditioning for such systems essentially depends on a sufficient fill-in, which apparently sacrifices the sparsity of the near-field matrix. In the light of this experience we propose a multilevel near-field matrix and its corresponding HSS representation as a hierarchical preconditioner in order to substantially reduce the number of iterations in the solution of the resulting system of equations.     

Analysis of a finite element PML approximation for the three dimensional time-harmonic Maxwell problem

In our paper [Math. Comp. 76, 2007, 597-614] we considered the acoustic and electromagnetic scattering problems in three spatial dimensions. In particular, we studied a perfectly matched layer (PML) approximation to an electromagnetic scattering problem. We demonstrated both the solvability of the continuous PML approximations and the exponential convergence of the resulting solution to the solution of the original acoustic or electromagnetic problem as the layer increased.

In this paper, we consider finite element approximation of the truncated PML electromagnetic scattering problem. Specifically, we consider approximations which result from the use of Nédélec (edge) finite elements. We show that the resulting finite element problem is stable and gives rise to quasi-optimal convergence when the mesh size is sufficiently small.

     

The Two-dimensional Electromagnetic Scattering from Periodic Chiral Structures and Its Finite Element Approximation

In this paper, we consider the electromagnetic scattering from periodic chiral structures. The structure is periodic in one direction and invariant in another direction. The electromagnetic fields in the chiral medium are governed by the Maxwell equations together with the Drude-Born-Fedorov equations. We simplify the problem to a two-dimensional scattering problem and we show that for all but possibly a discrete set of wave numbers, there is a unique quasi-periodic weak solution to the diffraction problem. The diffraction problem can be solved by finite element method. We also establish uniform error estimates for the finite element method and the error estimates when the truncation of the nonlocal transparent boundary operators takes place.     

A boundary element method for the Dirichlet eigenvalue problem of the Laplace operator

The solution of eigenvalue problems for partial differential operators by using boundary integral equation methods usually involves some Newton potentials which may be resolved by using a multiple reciprocity approach. Here we propose an alternative approach which is in some sense equivalent to the above. Instead of a linear eigenvalue problem for the partial differential operator we consider a nonlinear eigenvalue problem for an associated boundary integral operator. This nonlinear eigenvalue problem can be solved by using some appropriate iterative scheme, here we will consider a Newton scheme. We will discuss the convergence and the boundary element discretization of this algorithm, and give some numerical results.     

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.     

On the convergence of a numerical method for a hypersingular integral equation on a closed surface

We consider a linear integral equation with a hypersingular integral treated in the sense of the Hadamard finite value. This equation arises in the solution of the Neumann boundary value problem for the Laplace equation with a representation of a solution in the form of a double-layer potential. We consider the case in which the interior or exterior boundary value problem is solved in a domain; whose boundary is a smooth closed surface, and an integral equation is written out on that surface. For the integral operator in that equation, we suggest quadrature formulas like the method of vortical frames with a regularization, which provides its approximation on the entire surface for the use of a nonstructured partition. We construct a numerical scheme for the integral equation on the basis of suggested quadrature formulas, prove an estimate for the norm of the inverse matrix of the related system of linear equations and the uniform convergence of numerical solutions to the exact solution of the hypersingular integral equation on the grid.     

Numerical solution of 3D problems of electromagnetic wave diffraction on a system of ideally conducting surfaces by the method of hypersingular integral equations

We construct a numerical method for solving problems of electromagnetic wave diffraction on a system of solid and thin objects based on the reduction of the problem to a boundary integral equation treated in the sense of the Hadamard finite value. For the construction of such an equation, we construct a numerical scheme on the basis of the method of piecewise continuous approximations and collocations. Unlike earlier known schemes, by using the below-suggested scheme, we have found approximate analytic expressions for the coefficients of the arising system of linear equations by isolating the leading part of the kernel of the integral operator. We present examples of solution of a number of model problems of the diffraction of electromagnetic waves by the suggested method.     

Convergence of a numerical scheme of the type of the vortex loop method on a closed surface with an approximation to the surface shape

We consider a linear integral equation with a hypersingular integral treated in the sense of the Hadamard finite value. This equation arises when solving the Neumann boundary value problem for the Laplace equation with the use of the representation of the solution in the form of a double layer potential. We study the case in which an exterior or interior boundary value problem is solved in a domain whose boundary is a smooth closed surface and the integral equation is written out on that surface. For the numerical solution of the integral equation, the surface is approximated by spatial polygons whose vertices lie on the surface. We construct a numerical scheme for solving the integral equation on the basis of such an approximation to the surface with the use of quadrature formulas of the type of the method of discrete singularities with regularization. We prove that the numerical solutions converge to the exact solution of the hypersingular integral equation uniformly on the grid.     

On the convergence of the vortex loop method with regularization for the Neumann boundary value problem on a plane screen

We consider a linear integral equation with a supersingular integral treated in the sense of the Hadamard finite value, which arises in the solution of the Neumann boundary value problem for the Laplace equation with the representation of the solution in the form of a doublelayer potential. We consider the case in which the exterior boundary value problem is solved outside a plane surface (a screen). For the integral operator in the above-mentioned equation, we suggest quadrature formulas of the vortex loop method with regularization, which provide its approximation on the entire surface when using an unstructured partition. In the problem in question, the derivative of the unknown density of the double-layer potential, as well as the errors of quadrature formulas, has singularities in a neighborhood of the screen edge. We construct a numerical scheme for the integral equation on the basis of the suggested quadrature formulas and prove an estimate for the norm of the inverse matrix of the resulting system of linear equations and the uniform convergence of the numerical solutions to the exact solution of the supersingular integral equation on the grid.     

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.     

Integral equations in the wave scattering problem at an irregular interface of domains

Modeling the scattering of electromagnetic waves at an interface of media with different characteristics, one encounters the conjugation problem. Using the method of boundary integral equations and the theory of generalized potentials, we prove the classical resolvability of this problem. The boundary is assumed to be irregular. This means that the plane is divided into two domains by a curve which coincides with a straight line, except for a finite part, producing the irregularity. We propose algorithms for the approximate solution of the conjugation problem based on the spline methods for the solution of integral equations. We theoretically substantiate the computational scheme, namely, we prove the convergence and estimate the convergence rate.     

Scattering by Rough Surfaces: the Dirichlet Problem for the Helmholtz Equation in a Non-locally Perturbed Half-plane

We consider the Dirichlet boundary value problem for the Helmholtz equation in a non-locally perturbed half-plane, this problem arising in electromagnetic scattering by one-dimensional rough, perfectly conducting surfaces. We propose a new boundary integral equation formulation for this problem, utilizing the Green's function for an impedance half-plane in place of the standard fundamental solution. We show, at least for surfaces not differing too much from the flat boundary, that the integral equation is uniquely solvable in the space of bounded and continuous functions, and hence that, for a variety of incident fields including an incident plane wave, the boundary value problem for the scattered field has a unique solution satisfying the limiting absorption principle. Finally, a result of continuous dependence of the solution on the boundary shape is obtained.     

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.     

Furtivity and Masking Problems in Time-Dependent Electromagnetic Obstacle Scattering

In this paper, we consider furtivity and masking problems in time-dependent three-dimensional electromagnetic obstacle scattering. That is, we propose a criterion based on a merit function to minimize or to mask the electromagnetic field scattered by a bounded obstacle when hit by an incoming electromagnetic field and, with respect to this criterion, we drive the optimal strategy. These problems are natural generalizations to the context of electromagnetic scattering of the furtivity problem in time-dependent acoustic obstacle scattering presented in Ref. 1. We propose mathematical models of the furtivity and masking time-dependent three-dimensional electromagnetic scattering problems that consist in optimal control problems for systems of partial differential equations derived from the Maxwell equations. These control problems are approached using the Pontryagin maximum principle. We formulate the first-order optimality conditions for the control problems considered as exterior problems defined outside the obstacle for systems of partial differential equations. Moreover, the first-order optimality conditions derived are solved numerically with a highly parallelizable numerical method based on a perturbative series of the type considered in Refs. 2–3. Finally, we assess and validate the mathematical models and the numerical method proposed analyzing the numerical results obtained with a parallel implementation of the numerical method in several experiments on test problems. Impressive speedup factors are obtained executing the algorithms on a parallel machine when the number of processors used in the computation ranges between 1 and 100. Some virtual reality applications and some animations relative to the numerical experiments can be found in the website http://www.econ.unian.it/recchioni/w10/.     

Aspects of electromagnetic scattering in chiral media

In this work, we study two operators that arise in electromagnetic scattering in chiral media. We first consider electromagnetic scattering by a chiral dielectric with a perfectly conducting core. We define a chiral Calderon‐type surface operator in order to solve the direct electromagnetic scattering problem. For this operator, we state coercivity and prove compactness properties. In order to prove existence and uniqueness of the problem, we define some other operators that are also related to the chiral Calderon‐type operator, and we state some of their properties that they and their linear combinations satisfy. Then we sketch how to use these operators in order to prove the existence of the solution of the direct scattering problem. Furthermore, we focus on the electromagnetic scattering problem by a perfect conductor in a chiral environment. For this problem, we study the chiral far‐field operator that is defined on a unit sphere and contains the far‐field data, and we state and prove some of its properties that are preliminaries properties for solving the inverse scattering problem. Copyright © 2016 John Wiley & Sons, Ltd.     

A volume integral equation method for periodic scattering problems for anisotropic Maxwell's equations

This paper presents a volume integral equation method for an electromagnetic scattering problem for three-dimensional Maxwell's equations in the presence of a biperiodic, anisotropic, and possibly discontinuous dielectric scatterer. Such scattering problem can be reformulated as a strongly singular volume integral equation (i.e., integral operators that fail to be weakly singular). In this paper, we firstly prove that the strongly singular volume integral equation satisfies a Gårding-type estimate in standard Sobolev spaces. Secondly, we rigorously analyze a spectral Galerkin method for solving the scattering problem. This method relies on the periodization technique of Gennadi Vainikko that allows us to efficiently evaluate the periodized integral operators on trigonometric polynomials using the fast Fourier transform (FFT). The main advantage of the method is its simple implementation that avoids for instance the need to compute quasiperiodic Green's functions. We prove that the numerical solution of the spectral Galerkin method applied to the periodized integral equation converges quasioptimally to the solution of the scattering problem. Some numerical examples are provided for examining the performance of the method.     

On the numerical method for solving a hypersingular integral equation with the computation of the solution gradient

We consider a linear integral equation, which arises when solving the Neumann boundary value problem for the Laplace equation with the representation of the solution in the form of a double layer potential, with a hypersingular integral treated in the sense of Hadamard finite value. We consider the case in which the exterior or interior problem is solved in a domain whose boundary is a closed smooth surface and the integral equation is written over that surface. A numerical scheme for solving the integral equation is constructed with the use of quadrature formulas of the type of the method of discrete singularities with a regularization for the use of an irregular grid. We prove the convergence, uniform over the grid points, of the numerical solutions to the exact solution of the hypersingular equation and, in addition, the uniform convergence of the values of the approximate finite-difference derivative operator on the numerical solution to the values on the projection of the exact solution onto the subspace of grid functions with nodes at the collocation points.