The effect of angle on wall shear stresses in a LIMA to LAD anastomosis: numerical modelling of pulsatile flow

Hemodynamic Impact of Stenting on Carotid Bifurcation: A Potential Role of the Stented Segment and External Carotid Artery

Carotid stenting near the bifurcation carina is associated with adverse events, especially in-stent restenosis, thrombosis, and side branch occlusion in clinical data. This study is aimed at determining the potential biomechanical mechanisms for these adverse events after carotid stenting. The patient-specific carotid models were constructed with different stenting scenarios to study the flow distribution and hemodynamic parameters, such as wall shear stress (WSS), flow velocity, relative residence time (RRT), and oscillating shear index (OSI) in the carotid bifurcation. The results suggested that the existing stents surely reduced blood flow to the external carotid artery (ECA) but enhanced local flow disturbance both in ECA and stented internal carotid artery (ICA), and the inner posterior wall of the stented ICA and the outer posterior wall of ECA might endure a relatively low level of WSS and remarkably elevated OSI and RRT. In addition, the implanted stent leads to more ECA adverse flow than ICA after stenting. While disturbed flow near the strut increased as stent length increased, blood flow and areas of local flow disturbance in ECA slightly decreased as stent length increased. In conclusion, the results revealed that ECA might be in relatively high levels of abnormal local hemodynamics after stenting, followed by stented ICA, leading to potential adverse events after intervention.

Endothelial nitric oxide synthase (eNOS) mediates neointimal thickness in arteriovenous fistulae with different anastomotic angles in rats

BACKGROUND

It is known that the anastomotic angle can influence neointimal hyperplasia and patency in arteriovenous fistulae (AVF). Endothelial nitric oxide synthase (eNOS) is released from the vascular endothelium and can inhibit neointimal hyperplasia. Therefore, here, we aimed to test the hypothesis that the manipulation of eNOS expression could influence neointimal thickness in a rat AVF model with different anastomosis angles.

METHODS

Rat carotid artery (inflow, CA) and jugular vein (outflow, JV) AVF were created with acute, blunt, or end-to-end (ETE) anastomosis angles. Aspirin was used to increase eNOS expression in the acute angle group, while N(G)-nitro-L-arginine methyl ester (L-name) was used to decrease eNOS expression in the obtuse angle group. The rats were sacrificed on day 21, and tissues were harvested and analyzed histologically and with immunostaining.

RESULTS

A larger anastomosis diameter (p < 0.016) and smaller neointimal area (p < 0.01) were observed in the obtuse and end-to-end (ETE) groups compared to in the acute group. In the acute angle group, there were more proliferating cell nuclear antigen (PCNA) and α-actin dual-positive cells (p < 0.0001) and fewer phospho (p)-eNOS-positive endothelial cells (p < 0.0001) in the neointima than in the obtuse and ETE angle groups. On treating the acute angle and blunt angle groups with aspirin and L-name, respectively, no significant differences in the neointima/lumen rate were observed (p = 0.6526) between the groups; however, there were fewer von Willebrand factor (vWF) and p-eNOS dual-positive cells in the obtuse angle group treated with L-name (p = 0.0045).

CONCLUSIONS

We demonstrated that eNOS plays an important role in neointimal hyperplasia in AVF with different anastomosis angles; further, eNOS could potentially be used as a therapeutic target in patients with AVF in the future.

Can the Wall Shear Stress Values of Left Internal Mammary Artery Grafts during the Perioperative Period Reflect the One-Year Patency?

PURPOSE

The left internal mammary artery (LIMA) is the preferred graft for coronary artery bypass grafting, but the reasoning for LIMA occlusion is unclear. We sought to examine whether the wall shear stress (WSS) values of LIMA grafts during the perioperative period reflected the 1-year patency by using combining computational fluid dynamics (CFD) and coronary computed tomography angiography (CCTA) images.

METHODS

CCTA was performed in 233 patients with LIMA graft perioperatively and 1 year later from October 2014 to May 2017. LIMA occlusion was detected in six patients at the 1-year follow-up CCTA. Two patients were excluded due to poor imaging quality. The remaining four patients were enrolled as occlusive (OCC) group, and eight patients with patent LIMA were recruited as patent (PAT) group. The WSS values of LIMA during perioperative period were calculated. LIMA graft was artificially divided into three even segments, proximal (pLIMA), middle (mLIMA) and distal (dLIMA) segments. The independent samples t-test and the Student-Newman-Keuls test were used.

RESULTS

The WSS values of dLIMA were significantly higher in the PAT group than in the OCC group (4.43 vs. 2.56, p < 0.05). The WSS values of dLIMA in the PAT group were significantly higher than pLIMA, which was absent in the OCC group.

CONCLUSIONS

A higher WSS value of the distal segment of LIMA and a higher WSS value of the distal segment compared with the proximal segment of LIMA in the PAT were observed; this tendency might be helpful in predicting the 1-year patency of LIMA.

Hemodynamic and structural effects on bypass graft for different levels of stenosis using fluid structure interaction: A prospective analysis

Although the flow dynamics have been investigated using fluid-structure interaction scheme for the internal thoracic artery-left anterior descending coronary artery (ITA-LAD) bypass graft in different cases, a detailed analysis associating different degrees of LAD stenosis and its effects on hemodynamics and structural displacement are not comprehensively studied. Therefore, the primary focus of this work is to examine and determine the correlation between the hemodynamic effects and structural variations inside the bypass graft with different degrees of LAD stenosis (0%, 30%, 50%, and 75%). Navier-Stokes equation, arbitrary Lagrangian-Eulerian, and elasticity in the solid region are implemented by coupling incompressible viscous fluid, nonlinear viscous fluid, and the stress tensor equations, respectively. Using fluid-structure interaction, variations in the hemodynamic property and changes in wall shear stress (WSS) including the spatial WSS distribution and the changes in the displacement of different degree of LAD stenosis are determined.Maximum WSS is found to be around 1.58E1 Pa near the anastomosis region and maximum magnitude for the structural displacement is found to be approximately 1.25E-5 m close to the heel. The results demonstrate that the disturbance in the flow pattern is evident mainly in the anastomosis region. Consequently, higher WSS is observed near the toe and on the artery wall, near the anastomosis region.

Transfer of Low-Density Lipoproteins in Coronary Artery Bifurcation Lesions with Stenosed Side Branch: Numerical Study

Evidence from clinical data suggests that the stenotic side branch (SB) is one of the key predictors for SB occlusion-based adverse events. In this study, we hypothesized that coronary bifurcations with stenotic SB might lead to severe concentration polarization of atherogenic lipids, such as the low-density lipoproteins (LDL), motivating the adverse events in the clinic. To confirm this hypothesis, this work numerically investigated the transport of LDL in different bifurcation lesions based on the Medina classification with various location and stenosis severities. The results showed that the coronary bifurcations with stenotic SB might be suffering more serious concentration polarization of LDL on the luminal surface of the SB due to higher level of LDL concentrations. Moreover, compared to the other bifurcation lesion types, the type (1,0,1) had the highest luminal surface LDL concentration along the SB and the highest degree of risk to enhance the process of atherosclerosis. In addition, this study also showed that the luminal surface LDL concentration increased with elevated stenosis severity. The type (1,0,1) with the severe stenosis (75% diameter reduction) had the highest concentration at the SB. In conclusion, these results suggested that both location of lesions and stenosis severities had great influence on the distribution of LDL on the luminal surface of the SB. Therefore, the estimation of disease severity and the interventional therapy should be carried out not only according to the stenosis severities in clinic. Moreover, compared to the other bifurcation lesion types, the type (1,0,1), rather than the type (1,1,1) as usually considered, had the highest luminal surface LDL concentration along the SB and the highest degree of risk to enhance the process of atherosclerosis.

EFFECT OF COUPLE STRESSES ON ELECTROKINETIC OSCILLATORY FLOW OF BLOOD IN THE MICROCIRCULATORY SYSTEM

Of concern in the paper is the development of a mathematical model of blood flow in the microcirculatory system. The study pertains to a situation, where the system is subject to the action of an external AC electric field and the flow is oscillatory. The general case when the frequency of the electric field is different from that of the flow has been studied. Blood is treated as a couple stress fluid. The flow is supposed to take place between two oscillatory plates and is considered to be of electro-osmotic type. Our study is aimed at examining the effect of couple stress at different frequencies of the capillary wall and the applied AC electric field on electroosmotic flow of blood. The computational results indicate very clearly that the electrokinetic velocity of blood reduces as the couple stress effect increases and that the velocity reduces, when the Reynolds number rises. The results of the present study will serve as a reasonably good estimate for electro-osmotic transport of blood in small blood vessels.

A multiscale approach for determining the morphology of endothelial cells at a coronary artery

The morphology of endothelial cells (ECs) may be an indication for determining atheroprone sites. Until now, there has been no clinical imaging technique to visualize the morphology of ECs in the arteries. The present study introduces a computational technique for determining the morphology of ECs. This technique is a multiscale simulation consisting of the artery scale and the cell scale. The artery scale is a fluid-structure interaction simulation. The input for the artery scale is the geometry of the coronary artery, that is, the dynamic curvature of the artery due to the cardiac motion, blood flow, blood pressure, heart rate, and the mechanical properties of the blood and the arterial wall, the measurements of which can be obtained for a specific patient. The results of the artery scale are wall shear stress (WSS) and cyclic strains as the mechanical stimuli of ECs. The cell scale is an inventive mass-and-spring model that is able to determine the morphological response of ECs to any combination of mechanical stimuli. The results of the multiscale simulation show the morphology of ECs at different locations of the coronary artery. The results indicate that the atheroprone sites have at least 1 of 3 factors: low time-averaged WSS, high angle of WSS, and high longitudinal strain. The most probable sites for atherosclerosis are located at the bifurcation region and lie on the myocardial side of the artery. The results also indicated that a higher dynamic curvature is a negative factor and a higher pulse pressure is a positive factor for protection against atherosclerosis.

High shear stress induces atherosclerotic vulnerable plaque formation through angiogenesis

Rupture of atherosclerotic plaques causing thrombosis is the main cause of acute coronary syndrome and ischemic strokes. Inhibition of thrombosis is one of the important tasks developing biomedical materials such as intravascular stents and vascular grafts. Shear stress (SS) influences the formation and development of atherosclerosis. The current review focuses on the vulnerable plaques observed in the high shear stress (HSS) regions, which localizes at the proximal region of the plaque intruding into the lumen. The vascular outward remodelling occurs in the HSS region for vascular compensation and that angiogenesis is a critical factor for HSS which induces atherosclerotic vulnerable plaque formation. These results greatly challenge the established belief that low shear stress is important for expansive remodelling, which provides a new perspective for preventing the transition of stable plaques to high-risk atherosclerotic lesions.

DIRECT CORONARY COUPLING APPROACH FOR COMPUTING FFRCT

With the advances in computational fluid dynamics (CFD) and image-based modeling techniques, fractional flow reserve (FFR) can be computed from coronary computed tomography angiography (CTA) scans (FFRCT). However, this non-invasive approach requires large-scale computational resources, which limits its application in routine clinical setting. A 3D–0D coupling approach is proposed to improve the coupling efficiency of FFRCT. Aortic–root is modeled by a lumped parameter model and connected with the models of left ventricle and systemic circulation. With this approach, the interested coronary regions can be directly coupled to the lumped parameter model, resulting in a significant reduction (up to 20 times reduction) in the volume of the CFD computing domain. The proposed approach is applied to a patient-specific model and compared with previous non-reduced approach. Results show that the computed coronary flow rates, pressure waveforms and FFRCT contours by the proposed approach are consistent well with that of the non-reduced approach. These results demonstrate that the proposed approach can reduce the CFD computing domain of FFRCT significantly while maintaining the similar accuracy as compared with the non-reduced approach, and it can be further employed to promote FFRCT in routine clinical setting.

The numerical analysis of non-Newtonian blood flow in human patient-specific left ventricle

Recently, various non-invasive tools such as the magnetic resonance image (MRI), ultrasound imaging (USI), computed tomography (CT), and the computational fluid dynamics (CFD) have been widely utilized to enhance our current understanding of the physiological parameters that affect the initiation and the progression of the cardiovascular diseases (CVDs) associated with heart failure (HF). In particular, the hemodynamics of left ventricle (LV) has attracted the attention of the researchers due to its significant role in the heart functionality. In this study, CFD owing its capability of predicting detailed flow field was adopted to model the blood flow in images-based patient-specific LV over cardiac cycle. In most published studies, the blood is modeled as Newtonian that is not entirely accurate as the blood viscosity varies with the shear rate in non-linear manner. In this paper, we studied the effect of Newtonian assumption on the degree of accuracy of intraventricular hemodynamics. In doing so, various non-Newtonian models and Newtonian model are used in the analysis of the intraventricular flow and the viscosity of the blood. Initially, we used the cardiac MRI images to reconstruct the time-resolved geometry of the patient-specific LV. After the unstructured mesh generation, the simulations were conducted in the CFD commercial solver FLUENT to analyze the intraventricular hemodynamic parameters. The findings indicate that the Newtonian assumption cannot adequately simulate the flow dynamic within the LV over the cardiac cycle, which can be attributed to the pulsatile and recirculation nature of the flow and the low blood shear rate.

WAVE PROPAGATION IN A 1D FLUID DYNAMICS MODEL USING PRESSURE-AREA MEASUREMENTS FROM OVINE ARTERIES

This study considers a 1D fluid dynamics arterial network model with 14 vessels developed to assimilate ex vivo 0D temporal data for pressure-area dynamics in individual vessel segments from 11 male Merino sheep. A 0D model was used to estimate vessel wall parameters in a two-parameter elastic model and a four-parameter Kelvin viscoelastic model. This was done using nonlinear optimization minimizing the least squares error between model predictions and measured cross-sectional areas. Subsequently, estimated values for elastic stiffness and unstressed area were related to construct a nonlinear relationship. This relation was used in the network model. A 1D single vessel model of the aorta was then developed and used to estimate the inflow profile and parameters for total resistance and compliance for the downstream network and to demonstrate effects of incorporating viscoelasticity in the arterial wall. Lastly, the extent to which vessel wall parameters estimated from ex vivo data can be used to realistically simulate pressure and area in a vessel network was evaluated. Elastic wall parameters in the network simulations were found to yield pressure-area relationships across all vessel locations and sheep that were in ranges comparable to those in the ex vivo data.

FSI SIMULATION OF A HEALTHY CORONARY BIFURCATION FOR STUDYING THE MECHANICAL STIMULI OF ENDOTHELIAL CELLS UNDER DIFFERENT PHYSIOLOGICAL CONDITIONS

Atherosclerosis is a world-spread and well-known disease. This disease strongly relates to the endothelial cells (ECs) function. Normally, the endothelial cells align in the flow direction in the atheroprotected sites; however, in the case of atheroprone sites these cells orient randomly. The mechanical stimuli such as wall shear stress and strains could determine the morphology and function of the endothelial cells. In the present study, we numerically simulated the left main coronary artery (LCA) and its branches to left anterior descending (LAD) and left circumflex coronary (LCX) artery using fluid–structure interaction (FSI) modeling. The results were presented as longitudinal and circumferential strains of ECs as well as wall shear stress. Wide ranges of heart rate, cardiac motion, systolic and diastolic pressures were considered and their effects on mechanical stimuli were described in detail. The results showed that these factors could greatly influence the risk of atherosclerosis and the location of atherosclerotic lesions.

EFFECT OF CEREBROSPINAL FLUID MODELED WITH DIFFERENT MATERIAL PROPERTIES ON A HUMAN FINITE ELEMENT HEAD MODEL

The aim of this study was to enhance head-injury prediction, this paper investigated the behavior of cerebrospinal fluid (CSF) in finite element (FE) modeling. Nine different material properties selected according to material definitions and property values were used to represent CSF in FE head models. To evaluate the influence of CSF material parameters on brain mechanical responses, the models were validated against available cadaver experiment data. Results showed that coup pressure increased whereas contrecoup pressure decreased when the head sustained purely translational impact with increased bulk modulus when CSF was modeled as fluid. However, with increased bulk modulus, coup pressure, contrecoup pressure and relative skull-brain motions decreased under rotational impact. When CSF was assumed to be an elastic material, coup pressure increased whereas contrecoup pressure decreased with increased elastic modulus when the head was subjected to purely translational impact. However, the variation trend was not obvious during head rotation. Results also indicated that when subjected to brain strain and von Mises stress, the model was prone to underestimate brain injury when CSF was modeled as an elastic material, especially during purely translational impact to the head. The model with CSF as fluid reduced the strain rate of brain, which was more likely to be realistic than the model with CSF as a viscoelastic material. These findings suggested that significantly higher values of the bulk modulus of CSF modeled as fluid were needed to predict intracranial dynamic responses and brain injury during head impact.

Novel modular anastomotic valve device for hemodialysis vascular access: preliminary computational hemodynamic assessment

PURPOSE

Arteriovenous graft patency is limited by terminal occlusion caused by intimal hyperplasia (IH). Motivated by evidence that flow disturbances promote IH progression, a modular anastomotic valve device (MAVD) was designed to isolate the graft from the circulation between dialysis periods (closed position) and enable vascular access during dialysis (open position). The objective of this study was to perform a preliminary computational assessment of the device ability to normalize venous flow between dialysis periods and potentially limit IH development and thrombogenesis.

METHODS

Computational fluid dynamics simulations were performed to compare flow and wall shear stress (WSS) in a native vein and MAVD prototypes featuring anastomotic angles of 90° and 30°. Low WSS (LWSS) regions prone to IH development were characterized in terms of temporal shear magnitude (TSM), oscillatory shear index (OSI), and relative residence time (RRT). Thrombogenic potential was assessed by investigating the loading history of fluid particles traveling through the device.

RESULTS

The closed MAVD exhibited the same flow characteristics as the native vein (0.3% difference in pressure drop, 3.5% difference in surface-averaged WSS). The open MAVD generated five LWSS regions (TSM <0.5 Pa) exhibiting different degrees of flow reversal (surface-averaged OSI: 0.03-0.36) and stagnation (max RRT: 2.50-37.16). Reduction in anastomotic angle resulted in the suppression of three LWSS regions and overall reductions in flow reversal (surface-averaged OSI <0.21) and stagnation (max RRT <18.05).

CONCLUSIONS

This study suggests the ability of the MAVD to normalize venous flow between dialysis periods while generating the typical hemodynamics of end-to-side vein-graft anastomoses during dialysis.

NUMERICAL STUDY OF BIDIRECTIONAL GLENN WITH UNILATERAL PULMONARY ARTERY STENOSIS

Purpose: Hypoplastic left heart syndrome (HLHS) is a congenital heart disease and is usually associated with pulmonary artery stenosis. The superior vena cava-to-pulmonary artery (bidirectional Glenn) shunt is used primarily as a staging procedure to the total cava-to-pulmonary connection for single-ventricle complex. When HLHS coexists with pulmonary artery stenosis, the surgeons then face a multiple problem. This leads to high demand of optimized structure of Glenn surgery. The objective of this article is to investigate the influence of various anastomotic structures and the direction of superior vena cava (SVC) in Glenn on hemodynamics under pulse inflow conditions and try to find an optimal structure of SVC in Glenn surgery with unilateral pulmonary artery stenosis.

Method: First, 3D patient-specific models were constructed from medical images of a HLHS patient before any surgery by using the commercial software Mimics, and another software Free-form was used to deform the reconstructed models in the computer. Four 3D patient-specific Glenn models were constructed: model-1 (normal Glenn), model-2 (lean the SVC back to the stenotic pulmonary artery), model-3 (lean the SVC towards the stenotic pulmonary artery), model-4 (add patch at junction of the SVC toward stenosis at pulmonary artery). Second, a lumped parameter model (LPM) was established to predict boundary conditions for computational fluid dynamics (CFD). In addition, numerical simulations were conducted using CFD through the finite volume method. Finally, hemodynamic parameters were obtained and evaluated.

Results: It was showed that model-4 have relatively balanced vena cava blood perfusion into the left pulmonary artery (LPA) and right pulmonary artery (RPA), this may be due to less helical flow and the patch at junction of the SVC. Near stenosis of pulmonary artery, model-4 performed with the higher wall shear stress (WSS), which would benefit endothelial cell function and gene expression. In addition, results showed that model-4 performed with the lower oscillatory shear index (OSI) and wall shear stress gradient (WSSG), which would decrease the opportunity of vascular intimal hyperplasia.

Conclusion: It is benefited that surgeons adds patch at junction of the SVC towards stenosis at pulmonary artery. These results can impact the surgical design and planning of the Glenn surgery with unilateral pulmonary artery stenosis.

INTIMAL HYPERPLASIA AND WALL SHEAR IN ARTERIAL BYPASS Y-GRAFTING AND CONSEQUENCE GRAFTING: A NUMERICAL STUDY

The progression of intimal hyperplasia is considered to be the main cause of bypass failure and is directly related to the individual blood rheology, local arterial geometry and placement of the junctions, graft diameter and graft surface characteristics as well as the degree of compliance. In this paper we use commercial computational fluid dynamics (CFD) ANSYS to examine under the correct physiological flow conditions the hemodynamic forces of composite bypass with internal mammary artery in Y-grafting and consequence grafting which is known to achieve high patency rate and highly recommended by clinicians. Particular emphasis is given here on the parameters that could initiate the development of intimal hyperplasia within these bypass configurations. The hemodynamic flow patterns between the consequence grafting and the composite Y-grafting are observed here to be different. Moreover, on both end-to-side and side-to-side configurations, the circulating flows are detected in the vicinity of the junction area, while the Dean flow vortexes are only observed on the end-to-side configuration. Likewise, the hemodynamic flow on the end-to-side configuration on the LCX of both 45° and 90° Y-grafting is found to be smoother than that of the junction on the LCA, regardless of the changing of anastomosis angles. The high WSS gradients are observed at the vicinity of the toe and on the bed of the junction, while the low WSS are presented at the distal of the stenosis and at the stagnation point. The clinical relevance of the results are presented and discussed with particular focus on the factors and the flow patterns that trigger the development of intimal hyperplasia.

Fluid structure interaction (FSI) simulation of the left ventricle (LV) during the early filling wave (E-wave), diastasis and atrial contraction wave (A-wave)

In this paper, the hemodynamic characteristics inside a physiologically correct three-dimensional LV model using fluid structure interaction scheme are examined under various heart beat conditions during early filling wave (E-wave), diastasis and atrial contraction wave (A-wave). The time dependent and incompressible viscous fluid, nonlinear viscous fluid and the stress tensor equations are coupled with the full Navier-Stoke's equations together with the Arbitrary Lagrangian-Eulerian and elasticity in the solid domain are used in the analysis. The results are discussed in terms of the variation in the intraventricular pressure, wall shear stress (WSS) and the fluid flow patterns inside the LV model. Moreover, changes in the magnitude of displacements on the LV are also observed during diastole period. The results obtained demonstrate that the magnitude of the intraventricle pressure is found higher in the basal region of the LV during the beginning of the E-wave and A-wave, whereas the Ip is found much higher in the apical region when the flow propagation is in peak E-wave, peak A-wave and diastasis. The magnitude of the pressure is found to be 5.4E2 Pa during the peak E-wave. Additionally, WSS elevates with the rise in the E-wave and A-wave but the magnitude decreases during the diastasis phase. During the peak E-wave, maximum WSS is found to be 5.7 Pa. Subsequent developments, merging and shifting of the vortices are observed throughout the filling wave. Formations of clockwise vortices are evident during the peak E-wave and at the onset of the A-wave, but counter clockwise vortices are found at the end of the diastasis and at the beginning of the A-wave. Moreover, the maximum magnitude of the structural displacement is seen in the ventricle apex with the value of 3.7E-5 m.

The accelerated atherogenesis of venous grafts might be attributed to aggravated concentration polarization of low density lipoproteins: a numerical study

We hypothesize that after implantation the much elevated water filtration rate of venous grafts may cause aggravated concentration polarization of low density lipoproteins (LDLs), in turn lead to the accelerated atherogenesis of the grafts. To verify the hypothesis, we numerically simulated the transport of LDLs in various models of arterial bypasses with different grafts (veins or arteries) and geometrical configurations. The results showed that the venous grafts might endure abnormally high lipid infiltration/accumulation within the vessel wall due to severely elevated luminal surface LDL concentration. When compared to the conventional bypass models, the S-type bypass had the lowest luminal surface LDL concentration along its host artery floor, but the highest degree of risk to develop atherosclerotic lesions in its venous graft. Among the three conventional bypass models, the one with 30° anastomosis had the lowest risk to develop atherosclerosis in the venous graft. In conclusion, when compared with the bypass models with arterial grafts, the venous bypass models had rather high levels of LDL concentration polarization (cw) in the vein grafts, especially at the early stages of implantation. This might result in high infiltration/accumulation of LDLs within the walls of the venous grafts, leading to a fast genesis/development of atherosclerosis there.

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 anastomosis angle on the localization of disturbed flow in 'side-to-end' fistulae for haemodialysis access

BACKGROUND

Early failure of the vascular access for haemodialysis (HD) after the surgical creation of a radial-cephalic arteriovenous fistula (AVF) occurs mainly due to a juxta-anastomotic stenosis. Even if elevated blood flow induces high wall shear stress, we have recently shown that disturbed flow, characterized by low and reciprocating flow, may develop in zones of the AVF where it can provide a good indication of the sites of future stenoses. The present study was aimed at investigating whether the anastomosis angle influences disturbed flow in radial-cephalic 'side-to-end' AVF.

METHODS

By means of a parametric AVF model we created four equivalent meshes with anastomosis angles of 30°, 45°, 60° and 90°, respectively. We then performed transient, non-Newtonian computational fluid dynamics simulations using, as boundary conditions, previously measured blood volume flow and division ratio in subjects requiring primary access. The relative residence time (RRT), a robust indicator of disturbed flow, was calculated for the overall wall surface and disturbed flow was localized as areas having RRT > 1. Quantitative characterization and statistical tests were employed to assess the difference in RRT medians between the four anastomosis angle cases.

RESULTS

Disturbed flow was located in all AVF models in the same areas where flow recirculation and stagnation occurred, on the inner wall of the swing segment (SS) and on the arterial wall at the anastomosis floor (AF). A smaller angle AVF had smaller disturbed flow areas with lower RRT peak values, either on the venous or the arterial limb. There were significant differences in the RRT medians on the SS and on the AF between sharper (30° and 45°) and wider (60° or 90°) angles.

CONCLUSIONS

We have found that in 'side-to-end' radial-cephalic AVFs for HD, the anastomosis angle does impact on the local disturbed flow patterns. Among the four geometries we considered in this study, the smaller angle (30°) would be the preferred choice that minimizes the development of neointima. Clinicians should consider this at the time of AVF creation because the anastomosis angle is in part amenable to surgical manipulation.

Effect of the degree of LAD stenosis on "competitive flow" and flow field characteristics in LIMA-to-LAD bypass surgery

The long-term patency of the left internal mammary artery (LIMA) in left anterior descending (LAD) coronary stenosis bypass surgery is believed to be related to the degree of competitive flow between the LAD and LIMA. To investigate the effect of the LAD stenosis severity on this phenomenon and on haemodynamics in the LIMA and anastomosis region, a numerical LIMA-LAD model was developed based on 3D geometric (obtained from a cast) and hemodynamic data from an experimental pig study. Proximal LAD pressure was used as upstream boundary condition. The model counted 13 outlets (12 septal arteries and the distal LAD) where flow velocities were imposed in systole, while myocardial conductance was imposed in diastole via an implicit scheme. LAD stenoses of 100 (total occlusion), 90, 75 and 0 % area reduction were constructed. Low degree of LAD stenosis was associated with highly competitive flow and low wall shear stress (WSS) in the LIMA, an unfavourable hemodynamic regime which might contribute to WSS-related remodelling of the LIMA and suboptimal long-term LIMA bypass performance.

PIV MEASUREMENTS AND NUMERICAL VALIDATION OF END-TO-SIDE ANASTOMOSIS

In this article, particle image velocimetry (PIV) technique was used to determine the instantaneous velocity fields inside a model of end-to-side anastomosis under various physiological flow conditions. Using ANSYS software, a three-dimensional (3D) computational model at the peak systolic blood flow was simulated. The numerical and experimental results were presented and discussed in terms of velocity fields at various locations along the graft and the host artery. The numerical results were then compared with the experimental data and a large difference was found, which was attributed to the imperfection of manufacturing the glass model and measurements error associated with PIV. The findings indicated in general that the analysis at peak systole, steady flow could help in providing essential quantitative information of the hemodynamics in anastomotic artery.

Influences of Geometric Configurations of Bypass Grafts on Hemodynamics in End-to-Side Anastomosis

Background

Although considerable efforts have been made to improve the graft patency in coronary artery bypass surgery, the role of biomechanical factors remains underrecognized. The aim of this study is to investigate the influences of geometric configurations of the bypass graft on hemodynamic characteristics in relation to anastomosis.

Materials and Methods

The Numerical analysis focuses on understanding the flow patterns for different values of inlet and distal diameters and graft angles. The Blood flow field is treated as a two-dimensional incompressible laminar flow. A finite volume method is adopted for discretization of the governing equations. The Carreau model is employed as a constitutive equation for blood. In an attempt to obtain the optimal aorto-coronary bypass conditions, the blood flow characteristics are analyzed using in vitro models of the end-to-side anastomotic angles of 45°, 60° and 90°. To find the optimal graft configurations, the mass flow rates at the outlets of the four models are compared quantitatively.

Results

This study finds that Model 3, whose bypass diameter is the same as the inlet diameter of the stenosed coronary artery, delivers the largest amount of blood and the least pressure drop along the arteries.

Conclusion

Biomechanical factors are speculated to contribute to the graft patency in coronary artery bypass grafting.

Graft-artery junctions: design optimization and CAD development

Designing and manufacturing of vascular prosthesis for arterial bypass grafts is a very complex problem. The process involves the selection of suitable geometry, materials of appropriate characteristics, and manufacturing technique capable of constructing prosthesis in a cost-effective manner. In this chapter, all engineering aspects related to the design and optimization of an artificial graft are presented and discussed. These aspects include CAD design of the graft, in vitro hemodynamic analysis to ensure good mechanical integrity and functionality, and optimization of the manufacturing techniques. Brief discussion is also given on the endothelization and vascularization of the artificial vessels and the future directions of the development of synthetic vessels for human implementation.

Computational fluid dynamics in cardiovascular disease

Computational fluid dynamics (CFD) is a mechanical engineering field for analyzing fluid flow, heat transfer, and associated phenomena, using computer-based simulation. CFD is a widely adopted methodology for solving complex problems in many modern engineering fields. The merit of CFD is developing new and improved devices and system designs, and optimization is conducted on existing equipment through computational simulations, resulting in enhanced efficiency and lower operating costs. However, in the biomedical field, CFD is still emerging. The main reason why CFD in the biomedical field has lagged behind is the tremendous complexity of human body fluid behavior. Recently, CFD biomedical research is more accessible, because high performance hardware and software are easily available with advances in computer science. All CFD processes contain three main components to provide useful information, such as pre-processing, solving mathematical equations, and post-processing. Initial accurate geometric modeling and boundary conditions are essential to achieve adequate results. Medical imaging, such as ultrasound imaging, computed tomography, and magnetic resonance imaging can be used for modeling, and Doppler ultrasound, pressure wire, and non-invasive pressure measurements are used for flow velocity and pressure as a boundary condition. Many simulations and clinical results have been used to study congenital heart disease, heart failure, ventricle function, aortic disease, and carotid and intra-cranial cerebrovascular diseases. With decreasing hardware costs and rapid computing times, researchers and medical scientists may increasingly use this reliable CFD tool to deliver accurate results. A realistic, multidisciplinary approach is essential to accomplish these tasks. Indefinite collaborations between mechanical engineers and clinical and medical scientists are essential. CFD may be an important methodology to understand the pathophysiology of the development and progression of disease and for establishing and creating treatment modalities in the cardiovascular field.

Study of flow-induced hemolysis using novel Couette-type blood-shearing devices

To assist the development and application of blood-contacting medical devices, two novel flow-through Couette-type blood-shearing devices have been developed to study the quantitative relationship between blood damage indexes and flow-dependent parameters. One device is an axial flow-through Couette-type device supported by a pair of pin bearings adapted from the adult Jarvik 2000 blood pump. The other is a centrifugal flow-through Couette-type device supported with magnetic bearings adapted from the CentriMag blood pump. In both devices, a rotor spindle was used to replace the original impeller blades so that a small gap was created between the housing and the rotating spindle surface. Computational fluid dynamics simulations have shown that a uniform, high shear stress region can be generated inside the small gap while the shear stresses elsewhere are relatively low. The possibility of secondary blood damage caused by mechanical seals was eliminated due to the use of a magnetic rotor system. Blood flow through the gap was driven by an externally pressurized reservoir. By adjusting the rotational speed and blood flow rate, shear-induced hemolysis was quantified at a matrix of exposure time (0.039 to 1.48 s) and shear stress (50 to 320 Pa). All of the experiments were conducted at room temperature using heparinized ovine blood with a hematocrit value of 30%. The measured hemolysis levels were much lower than those published in the literature, and the overestimation of those earlier studies may be attributable to device-related secondary blood-damaging effects. A new set of coefficients for the power law model was derived from the regression of the experimental data.

Biomechanics of coronary artery and bypass graft disease: potential new approaches

The contribution of biomechanical factors to the incidence and distribution of coronary artery and bypass graft disease is underrecognized. This review examined the literature to determine which factors were relevant and the evidence for their importance. It identified two primary biomechanical factors that predispose to disease: (1) low-wall shear stress and (2) high-wall mechanical stress or strain. A range of secondary biomechanical factors have also been identified and include: vessel geometry; vessel movement; vessel wall characteristics and the presence of reflection waves. Potential surgical approaches for minimizing these effects are discussed.

Haemodynamic analysis of coronary artery bypass grafting in a non-linear deformable artery and Newtonian pulsatile blood flow

A three-dimensional (3D) computational model of stenotic coronary artery bypass grafting (CABG) system with fluid—structure interaction (FSI) using realistic physiological conditions is introduced. Unsteady pulsatile blood flow is applied to the wall of non-linear deformable arteries over the systolic period. In the analysis, the arbitrarily Lagrangian—Eulerian (ALE) formulation is used to couple the fluid region and solid domain. The method couples the equations of the deformation of the artery wall and applies them as the fluid domain boundary condition. The flow distribution and haemodynamic forces are presented in terms of velocity profiles and temporal and spatial wall shear stresses (WSSs) at the distal area. Rapid changes in the flow fields are observed in the early stages of the cardiac cycle, which alters the location of the recirculation zone from the toe to the host bed and then to the heel. The migration of the recirculation zone, considering the effect of deformability of the artery wall, indicates the same trend as the rigid wall model according to the location of low and high WSSs. However, the WSSs in the critical areas such as toe, heel, and suture lines are found to have dramatic drops in magnitudes in comparison with those of the rigid wall model. This could initiate the promotion of intimal hyperplasia (IH) and may cause an early graft failure in CABG.

Acute myocardial infarction of the right coronary artery originating from the distal left circumflex artery

A case of acute myocardial infarction (AMI) of the distal left circumflex, near the origin of an aberrant right coronary artery is presented. Coronary stenting was successfully performed. According to several reports, this anomaly is a common site for coronary atherosclerosis, but this is the first report of AMI.

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.