Computational model of blood flow in the aorto-coronary bypass graft

Abstract

BackgroundCoronary artery bypass grafting surgery is an effective treatment modality for patients with severe coronary artery disease. The conduits used during the surgery include both the arterial and venous conduits. Long- term graft patency rate for the internal mammary arterial graft is superior, but the same is not true for the saphenous vein grafts. At 10 years, more than 50% of the vein grafts would have occluded and many of them are diseased. Why do the saphenous vein grafts fail the test of time? Many causes have been proposed for saphenous graft failure. Some are non-modifiable and the rest are modifiable. Non-modifiable causes include different histological structure of the vein compared to artery, size disparity between coronary artery and saphenous vein. However, researches are more interested in the modifiable causes, such as graft flow dynamics and wall shear stress distribution at the anastomotic sites. Formation of intimal hyperplasia at the anastomotic junction has been implicated as the root cause of long- term graft failure.Many researchers have analyzed the complex flow patterns in the distal sapheno-coronary anastomotic region, using various simulated model in an attempt to explain the site of preferential intimal hyperplasia based on the flow disturbances and differential wall stress distribution. In this paper, the geometrical bypass models (aorto-left coronary bypass graft model and aorto-right coronary bypass graft model) are based on real-life situations. In our models, the dimensions of the aorta, saphenous vein and the coronary artery simulate the actual dimensions at surgery. Both the proximal and distal anastomoses are considered at the same time, and we also take into the consideration the cross-sectional shape change of the venous conduit from circular to elliptical. Contrary to previous works, we have carried out computational fluid dynamics (CFD) study in the entire aorta-graft-perfused artery domain. The results reported here focus on (i) the complex flow patterns both at the proximal and distal anastomotic sites, and (ii) the wall shear stress distribution, which is an important factor that contributes to graft patency.MethodsThe three-dimensional coronary bypass models of the aorto-right coronary bypass and the aorto-left coronary bypass systems are constructed using computational fluid-dynamics software (Fluent 6.0.1). To have a better understanding of the flow dynamics at specific time instants of the cardiac cycle, quasi-steady flow simulations are performed, using a finite-volume approach. The data input to the models are the physiological measurements of flow-rates at (i) the aortic entrance, (ii) the ascending aorta, (iii) the left coronary artery, and (iv) the right coronary artery.ResultsThe flow field and the wall shear stress are calculated throughout the cycle, but reported in this paper at two different instants of the cardiac cycle, one at the onset of ejection and the other during mid-diastole for both the right and left aorto-coronary bypass graft models. Plots...

c-Jun Regulates Shear- and Injury-inducible Egr-1 Expression, Vein Graft Stenosis after Autologous End-to-Side Transplantation in Rabbits, and Intimal Hyperplasia in Human Saphenous Veins

“…The process of intimal thickening as a prelude to bypass graft stenosis may be initiated by changes in fluid biomechanical forces sensed by the cells in the graft (18). For example, conduits undergo an acute static-to-shear transition after the reperfu- Ϫ ␣-SM actin ϩ cells reside at the inner surface of human saphenous vein conduits making contact with flowing blood (magnification, 200ϫ).…”

Physiological insights of recent clinical diagnostic and therapeutic technologies for cardiovascular diseases

“…To evaluate shear stress more accurately, it is necessary to solve the Navier–Stokes equations for the velocity vector ( u ) [ 67 ]. When the blood is assumed to be incompressible ( ρ = constant), the Navier–Stokes equations are described as below: where ρ is the blood density, ν is the kinematic viscosity ( ν = μ / ρ ), and p is the pressure.…”

Simulations of Magnetohemodynamics in Stenosed Arteries in Diabetic or Anemic Models

“…The potential for obtaining information helpful in therapeutic decision making was demonstrated by analyzing selected hemodynamic parameters such as blood flow pathlines, wall shear stresses, dynamics pressures, blood flow velocities, and flow particle residence times. Sankaranarayanan et al [ 7 ] transient model of aorta-coronary bypass graft showed that maximum perfusion of the occluded artery occurs during middiastole, and the maximum wall shear stress variation is observed around the distal anastomotic region. Pulsatile flow through a stenosed artery using Casson model for blood was investigated by Sankar and Lee [ 8 ].…”