Citations

Elasticity imaging methods have been used to study tissue mechanical properties and have demonstrated that tissue elasticity changes with disease state. Quantitative mechanical properties can be measured in a model independent manner if both shear wave speed and attenuation are known. However, measuring shear wave speed attenuation is challenging in the field of elasticity imaging. Typically, only shear wave speed is measured and rheological models, such as Kelvin-Voigt, Maxwell and Standard Linear Solid, are used to solve for shear viscoelastic complex modulus. Acoustic radiation force has been used to study quasi-static viscoelastic properties of tissue during creep and relaxation conditions, however, as with shear wave propagation methods, a rheological model needs to be fit to the creep or relaxation experimental data to solve for viscoelastic parameters. This paper presents a method to quantify viscoelastic properties in a model-independent way by estimating complex shear elastic modulus over a wide frequency range using time-dependent creep response induced by acoustic radiation force. The acoustic radiation force induced creep (RFIC) method uses a conversion formula that is the analytic solution of the constitutive equation relating time dependent stress and time dependent strain. The RFIC method in combination with shear wave propagation is used to measure the complex shear modulus so that knowledge of the applied radiation force magnitude is not necessary. Numerical simulation of creep strain and compliance using the Kelvin-Voigt model shown that the conversion formula is sensitive to sampling frequency, the first reliable measure in time and the long term viscosity approximation. Experimental data are obtained in homogeneous tissue mimicking phantoms and excised swine kidneys.

This study concerns the development and preliminary experimental verification of a topological sensitivity–based platform for the material characterization of tissue anomalies exposed by vibroacoustography-type imaging techniques. Vibroacoustography (VA) is a high-resolution imaging method that has been applied to the detection of pathological changes in soft tissues. Although the data provided by this method is related to the mechanical properties of tissue, the viscoelastic parameters of the object cannot be estimated by this imaging method itself. Topological sensitivity (TS) method is a data processing methodology that can be used to estimate the viscoelastic parameters of an object from vibration data. In this study, the concept of topological sensitivity is applied to interpret the vibroacoustography measurements for the purpose of lesion characterization. In the proposed approach, the topological sensitivity function, which signifies the variation of a given cost functional when an infinitesimal inclusion with trial material parameters is placed at the location of a point force, is formulated in terms of the adjoint field. The effectiveness of the resulting formula as a material indicator for lesion characterization is demonstrated by estimating the relative elastic parameters of a well-controlled neoprene sphere embedded in a tissue-mimicking phantom specimen.

Ultrasonic contrast agents consisting of microbubbles are used to assess tissue perfusion. The microbubbles are highly reflective and nonlinear and thus produce harmonics that are stronger than those from tissues. However, the magnitude of harmonic signals resulting from a region with microbubbles also depends on the acoustic pressure of incident ultrasound and the attenuation of intervening tissues in the ultrasound path. Therefore, the harmonic magnitude, as used in traditional harmonic imaging, may not be a reliable indicator of the presence or absence of microbubbles, and hence, tissue perfusion. To compensate for these effects, we present two parameters defined as the ratio of the harmonic to the fundamental component ( HFR) and the ratio of the harmonic to squared fundamental ( HSFR). A simplified model is used to illustrate the usefulness of these two parameters. Experiments show that both parameters improve detection of microbubbles and that HSFR performs better than HFR.

Vibro-acoustography is a method for imaging the acoustic energy emitted from objects in response to an oscillatory radiation force produced by two interfering focused beams of ultrasound [M. Fatemi and J. F. Greenleaf, Science 280, 82–85 (1998)]. To facilitate quantitative study, a model is presented that describes the behavior of an elastic sphere in two plane ultrasound waves of slightly different frequencies. The vibrating velocity and the resulting acoustic emission of the sphere at the difference frequency are derived from the radiation force and the impedance of the sphere. For small steel spheres in water, the velocity and emission are predicted to have a simple relationship with the difference frequency and the size of the sphere. Experiments on small stainless steel spheres of different size are carried out for a number of difference frequencies, with a laser interferometer to measure the velocity of the sphere and a calibrated hydrophone to measure the acoustic emission from the sphere. The experimental results fit the model description well. This means that vibro-acoustography may be useful in quantitative measurement.

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.

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.

The acoustic radiation force (ARF), which signifies the average transfer of momentum from the sound wave to the propagating medium, is nowadays frequently used to compute the mean displacement field in soft tissues due to the action of a focused ultrasound beam. In situations when the sound field used to generate the ARF is modulated, the latter depends on both spatial and temporal derivatives of the acoustic intensity. Unfortunately, the numerical computation of the nonlinear acoustic solution due to modulated focused ultrasound beam is complicated by the presence of two disparate time scales in the formulation. To deal with the problem, the KZK-type spatial scaling is complemented by its dual-time-scale companion, where the temporal coordinate is split into its “fast” and “slow” components, allowing one to track the ultrasound-scale oscillations and modulation-scale variations separately in the solution. In this case, the nonlinear acoustic solution (and consequently the ARF) can be effectively computed in the “fast” frequency domain and “slow” time domain for any given modulation envelope, transient, or steady-state. The proposed developments are both validated by the experimental results and illustrated via numerical examples that show the effectiveness of the computational scheme.

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.