Time reversal of ultrasonic fields. I. Basic principles

Time reversal of ultrasonic fields represents a way to focus through an inhomogeneous medium. This may be accomplished by a time-reversal mirror (TRM) made from an array of transmit-receive transducers that respond linearly and allow the incident acoustic pressure to be sampled. The pressure field is then time-reversed and re-emitted. This process can be used to focus through inhomogeneous media on a reflective target that behaves as an acoustic source after being insonified. The time-reversal approach is introduced in a discussion of the classical techniques used for focusing pulsed waves through inhomogeneous media (adaptive time-delay techniques). Pulsed wave time-reversal focusing is shown using reciprocity valid in inhomogeneous medium to be optimal in the sense that it realizes the spatial-temporal matched filter to the inhomogeneous propagation transfer function between the array and the target. The research on time-reversed wave fields has also led to the development of new concepts that are described: time-reversal cavity that extends the concept of the TRM, and iterative time-reversal processing for automatic sorting of targets according to their reflectivity and resonating of extended targets.     

Narrowband shear wave generation by a Finite-Amplitude radiation force: The fundamental component

A highly localized source of low-frequency shear waves can be created by the modulated radiation force resulting from two intersecting quasi-continuous-wave ultrasound beams of slightly different frequencies. In contrast to most other radiation force-based methods, these shear waves can be narrowband. Consequently, different frequency-dependent effects will not significantly affect their spectrum as they propagate within a viscoelastic medium, thereby enabling the viscoelastic shear properties of the medium to be determined at any given modulation frequency. This can be achieved by tracking the shear wave phase delay and change in amplitude over a specific distance. In this paper we explore the properties of short duration (dynamic) low-frequency shear wave propagation and study how the shear displacement field depends on the excitation conditions. Our investigations make use of the approximate Green's functions for viscoelastic media, and the evolution of such waves is studied in the spatiotemporal domain from a theoretical perspective. Although nonlinearities are included in our confocal source model, just the properties of the fundamental shear component are examined in this paper. We examine how the shear wave propagation is affected by the shear viscosity, the coupling wave, the spatial distribution of the force, the shear speed, and the duration of the modulated wave. A method is proposed for estimating the shear viscosity of a viscoelastic medium. In addition, it is shown how the Voigt model paremeters can be extracted from the frequency-dependent speed and attenuation.     

Ultrasonic backscatter imaging by shear-wave-induced echo phase encoding of target locations

We present a novel method for ultrasound backscatter image formation wherein lateral resolution of the target is obtained by using traveling shear waves to encode the lateral position of targets in the phase of the received echo. We demonstrate that the phase modulation as a function of shear wavenumber can be expressed in terms of a Fourier transform of the lateral component of the target echogenicity. The inverse transform, obtained by measurements of the phase modulation over a range of shear wave spatial frequencies, yields the lateral scatterer distribution. Range data are recovered from time of flight as in conventional ultrasound, yielding a B-mode-like image. In contrast to conventional ultrasound imaging, where mechanical or electronic focusing is used and lateral resolution is determined by aperture size and wavelength, we demonstrate that lateral resolution using the proposed method is independent of the properties of the aperture. Lateral resolution of the target is achieved using a stationary, unfocused, single-element transducer. We present simulated images of targets of uniform and non-uniform shear modulus. Compounding for speckle reduction is demonstrated. Finally, we demonstrate image formation with an unfocused transducer in gelatin phantoms of uniform shear modulus.     

Time-reversed Lamb waves

Lamb waves are extensively involved in plate structure inspection because of their guided nature. However, their dispersive nature often limits their use in flaw detection. In this paper we show that the use of a time-reversal mirror (TRM) allows to automatically compensate for the dispersive nature of Lamb waves. Experiments showing the spatial and temporal behavior of time-reversed Lamb waves, demonstrate the ability of TRMs to self-focus and to recompress dispersive pulses. This is demonstrated in a set of experiments in which a broadband ultrasonic laser source is used to simulate a point Lamb wave source and an optical interferometer is used to map the time reversed elastic field. We also show that TRM may work in pulse echo mode and allows to detect and to focus along large 2-D plates on any flaws located in the inspected area.     

3-D FDTD simulation of shear waves for evaluation of complex modulus imaging

The Navier equation describing shear wave propagation in 3-D viscoelastic media is solved numerically with a finite differences time domain (FDTD) method. Solutions are formed in terms of transverse scatterer velocity waves and then verified via comparison to measured wave fields in heterogeneous hydrogel phantoms. The numerical algorithm is used as a tool to study the effects on complex shear modulus estimation from wave propagation in heterogeneous viscoelastic media. We used an algebraic Helmholtz inversion (AHI) technique to solve for the complex shear modulus from simulated and experimental velocity data acquired in 2-D and 3-D. Although 3-D velocity estimates are required in general, there are object geometries for which 2-D inversions provide accurate estimations of the material properties. Through simulations and experiments, we explored artifacts generated in elastic and dynamic-viscous shear modulus images related to the shear wavelength and average viscosity.     

Some exact solutions for the propagation of transient electroacoustic waves. I: Piezoelectric half-space

In piezoelectric materials, coupled electromagnetic and horizontally polarized shear surface wave disturbances can exist which have no purely elastic counterpart. The properties of these electroacoustic surface waves and similarly coupled body wave components have to date only been described by far-field approximations in the frequency domain. This article obtains exact transient solutions for the electroacoustic surface and body waves generated by a dipole source on a piezoelectric half-space. These solutions are obtained for both conducting and nonconducting surface boundary conditions using a modification of the Lamb-Cagniard-Pekeris technique previously applied to similar elastic and acoustic wave propagation problems. Explicit results for the separated surface and body wave contributions at the surface of the piezoelectric are given and discussed.     

基于时间反转镜的宽带匹配场处理

声场的互易性和线性波动方程的时反不变性能用时间反转来代替频域的相位共轭。因此,时间反转技术可以应用于宽带的匹配场处理中。时间反转匹配场处理时,仍然需要知道海洋环境的先验知识,例如声场的声速分布,海洋边界条件等等。文章对时间反转应用于匹配场处理进行了仿真并分析海上试验数据。仿真和试验数据分析结果表明:该套时间的反转处理方法可以实现匹配场定位,定位精度较高,而且多元阵信号处理的增益也十分明显。     

Single Mode Tuning Effects on Lamb Wave Time Reversal with Piezoelectric Wafer Active Sensors for Structural Health Monitoring

Lamb wave time reversal method is a new and tempting baseline-free damage detection technique for structural health monitoring. With this method, certain types of damage can be detected without baseline data. However, the application of this method using piezoelectric wafer active sensors (PWAS) is complicated by the existence of at least two Lamb wave modes at any given frequency, and by the dispersion nature of the Lamb wave modes existing in thin-wall structures. The theory of PWAS-related Lamb wave time reversal has not yet been fully studied. This paper addresses this problem by developing a theoretical model for the analysis of PWAS-related Lamb wave time reversal based on the exact solutions of the Rayleigh-Lamb wave equation. The theoretical model is first used to predict the existence of single-mode Lamb waves. Then the time reversal behavior of single-mode and two-mode Lamb waves is studied numerically. The advantages of single-mode tuning in the application of time reversal damage detection are highlighted. The validity of the proposed theoretical model is verified through experimental studies. In addition, a similarity metric for judging time invariance of Lamb wave time reversal is presented. It is shown that, under certain condition, the use of PWAS-tuned single-mode Lamb waves can greatly improve the effectiveness of the time-reversal damage detection procedure.     

管桩屏障在柱面波源下近场隔振效应的解析求解

基于波函数展开法表示柱面波源自由场,首次提出了一种求解柱面波源下管桩屏障对弹性波隔振效应的解析方法。方法首先考虑柱面波源与管桩屏障的位置关系,采用波函数展开法对入射柱面波进行了0阶Hankel展开表示,并采用任意坐标系间变换的Graf加法定理在任意桩体坐标系中表示,将入射波场和散射波场叠加后通过满足所有桩体的边界条件以求解所有桩体的散射波场。该解析方法分析了柱面波源入射下群桩屏障的隔振效果,首次在排桩隔振问题中考虑了入射波曲率的影响,为柱面波的散射问题提供了理论解答。该文重点讨论了入射波曲率、管桩个数和桩排数等因素对管桩隔振效果的影响,结果表明：整体上管桩屏障对柱面P波的隔振效果优于柱面SV波;柱面波源距离管桩屏障更近时排桩后场地的位移反应显著增大;相比平面波,柱面波源作用下排桩数量的提升对隔振效果的影响较小;三排管桩屏障比两排管桩屏障隔振效果更强,更宜采用三排管桩屏障进行柱面波源的隔振。     

Supersonic shear imaging: a new technique for soft tissue elasticity mapping

Supersonic shear imaging (SSI) is a new ultrasound-based technique for real-time visualization of soft tissue viscoelastic properties. Using ultrasonic focused beams, it is possible to remotely generate mechanical vibration sources radiating low-frequency, shear waves inside tissues. Relying on this concept, SSI proposes to create such a source and make it move at a supersonic speed. In analogy with the "sonic boom" created by a supersonic aircraft, the resulting shear waves will interfere constructively along a Mach cone, creating two intense plane shear waves. These waves propagate through the medium and are progressively distorted by tissue heterogeneities. An ultrafast scanner prototype is able to both generate this supersonic source and image (5000 frames/s) the propagation of the resulting shear waves. Using inversion algorithms, the shear elasticity of medium can be mapped quantitatively from this propagation movie. The SSI enables tissue elasticity mapping in less than 20 ms, even in strongly viscous medium like breast. Modalities such as shear compounding are implementable by tilting shear waves in different directions and improving the elasticity estimation. Results validating SSI in heterogeneous phantoms are presented. The first in vivo investigations made on healthy volunteers emphasize the potential clinical applicability of SSI for breast cancer detection.     

Dynamic Stress around Two Interacting Cylindrical Nano-Inhomogeneities with Surface/Interface Effects

On the basis of continuum surface elasticity, two interacting cylindrical nano-inhomogeneities with surface/interface effect in a small-sized solid under anti-plane shear waves are investigated, and the dynamic stress around the nano-inhomogeneities is analyzed. The wave function expansion method is used to expressed the wave field around the two nano-inhomogeneities. The total wave field is obtained by the addition theorem for cylindrical wave function. Through analysis, it is found that the distance between the two nano-inhomogeneities shows great effect on the dynamic stress in nano composites. The effect of the distance is also related to the properties of the nano-inhomogeneities and the interface, the wave frequency, and the incident angle of shear waves. To show the accuracy of the results for certain given parameters, comparison with the existing results is also given.     

Ultrasonic measurements on hydrating cement slurries: Onset of shear wave propagation

The ultrasonic properties of three oil field cement slurries are studied during the early stages of the hydration process. As a percolating solid framework is established, the slurries develop mechanical integrity and shear waves begin to propagate. This transition is also apparent in the behavior of the compressional wave, but is considerably more clear-cut in the shear signals. For the three samples examined here, the ratios of the shear wave onset times are in good agreement with the corresponding ratios of the American Petroleum Institute (API) thickening times.     

Monitoring Early Age Properties of Cementitious Material Using Ultrasonic Guided Waves in Embedded Rebar

The evaluation of early age properties of concrete is critical for ensuring construction quality. This paper presents a sensing method to use ultrasonic guided waves in a rebar for monitoring the early age properties of cementitious material. An EMAT sensor was used to excite the longitudinal mode L(0,1) wave in a rebar embedded in cement/mortar, and an ultrasonic transducer was used for receiving the echo signals. Guided wave dispersion curves were developed to select appropriate frequency range. The leakage attenuations of the L(0,1) mode wave from the rebar to the surrounding cement materials were continuously monitored for the first 10 h. The evolution of the shear wave velocity was also monitored simultaneously. The leakage attenuation from experimental measurements was compared with the theory-predicted attenuation in both time and frequency domains, and showed good agreement. Experiments were performed on three cement paste samples and three mortar samples. The results indicated that attenuation is nearly linearly related to the shear wave velocity, and shear wave velocity is linearly related to the penetration resistance (ASTM C403) in logarithmic scale. These results suggest that mechanical properties and hardening process of cement materials can be monitored by using the ultrasonic guided waves in a rebar.     

The role of viscosity in the impulse diffraction field of elastic waves induced by the acoustic radiation force

Several ultrasound-based techniques for the estimation of soft tissue elasticity are currently being investigated. Most of them study the medium response to dynamic excitations. Such responses are usually modeled in a purely elastic medium using a Green's function solution of the motion equation. However, elasticity by itself is not necessarily a discriminant parameter for malignancy diagnosis. Modeling viscous properties of tissues could also be of great interest for tumor characterization. We report in this paper an explicit derivation of the Green's function in a viscous and elastic medium taking into account shear, bulk, and coupling waves. From this theoretical calculation, 3D simulations of mechanical waves in viscoelastic soft tissues are presented. The relevance of the viscoelastic Green's function is validated by comparing simulations with experimental data. The experiments were conducted using the supersonic shear imaging (SSI) technique which dynamically and remotely excites tissues using acoustic radiation force. We show that transient shear waves generated with SSI are modeled very precisely by the Green's function formalism. The combined influences of out-of-plane diffraction, beam shape, and shear viscosity on the shape of transient waves are carefully studied as they represent a major issue in ultrasound-based viscoelasticity imaging techniques.     

Propagation of shear waves generated by a modulated finite amplitude radiation force in a viscoelastic medium

An effective way to generate localized narrow-band low-frequency shear waves within tissue noninvasively, is by the modulated radiation force, resulting from the interference of two confocal quasi-CW ultrasound beams of slightly different frequencies. By using approximate viscoelastic Green's functions, investigations of the properties of the propagated shear-field component at the fundamental modulation frequency were previously reported by our group. However, high-amplitude source excitations may be needed to increase the signal-to-noise-ratio for shear-wave detection in tissue. This paper reports a study of the generation and propagation of dynamic radiation force components at harmonics of the modulation frequency for conditions that generally correspond to diagnostic safety standards. We describe the propagation characteristics of the resulting harmonic shear waves and discuss how they depend on the parameters of nonlinearity, focusing gain, and absorption. For conditions of high viscosity (believed to be characteristic of soft tissue) and higher modulation frequencies, the approximate shear wave Green's function is inappropriate. A more exact viscoelastic Green's function is derived in k-space, and using this, it is shown that the lowpass and dispersive effects, associated with a Voigt model of tissue, are more accurately represented. Finally, it is shown how the viscoelastic properties of the propagating medium can be estimated, based on several spectral components of the shear wave spectrum.     

Diffraction field of a low frequency vibrator in soft tissues using transient elastography

For the last 10 years, interest has grown in low frequency shear waves that propagate in the human body. However, the generation of shear waves by acoustic vibrators is a relatively complex problem, and the directivity patterns of shear waves produced by the usual vibrators are more complicated than those obtained for longitudinal ultrasonic transducers. To extract shear modulus parameters from the shear wave propagation in soft tissues, it is important to understand and to optimize the directivity pattern of shear wave vibrators. This paper is devoted to a careful study of the theoretical and the experimental directivity pattern produced by a point source in soft tissues. Both theoretical and experimental measurements show that the directivity pattern of a point source vibrator presents two very strong lobes for an angle around 35 degrees . This paper also points out the impact of the near field in the problem of shear wave generation.     

Experimental and theoretical study of Rayleigh-Lamb waves in a plate containing a surface-breaking crack

Scattering of Rayleigh-Lamb waves by a normal surface-breaking crack in a plate has been studied both theoretically and experimentally. The two-dimensionality of the far field, generated by a ball impact source, is exploited to characterize the source function using a direct integration technique. The scattering of waves generated by this impact source by the crack is subsequently solved by employing a Green's function integral expression for the scattered field coupled with a finite element representation of the near field. It is shown that theoretical results of plate response, both in frequency and time, are similar to those obtained experimentally. Additionally, implications for practical applications are discussed.     

Propagation of QSH (quasi shear horizontal) acoustic waves in piezoelectric plates

The characteristics of QSH (quasi shear horizontal) acoustic waves propagating in thin plates of Y-cut, X-propagation lithium niobate are investigated theoretically and experimentally. The fractional velocity change (Deltanu/nu) produced by electrical shorting of the surface is calculated as a function of the normalized plate thickness h/lambda (h=plate thickness, lambda=acoustic wavelength). It was found that values of Deltanu/nu as high as 0.18 could be obtained. Experimental measurements show good agreement with theory. The properties of QSH waves propagating in the presence of a perfectly conducting electrode separated from the piezoelectric plate by a small air gap have been studied theoretically and experimentally. It was found that by varying the height of the gap, the phase shift through a 3.2-MHz QSH wave delay line can be varied by more than 230 degrees . We have also theoretically investigated the influence of a thin layer of arbitrary conductivity on the velocity and attenuation of the QSH wave. Calculations show that the variations in these parameters can be as high as 18% and 5 dB per wavelength for a change in layer surface conductance from 10(-7) to 10(-5) S. Results obtained in this paper confirm the attractive properties of QSH waves for a variety of sensing and signal processing applications.     

Characterization of roll-drawn polypropylene using shear wave birefringence

Ultrasonic shear waves have been used to estimate the degree of anisotropy present within roll-drawn polypropylene samples, using shear wave birefringence. These measurements have been correlated with the physical properties that result from the draw ratio used during manufacture. It is shown that the technique has promise as a monitoring tool for this material.     

Dynamic mechanical response of elastic spherical inclusions to impulsive acoustic radiation force excitation

Acoustic radiation force impulse imaging has been used clinically to study the dynamic response of lesions relative to their background material to focused, impulsive acoustic radiation force excitations through the generation of dynamic displacement field images. Dynamic displacement data are typically displayed as a set of parametric images, including displacement immediately after excitation, maximum displacement, time to peak displacement, and recovery time from peak displacement. To date, however, no definitive trends have been established between these parametric images and the tissues' mechanical properties. This work demonstrates that displacement magnitude, time to peak displacement, and recovery time are all inversely related to the Young's modulus in homogeneous elastic media. Experimentally, pulse repetition frequency during displacement tracking limits stiffness resolution using the time to peak displacement parameter. The excitation pulse duration also impacts the time to peak parameter, with longer pulses reducing the inertial effects present during impulsive excitations. Material density affects tissue dynamics, but is not expected to play a significant role in biological tissues. The presence of an elastic spherical inclusion in the imaged medium significantly alters the tissue dynamics in response to impulsive, focused acoustic radiation force excitations. Times to peak displacement for excitations within and outside an elastic inclusion are still indicative of local material stiffness; however, recovery times are altered due to the reflection and transmission of shear waves at the inclusion boundaries. These shear wave interactions cause stiffer inclusions to appear to be displaced longer than the more compliant background material. The magnitude of shear waves reflected at elastic lesion boundaries is dependent on the stiffness contrast between the inclusion and the background material, and the stiffness and size of the inclusion dictate when shear wave reflections within the lesion will interfere with one another. Jitter and bias associated with the ultrasonic displacement tracking also impact the estimation of a tissue's dynamic response to acoustic radiation force excitation.