Fontan patients require a balanced hepatic blood flow distribution (HFD) to prevent pulmonary arteriovenous malformations. Currently, HFD is quantified by tracking Fontan conduit flow, assuming hepatic venous (HV) flow to be uniformly distributed within the Fontan conduit. However, this assumption may be unvalid leading to inaccuracies in HFD quantification with potential clinical impact. The aim of this study was to (i) assess the mixing of HV flow and inferior vena caval (IVC) flow within the Fontan conduit and (ii) quantify HFD by directly tracking HV flow and quantitatively comparing results with the conventional approach. Patient-specific, time-resolved computational fluid dynamic models of 15 total cavopulmonary connections were generated, including the HV and subhepatic IVC. Mixing of HV and IVC flow, on a scale between 0 (no mixing) and 1 (perfect mixing), was assessed at the caudal and cranial Fontan conduit. HFD was quantified by tracking particles from the caudal (HFDcaudal conduit) and cranial (HFDcranial conduit) conduit and from the hepatic veins (HFDHV). HV flow was non-uniformly distributed at both the caudal (mean mixing 0.66 ± 0.13) and cranial (mean 0.79 ± 0.11) level within the Fontan conduit. On a cohort level, differences in HFD between methods were significant but small; HFDHV (51.0 ± 20.6%) versus HFDcaudal conduit (48.2 ± 21.9%, p = 0.033) or HFDcranial conduit (48.0 ± 21.9%, p = 0.044). However, individual absolute differences of 8.2-14.9% in HFD were observed in 4/15 patients. HV flow is non-uniformly distributed within the Fontan conduit. Substantial individual inaccuracies in HFD quantification were observed in a subset of patients with potential clinical impact.
Keywords: Fontan; computational fluid dynamics; flow; hepatic flow distribution; hepatic venous; mixing.