Objective: Our objective was to analyze 3-dimensional (3D) blood flow patterns within the total cavopulmonary connection (TCPC) using in vivo phase contrast magnetic resonance imaging (PC MRI). Methods: Sixteen single-ventricle patients were prospectively recruited at 2 leading pediatric institutions for PC MRI evaluation of their Fontan pathway. Patients were divided into 2 groups. Group 1 comprised 8 patients with an extracardiac (EC) TCPC, and group 2 comprised 8 patients with a lateral tunnel (LT) TCPC. A coronal stack of 5 to 10 contiguous PC MRI slices with 3D velocity encoding (5-9 ms resolution) was acquired and a volumetric flow field was reconstructed. Results: Analysis revealed large vortices in LT TCPCs and helical flow structures in EC TCPCs. On average, there was no difference between LT and EC TCPCs in the proportion of inferior vena cava flow going to the left pulmonary artery (43% ± 7% vs 46% ± 5%; P = .34). However, for EC TCPCs, the presence of a caval offset was a primary determinant of inferior vena caval flow distribution to the pulmonary arteries with a significant bias to the offset side. Conclusions: 3D flow structures within LT and EC TCPCs were reconstructed and analyzed for the first time using PC MRI. TCPC flow patterns were shown to be different, not only on the basis of LT or EC considerations, but with significant influence from the superior vena cava connection as well. This work adds to the ongoing body of research demonstrating the impact of TCPC geometry on the overall hemodynamic profile.

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.