Echo Particle Image Velocimetry for Estimation of Carotid Artery Wall Shear Stress: Repeatability, Reproducibility and Comparison with Phase-Contrast Magnetic Resonance Imaging

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.

A comparison of intima media thickness in the common carotid artery, the bulb and plaque area as predictions of incident atherosclerotic events

BACKGROUND AND AIMS

There is a debate on how to evaluate carotid artery intima-media thickness (IMT). We here compared IMT of the common carotid artery (CCA) and bulb with plaque area regarding incident atherosclerotic disease.

METHODS

In the PIVUS study (age 70 at baseline, 53% women, n = 856), IMT-CCA, IMT-bulb and plaque area were measured at ages 70, 75 and 80 years and these three measurements were used in updated Cox proportional hazard analysis.

RESULTS

Over 15 years follow-up, 135 individuals experienced a first-time atherosclerotic disease (myocardial infarction or ischemic stroke). IMT-CCA was not significantly related to this composite endpoint (p = 0.10). IMT-bulb was significantly related to the endpoint (p = 0.003), but this relationship was attenuated following adjustment for CVD risk factors (p = 0.02). On the contrary, plaque area was consistently related to incident atherosclerotic disease also following adjustment for CVD risk factors (p<0.001). When added on top of traditional risk factors, both IMT-bulb and plaque area, but not IMT-CCA, improved the discrimination compared to the traditional risk factors (+5.2%, p = 0.0026 for IMT-bulb, +3.8%, p = 0.013 for plaque area and 0.0% for IMT-CCA).

CONCLUSION

In elderly subjects, both IMT-bulb and plaque area improved the discrimination regarding incident atherosclerotic disease when added to traditional risk factors. This was not seen for IMT-CCA. IMT-CCA was therefore inferior compared to the other two carotid artery ultrasonographic measurements in this sample of elderly subjects.

Imperatives of Mathematical Model of Arterial Blood Dynamics for Interpretation of Doppler Velocimetry: A Narrative Review

Clinicians frequently study arterial Doppler velocimetric waveforms depicted by Doppler sonography of the kidneys, the heart, the brain, and the feto-maternal circulation to assess the well-being of the aforementioned vital organs. The waveform interpretation of the Doppler indices can be studied using a mathematical model. The developed models serve as teaching tools and for easy comprehension of the regulatory mechanism of the organs. It will also obtain accurate wall shear stress (WSS) and likely atherosclerotic sites can be predicted early. The aim of this review is to reveal the imperatives of mathematical models in the study of the physical interpretation of Doppler velocimetry. The models will explore sonographic Doppler velocimetry and computational fluid dynamics (CFD) in determining the segments of the arteries that are prone to the development of atheromatous plaque. It will be achieved by comparing and computing the measurement differences of the WSS. A thorough literature review was carried out between 1971 and 2021 on the mathematical modeling of blood dynamics and Doppler velocimetry of different blood vessels, across various electronic databases including NC AHEC Digital Library, PUBMED, ERIC, MEDLINE, Free Medical Journals, and EMBASE. The results of the literature search were presented using the PRISMA flow chat. The narrative review of the mathematical models of arterial blood dynamics is based on incompressible Navier-Stokes equations, the Windkessel model, and CFD. It was deduced that the blood flow velocity decreased with time across the varying frequency from 0.2Hz to 0.50Hz in the interlobar arterial channels. The review also revealed that adult humans' Doppler indices of the renal-interlobar artery agree with developed models of renal interlobar arterial blood dynamics. The mathematical model measurements of the great vessels matched the sonographic Doppler velocimetry with <15% variation. In our fast-paced world of epidemiological transition, the imperatives of mathematical modeling of arterial flow dynamics based on the Navier-Stokes equations to represent various physiologic and pathologic situations cannot be overstated. The practical consequences include the possibility of mathematical models to acquire precise WSS distribution and early detection of potential atherosclerotic sites during cardiovascular Doppler sonography.

Dealiasing High-Frame-Rate Color Doppler Using Dual-Wavelength Processing

Doppler ultrasound is the premier modality to analyze blood flow dynamics in clinical practice. With conventional systems, Doppler can either provide a time-resolved quantification of the flow dynamics in sample volumes (spectral Doppler) or an average Doppler velocity/power [color flow imaging (CFI)] in a wide field of view (FOV) but with a limited frame rate. The recent development of ultrafast parallel systems made it possible to evaluate simultaneously color, power, and spectral Doppler in a wide FOV and at high-frame rates but at the expense of signal-to-noise ratio (SNR). However, like conventional Doppler, ultrafast Doppler is subject to aliasing for large velocities and/or large depths. In a recent study, staggered multi-pulse repetition frequency (PRF) sequences were investigated to dealias color-Doppler images. In this work, we exploit the broadband nature of pulse-echo ultrasound and propose a dual-wavelength approach for CFI dealiasing with a constant PRF. We tested the dual-wavelength bandpass processing, in silico, in laminar flow phantom and validated it in vivo in human carotid arteries ( n = 25 ). The in silico results showed that the Nyquist velocity could be extended up to four times the theoretical limit. In vivo, dealiased CFI were highly consistent with unfolded Spectral Doppler ( r²=0.83 , y=1.1x+0.1 , N=25 ) and provided consistent vector flow images. Our results demonstrate that dual-wavelength processing is an efficient method for high-velocity CFI.

Blood Flow Quantification in Peripheral Arterial Disease: Emerging Diagnostic Techniques in Vascular Surgery

The assessment of local blood flow patterns in patients with peripheral arterial disease is clinically relevant, since these patterns are related to atherosclerotic disease progression and loss of patency in stents placed in peripheral arteries, through mechanisms such as recirculating flow and low wall shear stress (WSS). However, imaging of vascular flow in these patients is technically challenging due to the often complex flow patterns that occur near atherosclerotic lesions. While several flow quantification techniques have been developed that could improve the outcomes of vascular interventions, accurate 2D or 3D blood flow quantification is not yet used in clinical practice. This article provides an overview of several important topics that concern the quantification of blood flow in patients with peripheral arterial disease. The hemodynamic mechanisms involved in the development of atherosclerosis and the current clinical practice in the diagnosis of this disease are discussed, showing the unmet need for improved and validated flow quantification techniques in daily clinical practice. This discussion is followed by a showcase of state-of-the-art blood flow quantification techniques and how these could be used before, during and after treatment of stenotic lesions to improve clinical outcomes. These techniques include novel ultrasound-based methods, Phase-Contrast Magnetic Resonance Imaging (PC-MRI) and Computational Fluid Dynamics (CFD). The last section discusses future perspectives, with advanced (hybrid) imaging techniques and artificial intelligence, including the implementation of these techniques in clinical practice.

Blood Flow Quantification in Peripheral Arterial Disease: Emerging Diagnostic Techniques in Vascular Surgery

The assessment of local blood flow patterns in patients with peripheral arterial disease is clinically relevant, since these patterns are related to atherosclerotic disease progression and loss of patency in stents placed in peripheral arteries, through mechanisms such as recirculating flow and low wall shear stress (WSS). However, imaging of vascular flow in these patients is technically challenging due to the often complex flow patterns that occur near atherosclerotic lesions. While several flow quantification techniques have been developed that could improve the outcomes of vascular interventions, accurate 2D or 3D blood flow quantification is not yet used in clinical practice. This article provides an overview of several important topics that concern the quantification of blood flow in patients with peripheral arterial disease. The hemodynamic mechanisms involved in the development of atherosclerosis and the current clinical practice in the diagnosis of this disease are discussed, showing the unmet need for improved and validated flow quantification techniques in daily clinical practice. This discussion is followed by a showcase of state-of-the-art blood flow quantification techniques and how these could be used before, during and after treatment of stenotic lesions to improve clinical outcomes. These techniques include novel ultrasound-based methods, Phase-Contrast Magnetic Resonance Imaging (PC-MRI) and Computational Fluid Dynamics (CFD). The last section discusses future perspectives, with advanced (hybrid) imaging techniques and artificial intelligence, including the implementation of these techniques in clinical practice.

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.

Measuring wall shear stress distribution in the carotid artery in an African population: Computational fluid dynamics versus ultrasound doppler velocimetry

INTRODUCTION

Computational fluid dynamics (CFD) and ultrasound Doppler velocimetry are diagnostic tools useful for determining carotid artery segments susceptible to atheromatous plaque development. This study computes and compares the difference in Wall Shear Stress (WSS) measurements between these two methods.

METHODS

The carotid artery of 204 volunteers selected using simple random sampling were scanned using standard carotid doppler protocols. Four segments of the carotid artery - the common, internal, external carotid, and the carotid bulb were sonographically assessed. The intima-media thickness, diameter, peak systolic velocity, and end-diastolic velocity were measured at a point 2 cm away from the carotid bifurcation for the three segments, while the carotid bulb was measured at the bifurcation. A 2D incompressible Navier-Stokes Equation for modelling Newtonian, pulsatile, and laminar flow in a viscoelastic pipe was applied to model velocity flow across the carotid artery using COMSOL software. WSS values were computed for experimental and CFD measurements and the results were compared.

RESULTS

The WSS values generated by the model had respectively peak and average values of 19.81 N/cm² and 15.76 ± 1.81 N/cm² for the common carotid, 10.77 N/cm² and 7.57 ± 1.66 N/cm² for the internal carotid, 11.51 N/cm² and 8.05 ± 1.65 N/cm² for the external carotid, 37.55 N/cm² and 26.55 ± 6.62 N/cm² for the carotid bifurcation, 1.39 N/cm² and 3.13 ± 1.34 N/cm² for the carotid bulb. The model measurements matched doppler velocimetry measurements with <15% variation.

CONCLUSION

Model based WSS values were higher but comparable with doppler velocimetry measurements. The carotid bulb had low WSS and is therefore the segment highly disposed to atheromatous plaque formation.

IMPLICATIONS FOR PRACTICE

Subject-specific mathematical models could be incorporated during cardiovascular scan work up for accurate WSS distribution and early prediction of possible atherosclerotic sites.

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.

Influence of malformation of right coronary artery originating from the left sinus in hemodynamic environment

Background

The anomalous origin of the right coronary artery (RCA) from the left coronary artery sinus (AORL) is one of the abnormal origins of the coronary arteries. Most of these issues rarely have any effects on human health, but some individuals may exhibit symptoms, such as myocardial ischemia or even sudden death. Recently, researchers have investigated the AORL through clinical cases, but studies based on computational fluid dynamics (CFD) have rarely been reported. In this study, the hemodynamic changes between the normal origin of the RCA and the AORL are compared based on numerical simulation results.

Methods

Realistic three-dimensional (3D) models of the 16 normal right coronary arteries and 26 abnormal origins of the RCAs were constructed, respectively. The blood flow was numerically simulated using the ANSYS software. This study used a one-way fluid–solid coupling finite element model, wherein the blood is assumed to be an incompressible Newtonian fluid, and the vessel is assumed to be made of an isotropic linear elastic material.

Results

The cross-sectional area differences between the inlet of the normal group and that of the abnormal group were significant (P < 0.0001). Moreover, there were significant differences in the volumetric flow (P = 0.0001) and pressure (P = 0.0002). Positive correlation exists for the ratio of the cross-sectional area of the RCA to the inlet area of the ascending aorta (AAO), and the ratio of the inlet volumetric flow of the RCA to the volumetric flow of the AAO, in the normal (P = 0.0001, r = 0.8178) and abnormal (P = 0.0033, r = 0.6107) groups.

Conclusion

This study demonstrates that the cross-sectional area of the AORL inlet may cause ischemia symptoms. The results obtained by this study may contribute to the further understanding of the clinical symptoms of the AORL based on the hemodynamics.

Determining Haemodynamic Wall Shear Stress in the Rabbit Aorta In Vivo Using Contrast-Enhanced Ultrasound Image Velocimetry

Abnormal blood flow and wall shear stress (WSS) can cause and be caused by cardiovascular disease. To date, however, no standard method has been established for mapping WSS in vivo. Here we demonstrate wide-field assessment of WSS in the rabbit abdominal aorta using contrast-enhanced ultrasound image velocimetry (UIV). Flow and WSS measurements were made independent of beam angle, curvature or branching. Measurements were validated in an in silico model of the rabbit thoracic aorta with moving walls and pulsatile flow. Mean errors over a cardiac cycle for velocity and WSS were 0.34 and 1.69%, respectively. In vivo time average WSS in a straight segment of the suprarenal aorta correlated highly with simulations (PC = 0.99) with a mean deviation of 0.29 Pa or 5.16%. To assess fundamental plausibility of the measurement, UIV WSS was compared to an analytic approximation derived from the Poiseuille equation; the discrepancy was 17%. Mapping of WSS was also demonstrated in regions of arterial branching. High time average WSS (TAWSS_xz = 3.4 Pa) and oscillatory flow (OSI_xz = 0.3) were observed near the origin of conduit arteries. In conclusion, we have demonstrated that contrast-enhanced UIV is capable of measuring spatiotemporal variation in flow velocity, arterial wall location and hence WSS in vivo with high accuracy over a large field of view.

Ultrasound Molecular Imaging of Atherosclerosis With Nanobodies: Translatable Microbubble Targeting Murine and Human VCAM (Vascular Cell Adhesion Molecule) 1

OBJECTIVE

Contrast-enhanced ultrasound molecular imaging (CEUMI) of endothelial expression of VCAM (vascular cell adhesion molecule)-1 could improve risk stratification for atherosclerosis. The microbubble contrast agents developed for preclinical studies are not suitable for clinical translation. Our aim was to characterize and validate a microbubble contrast agent using a clinically translatable single-variable domain immunoglobulin (nanobody) ligand. Approach and Results: Microbubble with a nanobody targeting VCAM-1 (MB_cAbVcam1-5) and microbubble with a control nanobody (MB_VHH2E7) were prepared and characterized in vitro. Attachment efficiency to VCAM-1 under continuous and pulsatile flow was investigated using activated murine endothelial cells. In vivo CEUMI of the aorta was performed in atherosclerotic double knockout and wild-type mice after injection of MB_cAbVcam1-5 and MB_VHH2E7. Ex vivo CEUMI of human endarterectomy specimens was performed in a closed-loop circulation model. The surface density of the nanobody ligand was 3.5×10⁵ per microbubble. Compared with MB_VHH2E7, MB_cAbVcam1-5 showed increased attachment under continuous flow with increasing shear stress of 1-8 dynes/cm² while under pulsatile flow attachment occurred at higher shear stress. CEUMI in double knockout mice showed signal enhancement for MB_cAbVcam1-5 in early (P=0.0003 versus MB_VHH2E7) and late atherosclerosis (P=0.007 versus MB_VHH2E7); in wild-type mice, there were no differences between MB_cAbVcam1-5 and MB_VHH2E7. CEUMI in human endarterectomy specimens showed a 100% increase in signal for MB_cAbVcam1-5versus MB_VHH2E7 (20.6±27.7 versus 9.6±14.7, P=0.0156).

CONCLUSIONS

CEUMI of the expression of VCAM-1 is feasible in murine models of atherosclerosis and on human tissue using a clinically translatable microbubble bearing a VCAM-1 targeted nanobody.

Arterial wall shear rate response to reactive hyperaemia is markedly different between young and older humans

KEY POINTS

The vasodilatory response to reactive hyperaemia is impaired with advancing age, but it is unclear whether this is because of an altered wall shear rate (WSR) stimulus or an altered flow-mediated dilatation (FMD) response. Using new technology that allows detailed WSR measurement, we assessed the WSR-FMD response in healthy older people. Our data show that older people have a markedly altered and diminished WSR response to reactive hyperaemia compared to young people, but reduced WSR alone does not fully explain reduced FMD. In young people, WSR appears to be coupled to FMD but, by age ∼65 years, the arterial vasodilatory response has begun to uncouple from the WSR stimulus. These findings point to the importance and utility of comprehensively characterizing the WSR-FMD response when using reactive hyperaemia to assess vascular function, as well as giving new insight into the age-related alteration in vascular function.

ABSTRACT

The vasodilatory response to reactive hyperaemia is impaired with age, but it is unknown whether this is because of an altered wall shear rate (WSR) stimulus or an altered flow-mediated dilatation (FMD) response to the WSR stimulus. Inherent difficulties in measuring blood flow velocity close to the arterial wall have prevented detailed assessment of the WSR-FMD response. Using an enhanced multigate spectral Doppler ultrasound system (ultrasound advanced open platform), we aimed to produce new data on the WSR-FMD relationship in healthy older adults. Sixty healthy people, comprising 28 young (27.5 ± 5.5 years) and 32 older (64.9 ± 3.7 years) individuals, underwent FMD assessment. Raw data were post-processed using custom-designed software to obtain WSR and diameter parameters. The data revealed that older people have a much altered and diminished WSR response to reactive hyperaemia compared to younger people [e.g. WSR peak: 622 (571-673) vs. 443 (396-491) 1/s in young and older respectively; P < 0.05]. However, reduced WSR alone does not appear to fully explain the reduced FMD response in older people because associations between WSR and FMD were few and weak. This was in contrast to young adults, where associations were strong. We conclude that WSR during FMD is much altered and diminished in older people, and there appears to be an 'uncoupling' of WSR from FMD in older people that may reflect a loss of precision in the reactive hyperaemia stimulus-response relationship. These findings also point to the importance and utility of comprehensively characterizing the WSR-FMD response when using reactive hyperaemia to assess vascular function.

Accuracy and Precision of a Plane Wave Vector Flow Imaging Method in the Healthy Carotid Artery

The objective of the study described here was to investigate the accuracy and precision of a plane wave 2-D vector flow imaging (VFI) method in laminar and complex blood flow conditions in the healthy carotid artery. The approach was to study (i) the accuracy for complex flow by comparing the velocity field from a computational fluid dynamics (CFD) simulation to VFI estimates obtained from the scan of an anthropomorphic flow phantom and from an in vivo scan; (ii) the accuracy for laminar unidirectional flow in vivo by comparing peak systolic velocities from VFI with magnetic resonance angiography (MRA); (iii) the precision of VFI estimation in vivo at several evaluation points in the vessels. The carotid artery at the bifurcation was scanned using both fast plane wave ultrasound and MRA in 10 healthy volunteers. The MRA geometry acquired from one of the volunteers was used to fabricate an anthropomorphic flow phantom, which was also scanned using the fast plane wave sequence. The same geometry was used in a CFD simulation to calculate the velocity field. Results indicated that similar flow patterns and vortices were estimated with CFD and VFI in the phantom for the carotid bifurcation. The root-mean-square difference between CFD and VFI was within 0.12 m/s for velocity estimates in the common carotid artery and the internal branch. The root-mean-square difference was 0.17 m/s in the external branch. For the 10 volunteers, the mean difference between VFI and MRA was -0.17 m/s for peak systolic velocities of laminar flow in vivo. The precision in vivo was calculated as the mean standard deviation (SD) of estimates aligned to the heart cycle and was highest in the center of the common carotid artery (SD = 3.6% for velocity magnitudes and 4.5° for angles) and lowest in the external branch and for vortices (SD = 10.2% for velocity magnitudes and 39° for angles). The results indicate that plane wave VFI measures flow precisely and that estimates are in good agreement with a CFD simulation and MRA.

Measurement of Wall Shear Stress Exerted by Flowing Blood in the Human Carotid Artery: Ultrasound Doppler Velocimetry and Echo Particle Image Velocimetry

Vascular endothelial cells lining the arteries are sensitive to wall shear stress (WSS) exerted by flowing blood. An important component of the pathophysiology of vascular diseases, WSS is commonly estimated by centerline ultrasound Doppler velocimetry (UDV). However, the accuracy of this method is uncertain. We have previously validated the use of a novel, ultrasound-based, particle image velocimetry technique (echo PIV) to compute 2-D velocity vector fields, which can easily be converted into WSS data. We compared WSS data derived from UDV and echo PIV in the common carotid artery of 27 healthy participants. Compared with echo PIV, time-averaged WSS was lower using UDV (28 ± 35%). Echo PIV revealed that this was due to considerable spatiotemporal variation in the flow velocity profile, contrary to the assumption that flow is steady and the velocity profile is parabolic throughout the cardiac cycle. The largest WSS underestimation by UDV was found during peak systole (118 ± 16%) and the smallest during mid-diastole (4.3± 46%). The UDV method underestimated WSS for the accelerating and decelerating systolic measurements (68 ± 30% and 24 ± 51%), whereas WSS was overestimated for end-diastolic measurements (-44 ± 55%). Our data indicate that UDV estimates of WSS provided limited and largely inaccurate information about WSS and that the complex spatiotemporal flow patterns do not fit well with traditional assumptions about blood flow in arteries. Echo PIV-derived WSS provides detailed information about this important but poorly understood stimulus that influences vascular endothelial pathophysiology.

Advanced Pediatric Neurosonography Techniques: Contrast-Enhanced Ultrasonography, Elastography, and Beyond

Recent technical advances in neurosonography continue broadening the diagnostic utility, sensitivity, and specificity of ultrasound for detecting intracranial abnormalities bed side. The clinical and functional applications of neurosonography have significantly expanded since the 1980s when transcranial Doppler sonography first allowed anatomic and hemodynamic delineation of the intracranial vessels through the thin temporal skull. In the past few years, contrast-enhanced ultrasonography, elastography, 3D/4D reconstruction tools, and high-resolution microvessel imaging techniques have further enhanced the diagnostic significance of neurosonography. Given these advances, a thorough familiarity with these new techniques and devices is crucial for a successful clinical application allowing improved patient care. It is essential that future neurosonography studies compare these advanced techniques against the current "gold standard" computed tomography and magnetic resonance imaging to assure the accuracy of their diagnostic potential. This review will provide a comprehensive update on currently available advanced neurosonography techniques.