冠状动脉移植管的血流动力学数值模拟

采用有限元数值计算的方法，模拟了冠状动脉搭桥术中移植管内的生理流动。计算模型包括了冠状动脉狭窄，并考虑了由于移植管直径大于冠状动脉直径而在两者缝合时移植管的变形。计算结果分析了缝合区附近的流场、二次流、壁面切应力在心动周期内的时空分布情况。计算结果表明，在缝合前端下游，存在一个低切应力、高切应力梯度的区域，在缝合区底部存在一个高切应力、高切应力梯度的区域。这两个区域都是内皮细胞增生并造成移植管术后再阻塞的危险部位。     

Computational haemodynamics analysis and comparison study of arterio-venous grafts

Haemodialysis arterio-venous graft failure is related to the development of stenotic lesions most commonly located near the venous anastomosis, especially in the toe region. 'Disturbed flow' interaction with the vessel wall surface, characterized by haemodynamic parameters based on the local wall shear stress or the radial pressure gradient, have been widely recognized as the trigger mechanism of a cascade of abnormal biological events leading to occlusive developments. Assuming incompressible laminar flow and rigid, in-plane vessel walls, validated haemodynamics are numerically simulated for a constant-diameter end-to-side base case, the Venaflo graft, and an improved graft-end configuration. The geometric design of the new graft-end was based on the reduction of three time- and area-averaged haemodynamic parameters, i.e. the wall shear stress gradient, wall shear stress angle gradient, and radial pressure gradient. Considering the critical toe region, the Venaflo graft has demonstrated measurable improvements over the base case configuration in predictive computer simulations as well as in clinical trials. The performance improvement should be further enhanced with the modifications illustrated by the new design.     

Numerical investigation and identification of susceptible sites of atherosclerotic lesion formation in a complete coronary artery bypass model

As hemodynamics is widely believed to correlate with anastomotic stenosis in coronary bypass surgery, this paper investigates the flow characteristics and distributions of the hemodynamic parameters (HPs) in a coronary bypass model (which includes both proximal and distal anastomoses), under physiological flow conditions. Disturbed flows (flow separation/reattachment, vortical and secondary flows) as well as regions of high oscillatory shear index (OSI) with low wall shear stress (WSS), i.e., high-OSI-and-low-WSS and low-OSI-and-high-WSS were found in the proximal and distal anastomoses, especially at the toe and heel regions of distal anastomosis, which indicate highly suspected sites for the onset of the atherosclerotic lesions. The flow patterns found in the graft and distal anastomoses of our model at deceleration phases are different from those of the isolated distal anastomosis model. In addition, a huge significant difference in segmental averages of HPs was found between the distal and proximal anastomoses. These findings further suggest that intimal hyperplasia would be more prone to form in the distal anastomosis than in the proximal anastomosis, particularly along the suture line at the toe and heel of distal anastomosis.     

Computational haemodynamics analysis and comparison study of arterio-venous grafts

Haemodialysis arterio-venous graft failure is related to the development of stenotic lesions most commonly located near the venous anastomosis, especially in the toe region. 'Disturbed flow' interaction with the vessel wall surface, characterized by haemodynamic parameters based on the local wall shear stress or the radial pressure gradient, have been widely recognized as the trigger mechanism of a cascade of abnormal biological events leading to occlusive developments. Assuming incompressible laminar flow and rigid, in-plane vessel walls, validated haemodynamics are numerically simulated for a constant-diameter end-to-side base case, the Venaflo^TM graft, and an improved graft-end configuration. The geometric design of the new graft-end was based on the reduction of three time- and area-averaged haemodynamic parameters, i.e. the wall shear stress gradient, wall shear stress angle gradient, and radial pressuregradient. Considering the critical toeregion, the Venaflo^TM graft has demonstrated measurable improvements over the base case configuration in predictive computer simulations as well as in clinical trials. The performance improvement should be further enhanced with the modifications illustrated by the new design.     

Influences of graft diameter on the blood flow in 2-way bypassing surgery

The graft diameter plays a critically important role in the long-term patency rates of bypass surgery. To clarify the influence of graft diameter on the blood flows in the femoral 2-way bypass surgery, the physiologically pulsatile flows in two femoral bypass models were simulated with numerical methods. For the sake of comparison, the models were constructed with identical geometry parameters except the different diameters of grafts. Two models with small and large grafts were studied. The boundary conditions for the simulation of blood flow were constant for both models. The maximum Reynolds number was 832.8, and the Womersley number was 6.14. The emphases of results were on the analysis of flow fields in the vicinity of the distal anastomosis. The temporal-spatial distributions of velocity vectors, pressure drop between the proximal and distal toe, wall shear stresses, wall shear stress gradients and oscillating shear index were compared. The present study indicated that femoral artery bypassed with a large graft demonstrated disturbed axial flow and secondary flow at the distal anastomosis while the axial flow at its downstream of toe was featured with larger and more uniform longitudinal velocities. Meanwhile, the large model exhibits less refluences, relatively uniform wall shear stresses, lower pressure and smaller wall shear stress gradients, whereas it does not have any advantages in the distributions of secondary flow and the oscillating shear index. In general, the large model exhibits better and more uniform hemodynamic phenomena near the vessel wall and may be effective in preventing the initiation and development of postoperative intimal hyperplasia and restenosis.     

Numerical simulations of unsteady flows in a stenosed coronary bypass graft

Using the finite element method, physiological blood flows through a three-dimensional model of a coronary graft are numerically analysed. The model includes a stenosis shape in the host artery upstream from the anastomosis. Recirculating areas, secondary flows, wall shear stress (WSS) and spatial wall shear stress gradients (WSSGs) are studied for different flow repartitions and at different times in the cycle. The temporal and spatial evolutions of the recirculating areas downstream from the stenosis, their interactions with the flow issued from the graft and their associated WSSs highlight that the presence of the stenosis in the recipient artery is essential for prediction of the evolution of a grafting at the beginning of its implantation. The areas downstream from the stenosis expansion, non-existent for a host artery without stenosis, are submitted to low and oscillating WSS between −0.5 and 0 Nm⁻². The stagnation point on the recipient artery floor is subjected to high positive and negative WSSG_nd values, and its location is dependent on the residual flow through the stenosis.     

Intimal hyperplasia and hemodynamic factors in arterial bypass and arteriovenous grafts: a review

Stenosis at the graft–vein junction caused by intimal hyperplasia (IH) is the major cause of failure of vascular access grafts used for hemodialysis. There is a strong relationship between hemodynamic factors and formation of IH. The hemodynamic pattern and the location of IH are different in arterial bypass grafts (ABGs) compared with arteriovenous grafts (AVGs). In an ABG, end-to-side anastomosis of the expanded polytetrafluoroethylene graft and artery produces hemodynamic changes around the junction. IH develops at the arterial floor and the toe and heel of the distal anastomosis. Low shear stress and oscillating shear forces at the arterial floor and the heel plus a high wall sheer stress (WSS) gradient at the toe probably promote IH development. Compliance mismatch between the graft and artery causes turbulence that may contribute to IH formation. The blood flow rate in AVGs is 5–10 times greater than that in ABGs. High flow causes turbulence that injures endothelial cells and eventually results in IH. The peak WSS in AVGs is about 6N/m², much higher than that in ABGs. Excessively high WSS may effect IH formation in AVGs. Several venous cuff or patch anastomotic designs have been used in attempts to regulate hemodynamic factors in grafts. In ABGs, these designs appear to help decrease IH formation. In AVGs, however, they generally have not improved patency rates. In a high-flow system such as an AVG, more drastic changes in anastomotic design may be required.     

Influence of Junction Angle Changed Models on Hemodynamic of Coronary Artery Bypass Graft

The restenosis after coronary artery bypass graft(CABG) is attributed to the formation of intimal hyperplasia(IH) at the anastomosis,which is closely related to hemodynamic depend on the geometric model. In order to give a reasonable assessment of the surgery effect and judge the long-term patency rate,the hemodynamic of CABG surgery program is compared with that of surgery design of the junction angle changed.Based on in-vivo CT coronary angiography datasets,the individual geometric model of CABG reconstructed instead of idealized geometric models are applied to simulate the real physiological blood flow utilizing pulsatile flow boundary waveforms in the present study. The simulation results show that the maximum wall shear rate(WSS) value is at the bottom of anastomosis. Moreover,the stagnation zone growing gradually with the greater angle downstream the anastomosis is prone to form the IH,which is consistent with clinical observation. It is proved that the surgery being better suited to maintain graft patency is successful.     

管径比对全阻塞动脉旁路移植流场影响的数值研究

目的 研究管径比对全阻塞动脉旁路移植流场的影响,为指导动脉旁路移植手术,减少术后再狭窄提供理论依据。方法 采用数值方法研究5种不同移植管与主血管的管径比对全阻塞情况下动脉旁路移植流场的影响,分析速度、二次流、壁面切应力和壁面切应力梯度等血流动力学参数的分布及其随管径比增大的改变。同时,为表明本文所采用模型的合理性,针对目前常使用的两类模型,比较在管径比1.0情况下全阻塞完整模型（Model A）、全阻塞局部模型（Model B）和75％狭窄完整模型（Model C）之间的血流动力学差异。结果 Model A和Model C的血流动力学特性是完全不同的;移植管顶部截面内的速度分布对下游吻合处的主血管底部壁面切应力的影响是显著的,最大相差达79％。大管径比时,主血管底部的壁面低切应力区较大,但壁面切应力分布均匀,壁面切应力梯度较小。而小管径比时,主血管底部的壁面低切应力区较小,但壁面切应力梯度较大。结论 采用整体模型单独研究全阻塞情况下的管径比对流场的影响是有必要的。管径比对全阻塞动脉旁路移植的流场具有显著影响,采用大管径比进行动脉旁路移植将有助于缓解吻合口处由于再狭窄而产生的阻塞。     

A Novel Coronary Artery Bypass Graft Design of Sequential Anastomoses

In this paper, the hemodynamics in a three-dimensional out-of-plane sequential bypass graft model is first investigated. Based on the advantageous flow characteristics observed within the side-to-side (STS) anastomosis in the sequential bypass graft simulation, a new CABG coupled-sequential anastomosis configuration is designed, entailing coupled STS and end-to-side (ETS) anastomotic components. In this new CABG design, the flow fields and distributions of various wall shear stress parameters within the STS and ETS anastomotic regions are studied, and compared to those of the conventional distal anastomosis, by means of computational fluid dynamics simulation of pulsatile Newtonian blood flow. Simulation results demonstrate that the new sequential anastomoses model provides: (i) a more uniform and smooth flow at the ETS anastomosis, without any stagnation point on the artery bed and vortex formation in the heel region of the ETS anastomosis within the coronary artery; (ii) a spare route for the blood flow to the coronary artery, to avoid re-operation in case of re-stenosis in either of the anastomoses; and (iii) improved distribution of hemodynamic parameters at the coronary artery bed and in the heel region of the ETS anastomosis, with more moderate shear stress indices. These advantages of the new design over the conventional ETS anastomosis are influenced by the occlusion ratio of the native coronary artery, and are most prominent when the proximal segment of the coronary artery is fully occluded. By varying the design parameters of the anastomotic angle and distance between the two anastomoses, the superior coupled STS–ETS anastomoses design is found to have the anastomotic angle of 30° and 30 mm distance between the two (STS and ETS) components.     

Human Saphenous Vein Coronary Artery Bypass Graft Morphology,Geometry and Hemodynamics

Coronary artery bypass graft (CABG) failure has been linked to graft hemodynamics, in particular wall shear stress. This study characterizes the morphology, geometry and wall shear stress patterns in human CABGs. The intimal thickness (IT) in 49 human saphenous vein CABGs was measured by digital light microscopy. The geometry of six saphenous vein CABGs was replicated by post-mortem infusion of Batsons #17 anatomical corrosion casting compound. Graft hemodynamics were evaluated in two flow models, fabricated from the casts, under steady (Re = 110) and pulsatile flow (Re = 110, = 2) conditions. Saphenous vein CABGs in situ for more than 2 months had, on average, the greatest IT on the hood and suture sites of the distal anastomosis. Floor thickening was highly variable and significantly less than IT at the hood, suture site and graft body. All casts showed an indentation along the floor and 5/6 casts displayed a sharp local curvature on the hood. In both flow models, a large increase in wall shear rate occurred on the hood, just proximal to the toe. The local geometry of the hood created this large spatial gradient in wall shear stress which is a likely factor in hood intimal hyperplasia.     

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.     

In vitro measurements of velocity and wall shear stress in a novel sequential anastomotic graft design model under pulsatile flow conditions

This study documents the superior hemodynamics of a novel coupled sequential anastomoses (SQA) graft design in comparison with the routine conventional end-to-side (ETS) anastomoses in coronary artery bypass grafts (CABG). The flow fields inside three polydimethylsiloxane (PDMS) models of coronary artery bypass grafts, including the coupled SQA graft design, a conventional ETS anastomosis, and a parallel side-to-side (STS) anastomosis, are investigated under pulsatile flow conditions using particle image velocimetry (PIV). The velocity field and distributions of wall shear stress (WSS) in the models are studied and compared with each other. The measurement results and WSS distributions, computed from the near wall velocity gradients reveal that the novel coupled SQA design provides: (i) a uniform and smooth flow at its ETS anastomosis, without any stagnation point on the artery bed and vortex formation in the heel region of the ETS anastomosis within the coronary artery; (ii) more favorable WSS distribution; and (iii) a spare route for the blood flow to the coronary artery, to avoid re-operation in case of re-stenosis in either of the anastomoses. This in vitro investigation complements the previous computational studies of blood flow in this coupled SQA design, and is another necessary step taken toward the clinical application of this novel design. At this point and prior to the clinical adoption of this novel design, in vivo animal trials are warranted, in order to investigate the biological effects and overall performance of this anastomotic configuration in vivo.     

一种可以避免血管再狭窄的双移植管搭桥方式的数值模拟

目的为改善冠状动脉旁路移植管的局部血流动力学,降低血管再狭窄的发生机率,研究一种可以避免血管再狭窄的双移植管搭桥方式。方法利用有限元分析方法,对传统模型和双移植管搭桥模型进行血流动力学模拟仿真,计算缝合区附近的流场、壁面切应力等血流动力学因素的分布情况。结果该双移植管搭桥具有较好的血流动力学分布,明显改善了主搭桥血管与冠状动脉缝合处的血流动力学参数,消除了该部位的涡流和流动停滞点,提高了底面的壁面切应力数值。在辅助搭桥血管与冠状动脉缝合处涡流区长度仅3 mm,与原主搭桥血管缝合处的涡流长度4.5 mm相比明显减小。辅助搭桥管分流了约36%的血液,只有约64%的血液流过了主搭桥管。结论该双移植管搭桥有助于减小内膜增生的发生机率。     

Evaluation of the hemodynamics in straight 6-mm and tapered 6- to 8-mm grafts as upper arm hemodialysis vascular access

The present study is intended to investigate and compare the hemodynamics in two different sizes of hemodialysis arteriovenous grafts for upper arm hemodialysis vascular access: 8-mm tapered to 6-mm at the arterial side and straight 6 mm. A computational simulation approach is presented for this study, which is validated against the available experimental and numerical pressure measurements in the literature. The imposed boundary conditions at the arterial inlet and venous outlet boundaries of the models are physiological velocity and pressure waveforms, respectively. Blood flow fields and distribution patterns of the hemodynamic indices including wall shear stress (WSS) as one of the major hemodynamic parameters of the cardiovascular system and spatial wall shear stress gradient (SWSSG) as an indicator of disturbed flow patterns and hence susceptible sites of lesion developments are analyzed and compared between the two grafts. The tapered 6- to 8-mm graft seemingly is associated with less disturbed flow patterns within the venous anastomosis (VA) and the vein downstream while benefiting from higher blood flow rates within. Also, it shows a definitive advantage in terms of WSS and SWSSG distribution patterns around the VA and throughout the vein downstream with significantly lower values, which reduce the risk of thrombosis formation and stenotic lesion developments. The only disadvantage encountered in using 6- to 8-mm tapered graft is higher values of hemodynamic parameters at the arterial junction attributable to its significantly higher mean blood flow rate within. The results clearly indicate that the tapered 6- to 8-mm graft entirely outperforms straight 6-mm graft hemodynamically as an upper arm hemodialysis vascular access graft and confirms clinical data in the literature, which suggests advantageous use of tapered 6- to 8-mm grafts in the creation of upper arm brachioaxillary hemodialysis vascular access grafts in selected groups of patients with expectably higher patency rates and lower complications.     

Effects of an Artery/Vascular Graft Compliance Mismatch on Protein Transport: A Numerical Study

Small-diameter vascular graft failure by intimal hyperplasia and thrombosis may result from flow disturbances and disruption of chemical transport in the fluid at the distal anastomosis, because of compliance mismatch between the graft and host artery. In previous studies. lower-than-normal wall shear stress (WSS), particle trapping, and high particle residence times were observed at the distal anastomosis due to a pulsatile tubular expansion effect caused by nonuniform radial deformations. This study was undertaken to examine effects of compliance and radius mismatch on the distribution of a model protein released at the graft-fluid interface. Finite element simulations of end-to-end vascular grafting were performed under pulsatile flow, using fluid-structure coupling to give physiologic wall displacements. Results showed that protein is convected smoothly downstream in a uniform compliant tube. A compliance mismatch disturbed the transport, causing positive and negative gradients in the concentration profile at the distal anastomosis. This was seen when the graft and artery radii were matched at zero pressure and at mean arterial pressure; low WSSs were only observed in the former case. Thus the distal intimal hypertrophy seen in noncompliant grafts may be caused partly by decreased WSS, and partly by concentration gradients of dissolved chemicals affecting chemotaxis of cells.     

Flow dynamics across end-to-end vascular bypass graft anastomoses

In this article, a numerical simulation of steady flow across an end-to-end vascular bypass graft anastomosis is presented.In vitro experiments were performed to determine the variations in the conduit cross section at the anastomosis. Penrose surgical drainage tubing was used to simulate an artery and was anastomosed with Polytetrafluoroethylene (PTFE) vascular grafts using a continuous suturing technique. Artery to artery anastomosis was simulated by suturing two Penrose tubing segments. The anastomotic specimens were subject to static transmural pressure in the physiologic range to determine the instantaneous diameter and compliance as a function of the distance from the anastomotic site. The experimentally determined geometries were used to simulate steady flow through an end-to-end anastomosis using the finite analytic (FA) numerical solution technique. The results demonstrated a region of flow separation 2 mm distal to the Penrose tubing-Penrose tubing anastomosis (simulating an artery-artery anastomosis) at higher transmural pressures. Moreover, wall shear stresses increased proximal to the anastomosis in flow from the Penrose tubing to the graft. In flow from the graft to the Penrose tubing, low wall shear stresses were observed distal to the anastomosis. Flow separation was observed distal to the anastomosis at higher transmural pressures with uniform inlet velocity condition. The region of low shear stress in flow from PTFE graft to the Penrose tubing was located nearer to the anastomosis with thin wall grafts than that with standard wall thickness grafts. Our steady flow model studies suggest a correlation between regions of low wall shear stress and the development of anastomotic neointimal fibrous hyperplasia in end-to-end anastomoses.     

单路和双路CABG中血流动力学的比较

为了改善冠状动脉搭桥术后的血流动力学，提出了对称双路搭桥的改进措施。利用有限元分析方法，对冠状动脉搭桥术中单路移植管和对称双路移植管内的生理流动进行了数值模拟，并对两种情况下的血流动力学计算结果进行了比较。计算结果分析了缝合区附近的流场、壁面剪应力等血流动力学因素在心动周期内的时空分布情况。研究结果表明，对称双路搭桥比单路搭桥具有更合理的血流动力学，可以避免动脉粥样硬化的危险性血流动力学因素，从而减少手术再狭窄的发生。     

Hemodynamic parameters and early intimal thickening in branching blood vessels

Intimal thickening due to atherosclerotic lesions or intimal hyperplasia in medium to large blood vessels is a major contributor to heart disease, the leading cause of death in the Western World. Balloon angioplasty with stenting, bypass surgery, and endarterectomy (with or without patch reconstruction) are some of the techniques currently applied to occluded blood vessels. On the basis of the preponderance of clinical evidence that disturbed flow patterns play a key role in the onset and progression of atherosclerosis and intimal hyperplasia, it is of interest to analyze suitable hemodynamic wall parameters that indicate susceptible sites of intimal thickening and/or favorable conditions for thrombi formation. These parameters, based on the wall shear stress, wall pressure, or particle deposition, are applied to interpret experimental/clinical observations of intimal thickening. Utilizing the parameters as "indicator" functions, internal branching blood vessel geometries are analyzed and possibly altered for different purposes: early detection of possibly highly stenosed vessel segments, prediction of future disease progression, and vessel redesign to potentially improve long-term patency rates. At the present time, the focus is on the identification of susceptible sites in branching blood vessels and their subsequent redesign, employing hemodynamic wall parameters. Specifically, the time-averaged wall shear stress (WSS), its spatial gradient (WSSG), the oscillatory shear index (OSI), and the wall shear stress angle gradient (WSSAG) are compared with experimental data for an aortoceliac junction. Then, the OSI, wall particle density (WPD), and WSSAG are segmentally averaged for different carotid artery bifurcations and compared with clinical data of intimal thickening. The third branching blood vessel under consideration is the graft-to-vein anastomosis of a vascular access graft. Suggested redesigns reduce several hemodynamic parameters (i.e., the WSSG, WSSAG, and normal pressure gradient [NPG]), thereby reducing the likelihood of restenosis, especially near the critical toe region.