Evaluation of a novel Y-shaped extracardiac Fontan baffle using computational fluid dynamics

Computational fluid dynamics models and congenital heart diseases

Mathematical modeling is a powerful tool to investigate hemodynamics of the circulatory system. With improving imaging techniques and detailed clinical investigations, it is now possible to construct patient-specific models of reconstructive surgeries for the treatment of congenital heart diseases. These models can help clinicians to better understand the hemodynamic behavior of different surgical options for a treated patient. This review outlines recent advances in mathematical modeling in congenital heart diseases, the discoveries and limitations these models present, and future directions that are on the horizon.

An integrated approach to patient-specific predictive modeling for single ventricle heart palliation

In patients with congenital heart disease and a single ventricle (SV), ventricular support of the circulation is inadequate, and staged palliative surgery (usually 3 stages) is needed for treatment. In the various palliative surgical stages individual differences in the circulation are important and patient-specific surgical planning is ideal. In this study, an integrated approach between clinicians and engineers has been developed, based on patient-specific multi-scale models, and is here applied to predict stage 2 surgical outcomes. This approach involves four distinct steps: (1) collection of pre-operative clinical data from a patient presenting for SV palliation, (2) construction of the pre-operative model, (3) creation of feasible virtual surgical options which couple a three-dimensional model of the surgical anatomy with a lumped parameter model (LPM) of the remainder of the circulation and (4) performance of post-operative simulations to aid clinical decision making. The pre-operative model is described, agreeing well with clinical flow tracings and mean pressures. Two surgical options (bi-directional Glenn and hemi-Fontan operations) are virtually performed and coupled to the pre-operative LPM, with the hemodynamics of both options reported. Results are validated against postoperative clinical data. Ultimately, this work represents the first patient-specific predictive modeling of stage 2 palliation using virtual surgery and closed-loop multi-scale modeling.

Predictive modeling of the virtual Hemi-Fontan operation for second stage single ventricle palliation: two patient-specific cases

Single ventricle hearts are congenital cardiovascular defects in which the heart has only one functional pumping chamber. The treatment for these conditions typically requires a three-staged operative process where Stage 1 is typically achieved by a shunt between the systemic and pulmonary arteries, and Stage 2 by connecting the superior venous return to the pulmonary circulation. Surgically, the Stage 2 circulation can be achieved through a procedure called the Hemi-Fontan, which reconstructs the right atrium and pulmonary artery to allow for an enlarged confluence with the superior vena cava. Based on pre-operative data obtained from two patients prior to Stage 2 surgery, we developed two patient-specific multi-scale computational models, each including the 3D geometrical model of the surgical junction constructed from magnetic resonance imaging, and a closed-loop systemic lumped-parameter network derived from clinical measurements. "Virtual" Hemi-Fontan surgery was performed on the 3D model with guidance from clinical surgeons, and a corresponding multi-scale simulation predicts the patient's post-operative hemodynamic and physiologic conditions. For each patient, a post-operative active scenario with an increase in the heart rate (HR) and a decrease in the pulmonary and systemic vascular resistance (PVR and SVR) was also performed. Results between the baseline and this "active" state were compared to evaluate the hemodynamic and physiologic implications of changing conditions. Simulation results revealed a characteristic swirling vortex in the Hemi-Fontan in both patients, with flow hugging the wall along the SVC to Hemi-Fontan confluence. One patient model had higher levels of swirling, recirculation, and flow stagnation. However, in both models, the power loss within the surgical junction was less than 13% of the total power loss in the pulmonary circulation, and less than 2% of the total ventricular power. This implies little impact of the surgical junction geometry on the SVC pressure, cardiac output, and other systemic parameters. In contrast, varying HR, PVR, and SVR led to significant changes in theses clinically relevant global parameters. Adopting a work-flow of customized virtual planning of the Hemi-Fontan procedure with patient-specific data, this study demonstrates the ability of multi-scale modeling to reproduce patient specific flow conditions under differing physiological states. Results demonstrate that the same operation performed in two different patients can lead to different hemodynamic characteristics, and that modeling can be used to uncover physiologic changes associated with different clinical conditions.

Caval blood flow distribution in patients with Fontan circulation: quantification by using particle traces from 4D flow MR imaging

PURPOSE

To validate the use of particle traces derived from four-dimensional (4D) flow magnetic resonance (MR) imaging to quantify in vivo the caval flow contribution to the pulmonary arteries (PAs) in patients who had been treated with the Fontan procedure.

MATERIALS AND METHODS

The institutional review boards approved this study, and informed consent was obtained. Twelve healthy volunteers and 10 patients with Fontan circulation were evaluated. The particle trace method consists of creating a region of interest (ROI) on a blood vessel, which is used to emit particles with a temporal resolution of approximately 40 msec. The flow distribution, as a percentage, is then estimated by counting the particles arriving to different ROIs. To validate this method, two independent observers used particle traces to calculate the flow contribution of the PA to its branches in volunteers and compared it with the contribution estimated by measuring net forward flow volume (reference method). After the method was validated, caval flow contributions were quantified in patients. Statistical analysis was performed with nonparametric tests and Bland-Altman plots. P < .05 was considered to indicate a significant difference.

RESULTS

Estimation of flow contributions by using particle traces was equivalent to estimation by using the reference method. Mean flow contribution of the PA to the right PA in volunteers was 54% ± 3 (standard deviation) with the reference method versus 54% ± 3 with the particle trace method for observer 1 (P = .4) and 54% ± 4 versus 54% ± 4 for observer 2 (P = .6). In patients with Fontan circulation, 87% ± 13 of the superior vena cava blood flowed to the right PA (range, 63%-100%), whereas 55% ± 19 of the inferior vena cava blood flowed to the left PA (range, 22%-82%).

CONCLUSION

Particle traces derived from 4D flow MR imaging enable in vivo quantification of the caval flow distribution to the PAs in patients with Fontan circulation. This method might allow the identification of patients at risk of developing complications secondary to uneven flow distribution.

SUPPLEMENTAL MATERIAL

http://radiology.rsna.org/lookup/suppl/doi:10.1148/radiol.12120778/-/DC1.

Optimization of a Y-Graft Design for Improved Hepatic Flow Distribution in the Fontan Circulation

Single ventricle heart defects are among the most serious congenital heart diseases, and are uniformly fatal if left untreated. Typically, a three-staged surgical course, consisting of the Norwood, Glenn, and Fontan surgeries is performed, after which the superior vena cava (SVC) and inferior vena cava (IVC) are directly connected to the pulmonary arteries (PA). In an attempt to improve hemodynamic performance and hepatic flow distribution (HFD) of Fontan patients, a novel Y-shaped graft has recently been proposed to replace the traditional tube-shaped extracardiac grafts. Previous studies have demonstrated that the Y-graft is a promising design with the potential to reduce energy loss and improve HFD. However these studies also found suboptimal Y-graft performance in some patient models. The goal of this work is to determine whether performance can be improved in these models through further design optimization. Geometric and hemodynamic factors that influence the HFD have not been sufficiently investigated in previous work, particularly for the Y-graft. In this work, we couple Lagrangian particle tracking to an optimal design framework to study the effects of boundary conditions and geometry on HFD. Specifically, we investigate the potential of using a Y-graft design with unequal branch diameters to improve hepatic distribution under a highly uneven RPA/LPA flow split. As expected, the resulting optimal Y-graft geometry largely depends on the pulmonary flow split for a particular patient. The unequal branch design is demonstrated to be unnecessary under most conditions, as it is possible to achieve the same or better performance with equal-sized branches. Two patient-specific examples show that optimization-derived Y-grafts effectively improve the HFD, compared to initial nonoptimized designs using equal branch diameters. An instance of constrained optimization shows that energy efficiency slightly increases with increasing branch size for the Y-graft, but that a smaller branch size is preferred when a proximal anastomosis is needed to achieve optimal HFD.

Computational fluid dynamics in congenital heart disease

Computational fluid dynamics has been applied to the design, refinement, and assessment of surgical procedures and medical devices. This tool calculates flow patterns and pressure changes within a virtual model of the cardiovascular system. In the field of paediatric cardiac surgery, computational fluid dynamics is being used to elucidate the optimal approach to staged reconstruction of specific defects and study the haemodynamics of the resulting anatomical configurations after reconstructive or palliative surgery. In this paper, we review the techniques and principal findings of computational fluid dynamics studies as applied to a few representative forms of congenital heart disease.

Large eddy simulation of powered Fontan hemodynamics

Children born with univentricular heart disease typically must undergo three open heart surgeries within the first 2-3 years of life to eventually establish the Fontan circulation. In that case the single working ventricle pumps oxygenated blood to the body and blood returns to the lungs flowing passively through the Total Cavopulmonary Connection (TCPC) rather than being actively pumped by a subpulmonary ventricle. The TCPC is a direct surgical connection between the superior and inferior vena cava and the left and right pulmonary arteries. We have postulated that a mechanical pump inserted into this circulation providing a 3-5 mmHg pressure augmentation will reestablish bi-ventricular physiology serving as a bridge-to-recovery, bridge-to-transplant or destination therapy as a "biventricular Fontan" circulation. The Viscous Impeller Pump (VIP) has been proposed by our group as such an assist device. It is situated in the center of the 4-way TCPC intersection and spins pulling blood from the vena cavae and pushing it into the pulmonary arteries. We hypothesized that Large Eddy Simulation (LES) using high-order numerical methods are needed to capture unsteady powered and unpowered Fontan hemodynamics. Inclusion of a mechanical pump into the CFD further complicates matters due to the need to account for rotating machinery. In this study, we focus on predictions from an in-house high-order LES code (WenoHemo(TM)) for unpowered and VIP-powered idealized TCPC hemodynamics with quantitative comparisons to Stereoscopic Particle Imaging Velocimetry (SPIV) measurements. Results are presented for both instantaneous flow structures and statistical data. Simulations show good qualitative and quantitative agreement with measured data.

Preliminary clinical experience with a bifurcated Y-graft Fontan procedure--a feasibility study

OBJECTIVE

Optimizing flow and diminishing power loss in the Fontan circuit can improve hemodynamic efficiency, potentially improving the long-term outcomes. Computerized modeling has predicted improved energetics with a Y-graft Fontan.

METHODS

From August to December 2010, 6 consecutive children underwent completion Fontan (n=3) or Fontan revision (n=3) using a bifurcated polytetrafluoroethylene Y-graft (18×9×9 mm in 2, 20×10×10 mm in 4) connecting the inferior vena cava to the right and left pulmonary arteries with separate graft limbs. The patents underwent magnetic resonance imaging (n=5) or computed tomography (n=1). Computational fluid dynamics assessed Fontan hemodynamics, power loss, and inferior vena cava flow splits to the branch pulmonary arteries. The clinical parameters were compared with those from 12 patients immediately preceding the present series who had undergone a lateral Fontan procedure.

RESULTS

Despite longer crossclamp and bypass times (not statistically significant), the Y-graft Fontan patients had postoperative courses similar to those of the conventional Fontan patients. Other than 2 early readmissions for pleural effusions managed with diuretics, at 6 to 12 months of follow-up (mean, 8 months), all 6 patients had done well. Postoperative flow modeling demonstrated a balanced distribution of inferior vena cava flow to both pulmonary arteries with minimal flow disturbance. Improvements in hemodynamics and efficiency were noted when the Y-graft branches were anastomosed distally and aligned tangentially with the branch pulmonary arteries.

CONCLUSIONS

The present preliminary surgical experience has demonstrated the clinical feasibility of the bifurcated Y-graft Fontan. Computational fluid dynamics showed acceptable hemodynamics with low calculated power losses and a balanced distribution of inferior vena cava flow to the pulmonary arteries as long as the branch grafts were anastomosed distally.

Comparing pre- and post-operative Fontan hemodynamic simulations: implications for the reliability of surgical planning

Virtual modeling of cardiothoracic surgery is a new paradigm that allows for systematic exploration of various operative strategies and uses engineering principles to predict the optimal patient-specific plan. This study investigates the predictive accuracy of such methods for the surgical palliation of single ventricle heart defects. Computational fluid dynamics (CFD)-based surgical planning was used to model the Fontan procedure for four patients prior to surgery. The objective for each was to identify the operative strategy that best distributed hepatic blood flow to the pulmonary arteries. Post-operative magnetic resonance data were acquired to compare (via CFD) the post-operative hemodynamics with predictions. Despite variations in physiologic boundary conditions (e.g., cardiac output, venous flows) and the exact geometry of the surgical baffle, sufficient agreement was observed with respect to hepatic flow distribution (90% confidence interval-14 ± 4.3% difference). There was also good agreement of flow-normalized energetic efficiency predictions (19 ± 4.8% error). The hemodynamic outcomes of prospective patient-specific surgical planning of the Fontan procedure are described for the first time with good quantitative comparisons between preoperatively predicted and postoperative simulations. These results demonstrate that surgical planning can be a useful tool for single ventricle cardiothoracic surgery with the ability to deliver significant clinical impact.

Pulsatile venous waveform quality affects the conduit performance in functional and "failing" Fontan circulations

OBJECTIVE

To investigate the effect of pulsatility of venous flow waveform in the inferior and superior caval vessels on the performance of functional and "failing" Fontan patients based on two primary performance measures - the conduit power loss and the distribution of inferior caval flow (hepatic factors) to the lungs.

METHODS

Doppler angiography flows were acquired from two typical extra-cardiac conduit "failing" Fontan patients, aged 13 and 25 years, with ventricle dysfunction. Using computational fluid dynamics, haemodynamic efficiencies of "failing", functional, and in vitro-generated mechanically assisted venous flow waveforms were evaluated inside an idealised total cavopulmonary connection with a caval offset. To investigate the effect of venous pulsatility alone, cardiac output was normalised to 3 litres per minute in all cases. To quantify the pulsatile behaviour of venous flows, two new performance indices were suggested.

RESULTS

Variations in the pulsatile content of venous waveforms altered the conduit efficiency notably. High-frequency and low-amplitude oscillations lowered the pulsatile component of the power losses in "failing" Fontan flow waveforms. Owing to the offset geometry, hepatic flow distribution depended strongly on the ratio of time-dependent caval flows and the pulsatility content rather than mixing at the junction. "Failing" Fontan flow waveforms exhibited less balanced hepatic flow distribution to lungs.

CONCLUSIONS

The haemodynamic efficiency of single-ventricle circulation depends strongly on the pulsatility of venous flow waveforms. The proposed performance indices can be calculated easily in the clinical setting in efforts to better quantify the energy efficiency of Fontan venous waveforms in pulsatile settings.

The Fontan Connections: Past, Present, and Future

The Fontan procedure is now considered the final common pathway for patients with anatomical or functional single ventricle. These patients initially have their systemic and pulmonary circulations in parallel, supported by one functional ventricular chamber. The ultimate goal with this procedure is to separate the two circulations, to prevent mixing of venous and arterial blood, and to provide adequate tissue oxygenation. The objective of this article is to review the Fontan procedure with its various modifications and refinements since its introduction to clinical practice in 1971, by Fontan and Baudet.

Three-dimensional simulations in Glenn patients: clinically based boundary conditions, hemodynamic results and sensitivity to input data

While many congenital heart defects can be treated without significant long term sequelae, some achieve successful palliation as their definitive endpoints. The single-ventricle defect is one such defect and leaves the child with only one operational ventricle, requiring the systemic and the pulmonary circulations to be placed in series through several operations performed during early childhood. Numerical simulations may be used to investigate these hemodynamic conditions and their relation to post-operative sequelae; however, they rely heavily on boundary condition prescription. In this study, we investigate the impact of hemodynamic input data uncertainties on simulation results. Imaged-based patient-specific models of the multi-branched pulmonary arteries and superior vena cava were built for five cavopulmonary connection (i.e. Glenn) patients. Magnetic resonance imaging and catheterization data were acquired for each patient prior to their Fontan surgery. Inflow and outflow boundary conditions were constructed to match available clinical data and resulted in the development of a framework to incorporate these types of clinical data into patient-specific simulations. Three-dimensional computational fluid dynamics simulations were run and hemodynamic indicators were computed. Power loss was low (and efficiency very high) and a linear correlation was found between power loss and cardiac index among the five patients. Other indicators such as low wall shear stress were considered to better characterize these patients. Flow was complex and oscillatory near the anastomosis, and laminar in the smaller branches. While common trends were seen among patients, results showed differences among patients, especially in the 3D maps, strengthening the importance of patient-specific simulations. A sensitivity analysis was performed to investigate the impact of input data (clinical and modeling) to construct boundary conditions on several indicators. Overall, the sensitivity of the output indicators to the input data was small but non-negligible. The sensitivity of commonly used hemodynamic indicators to compare patients is discussed in this context. Power efficiency was much more sensitive to pressure variation than power loss. To increase the precision of such indicators, mean flow split between right and left lungs needs to be measured with more accuracy with higher priority than refining the model of how the flow is distributed on average among the smaller branches. Although ± 10% flow split imprecision seemed reasonable in terms of patient comparison, this study suggests that the common practice of imposing a right pulmonary artery/left pulmonary artery flow split of 55%/45% when performing patient specific simulations should be avoided. This study constitutes a first step towards understanding the hemodynamic differences between pre- and post Fontan surgery, predicting these differences, and evaluating surgical outcomes based on preoperative data.

Hypoplastic left heart syndrome: current considerations and expectations

In the recent era, no congenital heart defect has undergone a more dramatic change in diagnostic approach, management, and outcomes than hypoplastic left heart syndrome (HLHS). During this time, survival to the age of 5 years (including Fontan) has ranged from 50% to 69%, but current expectations are that 70% of newborns born today with HLHS may reach adulthood. Although the 3-stage treatment approach to HLHS is now well founded, there is significant variation among centers. In this white paper, we present the current state of the art in our understanding and treatment of HLHS during the stages of care: 1) pre-Stage I: fetal and neonatal assessment and management; 2) Stage I: perioperative care, interstage monitoring, and management strategies; 3) Stage II: surgeries; 4) Stage III: Fontan surgery; and 5) long-term follow-up. Issues surrounding the genetics of HLHS, developmental outcomes, and quality of life are addressed in addition to the many other considerations for caring for this group of complex patients.

Virtual Interventions for Image-based Blood Flow Computation

Image-based blood flow computation provides great promise for evaluation of vascular devices and assessment of surgical procedures. However, many previous studies employ idealized arterial and device models or only patient-specific models from the image data after device deployment, since the tools for model construction are unavailable or limited and tedious to use. Moreover, in contrast to retrospective studies from existing data, there is a pressing need for prospective analysis with the goal of surgical planning. Therefore, it is necessary to construct models with deployed devices in a fast, virtual and interactive fashion. The goal of this paper is to develop new geometric methods to deploy stents or stent grafts virtually to patient-specific geometric models constructed from a 3D segmentation of medical images. A triangular surface representing the vessel lumen boundary is extracted from the segmentation. The diseased portion is either clipped and replaced by the surface of a deployed device or rerouted in the case of a bypass graft. For diseased arteries close to bifurcations, bifurcated device models are generated. A method to map a 2D strut pattern on the surface of a device is also presented. We demonstrate three applications of our methods in personalized surgical planning for aortic aneurysms, aortic coarctation, and coronary artery stenosis using blood flow computation. Our approach enables prospective model construction and may help to expand the throughput required by routine clinical uses in the future.

Image-based modeling of hemodynamics in coronary artery aneurysms caused by Kawasaki disease

Kawasaki Disease (KD) is the leading cause of acquired pediatric heart disease. A subset of KD patients develops aneurysms in the coronary arteries, leading to increased risk of thrombosis and myocardial infarction. Currently, there are limited clinical data to guide the management of these patients, and the hemodynamic effects of these aneurysms are unknown. We applied patient-specific modeling to systematically quantify hemodynamics and wall shear stress in coronary arteries with aneurysms caused by KD. We modeled the hemodynamics in the aneurysms using anatomic data obtained by multi-detector computed tomography (CT) in a 10-year-old male subject who suffered KD at age 3 years. The altered hemodynamics were compared to that of a reconstructed normal coronary anatomy using our subject as the model. Computer simulations using a robust finite element framework were used to quantify time-varying shear stresses and particle trajectories in the coronary arteries. We accounted for the cardiac contractility and the microcirculation using physiologic downstream boundary conditions. The presence of aneurysms in the proximal coronary artery leads to flow recirculation, reduced wall shear stress within the aneurysm, and high wall shear stress gradients at the neck of the aneurysm. The wall shear stress in the KD subject (2.95-3.81 dynes/sq cm) was an order of magnitude lower than the normal control model (17.10-27.15 dynes/sq cm). Particle residence times were significantly higher, taking 5 cardiac cycles to fully clear from the aneurysmal regions in the KD subject compared to only 1.3 cardiac cycles from the corresponding regions of the normal model. In this novel quantitative study of hemodynamics in coronary aneurysms caused by KD, we documented markedly abnormal flow patterns that are associated with increased risk of thrombosis. This methodology has the potential to provide further insights into the effects of aneurysms in KD and to help risk stratify patients for appropriate medical and surgical interventions.

Virtual surgeries in patients with congenital heart disease: a multi-scale modelling test case

The objective of this work is to perform a virtual planning of surgical repairs in patients with congenital heart diseases--to test the predictive capability of a closed-loop multi-scale model. As a first step, we reproduced the pre-operative state of a specific patient with a univentricular circulation and a bidirectional cavopulmonary anastomosis (BCPA), starting from the patient's clinical data. Namely, by adopting a closed-loop multi-scale approach, the boundary conditions at the inlet and outlet sections of the three-dimensional model were automatically calculated by a lumped parameter network. Successively, we simulated three alternative surgical designs of the total cavopulmonary connection (TCPC). In particular, a T-junction of the venae cavae to the pulmonary arteries (T-TCPC), a design with an offset between the venae cavae (O-TCPC) and a Y-graft design (Y-TCPC) were compared. A multi-scale closed-loop model consisting of a lumped parameter network representing the whole circulation and a patient-specific three-dimensional finite volume model of the BCPA with detailed pulmonary anatomy was built. The three TCPC alternatives were investigated in terms of energetics and haemodynamics. Effects of exercise were also investigated. Results showed that the pre-operative caval flows should not be used as boundary conditions in post-operative simulations owing to changes in the flow waveforms post-operatively. The multi-scale approach is a possible solution to overcome this incongruence. Power losses of the Y-TCPC were lower than all other TCPC models both at rest and under exercise conditions and it distributed the inferior vena cava flow evenly to both lungs. Further work is needed to correlate results from these simulations with clinical outcomes.

Cavopulmonary assist: (em)powering the univentricular fontan circulation

Since the Fontan/Kreutzer procedure was introduced, evolutionary clinical advances via a staged surgical reconstructive approach have markedly improved outcomes for patients with functional single ventricle. However, significant challenges remain. Early stage mortality risk seems impenetrable. Serious morbidities - construed as immutable consequences of palliation - have hardly been addressed. Late functional status is increasingly linked to pathophysiologic consequences of prior staged procedures. As more single-ventricle patients survive into adulthood, Fontan failure is emerging as an intractable problem for which there is no targeted therapy. Incremental solutions to address these ongoing problems have not had a measurable impact. Therefore, a fundamental reconsideration of the overall approach is reasonable and warranted. The ability to provide a modest pressure boost (2 to 6 mmHg) to existing blood flow at the total cavopulmonary connection can effectively restore more stable biventricular status. This would impact not only treatment of late Fontan failure, but also facilitate early surgical repair. A realistic means to provide such a pressure boost has never been apparent. Recent advances are beginning to unravel the unique challenges that must be addressed to realize this goal, with promise to open single-ventricle palliation to new therapeutic vistas.

A stochastic collocation method for uncertainty quantification and propagation in cardiovascular simulations

Simulations of blood flow in both healthy and diseased vascular models can be used to compute a range of hemodynamic parameters including velocities, time varying wall shear stress, pressure drops, and energy losses. The confidence in the data output from cardiovascular simulations depends directly on our level of certainty in simulation input parameters. In this work, we develop a general set of tools to evaluate the sensitivity of output parameters to input uncertainties in cardiovascular simulations. Uncertainties can arise from boundary conditions, geometrical parameters, or clinical data. These uncertainties result in a range of possible outputs which are quantified using probability density functions (PDFs). The objective is to systemically model the input uncertainties and quantify the confidence in the output of hemodynamic simulations. Input uncertainties are quantified and mapped to the stochastic space using the stochastic collocation technique. We develop an adaptive collocation algorithm for Gauss-Lobatto-Chebyshev grid points that significantly reduces computational cost. This analysis is performed on two idealized problems--an abdominal aortic aneurysm and a carotid artery bifurcation, and one patient specific problem--a Fontan procedure for congenital heart defects. In each case, relevant hemodynamic features are extracted and their uncertainty is quantified. Uncertainty quantification of the hemodynamic simulations is done using (a) stochastic space representations, (b) PDFs, and (c) the confidence intervals for a specified level of confidence in each problem.

Boundary conditions of patient-specific fluid dynamics modelling of cavopulmonary connections: possible adaptation of pulmonary resistances results in a critical issue for a virtual surgical planning

Cavopulmonary connections are surgical procedures used to treat a variety of complex congenital cardiac defects. Virtual pre-operative planning based on in silico patient-specific modelling might become a powerful tool in the surgical decision-making process. For this purpose, three-dimensional models can be easily developed from medical imaging data to investigate individual haemodynamics. However, the definition of patient-specific boundary conditions is still a crucial issue. The present study describes an approach to evaluate the vascular impedance of the right and left lungs on the basis of pre-operative clinical data and numerical simulations. Computational fluid dynamics techniques are applied to a patient with a bidirectional cavopulmonary anastomosis, who later underwent a total cavopulmonary connection (TCPC). Multi-scale models describing the surgical region and the lungs are adopted, while the flow rates measured in the venae cavae are used at the model inlets. Pre-operative and post-operative conditions are investigated; namely, TCPC haemodynamics, which are predicted using patient-specific pre-operative boundary conditions, indicates that the pre-operative balanced lung resistances are not compatible with the TCPC measured flows, suggesting that the pulmonary vascular impedances changed individually after the surgery. These modifications might be the consequence of adaptation to the altered pulmonary blood flows.

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

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

New insights into pacemaker lead-induced venous occlusion: simulation-based investigation of alterations in venous biomechanics

Venous obstruction is a major complication of transvenous pacemaker placement. Despite the increasing use of pacemakers and implantable cardiac defibrillators, a lack of understanding remains with regard to risk factors for the development of device-associated venous obstruction. We hypothesize that computational fluid dynamics simulations can reveal prothrombogenic locations and define thrombosis risk based on patient-specific anatomies. Using anatomic data derived from computed tomography, computer models of the superior vena cava, subclavian, innominate, and internal jugular veins were constructed for three adult patients with transvenous pacemakers. These models were used to perform patient-specific simulations examining blood flow velocity, wall shear stress, and blood pressure, both with and without the presence of the pacing leads. To better quantify stasis, mean exposure time fields were computed from the venous blood flow data. In comparing simulations with leads to those without, evident increases in stasis at locations between the leads and along the surface of the vessels closest to the leads were found. These locations correspond to regions at known risk for thrombosis. This work presents a novel application of computational methods to study blood flow changes induced by pacemaker leads and possible complications such as venous occlusion and thrombosis. This methodology may add to our understanding of the development of lead-induced thrombosis and occlusion in the clinical arena, and enable the development of new strategies to avoid such complications.

Cavopulmonary assist for the univentricular Fontan circulation: von Kármán viscous impeller pump

OBJECTIVE

In a univentricular Fontan circulation, modest augmentation of existing cavopulmonary pressure head (2-5 mm Hg) would reduce systemic venous pressure, increase ventricular filling, and thus substantially improve circulatory status. An ideal means of providing mechanical cavopulmonary support does not exist. We hypothesized that a viscous impeller pump, based on the von Kármán viscous pump principle, is optimal for this role.

METHODS

A 3-dimensional computational model of the total cavopulmonary connection was created. The impeller was represented as a smooth 2-sided conical actuator disk with rotation in the vena caval axis. Flow was modeled under 3 conditions: (1) passive flow with no disc; (2) passive flow with a nonrotating disk, and (3) induced flow with disc rotation (0-5K rpm). Flow patterns and hydraulic performance were examined for each case. Hydraulic performance for a vaned impeller was assessed by measuring pressure increase and induced flow over 0 to 7K rpm in a laboratory mock loop.

RESULTS

A nonrotating actuator disc stabilized cavopulmonary flow, reducing power loss by 88%. Disk rotation (from baseline dynamic flow of 4.4 L/min) resulted in a pressure increase of 0.03 mm Hg. A further increase in pressure of 5 to 20 mm Hg and 0 to 5 L/min flow was obtained with a vaned impeller at 0 to 7K rpm in a laboratory mock loop.

CONCLUSIONS

A single viscous impeller pump stabilizes and augments cavopulmonary flow in 4 directions, in the desired pressure range, without venous pathway obstruction. A viscous impeller pump applies to the existing staged protocol as a temporary bridge-to-recovery or -transplant in established univentricular Fontan circulations and may enable compressed palliation of single ventricle without the need for intermediary surgical staging or use of a systemic-to-pulmonary arterial shunt.

Optimization of inflow waveform phase-difference for minimized total cavopulmonary power loss

The Fontan operation is a palliative surgical procedure performed on children, born with congenital heart defects that have yielded only a single functioning ventricle. The total cavo-pulmonary connection (TCPC) is a common variant of the Fontan procedure, where the superior vena cava (SVC) and inferior vena cava (IVC) are routed directly into the pulmonary arteries (PA). Due to the limited pumping energy available, optimized hemodynamics, in turn, minimized power loss, inside the TCPC pathway is required for the best optimal surgical outcomes. To complement ongoing efforts to optimize the anatomical geometric design of the surgical Fontan templates, here, we focused on the characterization of power loss changes due to the temporal variations in between SVC and IVC flow waveforms. An experimentally validated pulsatile computational fluid dynamics solver is used to quantify the effect of phase-shift between SVC and IVC inflow waveforms and amplitudes on internal energy dissipation. The unsteady hemodynamics of two standard idealized TCPC geometries are presented, incorporating patient-specific real-time PC-MRI flow waveforms of "functional" Fontan patients. The effects of respiration and pulsatility on the internal energy dissipation of the TCPC pathway are analyzed. Optimization of phase-shift between caval flows is shown to lead to lower energy dissipation up to 30% in these idealized models. For physiological patient-specific caval waveforms, the power loss is reduced significantly (up to 11%) by the optimization of all three major harmonics at the same mean pathway flow (3 L/min). Thus, the hemodynamic efficiency of single ventricle circuits is influenced strongly by the caval flow waveform quality, which is regulated through respiratory dependent physiological pathways. The proposed patient-specific waveform optimization protocol may potentially inspire new therapeutic applications to aid postoperative hemodynamics and improve the well being of the Fontan patients.

A new multiparameter approach to computational simulation for Fontan assessment and redesign

INTRODUCTION

Despite an abundance of prior Fontan simulation articles, there have been relatively few clinical advances that are a direct result of computational methods. We address a few key limitations of previous Fontan simulations as a step towards increasing clinical relevance. Previous simulations have been limited in scope because they have primarily focused on a single energy loss parameter. We present a multi-parameter approach to Fontan modeling that establishes a platform for clinical decision making and comprehensive evaluation of proposed interventions.

METHODS

Time-dependent, 3-D blood flow simulations were performed on six patient-specific Fontan models. Key modeling advances include detailed pulmonary anatomy, catheterization-derived pressures, and MRI-derived flow with respiration. The following performance parameters were used to rank patients at rest and simulated exercise from best to worst performing: energy efficiency, inferior and superior vena cava (IVC, SVC) pressures, wall shear stress, and IVC flow distribution.

RESULTS

Simulated pressures were well matched to catheterization data, but low Fontan pressure did not correlate with high efficiency. Efficiency varied from 74% to 96% at rest, and from 63% to 91% with exercise. Distribution of IVC flow ranged from 88%/12% (LPA/RPA) to 53%/47%. A "transcatheter" virtual intervention demonstrates the utility of computation in evaluating interventional strategies, and is shown to result in increased energy efficiency.

CONCLUSIONS

A multiparameter approach demonstrates that each parameter results in a different ranking of Fontan performance. Ranking patients using energy efficiency does not correlate with the ranking using other parameters of presumed clinical importance. As such, current simulation methods that evaluate energy dissipation alone are not sufficient for a comprehensive evaluation of new Fontan designs.

Assisted Fontan procedure: animal and in vitro models and computational fluid dynamics study

Fontan connection with intermittent compression by wrapped latissimus dorsi (LD) was tested in vivo, in vitro and by means of computational fluid dynamics (CFD). Experimental study: LD was conditioned in four pigs for three weeks before Fontan connection by valved conduit wrapped with LD. Mock circuit: Inflatable cuff wrapped around valved conduit provided intermittent external compression, with pressure and flow measured at driving pressure of 8 or 16 mmHg. CFD study: A circuit was tested for possible increase above basal flow (4 l/min) with intermittent external compression. Experimental study: Intermittent conduit compression by LD provided mean 7% decrease of baseline PA pressure, with simultaneous flow increase of 2%. Mock circuit: By raising the driving pressure from 8 to 16 mmHg, the flow increased with baseline PVR (56%) and with elevated PVR (80%). Total pulmonary flow was reduced during intermittent external compression with both baseline and elevated PVR. CFD study: Compression with 13.0 mmHg provided 4.9% increase of total pulmonary flow with substantial increase of the peak flow (92%). In vivo and in vitro, the increased flow produced by compressing a conduit was confounded by the inevitable intermittent flow restriction. Mathematical model using lower pressure for intermittent external compression showed potential for increase in pulmonary flow.

Current progress in patient-specific modeling

We present a survey of recent advancements in the emerging field of patient-specific modeling (PSM). Researchers in this field are currently simulating a wide variety of tissue and organ dynamics to address challenges in various clinical domains. The majority of this research employs three-dimensional, image-based modeling techniques. Recent PSM publications mostly represent feasibility or preliminary validation studies on modeling technologies, and these systems will require further clinical validation and usability testing before they can become a standard of care. We anticipate that with further testing and research, PSM-derived technologies will eventually become valuable, versatile clinical tools.

Recent advances in the application of computational mechanics to the diagnosis and treatment of cardiovascular disease

During the last 30 years, research into the pathogenesis and progression of cardiovascular disease has had to employ a multidisciplinary approach involving a wide range of subject areas, from molecular and cell biology to computational mechanics and experimental solid and fluid mechanics. In general, research was driven by the need to provide answers to questions of critical importance for disease management. Ongoing improvements in the spatial resolution of medical imaging equipment coupled to an exponential growth in the capacity, flexibility and speed of computational techniques have provided a valuable opportunity for numerical simulations and complex experimental techniques to make a contribution to improving the diagnosis and clinical management of many forms of cardiovascular disease. This paper contains a review of recent progress in the numerical simulation of cardiovascular mechanics, focusing on three particular areas: patient-specific modeling and the optimization of surgery in pediatric cardiology, evaluating the risk of rupture in aortic aneurysms, and noninvasive characterization of intraventricular flow in the management of heart failure.

Patient-specific modeling of cardiovascular mechanics

Advances in numerical methods and three-dimensional imaging techniques have enabled the quantification of cardiovascular mechanics in subject-specific anatomic and physiologic models. Patient-specific models are being used to guide cell culture and animal experiments and test hypotheses related to the role of biomechanical factors in vascular diseases. Furthermore, biomechanical models based on noninvasive medical imaging could provide invaluable data on the in vivo service environment where cardiovascular devices are employed and on the effect of the devices on physiologic function. Finally, patient-specific modeling has enabled an entirely new application of cardiovascular mechanics, namely predicting outcomes of alternate therapeutic interventions for individual patients. We review methods to create anatomic and physiologic models, obtain properties, assign boundary conditions, and solve the equations governing blood flow and vessel wall dynamics. Applications of patient-specific models of cardiovascular mechanics are presented, followed by a discussion of the challenges and opportunities that lie ahead.

Characterization of coherent structures in the cardiovascular system

Recent advances in blood flow modeling have provided highly resolved, four-dimensional data of fluid mechanics in large vessels. The motivation for such modeling is often to better understand how flow conditions relate to health and disease, or to evaluate interventions that affect, or are affected by, blood flow mechanics. Vessel geometry and the pulsatile pumping of blood leads to complex flow, which is often difficult to characterize. This article discusses a computational method to better characterize blood flow kinematics. In particular, we compute Lagrangian coherent structures (LCS) to study flow in large vessels. We demonstrate that LCS can be used to characterize flow stagnation, flow separation, partitioning of fluid to downstream vasculature, and mechanisms governing stirring and mixing in vascular models. This perspective allows valuable understanding of flow features in large vessels beyond methods traditionally considered.