Locating structural changes in a multiple scattering domain with an irregular shape

Locadiff is a method for imaging local structural changes in a random, heterogeneous medium. It relies on the combination of a forward model to calculate the sensitivity kernel of the source-receiver pairs, with an inversion method to determine the position of the changes. So far, the sensitivity kernel has been evaluated based on an analytical solution of the diffusion equation, which lacks the flexibility to handle problems where the domain has boundaries with an irregular shape. Moreover, the accuracy of the previous inversion method, based on linear algebra tools, was very sensitive to the values of the inversion parameters. This paper introduces a more generic approach to solve both these issues. The first problem is tackled by the implementation of a numerical method as an alternative for solving the diffusion equation. The second problem is tackled by the introduction of enhanced optimization algorithms to improve the stability of the inversion. This improved version of Locadiff is validated via both numerical examples and experimental data from an actual civil engineering problem.