Finite element simulation of pulsatile flow through arterial stenosis.

The problem of blood flow through a stenosis is solved using the incompressible Navier-Stokes equations in a rigid circular tube presenting a partial occlusion. Calculations are based on a Galerkin finite element method. The time marching scheme employs a predictor-corrector technique using a variable time step. Results are obtained for steady and physiological pulsatile flows. Computational experiments analyse the effect of varying the degree of stenosis, the stricture length, the Reynolds number and Womersley number. The method gives results which agree well with previous computations for steady flows and experimental findings for steady and pulsatile flows.