High Frame Rate Ultrasound Particle Image Velocimetry for Estimating High Velocity Flow Patterns in the Left Ventricle

Ex-Gaussian vector metric and similarity index applied to reaction time analysis

In psychology and cognitive neuroscience, reaction time (RT) series and their ex-Gaussian distributions are commonly used as scalar quantities to explore the time course of attentional processes. However, we propose that such attentional processes can also be analyzed using an "ex-Gaussian vector", defined by successive triads of ex-Gaussian sigma, tau, and mu parameters from RT series. This geometrical object may help characterize interindividual differences between congruent and incongruent stimuli in the attentional Stroop task within a group of participants. To test these hypotheses, we calculated the similarity index of these geometrical objects in young adults without detectable neurological disorders. Our findings show that during two weeks of continuous Stroop task application, each participant displayed distinct ex-Gaussian RT vectors in a Cartesian 3D plot. Furthermore, our study found that the similarity index between ex-Gaussian RT vectors was significantly higher for incongruent stimuli than for congruent stimuli.

Validation of Left Ventricular High Frame Rate Echo-Particle Image Velocimetry against 4D Flow MRI in Patients

OBJECTIVE

Accurately measuring intracardiac flow patterns could provide insights into cardiac disease pathophysiology, potentially enhancing diagnostic and prognostic capabilities. This study aims to validate Echo-Particle Image Velocimetry (echoPIV) for in vivo left ventricular intracardiac flow imaging against 4D flow MRI.

METHODS

We acquired high frame rate contrast-enhanced ultrasound images from three standard apical views of 26 patients who required cardiac MRI. 4D flow MRI was obtained for each patient. Only echo image planes with sufficient quality and alignment with MRI were included for validation. Regional velocity, kinetic energy (KE) and viscous energy loss (EL˙) were compared between modalities using normalized mean absolute error (NMAE), cosine similarity and Bland-Altman analysis.

RESULTS

Among 24 included apical view acquisitions, we observed good correspondence between echoPIV and MRI regarding spatial flow patterns and vortex traces. The velocity profile at base-level (mitral valve) cross-section had cosine similarity of 0.92 ± 0.06 and NMAE of (14 ± 5)%. Peak spatial mean velocity differed by (3 ± 6) cm/s in systole and (6 ± 10) cm/s in diastole. The KE and rate of EL˙ also revealed a high level of cosine similarity (0.89 ± 0.09 and 0.91 ± 0.06) with NMAE of (23 ± 7)% and (52 ± 16)%.

CONCLUSION

Given good B-mode image quality, echoPIV provides a reliable estimation of left ventricular flow, exhibiting spatial-temporal velocity distributions comparable to 4D flow MRI. Both modalities present respective strengths and limitations: echoPIV captured inter-beat variability and had higher temporal resolution, while MRI was more robust to patient BMI and anatomy.

Ultrasound Speckle Decorrelation Analysis-Based Velocimetry for 3D-Velocity-Components Measurement Using a 1D Transducer Array

Ultrasound velocimetry has been widely used for blood flow imaging. However, the flow measurements are constrained to resolve the in-plane 2D flow components when using a 1D transducer array. In this work, an ultrasound speckle decorrelation analysis-based velocimetry (3C-vUS) is proposed for 3D velocity components measurement using a 1D transducer array. The 3C-vUS theory is first derived and validated with numerical simulations and phantom experiments. The in vivo testing results show that 3C-vUS can accurately measure the blood flow 3D-velocity-components of the human carotid artery at arbitrary probe-to-vessel angles throughout the cardiac cycle. With such capability, the 3C-vUS will alleviate the requirement of operators and promote disease screening for blood flow-related disorders.

Development of an in vitro setup for flow studies in a stented carotid artery bifurcation

To investigate flow conditions in a double-layered carotid artery stent, a bench-top in vitro flow setup including a bifurcation phantom was designed and fabricated. The geometry of the tissue-mimicking phantom was based on healthy individuals. Two identical phantoms were created using 3D-printing techniques and molding with PVA-gel. In one of them, a clinically available CGuard double-layer stent was inserted. Measurements were performed using both continuous and pulsatile flow conditions. Blood flow studies were performed using echoPIV: a novel ultrasound-based technique combined with particle image velocimetry. A maximum deviation of 3% was visible between desired and measured flow patterns. The echoPIV measurements showed promising results on visualization and quantification of blood flow in and downstream the stent. Further research could demonstrate the effects of a double-layered stent on blood flow patterns in a carotid bifurcation in detail.

The Effect of Spatial Velocity Gradients on Block-Matching Accuracy for Ultrasound Velocimetry

OBJECTIVE

Block matching serves as the foundation for ultrasound velocimetry techniques such as blood speckle tracking and echo-particle image velocimetry. Any spatial velocity gradients (SVGs) inside a block-matching pair will result in tracking error, due to both the finite block size and the ultrasound point-spread-function. We assess, using an in silico sinusoidal flow phantom, the effect of SVG magnitude and beam-to-flow angle on block-matching bias and precision. Secondarily we assess the effect that SVGs have on velocimetry bias when using angled plane-wave compounding.

METHODS

The magnitude and angle of SVGs were varied by adjusting the wavelength and direction of a sinusoidal flow profile. Scatterers displaced by this flow profile were used for simulating ultrasound radio frequency data at discrete time points. After beamforming, the 2-D flow field was estimated using block matching. Two imaging sequences were tested, a single plane-wave and a three-angled plane-wave.

RESULTS

Smaller sinusoidal flow wavelengths resulted in increased bias and reduced precision, revealing an inverse relationship between sinusoidal flow wavelength and tracking error, with median errors ranging from 69%-90% for the smallest flow wavelengths (highest SVGs) down to 3%-5% for the largest (lowest SVGs). The SVG angle was also important, in which lateral SVGs (with axially oriented flows) resulted in significant speckle decorrelation and high tracking errors in regions with high SVGs. Conversely, axial SVGs (laterally oriented flow) experienced higher bias in the peak velocity regions of the flow profile. Coherent compounding resulted in higher velocity errors than using a single transmission for lateral SVGs but not for axial SVGs.

CONCLUSION

The highest SVGs that could be measured with ≤10% error was when the sinusoidal flow wavelength was less than 20 times the ultrasound pulse wavelength. The clinical significance is that the high SVGs present in high kinetic energy flows, such as severe carotid stenosis and aortic regurgitation, will limit the ability to accurately quantify the velocities in these flow structures.

Blood Flow Quantification with High-Frame-Rate, Contrast-Enhanced Ultrasound Velocimetry in Stented Aortoiliac Arteries: In Vivo Feasibility

Local flow patterns influence stent patency, while blood flow quantification in stents is challenging. The aim of this study was to investigate the feasibility of 2-D blood flow quantification using high-frame-rate, contrast-enhanced ultrasound (HFR-CEUS) and particle image velocimetry (PIV), or echoPIV, in patients with aortoiliac stents. HFR-CEUS measurements were performed at 129 locations in 62 patients. Two-dimensional blood flow velocity fields were obtained using echoPIV. Visual inspection was performed by five observers to evaluate feasibility. The contrast-to-background ratio and average vector correlation were calculated and compared between stented and native vessel segments. Flow quantification with echoPIV was feasible in 128 of 129 locations (99%), with optimal quantification in 40 of 129 locations (31%). Partial quantification was achieved in 88 of 129 locations (68%), where one or multiple limiting issues occurred (not related to the stent) including loss of correlation during systole (57/129), short vessel segments (20/129), loss of contrast during diastole (20/129) and shadow regions (20/129). The contrast-to-background ratio and vector correlation were lower downstream in the imaged blood vessel, independent of the location of the stent. In conclusion, echoPIV was feasible in stents placed in the aortoiliac region, and the stents did not adversely affect flow tracking.

Left ventricular high frame rate echo-particle image velocimetry: clinical application and comparison with conventional imaging

Background

Echo-Particle Image Velocimetry (echoPIV) tracks speckle patterns from ultrasound contrast agent(UCA), being less angle-sensitive than colour Doppler. High frame rate (HFR) echoPIV enables tracking of high velocity flow in the left ventricle (LV). We aimed to demonstrate the potential clinical use of HFR echoPIV and investigate the feasibility and accuracy in patients.

Methods

Nineteen patients admitted for heart failure were included. HFR contrast images were acquired from an apical long axis view (ALAX), using a fully-programmable ultrasound system. A clinical UCA was continuously infused with a dedicated pump. Additionally, echocardiographic images were obtained using a clinical system, including LV contrast-enhanced images and pulsed-wave (PW) Doppler of the LV inflow and outflow in ALAX. 11 patients underwent CMR and 4 cardiac CT as clinically indicated. These CMR and CT images were used as reference. In 10 patients with good echoPIV tracking and reference imaging, the intracavitary flow was compared between echoPIV, conventional and UCA echocardiography.

Results

EchoPIV tracking quality was good in 12/19 (63%), moderate in 2/19 (10%) and poor in 5/19 (26%) subjects. EchoPIV could determine inflow velocity in 17/19 (89%), and outflow in 14/19 (74%) patients. The correlation of echoPIV and PW Doppler was good for the inflow (R² = 0.77 to PW peak; R² = 0.80 PW mean velocity) and moderate for the outflow (R² = 0.54 to PW peak; R² = 0.44 to PW mean velocity), with a tendency for echoPIV to underestimate PW velocities. In selected patients, echoPIV was able in a single acquisition to demonstrate flow patterns which required multiple interrogations with classical echocardiography. Those flow patterns could also be linked to anatomical abnormalities as seen in CMR or CT.

Conclusion

HFR echoPIV tracks multidirectional and complex flow patterns which are unapparent with conventional echocardiography, while having comparable feasibility. EchoPIV tends to underestimate flow velocities as compared to PW Doppler. It has the potential to provide in one acquisition all the functional information obtained by conventional imaging, overcoming the angle dependency of Doppler and low frame rate of classical contrast imaging.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12947-022-00283-4.

Ultrasonic biomechanics method for vortex and wall motion of left ventricle: a phantom and in vivo study

Background

The non-invasive quantitative evaluation of left ventricle (LV) function plays a critical role in clinical cardiology. This study proposes a novel ultrasonic biomechanics method by integrating both LV vortex and wall motion to fully assess and understand the LV structure and function. The purpose of this study was to validate the ultrasonic biomechanics method as a quantifiable approach to evaluate LV function.

Methods

Firstly, B-mode ultrasound images were acquired and processed, which were utilized to implement parameters for quantifying the LV vortex and wall motion respectively. Next, the parameters were compared in polyvinyl alcohol cryogen (PVA) phantoms with different degree of stiffness corresponding to different freezing and thawing cycles in vitro. Finally, the parameters were computed in vivo during one cardiac cycle to assess the LV function in normal and abnormal subjects in vivo.

Results

In vitro study, the velocity field of PVA phantom differed with stiffness (varied elasticity modulus). The peak of strain for wall motion decreases with the increase of elasticity modulus, and periodically changed values. Statistical analysis for parameters of vortex dynamics (energy dissipation index, DI; kinetic energy fluctuations, KEF; relative strength, RS; and vorticity, W) based on different elasticity (E) of phantom depicted the good viability of this algorithm. In vivo study, the results confirmed that subjects with LV dysfunction had lower vorticity and strain (S) compared to the normal group.

Conclusion

Ultrasonic biomechanics method can obtain the vortex and wall motion of left ventricle. The method may have potential clinical value in evaluation of LV dysfunction.

Optimization of Microbubble Concentration and Acoustic Pressure for Left Ventricular High-Frame-Rate EchoPIV in Patients

High-frame-rate (HFR) echo-particle image velocimetry (echoPIV) is a promising tool for measuring intracardiac blood flow dynamics. In this study, we investigate the optimal ultrasound contrast agent (UCA: SonoVue) infusion rate and acoustic output to use for HFR echoPIV (PRF = 4900 Hz) in the left ventricle (LV) of patients. Three infusion rates (0.3, 0.6, and 1.2 ml/min) and five acoustic output amplitudes (by varying transmit voltage: 5, 10, 15, 20, and 30 V-corresponding to mechanical indices of 0.01, 0.02, 0.03, 0.04, and 0.06 at 60-mm depth) were tested in 20 patients admitted for symptoms of heart failure. We assess the accuracy of HFR echoPIV against pulsed-wave Doppler acquisitions obtained for mitral inflow and aortic outflow. In terms of image quality, the 1.2-ml/min infusion rate provided the highest contrast-to-background ratio (CBR) (3-dB improvement over 0.3 ml/min). The highest acoustic output tested resulted in the lowest CBR. Increased acoustic output also resulted in increased microbubble disruption. For the echoPIV results, the 1.2-ml/min infusion rate provided the best vector quality and accuracy; mid-range acoustic outputs (corresponding to 15-20-V transmit voltages) provided the best agreement with the pulsed-wave Doppler. Overall, the highest infusion rate (1.2 ml/min) and mid-range acoustic output amplitudes provided the best image quality and echoPIV results.

In vitro performance of echoPIV for assessment of laminar flow profiles in a carotid artery stent

Purpose: Detailed blood flow studies may contribute to improvements in carotid artery stenting. High-frame-rate contrast-enhanced ultrasound followed by particle image velocimetry (PIV), also called echoPIV, is a technique to study blood flow patterns in detail. The performance of echoPIV in presence of a stent has not yet been studied extensively. We compared the performance of echoPIV in stented and nonstented regions in an in vitro flow setup.

Approach: A carotid artery stent was deployed in a vessel-mimicking phantom. High-frame-rate contrast-enhanced ultrasound images were acquired with various settings. Signal intensities of the contrast agent, velocity values, and flow profiles were calculated.

Results: The results showed decreased signal intensities and correlation coefficients inside the stent, however, PIV analysis in the stent still resulted in plausible flow vectors.

Conclusions: Velocity values and laminar flow profiles can be measured in vitro in stented arteries using echoPIV.

Numerical simulations of flow patterns in the human left ventricle model with a novel dynamic mesh morphing approach based on radial basis function

We present a new numerical simulation framework for prediction of flow patterns in the human left ventricle model. In this study, a radial basis function (RBF) mesh morphing method is developed and applied within the finite-volume computational fluid dynamics (CFD) approach. The numerical simulations are designed to closely mimic details of recent tomographic particle image velocimetry (TomoPIV) experiments. The numerically simulated dynamic motions of the left ventricle and tri-leaflet biological mitral valve are emulated through the RBF morphing method. The arbitrary Lagrangian-Eulerian (ALE) based CFD is performed with the RBF-defined deforming wall boundaries. The results obtained show a good agreement with experiments, confirming the reliability and accuracy of the developed simulation framework.

Contrast-Enhanced High-Frame-Rate Ultrasound Imaging of Flow Patterns in Cardiac Chambers and Deep Vessels

Cardiac function and vascular function are closely related to the flow of blood within. The flow velocities in these larger cavities easily reach 1 m/s, and generally complex spatiotemporal flow patterns are involved, especially in a non-physiologic state. Visualization of such flow patterns using ultrasound can be greatly enhanced by administration of contrast agents. Tracking the high-velocity complex flows is challenging with current clinical echographic tools, mostly because of limitations in signal-to-noise ratio; estimation of lateral velocities; and/or frame rate of the contrast-enhanced imaging mode. This review addresses the state of the art in 2-D high-frame-rate contrast-enhanced echography of ventricular and deep-vessel flow, from both technological and clinical perspectives. It concludes that current advanced ultrasound equipment is technologically ready for use in human contrast-enhanced studies, thus potentially leading to identification of the most clinically relevant flow parameters for quantifying cardiac and vascular function.

A Systematic Review for the Design of In Vitro Flow Studies of the Carotid Artery Bifurcation

PURPOSE

In vitro blood flow studies in carotid artery bifurcation models may contribute to understanding the influence of hemodynamics on carotid artery disease. However, the design of in vitro blood flow studies involves many steps and selection of imaging techniques, model materials, model design, and flow visualization parameters. Therefore, an overview of the possibilities and guidance for the design process is beneficial for researchers with less experience in flow studies.

METHODS

A systematic search to in vitro flow studies in carotid artery bifurcation models aiming at quantification and detailed flow visualization of blood flow dynamics results in inclusion of 42 articles.

RESULTS

Four categories of imaging techniques are distinguished: MRI, optical particle image velocimetry (PIV), ultrasound and miscellaneous techniques. Parameters for flow visualization are categorized into velocity, flow, shear-related, turbulent/disordered flow and other parameters. Model materials and design characteristics vary between study type.

CONCLUSIONS

A simplified three-step design process is proposed for better fitting and adequate match with the pertinent research question at hand and as guidance for less experienced flow study researchers. The three consecutive selection steps are: flow parameters, image modality, and model materials and designs. Model materials depend on the chosen imaging technique, whereas choice of flow parameters is independent from imaging technique and is therefore only determined by the goal of the study.

4-D Echo-Particle Image Velocimetry in a Left Ventricular Phantom

Left ventricular (LV) blood flow is an inherently complex time-varying 3-D phenomenon, where 2-D quantification often ignores the effect of out-of-plane motion. In this study, we describe high frame rate 4-D echocardiographic particle image velocimetry (echo-PIV) using a prototype matrix transesophageal transducer and a dynamic LV phantom for testing the accuracy of echo-PIV in the presence of complex flow patterns. Optical time-resolved tomographic PIV (tomo-PIV) was used as a reference standard for comparison. Echo-PIV and tomo-PIV agreed on the general profile of the LV flow patterns, but echo-PIV smoothed out the smaller flow structures. Echo-PIV also underestimated the flow rates at greater imaging depths, where the PIV kernel size and transducer point spread function were large relative to the velocity gradients. We demonstrate that 4-D echo-PIV could be performed in just four heart cycles, which would require only a short breath-hold, providing promising results. However, methods for resolving high velocity gradients in regions of poor spatial resolution are required before clinical translation.

High-Frame-Rate Contrast-Enhanced Echocardiography Using Diverging Waves: 2-D Motion Estimation and Compensation

Combining diverging ultrasound waves and microbubbles could improve contrast-enhanced echocardiography (CEE), by providing enhanced temporal resolution for cardiac function assessment over a large imaging field of view. However, current image formation techniques using coherent summation of echoes from multiple steered diverging waves (DWs) are susceptible to tissue and microbubble motion artifacts, resulting in poor image quality. In this study, we used correlation-based 2-D motion estimation to perform motion compensation for CEE using DWs. The accuracy of this motion estimation method was evaluated with Field II simulations. The root-mean-square velocity errors were 5.9% ± 0.2% and 19.5% ± 0.4% in the axial and lateral directions, when normalized to the maximum value of 62.8 cm/s which is comparable to the highest speed of blood flow in the left ventricle (LV). The effects of this method on image contrast ratio (CR) and contrast-to-noise ratio (CNR) were tested in vitro using a tissue mimicking rotating disk with a diameter of 10 cm. Compared against the control without motion compensation, a mean increase of 12 dB in CR and 7 dB in CNR were demonstrated when using this motion compensation method. The motion correction algorithm was tested in vivo on a CEE data set acquired with the Ultrasound Array Research Platform II performing coherent DW imaging. Improvement of the B-mode and contrast-mode image quality with cardiac motion and blood flow-induced microbubble motion was achieved. The results of motion estimation were further processed to interpret blood flow in the LV. This allowed for a triplex cardiac imaging technique, consisting of B mode, contrast mode, and 2-D vector flow imaging with a high frame rate of 250 Hz.