Computational modeling of vascular anastomoses

In vivo evaluation of compliance mismatch on intimal hyperplasia formation in small diameter vascular grafts

Small diameter synthetic vascular grafts have high failure rate due to the thrombosis and intimal hyperplasia formation. Compliance mismatch between the synthetic graft and native artery has been speculated to be one of the main causes of intimal hyperplasia. However, changing the compliance of synthetic materials without altering material chemistry remains a challenge. Here, we used poly(vinyl alcohol) (PVA) hydrogel as a graft material due to its biocompatibility and tunable mechanical properties to investigate the role of graft compliance in the development of intimal hyperplasia and in vivo patency. Two groups of PVA small diameter grafts with low compliance and high compliance were fabricated by dip casting method and implanted in a rabbit carotid artery end-to-side anastomosis model for 4 weeks. We demonstrated that the grafts with compliance that more closely matched with rabbit carotid artery had lower anastomotic intimal hyperplasia formation and higher graft patency compared to low compliance grafts. Overall, this study suggested that reducing the compliance mismatch between the native artery and vascular grafts is beneficial for reducing intimal hyperplasia formation.

New Computational Approach to Shunt Design in Congenital Heart Palliation

Shunts are commonly used to redirect blood to pulmonary arteries in procedures that palliate congenital cardiovascular defects. Previous clinical studies and hemodynamic simulations reveal a critical role of shunt diameter in balancing flow to pulmonary versus systemic vessels, but the biomechanical process of creating the requisite anastomosis between the shunt and host vessel has received little attention. Here, we report a new Lagrange multiplier-based finite element approach that represents the shunt and host vessels as individual structures and predicts the anastomosis geometry and attachment force that result when the shunt is sutured at an incision in the host, followed by pressurization. Simulations suggest that anastomosis orifice opening increases markedly with increasing length of the host incision and moderately with increasing blood pressure. The host artery is further predicted to conform to common stiff synthetic shunts, whereas more compliant umbilical vessel shunts should conform to the host, with orifice area transitioning between these two extremes via a Hill-type function of shunt stiffness. Moreover, a direct relationship is expected between attachment forces and shunt stiffness. This new computational approach promises to aid in surgical planning for diverse vascular shunts by predicting in vivo pressurized geometries.

Functional analysis of arteriovenous fistulae in non-contrast magnetic resonance images

BACKGROUND AND OBJECTIVE

Arteriovenous fistulae (AVF) are the preferred mode of hemodialysis vascular access and their successful maturation is critical to reduce patient morbidity, mortality, cost, and improve quality of life. Peri-anastomotic venous segment stenosis is the primary cause of AVF maturation failure. The objective is to develop a software protocol for the functional analysis of arteriovenous fistula.

METHOD

We have developed a standard protocol for the anatomical analysis of the AVF to better understand the mechanisms involved in AVF stenosis and to identify future imaging biomarkers for AVF success or failure using non-contrast magnetic resonance imaging (MRI). The 3D model of the AVF is created using a polar dynamic programming technique. Analysis has been performed on six Yorkshire cross domestic swine, but techniques can be applied into clinical settings.

RESULTS

Differences in AVF angles and vein curvature are associated with significant variability of venous cross-sectional area. This suggests that the pattern of stenosis is likely to be dependent upon hemodynamic profiles which are largely determined by AVF anatomical features and could play an important role in AVF maturation.

CONCLUSIONS

This protocol enables us to visualize and study the hemodynamic profiles indirectly allowing early stratification of patients into high and low risk groups for AVF maturation failure. High risk patients could then be targeted with an enhanced process of care or future maturation enhancing therapies resulting in a much-needed precision-medicine approach to dialysis vascular access.

Arteriovenous conduits for hemodialysis: how to better modulate the pathophysiological vascular response to optimize vascular access durability

Vascular access is the lifeline for patients on hemodialysis. Arteriovenous fistulas (AVFs) are the preferred vascular access, but AVF maturation failure remains a significant clinical problem. Currently, there are no effective therapies available to prevent or treat AVF maturation failure. AVF maturation failure frequently results from venous stenosis at the AVF anastomosis, which is secondary to poor outward vascular remodeling and excessive venous intimal hyperplasia that narrows the AVF lumen. Arteriovenous grafts (AVGs) are the next preferred vascular access when an AVF creation is not possible. AVG failure is primarily the result of venous stenosis at the vein-graft anastomosis, which originates from intimal hyperplasia development. Although there has been advancement in our knowledge of the pathophysiology of AVF maturation and AVG failure, this has not translated into effective therapies for these two important clinical problems. Further work will be required to dissect out the mechanisms of AVF maturation failure and AVG failure to develop more specific therapies. This review highlights the major recent advancements in AVF and AVG biology, reviews major clinical trials, and discusses new areas for future research.

Influence of microvascular sutures on shear strain rate in realistic pulsatile flow

Arterial thrombus formation is directly related to the mechanical shear experienced by platelets within flow. High shear strain rates (SSRs) and large shear gradients cause platelet activation, aggregation and production of thrombus. This study, for the first time, investigates the influence of pulsatile flow on local haemodynamics within sutured microarterial anastomoses. We measured physiological arterial waveform velocities experimentally using Doppler ultrasound velocimetry, and a representative example was applied to a realistic sutured microarterial geometry. Computational geometries were created using measurements taken from sutured chicken femoral arteries. Arterial SSRs were predicted using computational fluid dynamics (CFD) software, to indicate the potential for platelet activation, deposition and thrombus formation. Predictions of steady and sinusoidal inputs were compared to analyse whether the addition of physiological pulse characteristics affects local intravascular flow characteristics. Simulations were designed to evaluate flow in pristine and hand-sutured microarterial anastomoses, each with a steady-state and sinusoidal pulse component. The presence of sutures increased SSR_max in the anastomotic region by factors of 2.1 and 2.3 in steady-state and pulsatile flows respectively, when compared to a pristine vessel. SSR values seen in these simulations are analogous to the presence of moderate arterial stenosis. Steady-state simulations, driven by a constant inflow velocity equal to the peak systolic velocity (PSV) of the measured pulsatile flow, underestimated SSRs by ∼ 9% in pristine, and ∼ 19% in sutured vessels compared with a realistic pulse. Sinusoidal flows, with equivalent frequency and amplitude to a measured arterial waveform, represent a slight improvement on steady-state simulations, but still SSRs are underestimated by 1-2%. We recommend using a measured arterial waveform, of the form presented here, for simulating pulsatile flows in vessels of this nature. Under realistic pulsatile flow, shear gradients across microvascular sutures are high, of the order ∼ 7.9 × 10⁶ m^-1 s^-1, which may also be associated with activation of platelets and formation of aggregates.

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.

Dynamic behavior of suture-anastomosed arteries and implications to vascular surgery operations

BACKGROUND

Routine vascular surgery operations involve stitching of disconnected human arteries with themselves or with artificial grafts (arterial anastomosis). This study aims to extend current knowledge and provide better-substantiated understanding of the mechanics of end-to-end anastomosis through the development of an analytical model governing the dynamic behavior of the anastomotic region of two initially separated arteries.

METHODS

The formulation accounts for the arterial axial-circumferential deformation coupling and suture-artery interaction. The proposed model captures the effects of the most important parameters, including the geometric and mechanical properties of artery and sutures, number of sutures, loading characteristics, longitudinal residual stresses, and suture pre-tensioning.

RESULTS

Closed-form expressions are derived for the system response in terms of arterial radial displacement, anastomotic gap, suture tensile force, and embedding stress due to suture-artery contact interaction. Explicit objective functionalities are established to prevent failure at the anastomotic interface.

CONCLUSIONS

The mathematical formulation reveals useful interrelations among the problem parameters, thus making the proposed model a valuable tool for the optimal selection of materials and improved functionality of the sutures. By virtue of their generality and directness of application, the findings of this study can ultimately form the basis for the development of vascular anastomosis guidelines pertaining to the prevention of post-surgery implications.

Metric field construction for anisotropic mesh adaptation with application to blood flow simulations

The goal of this paper is to generate an anisotropic metric field suitable for cardiovascular geometries before a fluid simulation. Starting from a curvature map, an initial surface metric field is computed. This metric is used for anisotropic surface mesh adaptation and consecutively extended inside the volume in a frontal manner. The algorithm is based on the method proposed by Alauzet but replaces the metric intersection steps by an original metric 'blending'. This allows for a graded anisotropic volume mesh with a refinement layer close to the walls. The benefits of the resulting mesh are multiple: a reduced number of degrees of freedom, a priori refinement in areas with strong gradients of velocity and automatically increased resolution in regions with high surface curvature. The primal application of this method is in the domain of cardiovascular flows. Flow fields and derived quantities (wall shear stress) through a model bypass around a stenosed artery obtained on an adapted and standard isotropic mesh are compared. In addition, the mesh generation method is tested on a more complex patient-specific geometry. Values of computed wall shear stress are shown to be close to values obtained on anisotropic Hessian-adapted mesh, demonstrating the computational efficiency of the approach in comparison with adaptation based on error indicators derived from flow field.

A PARAMETRIC NUMERICAL INVESTIGATION OF LOCAL HAEMODYNAMICS IN THE END-TO-SIDE ANASTOMOSIS OF CERVICAL-TO-PETROUS BYPASS BASED ON REAL GEOMETRY OF INTERNAL CAROTID ARTERY

Bypass reconstructed from the cervical segment of internal carotid artery (ICA) to its petrous segment allows high-volume bypass flow without any risk of intracranial dissection. The purpose of this study was to investigate the geometric effect on the end-to-side anastomosis of cervical-to-petrous bypass, for its close relationship to local haemodynamic and the long-term performance of bypass. We focused on two controllable geometric parameters: diameter ratio (Φ) and angle (α) between the graft and host arteries. Different models covering a range of Φ (0.75, 1 and 1.25) and α (30°, 45°, 60° and 90°) were constructed based on real geometry of human ICA. Numerical simulations of blood flow were performed in physiological flow condition. The flow patterns, flow distributions, time-average wall shear stress (TAWSS) and oscillatory shear index (OSI) in different models were compared. Our results showed geometric factors have influence on both the local haemodynamic parameters and the flow velocity through downstream branches. Of models with different geometric parameters, the model with Φ ≥ 1 or α = 45° were the most optimized considering haemodynamic performance.

Flexible composition and execution of high performance, high fidelity multiscale biomedical simulations

Multiscale simulations are essential in the biomedical domain to accurately model human physiology. We present a modular approach for designing, constructing and executing multiscale simulations on a wide range of resources, from laptops to petascale supercomputers, including combinations of these. Our work features two multiscale applications, in-stent restenosis and cerebrovascular bloodflow, which combine multiple existing single-scale applications to create a multiscale simulation. These applications can be efficiently coupled, deployed and executed on computers up to the largest (peta) scale, incurring a coupling overhead of 1-10% of the total execution time.

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.

In silico vascular modeling for personalized nanoparticle delivery

AIMS

To predict the deposition of nanoparticles in a patient-specific arterial tree as a function of the vascular architecture, flow conditions, receptor surface density and nanoparticle properties.

MATERIALS & METHODS

The patient-specific vascular geometry is reconstructed from computed tomography angiography images. The isogeometric analysis framework integrated with a special boundary condition for the firm wall adhesion of nanoparticles is implemented. A parallel plate flow chamber system is used to validate the computational model in vitro.

RESULTS

Particle adhesion is dramatically affected by changes in patient-specific attributes, such as branching angle and receptor density. The adhesion pattern correlates well with the spatial and temporal distribution of the wall shear rates. For the case considered, the larger (2.0 µm) particles adhere two-times more in the lower branches of the arterial tree, whereas the smaller (0.5 µm) particles deposit more in the upper branches.

CONCLUSION

Our computational framework in conjunction with patient-specific attributes can be used to rationally select nanoparticle properties to personalize, and thus optimize, therapeutic interventions.

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.

Comparison between mechanical properties of human saphenous vein and umbilical vein

Background

As a main cause of mortality in developed countries, Coronary Artery Disease (CAD) is known as silent killer with a considerable cost to be dedicated for its treatment. Coronary Artery Bypass Graft (CABG) is a common remedy for CAD for which different blood vessels are used as a detour. There is a lack of knowledge about mechanical properties of human blood vessels used for CABG, and while these properties have a great impact on long-term patency of a CABG. Thus, studying these properties, especially those of human umbilical veins which have not been considered yet, looks utterly necessary.

Methods

Umbilical vein, as well as human Saphenous vein, are respectively obtained after cesarean and CABG. First, histological tests were performed to investigate different fiber contents of the samples. Having prepared samples carefully, force-displacement results of samples were rendered to real stress–strain measurements and then a fourth-order polynomial was used to prove the non-linear behavior of these two vessels.

Results

Results were analyzed in two directions, i.e. circumferentially and longitudinally, which then were compared with each other. The comparison between stiffness and elasticity of these veins showed that Saphenous vein’s stiffness is much higher than that of umbilical vein and also, it is less stretchable. Furthermore, for both vessels, longitudinal stiffness was higher than that of circumferential and in stark contrast, stretch ratio in circumferential direction came much higher than longitudinal orientation.

Conclusion

Blood pressure is very high in the region of aorta, so there should be a stiff blood vessel in this area and previous investigations showed that stiffer vessels would have a better influence on the flow of bypass. To this end, the current study has made an attempt to compare these two blood vessels’ stiffness, finding that Saphenous vein is stiffer than umbilical vein which is somehow as stiff as rat aortic vessels. As blood vessel’s stiffness is directly related to elastin and mainly collagen content, results showed the lower amount of these two contents in umbilical vein regarding Saphenous vein.

Longitudinal assessment of hemodynamic endpoints in predicting arteriovenous fistula maturation

Arteriovenous fistula (AVF) nonmaturation is currently a significant clinical problem; however, the mechanisms responsible for this have remained unanswered. Previous work by our group and others has suggested that anatomical configuration and the corresponding hemodynamic endpoints could have an important role in AVF remodeling. Thus, our goal was to assess the longitudinal (temporal) effect of wall shear stress (WSS) on remodeling process of AVFs with two different configurations. The hypothesis is that early assessment of hemodynamic endpoints such as temporal gradient of WSS will predict the maturation status of AVF at later time points. Two AVFs with curved (C-AVF) and straight (S-AVF) configurations were created between the femoral artery and vein of each pig. Three pigs were considered in this study and in total six AVFs (three C-AVF and three S-AVF) were created. The CT scan and ultrasound were utilized to numerically evaluate local WSS at 20 cross-sections along the venous segment of AVFs at 2D (D: days), 7D, and 28D postsurgery. These cross-sections were located at 1.5 mm increments from the anastomosis junction. Local WSS values at these cross-sections were correlated with their corresponding luminal area over time. The WSS in C-AVF decreased from 22.3 ± 4.8 dyn/cm(2) at 2D to 4.1 ± 5.1 dyn/cm(2) at 28D, while WSS increased in S-AVF from 13.0 ± 5.0 dyn/cm(2) at 2D to 36.7 ± 5.3 dyn/cm(2) at 28D. Corresponding to these changes in WSS levels, luminal area of C-AVF dilated (0.23 ± 0.14 cm(2) at 2D to 0.87 ± 0.14 cm(2) at 28D) with attendant increase in flow rate. However, S-AVF had minimal changes in area (0.26 ± 0.02 cm(2) at 2D to 0.27 ± 0.03 cm(2) at 28D) despite some increase in flow rate. Our results suggest that the temporal changes of WSS could have significant effects on AVF maturation. Reduction in WSS over time (regardless of initial values) may result in dilation (p < 0.05), while increase in WSS may be detrimental to maturation. Thus, creation of AVFs in a specific configuration which results in a decline in WSS over time may reduce AVF maturation failure.

Patient-specific multiscale modeling of blood flow for coronary artery bypass graft surgery

We present a computational framework for multiscale modeling and simulation of blood flow in coronary artery bypass graft (CABG) patients. Using this framework, only CT and non-invasive clinical measurements are required without the need to assume pressure and/or flow waveforms in the coronaries and we can capture global circulatory dynamics. We demonstrate this methodology in a case study of a patient with multiple CABGs. A patient-specific model of the blood vessels is constructed from CT image data to include the aorta, aortic branch vessels (brachiocephalic artery and carotids), the coronary arteries and multiple bypass grafts. The rest of the circulatory system is modeled using a lumped parameter network (LPN) 0 dimensional (0D) system comprised of resistances, capacitors (compliance), inductors (inertance), elastance and diodes (valves) that are tuned to match patient-specific clinical data. A finite element solver is used to compute blood flow and pressure in the 3D (3 dimensional) model, and this solver is implicitly coupled to the 0D LPN code at all inlets and outlets. By systematically parameterizing the graft geometry, we evaluate the influence of graft shape on the local hemodynamics, and global circulatory dynamics. Virtual manipulation of graft geometry is automated using Bezier splines and control points along the pathlines. Using this framework, we quantify wall shear stress, wall shear stress gradients and oscillatory shear index for different surgical geometries. We also compare pressures, flow rates and ventricular pressure-volume loops pre- and post-bypass graft surgery. We observe that PV loops do not change significantly after CABG but that both coronary perfusion and local hemodynamic parameters near the anastomosis region change substantially. Implications for future patient-specific optimization of CABG are discussed.

Comparative Study of Viscoelastic Arterial Wall Models in Nonlinear One-Dimensional Finite Element Simulations of Blood Flow

It is well known that blood vessels exhibit viscoelastic properties, which are modeled in the literature with different mathematical forms and experimental bases. The wide range of existing viscoelastic wall models may produce significantly different blood flow, pressure, and vessel deformation solutions in cardiovascular simulations. In this paper, we present a novel comparative study of two different viscoelastic wall models in nonlinear one-dimensional (1D) simulations of blood flow. The viscoelastic models are from papers by Holenstein et al. in 1980 (model V1) and Valdez-Jasso et al. in 2009 (model V2). The static elastic or zero-frequency responses of both models are chosen to be identical. The nonlinear 1D blood flow equations incorporating wall viscoelasticity are solved using a space-time finite element method and the implementation is verified with the Method of Manufactured Solutions. Simulation results using models V1, V2 and the common static elastic model are compared in three application examples: (i) wave propagation study in an idealized vessel with reflection-free outflow boundary condition; (ii) carotid artery model with nonperiodic boundary conditions; and (iii) subject-specific abdominal aorta model under rest and simulated lower limb exercise conditions. In the wave propagation study the damping and wave speed were largest for model V2 and lowest for the elastic model. In the carotid and abdominal aorta studies the most significant differences between wall models were observed in the hysteresis (pressure-area) loops, which were larger for V2 than V1, indicating that V2 is a more dissipative model. The cross-sectional area oscillations over the cardiac cycle were smaller for the viscoelastic models compared to the elastic model. In the abdominal aorta study, differences between constitutive models were more pronounced under exercise conditions than at rest. Inlet pressure pulse for model V1 was larger than the pulse for V2 and the elastic model in the exercise case. In this paper, we have successfully implemented and verified two viscoelastic wall models in a nonlinear 1D finite element blood flow solver and analyzed differences between these models in various idealized and physiological simulations, including exercise. The computational model of blood flow presented here can be utilized in further studies of the cardiovascular system incorporating viscoelastic wall properties.

Quantitative myocardial perfusion magnetic resonance imaging: the impact of pulsatile flow on contrast agent bolus dispersion

Myocardial blood flow (MBF) can be quantified using T1-weighted first-pass magnetic resonance imaging (MRI) in combination with a tracer-kinetic model, like MMID4. This procedure requires the knowledge of an arterial input function which is usually estimated from the left ventricle (LV). Dispersion of the contrast agent bolus may occur between the LV and the tissue of interest. The aim of this study was to investigate the dispersion under conditions of physiological pulsatile blood flow, and to simulate its effect on MBF quantification. The dispersion was simulated in coronary arteries using a computational fluid dynamics (CFD) approach. Simulations were accomplished on straight vessels with stenosis of different degrees and shapes. The results show that dispersion is more pronounced under resting conditions than during hyperemia. Stenosis leads to a reduction of dispersion. In consequence, dispersion results in a systematic MBF underestimation between -0.4% and -9.3%. The relative MBF error depends not only on the dispersion but also on the actual MBF itself. Since MBF under rest is more underestimated than under stress, myocardial perfusion reserve is overestimated between 0.1% and 4.5%. Considering other sources of errors in myocardial perfusion MRI, systematic errors of MBF by bolus dispersion are relatively small.

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.

Clinical Study Protocol for the ARCH Project Computational Modeling for Improvement of Outcome after Vascular Access Creation

Despite clinical guidelines and the possibility of diagnostic vascular imaging, creation and maintenance of a vascular access (VA) remains problematic: avoiding short- and long-term VA dysfunction is challenging. Although prognostic factors for VA dysfunction have been identified in previous studies, their potential interplay at a systemic level is disregarded. Consideration of multiple prognostic patient specific factors and their complex interaction using dedicated computational modeling tools might improve outcome after VA creation by enabling a better selection of VA configuration. These computational modeling tools are developed and validated in the ARCH project: a joint initiative of four medical centers and three industrial partners (FP7-ICT-224390). This paper reports the rationale behind computational modeling and presents the clinical study protocol designed for calibrating and validating these modeling tools. The clinical study is based on the pre-operative collection of structural and functional data at a vascular level, as well as a VA functional evaluation during the follow-up period. The strategy adopted to perform the study and for data collection is also described here.

A simple fluid-structure coupling algorithm for the study of the anastomotic mechanics of vascular grafts

Vascular anastomoses constitute a main factor in poor graft performance due to mismatches in distensibility between the host artery and the graft. This work aims at computational fluid-structure investigations of proximal and distal anastomoses of vein grafts and synthetic grafts. Finite element and finite volume models were developed and coupled with a user-defined algorithm. Emphasis was placed on the simplicity of the coupling algorithm. An artery and vein graft showed a larger dilation mismatch than an artery and synthetic graft. The vein graft distended nearly twice as much as the artery while the synthetic graft displayed only approximately half the arterial dilation. For the vein graft, luminal mismatching was aggravated by development of an anastomotic pseudo-stenosis. While this study focused on end-to-end anastomoses as a vehicle for developing the coupling algorithm, it may serve as useful point of departure for further investigations such as other anastomotic configurations, refined modelling of sutures and fully transient behaviour.

Constitutive modelling of arteries

This review article is concerned with the mathematical modelling of the mechanical properties of the soft biological tissues that constitute the walls of arteries. Many important aspects of the mechanical behaviour of arterial tissue can be treated on the basis of elasticity theory, and the focus of the article is therefore on the constitutive modelling of the anisotropic and highly nonlinear elastic properties of the artery wall. The discussion focuses primarily on developments over the last decade based on the theory of deformation invariants, in particular invariants that in part capture structural aspects of the tissue, specifically the orientation of collagen fibres, the dispersion in the orientation, and the associated anisotropy of the material properties. The main features of the relevant theory are summarized briefly and particular forms of the elastic strain-energy function are discussed and then applied to an artery considered as a thick-walled circular cylindrical tube in order to illustrate its extension–inflation behaviour. The wide range of applications of the constitutive modelling framework to artery walls in both health and disease and to the other fibrous soft tissues is discussed in detail. Since the main modelling effort in the literature has been on the passive response of arteries, this is also the concern of the major part of this article. A section is nevertheless devoted to reviewing the limited literature within the continuum mechanics framework on the active response of artery walls, i.e. the mechanical behaviour associated with the activation of smooth muscle, a very important but also very challenging topic that requires substantial further development. A final section provides a brief summary of the current state of arterial wall mechanical modelling and points to key areas that need further modelling effort in order to improve understanding of the biomechanics and mechanobiology of arteries and other soft tissues, from the molecular, to the cellular, tissue and organ levels.

ECM-based materials in cardiovascular applications: Inherent healing potential and augmentation of native regenerative processes

The in vivo healing process of vascular grafts involves the interaction of many contributing factors. The ability of vascular grafts to provide an environment which allows successful accomplishment of this process is extremely difficult. Poor endothelisation, inflammation, infection, occlusion, thrombosis, hyperplasia and pseudoaneurysms are common issues with synthetic grafts in vivo. Advanced materials composed of decellularised extracellular matrices (ECM) have been shown to promote the healing process via modulation of the host immune response, resistance to bacterial infections, allowing re-innervation and reestablishing homeostasis in the healing region. The physiological balance within the newly developed vascular tissue is maintained via the recreation of correct biorheology and mechanotransduction factors including host immune response, infection control, homing and the attraction of progenitor cells and infiltration by host tissue. Here, we review the progress in this tissue engineering approach, the enhancement potential of ECM materials and future prospects to reach the clinical environment.

Numerical investigation of regurgitation phenomena in pulmonary arteries of Tetralogy of Fallot patients after repair

Pulmonary regurgitation is a very common phenomenon in pulmonary arteries after repair of patients of Tetralogy of Fallot (TOF) which is the most common complex congenital heart diseases. The aim of this study is to use numerical approaches to simulate flow variations in pulmonary artery after repair of patients of TOF. We analyze the flow patterns in an in-vitro bifurcation pulmonary artery and consider effects of various regurgitation fractions (RF or b/f) in left pulmonary artery (LPA) and right pulmonary artery (RPA). We not only observe the variation of flow patterns, but also analyze the results of b/f and net volumetric flow rates in LPA and RPA. In general, the b/f of LPA is higher than RPA in the measured data provided by phase-contrast magnetic resonance imaging (PC-MRI). We validate the result using numerical approaches to analyze the flow patterns in pulmonary artery in this study. The results will be useful for medical doctors when they perform operations for TOF patients.

Neointimal hyperplasia associated with synthetic hemodialysis grafts

Stenosis is a major cause of failure of hemodialysis vascular grafts and is primarily caused by neointimal hyperplasia (NH) at the anastomoses. The objective of this article is to provide a scientific review of the biology underlying this disorder and a critical review of the state-of-the-art investigational preventive strategies in order to stimulate further research in this exciting area. The histology of the NH shows myofibroblasts (that are probably derived from adventitial fibroblasts), extracellular matrices, pro-inflammatory cells including foreign-body giant cells, a variety of growth factors and cytokines, and neovasculature. The contributing factors of the pathogenesis of NH include surgical trauma, bioincompatibility of the synthetic graft, and the various mechanical stresses that result from luminal hypertension and compliance mismatch between the vessel wall and graft. These mechanical stimuli are focal in nature and may have a significant influence on the preferential localization of the NH. Novel mechanical graft designs and local drug delivery strategies show promise in animal models in preventing graft NH development. Successful prevention of graft stenosis would provide a superior alternative to the native fistula as hemodialysis vascular access.

Computational modeling of microarterial anastomoses with size discrepancy (small-to-large)

OBJECTIVE

Among causes of anastomotic failure in microvascular surgery is vessel size mismatch. Where the option of an end-to-side anastomosis is unavailable, an end-to-end construct must be used. Several end-to-end techniques are described to deal with size mismatch. The aim of this study was to numerically model arterial flow patterns and wall shear stresses in four idealized end-to-end anastomoses, where the upstream or recipient artery is smaller. The four techniques modeled were: an invaginating anastomosis; a fish-mouth incision of the smaller vessel; an oblique section of the smaller vessel; and a wedge excision of the larger vessel.

MATERIALS AND METHODS

Flow rate in the right femoral artery of a single outbred male Wistar rat was recorded by transit time ultrasound. Initially, upstream vessel diameter in the models was set at 1 mm, and downstream at 2 mm. The wedge technique was further modeled using a shorter wedge, and using a downstream vessel diameter of 3 mm. Walls were deemed noncompliant. Flow was modeled by the finite volume method using the commercially available computational fluid dynamics code Fluent (Fluent Inc., Lebanon, NH; http://www.fluent.com).

RESULTS

Ring vortices were seen in the invagination and fish-mouth models and showed similar characteristics, although they were less pronounced in the fish-mouth model. The oblique section model demonstrated complex, spiral, counter-rotating vortices that dissipated downstream. Flow separation was least in the first wedge model, with centralization of flow during high but decelerating flow rate. Shortening the wedge length or increasing the downstream vessel diameter to 3 mm led to flow separation. Wall shear stresses were broadly similar for all constructs.

CONCLUSION

Of those modeled, excision of a wedge of the larger vessel proved the best construct. Where a vessel diameter ratio is 1:2, wedge length should be twice the diameter of the larger vessel. A vessel ratio of 1:3 leads to flow separation when using the wedge technique.

Numerical modelling of simulated blood flow in idealized composite arterial coronary grafts: Transient flow

In composite arterial coronary grafts (CACGs), transport phenomena and geometry may considerably alter blood flow dynamics. CACGs aim at revascularizing pathological arteries according to the human anatomy. However, the exact mechanisms causing the failure of coronary bypass grafting are not yet well elucidated. In the present study, computational fluid dynamics (CFD) techniques are applied for the simulation of multi-branched CACGs under physiologically realistic inflow waveforms. The numerical solution is obtained by a finite-volume method formulated in non-orthogonal, curvilinear coordinates and a multi-grid approach. The geometrical models, consisting of idealized and rigid vessels, include the typical T- and a rather new pi-graft configuration. The stenotic effect is also investigated by comparing computational results for three different degrees of area constriction, namely 25%, 50% and 75%, as well as the case without stenosis. Different grafting distances and various inflow rate ratios are imposed, to give an insight into haemodynamical alterations of CACGs and to study the process of restenosis. The results focus on the interaction between the grafts and coronary flows in terms of spatial and temporal variations of velocity and wall shear stress (WSS) distribution. Prominent variations among the different geometries, concerning the velocity profiles and secondary flow motion, are shown. Moreover, the residual flow emerging from different degrees of area constriction shows that low and oscillating shear stresses may arise for even moderate stenotic fields.

Doppler derived quantitative flow estimate in coronary artery bypass graft: a computational multiscale model for the evaluation of the current clinical procedure

In order to investigate the reliability of the so called mean velocity/vessel area formula adopted in clinical practice for the estimation of the flow rate using an intravascular Doppler guide wire instrumentation, a multiscale computational model was used to give detailed predictions on flow profiles within Y-shaped coronary artery bypass graft (CABG) models. At this purpose three CABG models were built from clinical patient's data and used to evaluate and compare, in each model, the computed flow rate and the flow rate estimated according to the assumption of parabolic velocity profile. A consistent difference between the exact and the estimated value of the flow rate was found in every branch of all the graft models. In this study we showed that this discrepancy in the flow rate estimation is coherent to the theory of Womersley regarding spatial velocity profiles in unsteady flow conditions. In particular this work put in evidence that the error in flow rate estimation can be reduced by using the estimation formula recently proposed by Ponzini et al. [Ponzini R, Vergara C, Redaelli A, Veneziani A. Reliable CFD-based estimation of flow rate in haemodynamics measures. Ultrasound Med Biol 2006;32(10):1545-55], accounting for the unsteady nature of blood, applicable in the clinical practice without resorting to further measurements.

A procedure to simulate coronary artery bypass graft surgery

In coronary artery bypass graft (CABG) surgery the involved tissues are overstretched, which may lead to intimal hyperplasia and graft failure. We propose a computational methodology for the simulation of traditional CABG surgery, and analyze the effect of two clinically relevant parameters on the artery and graft responses, i.e., incision length and insertion angle for a given graft diameter. The computational structural analyses are based on actual three-dimensional vessel dimensions of a human coronary artery and a human saphenous vein. The analyses consider the structure of the end-to-side anastomosis, the residual stresses and the typical anisotropic and nonlinear vessel behaviors. The coronary artery is modeled as a three-layer thick-walled tube. The finite element method is employed to predict deformation and stress distribution at various stages of CABG surgery. Small variations of the arterial incision have relatively big effects on the size of the arterial opening, which depends solely on the residual stress state. The incision length has a critical influence on the graft shape and the stress in the graft wall. Stresses at the heel region are higher than those at the toe region. The changes in the mechanical environment are severe along all transitions between the venous tissue and the host artery. Particular stress concentrations occur at the incision ends. The proposed computational methodology may be useful in designing a coronary anastomotic device for reducing surgical trauma. It may improve the quantitative knowledge of vessel diseases and serve as a tool for virtual planning of vascular surgery.

Efficiency differences in computational simulations of the total cavo-pulmonary circulation with and without compliant vessel walls

The Fontan operation is a palliative surgical procedure performed on children born with congenital defects of the heart that have yielded only a single functioning ventricle. The total cavo-pulmonary connection (TCPC) is the most popular variant of the Fontan procedure. The objective of the study was to quantify and compare the efficiency of numerical models of the TCPC with rigid versus elastic vessel wall models. The pressure drop and power loss through both type TCPC models was measured. Significant differences in efficiencies exist between rigid versus elastic numerical models. We have shown incorporating elasticity into numerical models of the total cavo-pulmonary connection is important when determining circuit efficiencies.