Shear Wave Motion Research Articles

Characterizing brain mechanics through 7 tesla magnetic resonance elastography

Magnetic resonance elastography (MRE) is a non-invasive method for determining the mechanical response of tissues using applied harmonic deformation and motion-sensitive MRI. MRE studies of the human brain are typically performed at conventional field strengths, with a few attempts at the ultra-high field strength, 7T, reporting increased spatial resolution with partial brain coverage. Achieving high-resolution human brain scans using 7T MRE presents unique challenges of decreased octahedral shear strain-based signal-to-noise ratio (OSS-SNR) and lower shear wave motion sensitivity. In this study, we establish high resolution MRE at 7T with a custom 2D multi-slice single-shot spin-echo echo-planar imaging sequence, using the Gadgetron advanced image reconstruction framework, applying Marchenko-Pastur Principal component analysis denoising, and using nonlinear viscoelastic inversion. These techniques allowed us to calculate the viscoelastic properties of the whole human brain at 1.1 mm isotropic imaging resolution with high OSS-SNR and repeatability. Using phantom models and 7T MRE data of eighteen healthy volunteers, we demonstrate the robustness and accuracy of our method at high-resolution while quantifying the feasible tradeoff between resolution, OSS-SNR, and scan time. Using these post-processing techniques, we significantly increased OSS-SNR at 1.1 mm resolution with whole-brain coverage by approximately 4-fold and generated elastograms with high anatomical detail. Performing high-resolution MRE at 7T on the human brain can provide information on different substructures within brain tissue based on their mechanical properties, which can then be used to diagnose pathologies (e.g. Alzheimer's disease), indicate disease progression, or better investigate neurodegeneration effects or other relevant brain disorders,in vivo.

Three-Dimensional Shear Wave Elastography Using Acoustic Radiation Force and a 2-D Row-Column Addressing (RCA) Array.

Acoustic radiation force (ARF)-based shear wave elastography (SWE) is a clinically available ultrasound imaging mode that noninvasively and quantitatively measures tissue stiffness. Current implementations of ARF-SWE are largely limited to 2-D imaging, which does not provide a robust estimation of heterogeneous tissue mechanical properties. Existing 3-D ARF-SWE solutions that are clinically available are based on wobbler probes, which cannot provide true 3-D shear wave motion detection. Although 3-D ARF-SWE based on 2-D matrix arrays have been previously demonstrated, they do not provide a practical solution because of the need for a high channel-count ultrasound system (e.g., 1024-channel) to provide adequate volume rates and the delicate circuitries (e.g., multiplexers) that are vulnerable to the long-duration "push" pulses. To address these issues, here we propose a new 3-D ARF-SWE method based on the 2-D row-column addressing (RCA) array which has a much lower element count (e.g., 256), provides an ultrafast imaging volume rate (e.g., 2000 Hz), and can withstand the push pulses. In this study, we combined the comb-push shear elastography (CUSE) technique with 2-D RCA for enhanced SWE imaging field-of-view (FOV). In vitro phantom studies demonstrated that the proposed method had robust 3-D SWE performance in both homogenous and inclusion phantoms. An in vivo study on a breast cancer patient showed that the proposed method could reconstruct 3-D elasticity maps of the breast lesion, which was validated using a commercial ultrasound scanner. These results demonstrate strong potential for the proposed method to provide a viable and practical solution for clinical 3-D ARF-SWE.

Coupling quantitative ultrasound using echo envelop statistics and shear wave propagation to provide new image contrast of mimicked liver lesions

In quantitative ultrasound, biomarkers are sensitive to the scatterers’ positioning and number density, and therefore, shear wave (SW) propagation may modify their time-varying behaviors and provide new image contrast. Herein, we considered echo envelop statistics by using one parameter of the homodyned-K distribution (HKD): the diffuse-to-total signal power ratio 1/(κ+1). A liver lesion-mimicking phantom was simulated with an inclusion and a background of elastic moduli of 40 and 4 kPa, and scatterers’ number densities per resolution cell of 15 and 5, respectively. Ten sets of spatially correlated scatterers were created, and a plane SW was simulated to reposition time and spatially varying scatterers according to the SW motion. Ultrasound images were simulated using MUST with added Gaussian noise to attain 20-dB SNR. Dynamic and static analyses were made over 30 frames with and without SW motion, respectively. Validation was based on the contrast of 1/(κ+1) between the inclusion and the background. HKD image contrast increased with SW motion with values from 0.12 ± 0.10 to 0.48 ± 0.16 (no units) (p < 0.001). Results suggest that the variation in scatterers’ organization under SW propagation may provide new information over its static counterpart to enhance the contrast of liver lesions.

Rayleigh-Lamb wave propagation in a prestressed transversely isotropic viscoelastic waveguide: Inverse modeling challenges

The functional role and structure of skeletal muscle results in anisotropy in both material properties and imposed stresses, as well as waveguide effects. Dynamic elastography reconstruction methods for estimating muscle tissue viscoelastic properties that are rooted in assumptions of isotropy and bulk wave motion may produce inaccurate estimates. The superposition of axially-aligned orthotropy (transverse isotropy) in material properties and axially-aligned prestress conditions due to passive stretch or muscle activation makes it difficult to independently discern how much of the apparent anisotropy is due to the muscle material or the imposed stress field. Furthermore, this stress field may result in large strain conditions that require use of higher order terms in the stress-strain relationship. The significance of this confounding condition and strategies for decoupling material and stress-based anisotropy are investigated with a series of numerical finite element and experimental elastography studies using scanning laser Doppler vibrometry and magnetic resonance elastography. Shear and Rayleigh-Lamb wave motion is studied in a polymeric muscle phantom that is in the shape of a rectangular rod and has either isotropic or transversely isotropic material properties under zero stress conditions. [Funding support: NIH AR071162]

Evaluation of Robustness of S-Transform Based Phase Velocity Estimation in Viscoelastic Phantoms and Renal Transplants.

Ultrasound shear wave elastography (SWE) methods are being used to differentiate healthy versus diseased tissue on the basis of their viscoelastic mechanical properties. Tissue viscoelasticity is often studied by analyzing shear wave phase velocity dispersion curves, which is the variation of phase velocity with frequency or wavelength. Recently, a unique approach using a generalized Stockwell transformation (GST-SFK) was proposed for the calculation of dispersion curves in viscoelastic media over expanded frequency band. In this work, the method's robustness was evaluated on data from five custom-made viscoelastic tissue-mimicking phantoms and sixty in vivo renal transplants. For each phantom, 15 shear wave motion data acquisitions were taken, while 10-13 acquisitions were acquired for renal transplants measured in the renal cortex. For each data-set mean and standard deviation (SD) of estimated phase velocity dispersion curves were studied. In addition, the viscoelastic parameters of the Zener model were examined, which were preceded by a convergence analysis. For viscoelastic phantoms scanned with a research ultrasound scanner, and for the in vivo renal transplants scanned with a clinical scanner, the decisive advantage of the GST-SFK method over the standard two-dimensional Fourier transform (2D-FT) method was shown. The GST-SFK method provided dispersion curve estimates with lower SD over a wider frequency band in comparison to the 2D-FT method. These advantages are relevant to the analysis of the mechanical properties of tissues in clinical practice to discriminate healthy from diseased tissue.

Focused shear wave beam propagation through a 3D printed human rib cage

Transient elastography (TE) is used in clinical screenings for liver fibrosis. TE generates shear wave motion in the liver by vibration of a piston at the skin. The small, flat piston in current devices radiates shear waves that diverge from the liver, resulting in low SNR and high failure rates, especially in obese patients. Our group recently introduced focused shear wave beams for TE [Cormack et al., IEEE TBME (2024)], in which vibration of a concave piston generates shear waves that converge towards the focal region, thereby increasing SNR for liver stiffness estimation. However, because piston sizes needed for efficient shear wave focusing are larger than the typical intercostal space, the rib cage that lies between the skin and the liver may cause shear wave aberration and influence stiffness measurement. Here, we present measurements of broadband focused shear wave propagation in tissue-mimicking gel in which are embedded 3D printed human ribs. Both straight elliptical rods and anatomically realistic 3D printed ribs are used as aberrators. Measurements are compared to 3D simulations of shear wave beam propagation [Archer et al., JASA (2023)]. Effects of shear wave aberration by the rib cage depend on shear wavelength, piston-to-rib distance, and piston radius and curvature.

Quantifying uniaxial prestress and waveguide effects on dynamic elastography estimates for a cylindrical rod.

Dynamic elastography attempts to reconstruct quantitative maps of the viscoelastic properties of materials by noninvasively measuring mechanical wave motion in them. The target motion is typically transversely-polarized relative to the wave propagation direction, such as bulk shear wave motion. In addition to neglecting waveguide effects caused by small lengths in one dimension or more, many reconstruction strategies also ignore nonzero, non-isotropic static preloads. Significant anisotropic prestress is inherent to the functional role of some biological materials of interest, which also are small in size relative to shear wavelengths in one or more dimensions. A cylindrically shaped polymer structure with isotropic material properties is statically elongated along its axis while its response to circumferentially-, axially-, and radially-polarized vibratory excitation is measured using optical or magnetic resonance elastography. Computational finite element simulations augment and aid in the interpretation of experimental measurements. We examine the interplay between uniaxial prestress and waveguide effects. A coordinate transformation approach previously used to simplify the reconstruction of un-prestressed transversely isotropic material properties based on elastography measurements is adapted with partial success to estimate material viscoelastic properties and prestress conditions without requiring advanced knowledge of either.

Comb Detection for Measuring Shear Wave Propagation.

Plane wave compounding (PWC) is widely used to measure the propagation of shear waves. Implementing PWC on most commercial ultrasound scanners is challenging because all channel (>128) data must be processed or transferred to the host computing unit in real time. Comb detection transmits multiple focused beams simultaneously and results in a reduced number of receive lines to be processed in parallel. These comb beams are scanned laterally to acquire receive lines at different lateral positions in order to obtain data over a large region of interest (ROI). One of the potential issues with using multiple simultaneously transmitted beams is the issue of crosstalk between the beams. Crosstalk is analyzed through simulated beam patterns, simulated B-mode images, and motion data from shear wave elastography (SWE) experiments. Using a Hamming window on transmit and receive can suppress crosstalk to 1.2% root-mean-square error (RMSE, normalized RMSE to the peak magnitude of the reference signal) for shear wave motion signals. Four comb beams with three laterally scanned locations cover almost the entire field of view (FOV) and achieve the same frame rate as PWC with three angles. Phantom and in vivo studies demonstrate comparable motion data of comb detection to PWC in terms of motion signal quality and measured phase velocity. In addition, comb detection provides motion with lower noise and stronger signals than PWC, which is believed to be due to the advantages of transmitting focused beams rather than plane waves (PWs).

Elastography using torsional wave motion in transverse isotropic material

An experimental configuration for driving and measuring torsional wave motion is designed, fabricated and experimentally tested using magnetic resonance elastography (MRE) in both a ultra-high field preclinical MRI system and in a cryogen-free low-field tabletop MRI system. Cylindrically-shaped isotropic and transversely isotropic phantoms, as well as the excised cat soleus muscle are statically stretched along their main axis to specific prestrain levels while simultaneously conducting MRE measurements with torsional waves. In the case of the excised cat soleus, diffusion tensor imaging is also performed to enable a co-registered mapping of fiber structure. The potential advantage of torsional wave excitation to drive shear wave motion over a large volume with uniform amplitude, and hence SNR, is tested by comparison to excitation of shear waves using other methods.

Crosstalk analysis of comb detection for measuring shear wave propagation

Applications, such as shear wave elastography and Doppler imaging, require high frame rates. Plane wave compounding (PWC) is widely used for these applications. Comb detection was proposed to combine the high frame rate of PWC and the high signal quality of focused beam scanning. Comb detection transmits multiple focused beams simultaneously to increase the frame rate. These comb beams are scanned laterally to cover the whole region-of-interest. The simultaneous multiple transmission of focused beams causes crosstalk between neighboring beams. Transmit crosstalk is related to the pressure field of transmit beams, and receive crosstalk is determined by the beam profile of receive beamforming. In this study, we varied f-number (F/N) and apodization window and measured their effects on crosstalk based on the pressure field and an arterial wall simulation using Field II. For transmit design, transmit apodization with a Hamming window significantly reduced crosstalk compared to a rectangular window and transmit F/N had little impact on the crosstalk. Regarding receive beamforming, a Hamming window led to lower crosstalk than the rectangular window and F/N from 1 to 5 suppressed crosstalk to below 30 dB. The effects of crosstalk in shear wave motion data from phantom experiments were also analyzed. Comb detection showed less motion artifact than PWC in phantom experiments.

Material characterization and simulation for soft gels subjected to impulsive loading

For impact and blast experiments of traumatic brain injury (TBI), soft gel materials are used as surrogates to imitate the mechanical responses of brain tissue. To properly model a viscoelastic gel brain in a surrogate head using a finite element (FE) model, material parameters such as the shear moduli and relaxation time at high strain rates are required. However, such information is scarce in the literature and its applicability for a range of dynamic conditions is unclear. We used an integrated experiment and simulation approach to efficiently determine mechanical properties of soft gels at finite strains, as well as over a wide range of strain rates. A novel impact experiment using a gel block was developed to capture the high strain rate behavior by maximizing the inherent shear wave motion at different impact conditions. A corresponding computational model was used to simulate the gel dynamics of the impact. Parametric simulations utilizing optimization and correlation analyses were used to calibrate multiple material parameters in the nonlinear viscoelastic model to the experimental data. The optimal parameters for gels, including Sylgards 184, 3–6636, and 527, were found. We ascertained the initial shear stiffening effect in gels at high strain rate loadings experimentally and incorporated this effect in the simulation. We have verified the integrated approach by comparing the material properties of the gels with analytical results based on shear wave propagation. This study provides a new approach to calibrate the material behavior of soft gels under high strain rate loading conditions.

Homogenization of the wave equation with non-uniformly oscillating coefficients

The focus of our work is a dispersive, second-order effective model describing the low-frequency wave motion in heterogeneous (e.g., functionally graded) media endowed with periodic microstructure. For this class of quasi-periodic medium variations, we pursue homogenization of the scalar wave equation in , , within the framework of multiple scales expansion. When either or , this model problem bears direct relevance to the description of (anti-plane) shear waves in elastic solids. By adopting the lengthscale of microscopic medium fluctuations as the perturbation parameter, we synthesize the germane low-frequency behavior via a fourth-order differential equation (with smoothly varying coefficients) governing the mean wave motion in the medium, where the effect of microscopic heterogeneities is upscaled by way of the so-called cell functions. In an effort to demonstrate the relevance of our analysis toward solving boundary value problems (deemed to be the ultimate goal of most homogenization studies), we also develop effective boundary conditions, up to the second order of asymptotic approximation, applicable to one-dimensional (1D) shear wave motion in a macroscopically heterogeneous solid with periodic microstructure. We illustrate the analysis numerically in one dimension by considering (i) low-frequency wave dispersion, (ii) mean-field homogenized description of the shear waves propagating in a finite domain, and (iii) full-field homogenized description thereof. In contrast to (i) where the overall wave dispersion appears to be fairly well described by the leading-order model, the results in (ii) and (iii) demonstrate the critical role that higher-order corrections may have in approximating the actual waveforms in quasi-periodic media.

Three-Dimensional Shear Wave Elastography Using a 2D Row Column Addressing (RCA) Array.

Objective. To develop a 3D shear wave elastography (SWE) technique using a 2D row column addressing (RCA) array, with either external vibration or acoustic radiation force (ARF) as the shear wave source. Impact Statement. The proposed method paves the way for clinical translation of 3D SWE based on the 2D RCA, providing a low-cost and high volume rate solution that is compatible with existing clinical systems. Introduction. SWE is an established ultrasound imaging modality that provides a direct and quantitative assessment of tissue stiffness, which is significant for a wide range of clinical applications including cancer and liver fibrosis. SWE requires high frame rate imaging for robust shear wave tracking. Due to the technical challenges associated with high volume rate imaging in 3D, current SWE techniques are typically confined to 2D. Advancing SWE from 2D to 3D is significant because of the heterogeneous nature of tissue, which demands 3D imaging for accurate and comprehensive evaluation. Methods. A 3D SWE method using a RCA array was developed with a volume rate up to 2000 Hz. The performance of the proposed method was systematically evaluated on tissue-mimicking elasticity phantoms and in an in vivo case study. Results. 3D shear wave motion induced by either external vibration or ARF was successfully detected with the proposed method. Robust 3D shear wave speed maps were reconstructed for phantoms and in vivo. Conclusion. The high volume rate 3D imaging provided by the 2D RCA array provides a robust and practical solution for 3D SWE with a clear pathway for future clinical translation.

Damping formulations for finite difference linear dynamic analyses: Performance and practical recommendations

An in-depth study is presented on some numerical formulations used to reproduce damping of motion in solid continua. The focus is addressed in particular to those formulations that in the finite difference geotechnical code FLAC can be exploited to perform linear analyses. The numerical investigations are performed through a simple 1D model of a homogeneous soil deposit interested by a shear wave motion coming from the bedrock, and a gravity dam model subjected to seismic excitations in presence of dam-water-foundation interaction. The Rayleigh damping method is tested in particular to find an appropriate choice of the control frequency in order to minimize the problems of high-frequency overdamping. The analyses also investigate two other formulations (the local damping and its variant, the combined damping) that present the major advantages in providing a frequency independent action and in needing not significant computational resources. Limits and principal drawbacks of their possible use are highlighted.

Anti-plane shear wave motion in a composite layered structure with slit

Anti-plane shear wave propagation in a functionally graded piezoelectric substrate and a piezoelectric substrate guided by a slit is analyzed theoretically. The equations of motion are taken from the linear theory of piezoelectricity for polarized ceramics. The dispersion relation is obtained by applying appropriate boundary conditions on the derived solutions. PZT-5 and ZnO are taken as numerical examples to manifest the effect of different material parameters (gradient factor, slit width etc.) graphically on the phase and group velocity of the considered wave. LiI O 3 and BaTi O 3 materials are taken separately to compare the results. Also, graphs are plotted for different modes of the wave as well as for different frequency ranges and for the case when slit is absent. As a special case, it is remarked that the present study has resemblance with one of the existing works. The highlights of the study are to characterize the anti-plane wave propagation within a functionally graded piezoelectric substrate/gap/piezoelectric substrate coupled structure, to show the effect of different parameters, modes, frequency, etc. on the velocity profile of the considered wave and to show the effect of gap layer. The outcomes of this study may be used to enhance the sensitivity of acoustic wave devices and sensors.

A fast and quantitative measurement of shear wave speed and attenuation using (cross) modal assurance criterion for acoustic radiation force elasticity imaging

The measurement of shear wave speed (SWS) or attenuation in soft materials, especially in biological tissues, plays a crucial role in acoustic radiation force elasticity imaging which provides a quantitative measurement of elastic properties. In this study, with the first attempt to apply modal analysis theory to acoustic radiation force elasticity imaging, we propose a fast tracking algorithm based on the modal assurance criterion (MAC) including 2D SWS-MAC for SWS measurement and the signed cross modal assurance criterion (CrossMAC) for shear wave attenuation measurement. These approaches achieved high speedup with a framework of high-level parallelism. The FEM simulation results proved the accuracy of these two approaches. Phantom experiments showed that the 2D SWS-MAC approach had good consistency in SWS estimation compared to existing approaches and could accurately preserve the boundary contour of the inclusion with narrowed assurance intervals (quarter-width assurance intervals). An excised porcine liver experiment verified the effectiveness of shear wave attenuation estimation with Signed CrossMAC. These results demonstrate that the MAC based approaches improved both the shear wave motion measurement precision and the speedup rate of the measurement compared to conventional time domain cross-correlation or time-to-peak methods.

Time-Aligned Plane Wave Compounding Methods for High-Frame-Rate Shear Wave Elastography: Experimental Validation and Performance Assessment on Tissue Phantoms

Shear wave elastography (SWE) is an ultrasonic technique able to quantitatively assess the mechanical properties of tissues by combining acoustic radiation force and ultrafast imaging. While utilizing coherent plane wave compounding enhances echo and shear wave motion signal-to-noise ratio (SNR), it also reduces the effective pulse repetition frequency (PRFe), affecting the accuracy of the measurements of motion and, consequently, of material properties. It is important to maintain both high-motion SNR and PRFe, particularly for the characterization of (material and/or geometrical) dispersive tissues such as arteries. This work proposes a method for SWE measurements with high SNR, while maintaining a high PRFe, using conventional clinical ultrasound scanners. A time alignment process is applied after acquiring data from plane wave transmissions at different angles. The time alignment uses interpolation to obtain data points at higher frame rates, and the time-aligned data are compounded to increase the SNR. The method is used for SWE in tissue-mimicking phantoms of different stiffness and is compared with traditional plane wave compounding. Increases of 58% and 36% in spatial and temporal bandwidth compared with conventional plane wave compounding, respectively, can be achieved for SWE measurements of representative arterial stiffness values. Improvements in phase velocity accuracy and bandwidth in an arterial phantom are also described, to emphasize the beneficial advantage in dispersive cases.

Viscoelastic parameter estimation using simulated shear wave motion and convolutional neural networks

Ultrasound shear wave elastography (SWE) techniques have been very useful for the analysis of tissue rheological properties, but there are still obstacles for robust evaluation of viscoelastic tissue properties. In this proof-of-concept study, we investigate whether convolutional neural networks (CNN) are capable of retrieving the elasticity and viscosity parameters from simulated shear wave motion images. Staggered-grid finite difference simulations based on a Kelvin-Voigt rheological model were used to generate data for this study. The wave motion datasets were created using Kelvin-Voigt shear elasticity values ranging from 1 to 25 kPa, shear viscosities ranging from 0 to 10 Pa⋅s, and two different push profiles using f-numbers of 1 and 2. The CNN architectures, optimized using mean squared error loss, were then trained to retrieve a specific viscoelastic parameter. Both elasticity and viscosity values were successfully retrieved, with regression R2 values above 0.99 when correlating the estimated mechanical properties versus the true mechanical properties. The CNN performance was also compared to estimation of shear elasticity and viscosity from fitting dispersion curves estimated from two-dimensional Fourier transform analysis. The results demonstrated that the CNN models were robust to noise, vertical position and partially to f-number. The architecture was proven to be robust to multiple push profiles if trained properly. The CNN results showed higher accuracy over the full viscoelastic parameter range compared to the Fourier-based analysis. The overall results showed the CNNs’ potential to be an alternative to complex mathematical analyses such as Fourier analysis and dispersion curve estimation used currently for shear wave viscoelastic parameter estimation.

Closed-chest measurement of diastolic and systolic shear wave speed to assess myocardial stiffness

Abstract Funding Acknowledgements Type of funding sources: Private grant(s) and/or Sponsorship. Main funding source(s): Research Foundation Flanders (FWO): grant 1211620N TTW – Dutch Heart Foundation partnership program "Earlier recognition of cardiovascular diseases": project number 14740 Background Echocardiographic shear wave elastography (SWE) encompasses all ultrasound techniques tracking shear wave (SW) motion in the cardiac wall, of which the propagation speed is linked to the intrinsic mechanical properties. SWs can be induced naturally, for example by valve closure, or externally by using an acoustic radiation force (ARF). Although the latter is technically more demanding, it enables instantaneous stiffness assessment throughout the entire cardiac cycle (fig. a). However, it is unknown how factors such as cardiac loading and contractility, next to intrinsic mechanical properties, affect ARF-based SW speeds. Purpose We performed transthoracic SWE measurements in pigs to study the effects of hemodynamic alterations, inotropic state and myocardial infarction (MI) on diastolic and systolic SW speeds. Methods Different cardiac conditions were considered in three pigs: (i) baseline (BL), (ii) preload decrease (PD), (iii) afterload increase (AI), (iv) preload increase (PI), (v) administration of dobutamine (DOB), (vi) BL2, (vii) MI through 60-100 min. occlusion of the LAD and (viii) 40 min. reperfusion (REP). For each condition, transthoracic high frame rate ARF-based SWE acquisitions were taken in a parasternal long-axis view with a research ultrasound system. SWs were induced in the septum at 34 Hz during 1.5 s to track SW speeds throughout the cardiac cycle (fig. a&b). Systolic and diastolic SW speeds were determined from the 10% highest and lowest median values per condition, respectively. Left ventricular pressure-volume (PV) loops were recorded to estimate end-diastolic pressure (EDP), end-systolic pressure (ESP) and passive chamber stiffness (dPdV). dPdV was determined as the slope of the tangent to the fitted end-diastolic PV relationship at mean ED volume. Linear regressions and Pearson’s correlation coefficients were computed. Results Diastolic SW speed was correlated to EDP for conditions with changes in loading, and to dPdV for conditions with changes in chamber stiffness (fig. c). Both relationships were significant, with a moderate positive correlation for EDP (R = 0.48, p = 0.02) and a strong positive correlation for dPdV (R = 0.76, p < 0.01). Furthermore, the observed change in diastolic SW speed was smaller when altering EDP compared to dPdV (0.4 m/s vs. 1.0 m/s). For systolic SW speed, very strong positive correlations were found with ESP (R = 0.91, p < 0.01), and with dPdV (R = 0.81, p < 0.01) in fig. d. Conclusion This study shows that both diastolic and systolic SW speed are related to passive chamber stiffness. Moreover, loading also influenced systolic SW speed and, to a lesser extent, diastolic SW speed, presumably because of material nonlinearity. Systolic SW speed is linked to contractility as well. Thus, while SWs after valve closure occur at a certain moment in the cardiac cycle, the timing of ARF-based SWs can be chosen such to assess specific aspects of the cardiac (structural and functional) status. Abstract Figure.

Phase Velocity Estimation With Expanded Bandwidth in Viscoelastic Phantoms and Tissues.

Ultrasound shear wave elastography (SWE) is a technique used to measure mechanical properties to evaluate healthy and pathological soft tissues. SWE typically employs an acoustic radiation force (ARF) to generate laterally propagating shear waves that are tracked in the spatiotemporal domains, and algorithms are used to estimate the wave velocity. The tissue viscoelasticity is often examined through analyzing the shear wave phase velocity dispersion curves, which is the variation of phase velocity with frequency or wavelength. A number of available methods to estimate dispersion exist, which can differ in resolution and variance. Moreover, most of these techniques reconstruct dispersion curves for a limited frequency band. In this work, we propose a novel method used for dispersion curve calculation. Our unique approach uses a generalized Stockwell transformation combined with a slant frequency-wavenumber analysis (GST-SFK). We tested the GST-SFK method on numerical phantom data generated using a finite-difference-based method in tissue-mimicking viscoelastic media. In addition, we evaluated the method on numerical shear wave motion data with different amounts of white Gaussian noise added. Additionally, we performed tests on data from custom-made tissue-mimicking viscoelastic phantom experiments, ex vivo porcine liver measurements, and in vivo liver tissue experiments. We compared results from our method with two other techniques used for estimating shear wave phase velocity: the two-dimensional Fourier transform (2D-FT) and the eigenvector (EV) method. Tests carried out revealed that the GST-SFK method provides dispersion curve estimates with lower errors over a wider frequency band in comparison to the 2D-FT and EV methods. In addition, the GST-SFK provides expanded bandwidth by a factor of two or more to be used for phase velocity estimation, which is meaningful for a tissue dispersion analysis in vivo.