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

Automatic Laplacian-based shape optimization for patient-specific vascular grafts

Cognitional heart disease is one of the leading causes of mortality among newborns. Tissue-engineered vascular grafts offer the potential to help treat cognitional heart disease through patient-specific vascular grafts. However, current methods often rely on non-personalized designs or involve significant human intervention. This paper presents a computational framework for the automatic shape optimization of patient-specific tissue-engineered vascular grafts for repairing the aortic arch, aimed at reducing the need for manual input and improving current treatment outcomes, which either use non-patient-specific geometry or require extensive human intervention to design the vascular graft. The paper's core innovation lies in an automatic shape optimization pipeline that combines Bayesian optimization techniques with the open-source finite volume solver, OpenFOAM, and a novel graft deformation algorithm. Specifically, our framework begins with Laplacian mode computation and the approximation of a computationally low-cost Gaussian process surrogate model to capture the minimum weighted combination of inlet-outlet pressure drop (PD) and maximum wall shear stress (WSS). Bayesian Optimization then performs a limited number of OpenFOAM simulations to identify the optimal patient-specific shape. We use imaging and flow data obtained from six patients diagnosed with cognitional heart disease to evaluate our approach. Our results showcase the potential of online training and hemodynamic surrogate model optimization for providing optimal graft shapes. These results show how our framework successfully reduces inlet-outlet PD and maximum WSS compared to pre-lofted models that include both the native geometry and human-designed grafts. Furthermore, we compare how the performance of each design optimized under steady-state simulation compares to that design's performance under transient simulation, and to what extent the optimal design remains similar under both conditions. Our findings underscore that the automated designs achieve at least a 16% reduction in blood flow pressure drop in comparison to geometries optimized by humans.

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

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

svMorph: Interactive Geometry-Editing Tools for Virtual Patient-Specific Vascular Anatomies

We propose svMorph, a framework for interactive virtual sculpting of patient-specific vascular anatomic models. Our framework includes three tools for the creation of tortuosity, aneurysms, and stenoses in tubular vascular geometries. These shape edits are performed via geometric operations on the surface mesh and vessel centerline curves of the input model. The tortuosity tool also uses the physics-based Oriented Particles method, coupled with linear blend skinning, to achieve smooth, elastic-like deformations. Our tools can be applied separately or in combination to produce simulation-suitable morphed models. They are also compatible with popular vascular modeling software, such as simvascular. To illustrate our tools, we morph several image-based, patient-specific models to create a range of shape changes and simulate the resulting hemodynamics via three-dimensional, computational fluid dynamics. We also demonstrate the ability to quickly estimate the hemodynamic effects of the shape changes via the automated generation of associated zero-dimensional lumped-parameter models.

Computational Characterization of Mechanical, Hemodynamic, and Surface Interaction Conditions: Role of Protein Adsorption on the Regenerative Response of TEVGs

Currently available small diameter vascular grafts (<6 mm) present several long-term limitations, which has prevented their full clinical implementation. Computational modeling and simulation emerge as tools to study and optimize the rational design of small diameter tissue engineered vascular grafts (TEVG). This study aims to model the correlation between mechanical-hemodynamic-biochemical variables on protein adsorption over TEVG and their regenerative potential. To understand mechanical-hemodynamic variables, two-way Fluid-Structure Interaction (FSI) computational models of novel TEVGs were developed in ANSYS Fluent 2019R3^® and ANSYS Transient Structural^® software. Experimental pulsatile pressure was included as an UDF into the models. TEVG mechanical properties were obtained from tensile strength tests, under the ISO7198:2016, for novel TEVGs. Subsequently, a kinetic model, linked to previously obtained velocity profiles, of the protein-surface interaction between albumin and fibrinogen, and the intima layer of the TEVGs, was implemented in COMSOL Multiphysics 5.3^®. TEVG wall properties appear critical to understand flow and protein adsorption under hemodynamic stimuli. In addition, the kinetic model under flow conditions revealed that size and concentration are the main parameters to trigger protein adsorption on TEVGs. The computational models provide a robust platform to study multiparametrically the performance of TEVGs in terms of protein adsorption and their regenerative potential.

Low Wall Shear Stress Is Associated with Saphenous Vein Graft Stenosis in Patients with Coronary Artery Bypass Grafting

Biomechanical forces may play a key role in saphenous vein graft (SVG) disease after coronary artery bypass graft (CABG) surgery. Computed tomography angiography (CTA) of 430 post-CABG patients were evaluated and 15 patients were identified with both stenosed and healthy SVGs for paired analysis. The stenosis was virtually removed, and detailed 3D models were reconstructed to perform patient-specific computational fluid dynamic (CFD) simulations. Models were processed to compute anatomic parameters, and hemodynamic parameters such as local and vessel-averaged wall shear stress (WSS), normalized WSS (WSS*), low shear area (LSA), oscillatory shear index (OSI), and flow rate. WSS* was significantly lower in pre-diseased SVG segments compared to corresponding control segments without disease (1.22 vs. 1.73, p = 0.012) and the area under the ROC curve was 0.71. No differences were observed in vessel-averaged anatomic or hemodynamic parameters between pre-stenosed and control whole SVGs. There are currently no clinically available tools to predict SVG failure post-CABG. CFD modeling has the potential to identify high-risk CABG patients who may benefit from more aggressive medical therapy and closer surveillance. Graphical Abstract.

Patient-Specific Hemodynamics of New Coronary Artery Bypass Configurations

PURPOSE

This study aims to quantify the patient-specific hemodynamics of complex conduit routing configurations of coronary artery bypass grafting (CABG) operation which are specifically suitable for off-pump surgeries. Coronary perfusion efficacy and local hemodynamics of multiple left internal mammary artery (LIMA) with sequential and end-to-side anastomosis are investigated. Using a full anatomical model comprised of aortic arch and coronary artery branches the optimum perfusion configuration in multi-vessel coronary artery stenosis is desired.

METHODOLOGY

Two clinically relevant CABG configurations are created using a virtual surgical planning tool where for each configuration set, the stenosis level, anastomosis distance and angle were varied. A non-Newtonian computational fluid dynamics solver in OpenFOAM incorporated with resistance boundary conditions representing the coronary perfusion physiology was developed. The numerical accuracy is verified and results agreed well with a validated commercial cardiovascular flow solver and experiments. For segmental performance analysis, new coronary perfusion indices to quantify deviation from the healthy scenario were introduced.

RESULTS

The first simulation configuration set;-a CABG targeting two stenos sites on the left anterior descending artery (LAD), the LIMA graft was capable of 31 mL/min blood supply for all the parametric cases and uphold the healthy LAD perfusion in agreement with the clinical experience. In the second end-to-side anastomosed graft configuration set;-the radial artery graft anastomosed to LIMA, a maximum of 64 mL/min flow rate in LIMA was observed. However, except LAD, the obtuse marginal (OM) and second marginal artery (m2) suffered poor perfusion. In the first set, average wall shear stress (WSS) were in the range of 4 to 35 dyns/cm² for in LAD. Nevertheless, for second configuration sets the WSS values were higher as the LIMA could not supply enough blood to OM and m2.

CONCLUSION

The virtual surgical configurations have the potential to improve the quality of operation by providing quantitative surgical insight. The degree of stenosis is a critical factor in terms of coronary perfusion and WSS. The sequential anastomosis can be done safely if the anastomosis angle is less than 90 degrees regardless of degree of stenosis. The smaller proposed perfusion index value, O(0.04 - 0) × 10², enable us to quantify the post-op hemodynamic performance by comparing with the ideal healthy physiological flow.

Automatic Shape Optimization of Patient-Specific Tissue Engineered Vascular Grafts for Aortic Coarctation

This paper proposes a computational framework for automatically optimizing the shapes of patient-specific tissue engineered vascular grafts. We demonstrate a proof-of-concept design optimization for aortic coarctation repair. The computational framework consists of three main components including 1) a free-form deformation technique exploring graft geometries, 2) high-fidelity computational fluid dynamics simulations for collecting data on the effects of design parameters on objective function values like energy loss, and 3) employing machine learning methods (Gaussian Processes) to develop a surrogate model for predicting results of high-fidelity simulations. The globally optimal design parameters are then computed by multistart conjugate gradient optimization on the surrogate model. In the experiment, we investigate the correlation among the design parameters and the objective function values. Our results achieve a 30% reduction in blood flow energy loss compared to the original coarctation by optimizing the aortic geometry.

Design and Characterization of a Fluidic Device for the Evaluation of SIS-Based Vascular Grafts

Currently available small diameter vascular conduits present several long-term limitations, which has prevented their full clinical implementation. Commercially available vascular grafts show no regenerative capabilities and eventually require surgical replacement; therefore, it is of great interest to develop alternative regenerative vascular grafts (RVG). Decellularized Small Intestinal Submucosa (SIS) is an attractive material for RVG, however, the evaluation of the performance of these grafts is challenging due to the absence of devices that mimic the conditions found in vivo. Thereby, the objective of this study is to design, manufacture and validate in silico and in vitro, a novel fluidic system for the evaluation of human umbilical vein endothelial cells (HUVECs) proliferation on SIS-based RVG under dynamical conditions. Our perfusion and rotational fluidic system was designed in Autodesk Inventor 2018. In silico Computational Fluid Dynamics (CFD) validation of the system was carried out using Ansys Fluent software from ANSYS, Inc for dynamical conditions of a pulsatile pressure function measured experimentally over a rigid wall model. Mechanical and biological parameters such as flow regime, pressure gradient, wall shear stress (WSS), sterility and indirect cell viability (MTT assay) were also evaluated. Cell adhesion was confirmed by SEM imaging. The fluid flow regime within the system remains laminar. The system maintained sterility and showed low cytotoxicity levels. HUVECs were successfully cultured on SIS-based RVG under both perfusion and rotation conditions. In silico analysis agreed well with our experimental and theoretical results, and with recent in vitro and in vivo reports for WSS. The system presented is a tool for evaluating RVG and represents an alternative to develop new methods and protocols for a more comprehensive study of regenerative cardiovascular devices.

Effect of modified Blalock-Taussig shunt anastomosis angle and pulmonary artery diameter on pulmonary flow

OBJECTIVE

This study aimed to identify the best graft-to-pulmonary artery (PA) anastomosis angle measuring pulmonary blood flow, wall shear stress (WSS), and shunt flow.

METHODS

A tetralogy of Fallot with pulmonary atresia computer model was used to study three different modified Blalock-Taussig shunt (mBTS) anastomosis angle configurations with three different PA diameter configurations. Velocity and WSS were analyzed, and the flow rates at the right PA (RPA) and left PA (LPA) were calculated.

RESULTS

A 4-mm and 8-mm diameter of RPA and LPA, respectively with vertical shunt angle produces the highest total flow. In the RPA larger diameter than the LPA configutations, the left-leaning shunt produces the lowest total PA flow whereas in the LPA larger diameter than the RPA configuratios, the right-leaning shunt produces the lowest total PA flow. Therefore, the shunt anastomosis should not be leaned through the narrow side of PA to reach best flow. As the flow inside the shunt increased, WSS also increased due to enhanced velocity gradients.

CONCLUSION

The anastomosis angle between the conduit and PA affects the flow to PA. Vertical anastomosis configurations increase the total PA flow; thus, these configurations are preferable than the leaned configurations.

The effect of inlet and outlet boundary conditions in image-based CFD modeling of aortic flow

Background

Computational modeling of cardiovascular flow is a growing and useful field, but such simulations usually require the researcher to guess the flow’s inlet and outlet conditions since they are difficult and expensive to measure. It is critical to determine the amount of uncertainty introduced by these assumptions in order to evaluate the degree to which cardiovascular flow simulations are accurate. Our work begins to address this question by examining the sensitivity of flow to several different assumed velocity inlet and outlet conditions in a patient-specific aorta model.

Methods

We examined the differences between plug flow, parabolic flow, linear shear flows, skewed cubic flow profiles, and Womersley flow at the inlet. Only the shape of the inlet velocity profile was varied—all other parameters were identical among these simulations. Secondary flow in the form of a counter-rotating pair of vortices was also added to parabolic axial flow to study its effect on the solution. In addition, we examined the differences between two-element Windkessel, three element Windkessel and the outflow boundary conditions. In these simulations, only the outlet boundary condition was varied.

Results

The results show axial and in-plane velocities are considerably different close to the inlet for the cases with different inlet velocity profile shapes. However, the solutions are qualitatively similar beyond 1.75D, where D is the inlet diameter. This trend is also observed in other quantities such as pressure and wall shear stress. Normalized root-mean-square deviation, a measure of axial velocity magnitude differences between the different cases, generally decreases along the streamwise coordinate. The linear shear inlet velocity boundary condition and plug velocity boundary condition solution exhibit the highest time-averaged wall shear stress, approximately \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$8\%$$\end{document}8% higher than the parabolic inlet velocity boundary condition. Upstream of 1D from the inlet, adding secondary flow has a significant impact on temporal wall shear stress distributions. This is especially observable during diastole, when integrated wall shear stress magnitude varies about \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$26\%$$\end{document}26% between simulations with and without secondary flow. The results from the outlet boundary condition study show the Windkessel models differ from the outflow boundary condition by as much as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$18\%$$\end{document}18% in terms of time-averaged wall shear stress. Furthermore, normalized root-mean-square deviation of axial velocity magnitude, a measure of deviation between Windkessel and the outflow boundary condition, increases along the streamwise coordinate indicating larger variations near outlets.

Conclusion

It was found that the selection of inlet velocity conditions significantly affects only the flow region close to the inlet of the aorta. Beyond two diameters distal to the inlet, differences in flow solution are small. Although additional studies must be performed to verify this result, the data suggest that it is important to use patient-specific inlet conditions primarily if the researcher is concerned with the details of the flow very close to the inlet. Similarly, the selection of outlet conditions significantly affects the flow in the vicinity of the outlets. Upstream of five diameters proximal to the outlet, deviations between the outlet boundary conditions examined are insignificant. Although the inlet and outlet conditions only affect the flow significantly in their respective neighborhoods, our study indicates that outlet conditions influence a larger percentage of the solution domain.

Physiological insights of recent clinical diagnostic and therapeutic technologies for cardiovascular diseases

Diagnostic and therapeutic methods for cardiovascular diseases continue to be developed in the 21st century. Clinicians should consider the physiological characteristics of the cardiovascular system to ensure successful diagnosis and treatment. In this review, we focus on the roles of cardiovascular physiology in recent diagnostic and therapeutic technologies for cardiovascular diseases. In the first section, we discuss how to evaluate and utilize left ventricular arterial coupling in the clinical settings. In the second section, we review unique characteristics of pulmonary circulation in the diagnosis and treatment of pulmonary hypertension. In the third section, we discuss physiological and anatomical factors associated with graft patency after coronary artery bypass grafting. In the last section, we discuss the usefulness of mechanical ventricular unloading after acute myocardial infarction. Clinical development of diagnostic methods and therapies for cardiovascular diseases should be based on physiological insights of the cardiovascular system.

Computational Fluid Dynamics and Additive Manufacturing to Diagnose and Treat Cardiovascular Disease

Noninvasive engineering models are now being used for diagnosing and planning the treatment of cardiovascular disease. Techniques in computational modeling and additive manufacturing have matured concurrently, and results from simulations can inform and enable the design and optimization of therapeutic devices and treatment strategies. The emerging synergy between large-scale simulations and 3D printing is having a two-fold benefit: first, 3D printing can be used to validate the complex simulations, and second, the flow models can be used to improve treatment planning for cardiovascular disease. In this review, we summarize and discuss recent methods and findings for leveraging advances in both additive manufacturing and patient-specific computational modeling, with an emphasis on new directions in these fields and remaining open questions.

Tetralogy of Fallot Surgical Repair: Shunt Configurations, Ductus Arteriosus and the Circle of Willis

In this study, hemodynamic performance of three novel shunt configurations that are considered for the surgical repair of tetralogy of Fallot (TOF) disease are investigated in detail. Clinical experience suggests that the shunt location, connecting angle, and its diameter can influence the post-operative physiology and the neurodevelopment of the neonatal patient. An experimentally validated second order computational fluid dynamics (CFD) solver and a parametric neonatal diseased great artery model that incorporates the ductus arteriosus (DA) and the full patient-specific circle of Willis (CoW) are employed. Standard truncated resistance CFD boundary conditions are compared with the full cerebral arterial system, which resulted 21, -13, and 37% difference in flow rate at the brachiocephalic, left carotid, and subclavian arteries, respectively. Flow splits at the aortic arch and cerebral arteries are calculated and found to change with shunt configuration significantly for TOF disease. The central direct shunt (direct shunt) has pulmonary flow 5% higher than central oblique shunt (oblique shunt) and 23% higher than modified Blalock Taussig shunt (RPA shunt) while the DA is closed. Maximum wall shear stress (WSS) in the direct shunt configuration is 9 and 60% higher than that of the oblique and RPA shunts, respectively. Patent DA, significantly eliminated the pulmonary flow control function of the shunt repair. These results suggests that, due to the higher flow rates at the pulmonary arteries, the direct shunt, rather than the central oblique, or right pulmonary artery shunts could be preferred by the surgeon. This extended model introduced new hemodynamic performance indices for the cerebral circulation that can correlate with the post-operative neurodevelopment quality of the patient.

HEMODYNAMICS-BASED LONG-TERM PATENCY OF DIFFERENT SEQUENTIAL GRAFTING: A PATIENT-SPECIFIC MULTI-SCALE STUDY

Background and aims: Sequential grafting is one of the common coronary artery bypass grafting (CABG) surgery. But the influence of the sequential grafting position on hemodynamics and the graft patency is still unclear. Materials and methods: The zero-dimensional/three-dimensional (0D/3D) coupling method was used to finalize the multi-scale simulation of two different sequential grafting models. First, a patient-specific 3D model was reconstructed based on coronary computed tomography angiography (CCTA) images. Two different sequential grafts were implemented on this patient-specific 3D model by using virtual surgery. Thus, two different postoperative 3D models were built. Then, a lumped parameter model (LPM; 0D) was built based on the patient physiological data to simulate the cardiovascular system. Finally, the 0D/3D coupling method was used to perform the numerical simulation by coupling a 0D LPM of the cardiovascular system and the patient-specific 3D models. Moreover, the long-term patency of these two different sequential grafts was discussed in this paper. Results: The coronary flow rate and the graft flow were calculated and illustrated. The instantaneous wave-free ratio (iFR) were calculated. Postoperative iFR values increase to over 0.90 for both sequential grafts. Some hemodynamics parameters were also illustrated, such as wall shear stress (WSS), oscillatory shear index (OSI). The area of low WSS in Model 1 was much less than that in Model 2. Two regions of high OSI exist in Model 2, while only one in Model 1. Conclusions: No significant differences exist on the short-term outcomes of two models. But the long-term patency of Model 2 was worse. The Model 1 may enhance long-term patency of grafting and should be priority when the sequential grafting need to be carried out.

Hemodynamics of patient-specific aorta-pulmonary shunt configurations

Optimal hemodynamics in aorta-pulmonary shunt reconstruction is essential for improved post-operative recovery of the newborn congenital heart disease patient. However, prior to in vivo execution, the prediction of post-operative hemodynamics is extremely challenging due to the interplay of multiple confounding physiological factors. It is hypothesized that the post-operative performance of the surgical shunt can be predicted through computational blood flow simulations that consider patient size, shunt configuration, cardiac output and the complex three-dimensional disease anatomy. Utilizing only the routine patient-specific pre-surgery clinical data sets, we demonstrated an intelligent decision-making process for a real patient having pulmonary artery atresia and ventricular septal defect. For this patient, a total of 12 customized candidate shunt configurations are contemplated and reconstructed virtually using a sketch-based computer-aided anatomical editing tool. Candidate shunt configurations are evaluated based on the parameters that are computed from the flow simulations, which include 3D flow complexity, outlet flow splits, shunt patency, coronary perfusion and energy loss. Our results showed that the modified Blalock-Taussig (mBT) shunt has 12% higher right pulmonary artery (RPA) and 40% lower left pulmonary artery (LPA) flow compared to the central shunt configuration. Also, the RPA flow regime is distinct from the LPA, creating an uneven flow split at the pulmonary arteries. For all three shunt sizes, right mBT innominate and central configurations cause higher pulmonary artery (PA) flow and lower coronary artery pressure than right and left mBT subclavian configurations. While there is a trade-off between energy loss, flow split and coronary artery pressure, overall, the mBT shunts provide sufficient PA perfusion with higher coronary artery pressures and could be preferred for similar patients having PA overflow risk. Central shunts would be preferred otherwise particularly for cases with very low PA overflow risk.

Numerical analysis of Venous External Scaffolding Technology for Saphenous Vein Grafts

This paper presents a method for analyzing and comparing numerically Saphenous Vein Grafts (SVGs) following Coronary Artery Bypass Graft surgery (CABG). The method analyses the flow dynamics inside vein grafts with and without supporting using Venous External Scaffolding Technology (VEST). The numerical method uses patients׳ specific computational fluid dynamics (CFD) methods to characterize the relevant hemodynamic parameters of patients׳ SVGs. The method was used to compare the hemodynamics of six patient׳s specific model and flow conditions of stented and non-stented SVGs, 12 months post-transplantation. The flow parameters used to characterize the grafts׳ hemodynamics include Time Averaged Wall Shear Stress (TAWSS), Oscillatory Shear Index (OSI) and Relative Residence Time (RRT). The effect of stenting was clearly demonstrated by the chosen parameters. SVGs under constriction of VEST were associated with similar spatial average of TAWSS (10.73 vs 10.29 dyn/cm(2)), yet had fewer lesions with low TAWSS, lower OSI (0.041 vs 0.08) and RRT (0.12 vs 0.24), and more uniform flow with less flow discrepancies. In conclusion, the suggested method and parameters well demonstrated the advantage of VEST support. Stenting vein grafts with VEST improved hemodynamic factors which are correlated to graft failure following CABG procedure.

Left atrial ligation alters intracardiac flow patterns and the biomechanical landscape in the chick embryo

BACKGROUND

Hypoplastic left heart syndrome (HLHS) is a major human congenital heart defect that results in single ventricle physiology and high mortality. Clinical data indicate that intracardiac blood flow patterns during cardiac morphogenesis are a significant etiology. We used the left atrial ligation (LAL) model in the chick embryo to test the hypothesis that LAL immediately alters intracardiac flow streams and the biomechanical environment, preceding morphologic and structural defects observed in HLHS.

RESULTS

Using fluorescent dye injections, we found that intracardiac flow patterns from the right common cardinal vein, right vitelline vein, and left vitelline vein were altered immediately following LAL. Furthermore, we quantified a significant ventral shift of the right common cardinal and right vitelline vein flow streams. We developed an in silico model of LAL, which revealed that wall shear stress was reduced at the left atrioventricular canal and left side of the common ventricle.

CONCLUSIONS

Our results demonstrate that intracardiac flow patterns change immediately following LAL, supporting the role of hemodynamics in the progression of HLHS. Sites of reduced WSS revealed by computational modeling are commonly affected in HLHS, suggesting that changes in the biomechanical environment may lead to abnormal growth and remodeling of left heart structures.

Regional variation of wall shear stress in ascending thoracic aortic aneurysms

The development of an ascending thoracic aortic aneurysm is likely caused by excessive hemodynamic loads exerted on the aneurysmal wall. Computational fluid-dynamic analyses were performed on patient-specific ascending thoracic aortic aneurysms obtained from patients with either bicuspid aortic valve or tricuspid aortic valve to evaluate hemodynamic and wall shear parameters, imparting aneurysm enlargement. Results showed an accelerated flow along the outer aortic wall with helical flow in the aneurysm center for bicuspid aortic valve ascending thoracic aortic aneurysms. In a different way, tricuspid aortic valve ascending thoracic aortic aneurysms exhibited normal systolic flow without substantial secondary pattern. Analysis of wall shear parameters evinced a high and locally varying wall shear stress on the outer aortic wall and high temporal oscillations in wall shear stress (oscillatory shear index) on either left or right side of aneurysmal aorta. These findings may explain the asymmetric dilatation typically observed in ascending thoracic aortic aneurysms. Simulations of a hypertensive scenario revealed an increase in wall shear stress upon 44% compared to normal systemic pressure models. Computational fluid-dynamics–based analysis may allow identification of wall shear parameters portending aneurysm dilatation and hence guide preventative intervention.

Computational fluid dynamics simulation to evaluate aortic coarctation gradient with contrast-enhanced CT

Coarctation of aorta (CoA) is a narrowing of the aorta leading to a pressure gradient (ΔP) across the coarctation, increased afterload and reduced peripheral perfusion pressures. Indication to invasive treatment is based on values of maximal (systolic) trans-coarctation ΔP. A computational fluid dynamic (CFD) approach is herein presented for the non-invasive haemodynamic assessment of ΔP across CoA. Patient-specific CFD simulations were created from contrast-enhanced computed tomography (CT) and appropriate flow boundary conditions. Computed ΔP was validated with invasive intravascular trans-CoA pressure measurements. Haemodynamic indices, including pressure loss coefficient (PLc), time-averaged wall shear stress (TAWSS) and oscillatory shear index (OSI), were also quantified. CFD-estimated ΔP values were comparable to the invasive ones. Moreover, the aorta proximal to CoA was exposed to altered TAWSS and OSI suggesting hypertension. PLc was found as a further geometric marker of CoA severity. Finally, CFD-estimated ΔP confirmed a significant reduction after percutaneous balloon dilatation and stenting of the CoA in one patient (e.g. from ΔP∼52 mmHg to ΔP∼3 mmHg). The validation of the ΔP computations with catheterisation measurements suggests that CFD simulation, based on CT-derived anatomical data, is a useful tool to readily quantify CoA severity.

Fontan conversion templates: patient-specific hemodynamic performance of the lateral tunnel versus the intraatrial conduit with fenestration

Intraatrial-conduit Fontan is considered a modification of both extracardiac and lateral-tunnel Fontan. In this study, the patient-specific hemodynamic performance of intraatrial-conduit and lateral-tunnel Fontan with fenestration, considered as conversion templates, was investigated based on the authors' patient cohort. Pulsatile computational fluid dynamics simulations were performed using patient-specific models of intraatrial-conduit and lateral-tunnel Fontan patients. Real-time "simultaneous" inferior and superior vena cava, pulmonary artery, and fenestration flow waveforms were acquired from ultrasound. Multiple hemodynamic performance indices were investigated, with particular focus on evaluation of the pulsatile flow performance. Power loss inside the lateral-tunnel Fontan appeared to be significantly higher than with the intraatrial-conduit Fontan for patient-specific cardiac output and normalized connection size. Inclusion of the 4-mm fenestration at a 0.24 L/min mean flow resulted in a lower cavopulmonary pressure gradient and less time-averaged power loss for both Fontan connections. Flow structures within the intraatrial conduit were notability more uniform than within the lateral tunnel. Hepatic flow majorly favored the left lung in both surgical connections: conversion from lateral-tunnel to intraatrial-conduit Fontan resulted in better hemodynamics with less power loss, a lower pressure gradient, and fewer stagnant flow zones along the conduit. This patient-specific computational case study demonstrated superior hemodynamics of intraatrial-conduit Fontan over those of lateral-tunnel Fontan with or without fenestration and improved performance after conversion of the lateral tunnel to the intraatrial conduit. The geometry-specific effect of the nonuniform hepatic flow distribution may motivate new rationales for the surgical design.

Aortic outflow cannula tip design and orientation impacts cerebral perfusion during pediatric cardiopulmonary bypass procedures

Poor perfusion of the aortic arch is a suspected cause for peri- and post-operative neurological complications associated with cardiopulmonary bypass (CPB). High-speed jets from 8 to 10FR pediatric/neonatal cannulae delivering ~1 L/min of blood can accrue sub-lethal hemolytic damage while also subjecting the aorta to non-physiologic flow conditions that compromise cerebral perfusion. Therefore, we emphasize the importance of cannulation strategy and hypothesize engineering better CPB perfusion through a redesigned aortic cannula tip. This study employs computational fluid dynamics to investigate novel diffuser-tipped aortic cannulae for shape sensitivity to cerebral perfusion, in an in silico cross-clamped aortic arch model modeled with fixed outflow resistances. 17 parametrically altered configurations of an 8FR end-hole and several diffuser cone angled tips in combination with jet incidence angles toward or away from the head-neck vessels were studied. Experimental pressure-flow characterizations were also conducted on these cannula tip designs. An 8FR end-hole aortic cannula delivering 1 L/min along the transverse aortic arch was found to give rise to backflow from the brachicephalic artery (BCA), irrespective of angular orientation, for the chosen ascending aortic insertion location. Parametric alteration of the cannula tip to include a diffuser cone angle (tested up to 7°) eliminated BCA backflow for any tested angle of jet incidence. Experiments revealed that a 1 cm long 10° diffuser cone tip demonstrated the best pressure-flow performance improvement in contrast with either an end-hole tip or diffuser cone angles greater than 10°. Performance further improved when the diffuser was preceded by an expanded four-lobe swirl inducer attachment-a novel component. In conclusion, aortic cannula orientation is crucial in determining net head-neck perfusion but precise angulations and insertion-depths are difficult to achieve practically. Altering the cannula tip to include a diffuser cone angle has been shown for the first time to have potential in ensuring a net positive outflow at the BCA. Cannula insertion distanced from the BCA inlet may also avoid backflow owing to the Venturi effect, but the diffuser tipped cannula design presents a promising solution to mitigate this issue irrespective of in vivo cannula tip orientation.

Critical transitions in early embryonic aortic arch patterning and hemodynamics

Transformation from the bilaterally symmetric embryonic aortic arches to the mature great vessels is a complex morphogenetic process, requiring both vasculogenic and angiogenic mechanisms. Early aortic arch development occurs simultaneously with rapid changes in pulsatile blood flow, ventricular function, and downstream impedance in both invertebrate and vertebrate species. These dynamic biomechanical environmental landscapes provide critical epigenetic cues for vascular growth and remodeling. In our previous work, we examined hemodynamic loading and aortic arch growth in the chick embryo at Hamburger-Hamilton stages 18 and 24. We provided the first quantitative correlation between wall shear stress (WSS) and aortic arch diameter in the developing embryo, and observed that these two stages contained different aortic arch patterns with no inter-embryo variation. In the present study, we investigate these biomechanical events in the intermediate stage 21 to determine insights into this critical transition. We performed fluorescent dye microinjections to identify aortic arch patterns and measured diameters using both injection recordings and high-resolution optical coherence tomography. Flow and WSS were quantified with 3D computational fluid dynamics (CFD). Dye injections revealed that the transition in aortic arch pattern is not a uniform process and multiple configurations were documented at stage 21. CFD analysis showed that WSS is substantially elevated compared to both the previous (stage 18) and subsequent (stage 24) developmental time-points. These results demonstrate that acute increases in WSS are followed by a period of vascular remodeling to restore normative hemodynamic loading. Fluctuations in blood flow are one possible mechanism that impacts the timing of events such as aortic arch regression and generation, leading to the variable configurations at stage 21. Aortic arch variations noted during normal rapid vascular remodeling at stage 21 identify a temporal window of increased vulnerability to aberrant aortic arch morphogenesis with the potential for profound effects on subsequent cardiovascular morphogenesis.

Hemodynamic characteristics of the vertebrobasilar system analyzed using MRI-based models

The vertebrobasilar system (VBS) is unique in human anatomy in that two arteries merge into a single vessel, and it is especially important because it supplies the posterior circulation of the brain. Atherosclerosis develops in this region, and atherosclerotic plaques in the vertebrobasilar confluence can progress with catastrophic consequences, including artery occlusion. Quantitative assessments of the flow characteristics in the VBS could elucidate the factors that influence flow patterns in this confluence, and deviations from normal patterns might then be used to predict locations to monitor for potential pathological changes, to detect early signs of disease, and to evaluate treatment options and efficacy. In this study, high-field MRI was used in conjunction with computational fluid dynamics (CFD) modeling to investigate the hemodynamics of subject-specific confluence models (n = 5) and to identify different geometrical classes of vertebrobasilar systems (n = 12) of healthy adult subjects. The curvature of the vessels and their mutual orientation significantly affected flow parameters in the VBS. The basilar artery geometry strongly influenced both skewing of the velocity profiles and the wall shear stress distributions in the VBS. All five subjects modeled possessed varying degrees of vertebral asymmetry, and helical flow was observed in four cases, suggesting that factors other than vertebral asymmetry influence mixing of the vertebral artery flow contributions. These preliminary studies verify that quantitative, MR imaging techniques in conjunction with subject-specific CFD models of healthy adult subjects may be used to characterize VBS hemodynamics and to predict flow features that have been related to the initiation and development of atherosclerosis in large arteries. This work represents an important first step towards applying this approach to study disease initiation and progression in the VBS.

Presurgical evaluation of Fontan connection options for patients with apicocaval juxtaposition using computational fluid dynamics

Apicocaval juxtaposition (ACJ) is a rare congenital heart defect associated with single ventricle physiology where optimal positioning of the Fontan conduit for completion of total cavopulmonary connection (TCPC) is still controversial. In ACJ, the cardiac apex is ipsilateral with the inferior vena cava (IVC), risking kinking and collapse of the Fontan conduit at the apex of the heart. The purpose of this study is to evaluate two viable routes for Fontan conduit connection in patients with ACJ, using computational fluid dynamics. Internal energy loss evaluations were used to determine contribution of conduit curvature to the energy efficiency of each cavopulmonary anastomosis configuration. This percentage of energy loss contribution was found to be greater in the case of a curved extracardiac conduit connection (44%, 4.1 mW) traveling behind the ventricular apex, connecting the IVC to the left pulmonary artery, than the straighter lateral tunnel conduit (6%, 1.4 mW) installed through the ventricular apex. In contrast, net energy loss across the anastomosis was significantly lower with extracardiac TCPC (9.3 mW) in comparison with lateral tunnel TCPC (23.2 mW), highlighting that a curved Fontan conduit is favorable provided that it is traded off for a superior cavopulmonary connection efficiency. Therefore, a relatively longer and curved Fontan conduit has been demonstrated to be a suitable connection option independent of anatomical situations.

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.

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.

Mechanotransduction in embryonic vascular development

A plethora of biochemical signals provides spatial and temporal cues that carefully orchestrate the complex process of vertebrate embryonic development. The embryonic vasculature develops not only in the context of these biochemical cues, but also in the context of the biomechanical forces imparted by blood flow. In the mature vasculature, different blood flow regimes induce distinct genetic programs, and significant progress has been made toward understanding how these forces are perceived by endothelial cells and transduced into biochemical signals. However, it cannot be assumed that paradigms that govern the mature vasculature are pertinent to the developing embryonic vasculature. The embryonic vasculature can respond to the mechanical forces of blood flow, and these responses are critical in vascular remodeling, certain aspects of sprouting angiogenesis, and maintenance of arterial-venous identity. Here, we review data regarding mechanistic aspects of endothelial cell mechanotransduction, with a focus on the response to shear stress, and elaborate upon the multifarious effects of shear stress on the embryonic vasculature. In addition, we discuss emerging predictive vascular growth models and highlight the prospect of combining signaling pathway information with computational modeling. We assert that correlation of precise measurements of hemodynamic parameters with effects on endothelial cell gene expression and cell behavior is required for fully understanding how blood flow-induced loading governs normal vascular development and shapes congenital cardiovascular abnormalities.

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.

Computational hemodynamic optimization predicts dominant aortic arch selection is driven by embryonic outflow tract orientation in the chick embryo

In the early embryo, a series of symmetric, paired vessels, the aortic arches, surround the foregut and distribute cardiac output to the growing embryo and fetus. During embryonic development, the arch vessels undergo large-scale asymmetric morphogenesis to form species-specific adult great vessel patterns. These transformations occur within a dynamic biomechanical environment, which can play an important role in the development of normal arch configurations or the aberrant arch morphologies associated with congenital cardiac defects. Arrested migration and rotation of the embryonic outflow tract during late stages of cardiac looping has been shown to produce both outflow tract and several arch abnormalities. Here, we investigate how changes in flow distribution due to a perturbation in the angular orientation of the embryonic outflow tract impact the morphogenesis and growth of the aortic arches. Using a combination of in vivo arch morphometry with fluorescent dye injection and hemodynamics-driven bioengineering optimization-based vascular growth modeling, we demonstrate that outflow tract orientation significantly changes during development and that the associated changes in hemodynamic load can dramatically influence downstream aortic arch patterning. Optimization reveals that balancing energy expenditure with diffusive capacity leads to multiple arch vessel patterns as seen in the embryo, while minimizing energy alone led to the single arch configuration seen in the mature arch of aorta. Our model further shows the critical importance of the orientation of the outflow tract in dictating morphogenesis to the adult single arch and accurately predicts arch IV as the dominant mature arch of aorta. These results support the hypothesis that abnormal positioning of the outflow tract during early cardiac morphogenesis may lead to congenital defects of the great vessels due to altered hemodynamic loading.