Objectives: This study sought to quantify average hemodynamic metrics of the Fontan connection as reference for future investigations, compare connection types (intra-atrial vs extracardiac), and identify functional correlates using computational fluid dynamics in a large patient-specific cohort. Fontan hemodynamics, particularly power losses, are hypothesized to vary considerably among patients with a single ventricle and adversely affect systemic hemodynamics and ventricular function if suboptimal. Methods: Fontan connection models were created from cardiac magnetic resonance scans for 100 patients. Phase velocity cardiac magnetic resonance in the aorta, vena cavae, and pulmonary arteries was used to prescribe patient-specific time-averaged flow boundary conditions for computational fluid dynamics with a customized, validated solver. Comparison with 4-dimensional cardiac magnetic resonance velocity data from selected patients was used to provide additional verification of simulations. Indexed Fontan power loss, connection resistance, and hepatic flow distribution were quantified and correlated with systemic patient characteristics. Results: Indexed power loss varied by 2 orders of magnitude, whereas, on average, Fontan resistance was 15% to 20%of published values of pulmonary vascular resistance in single ventricles. A significant inverse relationship was observed between indexed power loss and both systemic venous flow and cardiac index. Comparison by connection type showed no differences between intra-atrial and extracardiac connections. Instead, the least efficient connections revealed adverse consequences from localized Fontan pathway stenosis. Conclusions: Fontan power loss varies from patient to patient, and elevated levels are correlated with lower systemic flow and cardiac index. Fontan connection type does not influence hemodynamic efficiency, but an undersized or stenosed Fontan pathway or pulmonary arteries can be highly dissipative.

Total cavopulmonary connection is the result of a series of palliative surgical repairs performed on patients with single ventricle heart defects. The resulting anatomy has complex and unsteady hemodynamics characterized by flow mixing and flow separation. Although varying degrees of flow pulsatility have been observed in vivo, non-pulsatile (time-averaged) boundary conditions have traditionally been assumed in hemodynamic modeling, and only recently have pulsatile conditions been incorporated without completely characterizing their effect or importance. In this study, 3D numerical simulations with both pulsatile and non-pulsatile boundary conditions were performed for 24 patients with different anatomies and flow boundary conditions from Georgia Tech database. Flow structures, energy dissipation rates and pressure drops were compared under rest and simulated exercise conditions. It was found that flow pulsatility is the primary factor in determining the appropriate choice of boundary conditions, whereas the anatomic configuration and cardiac output had secondary effects. Results show that the hemodynamics can be strongly influenced by the presence of pulsatile flow. However, there was a minimum pulsatility threshold, identified by defining a weighted pulsatility index (wPI), above which the influence was significant. It was shown that when wPI < 30%, the relative error in hemodynamic predictions using time-averaged boundary conditions was less than 10% compared to pulsatile simulations. In addition, when wPI < 50, the relative error was less than 20%. A correlation was introduced to relate wPI to the relative error in predicting the flow metrics with non-pulsatile flow conditions.