Weight of single and recurrent scattering in the reflection matrix of complex media.

In a heterogeneous medium, the wave field can be decomposed as an infinite series known as the Born expansion. Each term of the Born expansion corresponds to a scattering order, it is thus theoretically possible to discriminate single and multiple scattering contribution to the field. Experimentally, what is actually measured is the total field in which all scattering orders interfere. Conventional imaging methods usually rely on the assumption that the multiple scattering contribution can be disregarded. In a back-scattering configuration, this assumption is valid for small depths, and begins to fail for depths larger than the scattering mean-free path ℓ_{s}. It is therefore a key issue to estimate the relative amount of single and multiple scattering in experimental data. To this end, a single-scattering estimator ρ[over ̂] computed from the reflection matrix has been introduced in order to assess the weight of single scattering in the backscattered wave field. In this paper, the meaning of this estimator is investigated and a particular attention is given to recurrent scattering. In a diffraction-limited experiment, a multiple scattering sequence is said to be recurrent if the first and last scattering events occur in the same resolution cell. Recurrent scattering is shown to be responsible for correlations between single scattering and higher scattering orders of the Born expansion, inducing a bias to the estimator ρ[over ̂] that should rather be termed confocal scattering ratio. Interestingly, a more robust estimator is built by projecting the reflection matrix in a focused basis. The argument is sustained by numerical simulations as well as ultrasonic data obtained around 1.5 MHz in a model medium made of nylon rods immersed in water. From a more general perspective, this work raises fundamental questions about the impact of recurrent scattering on wave imaging.

Evaluation of a multiple scattering filter to enhance defect detection in heterogeneous media.

Ultrasonic evaluation of coarse-grain materials generates multiple scattering at high frequency and large depth. Recent academic experiments with array probes showed the ability of a random matrix method [multiple scattering filter (MSF)] to reduce multiple scattering, hence improving detection. Here, MSF is applied to an industrial nickel-based alloy with coarse-grain structure. Two samples with average grain sizes 90 ± 60 µm and 750 ± 400 µm are inspected with wide-band 64-element arrays at central frequencies 2, 3, and 5 MHz. They contain cylindrical through-holes (1-mm radius) at various depths. The array transfer matrix is recorded and post-processed both in the flawless area and for eleven positions above each defect, which allows for a statistical analysis. MSF is compared with two conventional imaging techniques: the total focusing method (TFM) and the decomposition of the time-reversal operator (DORT). Several parameters to assess the performance of detection techniques are proposed and discussed. The results show the benefit of MSF, especially at high frequencies and for deep defects: at 5 MHz and 70 mm depth, i.e., more than three scattering mean-free paths, the detection rate for MSF ranges between 55% and 100% while it is found to be 0% both for TFM and DORT.

Radiative transfer of acoustic waves in continuous complex media: Beyond the Helmholtz equation.

Heterogeneity can be accounted for by a random potential in the wave equation. For acoustic waves in a fluid with fluctuations of both density and compressibility (as well as for electromagnetic waves in a medium with fluctuation of both permittivity and permeability) the random potential entails a scalar and an operator contribution. For simplicity, the latter is usually overlooked in multiple scattering theory: whatever the type of waves, this simplification amounts to considering the Helmholtz equation with a sound speed c depending on position r. In this work, a radiative transfer equation is derived from the wave equation, in order to study energy transport through a multiple scattering medium. In particular, the influence of the operator term on various transport parameters is studied, based on the diagrammatic approach of multiple scattering. Analytical results are obtained for fundamental quantities of transport theory such as the transport mean-free path ℓ^{*}, scattering phase function f, and anisotropy factor g. Discarding the operator term in the wave equation is shown to have a significant impact on f and g, yet limited to the low-frequency regime, i.e., when the correlation length of the disorder ℓ_{c} is smaller than or comparable to the wavelength λ. More surprisingly, discarding the operator part has a significant impact on the transport mean-free path ℓ^{*} whatever the frequency regime. When the scalar and operator terms have identical amplitudes, the discrepancy on the transport mean-free path is around 300% in the low-frequency regime, and still above 30% for ℓ_{c}/λ=10^{3} no matter how weak fluctuations of the disorder are. Analytical results are supported by numerical simulations of the wave equation and Monte Carlo simulations.

Scattering mean free path in continuous complex media: beyond the Helmholtz equation.

We present theoretical calculations of the ensemble-averaged (or effective or coherent) wave field propagating in a heterogeneous medium considered as one realization of a random process. In the literature, it is usually assumed that heterogeneity can be accounted for by a random scalar function of the space coordinates, termed the potential. Physically, this amounts to replacing the constant wave speed in Helmholtz' equation by a space-dependent speed. In the case of acoustic waves, we show that this approach leads to incorrect results for the scattering mean free path, no matter how weak the fluctuations. The detailed calculation of the coherent wave field must take into account both a scalar and an operator part in the random potential. When both terms have identical amplitudes, the correct value for the scattering mean free paths is shown to be more than 4 times smaller (13/3, precisely) in the low-frequency limit, whatever the shape of the correlation function. Based on the diagrammatic approach of multiple scattering, theoretical results are obtained for the self-energy and mean free path within Bourret's and on-shell approximations. They are confirmed by numerical experiments.

Measurements of ultrasonic diffusivity and transport speed from coda waves in a resonant multiple scattering medium.

Frequency-resolved experimental measurements of ultrasonic diffusivity in the MHz range are presented. The samples under study are two-dimensional random arrangements of parallel steel rods immersed in water and exhibit high-order multiple scattering. Their physical characteristics, particularly the density and pair-correlation functions of the scatterers, are well controlled. These synthetic samples are used as phantoms for actual inhomogeneous materials. The resonant nature of the scatterers has a strong effect on diffusivity, which is shown to vary significantly with frequency. This may affect the result of broadband measurements of apparent diffusivity, which can be expected to depend on time and sample thickness, whereas diffusivity is intrinsically an intensive parameter. Moreover, the transport speed is shown to vary drastically with frequency, sometimes by more than 50%, due to a very narrow resonance that slows down transport. Interestingly, this sharp resonance could only be revealed by experiments performed with coda waves, and not with ballistic or coherent waves whose frequency resolution is intrinsically limited from an experimental point of view.

Full transmission and reflection of waves propagating through a maze of disorder.

Thirty years ago, theorists showed that a properly designed combination of incident waves could be fully transmitted through (or reflected by) a disordered medium, based on the existence of propagation channels which are essentially either closed or open (bimodal law). In this Letter, we study elastic waves in a disordered waveguide and present direct experimental evidence of the bimodal law. Full transmission and reflection are achieved. The wave field is monitored by laser interferometry and highlights the interference effects that take place within the scattering medium.

Measurements of ultrasound velocity and attenuation in numerical anisotropic porous media compared to Biot's and multiple scattering models.

This article quantitatively investigates ultrasound propagation in numerical anisotropic porous media with finite-difference simulations in 3D. The propagation media consist of clusters of ellipsoidal scatterers randomly distributed in water, mimicking the anisotropic structure of cancellous bone. Velocities and attenuation coefficients of the ensemble-averaged transmitted wave (also known as the coherent wave) are measured in various configurations. As in real cancellous bone, one or two longitudinal modes emerge, depending on the micro-structure. The results are confronted with two standard theoretical approaches: Biot's theory, usually invoked in porous media, and the Independent Scattering Approximation (ISA), a classical first-order approach of multiple scattering theory. On the one hand, when only one longitudinal wave is observed, it is found that at porosities higher than 90% the ISA successfully predicts the attenuation coefficient (unlike Biot's theory), as well as the existence of negative dispersion. On the other hand, the ISA is not well suited to study two-wave propagation, unlike Biot's model, at least as far as wave speeds are concerned. No free fitting parameters were used for the application of Biot's theory. Finally we investigate the phase-shift between waves in the fluid and the solid structure, and compare them to Biot's predictions of in-phase and out-of-phase motions.

Coherent transmission of an ultrasonic shock wave through a multiple scattering medium.

We report measurements of the transmitted coherent (ensemble-averaged) wave resulting from the interaction of an ultrasonic shock wave with a two-dimensional random medium. Despite multiple scattering, the coherent waveform clearly shows the steepening that is typical of nonlinear harmonic generation. This is taken advantage of to measure the elastic mean free path and group velocity over a broad frequency range (2-15 MHz) in only one experiment. Experimental results are found to be in good agreement with a linear theoretical model taking into account spatial correlations between scatterers. These results show that nonlinearity and multiple scattering are both present, yet uncoupled.

Simulations of ultrasound propagation in random arrangements of elliptic scatterers: occurrence of two longitudinal waves.

Ultrasound propagation in clusters of elliptic (two-dimensional) or ellipsoidal (three-dimensional) scatterers randomly distributed in a fluid is investigated numerically. The essential motivation for the present work is to gain a better understanding of ultrasound propagation in trabecular bone. Bone microstructure exhibits structural anisotropy and multiple wave scattering. Some phenomena remain partially unexplained, such as the propagation of two longitudinal waves. The objective of this study was to shed more light on the occurrence of these two waves, using finite-difference simulations on a model medium simpler than bone. Slabs of anisotropic, scattering media were randomly generated. The coherent wave was obtained through spatial and ensemble-averaging of the transmitted wavefields. When varying relevant medium parameters, four of them appeared to play a significant role for the observation of two waves: (i) the solid fraction, (ii) the direction of propagation relatively to the scatterers orientation, (iii) the ability of scatterers to support shear waves, and (iv) a continuity of the solid matrix along the propagation. These observations are consistent with the hypothesis that fast waves are guided by the locally plate/bar-like solid matrix. If confirmed, this interpretation could significantly help developing approaches for a better understanding of trabecular bone micro-architecture using ultrasound.

Multiple scattering of ultrasound in weakly inhomogeneous media: application to human soft tissues.

Waves scattered by a weakly inhomogeneous random medium contain a predominant single-scattering contribution as well as a multiple-scattering contribution which is usually neglected, especially for imaging purposes. A method based on random matrix theory is proposed to separate the single- and multiple-scattering contributions. The experimental setup uses an array of sources/receivers placed in front of the medium. The impulse responses between every couple of transducers are measured and form a matrix. Single-scattering contributions are shown to exhibit a deterministic coherence along the antidiagonals of the array response matrix, whatever the distribution of inhomogeneities. This property is taken advantage of to discriminate single- from multiple-scattered waves. This allows one to evaluate the absorption losses and the scattering losses separately, by comparing the multiple-scattering intensity with a radiative transfer model. Moreover, the relative contribution of multiple scattering in the backscattered wave can be estimated, which serves as a validity test for the Born approximation. Experimental results are presented with ultrasonic waves in the megahertz range, on a synthetic sample (agar-gelatine gel) as well as on breast tissues. Interestingly, the multiple-scattering contribution is found to be far from negligible in the breast around 4.3 MHz.

Random matrix theory applied to acoustic backscattering and imaging in complex media.

The singular values distribution of the propagation operator in a random medium is investigated in a backscattering configuration. Experiments are carried out with pulsed ultrasonic waves around 3 MHz, using an array of transducers. Coherent backscattering and field correlations are taken into account. Interestingly, the distribution of singular values shows a dramatically different behavior in the single and multiple-scattering regimes. Based on a matrix separation of single and multiple-scattered waves, an experimental illustration of imaging through a highly scattering slab is presented.

Temperature-dependent diffusing acoustic wave spectroscopy with resonant scatterers.

The influence of a slight temperature change on the correlation of multiply scattered acoustic waves is studied, and experimental results are discussed. The technique presented here, similar to diffusing-acoustic-wave spectroscopy, is based on the sensitivity of a multiply scattering medium to a slight change. Ultrasonic waves around 3 MHz are transmitted through a sample made of steel rods in water and recorded by an array of transducers at different temperatures. The cross correlations between highly scattered signals are computed. The main effect of the temperature change is a simple dilation of the times of arrival, due to a change of the sound velocity in water. But the scatterers also play a role in the progressive decorrelation of wave forms. An analysis resolved in both time and frequency shows that at some particular frequencies, the resonant behavior of the scatterers is responsible for a significantly larger decorrelation. Interestingly, the experimental results allow one to detect the presence of a small resonance that was not detected earlier on the same scatterers with classical measurement of the scattering mean free path. A simple model is proposed to interpret the experimental results.

Ultrasonic imaging of highly scattering media from local measurements of the diffusion constant: separation of coherent and incoherent intensities.

As classical imaging fails with diffusive media, one way to image a multiple-scattering medium is to achieve local measurements of the dynamic transport properties of a wave undergoing diffusion. This paper presents a method to obtain local measurements of the diffusion constant D in a multiple-scattering medium. The experimental setup consists in an array of programmable transducers placed in front of the multiple-scattering medium to be imaged. By achieving Gaussian beamforming both at emission and reception, an array of virtual sources and receivers located in the near field is constructed. The time evolution of the incoherent component of the intensity backscattered on this virtual array is shown to represent directly the growth of the diffusive halo as sqrt[Dt]. A matrix treatment is proposed to separate the incoherent intensity from the coherent backscattering peak. Once the incoherent contribution is isolated, a local measurement of the diffusion constant is possible. The technique is applied to image the long-scale variations of D in a random-scattering sample made of two parts with a different concentration of cylindrical scatterers. This experimental result is obtained with ultrasonic waves around 3 MHz. It illustrates the possibility of imaging diffusive media from local measurements of the diffusion constant, based on coherent Gaussian beamforming and a matrix "antisymmetrization," which creates a virtual antireciprocity.

Influence of correlations between scatterers on the attenuation of the coherent wave in a random medium.

Experimental measurements of the coherent wave transmission for ultrasonic waves propagating in water through a random set of scatterers (metallic rods) are presented. Though the densities are moderate (6% and 14%) the experimental results show that the mean-free path deviates from the classical first-order approximation due to the existence of correlations between scatterers. Theoretical results for the mean free path obtained from different approaches are compared to the experimental measurements. The best agreement is obtained with the second-order diagrammatic expansion of the self-energy.

Passive retrieval of Rayleigh waves in disordered elastic media.

When averaged over sources or disorder, cross correlation of diffuse fields yields the Green's function between two passive sensors. This technique is applied to elastic ultrasonic waves in an open scattering slab mimicking seismic waves in the Earth's crust. It appears that the Rayleigh wave reconstruction depends on the scattering properties of the elastic slab. Special attention is paid to the specific role of bulk to Rayleigh wave coupling, which may result in unexpected phenomena, such as a persistent time asymmetry in the diffuse regime.

Weak localization and time reversal of ultrasound in a rotational flow.

A one-channel time-reversal (TR) experiment is performed inside a rough reverberating cavity in the presence of a rotational flow. The amplitude of the TR wave is plotted versus the distance between the TR channel and the initial source: when they coincide, it exhibits a "time-reversal enhancement" (TRE). With no flow, the TRE is the same as the coherent backscattering enhancement (CBE). But contrary to CBE, the TRE peak is shown to be insensitive to the breaking down of reciprocity due to the flow. This new effect of weak localization is sustained by a diagrammatic approach.

Relation between time reversal focusing and coherent backscattering in multiple scattering media: a diagrammatic approach.

In this paper, we revisit one-channel time reversal (TR) experiments through multiple scattering media in the framework of the multiple scattering theory. The hyperresolution and the self-averaging property are retrieved. The developed formalism leads to a deeper understanding of the role of the ladder and most-crossed diagrams in a TR experiment and also establishes the link between TR and coherent backscattering (CBS). Especially, we show that when the initial source and the time reversal point are at the same location, the time-reversed amplitude is twice higher. Surprisingly, this enhancement is due to the ladder diagrams and not to the most-crossed ones, contrary to CBS. These theoretical predictions are confirmed by experimental results. The experiments are performed with ultrasonic waves propagating through a random collection of parallel steel rods.

Taking advantage of multiple scattering to communicate with time-reversal antennas.

We present an experimental demonstration showing that, contrary to first intuition, the more scattering a mesoscopic medium is, the more information can be conveyed through it. We used a multiple input-multiple output configuration: a multichannel ultrasonic time-reversal antenna is used to transmit random series of bits simultaneously to different receivers which were only a few wavelengths apart. Whereas the transmission is free of error when multiple scattering occurs in the propagation medium, the error rate is huge in a homogeneous medium.

Time reversal versus phase conjugation in a multiple scattering environment.

We present experimental results on the reversibility of ultrasound in a multiple scattering medium. An ultrasonic pulsed wave is transmitted from a point source to a 128-element receiving array through 2D samples with various thickness. The samples consist of random collections of parallel steel rods immersed in water. The scattered waves are recorded, time reversed and sent back into the medium. The time-reversed waves are converging back to their source and the quality of spatial and temporal focusing on the source is related to the second-order moments of the scattered wave (correlation) in time and in space. Experimental results show that it is possible to obtain a robust estimation of the correlations on a single realisation of disorder, taking advantage of the wide frequency bandwidth. The spatial resolution of the system is only limited by the correlation length of the scattered field, and no longer by diffraction. Moreover, successful time-reversal focusing using a single element instead of an array is possible, whereas a one-channel monochromatic phase conjugation fails. The efficiency of broad-band time reversal compared to monochromatic phase conjugation lies in the number of 'information grains' in the frequency bandwidth.