Abstract and Applied Analysis

Volume 2012 (2012), Article ID 810676, 13 pages

doi:10.1155/2012/810676

Research Article

Inverse Scattering from a Sound-Hard Crack via Two-Step Method

Kuo-Ming Lee

Department of Mathematics, National Cheng Kung University, Tainan City 70101, Taiwan

Received 14 February 2012; Accepted 18 April 2012

Academic Editor: Ravshan Ashurov

Copyright © 2012 Kuo-Ming Lee. 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.

We present a two-step method for recovering an unknown sound-hard crack in from the measured far-field pattern. This method, based on a two-by-two system of nonlinear integral equations, splits the reconstruction into two consecutive steps which consists of a forward and an inverse problems. In this spirit, only the latter needs to be regularized.

1. Introduction

In this paper we consider an acoustic scattering problem from a sound-hard crack. This problem is modeled by an exterior boundary value problem governed by Helmholtz equation for an open arc with Neumann boundary conditions on both sides of the arc. Our major concern is the inverse problem which aim is to reconstruct the crack from some measurements. This kind of problem is of fundamental importance, for example, in material investigation, nondestructive testing, or in seismic exploration.

The inverse scattering problem for an open arc was first investigated by Kress [1]. In his article, integral equation method was used to solve both the direct and inverse problems for a sound-soft crack. The scattering problem in the unbounded domain was then converted into a boundary integral equation. Mönch [2] extended this approach to a Neumann problem. Their integral-equation-based method was to solve the so-called far-field equation.

for the unknown crack from the measured far-field data. Because of the nonlinearity and the compactness of the far-field operator, both linearization and regularization are needed to keep the solving of this equation in the framework of linear regularization theory. One of the major drawbacks of the regularized Newton’s method used in [1, 2] is that at every iteration, a direct problem with different boundary data must be solved for the Fréchet derivative of this far-field operator which is an essential part of this method.

The nonlinear integral equations method, proposed by Kress and Rundell [3], can be seen as a remedy for the Newton’s method. They proposed a two-by-two system of nonlinear integral equations which was shown to be equivalent to the original inverse problem. The Fréchet derivatives were obtained simply through the solving of this system. This method was extended to other boundary conditions and cracks in [4, 5] and to obstacle scattering in [6].

Another group of method which is not iterative is the decomposition method (see [7–9]). The idea is to split the nonlinearity and the ill-posedness of the original inverse problem. The scattered field is computed in the first linear ill-posed step from the far-field pattern . In the second step, which is nonlinear, the unknown crack is reconstructed as the location where the boundary condition fits in the least square sense. This method is thus carried out without computation of the direct problem. One of its disadvantages is that the reconstruction is less accurate than that of the Newton’s method.

In [10], a method consisting of two steps was proposed for a sound-soft crack. It is a mixture of the above mentioned methods. In this paper, we’d like to extend this method to the case of a sound-hard crack and to the case of limited aperture. The plan of the paper is as follows. For the sake of completeness and also for the introduction of notations, we will briefly summarize the main results of the direct problem in Section 2. In Section 3 we will consider the inverse scattering problem in its equivalence form of a two-by-two system of integral equations. In Section 4 we will discuss an iteration scheme for numerical computation of the inverse problem. This will be followed by some numerical examples in the final section.

2. Direct Neumann Problem

Given a regular nonintersecting -smooth open arc which can be represented as

The two end points of the crack are denoted by , respectively. The left hand side and the right-hand side of the crack are written by and , respectively. The unit outward normal to is denoted by . Further we set .

The direct scattering problem for a sound-hard crack that we are considering is as follows.

Problem 1 (DP). Given an incident plane wave with a wave number and a unit vector giving the direction of propagation, we find a solution to the Helmholtz equation

which satisfies the Neumann boundary conditions:

on both sides of the crack and the Sommerfeld radiation condition:

uniformly for all directions .

In (2.3), the limits

are required to exist in the sense of locally uniform convergence.

Note that the boundary conditions (2.3) can be reformulated as homogeneous Neumann conditions for the total field , that is, .

Using boundary integral equations, this direct problem can be solved via the layer approach. We refer to the monograph [11] for the details. In terms of the fundamental solution to the Helmholtz equation in

the following theorem, which was proved in [12], ensures the unique solvability of the direct scattering problem.

Theorem 2.1. The direct Neumann Problem 1 has a unique solution given by

where is the (unique) solution to the following integral equation

The function space is defined by

for .

At this place, we introduce the far-field pattern of the scattered field . The far-field pattern describes the behavior of the scattered wave at the infinity

uniformly for all directions . The one-to-one correspondence between radiating wave and its far-field pattern is established by the Rellich’s lemma. In the case of a sound-hard crack, the far-field pattern is given by

with the density function from Theorem 2.1 and the constant . We note that from the viewpoint of the inverse scattering problems, the far-field pattern is of particular interest since usually one has no access to the close neighborhood of the crack.

It is convenient for our further treatment to rewrite our integral equation (2.8) in operator form. For this purpose, we introduce the following integral operators:

The operator , which is termed the far-field operator in the literatures, plays an important role in the scattering theory. In terms of these operators, the solving of the direct scattering problem amounts to the solving of the following system of two operator equations:

In words, given the crack, the unknown density function is firstly solved from the first equation. Inserting into the very definition of the far-field operator, one obtains the far-field pattern from the second equation of the system (2.13) straight forward. This reasoning forms the basic idea of our numerical scheme for the inverse problem. At this place, we want to reduce the hypersingularity of the operator . Using Maue’s identity, the hypersingular operator is splitted into two parts (Theorem 7.29 in [11]):

where is the unit tangent vector.

3. Inverse Scattering Problem

After introducing the notations in the last section, we consider the following inverse problem.

Problem 2 (IP). Determine the crack if the far-field pattern is known for one incident wave.

For the uniqueness of this inverse problem, that is, for the identifiability of the arc, we refer to [2] where infinitely many incident waves are used. Although it is not yet proven, it is widely believed that the scatterer can be uniquely determined from one incident wave. Since our numerical method is based on Newton’s method which has the advantage of dealing with only one incident wave well, this leads to our problem setting. For the general solution theory of the inverse scattering problem, we refer to the monograph [13].

Motivated by the method of nonlinear integral equations in [3], a system of nonlinear integral equations which solution is the solution to the inverse problem was presented in [5] for a Neumann crack. In [5], the following system was proposed:

where the operators and are no more than the operators and in (2.13), respectively. The additional parameter refers to the dependence on the crack. In [5, Theorem 2], the equivalence of the inverse problem and the solving of the system (3.1) was proven. The basic idea of this system is that both the direct and inverse problem share the same integral equations, with slightly different interpretations. This makes the inverse problem really inverse to the original physical problem.

In [5], (3.1) was treated as a coupled system with two nonlinear ill-posed equations. This means that the two equations in (3.1) have to be linearized and regularized at the same time, and at every iterative step. The Fréchet derivatives for every operator both w.r.t. the unknown boundary and the unknown density have to be calculated at every step. Besides, two regularization parameters have to be selected for the scheme.

Start the inverse problem with the same system (3.1), we suggest a different scheme in the next section in details. Before formulating our algorithm, we need to parametrize the system.

As in Section 2, we use the parametrization for the sound-hard crack. To incorporate the square root singularities of the solution at the crack tips, we use the cosine substitution , as suggested in [14]. Keeping the notations in [5], we write for the crack and for the unknown density.

The system of integral equation (3.1), using Maue’s identity, is now parametrized as follows:

where

with the kernels

and the parametrized outward normal .

4. A Two-Step Method

Now we want to introduce our two-step algorithm. The parametrized systems (3.2), (3.3), or the original system (3.1) can be interpreted differently as in [5]. The first equation of the system can be seen as the equation for the direct problem. To be more precise, given the crack , which is equivalent to fixing the variable , (3.2) is solved for . Hence this equation is well posed according to Theorem 2.1. Moreover, it is a linear integral equation in . Another advantage is numerically as it can be just solved by the direct solver. With this at hand, we can solve (3.3) for . This is a nonlinear ill-posed equation which we will linearize and regularize. Based on Newton’s method, the Fréchet derivative of w.r.t. is needed for the linearization. The Fréchet derivative of the integral operator is simply given by the Fréchet derivative of its kernel (cf. [15]). For brevity, we set if . We have

where

The two-step scheme reads

Because of the compactness of the operator , (4.4) still has to be regularized. To accomplish this, we apply the Tikhonov regularization, that is, instead of (4.4), the equation

has to be solved with a regularization parameter .

At this point, we want to discuss the uniqueness and the solvability of the system (4.3), (4.5). From the solution theory of the direct problem in Section 2, (4.3) is uniquely solvable. However, the Fréchet derivative is not injective as pointed out in [1]. The null space of this operator is given by

This reflects the fact that different parametrizations of the arc leading to the same set of points turn out to have the same far-field pattern. We can avoid this ambiguity by limiting our solution space to the set of arcs representable as the graph of a function as suggested in [1]. Once this restriction is made, the operator is then an injective linear compact operator. As a consequence of the regularization theory, (4.5) is uniquely solvable, see [11]. Again, (4.3) is as a direct problem uniquely solvable as stated in Theorem 2.1. Thus, we can have the following theorem.

Theorem 4.1. If the pairs and are solutions to the system (4.3), (4.5), then and .

Numerically, we solve the systems (4.3), (4.5) in two steps: At the start, given an initial guess for the unknown crack . The first equation is treated as the corresponding direct problem and is solved for the density . Then the second equation updates the crack by means of the regularized Newton’s method.

The algorithm for our method can be formulated as follows:

(i)given an initial guess for the unknown crack;(ii)iterative steps, for .

(1)Step : solve

for .(2)Step : solve

for the update of .

(iii)

At this place we comment that from the numerical point of view our scheme is very attractive. As a variant of the nonlinear integral equations method, the Fréchet derivative of the integral operator is computed directly by solving the integral equation (4.8). Besides, the derivative is very simple as compared to those in [5] since only the far-field operator which has a smooth kernel is differentiated. Numerically it can be easily solved via the rectangular rule, for example. Another advantage of this method is that it splits the problem into two smaller parts and thus makes the computation cheaper. However, being a Newton’s method, the convergence of our numerical scheme is still open.

5. Numerical Results

In this section we will demonstrate the applicability of our method via some examples. We reconstruct the unknown crack from the knowledge of the far-field pattern at a number of points resulted from just one incident wave. For the direct problem, the forward solver is applied with 63 collocation points. In order to avoid committing an inverse crime, the number of collocation points used in the inverse problem is chosen to be different from that of the forward solver. The point is, in the inverse problem, (4.3) is solved with the same forward solver as for the direct problem. We therefore choose 31 collocation points in the inverse algorithm. In all our examples, we take the incident wave coming from the direction . The far-field pattern is measured at 15 points evenly distributed on the unit circle.

The basis functions for the parametrization are taken from the space:

where are the kth Chebyshev monomials . The selection of the Chebyshev polynomials is based on the fact that they can take care of the square root singularity in their own right. Thus the updates for the crack can be represented in the form:

with the unknown coefficients , , which have to be determined from (4.5). We test our scheme for two different wave numbers, . In the case where , is taken to be . For , we choose . The starting curve, that is, the initial guess for the regularized Newton’s method, is taken to be the straight line in all examples. According to the discrepancy principle, the stopping criterion for the iterative scheme is given by the relative error:

which is taken to be in case of exact data and in the case where data error are present. In all our figures, the dotted line (blue) represents the initial guess. We denote by the dashed line (red) the true solution and by the solid line (black) the reconstruction.

Example 5.1. For the first example, we take the arc

which is representable of a graph of a function. The numerical results are given in the Figures 1, 2 for and , respectively.

From the numerical results, we see that the reconstructions for exact data are very good. The reconstructions in the case where random are present, the reconstructions are not bad.

(a) Exact data,

(b) 3% error,

Figure 1: .

(a) Exact data,

(b) 3% error,

Figure 2: .

Example 5.2. To demonstrate the benefit of our method, we choose a curve which does not belong to our solution space. To this end, we take the arc

which is not representable as a function.

The results are shown in Figures 3 and 4.

(a) Exact data,

(b) 3% error,

Figure 3: .

(a) Exact data,

(b) 3% error,

Figure 4: .

We see that the reconstructions are very good, even in the case of erroneous data.

Example 5.3. In this final example, we choose the same curve as in the last example for the case of limited aperture. We assume that the data are only measurable within a certain range apart from the incident angle. Figure 5 shows the results for uniformly distributed far-field measurements within 90 degrees from both sides of the incident wave. We see that in this case, the reconstructions are very accurate both for the exact and erroneous data. In Figure 6, the results are shown for measurements taken within 45 degrees from both sides of the incident wave. We see that the reconstructions are convincible for only 15 data points.

(a) Exact data,

(b) 3% error,

Figure 5: .

(a) Exact data,

(b) 3% error,

Figure 6: .

We conclude this paper with some remarks. First of all, our examples above show the feasibility of the proposed numerical method. On one hand, being a Newton-type method, our method is conceptually simple and numerically more accurate than the traditional decomposition method. On the other hand, being a variant of the nonlinear integral equations method, the derivative is directly computable from the algorithm itself which makes the method more easily accessible than the classical Newton’s method. The splitting of the problem into two smaller parts makes our method more competitive as the computational cost concerns. Finally we want to point out that our method can be carried over to other type of boundary conditions and also to other type of scattering problems.

Acknowledgment

This work is partially supported by the NSC Grant NSC-100-2115-M-006-003-MY2.

References

R. Kress, “Inverse scattering from an open arc,” Mathematical Methods in the Applied Sciences, vol. 18, no. 4, pp. 267–293, 1995.

L. Mönch, “On the inverse acoustic scattering problem by an open arc: the sound-hard case,” Inverse Problems, vol. 13, no. 5, pp. 1379–1392, 1997.

R. Kress and W. Rundell, “Nonlinear integral equations and the iterative solution for an inverse boundary value problem,” Inverse Problems, vol. 21, no. 4, pp. 1207–1223, 2005.

O. Ivanyshyn and R. Kress, “Nonlinear integral equations for solving inverse boundary value problems for inclusions and cracks,” Journal of Integral Equations and Applications, vol. 18, no. 1, pp. 13–38, 2006.

K.-M. Lee, “Inverse scattering via nonlinear integral equations for a Neumann crack,” Inverse Problems, vol. 22, no. 6, pp. 1989–2000, 2006.

O. Ivanyshyn and R. Kress, “Nonlinear integral equations in inverse obstacle scattering,” in Proceedings of the 7th International Workshop on Mathematical Methods in Scattering Theory and Biomedical Engineering, Nymphaio, Greece, 2005.

A. Kirsch and R. Kress, “An optimization method in inverse acoustic scattering,” in Boundary Elements IX, Vol. 3: Fluid Flow and Potential Applications, C. A. Brebbia, W. L. Wendland, and G. Kuhn, Eds., pp. 3–18, Springer, Berlin, Germany, 1987.

R. Kress and P. Serranho, “A hybrid method for two-dimensional crack reconstruction,” Inverse Problems, vol. 21, no. 2, pp. 773–784, 2005.

R. Kress and P. Serranho, “A hybrid method for sound-hard obstacle reconstruction,” Journal of Computational and Applied Mathematics, vol. 204, no. 2, pp. 418–427, 2007.

K.-M. Lee, “A two step method in inverse scattering problem for a crack,” Journal of Mathematical Physics, vol. 51, no. 2, Article ID 023529, 10 pages, 2010.

R. Kress, Linear Integral Equations, Springer, Berlin, Germany, 2nd edition, 1999.

L. Mönch, “On the numerical solution of the direct scattering problem for an open sound-hard arc,” Journal of Computational and Applied Mathematics, vol. 71, no. 2, pp. 343–356, 1996.

D. Colton and R. Kress, Inverse Acoustic and Electromagnetic Scattering Theory, Springer, Berlin, Germany, 2nd edition, 1998.

Y. Yan and I. H. Sloan, “On integral equations of the first kind with logarithmic kernels,” Journal of Integral Equations and Applications, vol. 1, no. 4, pp. 549–579, 1988.

R. Potthast, “Fréchet differentiability of boundary integral operators in inverse acoustic scattering,” Inverse Problems, vol. 10, no. 2, pp. 431–447, 1994.