Numerical analysis of blood flow through a stenosed artery using a coupled, multiscale simulation method

A global system model of the systemic circulation is combined with a local finite element solution to simulate blood flow in a stenosed coronary artery. Local fluid dynamic issues arise in connection with the detailed flow patterns within the stenosed coronary artery while the global system model is used to simulate the response of the rest of the circulation to the local perturbation. A PISO type finite element technique is employed to compute the local blood flow. The Navier-Stokes equations are solved with the assumption of viscous incompressible flow across the stenosed coronary artery. A detailed lumped parameter model simulates the characteristics of the coronary circulation and is imbedded in a coarse-grained lumped parameter model of the entire cardiovascular system. These two methods are coupled in that the lumped parameter calculations provide the time-dependent boundary conditions for the local finite element calculation. In turn, the local fluid dynamical computation provides estimates for the pressure drop across the stenosis, which is subsequently used to refine the lumped parameter calculation. Results are obtained for an axisymmetric coronary artery model with a stenosis of 90% area reduction over one cardiac cycle. Numerical results show that the flow rate and resistance are strongly coupled. Compared with the flow rate distribution computed from the global simulation with constant resistance, the coupled solution predicts a flow rate with only slight changes. The high flow rate during diastole increases the stenosis pressure drop and resistance. In turn, this increased resistance of the stenosis slightly reduces the flow rate computed in the lumped parameter simulation.