Pulsatile pressure and flow in arterial stenoses simulated in a mathematical model

Abstract

This paper presents the detailed pulsatile pressure and flow velocity patterns inside an axis symmetric stenosis model with 75% constriction. The pressure and velocities have been calculated by solving the Navier-Stokes equations by the finite element method, the velocity profile in a straight tube caused by a pulsating driving pressure has been calculated first and then used as a boundary condition for the stenosis calculations. The results of the mathematical simulations of the stenosis model have been obtained in terms of velocity vectors, streamlines and isobars at 16 different instances in time, each 15 degrees apart during a cardiac cycle. The calculated velocity field shows that a vortex is developed at the wall distal to the stenosis as the velocity decreases from the peak systolic value. At the site of the vortex, a local pressure minimum is found due to the conversion of pressure to kinetic energy. When the flow is reversed, the reversal occurs first along the wall, thus forcing the vortex toward the the centre of the tube. As the reverse flow velocity increases, a vortex is also developed at the proximal site of the stenosis.