Reconstruction of Local Perturbations in Periodic Surfaces
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.
In this paper, we will introduce the numerical method of the inverse scattering problem in a locally perturbed periodic structure. Both the direct and inverse scattering problems in periodic structures have been studied in the last few years, especially for the case that the incident fields are quasi-periodic, e.g. plain waves. A classical way is to reduce the problems defined in an infinite periodic domain into one periodic cell, then the finitely defined problems could be solved in normal methods. However, if the periodicity of the solutions is destroyed, i.e., the incident field is not (quasi-)periodic, or the periodic structure is perturbed, the classical methods are no longer available and new techniques are needed.
The Floque-Bloch transform has been applied to perturbed periodic structures in [Coa12]. The direct scattering problems with non-periodic incident fields (Herglotz wave functions) have been studied in [LN15] and [Lec17] theoretically, for numerical method see [LZ17a, LZ16]; problems with locally perturbed periodic surfaces have also been studied, for theoretical part see [Lec17] and for numerical method see [LZ17b]. The Bloch transform could also be applied to waveguide problems, see [HN16]. In this paper, the analysis and numerical solutions of the direct problems are based on these results.
The numerical method developed in this paper is a combination of an initialization and an iteration scheme. The initialization step is to locate the perturbation from a relatively larger area. As was introduced by Ito, Jin and Zou in [IJZ12a], a sampling method that only involves one incident field and a simple evaluation of an integration, could roughly locate the inhomogeneity embedded in free space. For more results for this method, see [IJZ12b, IJZ13, LZ13]. Followed by their idea, we will design an initialization algorithm to locate the perturbation, such that we could continue the iteration step in a finite domain. In the next step, we will approximate the perturbation by finite number of spline basis, and rewritten the inverse problem as an optimization problem. In this paper, we will apply Newton-CG method (see [EHN96]) to solve the associated optimization method.
The rest of the paper is organized as follows. In Section 2, a description of the direct problem is made, together with the Green’s function in periodic structures. The inverse problem is formulated in Section 3, and the Fréchet derivative and its adjoint operator are studied. In Section 4, we will introduce the numerical methods for the inverse problem. In Section 5, several numerical examples are shown to illustrate the efficiency of the numerical method.
2 Direct Scattering Problem
Suppose is a -periodic interface in defined by the function that is -periodic. is a local perturbation of , where the function has a compact support in . is a straight line above and , where . For some , define . Define the unbounded domains , , , as
Define , and the define periodic cell s
Let , , , denote the domains , , , restricted in one period .
For simplicity, we assume that the functions and are smooth enough, and the support of lies in one periodic cell, i.e., there is some such that .
Let the incident field be a downward propagating Herglotz wave function, i.e.,
where , is the closure of the space in the norm
Then the total field satisfies the equations
with the scattered field is the upward propogating field, satisfies
where is defined by
It is a bounded operator from to for all , see [CE10]. Thus the total field satisfies the following boundar condition on :
Recall the definition of the function space . The space is defined as the closure of with the norm
for any . Define the space and as
for all , where . From [CE10], the problem is uniquley soluable in for any and .
For and any incident field , the variational problem (5) is uniquely solvable in .
We can also define the total and scattered fields with the non-perturbed interface , and the unique solvability in Theorem 2 also holds for this situation. The total field is denoted by and the scattered field is .
2.2 Bloch Transform
Define the Bloch transform in the periodic domain for
Define the functions space by the closure of with the following norm when
The definition could be extended to all by interpolation and duality arguments. With this function space, the operator has the following properties, see [Lec17].
The Bloch transform extends to an isomorphism between and . Further, is an isometry for with inverse
and the inverse transform equals to the adjoint operator of .
Bloch transform only defined on periodic domains, so following [LZ17b], define the diffeomorphism mapping onto as
then . Define the transformed total field , then from (5) it satisfies the following variational problem
for all , where
It is easy to deduce that the support of and are all subsets of . From the solution , the transformed function .
Let and , then it satisfies
For any , for each fixed , the -quasi-periodic solution has the following representation, i.e., the Rayleigh expansion
From [LZ17b], the variational problem is uniquely solvable in certain conditions.
Suppose is the graph of a Liptshitz continuous function.
The variational problem (7) is uniquely solvable in .
If for and , then .
For the fields with non-perturbed surface , there is no need to do the transformation . We denote the Bloch transform of by , and Theorem 5 also holds for this problem.
2.3 Green’s Functions in the Periodic Domain
Suppose is the Green’s function located in , then it satisfies
For we note that is the incident half-space Green’s function, where , and . From Theorem 5, with the fact that , the Green’s function is well-defined in . Moreover, .
Theorem 7 (Theorem III.1, [Lec08]).
The Green’s function is symmetric, i.e., with .
Suppose and , then from the boundary condition (11), if has the following form of Fourier transform, i.e.,
then its normal derivative on
Suppose are two different points in . Then the Green’s function satisfies the following property
As , for any , define (and similar for ), then . Denote , then
With these properties,
From the asympototic behaviers of and ,
From the Green’s identity and Dirichlet boundary conditions, for small enough
where the right hand side tends to as , i.e.,
As the Green’s function has the integral representation,
Define the function
the following property is also easy to prove
When is large enough, from the definition of UPRC,
When and are close enough, , thus
As has a singularity at , the imaginary part of the Green’s function may possess a sigularity, or a local maximum at that point. From the approximate behavior of when when falls in a small enough neighbourhood of , it could be expected to obtain relatively large values.
3 The Inverse Problem
Suppose and is a -periodic knowns function, and define the space
In the following, we will always assume that and , and , thus from Theorem 5, . Then the functions , . Define the following scattering operator
Suppose data is measured on , which is the scattered field with some noise, denoted by , then the inverse problem is described as follows.
Inverse Problem: Given a measured data , to find such that
In the rest part of this section, we will describe some properties of the scattering operator.
The proof is similar to that in [Kir93] so we omit it here.
Now we will study the property of , before that, the following property of the DtN operator is needed.
is the duality pairing defined by the inner product of . The DtN operator satisfies .
Suppose , then from the definition of the DtN map,
Then we have the following calculations.
From similar procedure,
the proof is finished from the denseness of in the space . ∎
With these results, the adjoint operator of from to is given in the following theorem.
For any and , as is the space of real-valued functions,
Use the Helmholtz equation with the homogeneous Dirichlet boundary condition, i.e., on , and the condition (21),
then we have the final results