Open Access

Flow and wall shear stress in end-to-side and side-to-side anastomosis of venous coronary artery bypass grafts

  • Thomas Frauenfelder1Email author,
  • Evangelos Boutsianis2,
  • Thomas Schertler1,
  • Lars Husmann1,
  • Sebastian Leschka1,
  • Dimos Poulikakos2,
  • Borut Marincek1 and
  • Hatem Alkadhi1
BioMedical Engineering OnLine20076:35

DOI: 10.1186/1475-925X-6-35

Received: 02 April 2007

Accepted: 26 September 2007

Published: 26 September 2007

Abstract

Purpose

Coronary artery bypass graft (CABG) surgery represents the standard treatment of advanced coronary artery disease. Two major types of anastomosis exist to connect the graft to the coronary artery, i.e., by using an end-to-side or a side-to-side anastomosis. There is still controversy because of the differences in the patency rates of the two types of anastomosis. The purpose of this paper is to non-invasively quantify hemodynamic parameters, such as mass flow and wall shear stress (WSS), in end-to-side and side-to-side anastomoses of patients with CABG using computational fluid dynamics (CFD).

Methods

One patient with saphenous CABG and end-to-side anastomosis and one patient with saphenous CABG and side-to-side anastomosis underwent 16-detector row computed tomography (CT). Geometric models of coronary arteries and bypasses were reconstructed for CFD analysis. Blood flow was considered pulsatile, laminar, incompressible and Newtonian. Peri-anastomotic mass flow and WSS were quantified and flow patterns visualized.

Results

CFD analysis based on in-vivo CT coronary angiography data was feasible in both patients. For both types of CABG, flow patterns were characterized by a retrograde flow into the native coronary artery. WSS variations were found in both anastomoses types, with highest WSS values at the heel and lowest WSS values at the floor of the end-to-side anastomosis. In contrast, the highest WSS values of the side-to-side anastomosis configuration were found in stenotic vessel segments and not in the close vicinity of the anastomosis. Flow stagnation zones were found in end-to-side but not in side-to-side anastomosis, the latter also demonstrating a smoother stream division throughout the cardiac cycle.

Conclusion

CFD analysis of venous CABG based on in-vivo CT datasets in patients was feasible producing qualitative and quantitative information on mass flow and WSS. Differences were found between the two types of anastomosis warranting further systematic application of the presented methodology on multiple patient datasets.

Background

Medical investigations

Coronary artery bypass graft (CABG) surgery represents the standard treatment of advanced coronary artery disease (CAD). Since the pioneering work of Favaloro [1], various grafts and different surgical techniques have been investigated. Regarding anastomoses, two major types exist to connect the graft to the coronary artery, i.e. by using an end-to-side or a side-to-side anastomosis. Although the latter has been reported to provide some advantages over individual grafting with end-to-side anastomosis [24], the sequential grafting technique has often been criticized because of differences in the patency rates of the two types of anastomosis [2, 5, 6]. The long-term clinical outcome after myocardial revascularization is dependent on graft patency. In venous CABG, however, accelerated atherosclerosis has been repetitively reported [7]. During the first year after CABG surgery up to 15% of venous grafts occlude, between 1 and 6 years the graft attrition rate is 1% to 2% per year, and between 6 to 10 years it is 4% per year. By 10 years after surgery only 60% of vein grafts are patent and only 50% of patent vein grafts are free of significant stenosis [8].

Neointimal hyperplasia, defined as the accumulation of smooth muscle cells and extracellular matrix in the intimal compartment, is the major disease process in venous CABG in the first year after surgery and sets the foundations for later development of graft atheroma [7]. In support of this proposal, the American Heart Association Council on Arteriosclerosis has defined the localized areas of adaptive neointimal hyperplasia as atherosclerosis prone regions [9]. During the first month after CABG surgery vein graft attrition results from thrombotic occlusion [10], while later on the dominant process is atherosclerotic obstruction occurring on a foundation of neointimal hyperplasia [11]. Although the risk factors predisposing to vein graft atherosclerosis are broadly similar to those recognized for native CAD, the pathogenic effects of these risk factors are amplified by inherent deficiencies of the vein as a conduit when transposed into the coronary arterial circulation [7].

In attempts to prevent acute and late graft occlusion, much effort has been invested in identifying the etiology of anastomotic neointimal hyperplasia and the plausibility of its prevention. Hypotheses related to this subject include the concept of compliance mismatch between graft and host artery [12], high frequency flow and wall shear stress (WSS) [13] as well as abnormal flow dynamics at the distal anastomosis [14]. Sottiurai et al. [15] have compared in an animal model the development of neointimal hyperplasia in end-to-side versus side-to-side anastomoses. Their study revealed that neointimal hyperplasia was present at the heel, toe, and floor of the end-to-side but not in the side-to-side anastomoses. Since compliance mismatch was not at issue in their model by using autogenous femorofemoral bypass grafts, the geometry of the distal anastomosis was attributed to be the causal factor. The configuration of end-to-side anastomosis is not a common occurrence in the primate cardiovascular system except for a patent "ductus arteriosus Botalli" [16].

Hemodynamic patterns at distal CABG anastomoses are thought to exhibit flow separation, recirculation and moving stagnation zones. Such flow features correspond typically to low time averaged WSS, long residence times, shear oscillation and increased spatial WSS gradients. These are the main hemodynamic features that have been connected with atherogenesis and intimal hyperplasia [17, 18]. In a recent review Kassab and Navia [19] put forward the homeostasis hypothesis as the underlying mechanism of graft failure. It is postulated that there is a mechanical homeostatic state for each blood vessel. Growth and remodeling takes place to abate perturbations from this state, such as when a vein segment is used in CABG. The remodeling procedure can fail or lead to a non-physiological state when the incurred perturbations are sufficiently large. Hemodynamic flow assessment in coronary arteries and CABGs is commonly performed with intravascular Doppler ultrasound [20]. However, this approach is invasive and the introduction of an ultrasound-catheter leads to flow disturbances thus interfering with the measurements [21]. Lately, phase-contrast magnetic resonance imaging (PC-MRI) [22] provides a non-invasive means of in-vivo three dimensional instantaneous velocity measurements on selected arterial cross-sections at a sufficient resolution.

CFD investigations

An alternative to invasive or non-invasive flow measurements is the simulation of blood flow by using computational fluid dynamics (CFD) and/or ex-vivo experimental flow setups. Numerical studies on physiological coronary flow as well as on CABG flow have been published extensively during the last 15 years. Steinman [23] presented a transient two dimensional flow simulation within a rigid 45 degree end-to-side anastomosis model. His results indicated elevated instantaneous WSS values at the toe and heel of the graft-host junction as well as along the host artery bed. In 1998 Ethier [24] simulated three dimensional unsteady flows in a symmetrical 45 degree end-to-side configuration with solid walls. It was shown that the perianastomotic WSS distributions are influenced by a complex interplay between secondary flow effects and the unsteadiness of the graft inflow waveform. An authoritative numerical description of the flow separation patterns of 45 degree junctions was given in [25] under steady flow conditions for a wide range of inlet Reynolds (Re) numbers, 250 – 1650. This study utilized an adaptive mesh refinement scheme to ensure grid independency.

Research on the hemodynamics of anastomoses has largely focused on determining the influence of the geometric characteristics of end-to-side configurations on WSS distributions. Fei [26] presented a series of steady flow simulations for a set of anastomosis angles ranging from 20 to 70 degrees. More recently Freshwater [27] utilized pulsatile flow boundary waveforms typical of a left internal mammary artery graft to the left anterior descending coronary artery scenario. They showed that higher anastomotic angles result in higher WSS magnitude values and flow oscillation around the toe and along the bed of the host branch. In all of the studies referenced so far, the centerlines of the graft and the host vessel lie within the same plane therefore forming planar configurations. Sherwin [28] investigated the influence of out of plane geometry in stiff end-to-side anastomosis models under steady and Papaharilaou [29] under transient flow conditions. Such configurations are closer to reality and result in non-symmetric flow fields. In general, the anastomosis bed was affected the most with reduced peak WSS values in comparison to similar planar models, while the mean oscillatory WSS magnitude was also shown to decrease. In most of the aforementioned investigations the native coronary or host vessel was fully occluded. Several researchers have addressed the issue of prograde and/or retrograde flow through the native coronary artery, as in [30, 31] under steady flow conditions and in [32] for pulsatile flow. In order to assess the potential role of activated platelets in the pathogenesis of intimal hyperplasia and/or thrombogenesis, numerical models have tried to quantify shear exposure and near wall residence times [33, 34]. Finally, Sankaranarayanan [35, 36] presented two studies, in a planar and a non-planar three dimensional CABG model respectively, that included the proximal anastomosis to the aorta.

It is worth noticing that most of the published CFD studies on anastomosis hemodynamics have been evaluated versus equivalent flow experiments by utilizing a variety of flow visualization, velocity and WSS measurement techniques such as in [3740]. Combined investigations utilizing CFD and PC-MRI on coronary and bypass flows are of particular interest because of the potential of MRI in providing in-vivo velocity measurements as well as anatomical information [41]. The vast majority of the above mentioned investigations focused on the hemodynamics of end-to-side anastomosis configurations. Bonert [42] examined the hemodynamics of side-to-side CABG anastomosis in solid idealized geometric models. They concluded that the parallel form of side-to-side anastomosis, as opposed to a non-parallel one, is better suited to maintain graft patency. A comparison between side-to-side and end-to-side anastomosis showed increased hemodynamic risk in the former approach due to the presence of larger low WSS areas.

Purpose of the present investigation

Advances in the CFD software and its corresponding hardware along with medical imaging are leading to the generation of increasingly reliable computational models. We can retrace a similar development into the fluid dynamics investigations of physiological coronary flow. Intracoronary flow can now be addressed in anatomically accurate configurations ranging from stiff multi-branched models [43] to moving and compliant arterial sections [44]. The present study investigated hemodynamic features of CABG anastomoses using patient-specific data. Pulsatile blood flow in venous CABGs was simulated using CFD on solid patient-specific geometric models based on in-vivo CT coronary angiography datasets. We assessed differences in WSS, mass flow, and flow patterns in venous conduits with end-to-side compared to side-to-side anastomosis.

Methods

Patients

Two male patients, 64 and 55-years old respectively, underwent cardiac CT. Patient 1 suffered from CAD with recurrent episodes of angina pectoris and a recent myocardial infarction. Invasive coronary angiography revealed a significant stenosis of the proximal right coronary artery (RCA) and non-significant stenoses of the middle left anterior descending (LAD) and distal left circumflex artery (LCX). Subsequently, the patient underwent saphenous CABG surgery with an end-to-side anastomosis onto the distal RCA. Patient 2 suffered from dyspnoea and instable angina pectoris. Invasive coronary angiography showed serial significant stenoses of the proximal LAD and a significant stenosis of the distal LCX. Saphenous CABG surgery was performed with a sequential side-to-side anastomosis onto the middle LAD and an end-to-side anastomosis onto the distal LCX.

CT data acquisition

CT was performed 10 and 15 days after surgery, respectively, on a 16-detector row scanner (Sensation 16, Siemens Medical Solutions, Forchheim, Germany) using the following parameters: detector collimation 16 × 0.75 mm, gantry rotation time 0.37 sec, pitch 0.38, tube potential 120 kV, tube current time product 400 mAs. A bolus of 150 ml iodinated contrast material (iodixanol, Visipaque 320, 320 mg/ml, GE Healthcare, Buckinghamshire, UK) followed by 30 ml saline solution was continuously injected into a right antecubital vein via a 18-gauge catheter at a flow rate of 5 ml/sec. Bolus tracking was performed with a region of interest in the ascending aorta and image acquisition was automatically started 5 sec after signal attenuation reached a threshold of 140 HU. Synchronized to the electrocardiogram (ECG), CT data sets were retrospectively reconstructed throughout the cardiac cycle in 5% steps of the R-R interval with a slice thickness of 1 mm and an increment of 0.5 mm using a medium soft-tissue convolution kernel (B30f). The adaptive cardio volume approach was used for image reconstruction and ECG-pulsing was applied to reduce radiation exposure. The reconstruction phase providing best image quality with the lowest degree of motion artifacts was determined by two readers in consensus and was used for further post-processing. The local ethics committee approved the study protocol and written informed consent was obtained from both patients.

Geometric reconstruction

Axial CT images were digitally processed to extract geometrical contours representing the coronary arteries and the CABGs. The lumen of all coronary arteries and grafts of the two patients were semi-automatically segmented using a commercially available software package (Amira 3.1, TGS, Belgium). In regions of reduced arterial opacification, segmentation was manually complemented. The outflows (i.e., the end of the branches) and inflows (i.e., the ostia) of the vessels were separately marked to allow the imposition of boundary conditions. As a next step, an unstructured surface mesh of triangles was generated covering the segmented volume using the marching cube algorithm. Manual smoothing and low-pass spatial filtering was then applied to further reduce fine-scale surface irregularities. The final model depicted the real three dimensional (3D) geometry of the coronary arteries and bypass grafts (Figure 1). Four computational models, two for each case, were subsequently built with 750.000 – 900.000 tetrahedral cells. The finest meshes represented a spatial resolution of about 0.15 mm and the number of elements per cross-section ranged from 150 to 200 depending on the vessel diameter.
Figure 1

Models. Anatomical depiction of the two models including the aortocoronary bypasses (arrows), the coronary arteries (arrowheads), the anastomoses (asterisks) and their relation to the heart and the aorta.

Model assumptions and boundary conditions

The flow for the simulation was considered transient, 3D, incompressible, and laminar. Corresponding to standard values from the literature, blood was assumed Newtonian with a viscosity of 0.0037 Pa· sec and a density of 1060 kg/m3. The walls were modeled as solid and stiff and a zero-velocity boundary condition was adopted at the fluid-solid interface, corresponding to a no-slip condition. In contrast to other flow simulations, we chose not to elongate the inlet part for the coronary arteries to allow for a full flow development. The instantaneous velocity at the different inlets (coronary ostia and proximal bypass anastomoses) was based on standard data [45, 46] reflecting the physiologically pulsatile, biphasic blood velocity from the ascending aorta into the coronary arteries (Figure 2). We adopted a spatially uniform profile for these boundary cross-sections. The maximum inlet Re was calculated at 1230 (bypass inlet) lasting for a very short period of time. The mean Re ranged between 380 (right coronary artey inlet) and 570 (left coronary artery inlet). The duration of the cardiac cycle was normalized to 1 second for both patients. The corresponding Womersley (Wo) number was 3.5 (bypass inlet), 2.1 (right coronary inlet) and 2.45 (left coronary inlet). The stress-free boundary condition (zero normal and tangential stresses), which arises naturally from the application of the finite element method, was imposed on the velocity field at all outlets to facilitate a common comparison frame between these two highly different coronary circulations by removing the influence of their downstream impedances. In this way, the observed flow differences are mainly the result of the adopted anastomosis type.
Figure 2

Input functions. Flow velocities throughout one cardiac cycle used as input functions at the ostium of the left main artery (LMA), right coronary artery (RCA), and the coronary artery bypass graft.

Computational fluid dynamics

The finite element software FIDAP (Version 8.6.2, Fluent Corp., Darmstadt, Germany) was used to perform the CFD simulations by solving the Navier-Stokes equations with linear basis functions. Calculated flow variables were flow velocities and pressure following a segregated solution approach. A convergence criterion of four orders of magnitude was adopted for the residuals in velocity and pressure. The instantaneous flow field was acquired at 100 steps per cardiac cycle with a constant time step using backward Euler implicit time integration. All the results presented herein belong to the third computational cardiac cycle to allow ample time for the attenuation of the effects of the initial conditions. The sensitivity of the numerical results on the underlying grid was examined under steady flow conditions. The time averaged values of the inflow velocities, shown in Figure 2, were used for this purpose by adopting a spatially uniform profile. The resulting differences in the calculated velocity values averaged below 5% of the maximum prescribed inflow velocity. This was judged sufficient for the requirements of the present feasibility investigation. The presented CFD results for both patients came from the finest available meshes. The Fieldview software (Version 11.0, Intelligent Light, Lyndhurst, NJ) was used for visualization of flow patterns, quantification of WSS and volumetric fluxes at selected sites.

Results

Mean heart rate during CT scanning was 47 ± 9 bpm for patient 1 and 51 ± 6 bpm for patient 2. The percent phase providing best image quality was found between 50 and 60% of the RR-interval for the RCA and between 60 and 70% for the left main artery (LMA), LAD, and LCX.

End-to-side anastomosis

Mass flow

The mean volumetric flow through the right coronary artery was 3.07 ml/sec. The mean volumetric flow through the CABG was 1.81 ml/sec (Figure 3; cut 5). Quantifications of mass flow through the proximal and distal part of the right coronary artery and the bypass close to the end-to-side anastomosis are demonstrated in Figure 3. The distribution of pulsatile flow could be determined by using the flow curves at each time-step of the cardiac cycle. There was a significant backflow into the proximal segment of the RCA (Figure 3; cut 2), reaching a maximum of 0.48 ml/sec. Due to this high backflow, the next proximal branch of the RCA (Figure 3; cut 4) was filled by blood coming from the bypass for the largest part of the cardiac cycle. The amount of backflow depends on the coronary mass flow and the bypass mass flow, which are characterized by different mass flow curves at their inlets, the angle between the bypass and the RCA, and the degree of the stenosis in the bypassed artery. The retrograde flow into the RCA reached a distance of 2.8 cm from the end-to-side anastomosis (Figure 4), leading to an area of stagnation between the two more proximal branches (Figure 3; cut 1 and cut 2).
Figure 3

Flux in end-to-side anastomosis. Integrated volumetric flux (ml/sec) at indicated cross-sections (cuts) in the RCA and the end-to-side anastomosis.

Figure 4

Velocity and shear stress in end-to-side anastomosis. Velocity-coded streamlines (m/sec) and wall shear stress (Pa) in an end-to-side anastomosis at three different time steps of the cardiac cycle.

Flow patterns and WSS

The WSS characterizes the tangential fluid forces that act on the vessel wall. The changes of WSS throughout the cardiac cycle showed a correlation with flow velocities, with WSS forces being high when blood flow was fast. Near the end-to-side anastomosis the maximum WSS spatial variation was approximately 1.5 Pa (Figures 4 and 5). The WSS ranged from 0.01 Pa at minimum inflow to approximately 2.0 Pa at maximum inflow. In the perianastomotic region, the time averaged WSS was 0.36 Pa. An area of high WSS was found during systole at the heel of the anastomosis. The lowest WSS was found in the RCA at the location where blood from the CABG hit the wall and was confronted by the flow coming from the native coronary artery. This formed a stagnation point at the impact site around which the bypass stream splits, as depicted by the flow streamlines in Figure 4. It is also shown here that backflow into the RCA persisted during most of the cardiac cycle. The highest WSS values were found during mid-systole (timestep 20) and were located at the heel of the anastomosis and at the side of a nearby clip that was placed during surgery (Figure 5). At the opposite side of the anastomosis the WSS values were lower, due to the presence of the stagnation zone.
Figure 5

Shear stress distribution in end-to-side anastomosis. Color-coded wall shear stress (Pa) and corresponding values plotted against the distance (mm) from the green to the red point at timestep 20 on three different paths upon the walls of the end-to-side anastomosis.

Side-to-side anastomosis

Mass flow

The mean volumetric flow through the left coronary artery was 3.35 ml/sec. The mean volumetric flow through the CABG was 1.51 ml/sec. Quantifications of volumetric flow through the proximal and distal segment of the LAD and the proximal and distal bypass close to the side-to-side anastomosis are demonstrated in Figure 6. The distribution of pulsatile flow could be determined by using the flow curves at each time-step of the cardiac cycle. There was a small backflow into the proximal part of the left coronary artery (Figure 6; cut 3) lasting almost 40% of the cardiac cycle and reaching a maximum of 0.25 ml/sec. The maximum volumetric flux into the distal part of the artery (Figure 6; cut 4) was 1.65 ml/sec or 46.2 % of the entire flow. Throughout a cardiac cycle, 53% of the mass flow into the bypass remained in it and reached the next anastomosis (Figure 6; cut 2). In diastole additional blood circulated from the proximal into the distal part of the coronary artery (Figure 6; cut 3), partially compensating for the decrease in the amount of blood that came from the bypass (Figure 6; cut 4).
Figure 6

Flux in side-to-side anastomosis. Integrated volumetric flux (ml/sec) at indicated cross-sections (cuts) in the LAD and the side-to-side anastomosis.

Flow patterns and WSS

Comparing side-to-side to end-to-side anastomosis, the mean WSS values were lower in the side-to-side anastomosis with a time-averaged value of 0.29 Pa in the perianastomotic region (Figure 7 and 8). The WSS values ranged from 0.02 Pa up to approximately 2.1 Pa at mid-systole. There were some areas of elevated WSS close to the anastomosis, but the highest WSS values were found in the atherosclerotic distal coronary artery (Figure 8). The elevated WSS in the distal coronary artery was mainly caused by increased bypass mass flow during systole. Furthermore, there was a high WSS region in the proximal part of the bypass due to a mild stenosis (Figure 7 and 8). In contrast, in the end-to-side anastomosis highest WSS values were located in the distal bypass graft close to the anastomosis. Lowest WSS values were found in the native coronary artery proximal to the side-to-side anastomosis. The maximum WSS spatial variation was approximately 1.0 Pa, being smaller than the WSS variation in end-to-side anastomosis. The streamline patterns presented in Figure 7 indicate absence of stagnation zones in the perianastomotic region of side-to-side anastomosis and a smoother stream division throughout the cardiac cycle as compared to the end-to-side anastomosis.
Figure 7

Velocity and shear-stress in side-to-side anastomosis. Velocity-coded streamlines (m/sec) and wall shear stress (Pa) in a side-to-side anastomosis at three different time steps of the cardiac cycle.

Figure 8

Shear stress distribution in side-to-side anastomosis. Color-coded wall shear stress (Pa) and corresponding values plotted against the distance (mm) from the green to the red point at timestep 20 on four different paths upon the walls of the side-to-side anastomosis.

Discussion

The results from our study indicate that considerable variations of blood flow exist in the distal perianastomotic region of venous CABG that may coincide with areas of neointimal hyperplasia. We found WSS variations in both types of anastomoses, with highest WSS values at the heel and lowest WSS values at the floor of the end-to-side anastomosis case. In contrast, high WSS in the side-to-side anastomosis configuration was only found in stenotic vessel segments but not in the close vicinity of the anastomosis. Across both types of anastomoses, elevated WSS values did coincide with vessel stenosis, being either caused by atherosclerotic wall changes or by a surgical clip.

CFD provided detailed information on instantaneous mass flow. In the case of side-to-side anastomosis, this information was important to predict the mass flow reaching the next anastomosis. In addition, the amount of blood running retrograde in the coronary artery and eventually reaching a more proximal branch can be quantitatively measured. The numerical results in our patients revealed considerable retrograde flow into the host coronary artery in both types of anastomoses, a finding that has been previously described in experimental models [38]. On the other hand stagnation zones were only found in end-to-side as compared to side-to-side anastomosis. The latter configuration showed smoother stream divisions of mass flow throughout the cardiac cycle. Finally, it was clearly demonstrated that the instantaneous blood influx through the bypasses maximizes during systole. Conversely, under physiological coronary flow conditions the maximum flux into the system occurs during diastole. This is a major difference that leads to a more balanced total instantaneous flow in the cases reported herein. Possible pathophysiological implications of this fact warrant further research.

When investigating hemodynamic features of blood flow, numerical simulation has become an important tool. This study presents realistic, patient-specific models based on CT angiography datasets of coronary and CABG anatomy. This approach circumvents the need for idealization of the main geometrical features. Most of the CFD studies that were mentioned in the introduction had to rely on some sort of geometrical regularization by introducing pre-specified parameters such as the CABG angle, the host to graft diameter ratio and the planarity or non-planarity of the entire configuration. There are merits in the latter methodology. There is the ability to parametrically study the hemodynamic effects of these strictly defined morphological parameters. Furthermore the generation of the underlying surface and volumetric grids is much easier and can lead to numerical results of greater accuracy. On the other hand, a patient specific model is much more difficult to construct and to ensure sufficient numerical accuracy. It is a closer match to reality but at the same time it is harder to parameterize the results and extract generalized conclusions.

CT angiography currently provides the most accurate representation of real anatomy. In contrast to most computational studies, our geometric models did not only include a small part of a coronary artery or CABG bypass but entailed the coronary tree and bypass grafts from their origin to the distal anastomoses. This became possible with the improvement of multi-detector row CT scanner technology with fast gantry rotation times enabling accurate and reliable depiction of coronary and bypass graft anatomy in a short imaging time [47, 48]. With the advent of 64-slice CT [49] and dual-source CT [50] further improvements with regard to temporal and spatial resolution have been made allowing an even more precise depiction of anatomy that will further improve the accuracy of numerical flow simulations.

The following study limitations have to be acknowledged. Firstly, 16-detector row CT does not offer currently the best imaging for coronary arteries and CABG as compared to newer scanner technology. In addition, geometric data for CFD was obtained from a single reconstruction time-point in early to mid-diastole and vessel contours and diameters might be different at other reconstruction time-points. Secondly, the assumptions concerning inflow and outflow may not be valid under pathologic conditions and thus leading to exaggerated mass flow in the bypass grafts. The localization of the bypasses and native coronary inlets is extracted directly from the CT scans. It is expected that a realistic inflow velocity profile would not exhibit developed flow characteristics, as in Poisseuille or Womersley theories. The main discrepancy lies on the significant presence of in plane secondary flow patterns [51] due to the motion of the coronary sinuses. Typical PC-MRI resolution does not allow for the quantification of these in-plane velocity distributions. Conversely, PC-MRI or Doppler Ultrasound can be used to quantify the blood flux discharge ratios between the first few branches of the coronary and bypass configurations. It is our intention to enhance our investigation protocol with such in-vivo measurements in order to improve the imposed outflow boundary conditions.

Moreover, both the arterial and bypass walls were simplified as being stiff. Heart movement and changing pressure at the outer wall of the vessels due to the myocardium contraction were not included in the present calculations. The importance of the large scale coronary motion induced by the beating heart is hard to quantify in CABG configurations. Although the coronary arteries move considerably throughout the cardiac cycle, it has been shown that pulsatility is the main characterizing factor of WSS distributions [52]. Additionally, the bypass grafts are intentionally fixed with surgical clips and thus experience only minor displacements near the distal anastomosis sites. Their interaction is expected to affect the pressure and WSS distributions near the anastomosis. Results might differ in elastic models of coronary arterial and CABG walls that can also take into account the compliance mismatch between the host and graft sections [53, 54]. Recently, Ramaswamy et al. [44] performed a fluid dynamics analysis in a diseased section of a human coronary artery by taking into account both the large scale motion as well as compliance. This was accomplished by utilizing bi-plane angiography combined with intravascular ultrasound images. Their results showed that compliance caused substantial differences in the circumferential WSS distributions. The profound change in the distal coronary lumen dimensions during systole has major consequences as far as the outflow boundary conditions are concerned, by significantly increasing the flow resistance. Additionally, the concurrent increase of arterial pressure leads to the dilatation of the proximal coronary sections. Coronary hemodynamics and the corresponding WSS distributions will now depend on the interaction between pulsatile inflow, which maximizes during diastole, and the variation of the luminal diameter, which maximizes during systole.

Finally, the assumption that blood behaves as a Newtonian liquid must be examined further. Although it is assumed that blood exhibits shear thinning behavior in vessels of smaller diameter than the coronary arteries, there is a multitude of patient dependent parameters like the hematocrit and the intake of drugs that can alter blood viscosity. Numerical simulations of non-Newtonian flow in a two dimensional end-to-side anastomosis under pulsatile flow conditions have shown only minor effects on WSS distributions [55]. However such results are bound to depend on the particulars of the constitutive equations utilized for determining blood's viscosity.

Conclusion

Our study has investigated the influence of patient specific geometry of end-to-side and side-to-side anastomosis on perianastomotic hemodynamics to identify geometrically driven flow features that might increase the propensity of venous graft failure. CFD analysis allowed us to differentiate hemodynamics in the two types of anastomoses and indicated significant spatial WSS variations especially in the end-to-side anastomosis configuration and an absence of stagnation areas in side-to-side anastomosis. The CFD simulations were performed in only two patients, a fact that limits the significance of our findings, which might differ in patients with different geometry. The numerical method applied herein may provide the basis for future prospective investigations correlating hemodynamic features with neointimal hyperplasia after CABG surgery.

Declarations

Acknowledgements

Supported by the National Center of Competence in Research, Computer Aided and Image Guided Medical Interventions (NCCR-CO ME) of the Swiss National Science Foundation.

Authors’ Affiliations

(1)
Institute of Diagnostic Radiology, University Hospital Zurich
(2)
Laboratory of Thermodynamics in Emerging Technologies, ETH Zurich

References

  1. Favaloro RG: Saphenous vein autograft replacement of severe segmental coronary artery occlusion: operative technique. Ann Thorac Surg 1968, 5: 334–339.View ArticleGoogle Scholar
  2. Grondin CM, Campeau L, Lesperance J, Enjalbert M, Bourassa MG: Comparison of late changes in internal mammary artery and saphenous vein grafts in two consecutive series of patients 10 years after operation. Circulation 1984, 70: I208–212.Google Scholar
  3. Meeter K, Veldkamp R, Tijssen JG, van Herwerden LL, Bos E: Clinical outcome of single versus sequential grafts in coronary bypass operations at ten years' follow-up. J Thorac Cardiovasc Surg 1991, 101: 1076–1081.Google Scholar
  4. O'Neill MJ Jr, Wolf PD, O'Neill TK, Montesano RM, Waldhausen JA: A rationale for the use of sequential coronary artery bypass grafts. J Thorac Cardiovasc Surg 1981, 81: 686–690.Google Scholar
  5. Christenson JT, Simonet F, Schmuziger M: Sequential vein bypass grafting: tactics and long-term results. Cardiovasc Surg 1998, 6: 389–397. 10.1016/S0967-2109(98)00024-6View ArticleGoogle Scholar
  6. Kieser TM, FitzGibbon GM, Keon WJ: Sequential coronary bypass grafts. Long-term follow-up. J Thorac Cardiovasc Surg 1986, 91: 767–772.Google Scholar
  7. Motwani JG, Topol EJ: Aortocoronary saphenous vein graft disease: pathogenesis, predisposition, and prevention. Circulation 1998, 97: 916–931.View ArticleGoogle Scholar
  8. Campeau L, Enjalbert M, Lesperance J, Bourassa MG, Kwiterovich P Jr, Wacholder S, Sniderman A: The relation of risk factors to the development of atherosclerosis in saphenous-vein bypass grafts and the progression of disease in the native circulation. A study 10 years after aortocoronary bypass surgery. N Engl J Med 1984, 311: 1329–1332.View ArticleGoogle Scholar
  9. Stary HC, Blankenhorn DH, Chandler AB, Glagov S, Insull W Jr, Richardson M, Rosenfeld ME, Schaffer SA, Schwartz CJ, Wagner WD, et al.: A definition of the intima of human arteries and of its atherosclerosis-prone regions. A report from the Committee on Vascular Lesions of the Council on Arteriosclerosis, American Heart Association. Circulation 1992, 85: 391–405.View ArticleGoogle Scholar
  10. Bourassa MG: Fate of venous grafts: the past, the present and the future. J Am Coll Cardiol 1991, 17: 1081–1083.View ArticleGoogle Scholar
  11. Dilley RJ, McGeachie JK, Tennant M: Vein to artery grafts: a morphological and histochemical study of the histogenesis of intimal hyperplasia. Aust N Z J Surg 1992, 62: 297–303.View ArticleGoogle Scholar
  12. Abbott WM, Megerman J: Does compliance mismatch alone cause neointimal hyperplasia? J Vasc Surg 1989, 9: 507.View ArticleGoogle Scholar
  13. Imparato AM, Bracco A, Kim GE, Zeff R: Intimal and neointimal fibrous proliferation causing failure of arterial reconstructions. Surgery 1972, 72: 1007–1017.Google Scholar
  14. Faulkner SL, Fisher RD, Conkle DM, Page DL, Bender HW Jr: Effect of blood flow rate on subendothelial proliferation in venous autografts used as arterial substitutes. Circulation 1975, 52: I163–172.Google Scholar
  15. Sottiurai VS, Sue SL, Feinberg EL 2nd, Bringaze WL, Tran AT, Batson RC: Distal anastomotic intimal hyperplasia: biogenesis and etiology. Eur J Vasc Surg 1988, 2: 245–256. 10.1016/S0950-821X(88)80034-3View ArticleGoogle Scholar
  16. Sottiurai VS: Distal Anastomotic Intimal Hyperplasia: Histocytomorphology, Pathophysiology, Etiology, and Prevention. International Journal Of Angiology 1999, 8: 1–10.View ArticleGoogle Scholar
  17. Giddens DP, Zarins CK, Glagov S: The Role of Fluid-Mechanics in the Localization and Detection of Atherosclerosis. Journal of Biomechanical Engineering-Transactions of the Asme 1993,115(4B):588–594.View ArticleGoogle Scholar
  18. Nerem RM: Vascular Fluid-Mechanics, the Arterial-Wall, and Atherosclerosis. Journal of Biomechanical Engineering-Transactions of the Asme 1992,114(3):274–282.View ArticleGoogle Scholar
  19. Kassab GS, Navia JA: Biomechanical considerations in the design of graft: The homeostasis hypothesis. Annual Review of Biomedical Engineering 2006, 8: 499–535. 10.1146/annurev.bioeng.8.010506.105023View ArticleGoogle Scholar
  20. Pijls NHJ, De Bruyne B, Bech GJW, Liistro F, Heyndrickx GR, Bonnier HJRM, Koolen JJ: Coronary pressure measurement to assess the hemodynamic significance of serial stenoses within one coronary artery – Validation in humans. Circulation 2000, 102: 2371–2377.View ArticleGoogle Scholar
  21. Krams R, Wentzel JJ, Cespedes I, Vinke R, Carlier S, van der Steen AF, Lancee CT, Slager CJ: Effect of catheter placement on 3-D velocity profiles in curved tubes resembling the human coronary system. Ultrasound Med Biol 1999, 25: 803–810. 10.1016/S0301-5629(99)00018-6View ArticleGoogle Scholar
  22. Lotz J, Meier C, Leppert A, Galanski M: Cardiovascular flow measurement with phase-contrast MR imaging: basic facts and implementation. Radiographics 2002, 22: 651–671.View ArticleGoogle Scholar
  23. Steinman DA, Vinh B, Ethier CR, Ojha M, Cobbold RS, Johnston KW: A numerical simulation of flow in a two-dimensional end-to-side anastomosis model. J Biomech Eng 1993, 115: 112–118.View ArticleGoogle Scholar
  24. Ethier CR, Steinman DA, Zhang X, Karpik SR, Ojha M: Flow waveform effects on end-to-side anastomotic flow patterns. J Biomech 1998, 31: 609–617. 10.1016/S0021-9290(98)00059-1View ArticleGoogle Scholar
  25. Ethier CR, Prakash S, Steinman DA, Leask RL, Couch GG, Ojha M: Steady flow separation patterns in a 45 degree junction. J Fluid Mech 2000, 411: 1–38. 10.1017/S0022112099008022View ArticleMATHGoogle Scholar
  26. Fei DY, Thomas JD, Rittgers SE: The effect of angle and flow rate upon hemodynamics in distal vascular graft anastomoses: a numerical model study. J Biomech Eng 1994, 116: 331–336.View ArticleGoogle Scholar
  27. Freshwater IJ, Morsi YS, Lai T: The effect of angle on wall shear stresses in a LIMA to LAD anastomosis: numerical modelling of pulsatile flow. Proc Inst Mech Eng [H] 2006, 220: 743–757.View ArticleGoogle Scholar
  28. Sherwin SJ, Shah O, Doorly DJ, Peiro J, Papaharilaou Y, Watkins N, Caro CG, Dumoulin CL: The influence of out-of-plane geometry on the flow within a distal end-to-side anastomosis. J Biomech Eng 2000, 122: 86–95. 10.1115/1.429630View ArticleGoogle Scholar
  29. Papaharilaou Y, Doorly DJ, Sherwin SJ: The influence of out-of-plane geometry on pulsatile flow within a distal end-to-side anastomosis. Journal of Biomechanics 2002, 35: 1225–1239. 10.1016/S0021-9290(02)00072-6View ArticleGoogle Scholar
  30. Bertolotti C, Deplano V: Three-dimensional numerical simulations of flow through a stenosed coronary bypass. J Biomech 2000, 33: 1011–1022. 10.1016/S0021-9290(00)00012-9View ArticleGoogle Scholar
  31. Kute SM, Vorp DA: The effect of proximal artery flow on the hemodynamics at the distal anastomosis of a vascular bypass graft: computational study. J Biomech Eng 2001, 123: 277–283. 10.1115/1.1374203View ArticleGoogle Scholar
  32. Bertolotti C, Deplano V, Fuseri J, Dupouy P: Numerical and experimental models of post-operative realistic flows in stenosed coronary bypasses. J Biomech 2001, 34: 1049–1064. 10.1016/S0021-9290(01)00027-6View ArticleGoogle Scholar
  33. Longest PW, Kleinstreuer C: Numerical simulation of wall shear stress conditions and platelet localization in realistic end-to-side arterial anastomoses. J Biomech Eng 2003, 125: 671–681. 10.1115/1.1613298View ArticleGoogle Scholar
  34. Sherwin SJ, Doorly DJ, Franke P, Peiro J: Unsteady near wall residence times and shear exposure in model distal arterial bypass grafts. Biorheology 2002, 39: 365–371.Google Scholar
  35. Sankaranarayanan M, Chua LP, Ghista DN, Tan YS: Computational model of blood flow in the aorto-coronary bypass graft. Biomed Eng Online 2005, 4: 14. 10.1186/1475-925X-4-14View ArticleGoogle Scholar
  36. Sankaranarayanan M, Ghista DN, Poh CL, Seng TY, Kassab GS: Analysis of blood flow in an out-of-plane CABG model. Am J Physiol Heart Circ Physiol 2006, 291: H283–295. 10.1152/ajpheart.01347.2005View ArticleGoogle Scholar
  37. Heise M, Schmidt S, Kruger U, Ruckert R, Rosler S, Neuhaus P, Settmacher U: Flow pattern and shear stress distribution of distal end-to-side anastomoses. A comparison of the instantaneous velocity fields obtained by particle image velocimetry. J Biomech 2004, 37: 1043–1051. 10.1016/j.jbiomech.2003.11.030View ArticleGoogle Scholar
  38. Hughes PE, How TV: Effects of geometry and flow division on flow structures in models of the distal end-to-side anastomosis. J Biomech 1996, 29: 855–872. 10.1016/0021-9290(95)00168-9View ArticleGoogle Scholar
  39. Li XM, Rittgers SE: Hemodynamic factors at the distal end-to-side anastomosis of a bypass graft with different POS:DOS flow ratios. J Biomech Eng 2001, 123: 270–276. 10.1115/1.1372323View ArticleGoogle Scholar
  40. Ojha M: Spatial and temporal variations of wall shear stress within an end-to-side arterial anastomosis model. J Biomech 1993, 26: 1377–1388. 10.1016/0021-9290(93)90089-WView ArticleGoogle Scholar
  41. Papaharilaou Y, Doorly DJ, Sherwin SJ, Peiro J, Griffith C, Cheshire N, Zervas V, Anderson J, Sanghera B, Watkins N, Caro CG: Combined MR imaging and numerical simulation of flow in realistic arterial bypass graft models. Biorheology 2002, 39: 525–531.Google Scholar
  42. Bonert M, Myers JG, Fremes S, Williams J, Ethier CR: A numerical study of blood flow in coronary artery bypass graft side-to-side anastomoses. Annals of Biomedical Engineering 2002, 30: 599–611. 10.1114/1.1481052View ArticleGoogle Scholar
  43. Boutsianis E, Dave H, Frauenfelder T, Poulikakos D, Wildermuth S, Turina M, Ventikos Y, Zund G: Computational simulation of intracoronary flow based on real coronary geometry. Eur J Cardiothorac Surg 2004, 26: 248–256. 10.1016/j.ejcts.2004.02.041View ArticleGoogle Scholar
  44. Ramaswamy SD, Vigmostad SC, Wahle A, Lai YG, Olszewski ME, Braddy KC, Brennan TM, Rossen JD, Sonka M, Chandran KB: Fluid dynamic analysis in a human left anterior descending coronary artery with arterial motion. Ann Biomed Eng 2004, 32: 1628–1641. 10.1007/s10439-004-7816-3View ArticleGoogle Scholar
  45. Canver CC, Dame NA: Ultrasonic Assessment of Internal Thoracic Artery Graft Flow in the Revascularized Heart. Annals of Thoracic Surgery 1994, 58: 135–138.View ArticleGoogle Scholar
  46. Drost C: On/Off Pump Graft Patency Assessment. United States of America: Transonic Systems, Inc; 2002.Google Scholar
  47. Nieman K, Pattynama PM, Rensing BJ, Van Geuns RJ, De Feyter PJ: Evaluation of patients after coronary artery bypass surgery: CT angiographic assessment of grafts and coronary arteries. Radiology 2003, 229: 749–756. 10.1148/radiol.2293020856View ArticleGoogle Scholar
  48. Pache G, Saueressig U, Frydrychowicz A, Foell D, Ghanem N, Kotter E, Geibel-Zehender A, Bode C, Langer M, Bley T: Initial experience with 64-slice cardiac CT: non-invasive visualization of coronary artery bypass grafts. European Heart Journal 2006, 27: 976–980. 10.1093/eurheartj/ehi824View ArticleGoogle Scholar
  49. Flohr TG, Stierstorfer K, Ulzheimer S, Bruder H, Primak AN, McCollough CH: Image reconstruction and image quality evaluation for a 64-slice CT scanner with z-flying focal spot. Med Phys 2005, 32: 2536–2547. 10.1118/1.1949787View ArticleGoogle Scholar
  50. Scheffel H, Alkadhi H, Plass A, Vachenauer R, Desbiolles L, Gaemperli O, Schepis T, Frauenfelder T, Schertler T, Husmann L, et al.: Accuracy of dual-source CT coronary angiography: First experience in a high pre-test probability population without heart rate control. Eur Radiol 2006, 16: 2739–2747. 10.1007/s00330-006-0474-0View ArticleGoogle Scholar
  51. Suo J: Investigation of Blood Flow Patterns and Hemodynamics in the human ascending aorta and major trunks of right and left coronary arteries using magnetic resonance imaging and computational fluid dynamics. In PhD Thesis. Georgia Institute of Technology, Biomedical Engineering; 2005.Google Scholar
  52. Zeng DH, Ding ZH, Friedman MH, Ethier CR: Effects of cardiac motion on right coronary artery hemodynamics. Annals of Biomedical Engineering 2003, 31: 420–429. 10.1114/1.1560631View ArticleGoogle Scholar
  53. Hofer M, Rappitsch G, Perktold K, Trubel W, Schima H: Numerical study of wall mechanics and fluid dynamics in end-to-side anastomoses and correlation to intimal hyperplasia. J Biomech 1996, 29: 1297–1308. 10.1016/0021-9290(96)00036-XView ArticleGoogle Scholar
  54. Leuprecht A, Perktold K, Prosi M, Berk T, Trubel W, Schima H: Numerical study of hemodynamics and wall mechanics in distal end-to-side anastomoses of bypass grafts. J Biomech 2002, 35: 225–236. 10.1016/S0021-9290(01)00194-4View ArticleGoogle Scholar
  55. Ballyk PD, Steinman DA, Ethier CR: Simulation of non-Newtonian blood flow in an end-to-side anastomosis. Biorheology 1994, 31: 565–586.Google Scholar

Copyright

© Frauenfelder et al; licensee BioMed Central Ltd. 2007

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Advertisement