Finite element analysis of nonlinear pulsatile suspension flow dynamics in blood vessels with aneurysm

A nonlinear pulsatile suspension flow in a dilated vessel is numerically analysed. Two sets of highly coupled nonlinear partial differential equations governing the suspension flow are numerically solved, to simulate the suspension flow dynamics. A transient velocity-pressure (UVP) finite element method (FEM) and a stable time integration scheme, based on a predictor-corrector strategy, with constant error monitoring are employed in the flow analysis. The pulsatile suspension flow is characterized by analysing the flow, pressure and stress fields. Effects of the nonlinear particulate phase on the nonlinear suspending fluid phase are brought out by comparing the suspension flow results with those of homogeneous flow. Particles are seen to dampen the flow velocity, wall and central axis pressure, pressure gradient and wall shear stress. time-dependent recirculation regions which are sensitive to the presence of particles are seen in the dilated portion of the vessel. These recirculation regions favour thrombogenesis. The nonlinear effects due to the vessel geometry and those due to the convective terms dominate the dampening effect of the particles. These nonlinear effects are depicted through the transverse velocity and pressure plots. Wall shear stresses of suspension flow are not only high but also alternate in direction.