Numerical analysis of flow through a severely stenotic carotid artery bifurcation

The results of computational simulations may supplement MR and other in vivo diagnostic techniques to provide an accurate picture of the hemodynamics in particular vessels, which may help demonstrate the risks of embolism or plaque rupture posed by particular plaque deposits. In this study, a model based on an endarterectomy specimen of the plaque in a carotid bifurcation was examined. The flow conditions include steady flow at Reynolds numbers of 300, 600, and 900 as well as unsteady pulsatile flow. Both dynamic pressure and wall shear stress are very high, with shear values up to 70 N/m2, proximal to the stenosis throat in the internal carotid artery, and both vary significantly through the flow cycle. The wall shear stress gradient is also strong along the throat. Vortex shedding is observed downstream of the most severe occlusion. Two turbulence models, the Chien and Goldberg varieties of k-epsilon, are tested and evaluated for their relevance in this geometry. The Chien model better captures phenomena such as vortex shedding. The flow distal to stenosis is likely transitional, so a model that captures both laminar and turbulent behavior is needed.