The influence of out-of-plane geometry on pulsatile flow within a distal end-to-side anastomosis

Multi-Objective Optimisation of a Novel Bypass Graft with a Spiral Ridge

The low long-term patency of bypass grafts is a major concern for cardiovascular treatments. Unfavourable haemodynamic conditions in the proximity of distal anastomosis are closely related to thrombus creation and lumen lesions. Modern graft designs address this unfavourable haemodynamic environment with the introduction of a helical component in the flow field, either by means of out-of-plane helicity graft geometry or a spiral ridge. While the latter has been found to lack in performance when compared to the out-of-plane helicity designs, recent findings support the idea that the existing spiral ridge grafts can be further improved in performance through optimising relevant design parameters. In the current study, robust multi-objective optimisation techniques are implemented, covering a wide range of possible designs coupled with proven and well validated computational fluid dynamics (CFD) algorithms. It is shown that the final set of suggested design parameters could significantly improve haemodynamic performance and therefore could be used to enhance the design of spiral ridge bypass grafts.

Distinct Temporal Pattern of the Prediction of Lumen Remodeling of Lower Extremity Vein Bypass Grafts by Initial Local Hemodynamics

We predicted human lower extremity vein bypass graft remodeling by hemodynamics. Computed tomography and duplex ultrasound scans of 55 patients were performed at 1 week and 1, 6, and 12 months post-implantation to obtain wall shear stress (WSS) and oscillatory shear index (OSI) at 1-mm intervals via computational fluid dynamics simulations. Graft remodeling was quantified by computed tomography-measured lumen diameter changes in the early (1 week-1 month), intermediate (1-6 months), and late (6-12 months) periods. Linear mixed-effect models were constructed to examine the overall relationship between remodeling and initial hemodynamics using the average data of all cross sections within the same graft. A significant association of graft remodeling with WSS (p < 0.001) and time (p = 0.001) was found; however, the effect size decreased with time (every 2.7 dyne/cm² increase of WSS was associated with a 0.39, 0.35, 0.002 mm diameter increase in the three periods, respectively). The association of remodeling with OSI was significant only in the intermediate period (every 0.1 increase of OSI was associated with a 0.25 mm lumen diameter decrease, p = 0.004). Therefore, the association of graft lumen remodeling with local hemodynamics has a distinct temporal pattern; WSS and OSI are predictive of remodeling only in certain postoperative periods.

Secondary flow in bifurcations - Important effects of curvature, bifurcation angle and stents

Coronary bifurcations have complex flow patterns including secondary flow zones and helical flow, which directly affect pathophysiological mechanisms such as the development of atherosclerosis. The objective of this study was to generate insights into the effects of curvature, bifurcation angle and the presence of stents on flow patterns and resulting haemodynamics in coronary left main bifurcations. The blood flow and associated metrics were modelled in both idealised and patient-specific bifurcations with varying curvature and bifurcation angles with and without stents, resulting in a total of 128 geometries considered. The results showed that larger curvature of bifurcating vessels has a significant influence on secondary flow, especially with distance to the bifurcation region, causing a skew, spin and asymmetry of Dean vortices, an increase in helical flow intensity with symmetry loss, and a decrease in adversely low time-average wall shear stress (TAWSS). Generally, asymmetric flow patterns coincided with adversely low TAWSS regions. In identical stented geometries, the presence of the stents induced local recirculation immediately adjacent to the stent struts, thus generating adversely low TAWSS in these areas, with some effect on the overall secondary flow. Overall, the effect of stents outweighed the effect of curvature and BA. This new knowledge contributes to a better understanding of the joint effects of curvature, bifurcation angle, and stents on flow patterns and haemodynamics in coronary bifurcations.

Flow analysis of aortic dissection: comparison of inflow boundary conditions for computational models based on 4D PCMRI and Doppler ultrasound

Computational hemodynamics quantifying the flow environment is an important tool in understanding aortic dissection. In this study, various inflow boundaries were applied on a patient-specific model and compared to the individualized velocimetry. The results indicated that the computations generally overestimated the flow volume and underestimated the wall shear stress. By quantifying the accuracy of the simulation results, two inflow settings were suggested. One was individualized, the PCMRI-extracted 4D flow information, and the other was averaged by healthy data, the ultrasound-extracted averaged flow waveform with parabolic velocity profile. This study might contribute to improving the precise computation of aortic dissection hemodynamics.

Retrograde branched extension limb assembling stent of pararenal abdominal aortic aneurysm: A longitudinal hemodynamic analysis for stent graft migration

PURPOSE

Pararenal abdominal aortic aneurysms (PRAAAs) are a life-threatening disease, and hemodynamic analysis may provide greater insight into the effectiveness and long-term outcomes of endovascular aneurysm repair (EVAR). However, the lack of patient-specific boundary conditions on the periphery compromises the accuracy. Windkessel (WK) boundary conditions coupled to hemodynamic follow-up models of a PRAAA patient, aims to provide insights into the link between hemodynamics and poor prognosis.

METHOD

One PRAAA patient underwent EVAR and reintervention after one branch of stent-graft (SG) had migrated. Totally five computational follow-up models were studied. Patient-specific flow data acquired via ultrasound were used to define the boundary conditions in the ascending aorta and the following three branches. Coupled zero-dimensional WK models representing the distal vasculature were used to define the outlet boundary conditions under the abdomen.

RESULTS

Flow divisions of the main SG branches were 40.7% and 24.7%, respectively. Time-averaged wall shear stress and oscillatory shear index (OSI) increased at the junction connected the SG branch and the stent leading to the right common iliac artery (RCIA) where the stent migrated. The OSI and relative residence time (RRT) value in superior mesenteric artery increased notably after the migration, the RRT continuously increased following the reintervention.

CONCLUSION

Unbalanced flow, resulting in locally high-speed flow, high WSS and OSI might significantly affect stent stability. Results suggest that diameters and interconnection design of stents in complex cases should take the flow division into consideration and computational simulations might be considered as a tool for intervention protocol design.

Hemodynamic parameters that may predict false-lumen growth in type-B aortic dissection after endovascular repair: A preliminary study on long-term multiple follow-ups

Thoracic endovascular aortic repair (TEVAR) is commonly applied in type-B aortic dissection. For patients with dissection affects descending aorta and extends downward to involve abdominal aorta and possibly iliac arteries, false lumen (FL) expansion might occur post-TEVAR. Predictions of dissection development may assist in medical decision on re-intervention or surgery. In this study, two patients are selected with similar morphological features at initial presentation but with different long-term FL development post-TEVAR (stable and enlarged FL). Patient-specific models are established for each of the follow-ups. Flow boundaries and computational validations are obtained from Doppler ultrasound velocimetry. By analyzing the hemodynamic parameters, the false-to-true luminal pressure difference (PDiff) and particle relative residence time (RRT) are found related to FL remodeling. It is found that (i) the position of the first FL flow entry is the watershed of negative-and-positive PDiff and, in long-term follow-ups, and the position of largest PDiff is consistent with that of the greatest increase of FL width; (ii) high RRT occurs at the FL proximal tip and similar magnitude of RRT is found in both stable and enlarged cases; (iii) comparing to the RRT at 7days post-TEVAR, an increase of RRT afterwards in short-term is found in the stable case while a slight decrease of this parameter is found in the enlarged case, indicating that the variation of RRT in short-term post-TEVAR might be potential to predict long-term FL remodeling.

Optimisation of a Novel Spiral-Inducing Bypass Graft Using Computational Fluid Dynamics

Graft failure is currently a major concern for medical practitioners in treating Peripheral Vascular Disease (PVD) and Coronary Artery Disease (CAD). It is now widely accepted that unfavourable haemodynamic conditions play an essential role in the formation and development of intimal hyperplasia, which is the main cause of graft failure. This paper uses Computational Fluid Dynamics (CFD) to conduct a parametric study to enhance the design and performance of a novel prosthetic graft, which utilises internal ridge(s) to induce spiral flow. This design is primarily based on the identification of the blood flow as spiral in the whole arterial system and is believed to improve the graft longevity and patency rates at distal graft anastomoses. Four different design parameters were assessed in this work and the trailing edge orientation of the ridge was identified as the most important parameter to induce physiological swirling flow, while the height of the ridge also significantly contributed to the enhanced performance of this type of graft. Building on these conclusions, an enhanced configuration of spiral graft is proposed and compared against conventional and spiral grafts to reaffirm its potential benefits.

Numerical Assessment of Novel Helical/Spiral Grafts with Improved Hemodynamics for Distal Graft Anastomoses

In the present work, numerical simulations were conducted for a typical end-to-side distal graft anastomosis to assess the effects of inducing secondary flow, which is believed to remove unfavourable flow environment. Simulations were carried out for four models, generated based on two main features of 'out-of-plane helicity' and 'spiral ridge' in the grafts as well as their combination. Following a qualitative comparison against in vitro data, various mean flow and hemodynamic parameters were compared and the results showed that helicity is significantly more effective in inducing swirling flow in comparison to a spiral ridge, while their combination could be even more effective. In addition, the induced swirling flow was generally found to be increasing the wall shear stress and reducing the flow stagnation and particle residence time within the anastomotic region and the host artery, which may be beneficial to the graft longevity and patency rates. Finally, a parametric study on the spiral ridge geometrical features was conducted, which showed that the ridge height and the number of spiral ridges have significant effects on inducing swirling flow, and revealed the potential of improving the efficiency of such designs.

Haemodynamics and Flow Modification Stents for Peripheral Arterial Disease: A Review

Endovascular stents are widely used for the treatment of peripheral arterial disease (PAD). However, the development of in-stent restenosis and downstream PAD progression remain a challenge. Stent revascularisation of PAD causes arterial trauma and introduces abnormal haemodynamics, which initiate complicated biological processes detrimental to the arterial wall. The interaction between stent struts and arterial cells in contact, and the blood flow field created in a stented region, are highly affected by stent design. Spiral flow is known as a normal physiologic characteristic of arterial circulation and is believed to prevent the development of flow disturbances. This secondary flow motion is lost in atheromatous disease and its re-introduction after endovascular treatment of PAD has been suggested as a method to induce stabilised and coherent haemodynamics. Stent designs able to generate spiral flow may support endothelial function and therefore increase patency rates. This review is focused on secondary flow phenomena in arteries and the development of flow modification stent technologies for the treatment of PAD.

The Tissue-Engineered Vascular Graft-Past, Present, and Future

Cardiovascular disease is the leading cause of death worldwide, with this trend predicted to continue for the foreseeable future. Common disorders are associated with the stenosis or occlusion of blood vessels. The preferred treatment for the long-term revascularization of occluded vessels is surgery utilizing vascular grafts, such as coronary artery bypass grafting and peripheral artery bypass grafting. Currently, autologous vessels such as the saphenous vein and internal thoracic artery represent the gold standard grafts for small-diameter vessels (<6 mm), outperforming synthetic alternatives. However, these vessels are of limited availability, require invasive harvest, and are often unsuitable for use. To address this, the development of a tissue-engineered vascular graft (TEVG) has been rigorously pursued. This article reviews the current state of the art of TEVGs. The various approaches being explored to generate TEVGs are described, including scaffold-based methods (using synthetic and natural polymers), the use of decellularized natural matrices, and tissue self-assembly processes, with the results of various in vivo studies, including clinical trials, highlighted. A discussion of the key areas for further investigation, including graft cell source, mechanical properties, hemodynamics, integration, and assessment in animal models, is then presented.

Effect of swirling blood flow on vortex formation at post-stenosis

Various clinical observations reported that swirling blood flow is a normal physiological flow pattern in various vasculatures. The swirling flow has beneficial effects on blood circulation through the blood vessels. It enhances oxygen transfer and reduces low-density lipoprotein concentration in the blood vessel by enhancing cross-plane mixing of the blood. However, the fluid-dynamic roles of the swirling flow are not yet fully understood. In this study, inhibition of material deposition at the post-stenosis region by the swirling flow was observed. To reveal the underlying fluid-dynamic characteristics, pathline flow visualization and time-resolved particle image velocimetry measurements were conducted. Results showed that the swirling inlet flow increased the development of vortices at near wall region of the post-stenosis, which can suppress further development of stenosis by enhancing transport and mixing of the blood flow. The fluid-dynamic characteristics obtained in this study would be useful for improving hemodynamic characteristics of vascular grafts and stents in which the stenosis frequently occurred. Moreover, the time-resolved particle image velocimetry measurement technique and vortex identification method employed in this study would be useful for investigating the fluid-dynamic effects of the swirling flow on various vascular environments.

Numerical investigation of blood flow in a deformable coronary bifurcation and non-planar branch

INTRODUCTION

Among cardiovascular diseases, arterials stenosis is recognized more commonly than the others. Hemodynamic characteristics of blood play a key role in the incidence of stenosis. This paper numerically investigates the pulsatile blood flow in a coronary bifurcation with a non-planar branch. To create a more realistic analysis, the wall is assumed to be compliant. Furthermore, the flow is considered to be three-dimensional, incompressible, and laminar.

METHODS

The effects of non-Newtonian blood, compliant walls and different angles of bifurcation on hemodynamic characteristics of flow were evaluated. Shear thinning of blood was simulated with the Carreau-Yasuda model. The current research was mainly focused on the flow characteristics in bifurcations since atherosclerosis occurs mostly in bifurcations. Moreover, as the areas with low shear stresses are prone to stenosis, these areas were identified.

RESULTS

Our findings indicated that the compliant model of the wall, bifurcation's angle, and other physical properties of flow have an impact on hemodynamics of blood flow. Lower wall shear stress was observed in the compliant wall than that in the rigid wall. The outer wall of bifurcation in all models had lower wall shear stress. In bifurcations with larger angles, wall shear stress was higher in outer walls, and lower in inner walls.

CONCLUSION

The non-Newtonian blood vessels and different angles of bifurcation on hemodynamic characteristics of flow evaluation confirmed a lower wall shear stress in the compliant wall than that in the rigid wall, while the wall shear stress was higher in outer walls but lower in inner walls in the bifurcation regions with larger angles.

A computational fluid dynamics study on hemodynamics for different locations of the distal anastomosis of a bypass nearby a collateral vessel in the femoropopliteal area

Revascularization of the femoropopliteal sector is often performed by the placement of a bypass. In this paper, we have studied the effects of hemodynamics on patency of the bypass for different positions of the distal anastomosis close to a collateral artery. Computational fluid dynamics (CFD) are used for this study. The cardiac cycle-averaged wall shear stress (WSS) and oscillation index (OSI) have been analyzed. Low WSS and high OSI may increase the risk of intimal hyperplasia (IH), which may reduce bypass patency. From the CFD simulations, spots of low WSS and high OSI are found within and near the entrance of the collateral artery, near the suture line, at the floor, toe, and heel. We regarded flow ratios of 20:80 and of 35:65. It is found that for the high flow ratio anastomosis located proximal to the collateral artery is clearly more advantageous. However for the low flow ratio anastomosis located distal to the collateral artery seems to be slightly more advantageous, the results are less conclusive. One of the studied flow geometries has been validated by in vitro experiments using a time resolved particle image velocimetry technique. Velocity fields from these experiments are in good agreement with the CFD results.

Fluid-dynamic optimal design of helical vascular graft for stenotic disturbed flow

Although a helical configuration of a prosthetic vascular graft appears to be clinically beneficial in suppressing thrombosis and intimal hyperplasia, an optimization of a helical design has yet to be achieved because of the lack of a detailed understanding on hemodynamic features in helical grafts and their fluid dynamic influences. In the present study, the swirling flow in a helical graft was hypothesized to have beneficial influences on a disturbed flow structure such as stenotic flow. The characteristics of swirling flows generated by helical tubes with various helical pitches and curvatures were investigated to prove the hypothesis. The fluid dynamic influences of these helical tubes on stenotic flow were quantitatively analysed by using a particle image velocimetry technique. Results showed that the swirling intensity and helicity of the swirling flow have a linear relation with a modified Germano number (Gn*) of the helical pipe. In addition, the swirling flow generated a beneficial flow structure at the stenosis by reducing the size of the recirculation flow under steady and pulsatile flow conditions. Therefore, the beneficial effects of a helical graft on the flow field can be estimated by using the magnitude of Gn*. Finally, an optimized helical design with a maximum Gn* was suggested for the future design of a vascular graft.

Coronary artery bypass grafting hemodynamics and anastomosis design: a biomedical engineering review

In this paper, coronary arterial bypass grafting hemodynamics and anastomosis designs are reviewed. The paper specifically addresses the biomechanical factors for enhancement of the patency of coronary artery bypass grafts (CABGs). Stenosis of distal anastomosis, caused by thrombosis and intimal hyperplasia (IH), is the major cause of failure of CABGs. Strong correlations have been established between the hemodynamics and vessel wall biomechanical factors and the initiation and development of IH and thrombus formation. Accordingly, several investigations have been conducted and numerous anastomotic geometries and devices have been designed to better regulate the blood flow fields and distribution of hemodynamic parameters and biomechanical factors at the distal anastomosis, in order to enhance the patency of CABGs. Enhancement of longevity and patency rate of CABGs can eliminate the need for re-operation and can significantly lower morbidity, and thereby reduces medical costs for patients suffering from coronary stenosis. This invited review focuses on various endeavors made thus far to design a patency-enhancing optimized anastomotic configuration for the distal junction of CABGs.

Computation in the rabbit aorta of a new metric - the transverse wall shear stress - to quantify the multidirectional character of disturbed blood flow

Spatial variation of the haemodynamic stresses acting on the arterial wall is commonly assumed to explain the focal development of atherosclerosis. Disturbed flow in particular is thought to play a key role. However, widely-used metrics developed to quantify its extent are unable to distinguish between uniaxial and multidirectional flows. We analysed pulsatile flow fields obtained in idealised and anatomically-realistic arterial geometries using computational fluid dynamics techniques, and in particular investigated the multidirectionality of the flow fields, capturing this aspect of near-wall blood flow with a new metric - the transverse wall shear stress (transWSS) - calculated as the time-average of wall shear stress components perpendicular to the mean flow direction. In the idealised branching geometry, multidirectional flow was observed downstream of the branch ostium, a region of flow stagnation, and to the sides of the ostium. The distribution of the transWSS was different from the pattern of traditional haemodynamic metrics and more dependent on the velocity waveform imposed at the branch outlet. In rabbit aortas, transWSS patterns were again different from patterns of traditional metrics. The near-branch pattern varied between intercostal ostia, as is the case for lesion distribution; for some branches there were striking resemblances to the age-dependent patterns of disease seen in rabbit and human aortas. The new metric may lead to improved understanding of atherogenesis.

Numerical analysis of non-Newtonian blood flow and wall shear stress in realistic single, double and triple aorto-coronary bypasses

Considering the fact that hemodynamics plays an important role in the patency and overall performance of implanted bypass grafts, this work presents a numerical investigation of pulsatile non-Newtonian blood flow in three different patient-specific aorto-coronary bypasses. The three bypass models are distinguished from each other by the number of distal side-to-side and end-to-side anastomoses and denoted as single, double and triple bypasses. The mathematical model in the form of time-dependent nonlinear system of incompressible Navier-Stokes equations is coupled with the Carreau-Yasuda model describing the shear-thinning property of human blood and numerically solved using the principle of the SIMPLE algorithm and cell-centred finite volume method formulated for hybrid unstructured tetrahedral grids. The numerical results computed for non-Newtonian and Newtonian blood flow in the three aorto-coronary bypasses are compared and analysed with emphasis placed on the distribution of cycle-averaged wall shear stress and oscillatory shear index. As shown in this study, the non-Newtonian blood flow in all of the considered bypass models does not significantly differ from the Newtonian one. Our observations further suggest that, especially in the case of sequential grafts, the resulting flow field and shear stimulation are strongly influenced by the diameter of the vessels involved in the bypassing.

Competitive flow in coronary bypass surgery: is it a problem?

PURPOSE OF REVIEW

Competition flow is a common finding in arterial grafting. This phenomenon can lead to graft occlusion.

RECENT FINDINGS

Internal thoracic arteries are the best equipped arterial conduit to withstand the competition flow thanks to their endothelial function. Radial as well as right gastroepiploic arteries support much less flow competition because of their different anatomy, histology and endothelial function than that of internal thoracic artery, leading to spasm and occlusion. Therefore, these two arterial conduits should be used only in case of critical lesion to avoid graft occlusion. Saphenous vein graft is the only conduit that is not significantly affected by flow competition, mainly because of its nonresistivity and common reimplantation in the aorta. Graft configuration is the second important factor influencing the equation between graft flow and native coronary flow. Therefore, composite grafting should be reserved for the case of severely stenotic coronary target, especially if multiple arterial sequential anastomoses have to be performed on the lateral-inferior wall of the heart. Finally, an accurate tool such as the fractional flow reserve to evaluate the stenosis severity should be the milestone of coronary surgery in order to decrease the rate of flow competition and improve arterial grafting functionality.

CONCLUSION

Competition flow plays a crucial role in arterial grafting functionality. Grafting strategy should address this by appropriate graft choice and configuration in order to avoid graft attrition.

Hemodynamic variations due to spiral blood flow through four patient-specific bifurcated stent graft configurations for the treatment of abdominal aortic aneurysms

Endovascular repair is now a recognised procedure for treating abdominal aortic aneurysms. However, post-operative complications such as stent graft migration and thrombus may still occur. To assess these complications numerically, the correct input boundary conditions, which include the full human aorta with associated branching, should be included. Four patient-specific computed tomography scanned bifurcated stent grafts (SGs) were modelled and attached onto a full human aorta, which included the ascending, aortic arch and descending aortas. Two of the SG geometries had a twisted leg configuration, while the other two had conventional nontwisted leg configurations. Computational fluid dynamics was completed for both geometries and the hemodynamics assessed. The complexity of the flow patterns and secondary flows were influenced by the inclusion of the full human aorta at the SG proximal section. During the decelerating phase significant recirculations occurred along the main body of all SG configurations. The inclusion of the full human aorta did not impact the velocity contours within the distal legs and there was no difference in drag forces with the SG containing the full human aorta and those without. A twisted leg configuration further promoted a spiral flow formation along its distal legs.

A computational exploration of helical arterio-venous graft designs

Although arterio-venous grafts (AVGs) are the second best option as long-term vascular access for hemodialysis, they suffer from complications caused by intimal hyperplasia, mainly located in vessel regions of low and oscillating wall shear stress. However, certain flow patterns in the bulk may reduce these unfavorable hemodynamic conditions. We therefore studied, with computational fluid dynamics (CFD), the impact of a helical AVG design on the occurrence of (un)favorable hemodynamic conditions at the venous anastomosis. Six CFD-models of an AVG in closed-loop configuration were constructed: one conventional straight graft, and five helical designed grafts with a pitch of 105 mm down to 35 mm. At the venous anastomosis, disturbed shear was assessed by quantifying the area with unfavorable conditions, and by analyzing averaged values in a case-specific patch. The bulk hemodynamics were assessed by analyzing the kinetic helicity in and the pressure drop over the graft. The most helical design scores best, being instrumental to suppress disturbed shear in the venous segment. There is, however, no trivial relationship between the number of helix turns of the graft and disturbed shear in the venous segment, when a realistic closed-loop AVG model is investigated. Bulk flow investigation showed a marked increase of helicity intensity in, and a moderate pressure drop over the AVG by introducing a lower pitch. At the venous anastomosis, unfavorable hemodynamic conditions can be reduced by introducing a helical design. However, due to the complex flow conditions, the optimal helical design for an AVG cannot be derived without studying case by case.

Effect of head posture on the healthy human carotid bifurcation hemodynamics

Head and neck postures may cause morphology changes to the geometry of the carotid bifurcation (CB) that alter the low and oscillating wall shear stress (WSS) regions previously reported as important in the development of atherosclerosis. Here the right and left CB were imaged by MRI in two healthy subjects in the neutral head posture with the subject in the supine position and in two other head postures with the subject in the prone position: (1) rightward rotation up to 80°, and (2) leftward rotation up to 80°. Image-based computational models were constructed to investigate the effect of posture on arterial geometry and local hemodynamics. The area exposure to unfavorable hemodynamics, based on thresholds set for oscillatory shear index (OSI), WSS and relative residence time, was used to quantify the hemodynamic impact on the wall. Torsion of the head was found to: (1) cause notable changes in the bifurcation and internal carotid artery angles and, in most cases, on cross-sectional area ratios for common, internal and external carotid artery, (2) change the spatial distribution of wall regions exposed to unfavorable hemodynamics, and (3) cause a marked change in the hemodynamic burden on the wall when the OSI was considered. These findings suggest that head posture may be associated with the genesis and development of atherosclerotic disease as well as complications in stenotic and stented vessels.

Numerical analysis of coronary artery bypass grafts: an over view

Arterial bypass grafts tend to fail after some years due to the development of intimal thickening (restenosis). Non-uniform hemodynamics following a bypass operation contributes to restenosis and bypass failure can occur due to the focal development of anastomotic intimal hyperplasia. Additionally, surgical injury aggravated by compliance mismatch between the graft and artery has been suggested as an initiating factor for progress of wall thickening along the suture line Vascular grafts that are small in diameter tend to occlude rapidly. Computational fluid dynamics (CFD) methods have been effectively used to simulate the physical and geometrical parameters characterizing the hemodynamics of various arteries and bypass configurations. The effects of such changes on the pressure and flow characteristics as well as the wall shear stress during a cardiac cycle can be simulated. Recently, utilization of fluid and structure interactions have been used to determine fluid flow parameters and structure forces including stress and strains relationships under steady and transient conditions. In parallel to this, experimental diagnostics techniques such as Laser Doppler Anemometry, Particle Image Velocimetry, Doppler Guide wire and Magnetic Resonance Imaging have been used to provide essential information and to validate the numerical results. Moreover, clinical imaging techniques such as magnetic resonance or computed tomography have assisted considerably in gaining a detailed patient-specific picture of the blood flow and structure dynamics. This paper gives a review of recent numerical investigations of various configurations of coronary artery bypass grafts (CABG). In addition, the paper ends with a summary of the findings and the future directions.

Effect of branchings on blood flow in the system of human coronary arteries

In this work, we investigate the behavior of the pulsatile blood flow in the system of human coronary arteries. Blood is modeled as an incompressible non-Newtonian fluid. The transient phenomena of blood flow through the coronary system are simulated by solving the three dimensional unsteady state Navier-Stokes equations and continuity equation. Distributions of velocity, pressure and wall shear stresses are determined in the system under pulsatile conditions on the boundaries. Effect of branching vessel on the flow problem is investigated. The numerical results show that blood pressure in the system with branching vessels of coronary arteries is lower than the one in the system with no branch. The magnitude of wall shear stresses rises at the bifurcation.

Disturbed flow in radial-cephalic arteriovenous fistulae for haemodialysis: low and oscillating shear stress locates the sites of stenosis

BACKGROUND

Despite recent clinical and technological advancements, the vascular access (VA) for haemodialysis still has significant early failure rates after arteriovenous fistula (AVF) creation. VA failure is mainly related to the haemodynamic conditions that trigger the phenomena of vascular wall disease such as intimal hyperplasia (IH) or atherosclerosis.

METHODS

We performed transient computational fluid dynamics simulations within idealized three-dimensional models of 'end-to-side' and 'end-to-end' radio-cephalic anastomosis, using non-Newtonian blood and previously measured flows and division ratio in subjects requiring primary access procedure as boundary conditions.

RESULTS

The numerical simulations allowed full characterization of blood flow inside the AVF and of patterns of haemodynamic shear stress, known to be the major determinant of vascular remodelling and disease. Wall shear stress was low and oscillating in zones where flow stagnation occurs on the artery floor and on the inner wall of the juxta-anastomotic vein.

CONCLUSIONS

Zones of low and oscillatory shear stress were located in the same sites where luminal reduction was documented in previous experimental studies on sites stenosis distribution in AVF. We conclude that even when exposed to high flow rates, there are spot regions along the AVF exposed to athero-prone shear stress that favour vessel stenosis by triggering IH.

Computer-Aided Patient-Specific Coronary Artery Graft Design Improvements Using CFD Coupled Shape Optimizer

This study aims to (i) demonstrate the efficacy of a new surgical planning framework for complex cardiovascular reconstructions, (ii) develop a computational fluid dynamics (CFD) coupled multi-dimensional shape optimization method to aid patient-specific coronary artery by-pass graft (CABG) design and, (iii) compare the hemodynamic efficiency of the sequential CABG, i.e., raising a daughter parallel branch from the parent CABG in patient-specific 3D settings. Hemodynamic efficiency of patient-specific complete revascularization scenarios for right coronary artery (RCA), left anterior descending artery (LAD), and left circumflex artery (LCX) bypasses were investigated in comparison to the stenosis condition. Multivariate 2D constraint optimization was applied on the left internal mammary artery (LIMA) graft, which was parameterized based on actual surgical settings extracted from 2D CT slices. The objective function was set to minimize the local variation of wall shear stress (WSS) and other hemodynamic indices (energy dissipation, flow deviation angle, average WSS, and vorticity) that correlate with performance of the graft and risk of re-stenosis at the anastomosis zone. Once the optimized 2D graft shape was obtained, it was translated to 3D using an in-house "sketch-based" interactive anatomical editing tool. The final graft design was evaluated using an experimentally validated second-order non-Newtonian CFD solver incorporating resistance based outlet boundary conditions. 3D patient-specific simulations for the healthy coronary anatomy produced realistic coronary flows. All revascularization techniques restored coronary perfusions to the healthy baseline. Multi-scale evaluation of the optimized LIMA graft enabled significant wall shear stress gradient (WSSG) relief (~34%). In comparison to original LIMA graft, sequential graft also lowered the WSSG by 15% proximal to LAD and diagonal bifurcation. The proposed sketch-based surgical planning paradigm evaluated the selected coronary bypass surgery procedures based on acute hemodynamic readjustments of aorta-CA flow. This methodology may provide a rational to aid surgical decision making in time-critical, patient-specific CA bypass operations before in vivo execution.

Vascular corrosion casting: analyzing wall shear stress in the portal vein and vascular abnormalities in portal hypertensive and cirrhotic rodents

Vascular corrosion casting is an established method of anatomical preparation that has recently been revived and has proven to be an excellent tool for detailed three-dimensional (3D) morphological examination of normal and pathological microcirculation. In addition, the geometry provided by vascular casts can be further used to calculate wall shear stress (WSS) in a vascular bed using computational techniques. In the first part of this study, the microvascular morphological changes associated with portal hypertension (PHT) and cirrhosis in vascular casts are described. The second part of this study consists of a quantitative analysis of the WSS in the portal vein in casts of different animal models of PHT and cirrhosis using computational fluid dynamics (CFD). Microvascular changes in the splanchnic, hepatic and pulmonary territory of portal hypertensive and cirrhotic mice are described in detail with stereomicroscopic examination and scanning electron microscopy. To our knowledge, our results are the first to report the vascular changes in the common bile duct ligation cirrhotic model. Calculating WSS using CFD methods is a feasible technique in PHT and cirrhosis, enabling the differentiation between different animal models. First, a dimensional analysis was performed, followed by a CFD calculation describing the spatial and temporal WSS distributions in the portal vein. WSS was significantly different between sham/cirrhotic/pure PHT animals with the highest values in the latter. Up till now, no techniques have been developed to quantify WSS in the portal vein in laboratory animals. This study showed for the first time that vascular casting has an important role not only in the morphological evaluation of animal models of PHT and cirrhosis, but also in defining the biological response of the portal vein wall to hemodynamic changes. CFD in 3D geometries can be used to describe the spatial and temporal variations in WSS in the portal vein and to better understand the forces affecting mechanotransduction in the endothelium.

Hemodynamics of an end-to-side anastomotic graft for a pulsatile pediatric ventricular assist device

Numerical simulations are performed to investigate the flow within the end-to-side proximal anastomosis of a pulsatile pediatric ventricular assist device (PVAD) to an aorta. The anastomotic model is constructed from a patient-specific pediatric aorta. The three great vessels originating from the aortic arch--brachiocephalic (innominate), left common carotid, and left subclavian arteries--are included. An implicit large eddy simulation method based on a finite volume approach is used to study the resulting turbulent flow. A resistance boundary condition is applied at each branch outlet to study flow splitting. The PVAD anastomosis is found to alter the aortic flow dramatically. More flow is diverted into the great vessels with the PVAD support. Turbulence is found in the jet impingement area at peak systole for 100% bypass, and a maximum principal normal Reynolds stress of 7081 dyn/cm(2) is estimated based on ten flow cycles. This may be high enough to cause hemolysis and platelet activation. Regions prone to intimal hyperplasia are identified by combining the time-averaged wall shear stress and oscillatory shear index. These regions are found to vary, depending on the percentage of the flow bypass.

Arterial geometry, flow pattern, wall shear and mass transport: potential physiological significance

We have studied numerically steady and unsteady flow in a straight and a helically stented common carotid artery, in order to model porcine experimental results that show reduced intimal hyperplasia (IH) in the helical case. The combination of flow pulsatility and three-dimensionality generates a sweeping motion of the Dean vortices, which overall reduced extremes of both oxygen flux to the vessel wall and wall shear stress (WSS). Since IH and atherosclerosis affect preferentially low WSS regions, these findings imply that vessel three-dimensionality and flow pulsatility can play important protective roles in respect of these diseases. The amplitude and frequency of the velocity waveform are important parameters of the system. Increase in amplitude increases WSS and oxygen flux to the vessel wall. Increase in frequency has a small effect; it increases WSS but has no effect on the oxygen flux to the vessel wall.

Effect of Reynolds number and flow division on patterns of haemodynamic wall shear stress near branch points in the descending thoracic aorta

Atherosclerotic lesions are non-uniformly distributed at arterial bends and branch sites, suggesting an important role for haemodynamic factors, particularly wall shear stress (WSS), in their development. The pattern of lesions at aortic branch sites depends on age and species. Using computational flow simulations in an idealized model of an intercostal artery emerging perpendicularly from the thoracic aorta, we studied the effects of Reynolds number and flow division under steady conditions. Patterns of flow and WSS were strikingly dependent on these haemodynamic parameters. With increasing Reynolds number, WSS, normalized by the fully developed aortic value, was lowered at the sides of the ostium and increased upstream and downstream of it. Increasing flow into the side branch exacerbated these patterns and gave rise to a reversing flow region downstream of the ostium. Incorporation of more realistic geometric features had only minor effects and patterns of mean WSS under pulsatile conditions were similar to the steady flow results. Aspects of the observed WSS patterns correlate with, and may explain, some but not all of the lesion patterns in human, rabbit and mouse aortas.

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

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.

A decoupled fluid structure approach for estimating wall stress in abdominal aortic aneurysms

Abdominal aortic aneurysm (AAA) is a localized dilatation of the aortic wall. The lack of an accurate AAA rupture risk index remains an important problem in the clinical management of the disease. To accurately estimate AAA rupture risk, detailed information on patient-specific wall stress distribution and aortic wall tissue yield stress is required. A complete fluid structure interaction (FSI) study is currently impractical and thus of limited clinical value. On the other hand, isolated static structural stress analysis based on a uniform wall loading is a widely used approach for AAA rupture risk estimation that, however, neglects the flow-induced wall stress variation. The aim of this study was to assess the merit of a decoupled fluid structure analysis of AAA wall stress. Anatomically correct, patient specific AAA wall models were created by 3D reconstruction of computed tomography images. Flow simulations were carried out with inflow and outflow boundary conditions obtained from patient extracted data. Static structural stress analysis was performed applying both a uniform pressure wall loading and a flow induced non-uniform pressure distribution obtained during early systolic deceleration. For the structural analysis, a hyperelastic arterial wall model and an elastic intraluminal thrombus model were assumed. The results of this study demonstrate that although the isolated static structural stress analysis approach captures the gross features of the stress distribution it underestimates the magnitude of the peak wall stress by as much as 12.5% compared to the proposed decoupled fluid structure approach. Furthermore, the decoupled approach provides potentially useful information on the nature of the aneurysmal sac flow.

Flow capacity of gastroepiploic artery versus vein grafts for intermediate coronary artery stenosis

BACKGROUND

Native flow competition is a significant factor affecting bypass graft patency. The objective of this study was to compare the effect of competitive flow on conduit flow dynamics in the gastroepiploic artery (GEA) and the saphenous vein graft (SVG).

METHODS

In 51 patents, 23 GEAs (in-situ grafts) and 28 SVGs (aortocoronary grafts) were examined using a Doppler-tipped guidewire during coronary angiography after coronary artery bypass. Graft flow volume at rest and maximum graft flow volume during hyperemia were calculated from graft diameter and average peak velocity at rest and maximum average peak velocity induced by papaverine hydrochloride injection. Grafts were classified according to the grade of native coronary artery stenosis; group S (14 GEAs and 16 SVGs) displayed over 75% stenosis and group M (9 GEAs and 12 SVGs) exhibited over 50% up to 75% stenosis.

RESULTS

In group S, no difference in flow volume was apparent between the GEA and the SVG at rest (36+/- 17 vs 42 +/- 16) and during hyperemia (78 +/- 30 vs 88 +/- 28). In group M, flow volume of the GEA was significantly lower than that of the SVG at rest (17 +/- 11 vs 38 +/- 12; p = 0.029) and during hyperemia (32 +/- 19 vs 94 +/- 46; p = 0.001).

CONCLUSIONS

These data suggest that in intermediate coronary stenosis, GEA flow is compromised by native flow competition, whereas the SVG flow dynamics is maintained. However, the GEA can provide comparable flow capacity to the SVG and will achieve good surgical results when target coronary artery selection is appropriate.

Effects of flat, parabolic and realistic steady flow inlet profiles on idealised and realistic stent graft fits through Abdominal Aortic Aneurysms (AAA)

Computational fluid dynamics modelling of bifurcated stent grafts positioned across Abdominal Aortic Aneursysms (AAAs) is influenced significantly by the boundary conditions and geometry. To assess these influencing factors, simulations using three different steady inlet profiles and two geometries were run. The steady velocity inlet profiles (flat, parabolic and realistic) were used as input boundary conditions with a constant pressure outlet for each model. The geometry of an out-of-plane realistic stent graft model, fitted to the internal geometry of an AAA (a realistic geometry generated from CT scans) was compared to an idealised geometry of an in-plane model. It was found that the realistic boundary condition when applied to the idealised model had a 5% difference in outlet flow rates when compared to flat and parabolic inlet profiles, while the realistic models had an 11-16% difference. The blood flowed in a parallel streamlined fashion for all the idealised models with slight recirculation for the realistic boundary condition. All realistic models had significant recirculation and flowed in a left-handed helix motion in the proximal end and left leg and in a right-handed helix motion in the right leg. The drag force was unidirectional in the idealised models, while it was found to act in all three directions for the realistic models with a 26% increase over the idealised models. It was concluded that geometry has the greatest influence on the outlet flow rates, flow patterns and drag forces. Input boundary condition cause variations in the skewness and recirculation of the flow throughout the models.

Swirling flow pattern in a non-planar model of an interposition vein cuff anastomosis

One of the main causes of long-term failure of ePTFE grafts is the development of anastomotic intimal hyperplasia which leads to graft thrombosis. Experimental studies with bypass grafts have shown an inverse relationship between mean wall shear stress and intimal hyperplasia. The geometry of the anastomosis has a strong influence on the flow pattern and wall shear stress distribution. The aim of this in vitro study was to investigate the influence of non-planarity in a model of a distal anastomosis with interposition vein cuff, an anastomosis configuration that is increasingly being used because of improved clinical results. Laser Doppler anemometer measurements were carried out in silicone rubber models of interposition vein cuff anastomoses with planar and non-planar inflow. The pulsatile flow waveforms were typical of those found in femoro-infrapopliteal bypass. Axial and radial velocities were measured in the proximal and distal outflow segments. As expected a symmetrical helical flow pattern (Dean flow) was evident in the planar model. The model with non-planar inflow, however, gave rise to swirling flow in both the distal and proximal artery outflow segments for during the systolic phase. In patients, the anastomosis is usually non-planar. Since the configuration depends in part upon the tunnelling of the graft, this may be altered to some extent. Non-planar anastomotic configurations induce a swirling flow pattern, which may normalise wall shear stress, thereby potentially reducing intimal hyperplasia.

3-D numerical simulation of blood flow through models of the human aorta

A Spiral Computerized Tomography (CT) scan of the aorta were obtained from a single subject and three model variations were examined. Computational fluid dynamics modeling of all three models showed variations in the velocity contours along the aortic arch with differences in the boundary layer growth and recirculation regions. Further down-stream, all three models showed very similar velocity profiles during maximum velocity with differences occurring in the decelerating part of the pulse. Flow patterns obtained from transient 3-D computational fluid dynamics are influenced by different reconstruction methods and the pulsatility of the flow. Caution is required when analyzing models based on CT scans.

On reducing abnormal hemodynamics in the femoral end-to-side anastomosis: the influence of mechanical factors

This study was concerned with investigating the influence of mechanical factors on the hemodynamics of the end-to-side anastomosis in an attempt to identify critical factors and establish if it is possible to re-engineer existing, patient-specific, by-pass grafts with a view to increasing their patency. The study chose the femoral artery as the principal subject of interest. Wall shear stresses (WSS) and wall shear stress gradients (WSSG) were taken as the primary quantities of interest. Angle, graft calibre, interposition cuffs, proximal outflow and inlet waveform were studied. The study found that the use of cuffs and patches can significantly reduce abnormal WSS and WSSG by up to 70% when compared to a benchmark 45 degrees conventional anastomosis. The Taylor patch was found to be more robust in reducing peak WSS magnitudes and gradients than the Miller cuff, where design variables proved to be more critical. On the addition of a Taylor patch to a realistic end-to-side femoral anastomosis, the peak WSS and WSSG were found to be reduced by 27% and 57%, respectively. In conclusion, it is possible to use idealised models to identify critical disease influencing factors and to use these findings to reduce the effects of abnormal hemodynamics in realistic, patient-dependant models.

Local and Global Geometric Influence on Steady Flow in Distal Anastomoses of Peripheral Bypass Grafts

We consider the effect of geometrical configuration on the steady flow field of representative geometries from an in vivo anatomical data set of end-to-side distal anastomoses constructed as part of a peripheral bypass graft. Using a geometrical classification technique, we select the anastomoses of three representative patients according to the angle between the graft and proximal host vessels (GPA) and the planarity of the anastomotic configuration. The geometries considered include two surgically tunneled grafts with shallow GPAs which are relatively planar but have different lumen characteristics, one case exhibiting a local restriction at the perianastomotic graft and proximal host whilst the other case has a relatively uniform cross section. The third case is nonplanar and characterized by a wide GPA resulting from the graft being constructed superficially from an in situ vein. In all three models the same peripheral resistance was imposed at the computational outflows of the distal and proximal host vessels and this condition, combined with the effect of the anastomotic geometry, has been observed to reasonably reproduce the in vivo flow split. By analyzing the flow fields we demonstrate how the local and global geometric characteristics influences the distribution of wall shear stress and the steady transport of fluid particles. Specifically, in vessels that have a global geometric characteristic we observe that the wall shear stress depends on large scale geometrical factors, e.g., the curvature and planarity of blood vessels. In contrast, the wall shear stress distribution and local mixing is significantly influenced by morphology and location of restrictions, particular when there is a shallow GPA. A combination of local and global effects are also possible as demonstrated in our third study of an anastomosis with a larger GPA. These relatively simple observations highlight the need to distinguish between local and global geometric influences for a given reconstruction. We further present the geometrical evolution of the anastomoses over a series of follow-up studies and observe how the lumen progresses towards the faster bulk flow of the velocity in the original geometry. This mechanism is consistent with the luminal changes in recirculation regions that experience low wall shear stress. In the shallow GPA anastomoses the proximal part of the native host vessel occludes or stenoses earlier than in the case with wide GPA. A potential contribution to this behavior is suggested by the stronger mixing that characterizes anastomoses with large GPA.

Computational modeling of vascular anastomoses

Recent development of computational technology allows a level of knowledge of biomechanical factors in the healthy or pathological cardiovascular system that was unthinkable a few years ago. In particular, computational fluid dynamics (CFD) and computational structural (CS) analyses have been used to evaluate specific quantities, such as fluid and wall stresses and strains, which are very difficult to measure in vivo. Indeed, CFD and CS offer much more variability and resolution than in vitro and in vivo methods, yet computations must be validated by careful comparison with experimental and clinical data. The enormous parallel development of clinical imaging such as magnetic resonance or computed tomography opens a new way toward a detailed patient-specific description of the actual hemodynamics and structural behavior of living tissues. Coupling of CFD/CS and clinical images is becoming a standard evaluation that is expected to become part of the clinical practice in the diagnosis and in the surgical planning in advanced medical centers. This review focuses on computational studies of fluid and structural dynamics of a number of vascular anastomoses: the coronary bypass graft anastomoses, the arterial peripheral anastomoses, the arterio-venous graft anastomoses and the vascular anastomoses performed in the correction of congenital heart diseases.

On the Existence of an Optimum End-to-side Junctional Geometry in Peripheral Bypass Surgery—A Computer Generated Study

BACKGROUND

To investigate hemodynamic flow changes associated with vein cuffs and patches that may be responsible for improved patency of prosthetic infrainguinal grafts.

METHODS

The role of the graft-artery junction angle was examined by computational fluid dynamics to assess the influence of anastomotic geometry on wall shear stress (WSS) distributions. Three geometrically different junction configurations were studied and the WSS and WSS gradient (WSSG) values were compared.

RESULTS

The inclusion of a patch or a cuff moves the bed stagnation point (BSP) distally, increasing the area on the bed of the junction which experiences a BSP and reducing the strength of the recirculation region opposite the heel of the junction by 54.8 and 50.8%, respectively. The patched geometry promotes earlier recovery of the flow in the distal outflow segment (DOS) than for the unpatched model. Also, the helical flow patterns in the DOS associated with the cuffed geometry are stronger. The net effect of these changes are that peak WSSG values for the patched and cuffed geometries are three times lower than those for the uncuffed geometry.

CONCLUSION

This study provides some additional insights into the hemodynamics of graft-artery junction geometry which may influence future clinical practice.