Contrast source inversion methods in elastodynamics

In this paper a nonlinear inversion method is presented for determining the mass density of an elastic inclusion from the knowledge of how the inclusion scatters known incident elastic waves. The algorithm employed is an extension of the multiplicative regularized contrast source inversion method (MR-CSI) to elasticity. This method involves alternate determination of the mass density contrast and the contrast sources (the product of the contrast and the fields) in each iterative step. The simple updating schemes of the method allow the introduction of an extra regularization term to the cost functional as a multiplicative constraint. This so-called MR-CSI method (MR-CSI) has been proven to be very effective for the acoustic and electromagnetic inverse scattering problems. Numerical examples demonstrate that the MR-CSI method shows excellent edge preserving properties by robustly handling noisy data very well, even for more complicated elastodynamic problems.