by Sebastian Laudenschlager, Samuel Schofield, Nicolas Drysdale, Matthew Stone, Jennifer Romanowicz, Benjamin Frank, Michael DiMaria, Vitaly O. Kheyfets, Mehdi Hedjazi-Moghari
Children with single ventricle heart disease typically require a series of three operations, (1) Norwood, (2) Glenn, and (3) Fontan, which ultimately results in complete separation of the pulmonary and systemic circuits to improve pulmonary/systemic circulation. In the last stage, the Fontan operation, the inferior vena cava (IVC) is connected to the pulmonary arteries (PAs), allowing the remainder of deoxygenated blood to passively flow to the pulmonary circuit. It is hypothesized that optimizing the Fontan anatomy would lead to decreased power loss and more balanced hepatic flow distribution. One approach to optimizing the geometry is to create a patient-specific digital twin to simulate various configurations of the Fontan conduit, which requires a computational model of the proximal PA anatomy and resistance, as well as the distal Pulmonary Vascular Resistance (PVR), at the Glenn stage. To that end, an optimization pipeline was developed using 3D computational fluid dynamics (CFD) and 0D lumped parameter (LP) simulations to iteratively refine the PVR of each lung by minimizing the simulated flow and pressure error relative to patients’ cardiac magnetic resonance (CMR) and catheterization (CATH) data. While the PVR can also be estimated directly by computing the ratio of pressure gradients and flow from CATH and CMR data, the computational approach can separately identify the different components of PVR along the Glenn pathway, allowing for a more detailed depiction of the Glenn vasculature. Results indicate good correlation between the optimized PVR of the CFD and LP models (n = 16), with an intraclass correlation coefficient (ICC) of 0.998 (p = 0.976) and 0.991 (p = 0.943) for the left and right lung, respectively. Furthermore, compared to CMR flow and CATH pressure data, the optimized PVR estimates result in mean outlet flow and pressure errors of less than 5%. The optimized PVR estimates also agree well with the computed PVR estimates from CATH pressure and CMR flow for both lungs, yielding a mean difference of less than 4%.