Common Carotid Artery Flow Measured by 3-D Ultrasonic Vector Flow Imaging and Validated with Magnetic Resonance Imaging

Intra-Cardiac Flow from Geometry Prescribed Computational Fluid Dynamics: Comparison with Ultrasound Vector Flow Imaging

PURPOSE

This paper investigates the accuracy of blood flow velocities simulated from a geometry prescribed computational fluid dynamics (CFD) pipeline by applying it to a dynamic heart phantom. The CFD flow patterns are compared to a direct flow measurement by ultrasound vector flow imaging (VFI). The hypothesis is that the simulated velocity magnitudes are within one standard deviation of the measured velocities.

METHODS

The CFD pipeline uses computed tomography angiography (CTA) images with 20 volumes per cardiac cycle as geometry input. Fluid domain movement is prescribed from volumetric image registration using the CTA image data. Inlet and outlet conditions are defined by the experimental setup. VFI is systematically measured in parallel planes, and compared to the corresponding planes in the simulated time dependent three dimensional fluid velocity field.

RESULTS

The measured VFI and simulated CFD have similar flow patterns when compared qualitatively. A quantitative comparison of the velocity magnitude is also performed at specific regions of interest. These are evaluated at 11 non-overlapping time bins and compared by linear regression giving R2 = 0.809, SD = 0.060 m/s, intercept = - 0.039 m/s, and slope = 1.09. Excluding an outlier at the inlet, the correspondence between CFD and VFI improves to: R2 = 0.823, SD = 0.048 m/s, intercept = -0.030 m/s, and slope = 1.01.

CONCLUSION

The direct comparison of flow patterns shows that the proposed CFD pipeline provide realistic flow patterns in a well-controlled experimental setup. The demanded accuracy is obtained close to the inlet and outlet, but not in locations far from these.

Precise Estimation of Intravascular Pressure Gradients

This study presents a method for noninvasive pressure gradient estimation, which allows the detection of small pressure differences with higher precision compared to invasive catheters. It combines a new method for estimating the temporal acceleration of the flowing blood with the Navier-Stokes equation. The acceleration estimation is based on a double cross-correlation approach, which is hypothesized to minimize the influence of noise. Data are acquired using a 256-element, 6.5-MHz GE L3-12-D linear array transducer connected to a Verasonics research scanner. A synthetic aperture (SA) interleaved sequence with 2 ×12 virtual sources evenly distributed over the aperture and permuted in emission order is used in combination with recursive imaging. This enables a temporal resolution between correlation frames equal to the pulse repetition time at a frame rate of half the pulse repetition frequency. The accuracy of the method is evaluated against a computational fluid dynamic simulation. Here, the estimated total pressure difference complies with the CFD reference pressure difference, which yields an R -square of 0.985 and an RMSE of 3.03 Pa. The precision of the method is tested on experimental data, measured on a carotid phantom of the common carotid artery. The volume profile used during measurement was set to mimic flow in the carotid artery with a peak flow rate of 12.9 mL/s. The experimental setup showed that the measured pressure difference changes from -59.4 to 31 Pa throughout a single pulse cycle. This was estimated with a precision of 5.44% (3.22 Pa) across ten pulse cycles. The method was also compared to invasive catheter measurements in a phantom with a 60% cross-sectional area reduction. The ultrasound method detected a maximum pressure difference of 72.3 Pa with a precision of 3.3% (2.22 Pa). The catheters measured a maximum pressure difference of 105 Pa with a precision of 11.2% (11.4 Pa). This was measured over the same constriction and with a peak flow rate of 12.9 mL/s. The double cross-correlation approach revealed no improvement compared to a normal differential operator. The method's strength, thus, lies primarily in the ultrasound sequence, which allows precise and accurate velocity estimations, at which acceleration and pressure differences can be acquired.

Anatomic and Functional Imaging Using Row-Column Arrays

Row-column (RC) arrays have the potential to yield full 3-D ultrasound imaging with a greatly reduced number of elements compared to fully populated arrays. They, however, have several challenges due to their special geometry. This review article summarizes the current literature for RC imaging and demonstrates that full anatomic and functional imaging can attain a high quality using synthetic aperture (SA) sequences and modified delay-and-sum beamforming. Resolution can approach the diffraction limit with an isotropic resolution of half a wavelength with low sidelobe levels, and the field of view can be expanded by using convex or lensed RC probes. GPU beamforming allows for three orthogonal planes to be beamformed at 30 Hz, providing near real-time imaging ideal for positioning the probe and improving the operator's workflow. Functional imaging is also attainable using transverse oscillation and dedicated SA sequence for tensor velocity imaging for revealing the full 3-D velocity vector as a function of spatial position and time for both blood velocity and tissue motion estimation. Using RC arrays with commercial contrast agents can reveal super-resolution imaging (SRI) with isotropic resolution below [Formula: see text]. RC arrays can, thus, yield full 3-D imaging at high resolution, contrast, and volumetric rates for both anatomic and functional imaging with the same number of receive channels as current commercial 1-D arrays.

Performance Assessment of Row-Column Transverse Oscillation Tensor Velocity Imaging Using Computational Fluid Dynamics Simulation of Carotid Bifurcation Flow

In this work, the accuracy of row-column tensor velocity imaging (TVI), i.e., 3-D vector flow imaging (VFI) in 3-D space over time, is quantified on a complex, clinically relevant flow. The quantification is achieved by transferring the flow simulated using computational fluid dynamics (CFD) to a Field II simulation environment, and this allows for a direct comparison between the actual and estimated velocities. The carotid bifurcation flow simulations were performed with a peak inlet velocity of 80 cm/s, nonrigid vessel walls, and a flow cycle duration of 1.2 s. The flow was simulated from two observation angles, and it was acquired using a 3-MHz 62+62 row-column addressed array (RCA) at a pulse repetition frequency ( f_prf ) of 10 and 20 kHz. The tensor velocities were obtained at a frame rate of 208.3 Hz, at f_prf = 10 kHz , and the results from two velocity estimators were compared. The two estimators were the directional transverse oscillation (TO) cross correlation estimator and the proposed autocorrelation estimator. Linear regression between the actual and estimated velocity components yielded, for the cross correlation estimator, an R ² value in the range of 0.89-0.91, 0.46-0.77, and 0.91-0.97 for the x -, y -, and z -components, and 0.87-0.89, 0.40-0.83, and 0.91-0.96 when using the autocorrelation estimator. The results demonstrate that an RCA can, with just 62 receive channels, measure complex 3-D flow fields at a high volume rate.

Quantitative Blood Flow Measurements in the Common Carotid Artery: A Comparative Study of High-Frame-Rate Ultrasound Vector Flow Imaging, Pulsed Wave Doppler, and Phase Contrast Magnetic Resonance Imaging

V Flow is commercially developed by high-frame-rate ultrasound vector flow imaging. Compared to conventional color Doppler, V Flow is angle-independent and is capable of measuring both the magnitude and the direction of blood flow velocities. This paper aims to investigate the differences between V Flow and pulsed wave Doppler (PW) relative to phase contrast magnetic resonance imaging (PC-MRI), for the quantitative measurements of blood flow in common carotid arteries (CCA) and, consequently, to evaluate the accuracy of the new technique, V Flow. Sixty-four CCAs were measured using V Flow, PW, and PC-MRI. The maximum velocities, time-averaged mean (TAMEAN) velocities, and volume flow were measured using different imaging technologies. The mean error with standard deviation (Std), the median of absolute errors, and the r-values between V Flow and PC-MRI results for the maximum velocity, the TAMEAN velocity, and the volume flow measurements are {9.40% ± 14.91%; 11.84%; 0.84}, {21.52% ± 14.46%; 19.28%; 0.86}, and {−2.80% ± 14.01%; 10.38%; 0.7}, respectively, and are {53.44% ± 29.68%; 49.79%; 0.74}, {27.83% ± 31.60%; 23.83; 0.71}, and {21.01% ± 29.64%; 25.48%; 0.34}, respectively, for those between PW and PC-MRI. The boxplot, linear regression and Bland–Altman plots were performed for each comparison, which illustrated that the results measured via V Flow rather than via PW agreed more closely with those measured via PC-MRI.

Pressure Difference Estimation in Non-stenotic Carotid Bifurcation Phantoms Using Vector Flow Imaging

Local pressure differences estimated using vector flow imaging (VFI) and direct catheterization in seven carotid bifurcation phantoms were compared with simulated pressure fields. VFI correlated strongly with simulated peak pressure differences (r = 0.99, p < 0.00001), with an average overestimation of 12.3 Pa (28.6%). The range between the lowest and highest pressure difference of VFI underestimated simulations by 4.6 Pa (8.06%; r = 0.99, p < 0.0001). The catheter method exhibited no correlation (r = -0.09, p = 0.85). Ten repeated measurements on one phantom revealed a small standard deviation (SD) for VFI (SD = 8.4%, mean estimated SD = 11.5%), but not for the catheter method (SD = 785.6%). An in vivo peak systolic pressure difference of 97.9 Pa (estimated SD = 30%) was measured using VFI in one healthy individual. This study indicates that VFI pressure difference estimation is feasible in phantoms and in vivo and realistic estimates of the SD can be attained from the data.

Transthoracic Vector Flow Imaging in Pediatric Patients with Valvular Stenosis - A Proof of Concept Study

Purpose Continuous wave Doppler ultrasound is routinely used to detect cardiac valve stenoses. Vector flow imaging (VFI) is an angle-independent real-time ultrasound method that can quantify flow complexity. We aimed to evaluate if quantification of flow complexity could reliably assess valvular stenosis in pediatric patients.

Materials and Methods Nine pediatric patients with echocardiographically confirmed valvular stenosis were included in the study. VFI and Doppler measurements were compared with transvalvular peak-to-peak pressure differences derived from invasive endovascular catheterization.

Results Vector concentration correlated with the catheter measurements before intervention after exclusion of one outlier (r=−0.83, p=0.01), whereas the Doppler method did not (r=0.49, p=0.22). The change in vector concentration after intervention correlated strongly with the change in the measured catheter pressure difference (r=−0.86, p=0.003), while Doppler showed a tendency for a moderate correlation (r=0.63, p=0.07).

Conclusion Transthoracic flow complexity quantification calculated from VFI data is feasible and may be useful for assessing valvular stenosis severity in pediatric patients.

Turbulence Intensity as an Indicator for Ischemic Stroke in the Carotid Web

OBJECTIVE

Carotid web (CaW) is reported as one of the causes of cryptogenic strokes. However, its pathophysiology is not known, which makes clinicians debate CaW treatment. The current study presents the mechanism of ischemic stroke in CaW and finds a potential indicator for ischemic stroke in CaW.

METHODS

Computational fluid dynamics (CFD) was applied to computer-aided artificial human CaW models to investigate the correlation between stenosis rate and CaW angle. The CFD results were compared with stroke or symptomatic CaW including adequate reported images for anatomic analysis of the stenosis rate and CaW angle.

RESULTS

The CFD result of the CaW models shows that the velocity decreases and the turbulence generation of blood flow increases in the area under CaW (AuCaW). In AuCaW, turbulence intensity (TI) is high in the high stenosis rate and small angle conditions. Moreover, TI remained high in the AuCaW throughout the cardiac cycles. The analysis of 67 anatomic figures appeared in 36 articles merged with CFD results, and a CaW scoring system was introduced for clinicians to determine patients with CaW with a high risk of ischemic stroke.

CONCLUSIONS

This study is believed to be the first to apply TI to CaW. Because most ischemic strokes in CaW occurred at a high stenosis rate and small angle, TI may be used as a potential indicator of ischemic stroke in CaW. Furthermore, TI is easily used in the daily clinical field by applying CaW scores.

Common Carotid Artery Volume Flow: A Comparison Study between Ultrasound Vector Flow Imaging and Phase Contrast Magnetic Resonance Imaging

Volume flow estimation in the common carotid artery (CCA) can assess the absolute hemodynamic effect of a carotid stenosis. The aim of this study was to compare a commercial vector flow imaging (VFI) setup against the reference method magnetic resonance phase contrast angiography (MRA) for volume flow estimation in the CCA. Ten healthy volunteers were scanned with VFI and MRA over the CCA. VFI had an improved precision of 19.2% compared to MRA of 31.9% (p = 0.061). VFI estimated significantly lower volume flow than MRA (mean difference: 63.2 mL/min, p = 0.017), whilst the correlation between VFI and MRA was strong (R² = 0.81, p < 0.0001). A Bland–Altman plot indicated a systematic bias. After bias correction, the percentage error was reduced from 41.0% to 25.2%. This study indicated that a VFI setup for volume flow estimation is precise and strongly correlated to MRA volume flow estimation, and after correcting for the systematic bias, VFI and MRA become interchangeable.

Helical toroid phantom for 3D flow imaging investigations

The medical physics community has hitherto lacked an effective calibration phantom to holistically evaluate the performance of three-dimensional (3D) flow imaging techniques. Here, we present the design of a new omnidirectional, three-component (3-C) flow phantom whose lumen is consisted of a helical toroid structure (4 mm lumen diameter; helically winded for 5 revolutions over a torus with 10 mm radius; 5 mm helix radius). This phantom's intraluminal flow trajectory embraces all combinations of x, y, and z directional components, as confirmed using computational fluid dynamics (CFD) simulations. The phantom was physically fabricated via lost-core casting with polyvinyl alcohol cryogel (PVA) as the tissue mimic. 3D ultrasound confirmed that the phantom lumen expectedly resembled a helical toroid geometry. Pulsed Doppler measurements showed that the phantom, when operating under steady flow conditions (3 ml s^-1 flow rate), yielded flow velocity magnitudes that agreed well with those derived from CFD at both the inner torus (-47.6 ± 5.7 versus -52.0 ± 2.2 cm s^-1; mean ± 1 S.D.) and the outer torus (49.5 ± 4.2 versus 48.0 ± 1.7 cm s^-1). Additionally, 3-C velocity vectors acquired from multi-angle pulsed Doppler showed good agreement with CFD-derived velocity vectors (<7% and 10° difference in magnitude and flow angle, respectively). Ultrasound color flow imaging further revealed that the phantom's axial flow pattern was aligned with the CFD-derived flow profile. Overall, the helical toroid phantom has strong potential as an investigative tool in 3D flow imaging innovation endeavors, such as the development of flow vector estimators and visualization algorithms.

Adaptive transverse blood velocity estimation in medical ultrasound: A simulation study

Undoubtedly, highly valuable information about vascular anomalies is attained by the examination of the blood flow profile. The chief drawback of the conventional medical ultrasound in preparation of the blood periodogram is the measurement system shortcoming at the beam to flow angles near 90°. Recently, a method based on transverse oscillation (TO) approach, known as "Fourth-order estimation", has been developed to directly estimate the transverse power spectral density (PSD) of the fully transverse blood flow. One of the basic requirements to accomplish acceptable PSDs by this technique is the sufficiently large observation window. In this paper, two adaptive approaches for efficient estimation of the velocity spectrum of a fully transverse flow by a limited observation window length are described. The first proposed adaptive approach is based on the minimum variance adaptive spectral estimation in combination with the well-known TO technique (TO-MV). Then, by exploiting the eigenspace separation of the observed data to eliminate the contribution of the undesired components, the second technique (TO-EIBMV) is developed. The approaches are validated using Field II simulations for pulsating flow. The proposed methods are tested and compared to the conventional TO transverse spectral estimator by metrics of relative standard deviation (RSD) and relative bias (RB). One of the main achievements is the decrement of the required data samples for spectrogram estimation, which leads to a better temporal resolution. Moreover, for the analyzed adaptive techniques, the robustness of the estimation results for the beam to flow angles of 60-90° and vessel depths ranging from 20 mm to 60 mm are investigated.

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.

Flow Complexity Estimation in Dysfunctional Arteriovenous Dialysis Fistulas using Vector Flow Imaging

Non-invasive assessment is preferred for monitoring arteriovenous dialysis fistulas (AVFs). Vector concentration assesses flow complexity, which may correlate with stenosis severity. We determined whether vector concentration could assess stenosis severity in dysfunctional AVFs. Vector concentration was estimated in four stenotic phantoms at different pulse repetition frequencies. Spectral Doppler peak velocity and vector concentration were measured in 12 patients with dysfunctional AVFs. Additionally, 5 patients underwent digital subtraction angiography (DSA). In phantoms, vector concentration exhibited an inverse relationship with stenosis severity and was less affected by aliasing in severe stenoses. In nine stenoses of 5 patients undergoing DSA, vector concentration correlated strongly with stenosis severity (first stenosis: r = -0.73, p = 0.04; other stenoses; r = -0.69, p = 0.02) and mid-stenotic diameter (first stenosis: r = 0.87, p = 0.006; other stenoses: r = 0.70, p = 0.02) as opposed to peak velocities (p > 0.05). Vector concentration is less affected by aliasing in severe stenoses and correlates with DSA in patients with dysfunctional AVF.

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.

Real-Time Transthoracic Vector Flow Imaging of the Heart in Pediatric Patients

In children with congenital heart defects, Doppler ultrasound is the standard, bedside imaging modality. However, precise characterization of blood flow is challenging due to angle-dependent and one-dimensional velocity estimation. Contrast agent free Vector Flow Imaging is a new ultrasound technology that enables angle-independent visualization of the detailed flow field. Two piglets, one with normal cardiac anatomy and one with congenital heart disease comprised of valvular pulmonary stenosis, a dilated main pulmonary artery, and an incomplete atrioventricular canal defect, were imaged transthoracically and epicardially using a BK Ultrasound bk5000 with built-in vector flow imaging and a 5MHz linear probe. Subsequently, two children, one with normal cardiac anatomy and one with congenital heart disease comprised of aortic valve stenosis and coarctation of the aorta were imaged transthoracically. Transthoracic two-dimensional echocardiography and vector flow imaging were readily performed in both animals and were limited only by the geometry of the porcine thorax. In addition, transthoracic vector flow imaging was successfully performed in both children, and abnormal flow secondary to cardiac anomalies was visible. Adequate penetration was obtained to a depth of 6.5 cm. Our group has previously demonstrated for the first time that transthoracic vector flow imaging echocardiography is feasible and practicable in pediatric-sized patients, and this paper describes examples of these concepts and in-depth comparisons with traditional imaging modalities. This paper demonstrates that commercially available vector flow imaging technology can be utilized in pediatric cardiac applications as a bedside transthoracic imaging modality, providing advanced detail of blood flow patterns within the cardiac chambers, across valves, and in the great arteries.

Non-Invasive Assessment of Intravascular Pressure Gradients: A Review of Current and Proposed Novel Methods

Invasive catheterization is associated with a low risk of serious complications. However, although it is the gold standard for measuring pressure gradients, it induces changes to blood flow and requires significant resources. Therefore, non-invasive alternatives are urgently needed. Pressure gradients are routinely estimated non-invasively in clinical settings using ultrasound and calculated with the simplified Bernoulli equation, a method with several limitations. A PubMed literature search on validation of non-invasive techniques was conducted, and studies were included if non-invasively estimated pressure gradients were compared with invasively measured pressure gradients in vivo. Pressure gradients were mainly estimated from velocities obtained with Doppler ultrasound or magnetic resonance imaging. Most studies used the simplified Bernoulli equation, but more recent studies have employed the expanded Bernoulli and Navier⁻Stokes equations. Overall, the studies reported good correlation between non-invasive estimation of pressure gradients and catheterization. Despite having strong correlations, several studies reported the non-invasive techniques to either overestimate or underestimate the invasive measurements, thus questioning the accuracy of the non-invasive methods. In conclusion, more advanced imaging techniques may be needed to overcome the shortcomings of current methods.

4-D Intracardiac Ultrasound Vector Flow Imaging-Feasibility and Comparison to Phase-Contrast MRI

In vivo characterization of intracardiac blood velocity vector fields may provide new clinical information but is currently not available for bedside evaluation. In this paper, 4-D vector flow imaging for intracardiac flow assessment is demonstrated using a clinical ultrasound (US) system and a matrix array transducer, without the use of contrast agent. Two acquisition schemes were developed, one for full volumetric coverage of the left ventricle (LA) at 50 vps and a 3-D thick-slice setup with continuous frame acquisition (4000 vps), both utilizing ECG-gating. The 3-D vector velocity estimates were obtained using a novel method combining phase and envelope information. In vitro validation in a rotating tissue-mimicking phantom revealed velocity estimates in compliance with the ground truth, with a linear regression slope of 0.80, 0.77, and 1.03 for the , , and velocity components, and with standard deviations of 2.53, 3.19, and 0.95 cm/s, respectively. In vivo measurements in a healthy LV showed good agreement with PC-MRI. Quantitative analysis of energy loss (EL) and kinetic energy (KE) further showed similar trends, with peak KE at 1.5 and 2.4 mJ during systole and 3.6 and 3.1 mJ for diastole for US and PC-MRI. Similar for EL, 0.15- 0.2 and 0.7 mW was found during systole and 0.6 and 0.7 mW during diastole, for US and PC-MRI, respectively. Overall, a potential for US as a future modality for 4D cardiac vector flow imaging was demonstrated, which will be further evaluated in clinical studies.

Real-Time 2-D Phased Array Vector Flow Imaging

Echocardiography examination of the blood flow is currently either restricted to 1-D techniques in real-time or experimental offline 2-D methods. This paper presents an implementation of transverse oscillation for real-time 2-D vector flow imaging (VFI) on a commercial BK Ultrasound scanner. A large field-of-view (FOV) sequence for studying flow dynamics at 11 frames per second (fps) and a sequence for studying peak systolic velocities (PSVs) with a narrow FOV at 36 fps were validated. The VFI sequences were validated in a flow rig with continuous laminar parabolic flow and in a pulsating flow pump system before being tested in vivo, where measurements were obtained on two healthy volunteers. Mean PSV from 11 cycles was 155 cms^-1 with a precision of ±9.0% for the pulsating flow pump. In vivo, PSV estimated in the ascending aorta was 135 cms^-1 ± 16.9% for eight cardiac cycles. Furthermore, in vivo flow dynamics of the left ventricle and in the ascending aorta were visualized. In conclusion, angle independent 2-D VFI on a phased array has been implemented in real time, and it is capable of providing quantitative and qualitative flow evaluations of both the complex and fully transverse flow.