Greenleaf, James F.

New vibration pulses are proposed to increase the power of shear waves induced by ultrasound radiation force in a tissue region with a preferred spectral distribution. The new pulses are sparsely sampled from an orthogonal frequency wave composed of several sinusoidal signals. Those sinusoidal signals have different frequencies and are orthogonal to each other. The phase and amplitude of each sinusoidal signal are adjusted to control the shape of the orthogonal frequency wave. Amplitude of the sinusoidal signal is increased as its frequency increases to compensate for higher loss at higher frequency in the tissue region. The new vibration pulses and detection pulses can be interleaved for array transducer applications. The experimental results show that the new vibration pulses significantly increases induced tissue vibration with the same peak ultrasound intensity, compared with the binary vibration pulses.

In the past several decades, the fields of ultrasound and magnetic resonance elastography have shown promising results in noninvasive estimates of mechanical properties of soft tissues. These techniques often rely on measuring shear wave velocity due to an external or internal source of force and relating the velocity to viscoelasticity of the tissue. The mathematical relationship between the measured velocity and material properties of the myocardial wall, arteries, and other organs with non-negligible boundary conditions is often complicated and computationally expensive. A simple relationship between the Lamb–Rayleigh dispersion and the shear wave dispersion is derived for both the velocity and attenuation. The relationship shows that the shear wave velocity is around 20% higher than the Lamb–Rayleigh velocity and that the shear wave attenuation is about 20% lower than the Lamb–Rayleigh attenuation. Results of numerical simulations in the frequency range 0–500 Hz are presented.

Tissue mechanical properties can be obtained by algebraic direct inversion (ADI) of the shear-wave motion differential equation, which is insensitive to reflection and geometry of the pushing beam. Shear waves normally generated by a focused ultrasound pushing beam have limited spatial extent in depth and are transient in time, leading to noisy and unstable ADI results. A comb-shape acoustic radiation force distribution can generate shear waves with longer spatial extent and time duration, facilitating more robust ADI. For 400 μs, a linear array transducer simultaneously transmitted four unfocused pushing beams (12 on elements for each pushing beam, 8 off elements between beams) into a tissue-mimicking phantom (∼1.5 kPa) with a cylindrical inclusion (∼4 kPa). The ultrasound system (Verasonics Inc.) then immediately switched to flash imaging mode (frame rate = 2 kHz, spatial resolution = 0.31 mm) to measure shear-wave displacements in a 38 mm by 30 mm domain that was used for ADI. The reconstructed 2D shear elasticity map provides accurate shear elasticity estimates (background = 1.5 kPa with 7% variance; inclusion = 3.9 kPa with 15% variance) and excellent contrast between the background and inclusion.

In the field of tissue engineering, monitoring the elasticity of developing tissue cultures is important for evaluating their growth and maturation. In this work, we used wave propagation methods to measure the elasticity of thin tissue-mimicking layers to test the feasibility of measuring the elasticity of developing tissue cultures. We developed finite element and analytical models for investigating this problem and found good agreement between the results from both models. We used ultrasound radiation force to induce propagating waves, and the waves are measured using high frame rate ultrasound imaging at 10 kHz. The wave modes were identified and compared with an analytic model for Rayleigh wave propagation in a thin elastic layer bound to a solid substrate. The same gelatin mixture was used to make phantoms 1, 4, and 50 mm thick, where the phantom with 50 mm thickness was used as control and evaluated with shear wave imaging methods. The analytic model was used to fit data from experiments in the large block and thin layers of 1 and 4 mm thick, and the measured shear moduli were 22.7, 21.8, and 22.5 kPa, respectively. These results provide a validation for measurement of elasticity in thin elastic layers.

Propagating shear waves in tissue can be measured with high frame rate Doppler or correlation methods. The measured characteristics of the shear waves, such as speed versus frequency, can be used to deduce material properties such as complex viscoelastic modulus using physics models appropriate to the geometry and properties of the tissue. This inverse problem is characterized by calculating the storage and loss modulus as a function of frequency and requires appropriate tissue motion, which in turn requires optimized dynamic radiation force distributions. We will discuss novel radiation force distributions that provide enhanced tissue motions appropriate to dealing with the inverse problem of determining tissue material properties from ultrasonically measured tissue motion.

The response of an embedded sphere in a viscoelastic medium excited by acoustic radiation force has been studied in both the time- and frequency-domains. This model is important because it can be used to characterize the viscoelastic properties of the medium by fitting the response to the theoretical model. The Kelvin–Voigt model has been used exclusively in these models. An extension to the previously reported models is described so that any viscoelastic rheological model can be used. This theoretical development describes the generalized embedded sphere response both in the time and frequency domains. Comparing the results from derivations in both domains showed very good agreement with a median absolute error (MAE) ranging from 0.0044 to 0.0072. Good agreement is demonstrated with finite element model simulations and the theory with a MAE of 0.006. Lastly, results for characterization of gelatin and rubber materials with the new theory are shown where the MAE values were used to determine which rheological model best describes the measured responses.

Characterization of viscoelastic soft tissues and materials has been a subject of active research for the past two decades. Techniques involving measuring the motion of embedded spheres in viscoelastic media have been employed to measure the viscoelastic material properties of the media surrounding the sphere. The published work using this technique has been restricted to using the Kelvin–Voigt model, which has worked for gelatin-based materials. We report an extension of this method using different viscoelastic rheological models such as the Maxwell, generalized Maxwell, Zener, and Kelvin–Voigt fractional derivative models. We performed experiments using a focused 3-MHz transducer to produce radiation force and a 7.5-MHz pulse-echo transducer to measure the sphere motion. Two phantoms were tested; one made of gelatin and another made of rubber with stainless steel spheres embedded in the materials. The response was fitted with the generalized theoretical model to find the material properties over bandwidths of 25–3000 Hz. The median absolute errors ranged from 0.0036–0.0421 and 0.0365–0.0879 for the gelatin and rubber phantoms, respectively. This generalized method offers added flexibility to characterize different types of viscoelastic materials using the most appropriate rheological model. [Work supported in part by NIH Grant nos. EB002167 and DK082408.]

Multiple methods have been proposed to characterize the viscoelastic material properties of soft tissues to differentiate normal from pathological tissue. Shearwave dispersion ultrasound vibrometry is a method that uses the dispersion of the shear wave velocity to characterize tissue shear elasticity and viscosity. It has recently been shown [M. W. Urban and J. F. Greenleaf, Phys. Med. Biol. 54, 5919–5933 (2009)] that the Kramers–Kronig relationships can be applied to shear wave measurements in ex vivo skeletal muscle over a frequency bandwidth of 100–600 Hz. We will describe this process where shear wave attenuation is modeled using a power law, and the shear wave velocity dispersion is predicted by applying the Kramers–Kronig relationships. In this study, we apply this method to shear wave propagation measurements made in different soft tissues such as ex vivo porcine skeletal muscle and kidney and in vivo porcine liver. The shear wave velocity predictions using the Kramers–Kronig relationships agree with the velocity measurements. The mean square error for these tissues ranged from 0.0024–0.0625 m2/s2. The Kramers–Kronig relationships can be used in shear wave propagation measurements in different types of soft tissues. [This work was supported in part by NIH Grants Nos. EB002167 and DK082408.]

Dynamic radiation force has been used in several types of applications, and is performed by modulating ultrasound with different methods. By modulating ultrasound, energy can be transmitted to tissue, in this case a dynamic force to elicit a low frequency cyclic displacement to inspect the material properties of the tissue. In this paper, different types of modulation are explored including amplitude modulation (AM), double sideband suppressed carrier amplitude modulation AM, linear frequency modulation, and frequency-shift keying. Generalized theory is presented for computing the radiation force through the short-term time average of the energy density for these various types of modulation. Examples of modulation with different types of signals including sine waves, square waves, and triangle waves are shown. Using different modulating signals, multifrequency radiation force with different numbers of frequency components can be created, and can be used to characterize tissue mimicking materials and soft tissue. Results for characterization of gelatin phantoms using a method of vibrating an embedded sphere are presented. Different degrees of accuracy were achieved using different modulation techniques and modulating signals. Modulating ultrasound is a very flexible technique to produce radiation force with multiple frequency components that can be used for various applications.