Computational Fluid Dynamics Turbulence Model and Experimental Study for a Fontan Cavopulmonary Assist Device.

Head-flow HQ curves for a Fontan cavopulmonary assist device (CPAD) were measured using a blood surrogate in a mock circulatory loop and simulated with various computational fluid dynamics (CFD) models. The tests benchmarked the CFD tools for further enhancement of the CPAD design. Recommended Reynolds-Averaged Navier-Stokes (RANS) CFD approaches for the development of conventional ventricular assist devices (VAD) were found to have shortcomings when applied to the Fontan CPAD, which is designed to neutralize off-condition obstruction risks that could contribute to a major adverse event. The no-obstruction condition is achieved with a von Karman pump, utilizing large clearances and small blade heights, which challenge conventional VAD RANS-based CFD hemodynamic simulations. High-fidelity large eddy simulation (LES) is always recommended; however, this may be cost-inhibitive for optimization studies in commercial settings, thus the reliance on RANS models. This study compares head and power predictions of various RANS turbulence models, employing experimental measurements and LES results as a basis for comparison. The models include standard k-ϵ, re-normalization group k-ϵ, realizable k-ϵ, shear stress transport (SST) k-ω, SST with transitional turbulence, and Generalized k-ω. For the pressure head predictions, it was observed that the standard k-ϵ model provided far better agreement with experiment. For the rotor torque, k-ϵ predictions were 30% lower than LES, while the SST and LES torque values were near identical. For the Fontan CPAD, the findings support using LES for the final design simulations, k-ϵ model for head and general flow simulation, and SST for power, shear stress, hemolysis, and thrombogenicity predictions.

Passive performance evaluation and validation of a viscous impeller pump for subpulmonary fontan circulatory support.

Patients with single ventricle defects undergoing the Fontan procedure eventually face Fontan failure. Long-term cavopulmonary assist devices using rotary pump technologies are currently being developed as a subpulmonary power source to prevent and treat Fontan failure. Low hydraulic resistance is a critical safety requirement in the event of pump failure (0 RPM) as a modest 2 mmHg cavopulmonary pressure drop can compromise patient hemodynamics. The goal of this study is therefore to assess the passive performance of a viscous impeller pump (VIP) we are developing for Fontan patients, and validate flow simulations against in-vitro data. Two different blade heights (1.09 mm vs 1.62 mm) and a blank housing model were tested using a mock circulatory loop (MCL) with cardiac output ranging from 3 to 11 L/min. Three-dimensional flow simulations were performed and compared against MCL data. In-silico and MCL results demonstrated a pressure drop of < 2 mmHg at a cardiac output of 7 L/min for both blade heights. There was good agreement between simulation and MCL results for pressure loss (mean difference - 0.23 mmHg 95% CI [0.24-0.71]). Compared to the blank housing model, low wall shear stress area and oscillatory shear index on the pump surface were low, and mean washout times were within 2 s. This study demonstrated the low resistance characteristic of current VIP designs in the failed condition that results in clinically acceptable minimal pressure loss without increased washout time as compared to a blank housing model under normal cardiac output in Fontan patients.

Passive Performance Evaluation and Validation of a Viscous Impeller Pump for Subpulmonary Fontan Circulatory Support.

Patients with single ventricle defects undergoing the Fontan procedure eventually face Fontan failure. Long-term cavopulmonary assist devices using rotary pump technologies are currently being developed as a subpulmonary power source to prevent and treat Fontan failure. Low hydraulic resistance is a critical safety requirement in the event of pump failure (0 RPM) as a modest 2 mmHg cavopulmonary pressure drop can compromise patient hemodynamics. The goal of this study is therefore to assess the passive performance for a viscous impeller pump (VIP) we are developing for Fontan patients, and validate flow simulations against in-vitro data. Two different blade heights (1.09 mm vs 1.62 mm) and a blank housing model were tested using a mock circulatory loop (MCL) with cardiac output ranging from 3 to 11 L/min. Three-dimensional flow simulations were performed and compared against MCL data. In-silico and MCL results demonstrated a clinically insignificant pressure drop of $<$ 2 mmHg at a cardiac output of 7 L/min for both blade heights. There was good agreement between simulation and MCL results for pressure loss (mean difference -0.23 mmHg 95% CI [0.24 -0.71]). Compared to the blank housing model, low wall shear stress area and oscillatory shear index on the pump surface were low, and mean washout times were within 2 seconds. This study demonstrated the low resistance characteristic of current VIP designs in the failed condition that results in clinically acceptable minimal pressure loss with low risk of thrombosis.

Computational Modeling to Support Surgical Decision Making in Single Ventricle Physiology.

Many of the advances in congenital heart surgery were built upon lessons and insights gained from model simulations. While animal and mock-circuit models have historically been the main arena to test new operative techniques and concepts, the recognition that complex cardiovascular anatomy and circulation can be modeled mathematically ushered a new era of collaboration between surgeons and engineers. In 1996, the computational age in congenital heart surgery began when investigators in London and Milan tapped the power of the computer to simulate the Fontan procedure and introduced operative improvements. Since then, computational modeling has led to numerous contributions in congenial heart surgery as continuing sophistication and advances in numerical and imaging methods furthered the ability to refine anatomic and physiologic details. Idealized generic models have given way to precise patient-specific simulations of the 3-dimensional anatomy, reconstructed circulation, affected hemodynamics, and altered physiology. Tools to perform virtual surgery, and predict flow dynamic and circulatory results, have been developed for some of the most complex defects, such as those requiring single ventricle palliation. In today's quest for personalized medicine and precision care, computational modeling's role to assist surgical planning in complex congenital heart surgery will continue to grow and evolve. With ever closer collaboration between surgeons and engineers, and clear understanding of modeling limitations, computational simulations can be a valuable adjunct to support preoperative surgical decision making.

Multiscale Modeling of Superior Cavopulmonary Circulation: Hemi-Fontan and Bidirectional Glenn Are Equivalent.

Superior cavopulmonary circulation (SCPC) can be achieved by either the Hemi-Fontan (hF) or Bidirectional Glenn (bG) connection. Debate remains as to which results in best hemodynamic results. Adopting patient-specific multiscale computational modeling, we examined both the local dynamics and global physiology to determine if surgical choice can lead to different hemodynamic outcomes. Six patients (age: 3-6 months) underwent cardiac magnetic resonance imaging and catheterization prior to SCPC surgery. For each patient: (1) a finite 3-dimensional (3D) volume model of the preoperative anatomy was constructed to include detailed definition of the distal branch pulmonary arteries, (2) virtual hF and bG operations were performed to create 2 SCPC 3D models, and (3) a specific lumped network representing each patient's entire cardiovascular circulation was developed from clinical data. Using a previously validated multiscale algorithm that couples the 3D models with lumped network, both local flow dynamics, that is, power loss, and global systemic physiology can be quantified. In 2 patients whose preoperative imaging demonstrated significant left pulmonary artery (LPA) stenosis, we performed virtual pulmonary arterioplasty to assess its effect. In one patient, the hF model showed higher power loss (107%) than the bG, while in 3, the power losses were higher in the bG models (18-35%). In the remaining 2 patients, the power loss differences were minor. Despite these variations, for all patients, there were no significant differences between the hF and bG models in hemodynamic or physiological outcomes, including cardiac output, superior vena cava pressure, right-left pulmonary flow distribution, and systemic oxygen delivery. In the 2 patients with LPA stenosis, arterioplasty led to better LPA flow (5-8%) while halving the power loss, but without important improvements in SVC pressure or cardiac output. Despite power loss differences, both hF and bG result in similar SCPC hemodynamics and physiology outcome. This suggests that for SCPC, the pre-existing patient-specific physiology and condition, such as pulmonary vascular resistance, are more deterministic in the hemodynamic performance than the type of surgical palliation. Multiscale modeling can be a decision-assist tool to assess whether an extensive LPA reconstruction is needed at the time of SCPC for LPA stenosis.

Cavopulmonary assist: Long-term reversal of the Fontan paradox.

OBJECTIVE: Fontan circulatory inefficiency can be addressed by replacing the missing subpulmonary power source to reverse the Fontan paradox. An implantable cavopulmonary assist device is described that will simultaneously reduce systemic venous pressure and increase pulmonary arterial pressure, improving preload and cardiac output, in a univentricular Fontan circulation on a long-term basis. METHODS: A rotary blood pump that was based on the von Karman viscous pump was designed for implantation into the total cavopulmonary connection (TCPC). It will impart modest pressure energy to augment Fontan flow without risk of obstruction. In the event of rotational failure, it is designed to default to a passive flow diverter. Pressure-flow performance was characterized in vitro in a Fontan mock circulatory loop with blood analog. RESULTS: The pump performed through the fully specified operating range, augmenting flow in all 4 directions of the TCPC. Pressure rise of 6 to 8 mm Hg was readily achieved, ranging to 14 mm Hg at highest speed (5600 rpm). Performance was consistent across a wide range of cardiac outputs. In stalled condition (0 rpm), there was no discernible pressure loss across the TCPC. CONCLUSIONS: A blood pump technology is described that can reverse the Fontan paradox and may permit a surgical strategy of long-term biventricular maintenance of a univentricular Fontan circulation. The technology is intended for Fontan failure in which right-sided circulatory inefficiencies predominate and ventricular systolic function is preserved. It may also apply before clinical Fontan failure as health maintenance to preempt the progression of Fontan disease.

Patient-Specific Multiscale Modeling of the Assisted Bidirectional Glenn.

BACKGROUND: First-stage palliation of neonates with single-ventricle physiology is associated with poor outcomes and challenging clinical management. Prior computational modeling and in vitro experiments introduced the assisted bidirectional Glenn (ABG), which increased pulmonary flow and oxygenation over the bidirectional Glenn (BDG) and the systemic-to-pulmonary shunt in idealized models. In this study, we demonstrate that the ABG achieves similar performance in patient-specific models and assess the influence of varying shunt geometry. METHODS: In a small cohort of single-ventricle prestage 2 patients, we constructed three-dimensional in silico models and tuned lumped parameter networks to match clinical measurements. Each model was modified to produce virtual BDG and ABG surgeries. We simulated the hemodynamics of the stage 1 procedure, BDG, and ABG by using multiscale computational modeling, coupling a finite-element flow solver to the lumped parameter network. Two levels of pulmonary vascular resistances (PVRs) were investigated: baseline (low) PVR of the patients and doubled (high) PVR. The shunt nozzle diameter, anastomosis location, and shape were also manipulated. RESULTS: The ABG increased the pulmonary flow rate and pressure by 15% to 20%, which was accompanied by a rise in superior vena caval pressure (2 to 3 mm Hg) at both PVR values. Pulmonary flow rate and superior vena caval pressures were most sensitive to the shunt nozzle diameter. CONCLUSIONS: Patient-specific ABG performance was similar to prior idealized simulations and experiments, with good performance at lower PVR values in the range of measured clinical data. Larger shunt outlet diameters and lower PVR led to improved ABG performance.

Optimization of the Assisted Bidirectional Glenn Procedure for First Stage Single Ventricle Repair.

BACKGROUND: First-stage single-ventricle palliation is challenging to manage, and significant interstage morbidity and mortality remain. Prior computational and in vitro studies of the assisted bidirectional Glenn (ABG), a novel first-stage procedure that has shown potential for early conversion to a more stable augmented Glenn physiology, demonstrated increased pulmonary flow and oxygen delivery while decreasing cardiac work, as compared to conventional stage-1 alternatives. This study aims to identify optimal shunt designs for the ABG to improve pulmonary flow while maintaining or decreasing superior vena caval (SVC) pressure. METHODS: A representative three-dimensional model of a neonatal bidirectional Glenn (BDG) was created, with a shunt connecting the innominate artery to the SVC. The shunt design was studied as a six-parameter constrained shape optimization problem. We simulated hemodynamics for each candidate designs using a multiscale finite element flow solver and compared performance against designs with taper-less shunts, the standalone BDG, and a simplified control volume model. Three values of pulmonary vascular resistance (PVR) of 2.3, 4.3, and 7.1 WUm2 were studied. RESULTS: Increases in pulmonary flow were generally accompanied by increases in SVC pressure, except at low PVR (2.3 WUm2), where the optimal shunt geometry achieved a 13% increase in pulmonary flow without incurring any increase in SVC pressure. Shunt outlet area was the most influential design parameter, while others had minimal effect. CONCLUSION: Assisted bidirectional Glenn performance is sensitive to PVR and shunt outlet diameter. An increase in pulmonary flow without a corresponding increase in SVC pressure is possible only when PVR is low.

An interactive simulation tool for patient-specific clinical decision support in single-ventricle physiology.

OBJECTIVE: Modeling of single-ventricle circulations has yielded important insights into their unique flow dynamics and physiology. Here we translated a state-of-the-art mathematical model into a patient-specific clinical decision support interactive Web-based simulation tool and show validation for all 3 stages of single-ventricular palliation. METHODS: Via the adoption a validated lumped parameter method, complete cardiovascular-pulmonary circulatory models of all 3 stages of single-ventricle physiology were created within a simulation tool. The closed-loop univentricular heart model includes scaling for growth and respiratory effects, and typical patient-specific parameters are entered through an intuitive user interface. The effects of medical or surgical interventions can be simulated and compared. To validate the simulator, patient parameters were collected from catheterization reports. Four simulator outputs were compared against catheterization findings: pulmonary to systemic flow ratio (Qp:Qs), systemic arterial saturation (SaO2), mean pulmonary arterial pressure (mPAp), and systemic-venous oxygen difference (SaO2-SvO2). RESULTS: Data from 60 reports were used. Compared with the clinical values, the simulator results were not significantly different in mean Qp:Qs, SaO2, or mPAp (P > .09). There was a statistical but clinically insignificant difference in average SaO-SvO2 (average difference 1%, P < .01). Linear regression analyses revealed a good prediction for each variable (Qp:Qs, R2 = 0.79; SaO2, R2 = 0.64; mPAp, R2 = 0.69; SaO2-SvO2, R2 = 0.93). CONCLUSIONS: This simulator responds quickly and predicts patient-specific hemodynamics with good clinical accuracy. By predicting postoperative and postintervention hemodynamics in all 3 stages of single-ventricle physiology, the simulator could assist in clinical decision-making, training, and consultation. Continuing model refinement and validation will further its application to the bedside.

In Vitro Validation of a Multiscale Patient-Specific Norwood Palliation Model.

In Norwood physiology, shunt size and the occurrence of coarctation can affect hemodynamics significantly. The aim of the study was to validate an in vitro model of the Norwood circulation against clinical measurements for patients presenting differing aortic morphologies. The mock circulatory system used coupled a lumped parameter network of the neonatal Norwood circulation with modified Blalock-Taussig (mBT) shunt with a three-dimensional aorta model. Five postoperative aortic arch anatomies of differing morphologies were reconstructed from imaging data, and the system was tuned to patient-specific clinical values. Experimentally measured flow rates and pressures were compared with clinical measurements. Time-based experimental and clinical pressure and flow signals within the aorta and pulmonary circulation branches agreed closely (0.72 < R < 0.95) for the five patients, whereas mean values within the systemic and pulmonary branches showed no significant differences (95% confidence interval). We validated an experimental multiscale model of the Norwood circulation with mBT shunt by showing it capable of reproducing clinical pressure and flow rates at various positions of the circulation with very good fidelity across a range of patient physiologies and morphologies. The multiscale aspect of the model provides a means to study variables in isolation with their effects both locally and at the system level. The model serves as a tool to further the understanding of the complex physiology of single-ventricle circulation.

In Vitro Assessment of the Assisted Bidirectional Glenn Procedure for Stage One Single Ventricle Repair.

This in vitro study compares the hemodynamic performance of the Norwood and the Glenn circulations to assess the performance of a novel assisted bidirectional Glenn (ABG) procedure for stage one single ventricle surgery. In the ABG, the flow in a bidirectional Glenn procedure is assisted by injection of a high-energy flow stream from the systemic circulation using an aorta-caval shunt with nozzle. The aim is to explore experimentally the potential of the ABG as a surgical alternative to current surgical practice. The experiments are directly compared against previously published numerical simulations. A multiscale mock circulatory system was used to measure the hemodynamic performance of the three circulations. For each circulation, the system was tested using both low and high values of pulmonary vascular resistance. Resulting parameters measured were: pressure and flow rate at left/right pulmonary artery and superior vena cava (SVC). Systemic oxygen delivery (OD) was calculated. A parametric study of the ratio of ABG nozzle to shunt diameter was done. We report time-based comparisons with numerical simulations for the three surgical variants tested. The ABG circulation demonstrated an increase of 30-38% in pulmonary flow with a 2-3.7 mmHg increase in SVC pressure compared to the Glenn and a 4-14% higher systemic OD than either the Norwood or the Glenn. The nozzle/shunt diameter ratio affected the local hemodynamics. These experimental results agreed with those of the numerical model: mean flow values were not significantly different (p > 0.05) while mean pressures were comparable within 1.2 mmHg. The results verify the approaches providing two tools to study this complicated circulation. Using a realistic experimental model we demonstrate the performance of a novel surgical procedure with potential to improve patient hemodynamics in early palliation of the univentricular circulation.

A multiscale model for the study of cardiac biomechanics in single-ventricle surgeries: a clinical case.

Complex congenital heart disease characterized by the underdevelopment of one ventricular chamber (single ventricle (SV) circulation) is normally treated with a three-stage surgical repair. This study aims at developing a multiscale computational framework able to couple a patient-specific three-dimensional finite-element model of the SV to a patient-specific lumped parameter (LP) model of the whole circulation, in a closed-loop fashion. A sequential approach was carried out: (i) cardiocirculatory parameters were estimated by using a fully LP model; (ii) ventricular material parameters and unloaded geometry were identified by means of the stand-alone, three-dimensional model of the SV; and (iii) the three-dimensional model of SV was coupled to the LP model of the circulation, thus closing the loop and creating a multiscale model. Once the patient-specific multiscale model was set using pre-operative clinical data, the virtual surgery was performed, and the post-operative conditions were simulated. This approach allows the analysis of local information on ventricular function as well as global parameters of the cardiovascular system. This methodology is generally applicable to patients suffering from SV disease for surgical planning at different stages of treatment. As an example, a clinical case from stage 1 to stage 2 is considered here.

Control of respiration-driven retrograde flow in the subdiaphragmatic venous return of the Fontan circulation.

Respiration influences the subdiaphragmatic venous return in the total cavopulmonary connection (TCPC) of the Fontan circulation whereby both the inferior vena cava (IVC) and hepatic vein flows can experience retrograde motion. Controlling retrograde flows could improve patient outcomes. Using a patient-specific model within a Fontan mock circulatory system with respiration, we inserted a valve into the IVC to examine its effects on local hemodynamics while varying retrograde volumes by changing vascular impedances. A bovine valved conduit reduced IVC retrograde flow to within 3% of antegrade flow in all cases. The valve closed only under conditions supporting retrograde flow and its effects on local hemodynamics increased with larger retrograde volume. Liver and TCPC pressures improved only when the valve leaflets were closed whereas cycle-averaged pressures improved only slightly (<1 mm Hg). Increased pulmonary vascular resistance raised mean circulation pressures, but the valve functioned and cardiac output improved and stabilized. Power loss across the TCPC improved by 12%-15% (p < 0.05) with a valve. The effectiveness of valve therapy is dependent on patient vascular impedance.

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.

Mock circulatory system of the Fontan circulation to study respiration effects on venous flow behavior.

We describe an in vitro model of the Fontan circulation with respiration to study subdiaphragmatic venous flow behavior. The venous and arterial connections of a total cavopulmonary connection (TCPC) test section were coupled with a physical lumped parameter (LP) model of the circulation. Intrathoracic and subdiaphragmatic pressure changes associated with normal breathing were applied. This system was tuned for two patients (5 years, 0.67 m2; 10 years, 1.2 m2) to physiological values. System function was verified by comparison to the analytical model on which it was based and by consistency with published clinical measurements. Overall, subdiaphragmatic venous flow was influenced by respiration. Flow within the arteries and veins increased during inspiration but decreased during expiration, with retrograde flow in the inferior venous territories. System pressures and flows showed close agreement with the analytical LP model (p < 0.05). The ratio of the flow rates occurring during inspiration to expiration were within the clinical range of values reported elsewhere. The approach used to set up and control the model was effective and provided reasonable comparisons with clinical data.

Modeling single ventricle physiology: review of engineering tools to study first stage palliation of hypoplastic left heart syndrome.

First stage palliation of hypoplastic left heart syndrome, i.e., the Norwood operation, results in a complex physiological arrangement, involving different shunting options (modified Blalock-Taussig, RV-PA conduit, central shunt from the ascending aorta) and enlargement of the hypoplastic ascending aorta. Engineering techniques, both computational and experimental, can aid in the understanding of the Norwood physiology and their correct implementation can potentially lead to refinement of the decision-making process, by means of patient-specific simulations. This paper presents some of the available tools that can corroborate clinical evidence by providing detailed insight into the fluid dynamics of the Norwood circulation as well as alternative surgical scenarios (i.e., virtual surgery). Patient-specific anatomies can be manufactured by means of rapid prototyping and such models can be inserted in experimental set-ups (mock circulatory loops) that can provide a valuable source of validation data as well as hydrodynamic information. Such models can be tuned to respond to differing the patient physiologies. Experimental set-ups can also be compatible with visualization techniques, like particle image velocimetry and cardiovascular magnetic resonance, further adding to the knowledge of the local fluid dynamics. Multi-scale computational models include detailed three-dimensional (3D) anatomical information coupled to a lumped parameter network representing the remainder of the circulation. These models output both overall hemodynamic parameters while also enabling to investigate the local fluid dynamics of the aortic arch or the shunt. As an alternative, pure lumped parameter models can also be employed to model Stage 1 palliation, taking advantage of a much lower computational cost, albeit missing the 3D anatomical component. Finally, analytical techniques, such as wave intensity analysis, can be employed to study the Norwood physiology, providing a mechanistic perspective on the ventriculo-arterial coupling for this specific surgical scenario.

Implementing the Sano modification in an experimental model of first-stage palliation of hypoplastic left heart syndrome.

This study describes the implementation of an experimental model of the "Sano" variant (right-ventricle to pulmonary-artery shunt) of the Norwood operation used to treat hypoplastic left-heart syndrome (HLHS). The Sano operation is an alternative to the modified Blalock-Taussig shunt (innominate to pulmonary artery shunt). In the experimental setup, the single ventricle is simulated using a Berlin Heart Excor ventricular assist device and the Sano shunt is constructed by attaching a Tygon tube (6 mm internal diameter, 40 mm long) to the perforated de-airing valve of the Berlin Heart at one end, and to a pulmonary compliance chamber at the other end. The feasibility of the setup was verified by testing two rapid-prototyped patient-specific anatomical models (one without and one with aortic coarctation) under pulsatile flow conditions typical of Norwood patients. Results showed physiological and repeatable pressure and flow signals, as well as physiological values for pulmonary and aortic flow. Shunt flow was regulated by shunt size, and diastolic runoff through the shunt was also observed, both being features of Sano physiology. This system allows for comparing variations of first stage palliation of HLHS in vitro, and it also represents a source of data for validation of computational models.