Various issues relating to computational fluid dynamics simulations of carotid bifurcation flow based on models reconstructed from three-dimensional ultrasound images.

Computational fluid dynamics (CFD) flow simulation techniques have the potential to enhance understanding of how haemodynamic factors are involved in atherosclerosis. Recently, three-dimensional ultrasound has emerged as an alternative to other three-dimensional imaging techniques, such as magnetic resonance angiography (MRA). The method can be used to generate accurate vascular geometry suitable for CFD simulations and can be coupled with Doppler ultrasound to provide physiologically realistic flow boundary conditions. However, there are various ways to utilize the flow data acquired, possibly leading to different results regarding both flow and wall shear stress patterns. A disadvantage of three-dimensional ultrasound for imaging the carotid bifurcation has been established as being the scanning limitation of the jawbone position. This may make artificial extensions of the internal and/or external carotid arteries necessary, which in turn may influence the predicted flow patterns. Flow simulations were carried out for three outflow calculation schemes as well as four geometries with different extensions to the carotid daughter vessels. It was found that variation of flow patterns was more strongly influenced by the outflow conditions than by the extensions of the daughter vessels. Consequently, it is recommended that for future CFD simulations of carotid flow using three-dimensional ultrasound data, the outflow boundary conditions should rely on the most accurate measurement available, and flow data recorded in the common and internal carotid are considered more reliable than data from the external carotid. Even though the extended lengths of the daughter vessels have insignificant effects on the predicted haemodynamic parameters, it would be a safer option to extend the internal carotid by approximately three times the diameter of the common carotid artery.

Analysis of complex flow and the relationship between blood pressure, wall shear stress, and intima-media thickness in the human carotid artery.

BACKGROUND: Previous clinical studies have observed relationships between increased intima-media thickness (IMT) in the carotid artery, elevated blood pressure, and low wall shear stress (WSS) calculated from the Poiseuille equation. This study used numerical methods to more accurately determine WSS in the carotid artery and to investigate possible determinants of increased IMT. METHODS: IMT [common carotid artery (CCA) and bulb], CCA flow velocity, brachial systolic (SBP) and diastolic blood pressure (DBP), and carotid systolic pressure (cSBP) were measured in 14 healthy subjects (aged 44 +/- 16 yr). Flow patterns in the carotid bifurcation were determined by computational fluid dynamics (CFD) based on three-dimensional ultrasound geometry. Instantaneous and time-averaged wall shear stress (WSS(av)), oscillatory shear index (OSI), and wall shear stress angle gradients (WSSAG) were calculated. RESULTS: IMT was positively related to SBP, DBP, cSBP, and WSSAG and inversely related to WSS(av) in the CCA. In the bulb, IMT was positively related to SBP and cSBP but was not significantly related to WSS(av) or WSSAG. IMT was unrelated to OSI in both the CCA and the bulb. CONCLUSION: Increased carotid artery IMT in healthy subjects with no evidence of focal plaques is primarily a response to elevated pressure.

Reproducibility study of magnetic resonance image-based computational fluid dynamics prediction of carotid bifurcation flow.

The importance of shear stress in the initiation and progression of atherosclerosis has been recognized for some time. A novel way to quantify wall shear stress under physiologically realistic conditions is to combine magnetic resonance imaging (MRI) and computational fluid dynamics. The present study aims to investigate the reproducibility of the simulated flow by using this combined approach. The right carotid bifurcations of eight healthy subjects were scanned twice with MRI within a few weeks. Three-dimensional geometries of the vessels were reconstructed for each scan and each subject. Pulsatile flows through these models were calculated to assess errors associated with the predicted flow parameters. This was done by comparing various wall shear stress indices, including the time-averaged wall shear stress (WSS), oscillating shear index (OSI), WSS Gradients (WSSG) and WSS Angle Deviation (WSSAD). Qualitatively, all the wall shear parameters proved to be highly reproducible. Quantitatively, the reproducibility was over 90% for OSI and WSSAD, but less impressive (60%) for other parameters. Our results indicated that WSS and WSSG values were extremely sensitive to subtle variations in local geometry and mesh design, particularly in regions around the bifurcation apex where WSS values were high and least reproducible.

Accuracy and reproducibility of CFD predicted wall shear stress using 3D ultrasound images.

Computational fluid dynamics (CFD) flow simulation techniques have the potential to enhance our understanding of how haemodynamic factors are involved in atherosclerosis. Recently, 3D ultrasound has emerged as an alternative to other 3D imaging techniques, such as magnetic resonance angiography (MRA). The method can be used to generate realistic vascular geometry suitable for CFD simulations. In order to assess accuracy and reproducibility of the procedure from image acquisition to reconstruction to CFD simulation, a human carotid artery bifurcation phantom was scanned three times using 3D ultrasound. The geometry was reconstructed and flow simulations were carried out on the three sets as well as on a model generated using computer aided design (CAD) from the geometric information given by the manufacturer. It was found that the three reconstructed sets showed good reproducibility as well as satisfactory quantitative agreement with the CAD model. Analyzing two selected locations probably representing the 'worst cases,' accuracy comparing ultrasound and CAD reconstructed models was estimated to be between 7.2% and 7.7% of the maximum instantaneous WSS and reproducibility comparing the three scans to be between 8.2% and 10.7% of their average maximum.

Influence of head position on carotid hemodynamics in young adults.

Studies in adults have shown marked changes in geometry and relative positions of the carotid arteries when rotating the head. The aim of this study was to quantify the change in geometry and analyze its effect on carotid hemodynamics as a result of head rotation. The right carotid arteries of nine young adult subjects were investigated in supine position with straight and left turned head positions, respectively. The three-dimensional (3D) carotid geometry was reconstructed by using 3D ultrasound (3D US), and the carotid hemodynamics were calculated by combining 3D US with computational fluid dynamics. It was observed that cross-sectional areas and shapes did not change markedly with head rotation, but carotid vessel center lines altered with planarification of the common carotid artery as a main feature (P < 0.05). Measured common carotid flow rates changed significantly at the individual level when the head was turned, but on the average, the change in mean common carotid flow rate was relatively small (0.37 +/- 1.11 ml/s). The effect of the altered center lines and flow rates on the atherogenic nature of the carotid bifurcation was evaluated by using calculated hemodynamic wall parameters, such as wall shear stress (WSS) and oscillatory shear index (OSI). It was found that WSS and OSI patterns changed significantly with head rotation, but the variations were very subject dependent and could not have been predicted without assessing the altered geometry and flow of the carotid bifurcation for individual cases. This study suggests that there is a need for standardization of the choice of head position in the 3D US scan protocol, and that carotid stents and emboli diverters should be studied in different head positions.