Mathematical model of crack diagnosis: inverse acoustic scattering problem and its high-precision numerical solution

Victor A. Kovtunenko1

1Institute for Mathematics and Scientific Computing, Karl-Franzens University of Graz, NAWI Graz, Heinrich str. 36, 8010 Graz, Austria

1Lavrent’ev Institute of Hydrodynamics, Siberian Division of the Russian Academy of Sciences, 630090, Novosibirsk, Russia

1Corresponding author

Vibroengineering PROCEDIA, Vol. 22, 2019, p. 31-35.
Received 10 January 2019; accepted 20 January 2019; published 15 March 2019

Copyright © 2019 Victor A. Kovtunenko. 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.
Creative Commons License

The inverse acoustic scattering model for crack diagnosis is described by Helmholtz problem within mathematic framework and investigated for the sake of scientific computing. Minimizing the misfit from given measurements leads to an optimality condition-based imaging function which is used for non-iterative identification of the center of an unknown crack put in a test domain. The numerical tests are presented for the cracks of T-junction shape and are carried out based on the Petrov-Galerkin generalized FEM using wavelets basis and level-sets. This shows high-precision identification result and stability to noisy data of the diagnosis, which is illustrated for sound-soft as well as moderately sound-hard cracks when varying the coefficient of surface impedance.

  • The inverse acoustic scattering problem is investigated for scientific computing of crack diagnosis
  • Minimizing the misfit from given measurements leads to an optimality condition-based imaging
  • Imaging function is used for non-iterative identification of the center of a crack put in test domain
  • Numerical tests are based on Petrov-Galerkin GFEM and wavelet basis
  • T-junction shaped cracks are sound-soft as well as sound-hard with moderate surface impedance
  • High-precision identification result and stability to noisy data of the diagnosis are reported

Keywords: crack, acoustic scattering, inverse Helmholtz problem, optimization, imaging, variational method, level set, Petrov-Galerkin generalized FEM, wavelet basis, noisy data.

1. Introduction

In this contribution, the Helmholtz problem describing scattering of homogeneous acoustic medium that contains the single defect described by a cluster of cracks like T-junctions and characterized by surface impedance (thus, inhomogeneous medium) is considered. The inverse acoustic scattering problem is aimed at crack diagnosis motivated by applications to non-destructive testing with acoustic, electromagnetic, and elastic waves in engineering sciences.

The classic approach to inverse scattering problems in the context of mathematical modeling was established based on the operator theory, which results in direct (non-iterative) solution methods such as factorization, sampling, enclosure, MUSIC-type algorithms, and alike. The classic analysis utilizes fundamental solutions and Green’s functions for respective boundary value problems. In comparison, iterative methods are mostly used for the reason of numerical computing of the solution.

Recently, the shape optimization approach to identification of small geometric objects was developed based on the concept of topological derivative (when relating inhomogeneous medium to the background homogeneous one). The corresponding asymptotic analysis realizes the methods of singular perturbations. For the optimization approach adapted to inverse scattering problems, see the author’s works [1] related to sound-soft and sound-hard obstacles, and [2] for diagnosis of inhomogeneities in the background medium. The first-order topological derivative specified for geometric objects with real-valued refractive index was found in [3], and the high-order topological expansions were derived in [4].

Following the shape optimization, the imaging function is determined by the incident background field and the known measurement. In [1] it was applied for identification of the center of impedance obstacles from boundary measurements, and in [2] for identification of the center of an inhomogeneity from the known far-field pattern. The latter work used an equivalent formulation of the scattering problem in the form of a weakly singular integral equation of Lippmann-Schwinger. The numerical tests were carried out in 2d and 3d and demonstrated high-precision of the identification result and its stability with respect to the a-posteriori errors due to discretization as well as noisy data.

In the current contribution we extend the optimality condition-based imaging technique and successfully apply numerical algorithm for diagnosis of cracks. For the variational theory of cracks, see [5].

2. Theory

The theoretical result is established in the following manner. A least square cost functional of the misfit from a given measurement is minimized. To define the set of feasible parameters we describe the unknown crack ω with the help of geometric variables of its shape, center, size, and surface impedance coefficient α [0, ]. The weak variational formulation of the minimization problem with respect to these variables provides us with the zero-order necessary optimality condition, which is derivative-free. Further applying asymptotic arguments as the size tends to zero, for the sound-soft crack as α the optimality condition with respect to the center x* was derived in [1, 2] as follows:

I l x * = 0 ,

and the imaging function Il is defined for spatial points x by the following imaginary part:

I l x = Im u l x   v l ¯ x .

In Eq. (2) the complex-valued function ulx yields the incident field determined in the background homogeneous medium with the wave number k that satisfies the Helmholtz equation:

2 u l x + k 2 u l x = 0 ,       x ,

where stands for the gradient vector and 2 for the Laplace operator. The other complex-valued function vlx solves the adjoint Helmholtz problem:

2 v l x + k 2 v l x = 0 ,       x ,         n v l x = u l * x - u l x ,       x Γ ,

for the measurement ul*x given at the observation boundary Γ with the normal vector n, and vl¯ in Eq. (2) is its complex conjugate.

3. Methods

Motivated by the special need of highly accurate and stable interpolation for the variational solutions of the acoustic scattering problem, in [6] we developed the concept of Petrov-Galerkin enrichment by a finite element method (FEM). To reduce the discretization error, we enriched the space of test functions based on necessary optimality conditions for interpolation properties. Using a Petrov-Galerkin approach, we suggested low order interpolation polynomials for the trial space, and we enriched the test space with high order shape functions. Henceforth, the resulting Petrov-Galerkin enrichment (PGE) has low computational costs, but it improves significantly the accuracy of interpolation which is of the seventh order. In [6] justification of PGE was provided by local wavelets with vanishing moments based on the Gegenbauer polynomial approximation. Practical formulas were derived for calculation of the system matrix for the reference Helmholtz equation given over uniform meshes in 2d and 3d. In the current work we apply the PGE algorithm for computation of the Helmholtz problems Eqs. (4) and (7).

From Eqs. (1) and (2) it follows that the true center x* belongs to zero-level set Zl of Il that is:

x * Z l = x   I l x = 0 } ,

hence forces the numerical algorithm finding approximate center xh* from multiple measurements u1*, u2*, … and zero-level sets Z1h, Z2h, … by the mean of intersection:

x h * = l = 1,2 , Z l h   ,

where h associates the mesh size appeared in numerical computation of problem Eq. (4).

4. Results

We present our numerical tests over the uniform quadrilateral grid of size h= 2-7 in 2d implying the degrees of freedom DOF = 16641 in the computational domain = (0, 1) × (0, 1). The crack of T-junction shape and of size 0.0194 posed at the center x*= (0.1467, 0.3750) is drawn in Fig. 1 and Fig. 2. The measurements u1* and u2* are synthesized by solving numerically the forward scattering problem for l= 1, 2:

2 u l * x + k 2 u l * x = 0 ,       x \ ω ,         n ( u l * - u l ) x = 0 ,       x Γ ,
n u l * x + α u l * x = 0 ,         x ω ,

where the incident field is prescribed by the plane wave:

u l x = e i k ( x 1 cos θ l + x 2 sin θ l ) ,         i 2 = - 1 ,

in the scattering directions θ1=π/8 and θ2=π/4, and the wave number k=π/2 fixed.

We distinguish the crack scattering by varying the surface impedance coefficient α  [0, ] in Eq. (7) thus approaching the sound-hard α= 0 and the sound-soft α= cases. The typical identification result of xh* due to Eq. (6) for large α is depicted in Fig. 1 as α= 1.

Fig. 1. Identification of the center of crack when impedance α= 1 (approaching sound-soft case)

 Identification of the center of crack when impedance α= 1 (approaching sound-soft case)

a) Wave direction π/8

 Identification of the center of crack when impedance α= 1 (approaching sound-soft case)

b) Wave direction π/4

 Identification of the center of crack when impedance α= 1 (approaching sound-soft case)

c) Zero level sets

In Fig. 2(a) and 2(b) respectively there are drawn the numerically computed imaging functions I1h and I2h. Its zero-level sets Z1h and Z2h are utilized numerically by the narrow band technique and form the straight lines crossing the crack. Both the zero-level sets are depicted together in the computational test domain in Fig. 2(c). Their intersection point xh* is very close to the true center x* as can be verified in Table 1.

For comparison, the cross-point xh* for small α is presented in Fig. 2 as α= 0.01. This picture shows the error of approximation of x* caused by curving the zero-level sets.

All the results of numerical tests are gathered together in Table 1 also quantifying the absolute error |x*-xh*| of the identification algorithm. In particular, in the last column we observe the relative error |x*-xh*|/h given with respect to the mesh-size h. This fact clearly demonstrates that for moderate α 1 the Euclidean distance between the centers  x* and xh* is significantly less than the mesh-size (in the range of 1-30 % of h), which definitively justifies the high-precision property of our algorithm.

Fig. 2. Approximation of the center of crack when impedance α= 0.01 (approaching sound-hard case)

 Approximation of the center of crack when impedance α= 0.01 (approaching sound-hard case)

a) Wave direction π/8

 Approximation of the center of crack when impedance α= 0.01 (approaching sound-hard case)

b) Wave direction π/4

 Approximation of the center of crack when impedance α= 0.01 (approaching sound-hard case)

c) Zero level sets

Table 1. Identification result of center x*= (0.1467, 0.3750)

Surface impedance α
Approximate center xh*
Absolute error
Error related to h
(0.1490, 0.3750)
(0.1489, 0.3750)
(0.1479, 0.3750)
(0.1459, 0.3750)
(0.1479, 0.3752)
(0.1592, 0.3906)
(0.2031, 0.4280)
(0.2422, 0.4438)

The error  |x*-xh*| versus α from Table 1 is also illustrated for convenience in the log-log plot in Fig. 1. This curve exhibits the bounded limit behavior as α 0 and α.

Fig. 3. Identification error when varying impedance (sound-hard α= 0 and sound-soft α=)

 Identification error when varying impedance (sound-hard α= 0 and sound-soft α=∞)

5. Conclusions

Although the inverse acoustic scattering theory of Section 2 is derived rigorously when getting to the limit of the size and the surface impedance coefficient, the computer simulation ensures that the imaging function in Eq. (2) is still acceptable for moderate values of the parameters. This conclusion holds true also for the diagnosis of cracks treated here.


The author is supported by the Austrian Science Fund (FWF) Project P26147-N26: ‘Object Identification Problems: Numerical Analysis’ (PION) and the Austrian Academy of Sciences (OeAW), the RFBR and JSPS research Project 19-51-50004.


  1. Kovtunenko V. A., Kunisch K. High precision identification of an object: optimality conditions based concept of imaging. SIAM Journal of Control and Optimization, Vol. 52, Issue 1, 2014, p. 773-796. [Publisher]
  2. Cakoni F., Kovtunenko V. A. Topological optimality condition for the identification of the center of an inhomogeneity. Inverse Problems, Vol. 34, Issue 3, 2018, p. 035009. [Publisher]
  3. Kovtunenko V. A. Two-parameter topological expansion of Helmholtz problems with inhomogeneity. Mathematical Analysis of Continuum Mechanics and Industrial Applications. Mathematics for Industry, Vol. 26, 2017, p. 51-81. [CrossRef]
  4. Kovtunenko V. A. High-order topological expansions for Helmholtz problems in 2d. Topological Optimization and Optimal Transport, Radon Series on Computational and Applied Mathematics, Vol. 17, 2017, p. 64-122. [CrossRef]
  5. Khludnev A. M., Kovtunenko V. A. Analysis of Cracks in Solids. WIT-Press, Southampton, Boston, 2000. [CrossRef]
  6. Kovtunenko V. A., Kunisch K. Revisiting generalized FEM: a Petrov-Galerkin enrichment based FEM interpolation for Helmholtz problem. Calcolo, Vol. 55, 2018, p. 38. [Publisher]