Computational simulation of intracoronary flow based on real coronary geometry

To assess the feasibility of computationally simulating intracoronary blood flow based on real coronary artery geometry and to graphically depict various mechanical characteristics of this flow.

Explanted fresh pig hearts were fixed using a continuous perfusion of 4% formaldehyde at physiological pressures. Omnipaque dye added to lead rubber solution was titrated to an optimum proportion of 1:25, to cast the coronary arterial tree. The heart was stabilized in a phantom model so as to suspend the base and the apex without causing external deformation. High resolution computerized tomography scans of this model were utilized to reconstruct the three-dimensional coronary artery geometry, which in turn was used to generate several volumetric tetrahedral meshes of sufficient density needed for numerical accuracy. The transient equations of momentum and mass conservation were numerically solved by employing methods of computational fluid dynamics under realistic pulsatile inflow boundary conditions.

The simulations have yielded graphic distributions of intracoronary flow stream lines, static pressure drop, wall shear stress, bifurcation mass flow ratios and velocity profiles. The variability of these quantities within the cardiac cycle has been investigated at a temporal resolution of 1/100th of a second and a spatial resolution of about 10 microm. The areas of amplified variations in wall shear stress, mostly evident in the neighborhoods of arterial branching, seem to correlate well with clinically observed increased atherogenesis. The intracoronary flow lines showed stasis and extreme vorticity during the phase of minimum coronary flow in contrast to streamlined undisturbed flow during the phase of maximum flow.

Computational tools of this kind along with a state-of-the-art multislice computerized tomography or magnetic resonance-based non-invasive coronary imaging, could enable realistic, repetitive, non-invasive and multidimensional quantifications of the effects of stenosis on distal hemodynamics, and thus help in precise surgical/interventional planning. It could also add insights into coronary and bypass graft atherogenesis.