An inverse problem for reconstructing arbitrary-shaped thin penetrable electromagnetic inclusions concealed in a homogeneous material is considered in this paper. For this purpose, the level-set evolution method is adopted. The topological derivative concept is incorporated in order to evaluate the evolution speed of the level-set functions. The results of the corresponding numerical simulations with and without noise are presented in this paper.

1. Introduction

In general, one of the main purposes of the inverse scattering problem is to determine the characteristics of an unknown target, for example, the shape, location, and internal constitution, based on the measured scattered field data. Therefore, there exist a considerable number of significant inverse scattering problems for identifying the information of arbitrary-shaped thin penetrable electromagnetic inhomogeneities or perfectly conducting cracks concealed in structures such as concrete walls of buildings and bridges (refer to [1] and the references therein).

In order to solve related inverse problems, various algorithms have been developed. Among them, an iterative algorithm based on level-set evolution has been established in [2] and many generalizations have been investigated (refer to [3] and the references therein). However, most research studies have focused on the shape reconstruction of volumetric targets. Therefore, in order to reconstruct crack-like targets via level-set evolution, a new concept is required.

In recent work [4], arbitrary-shaped thin penetrable conductivity inclusions have been successfully reconstructed via the level-set method by adopting two different level-set functions. Shape reconstruction of thin electromagnetic inclusions is considered in [5]. However, in [5], only the case of the purely permittivity contrast or permeability contrast with respect to the background material is considered. So, extension to both contrast cases is still remaining as a research topic.

Motivated by the mentioned above, the main objective of this paper is to reconstruct the shape of thin electromagnetic inclusions by adopting the level-set method in the case of both the dielectric permittivity and the magnetic permeability contrasts with respect to the background material. Note that, in general research papers, the evolution speed is evaluated via a complex calculation of the Fréchet derivatives using an adjoint technique. Instead of such evaluations, we incorporate a topological derivative to derive the evolution speed of two level-set functions introduced in [610].

This paper is organized as follows. In Section 2, we introduce two-dimensional direct scattering from an arbitrary-shaped thin penetrable electromagnetic inclusion and the level-set functions for describing such inclusions. In Section 3, the evolution speed of the level-set functions is derived via the topological derivative. We present some numerical simulation results in Section 4 to show the effectiveness of the level-set method. Finally, we provide short conclusions in Section 5.

2. Direct Scattering Problem and Two Level-Set Functions: A Brief Survey

2.1. Two-Dimensional Direct Scattering Problem

Let be a two-dimensional homogeneous domain with a smooth boundary . Throughout this paper, we assume that conceals an electromagnetic inclusion with thickness , which is characterized in the neighborhood of the smooth curve : where is a unit normal to at . Throughout this paper, we assume that does not touch ; that is, we assume that .

In this paper, the electric permittivities and magnetic permeabilities of and at a given nonzero frequency are assumed to be known; they are finite valued and they differ (either one or both of them) from the ones of the homogeneous embedding medium. For the sake of simplicity, let and denote the permittivity and permeability of , and let and denote those of , respectively. Then, we can define the following piecewise constants:

Let denote a time-harmonic total field that satisfies the following Helmholtz equation: with boundary condition where , denotes the unit outward normal to , and denotes a vector on the unit circle . In the same manner, we denote as a background solution of (3) with boundary condition (4).

2.2. Description of Thin Inclusion Using Two Level-Set Functions

Now, we describe the thin inclusion by two different level-set functions (see [4, 5] for a detailed description). From now on, we assume that the level-set function is a function. With this assumption, we introduce the zero level-set function and its boundary as respectively. Then, a thin region of thickness and inner boundary can be written as follows: respectively.

In order to describe a thin region of finite length (which is connected or disconnected) in , we assume that another level-set function is also a continuously differentiable function. Let us define a band structure by Then, we can describe as and the boundary of for evolving at inclusion tips as

Figure 1 shows the illustration of the thin inclusion using two level-set functions.

3. Evolution Speed of the Level-Set Functions and the Topological Derivative

The evolution of the level-set function satisfies the Hamilton-Jacobi type equation: where is or introduced in Section 2.2 and is given by (see, e.g., [4, 5]) the following. (i)For level-set function and , (ii)For level-set function and , Here, denotes the speed function of the level sets that we have to evaluate. In a recent work [5], has been derived for the purely permittivity contrast ( and ) or permeability contrast ( and ) cases by evaluating the Fréchet derivative. In order to evaluate for both the permittivity and permeability contrast cases, we introduce the following discrepancy functional: Then, we adopt the following relationship between the evolution speed of level sets and the topological derivative.

Theorem 1. Let denote the topological derivative of (13); then satisfies

Proof. A rigorous derivation can be found in [6, 8, 9].

Hence, we can evaluate the evolution speed by computing the topological derivative . The result is as follows.

Theorem 2. In both the permittivity and permeability contrast cases, can be written as where satisfies the following adjoint problem:

Proof. See the appendix.

Note that, if either or , we can obtain the Fréchet derivatives based on [5, Theorem 3.1] or [5, Theorem 3.2], respectively. This result is identical to that of [11, Theorem 3.2], but the proof is considerably easier.

4. Numerical Simulations

In this section, we choose the background material as a unit circle centered at origin with . For thin inclusions , the thickness is set to and the parameters of are chosen as , where and are the permittivity and permeability contrasts of , respectively. The supporting curve is selected as

For a given inclusion , let be the solution to (3) with boundary condition (4). In the same manner, for an inclusion at the iteration step (if , it is an initial guess), let denote the solution to (3) with boundary condition (4). For the least-square functional at a given frequency, we introduce a normalized error function: where In this section, we adopt and the locations are selected as

A finite-difference time discretization form of (10) is For the evolving level sets and , the time step is chosen as (see [5, 12] for a detailed discussion). With this configuration, the forward and adjoint problems are solved via the finite element method (FEM) to avoid an inverse crime, and a 20 dB white Gaussian random noise is added to the boundary measurements data in order to present the robustness.

Remark 3 (generating a good initial guess). It is worth mentioning that, for a successful performance of the iterative based algorithm, one must start the iteration procedure with a good initial guess that is close to the unknown target. Without this guess, one will encounter unexpected scenarios such as nonconvergence (refer to Figures  7 and  14 in [5]). There are various noniterative shape reconstruction algorithms such as the multiple signal classification (MUSIC)-type algorithm [1317], the end-point identification algorithm [1820], the multifrequency based subspace migration imaging technique [2127], the topological derivative strategy [10, 2832], and the gradient-for-the-initiation (GFI) method [33].

First, we consider the shape reconstruction of . In this example, we assume that the location of the end points of has been identified based on the location search algorithm in [1820]. Thus, by connecting them by a straight line segment, a good initial guess can be obtained.

Figure 2 shows the shape reconstruction of at the operated frequency . Based on this result, we can confirm that the shape of is successfully reconstructed and the residual considerably decreases accordingly.

Next, we consider the reconstruction of . In contrast to the previous example, we apply the topological derivative strategy introduced in [10, 30] to obtain a good initial guess. In Figure 3, the obtained shape of via the topological derivative with and is illustrated. In this result, we can construct a supporting curve via well-identified points as for (see Figure  9 in [10]). Here, denotes the Chebyshev polynomials of the first kind. Considering this curve as the supporting curve of the initial guess, the reconstruction of the complete shape of is illustrated in Figure 3. Note that, because of the considerably close initial guess, only 31 iterations are sufficient to obtain a good result.

5. Conclusion Remarks

In this paper, a level-set method to perform the shape reconstruction of thin electromagnetic inclusion using both the dielectric permittivity and magnetic permeability contrasts with respect to the homogeneous domain was considered. The relationship between the topological derivative and the evolution speed of level-set functions allowed us to reconstruct the complete shape of the thin inclusion with a good initial guess. It is discovered from the numerical experiments that the proposed technique is stable with respect to random noise. Further, incorporating the topological derivative guarantees a successful reconstruction procedure.

In this paper, we performed shape reconstruction with a priori information of unknown targets, for example, thickness, permittivity, and permeability values. Performing simultaneous reconstruction without a priori information will be a challenge (refer to [3, Section 10.3]).


Proof of Theorem 2. In this appendix, we prove Theorem 2. In order to derive of (13), we construct a small electromagnetic inclusion in with the same permittivity and permeability of and denote this domain by . Then, by virtue in [10], the following asymptotic expansion holds: where as . Let be a solution to (3) with boundary condition (4) in the existence of . Then, the following asymptotic expansion formula holds (see [13]): where denotes the Neumann function that is the solution to
Then, applying (A.2) to (13), we can compute where Then, applying boundary condition (16) and asymptotic formula (A.2) yields
First, applying integration by parts and (A.3), we can compute as
Similarly, becomes Hence, and by considering the real part (A.9) and comparing (A.1) and (A.4), we can observe that and This completes the proof.


The author would like to acknowledge the anonymous reviewer for his valuable comments. This work was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF), funded by the Ministry of Education, Science, and Technology (no. 2011-0007705), and the Research Program of the Kookmin University in Korea.