A Simulation Protocol for Exercise Physiology in Fontan Patients Using a Closed Loop Lumped-Parameter Model

Background: Reduced exercise capacity is nearly universal among Fontan patients, though its etiology is not yet fully understood. While previous computational studies have attempted to model Fontan exercise, they did not fully account for global physiologic mechanisms nor directly compare results against clinical and physiologic data. Methods: In this study, we developed a protocol to simulate Fontan lower-body exercise using a closed-loop lumped-parameter model describing the entire circulation. We analyzed clinical exercise data from a cohort of Fontan patients, incorporated previous clinical findings from literature, quantified a comprehensive list of physiological changes during exercise, translated them into a computational model of the Fontan circulation, and designed a general protocol to model Fontan exercise behavior. Using inputs of patient weight, height, and if available, patient-specific reference heart rate (HR) and oxygen consumption, this protocol enables the derivation of a full set of parameters necessary to model a typical Fontan patient of a given body-size over a range of physiologic exercise levels. Results: In light of previous literature data and clinical knowledge, the model successfully produced realistic trends in physiological parameters with exercise level. Applying this method retrospectively to a set of clinical Fontan exercise data, direct comparison between simulation results and clinical data demonstrated that the model successfully reproduced the average exercise response of a cohort of typical Fontan patients. Conclusion: This work is intended to offer a foundation for future advances in modeling Fontan exercise, highlight the needs in clinical data collection, and provide clinicians with quantitative reference exercise physiologies for Fontan patients

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

ObjectivesThe objective of this work is to evaluate the hemodynamic performance of a new Y-graft modification of the extracardiac conduit Fontan operation. The performance of the Y-graft design is compared to two designs used in current practice: a t-junction connection of the venae cavae and an offset between the inferior and superior venae cavae.MethodsThe proposed design replaces the current tube grafts used to connect the inferior vena cava to the pulmonary arteries with a Y-shaped graft. Y-graft hemodynamics were evaluated at rest and during exercise with a patient-specific model from magnetic resonance imaging data together with computational fluid dynamics. Four clinically motivated performance measures were examined: Fontan pressures, energy efficiency, inferior vena cava flow distribution, and wall shear stress. Two variants of the Y-graft were evaluated: an “off-the-shelf” graft with 9-mm branches and an “area-preserving” graft with 12-mm branches.ResultsEnergy efficiency of the 12-mm Y-graft was higher than all other models at rest and during exercise, and the reduction in efficiency from rest to exercise was improved by 38%. Both Y-graft designs reduced superior vena cava pressures during exercise by as much as 5 mm Hg. The Y-graft more equally distributed the inferior vena cava flow to both lungs, whereas the offset design skewed 70% of the flow to the left lung. The 12-mm graft resulted in slightly larger regions of low wall shear stress than other models; however, minimum shear stress values were similar.ConclusionsThe area-preserving 12-mm Y-graft is a promising modification of the Fontan procedure that should be clinically evaluated. Further work is needed to correlate our performance metrics with clinical outcomes, including exercise intolerance, incidence of protein-losing enteropathy, and thrombus formation

Outflow boundary conditions for 3D simulations of non-periodic blood flow and pressure fields in deformable arteries

The simulation of blood flow and pressure in arteries requires outflow boundary conditions that incorporate models of downstream domains. We previously described a coupled multidomain method to couple analytical models of the downstream domains with 3D numerical models of the upstream vasculature. This prior work either included pure resistance boundary conditions or impedance boundary conditions based on assumed periodicity of the solution. However, flow and pressure in arteries are not necessarily periodic in time due to heart rate variability, respiration, complex transitional flow or acute physiological changes. We present herein an approach for prescribing lumped parameter outflow boundary conditions that accommodate transient phenomena. We have applied this method to compute haemodynamic quantities in different physiologically relevant cardiovascular models, including patient-specific examples, to study non-periodic flow phenomena often observed in normal subjects and in patients with acquired or congenital cardiovascular disease. The relevance of using boundary conditions that accommodate transient phenomena compared with boundary conditions that assume periodicity of the solution is discussed

Computational Fluid Dynamics Proof of Concept and Analysis of a Self-Powered Fontan Circulation

The Fontan circulation is a result of the last (third stage) surgical procedure to correct a single ventricle congenital cardiac disorder in children. Although the Fontan circulation has been successfully established in surgeries over the years, it is flawed and can lead in certain cases to pre-mature death. The main cause of this failure is due to increased pulmonary vascular resistance due to loss pulse pressure and blood flow. In healthy circulations, the heart pumps directly to the lungs, where as Single Ventricle patients must use a single sided heart to supply blood to the rest of the body before the lungs. Improvements to the Fontan circulation have been proposed, but they require extensive care or external devices. We propose a Self-Powered Fontan circulation that will inject energy into the pulmonary system by adding an injection jet shunt (IJS) directly from the heart. The IJS will provide the pulse pressure, blood flow, and entrainment that the pulmonary vascular system needs to function at a healthy level. The difference between a healthy and sick Fontan circulation is 3-5[mmHg] in the IVC. The goal of the IJS is to cause this 3-5[mmHg] pressure drop in the IVC. In the analysis of the Fontan, ascertaining energy losses due to flow jet impingements and flow mixing is critical. Moreover, in order to better understand surgical alternatives is it important to have a robust multi-scale 0D-3D CFD analysis tool that permits investigation of surgical alternatives in a virtual physics-based environment. To this end, a lumped parameter model (LPM) is tightly coupled at the time step level with a full 3D computational fluid dynamics (CFD) model. Using this model scheme, the Fontan test section is no longer being modeled by the LPM. Therefore, it is not limited by the 0D nature of the vascular resistance, capacitance, and inertia bed model. The CFD can take over at the area of interest which accounts for flow directionality and momentum transfer that the LPM is unable to capture. To efficiently calculate optimal IJS configurations, a closed loop steady state model was created to solve a simplified Fontan circulation in 3D. Three models were created with several different optimized configurations, a synthetic model (average dimensions of 2-4 year-old Fontan patients), and two patient-specific models (10 and 24-year-old). The model configurations include changes in the IJS nozzle diameter and IJS placement along the pulmonary artery. These configurations are compared to a baseline model with no IJS. All three models suggest that the IJS helps to decrease IVC pressure while increasing pulse pressure and blood flow to the pulmonary system

Computational Simulations for Aortic Coarctation: Representative Results From a Sampling of Patients

Treatments for coarctation of the aorta (CoA) can alleviate blood pressure (BP) gradients(D), but long-term morbidity still exists that can be explained by altered indices of hemodynamics and biomechanics. We introduce a technique to increase our understanding of these indices for CoA under resting and nonresting conditions, quantify their contribution to morbidity, and evaluate treatment options. Patient-specific computational fluid dynamics (CFD) models were created from imaging and BP data for one normal and four CoA patients (moderate native CoA: D12 mmHg, severe native CoA: D25 mmHg and postoperative end-to-end and end-to-side patients: D0 mmHg). Simulations incorporated vessel deformation, downstream vascular resistance and compliance. Indices including cyclic strain, time-averaged wall shear stress (TAWSS), and oscillatory shear index (OSI) were quantified. Simulations replicated resting BP and blood flow data. BP during simulated exercise for the normal patient matched reported values. Greatest exercise-induced increases in systolic BP and mean and peak DBP occurred for the moderate native CoA patient (SBP: 115 to 154 mmHg; mean and peak DBP: 31 and 73 mmHg). Cyclic strain was elevated proximal to the coarctation for native CoA patients, but reduced throughout the aorta after treatment. A greater percentage of vessels was exposed to subnormal TAWSS or elevated OSI for CoA patients. Local patterns of these indices reported to correlate with atherosclerosis in normal patients were accentuated by CoA. These results apply CFD to a range of CoA patients for the first time and provide the foundation for future progress in this area

Cavopulmonary assist for the failing Fontan circulation: impact of ventricular function on mechanical support strategy

Mechanical circulatory support--either ventricular assist device (VAD, left-sided systemic support) or cavopulmonary assist device (CPAD, right-sided support)--has been suggested as treatment for Fontan failure. The selection of left- versus right-sided support for failing Fontan has not been previously defined. Computer simulation and mock circulation models of pediatric Fontan patients (15-25 kg) with diastolic, systolic, and combined systolic and diastolic dysfunction were developed. The global circulatory response to assisted Fontan flow using VAD (HeartWare HVAD, Miami Lakes, FL) support, CPAD (Viscous Impeller Pump, Indianapolis, IN) support, and combined VAD and CPAD support was evaluated. Cavopulmonary assist improves failing Fontan circulation during diastolic dysfunction but preserved systolic function. In the presence of systolic dysfunction and elevated ventricular end-diastolic pressure (VEDP), VAD support augments cardiac output and diminishes VEDP, while increased preload with cavopulmonary assist may worsen circulatory status. Fontan circulation can be stabilized to biventricular values with modest cavopulmonary assist during diastolic dysfunction. Systemic VAD support may be preferable to maintain systemic output during systolic dysfunction. Both systemic and cavopulmonary support may provide best outcome during combined systolic and diastolic dysfunction. These findings may be useful to guide clinical cavopulmonary assist strategies in failing Fontan circulations

In-Vitro and In-Silico Investigations of Alternative Surgical Techniques for Single Ventricular Disease

Single ventricle (SV) anomalies account for one-fourth of all cases of congenital Heart disease. The conventional second and third stage i.e. Comprehensive stage II and Fontan procedure of the existing three-staged surgical approach serving as a palliative treatment for this anomaly, entails multiple complications and achieves a survival rate of 50%. Hence, to reduce the morbidity and mortality rate associated with the second and third stages of the existing palliative procedure, the novel alternative techniques called “Hybrid Comprehensive Stage II” (HCSII), and a “Self-powered Fontan circulation” have been proposed. The goal of this research is to conduct in-vitro investigations to validate computational and clinical findings on these proposed novel surgical techniques. The research involves the development of a benchtop study of HCSII and self-powered Fontan circulation

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\u27s 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