The focus of this work is on rigorous mathematical analysis of the topological derivative based detection algorithms for the localization of an elastic inclusion of vanishing characteristic size. A filtered quadratic misfit is considered, and the performance of the topological derivative imaging functional resulting therefrom is analyzed. Our analysis reveals that the imaging functional may not attain its maximum at the location of the inclusion. Moreover, the resolution of the image is below the diffraction limit. Both phenomena are due to the coupling of pressure and shear waves propagating with different wave speeds and polarization directions. A novel imaging functional based on the weighted Helmholtz decomposition of the topological derivative is, therefore, introduced. It is thereby substantiated that the maximum of the imaging functional is attained at the location of the inclusion and the resolution is enhanced and proves to be the diffraction limit. Finally, we investigate the stability of the proposed imaging functionals with respect to measurement and medium noises.

We consider reflector imaging in a weakly random waveguide. We address the situation in which the source is farther from the reflector to be imaged than the energy equipartition distance, but the receiver array is closer to the reflector to be imaged than the energy equipartition distance. As a consequence, the reflector is illuminated by a partially coherent field and the signals recorded by the receiver array are noisy. This paper shows that migration of the recorded signals cannot give a good image, but an appropriate migration of the cross correlations of the recorded signals can give a very good image. The resolution and stability analysis of this original functional shows that the reflector can be localized with an accuracy of the order of the wavelength even when the receiver array has small aperture, and that broadband sources are necessary to ensure statistical stability, whatever the aperture of the array.

This paper analyzes the random fluctuations obtained by a heterogeneous multi-scale first-order finite element method applied to solve elliptic equations with a random potential. Several multiscale numerical algorithms have been shown to correctly capture the homogenized limit of solutions of elliptic equations with coefficients modeled as stationary and ergodic random fields. Because theoretical results are available in the continuum setting for such equations, we consider here the case of a secondorder elliptic equations with random potential in two dimensions of space. We show that the random fluctuations of such solutions are correctly estimated by the heterogeneous multi-scale algorithm when appropriate fine-scale problems are solved on subsets that cover the whole computational domain. However, when the fine-scale problems are solved over patches that do not cover the entire domain, the random fluctuations may or may not be estimated accurately. In the case of random potentials with short-range interactions, the variance of the random fluctuations is amplified as the inverse of the fraction of the medium covered by the patches. In the case of random potentials with long-range interactions, however, such an amplification does not occur and random fluctuations are correctly captured independent of the (macroscopic) size of the patches. These results are consistent with those obtained in [9] for more general equations in the one-dimensional setting and provide indications on the loss in accuracy that results from using coarser, and hence computationally less intensive, algorithms.

The aim of this paper is to develop a mathematical framework for opto-elastography. In opto-elastography, a mechanical perturbation of the medium produces a decorrelation of optical speckle patterns due to the displacements of optical scatterers. To model this, we consider two optically random media, with the second medium obtained by shifting the first medium in some local region. We derive the radiative transfer equation for the cross-correlation of the wave fields in the media. Then we derive its diffusion approximation. In both the radiative transfer and the diffusion regimes, we relate the correlation of speckle patterns to the solutions of the radiative transfer and the diffusion equations. We present numerical simulations based on our model which are in agreement with recent experimental measurements.

The aim of this paper is to provide a fast and efficient procedure for (real-time) target identification in imaging based on matching on a dictionary of precomputed generalized polarization tensors (GPTs). The approach is based on some important properties of the GPTs and new invariants. A new shape representation is given and numerically tested in the presence of measurement noise. The stability and resolution of the proposed identification algorithm is numerically quantified. We compare the proposed GPT-based shape representation with a moment-based one.