Impact of minimal lumen segmentation uncertainty on patient-specific coronary simulations: A look at FFRCT

We examined the effect of minimal lumen segmentation uncertainty on Fractional Flow Reserve obtained from Coronary Computed Tomography AngiographyFFR CT . A total of 14 patient-specific coronary models with different stenosis locations and degrees of severity were enrolled in this study. The optimal segmented coronary lumens were disturbed using intra± 6 % and inter-operator± 15 % variations on the segmentation threshold.FFR CT was evaluated in each case by 3D-OD CFD simulations. The findings suggest that the sensitivity ofFFR CT to this type of uncertainty increases distally and with the stenosis severity. Cases with moderate or severe distal coronary lesions should undergo either exact and thorough segmentation operations or invasive FFR measurements, particularly if theFFR CT is close to the cutoff (0.80). Therefore, we conclude that it is crucial to consider the lesion's location and degree of severity when evaluatingFFR CT results.

Stroke risk prediction by color Doppler ultrasound of carotid artery-based deep learning using Inception V3 and VGG-16

Purpose

This study aims to automatically classify color Doppler images into two categories for stroke risk prediction based on the carotid plaque. The first category is high-risk carotid vulnerable plaque, and the second is stable carotid plaque.

Method

In this research study, we used a deep learning framework based on transfer learning to classify color Doppler images into two categories: one is high-risk carotid vulnerable plaque, and the other is stable carotid plaque. The data were collected from the Second Affiliated Hospital of Fujian Medical University, including stable and vulnerable cases. A total of 87 patients with risk factors for atherosclerosis in our hospital were selected. We used 230 color Doppler ultrasound images for each category and further divided those into the training set and test set in a ratio of 70 and 30%, respectively. We have implemented Inception V3 and VGG-16 pre-trained models for this classification task.

Results

Using the proposed framework, we implemented two transfer deep learning models: Inception V3 and VGG-16. We achieved the highest accuracy of 93.81% by using fine-tuned and adjusted hyperparameters according to our classification problem.

Conclusion

In this research, we classified color Doppler ultrasound images into high-risk carotid vulnerable and stable carotid plaques. We fine-tuned pre-trained deep learning models to classify color Doppler ultrasound images according to our dataset. Our suggested framework helps prevent incorrect diagnoses caused by low image quality and individual experience, among other factors.

A detailed analysis in thoracic aorta by means of the entropy generation rate: Prediction of the atherosclerotic lesion

A detailed numerical analysis is carried out in a real human thoracic aorta by means of the Computational Fluid Dynamics (CFD) for the prediction of the atherosclerosis lesion. Common hemodynamics parameters, such as, the oscillatory shear index (OSI) and the time average wall shear stress (TAWSS) are used for the prediction of the atherosclerosis lesion. Furthermore, the entropy generation rate is considered to obtain the main irreversibilities that occurs inside the thoracic aorta for the prediction of the atherosclerosis lesion. The model considers the blood flow inside the thoracic aorta in an unsteady state. The results show contours of velocity, streams lines, velocity profiles and the comparison of the hemodynamics parameters OSI versus TAWSS. Moreover, contours of the entropy generation rate are showed inside the aorta. The time averaged entropy generation rate (TAEGR) is obtained as a result of the entropy generation analysis. Finally, TAEGR index is compared and discussed with the common hemodynamics parameters, OSI and TAWSS. The accuracy to detect prone locations to atherosclerotic development in the real aorta using the TAEGR in comparison to the OSI and the TAWSS is in good agreement.

Magnetic resonance analysis of deep cerebral venous vasospasm after subarachnoid hemorrhage in rabbits

Objective

Arterial spasm is proved to be an inducer of cerebral ischemia and cerebral infarction, while when a venous spasm occurs, cerebral edema is seen to be caused by a disturbance in cerebral blood flow. However, it is unclear and unproven whether venous spasm occurs after subarachnoid hemorrhage (SAH). To provide the theoretical basis for treating cerebral vasospasm after SAH, magnetic resonance imaging (MRI) was employed to observe the changes in the diameter of deep cerebral veins in rabbits after SAH.

Methods

Fourteen New Zealand rabbits were randomly divided into the SAH group (n = 10) and the normal saline group (NS group, n = 4). Specifically, the SAH models were established by the ultrasound-guided double injections of blood into cisterna magna. Moreover, the MRI was performed to observe the changes in the diameter of deep cerebral veins (internal cerebral vein, basilar vein, and great cerebral vein) and basilar artery before modeling (0 d) and 1, 3, 5, 7, 9, and 11 d after modeling.

Results

In the SAH group, the diameter of the basilar artery showed no evident change on the 1st d. However, it became narrower obviously on the 3rd d and 5th d, and the stenosis degree was more than 30%. The diameter gradually relieved from 7th to 9th d, and finally returned to normal on the 11th d. Moreover, the diameter of the internal cerebral vein significantly narrowed on the 1st d, the stenosis degree of which was 19%; the stenosis then relieved slightly on the 3rd d (13%), reached the peak (34%) on the 5th d, and gradually relieved from 7th d to 11th d. Moreover, the stenosis degree of the basilar vein was 18% on the 1st d, 24% on the 3rd d, and reached the peak (34%) on the 5th d.

Conclusion

After SAH in rabbits, the cerebral vasospasm was seen to occur in the basilar artery, and likewise, spasmodic changes took place in the deep cerebral vein. Furthermore, the time regularity of spasmodic changes between the cerebral vein and basilar artery was of significant difference, indicating that the venous vasospasm resulted in active contraction.

COMPUTATIONAL STUDY OF THE INFLUENCE OF BIFURCATION ANGLE ON HAEMODYNAMICS AND OXYGEN TRANSPORT IN THE CAROTID BIFURCATION

In this study, blood flow associated with oxygen transport in the human carotid bifurcation was investigated numerically to assess the effects of bifurcation geometry on distribution and magnitude of the wall shear stress (WSS) and Sherwood number (Sh: dimensionless oxygen wall flux) at the favourable site of atherosclerotic lesion. Three-dimensional average models of the rigid-walled carotid bifurcation were constructed to perform simulations of steady blood flow under the wall boundary condition of a constant oxygen tension. The results demonstrated that changes in the bifurcation angle significantly altered the distribution of both the WSS and the Sh, even though the pattern of the axial flow was not very sensitive to the change in bifurcation angle. Flow with large inertia bifurcated at the flow divider and created a flow recirculation zone with low WSS and Sh on the outer wall of the internal carotid artery (ICA) sinus, where atherosclerotic lesions tend to develop. A wider bifurcation angle made the area of low Sh in the ICA sinus smaller, but the level of Sh along the outer wall of the ICA sinus extremely low. Another finding was that low Sh was associated with high WSS at the region distal to the ICA sinus. The Sh distribution did not readjust as fast as the WSS in this region, as reflected by the different rates of recovery of the WSS and Sh, thus uncoupling the transport process of oxygen transport from WSS.

Numerical analysis of the pressure drop across highly-eccentric coronary stenoses: application to the calculation of the fractional flow reserve

Background

Fractional flow reverse (FFR) is the gold standard assessment of the hemodynamic significance of coronary stenoses. However, it requires the catheterization of the coronary artery to determine the pressure waveforms proximal and distal to the stenosis. On the contrary, computational fluid dynamics enables the calculation of the FFR value from relatively non-invasive computed tomography angiography (CTA).

Methods

We analyze the flow across idealized highly-eccentric coronary stenoses by solving the Navier–Stokes equations. We examine the influence of several aspects (approximations) of the simulation method on the calculation of the FFR value. We study the effects on the FFR value of errors made in the segmentation of clinical images. For this purpose, we compare the FFR value for the nominal geometry with that calculated for other shapes that slightly deviate from that geometry. This analysis is conducted for a range of stenosis severities and different inlet velocity and pressure waveforms.

Results and conclusions

The errors made in assuming a uniform velocity profile in front of the stenosis, as well as those due to the Newtonian and laminar approximations, are negligible for stenosis severities leading to FFR values around the threshold 0.8. The limited resolution of the stenosis geometry reconstruction is the major source of error when predicting the FFR value. Both systematic errors in the contour detection of just 1-pixel size in the CTA images and a low-quality representation of the stenosis surface (coarse faceted geometry) may yield wrong outcomes of the FFR assessment for an important set of eccentric stenoses. On the contrary, the spatial resolution of images acquired with optical coherence tomography may be sufficient to ensure accurate predictions for the FFR value.

The effect of inlet and outlet boundary conditions in image-based CFD modeling of aortic flow

Background

Computational modeling of cardiovascular flow is a growing and useful field, but such simulations usually require the researcher to guess the flow’s inlet and outlet conditions since they are difficult and expensive to measure. It is critical to determine the amount of uncertainty introduced by these assumptions in order to evaluate the degree to which cardiovascular flow simulations are accurate. Our work begins to address this question by examining the sensitivity of flow to several different assumed velocity inlet and outlet conditions in a patient-specific aorta model.

Methods

We examined the differences between plug flow, parabolic flow, linear shear flows, skewed cubic flow profiles, and Womersley flow at the inlet. Only the shape of the inlet velocity profile was varied—all other parameters were identical among these simulations. Secondary flow in the form of a counter-rotating pair of vortices was also added to parabolic axial flow to study its effect on the solution. In addition, we examined the differences between two-element Windkessel, three element Windkessel and the outflow boundary conditions. In these simulations, only the outlet boundary condition was varied.

Results

The results show axial and in-plane velocities are considerably different close to the inlet for the cases with different inlet velocity profile shapes. However, the solutions are qualitatively similar beyond 1.75D, where D is the inlet diameter. This trend is also observed in other quantities such as pressure and wall shear stress. Normalized root-mean-square deviation, a measure of axial velocity magnitude differences between the different cases, generally decreases along the streamwise coordinate. The linear shear inlet velocity boundary condition and plug velocity boundary condition solution exhibit the highest time-averaged wall shear stress, approximately \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$8\%$$\end{document}8% higher than the parabolic inlet velocity boundary condition. Upstream of 1D from the inlet, adding secondary flow has a significant impact on temporal wall shear stress distributions. This is especially observable during diastole, when integrated wall shear stress magnitude varies about \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$26\%$$\end{document}26% between simulations with and without secondary flow. The results from the outlet boundary condition study show the Windkessel models differ from the outflow boundary condition by as much as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$18\%$$\end{document}18% in terms of time-averaged wall shear stress. Furthermore, normalized root-mean-square deviation of axial velocity magnitude, a measure of deviation between Windkessel and the outflow boundary condition, increases along the streamwise coordinate indicating larger variations near outlets.

Conclusion

It was found that the selection of inlet velocity conditions significantly affects only the flow region close to the inlet of the aorta. Beyond two diameters distal to the inlet, differences in flow solution are small. Although additional studies must be performed to verify this result, the data suggest that it is important to use patient-specific inlet conditions primarily if the researcher is concerned with the details of the flow very close to the inlet. Similarly, the selection of outlet conditions significantly affects the flow in the vicinity of the outlets. Upstream of five diameters proximal to the outlet, deviations between the outlet boundary conditions examined are insignificant. Although the inlet and outlet conditions only affect the flow significantly in their respective neighborhoods, our study indicates that outlet conditions influence a larger percentage of the solution domain.

The effect of heart failure and left ventricular assist device treatment on right ventricular mechanics: a computational study

Background and aims

Although it is important to analyze the hemodynamic factors related to the right ventricle (RV) after left ventricular assist device (LVAD) implantation, previous studies have focused only on the alteration of the ventricular shape and lack quantitative analysis of the various hemodynamic parameters. Therefore, we quantitatively analyzed various hemodynamic parameters related to the RV under normal, heart failure (HF), and HF incorporated with continuous flow LVAD therapy by using a computational model.

Methods

In this study, we combined a three-dimensional finite element electromechanical model of ventricles, which is based on human ventricular morphology captured by magnetic resonance imaging (MRI) with a lumped model of the circulatory system and continuous flow LVAD function in order to construct an integrated model of an LVAD implanted-cardiovascular system. To induce systolic dysfunction, the magnitude of the calcium transient function under HF condition was reduced to 70% of the normal value, and the time constant was reduced by 30% of the normal value.

Results

Under the HF condition, the left ventricular end systolic pressure decreased, the left ventricular end diastolic pressure increased, and the pressure in the right atrium (RA), RV, and pulmonary artery (PA) increased compared with the normal condition. The LVAD therapy decreased the end-systolic pressure of the LV by 41%, RA by 29%, RV by 53%, and PA by 71%, but increased the right ventricular ejection fraction by 52% and cardiac output by 40%, while the stroke work was reduced by 67% compared with the HF condition without LVAD. The end-systolic ventricular tension and strain decreased with the LVAD treatment.

Conclusion

LVAD enhances CO and mechanical unloading of the LV as well as those of the RV and prevents pulmonary hypertension which can be induced by HF.

Hemodynamic effects of perfusion level of peripheral ECMO on cardiovascular system

Background

Peripheral ECMO is an effective cardiopulmonary support in clinical. The perfusion level could directly influence the performances and complications. However, there are few studies on the effects of the perfusion level on hemodynamics of peripheral ECMO.

Methods

The geometric model of cardiovascular system with peripheral ECMO was established. The blood assist index was used to classify the perfusion level of the ECMO. The flow pattern from the aorta to the femoral artery and their branches, blood flow rate from aorta to brain and limbs, flow interface, harmonic index of blood flow, wall shear stress and oscillatory shear index were chosen to evaluate the hemodynamic effects of peripheral ECMO.

Results

The results demonstrated that the flow rate of aorta outlets increased and perfusion condition had been improved. And the average flow to the upper limbs and brain has a positive correlation with BAI (r = 0.037, p < 0.05), while there is a negative correlation with lower limbs (r = − 0.054, p < 0.05). The HI has negative correlation with BAI (p < 0.05, r < 0). The blood interface is further from the heart with the BAI decrease. And the average WSS has negative correlation with BAI (p < 0.05, r = − 0.983) at the bifurcation of femoral aorta and has positive correlation with BAI (p < 0.05, r = 0.99) at the inner aorta. The OSI under different BAI is higher (reaching 0.4) at the inner wall of the aortic arch, the descending aorta and the femoral access.

Conclusions

The pathogenesis of peripheral ECMO with different perfusion levels varies; its further research will be thorough and extensive.

Clinical validation and assessment of aortic hemodynamics using computational fluid dynamics simulations from computed tomography angiography

Background

Hemodynamic information including peak systolic pressure (PSP) and peak systolic velocity (PSV) carry an important role in evaluation and diagnosis of congenital heart disease (CHD). Since MDCTA cannot evaluate hemodynamic information directly, the aim of this study is to provide a noninvasive method based on a computational fluid dynamics (CFD) model, derived from multi-detector computed tomography angiography (MDCTA) raw data, to analyze the aortic hemodynamics in infants with CHD, and validate these results against echocardiography and cardiac catheter measurements.

Methods

This study included 25 patients (17 males, and 8 females; a median age of 2 years, range: 4 months–4 years) with CHD. All patients underwent both transthoracic echocardiography (TTE) and MDCTA within 2 weeks prior to cardiac catheterization. CFD models were created from MDCTA raw data. Boundary conditions were confirmed by lumped parameter model and transthoracic echocardiography (TTE). Peak systolic velocity derived from CFD models (PSV_CFD) was compared to TTE measurements (PSV_TTE), while the peak systolic pressure derived from CFD (PSP_CFD) was compared to catheterization (PSP_CC). Regions with low and high peak systolic wall shear stress (PSWSS) were also evaluated.

Results

PSV_CFD and PSP_CFD showed good agreements between PSV_TTE (r = 0.968, p < 0.001; mean bias = − 7.68 cm/s) and PSP_CC (r = 0.918, p < 0.001; mean bias = 1.405 mmHg). Regions with low and high PSWSS) can also be visualized. Skewing of velocity or helical blood flow was also observed at aortic arch in patients.

Conclusions

Our result demonstrated that CFD scheme based on MDCTA raw data is an accurate and convenient method in obtaining the velocity and pressure from aorta and displaying the distribution of PSWSS and flow pattern of aorta. The preliminary results from our study demonstrate the capability in combining clinical imaging data and novel CFD tools in infants with CHD and provide a noninvasive approach for diagnose of CHD such as coarctation of aorta in future.

Should fluid dynamics be included in computer models of RF cardiac ablation by irrigated-tip electrodes?

Background

Although accurate modeling of the thermal performance of irrigated-tip electrodes in radiofrequency cardiac ablation requires the solution of a triple coupled problem involving simultaneous electrical conduction, heat transfer, and fluid dynamics, in certain cases it is difficult to combine the software with the expertise necessary to solve these coupled problems, so that reduced models have to be considered. We here focus on a reduced model which avoids the fluid dynamics problem by setting a constant temperature at the electrode tip. Our aim was to compare the reduced and full models in terms of predicting lesion dimensions and the temperatures reached in tissue and blood.

Results

The results showed that the reduced model overestimates the lesion surface width by up to 5 mm (i.e. 70%) for any electrode insertion depth and blood flow rate. Likewise, it drastically overestimates the maximum blood temperature by more than 15 °C in all cases. However, the reduced model is able to predict lesion depth reasonably well (within 0.1 mm of the full model), and also the maximum tissue temperature (difference always less than 3 °C). These results were valid throughout the entire ablation time (60 s) and regardless of blood flow rate and electrode insertion depth (ranging from 0.5 to 1.5 mm).

Conclusions

The findings suggest that the reduced model is not able to predict either the lesion surface width or the maximum temperature reached in the blood, and so would not be suitable for the study of issues related to blood temperature, such as the incidence of thrombus formation during ablation. However, it could be used to study issues related to maximum tissue temperature, such as the steam pop phenomenon.

A new CFD based non-invasive method for functional diagnosis of coronary stenosis

Background

Accurate functional diagnosis of coronary stenosis is vital for decision making in coronary revascularization. With recent advances in computational fluid dynamics (CFD), fractional flow reserve (FFR) can be derived non-invasively from coronary computed tomography angiography images (FFR_CT) for functional measurement of stenosis. However, the accuracy of FFR_CT is limited due to the approximate modeling approach of maximal hyperemia conditions. To overcome this problem, a new CFD based non-invasive method is proposed.

Methods

Instead of modeling maximal hyperemia condition, a series of boundary conditions are specified and those simulated results are combined to provide a pressure-flow curve for a stenosis. Then, functional diagnosis of stenosis is assessed based on parameters derived from the obtained pressure-flow curve.

Results

The proposed method is applied to both idealized and patient-specific models, and validated with invasive FFR in six patients. Results show that additional hemodynamic information about the flow resistances of a stenosis is provided, which cannot be directly obtained from anatomy information. Parameters derived from the simulated pressure-flow curve show a linear and significant correlations with invasive FFR (r > 0.95, P < 0.05).

Conclusion

The proposed method can assess flow resistances by the pressure-flow curve derived parameters without modeling of maximal hyperemia condition, which is a new promising approach for non-invasive functional assessment of coronary stenosis.

Intrinsic Frequency Analysis and Fast Algorithms

Intrinsic Frequency (IF) has recently been introduced as an ample signal processing method for analyzing carotid and aortic pulse pressure tracings. The IF method has also been introduced as an effective approach for the analysis of cardiovascular system dynamics. The physiological significance, convergence and accuracy of the IF algorithm has been established in prior works. In this paper, we show that the IF method could be derived by appropriate mathematical approximations from the Navier-Stokes and elasticity equations. We further introduce a fast algorithm for the IF method based on the mathematical analysis of this method. In particular, we demonstrate that the IF algorithm can be made faster, by a factor or more than 100 times, using a proper set of initial guesses based on the topology of the problem, fast analytical solution at each point iteration, and substituting the brute force algorithm with a pattern search method. Statistically, we observe that the algorithm presented in this article complies well with its brute-force counterpart. Furthermore, we will show that on a real dataset, the fast IF method can draw correlations between the extracted intrinsic frequency features and the infusion of certain drugs.

Influence of vascular geometry on local hemodynamic parameters: phantom and small rodent study

BACKGROUND

Many studies have demonstrated that the geometry of the carotid bifurcation enables prediction of blood flow variation associated with atherosclerotic plaque formation. The phase angle between the arterial wall circumferential strain and its instantaneous wall shear stress is known as stress phase angle (SPA). This parameter is used to evaluate hemodynamic factors of atherogenesis. Note that SPA can be numerically computed for the purpose of locating atherosclerosis in different artery geometries. However, there is no experimental data to verify its role in the location of atherosclerosis in different artery geometries. In this study, we use an ultrasonic biomechanical method to experimentally evaluate the role of SPA for locating atherosclerosis in carotid bifurcation.

RESULTS

For carotid anthropomorphic vascular phantom experiments, the SPAs of common carotid arteries (CCAs), external carotid arteries (ECAs) and internal carotid arteries (ICAs) are - 148.53 ± 6.92°, - 153.95 ± 5.11°, and - 238.69 ± 1.72°, respectively. The corresponding SPAs are - 173.47 ± 0.065°, - 115.57 ± 4.83° and - 233.9 ± 8.12° for the polyvinyl alcohol (PVA-c) phantoms. In vivo mouse experiments indicated that the wall shear stress and circumferential strain were out of phase in the ICAs (- 280.08 ± 13.12°) to a greater extent as compared to CCAs (- 141.97 ± 8.03°) and ECAs (- 170.07 ± 9.24°).

CONCLUSIONS

The results suggested that SPA may be a useful indicator to locate the atherosclerosis position in carotid bifurcation.

Investigating the performance of four specific types of material grafts and their effects on hemodynamic patterns as well as on von Mises stresses in a grafted three-layer aortic model using fluid-structure interaction analysis

One of the important parts of the cardiac system is aorta which is the fundamental channel and supply of oxygenated blood in the body. Diseases of the aorta represent critical cardiovascular bleakness and mortality around the world. This study aims at investigation of hemodynamic parameters in a two-dimensional axisymmetric model of three-layer grafted aorta using fluid-structure interaction (FSI). It assumes that a damaged part of aorta, which may happen as a result of some diseases like aneurysm, dissection and post-stenotic dilatation, is replaced with a biomaterial graft. Four types of grafts materials so-called Polyurethane, Silicone rubber, Polytetrafluoroethylene (PTFE) and Dacron are considered in the present study. The assumption of linear elastic and isotropic material is set for the both aorta's wall and aforementioned grafts. Blood is considered as an incompressible and Newtonian fluid. The results indicate higher displacement in Polyurethane and silicone rubber in comparison with other two. Furthermore, results reveal that blood flow velocity has slightly higher values in PTFE and Dacron grafted models compared to Polyurethane and Silicone rubber ones. Even though there are some differences in hemodynamic patterns in these grafted models, they are not considerable as much as von Mises stresses across the graft-aorta intersections are. This study shows that the types of material grafts play an important role in the amount of stresses particularly at intersections of aorta and graft.

Computational and experimental assessment of influences of hemodynamic shear stress on carotid plaque

Background

Studies have identified hemodynamic shear stress as an important determinant of endothelial function and atherosclerosis. In this study, we assess the influences of hemodynamic shear stress on carotid plaques.

Methods

Carotid stenosis phantoms with three severity (30, 50, 70%) were made from 10% polyvinyl alcohol (PVA) cryogel. The phantoms were placed in a pulsatile flow loop with the same systolic/diastolic phase (35/65) and inlet flow rate (16 L/h). Ultrasonic particle imaging velocimetry (Echo PIV) and computational fluid dynamics (CFD) were used to calculate the velocity profile and shear stress distribution in the carotid stenosis phantoms. Inlet/outlet boundary conditions used in CFD were extracted from Echo PIV experiments to make sure that the results were comparable.

Results

Echo PIV and CFD results showed that velocity was largest in 70% than those in 30 and 50% at peak systole. Echo PIV results indicated that shear stress was larger in the upper wall and the surface of plaque than in the center of vessel. CFD results demonstrated that wall shear stress in the upstream was larger than in downstream of plaque. There was no significant difference in average velocity obtained by CFD and Echo PIV in 30% (p = 0.25). Velocities measured by CFD in 50% (93.01 cm/s) and in 70% (115.07 cm/s) were larger than those by Echo PIV in 50% (60.26 ± 5.36 cm/s) and in 70% (89.11 ± 7.21 cm/s).

Conclusions

The results suggested that Echo PIV and CFD could obtain hemodynamic shear stress on carotid plaques. Higher WSS occurred in narrower arteries, and the shoulder of plaque bore higher WSS than in bottom part.