Physiological flow analysis in significant human coronary artery stenoses

To evaluate the local hemodynamics in flow limiting coronary lesions, computational hemodynamics was applied to a group of patients previously reported by Wilson et al. (1988) with representative pre-angioplasty stenosis geometry (minimal lesion size d(m)=0.95 mm; 68% mean diameter stenosis) and with measured values of coronary flow reserve (CFR) in the abnormal range (2.3+/-0.1). The computations were at mean flow rates (Q) of 50, 75 and 100 ml/min (the limit of our converged calculations). Computed mean pressure drops Deltap were approximately 9 mmHg for basal flow (50 ml/min), approximately 27 mmHg for elevated flow (100 ml/min) and increased to an extrapolated value of approximately 34 mmHg for hyperemic flow (115 ml/min), which led to a distal mean coronary pressure p(rh) of approximately 55 mmHg, a level known to cause ischemia in the subendocardium (Brown et al., 1984), and consistent with the occurrence of angina in the patients. Relatively high levels of wall shear stress were computed in the narrow throat region and ranged from about 600 to 1500 dyn/cm(2), with periodic (phase shifted) peak systolic values of about 3500 dyn/cm(2). In the distal vessel, the interaction between the separated shear layer wave, convected downstream by the core flow, and the wall shear layer flow, led to the formation of vortical flow cells along the distal vessel wall during the systolic phase where Reynolds numbers Re(e)(t) were higher. During the phasic vortical mode observed at both basal and elevated mean flow rates, wide variations in distal wall shear stress occurred, distal transmural pressures were depressed below throat levels, and pressure recovery was larger farther along the distal vessel. Along the constriction (convergent) and throat segments of the lesion the pulsatile flow field was principally quasi-steady before flow separation occurred. The flow regimes were complex in the narrow mean flow Reynolds number range Re(e)=100-230 and a frequency parameter of alphae=2.25. The shear layer flow disturbances diminished in strength due to viscous damping along the distal vessel at these relatively low values of Re(e), typical of flow through diseased epicardial coronary vessels. The distal hyperemic flow field was likely to be in an early stage of turbulent flow development during the peak systolic phase.