Ruming Zhang

Ruihao Huang, Allan A. Struthers, Jiguang Sun et al.

Recently, a new eigenvalue problem, called the transmission eigenvalue problem, has attracted many researchers. The problem arose in inverse scattering theory for inhomogeneous media and has important applications in a variety of inverse problems for target identification and nondestructive testing. The problem is numerically challenging because it is non-selfadjoint and nonlinear. In this paper, we propose a recursive integral method for computing transmission eigenvalues from a finite element discretization of the continuous problem. The method, which overcomes some difficulties of existing methods, is based on eigenprojectors of compact operators. It is self-correcting, can separate nearby eigenvalues, and does not require an initial approximation based on some a priori spectral information. These features make the method well suited for the transmission eigenvalue problem whose spectrum is complicated. Numerical examples show that the method is effective and robust.

Armin Lechleiter, Ruming Zhang

Scattering problems for periodic structures have been studied a lot in the past few years. A main idea for numerical solution methods is to reduce such problems to one periodicity cell. In contrast to periodic settings, scattering from locally perturbed periodic surfaces is way more challenging. In this paper, we introduce and analyze a new numerical method to simulate scattering from locally perturbed periodic structures based on the Bloch transform. As this transform is applied only in periodic domains, we firstly rewrite the scattering problem artificially in a periodic domain. With the help of the Bloch transform, we secondly transform this problem into a coupled family of quasiperiodic problems posed in the periodicity cell. A numerical scheme then approximates the family of quasiperiodic solutions (we rely on the finite element method) and backtransformation provides the solution to the original scattering problem. In this paper, we give convergence analysis and error bounds for a Galerkin discretization in the spatial and the quasiperiodicity's unit cells. We also provide a simple and efficient way for implementation that does not require numerical integration in the quasiperiodicity, together with numerical examples for scattering from locally perturbed periodic surfaces computed by this scheme.

Armin Lechleiter, Ruming Zhang

Periodic surface structures are nowadays standard building blocks of optical devices. If such structures are illuminated by aperiodic time-harmonic incident waves as, e.g., Gaussian beams, the resulting surface scattering problem must be formulated in an unbounded layer including the periodic surface structure. An obvious recipe to avoid the need to discretize this problem in an unbounded domain is to set up an equivalent system of quasiperiodic scattering problems in a single (bounded) periodicity cell via the Floquet-Bloch transform. The solution to the original surface scattering problem then equals the inverse Floquet-Bloch transform applied to the family of solutions to the quasiperiodic problems, which simply requires to integrate these solutions in the quasiperiodicity parameter. A numerical scheme derived from this representation hence completely avoids the need to tackle differential equations on unbounded domains. In this paper, we provide rigorous convergence analysis and error bounds for such a scheme when applied to a two-dimensional model problem, relying upon a quadrature-based approximation to the inverse Floquet-Bloch transform and finite element approximations to quasiperiodic scattering problems. Our analysis essentially relies upon regularity results for the family of solutions to the quasiperiodic scattering problems in suitable mixed Sobolev spaces. We illustrate our error bounds as well as efficiency of the numerical scheme via several numerical examples.

Armin Lechleiter, Ruming Zhang

This paper concerns the inverse scattering problem to reconstruct a local perturbation in a periodic structure. Unlike the periodic problems, the periodicity for the scattered field no longer holds, thus classical methods, which reduce quasi-periodic fields in one periodic cell, are no longer available. Based on the Floquet-Bloch transform, a numerical method has been developed to solve the direct problem, that leads to a possibility to design an algorithm for the inverse problem. The numerical method introduced in this paper contains two steps. The first step is initialization, that is to locate the support of the perturbation by a simple method. This step reduces the inverse problem in an infinite domain into one periodic cell. The second step is to apply Newton-CG method to solve the associated optimization problem. The perturbation is then approximated by a finite spline basis. Numerical examples are given at the end of this paper, shows the efficiency of the numerical method.

Ruming Zhang

In this paper, we will introduce a high order numerical method to solve the scattering problems with non-periodic incident fields and (locally perturbed) periodic surfaces. For the problems we are considering, the classical methods to treat quasi-periodic scattering problems no longer work, while a Bloch transform based numerical method was proposed in [LZ17b]. This numerical method, on one hand, is able to solve this kind of problems convergently; on the other hand, it takes up a lot of time and memory during the computation. The motivation of this paper is to improve this numerical method, from the regularity results of the Bloch transform of the total field, which have been studied in [Zha17]. As the set of the singularities of the total field is discrete in $\R$, and finite in one periodic cell, we are able to improve the numerical method by designing a proper integration contour with special conditions at the singularities. With a good choice of the transformation, we can prove that the new numerical method could possess a super algebraic convergence rate. \high{This new method improves the efficient significantly. At the end of this paper, several numerical results will be provided to show the fast convergence of the new method.} The method also provides a possibility to solve more complicated problems efficiently, e.g., three dimensional problems, or electromagnetic scattering problems.

Bo Zhang, Ruming Zhang

The integral equation method is widely used in numerical simulations of 2D/3D acoustic and electromagnetic scattering problems, which needs a large number of values of the Green's functions. A significant topic is the scattering problems in periodic domains, where the corresponding Green's functions are quasi-periodic. The quasi-periodic Green's functions are defined by series that converge too slowly to be used for calculations. Many mathematicians have developed several efficient numerical methods to calculate quasi-periodic Green's functions. In this paper, we will propose a new FFT-based fast algorithm to compute the 2D/3D quasi-periodic Green's functions for both the Helmholtz equations and Maxwell's equations. The convergence results and error estimates are also investigated in this paper. Further, the numerical examples are given to show that, when a large number of values are needed, the new algorithm is very competitive.

Ruming Zhang

In this paper, we will study the Bloch transformed rough surface scattering problems, and propose a numerical method based on the Bloch transformed problems. Based on the mathematical theory of the scattering problems from locally perturbed periodic surfaces, the same techniques will be applied to the rough surface scattering problems, and an equivalent coupled family of quasi-periodic scattering problems in one periodic cell will be established. The most important result obtained in this paper is on the finite Fourier series approximation of the Bloch transformed field with respect to the quasi-periodicity parameter. It will be proved that the finite series is exactly the Bloch transformed solution corresponds to truncated rough surfaces. Thus the truncation provides a reasonable approximation, and could be applied to the numerical solutions. Based on the approximation, a numerical method is proposed for the rough surface scattering problems. The convergence of the numerical method is proved and illustrated by the numerical experiments. The method provides a completely new perspective for the rough surface scattering problems. There is possibility that some high order method will be developed based on this new method.

Ruming Zhang

Waves scattering from unbounded structures are always complicated problems for numerical simulations. For the case that the non-periodic incident field scattered by (locally perturbed) periodic surfaces, with the help of the Bloch transform, the problem could be solved by some finite element methods, if the incident fields decay at certain rate at the infinity. For faster decaying incident fields, a high order numerical method is also available. However, in these cases, the plain waves, which belong to a very important family of incident fields but do not decay at the infinity, are not included. In this paper, we aim to develop the Bloch transform based standard finite method for this certain case, and then establish the high order method afterwards. Numerical experiments have been carried out for both the standard and high order numerical methods. Based on the algorithms for incident plain waves, we could also extend the numerical methods to more generalized cases when only the not so efficient standard method is available.

Ruming Zhang

In this paper, we develop a high order numerical method for the numerical solutions of scattering problems with slightly perturbed periodic surfaces in two dimensional spaces. Based on the regularity property introduced in Part I, the decaying rate of the incident field could be transferred directly to the total field for small perturbations. Thus the finite section method could reach a high accuracy rate. With the help of a modification of the truncated problem, the problem is solved by a finite element method. The convergence of the finite element method is proved and numerical examples have been carried out to show the efficiency of the numerical scheme.