Advances in Mathematical Physics

Volume 2018, Article ID 2096058, 8 pages

https://doi.org/10.1155/2018/2096058

## Topological Derivative for Imaging of Thin Electromagnetic Inhomogeneity: Least Condition of Incident Directions

Correspondence should be addressed to Won-Kwang Park; rk.ca.nimkook@kwkrap

Received 19 June 2017; Accepted 28 January 2018; Published 22 February 2018

Academic Editor: Antonio Scarfone

Copyright © 2018 Won-Kwang Park. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

It is well-known that using topological derivative is an effective noniterative technique for imaging of crack-like electromagnetic inhomogeneity with small thickness when small number of incident directions are applied. However, there is no theoretical investigation about the configuration of the range of incident directions. In this paper, we carefully explore the mathematical structure of topological derivative imaging functional by establishing a relationship with an infinite series of Bessel functions of integer order of the first kind. Based on this, we identify the condition of the range of incident directions and it is highly depending on the shape of unknown defect. Results of numerical simulations with noisy data support our identification.

#### 1. Introduction

In this paper, we consider topological derivative [1] based imaging technique for thin, curve-like penetrable electromagnetic inhomogeneity with small thickness. Originally, this has been considered for shape optimization problems [2–7] and was then successfully combined with the level-set method (see [8–13]) for various inverse scattering problems. Surprisingly, throughout some researches [14, 15], it has been confirmed that topological derivative is also one of noniterative imaging techniques and very effective algorithm. However, as we can see [15–21], the most work has considered full-aperture problems.

Throughout results of numerical simulations, it has been confirmed that topological derivatives can be applied in limited-aperture problems, and an analysis in limited-aperture problems has been performed in [22]. In this interesting research, a relationship between topological derivative imaging function and an infinite series of Bessel functions of first kind has been established and, correspondingly, a sufficient condition of the range of incident directions for application has been identified theoretically. However, a least condition of application still remains unknown. Motivated by this fact, we identify a least condition of the range of incident directions for a successful application in limited-aperture inverse scattering problem and confirm this condition is highly depending on the unknown shape of thin inhomogeneity.

The remainder of paper is organized as follows. In Section 2, we survey two-dimensional direct scattering problem and topological derivative based imaging technique. In Section 3, we investigate a least condition of the range of incident directions and discuss its properties. In Section 4, several results of numerical simulations with noisy data are presented in order to support our investigation. A brief conclusion is given in Section 5.

#### 2. Introduction to Direct Scattering Problem and Topological Derivative

Let be a homogeneous domain with smooth boundary that contains a homogeneous thin inhomogeneity with a small thickness . That is,where is the unit normal to at and denotes the a simple, smooth curve in which describes the supporting curve of . In this contribution, we assume that the applied angular frequency is of the form . We assume that all materials are characterized by their dielectric permittivity and magnetic permeability at frequency of operation ; we define the piecewise constant permittivity and permeability asrespectively. For this sake, we set and denote as the wavenumber, where is a given wavelength and satisfies .

Let be the time-harmonic total field that satisfies Helmholtz equationwith boundary conditionand with transmission conditions on the boundary of . Here, denotes a two-dimensional vector on the connected, proper subset of unit circle such thatSimilarly, let be the background solution of (3) with boundary condition (4).

Now, we introduce the basic concept of topological derivative operated at a fixed single frequency. The problem considered herein is the minimization of the tracking type functional depending on the solution :

Assume that an electromagnetic inclusion of small diameter is created at a certain position , and let denote this domain. Since the topology of the entire domain has changed, we can consider the corresponding topological derivative based on with respect to point aswhere as . From (7), we can obtain an asymptotic expansion:

In [21], the following normalized topological derivative imaging function has been introduced:Here, and satisfying (8) for purely dielectric permittivity contrast ( and ) and magnetic permeability contrast ( and ) cases, respectively, are explicitly expressed as (see [21])where satisfies the adjoint problem

#### 3. Least Condition of Incident Directions

In this section, we identify the least condition of incident directions for applying topological derivative. For this, we introduce the structure of and as follows.

Lemma 1 (see [20, 21]). *Suppose that and are sufficiently large; thenwhere and are unit vectors that are, respectively, tangent and normal to the supporting curve at , and*

Based on Lemma 1, the structure of (9) in limited-aperture problem can be represented as follows. This result plays an important role in identifying least condition of incident directions. For a detailed proof, we refer to [22].

Theorem 2. *Let and . If and are sufficiently large, the structure of (9) becomeswherewithand the term does not contribute to the imaging performance.*

In recent work [20], it has been confirmed that the application of multifrequencies guarantees better a imaging performance than the application of a single frequency. Therefore, we consider the following normalized multifrequency topological derivative: for several frequencies , define

Based on the structure of in Theorem 2, we can observe that the terms and contribute to and disturb the imaging performance, respectively, for . Hence, eliminating the term will guarantee a good result. This means that the least condition comes fromfor all and . Note that since satisfies the asymptotic propertyfor , good results can appear in the map of when . Unfortunately, this is an ideal condition. Furthermore, since is arbitrary, we cannot control the value of the term . So, we must find a condition of and :for all . A simple way is to select and such thatthat is, for any , a selection and will guarantee good results via topological derivative. Based on this, we can obtain the following theoretical result of the range of incident directions, which has been examined heuristically.

Theorem 3 (least condition of range). *The least range of incident directions for successful application of topological derivative is highly depending on the shape of unknown thin inhomogeneity and the range of directions must be wider than .*

#### 4. Simulations Results and Discussions

In this section, some results of numerical simulations are exhibited to support identified condition mentioned in previous section. The homogeneous domain is chosen as a unit circle centered at the origin, and two supporting curves of the are selected as

The thickness of all is set to , and parameters , are chosen as 1. Let and for denote the permittivity and permeability of , respectively. The applied frequency is selected as , with , , and different incident directions chosen. In order to show the robustness, a white Gaussian noise with 20 dB signal-to-noise ratio (SNR) was added to the unperturbed boundary data.

First, let us consider the imaging of straight line shaped thin inhomogeneity . In this case, since , when , and will be a good choice for any . Corresponding results for , , , and are exhibited in Figure 1. Based on these results, we can observe that, for any value of , the shape of was retrieved satisfactorily via with the range of incident directions. However, although one can recognize the existence of inhomogeneity, the shape of cannot be reconstructed satisfactorily when .