Non-Newtonian Effects on Patient-Specific Modeling of Fontan Hemodynamics

Advancement in computational simulation and validation of congenital heart disease: a review

The improvement in congenital heart disease (CHD) treatment and management has increased the life expectancy in infants. However, the long-term efficacy is difficult to assess and thus, computational modelling has been applied for evaluating this. Here, we provide an overview of the applications of computational modelling in CHD based on three categories; CHD involving large blood vessels only, heart chambers only, and CHD that occurs at multiple heart structures. We highlight the advancement of computational simulation of CHD that uses multiscale and multiphysics modelling to ensure a complete representation of the heart and circulation. We provide a brief future direction of computational modelling of CHD such as to include growth and remodelling, detailed conduction system, and occurrence of myocardial infarction. We also proposed validation technique using advanced three-dimensional (3D) printing and particle image velocimetry (PIV) technologies to improve the model accuracy.

Impact of Age-Related Change in Caval Flow Ratio on Hepatic Flow Distribution in the Fontan Circulation

BACKGROUND

The Fontan operation is a palliative technique for patients born with single ventricle heart disease. The superior vena cava (SVC), inferior vena cava (IVC), and hepatic veins are connected to the pulmonary arteries in a total cavopulmonary connection by an extracardiac conduit or a lateral tunnel connection. A balanced hepatic flow distribution (HFD) to both lungs is essential to prevent pulmonary arteriovenous malformations and cyanosis. HFD is highly dependent on the local hemodynamics. The effect of age-related changes in caval inflows on HFD was evaluated using cardiac magnetic resonance data and patient-specific computational fluid dynamics modeling.

METHODS

SVC and IVC flow from 414 patients with Fontan were collected to establish a relationship between SVC:IVC flow ratio and age. Computational fluid dynamics modeling was performed in 60 (30 extracardiac and 30 lateral tunnel) patient models to quantify the HFD that corresponded to patient ages of 3, 8, and 15 years, respectively.

RESULTS

SVC:IVC flow ratio inverted at ≈8 years of age, indicating a clear shift to lower body flow predominance. Our data showed that variation of HFD in response to age-related changes in caval inflows (SVC:IVC, 2, 1, and 0.5 corresponded to ages, 3, 8, and 15+, respectively) was not significant for extracardiac but statistically significant for lateral tunnel cohorts. For all 3 caval inflow ratios, a positive correlation existed between the IVC flow distribution to both the lungs and the HFD. However, as the SVC:IVC ratio changed from 2 to 0.5 (age, 3-15+) years, the correlation's strength decreased from 0.87 to 0.64, due to potential flow perturbation as IVC flow momentum increased.

CONCLUSIONS

Our analysis provided quantitative insights into the impact of the changing caval inflows on Fontan's long-term HFD, highlighting the importance of SVC:IVC variations over time on Fontan's long-term hemodynamics. These findings broaden our understanding of Fontan hemodynamics and patient outcomes.

Decreased erythrocyte aggregation in Glenn and Fontan: univentricular circulation as a rheologic disease model

BACKGROUND

In the Fontan palliation for single ventricle heart disease (SVHD), pulmonary blood flow is non-pulsatile/passive, low velocity, and low shear, making viscous power loss a critical determinant of cardiac output. The rheologic properties of blood in SVHD patients are essential for understanding and modulating their limited cardiac output and they have not been systematically studied. We hypothesize that viscosity is decreased in single ventricle circulation.

METHODS

We evaluated whole blood viscosity, red blood cell (RBC) aggregation, and RBC deformability to evaluate changes in healthy children and SVHD patients. We altered suspending media to understand cellular and plasma differences contributing to rheologic differences.

RESULTS

Whole blood viscosity was similar between SVHD and healthy at their native hematocrits, while viscosity was lower at equivalent hematocrits for SVHD patients. RBC deformability is increased, and RBC aggregation is decreased in SVHD patients. Suspending SVHD RBCs in healthy plasma resulted in increased RBC aggregation and suspending healthy RBCs in SVHD plasma resulted in lower RBC aggregation.

CONCLUSIONS

Hematocrit corrected blood viscosity is lower in SVHD vs. healthy due to decreased RBC aggregation and higher RBC deformability, a viscous adaptation of blood in patients whose cardiac output is dependent on minimizing viscous power loss.

IMPACT

Patients with single ventricle circulation have decreased red blood cell aggregation and increased red blood cell deformability, both of which result in a decrease in blood viscosity across a large shear rate range. Since the unique Fontan circulation has very low-shear and low velocity flow in the pulmonary arteries, blood viscosity plays an increased role in vascular resistance, therefore this work is the first to describe a novel mechanism to target pulmonary vascular resistance as a modifiable risk factor. This is a novel, modifiable risk factor in this patient population.

In vivo hepatic flow distribution by computational fluid dynamics can predict pulmonary flow distribution in patients with Fontan circulation

In Fontan patients, a lung deprived of hepatic blood may develop pulmonary arterio-venous malformations (PAVMs) resulting in shunting, reduced pulmonary vascular resistance (PVR) and decreased oxygenation. To provide guidance for corrective invasive interventions, we aimed to non-invasively determine how the hepatic to pulmonary blood flow balance correlates with pulmonary flow, PVR, and with oxygen saturation. Magnetic resonance imaging (MRI) data from eighteen Fontan patients (eight females, age 3-14 years) was used to construct patient-specific computational fluid dynamics (CFD) models to calculate the hepatic to pulmonary blood flow. This was correlated with pulmonary vein flow, simulated PVR and oxygen saturation. Clinical applicability of the findings was demonstrated with an interventional patient case. The hepatic to pulmonary blood flow balance correlated with right/left pulmonary vein flow (R2 = 0.50), left/right simulated PVR (R2 = 0.47), and oxygen saturation at rest (R2 = 0.56). In the interventional patient, CFD predictions agreed with post-interventional MRI measurements and with regressions in the cohort. The balance of hepatic blood to the lungs has a continuous effect on PVR and oxygen saturation, even without PAVM diagnosis. MRI combined with CFD may help in planning of surgical and interventional designs affecting the hepatic to pulmonary blood flow balance in Fontan patients.

Impact of Age-related change in Caval Flow Ratio on Hepatic Flow Distribution in Fontan

Background

The Fontan operation is a palliative technique for patients born with single ventricle heart disease. The superior vena cava (SVC), inferior vena cava (IVC), and hepatic veins are connected to the pulmonary arteries in a total cavopulmonary connection by an extracardiac (EC) conduit or a lateral tunnel (LT) connection. A balanced hepatic flow distribution (HFD) to both lungs is essential to prevent pulmonary arteriovenous malformations and cyanosis. HFD is highly dependent on the local hemodynamics.

Objective

The effect of age-related changes in caval inflows on HFD was evaluated using cardiac MRI (CMR) data and patient-specific computational fluid dynamics (CFD) modeling.

Methods

SVC and IVC flow from 414 Fontan patients were collected to establish a relationship between SVC:IVC flow ratio and age. CFD modeling was performed in 60 (30 EC and 30 LT) patient models to quantify the HFD that corresponded to patient ages of 3, 8, and 15 years, respectively.

Results

SVC:IVC flow ratio inverted at ∼8 years of age, indicating a clear shift to lower body flow predominance. Our data showed that variation of HFD in response to age-related changes in caval inflows (SVC:IVC = 2,1, and 0.5 corresponded to ages 3, 8, and 15+ respectively) was not significant for EC but statistically significant for LT cohorts. For all three caval inflow ratios, a positive correlation existed between the IVC flow distribution to both the lungs and the HFD. However, as the SVC:IVC ratio changed from 2→0.5 (age 3→15+), the correlation's strength decreased from 0.87→0.64, due to potential flow perturbation as IVC flow momentum increased.

Conclusion

Our analysis provided quantitative insights into the impact of the changing caval inflows on Fontan's long-term HFD, highlighting the importance of including SVC:IVC variations over time to understand Fontan's long-term hemodynamics. These findings broaden our understanding of Fontan hemodynamics and patient outcomes.

Clinical Perspective

With improvement in standard of care and management of single ventricle patients with Fontan physiology, the population of adults with Fontan circulation is increasing. Consequently, there is a clinical need to comprehend the impact of patient growth on Fontan hemodynamics. Using CMR data, we were able to quantify the relationship between changing caval inflows and somatic growth. We then used patient-specific computational flow modeling to quantify how this relationship affected the distribution of long-term hepatic flow in extracardiac and lateral tunnel Fontan types. Our findings demonstrated the significance of including SVC:IVC changes over time in CFD modeling to learn more about the long-term hemodynamics of Fontan. Fontan surgical approaches are increasingly planned and optimized using computational flow modeling. For a patient undergoing a Fontan procedure, the workflow presented in this study that takes into account the variations in Caval inflows over time can aid in predicting the long-term hemodynamics in a planned Fontan pathway.

Numerical simulation of the cavopulmonary connection flow with conduit stenoses of varying configurations

The circulation in the total cavopulmonary connection (TCPC) is a low-energy system which operation and efficiency are subjected to multiple factors. Some retrospective studies report that the abnormal narrowing of vessels in the system, i.e. stenosis, is one of the most dangerous geometric factors which can result in heart failure. In the present study, the effect of varying extracardiac conduit (ECC) stenosis on the hemodynamics in a surrogate TCPC model is investigated using high-fidelity numerical simulations. The efficiency of the surrogate TCPC model was quantified according to the power loss, relative perfusion in lungs and the percentage of conduit surface area with abnormally low and high wall shear stress for venous flow. Additionally, the impact of respiration and asymmetry in the stenosis geometry to the system was examined. The results show that the flow in the TCPC model exhibits pronounced unsteadiness even under the steady initial boundary conditions, while the uneven pulmonary flow distribution and the presence of the ECC stenosis amplify the chaotic nature of the flow. Energy efficiency of the system is shown to strongly correlate with amount of vortical structures in the model and their range of scales. Finally, the study demonstrates that the presence of respiration in the model adds to perturbations in the flow which causes increase in the power loss. Results obtained in the study provide valuable insights on how the ECC stenosis effect the flow in the surrogate TCPC model under different flow conditions.

Thermal and Postural Effects on Fluid Mixing and Irrigation Patterns for Intraventricular Hemorrhage Treatment

Intraventricular hemorrhage is characterized by blood leaking into the cerebral ventricles and mixing with cerebrospinal fluid. A standard treatment method involves inserting a passive drainage catheter, known as an external ventricular drain (EVD), into the ventricle. EVDs have common adverse complications, including the occlusion of the catheter, that may lead to permanent neural damage or even mortality. In order to prevent such complications, a novel dual-lumen catheter (IRRAflow®) utilizing an active fluid exchange mechanism has been recently developed. However, the fluid dynamics of the exchange system have not been investigated. In this study, convective flow in a three-dimensional cerebral lateral ventricle with an inserted catheter is evaluated using an in-house lattice-Boltzmann-based fluid-solid interaction solver. Different treatment conditions are simulated, including injection temperature and patient position. Thermal and gravitational effects on medication distribution are studied using a dye simulator based on a recently-introduced (pseudo)spectral convection-diffusion equation solver. The effects of injection temperature and patient position on catheter performance are presented and discussed in terms of hematoma irrigation, vortical structures, mixing, and medication volume distribution. Results suggest that cold-temperature injections can increase catheter efficacy in terms of dye distribution and irrigation potential, both of which can be further guided by patient positioning.

Haemodynamic performance of 16-20-mm extracardiac Goretex conduits in adolescent Fontan patients at rest and during simulated exercise

OBJECTIVES

To date, it is not known if 16-20-mm extracardiac conduits are outgrown during somatic growth from childhood to adolescence. This study aims to determine total cavopulmonary connection (TCPC) haemodynamics in adolescent Fontan patients at rest and during simulated exercise and to assess the relationship between conduit size and haemodynamics.

METHODS

Patient-specific, magnetic resonance imaging-based computational fluid dynamic models of the TCPC were performed in 51 extracardiac Fontan patients with 16-20-mm conduits. Power loss, pressure gradient and normalized resistance were quantified in rest and during simulated exercise. The cross-sectional area (CSA) (mean and minimum) of the vessels of the TCPC was determined and normalized for flow rate (mm2/l/min). Peak (predicted) oxygen uptake was assessed.

RESULTS

The median age was 16.2 years (Q1-Q3 14.0-18.2). The normalized mean conduit CSA was 35-73% smaller compared to the inferior and superior vena cava, hepatic veins and left/right pulmonary artery (all P < 0.001). The median TCPC pressure gradient was 0.7 mmHg (Q1-Q3 0.5-0.8) and 2.0 (Q1-Q3 1.4-2.6) during rest and simulated exercise, respectively. A moderate-strong inverse non-linear relationship was present between normalized mean conduit CSA and TCPC haemodynamics in rest and exercise. TCPC pressure gradients of ≥1.0 at rest and ≥3.0 mmHg during simulated exercise were observed in patients with a conduit CSA ≤ 45 mm2/l/min and favourable haemodynamics (<1 mmHg during both rest and exercise) in conduits ≥125 mm2/l/min. Normalized TCPC resistance correlated with (predicted) peak oxygen uptake.

CONCLUSIONS

Extracardiac conduits of 16-20 mm have become relatively undersized in most adolescent Fontan patients leading to suboptimal haemodynamics.

Hemodynamic characteristics in a cerebral aneurysm model using non-Newtonian blood analogues

This study aims to develop an experimentally validated computational fluid dynamics (CFD) model to estimate hemodynamic characteristics in cerebral aneurysms (CAs) using non-Newtonian blood analogues. Blood viscosities varying with shear rates were measured under four temperatures first, which serves as the reference for the generation of blood analogues. Using the blood analogue, particle image velocimetry (PIV) measurements were conducted to quantify flow characteristics in a CA model. Then, using the identical blood properties in the experiment, CFD simulations were executed to quantify the flow patterns, which were used to compare with the PIV counterpart. Additionally, hemodynamic characteristics in the simplified Newtonian and non-Newtonian models were quantified and compared using the experimentally validated CFD model. Results showed the proposed non-Newtonian viscosity model can predict blood shear-thinning properties accurately under varying temperatures and shear rates. Another developed viscosity model based on the blood analogue can well represent blood rheological properties. The comparisons in flow characteristics show good agreements between PIV and CFD, demonstrating the developed CFD model is qualified to investigate hemodynamic factors within CAs. Furthermore, results show the differences of absolute values were insignificant between Newtonian and non-Newtonian fluids in the distributions of wall shear stress (WSS) and oscillatory shear index (OSI) on arterial walls. However, not only does the simplified Newtonian model underestimate WSS and OSI in most regions of the aneurysmal sac, but it also makes mistakes in identifying the high OSI regions on the sac surface, which may mislead the hemodynamic assessment on the pathophysiology of CAs.

Computational Fluid Dynamics Support for Fontan Planning in Minutes, Not Hours: The Next Step in Clinical Pre-Interventional Simulations

Computational fluid dynamics (CFD) modeling may aid in planning of invasive interventions in Fontan patients. Clinical application of current CFD techniques is however limited by complexity and long computation times. Therefore, we validated a “lean” CFD method to magnetic resonance imaging (MRI) and an “established” CFD method, ultimately aiming to reduce complexity to enable predictive CFD during ongoing interventions. Fifteen Fontan patients underwent MRI for CFD modeling. The differences between lean and established approach, in hepatic and total flow percentage to the left pulmonary artery (%LPA), power loss and relative wall shear stress area were 1.5 ± 4.0%, -0.17 ± 1.1%, -0.055 ± 0.092 mW and 1.1 ± 1.4%. Compared with MRI, the lean and established method showed a bias in %LPA of -1.9 ± 3.4% and -1.8 ± 3.1%. Computation time was for the lean and established approach 3.0 ± 2.0 min and 7.0 ± 3.4 h, respectively. We conclude that the proposed lean method provides fast and reliable results for future CFD support during interventions.

Engineering Perspective on Cardiovascular Simulations of Fontan Hemodynamics: Where Do We Stand with a Look Towards Clinical Application

BACKGROUND

Cardiovascular simulations for patients with single ventricles undergoing the Fontan procedure can assess patient-specific hemodynamics, explore surgical advances, and develop personalized strategies for surgery and patient care. These simulations have not yet been broadly accepted as a routine clinical tool owing to a number of limitations. Numerous approaches have been explored to seek innovative solutions for improving methodologies and eliminating these limitations.

PURPOSE

This article first reviews the current state of cardiovascular simulations of Fontan hemodynamics. Then, it will discuss the technical progress of Fontan simulations with the emphasis of its clinical impact, noting that substantial improvements have been made in the considerations of patient-specific anatomy, flow, and blood rheology. The article concludes with insights into potential future directions involving clinical validation, uncertainty quantification, and computational efficiency. The advancements in these aspects could promote the clinical usage of Fontan simulations, facilitating its integration into routine clinical practice.

The Rheology of the Carotid Sinus: A Path Toward Bioinspired Intervention

The association between blood viscosity and pathological conditions involving a number of organ systems is well known. However, how the body measures and maintains appropriate blood viscosity is not well-described. The literature endorsing the function of the carotid sinus as a site of baroreception can be traced back to some of the earliest descriptions of digital pressure on the neck producing a drop in blood delivery to the brain. For the last 30 years, improved computational fluid dynamic (CFD) simulations of blood flow within the carotid sinus have demonstrated a more nuanced understanding of the changes in the region as it relates to changes in conventional metrics of cardiovascular function, including blood pressure. We suggest that the unique flow patterns within the carotid sinus may make it an ideal site to transduce flow data that can, in turn, enable real-time measurement of blood viscosity. The recent characterization of the PIEZO receptor family in the sinus vessel wall may provide a biological basis for this characterization. When coupled with other biomarkers of cardiovascular performance and descriptions of the blood rheology unique to the sinus region, this represents a novel venue for bioinspired design that may enable end-users to manipulate and optimize blood flow.

A Tribute to Ajit Yoganathan's Cardiovascular Fluid Mechanics Lab: A Survey of Its Contributions to Our Understanding of the Physiology and Management of Single-Ventricle Patients

INTRODUCTION

Among patients with congenital heart disease, those born with only a single working ventricle represent a particularly complex sub-population, typically requiring multiple surgeries and suffering from high levels of mortality and morbidity. Their cardiac care is complex and has evolved considerably since surgical palliation was first introduced more than 50 years ago. Improvements in treatment have been driven both by growing clinical experience and by knowledge gained through experimental and computational studies of blood flow in these patients. The Cardiovascular Fluid Mechanics Lab at the Georgia Institute of Technology, founded 30 years ago by Dr. Ajit Yoganathan, has pioneered work in this field.

METHODS

In this review, key contributions of Dr. Yoganathan's Cardiovascular Fluid Dynamics Lab are surveyed, including experimental flow loop studies as well as computational fluid dynamics analyses that address many of the critical challenges that cardiologists and surgeons face in treating these patients, including how to reconstruct cardiovascular anatomy to minimize power loss, balance blood flow distribution at key bifurcation points, and avoid other unfavorable hemodynamic conditions.

CONCLUSIONS

Among many contributions in this field, work from the Cardiovascular Fluid Mechanics Lab has led to novel medical devices and patient-specific computational modeling workflows and software tools. These key contributions from this group have enhanced our understanding of the physiology and management of single-ventricle patients.

Is Doppler Echocardiography Adequate for Surgical Planning of Single Ventricle Patients?

PURPOSE

Surgical planning has shown great potential for optimizing outcomes for patients affected by single ventricle (SV) malformations. Phase-contrast magnetic resonance imaging (PC-MRI) is the routine technique used for flow acquisition in the surgical planning paradigm. However, PC-MRI may suffer from possible artifacts in certain cases; furthermore, this technology may not be readily available for patients in low and lower-middle-income countries. Therefore, this study aims to investigate the effectiveness of using Doppler echocardiography (echo-Doppler) for flow acquisitions of SV surgical planning.

METHODS

This study included eight patients whose blood flow data was acquired by both PC-MRI and echo-Doppler. A virtual surgery platform was used to generate two surgical options for each patient: (1) a traditional Fontan conduit and (2) a Y-graft. Computational fluid dynamics (CFD) simulations were conducted using the two flow acquisitions to assess clinically relevant hemodynamic metrics: indexed power loss (iPL) and hepatic flow distribution (HFD).

RESULTS

Differences exist in flow data acquired by PC-MRI and echo-Doppler, but no statistical significance was obtained. Flow fields, therefore, exhibit discrepancies between simulations using flow acquisitions by PC-MRI and echo-Doppler. In virtual surgery, the two surgical options were ranked based on these metrics. No difference was observed in the ranking of surgical options between using different flow acquisitions.

CONCLUSION

Doppler echocardiography is an adequate alternative approach to acquire flow data for SV surgical planning. This finding encourages broader usage of SV surgical planning with echo-Doppler when MRI may present artifacts or is not available, especially in low and lower-middle-income countries.

Non-uniform mixing of hepatic venous flow and inferior vena cava flow in the Fontan conduit

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 (HFD_{caudal conduit}) and cranial (HFD_{cranial conduit}) conduit and from the hepatic veins (HFD_HV). 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; HFD_HV (51.0 ± 20.6%) versus HFD_{caudal conduit} (48.2 ± 21.9%, p = 0.033) or HFD_{cranial 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.

A Simplified In Silico Model of Left Ventricular Outflow in Patients After Transcatheter Mitral Valve Replacement with Anterior Leaflet Laceration

In silico modeling has been proposed as a tool to simulate left ventricular (LV) outflow tract (LVOT) obstruction in patients undergoing transcatheter mitral valve replacement (TMVR). This study validated a simplified approach to simulate LV outflow hemodynamics in the setting of TMVR with anterior leaflet laceration, a clinical technique used to mitigate the risk of LVOT obstruction. Personalized, 3-dimensional computational fluid dynamics models were developed from computed tomography images of six patients who underwent TMVR with anterior leaflet laceration. LV outflow hemodynamics were simulated using the patient-specific anatomy and the peak systolic flow rate as boundary conditions. The peak outflow velocity, a clinically relevant hemodynamic metric, was extracted from each simulation (v_sim-peak) and compared with the clinical measurement from Doppler echocardiography (v_clin-peak) for validation. In silico models were successfully developed and implemented for all patients. The pre-processing time was 2 h per model and the simulation could be completed within 3 h. In three patients, the lacerated anterior leaflet exposed open cells of the transcatheter valve to flow. Good agreement was obtained between v_sim-peak and v_clin-peak (r = 0.97, p < 0.01) with average discrepancies of 5 ± 2% and 14 ± 1% for patients with exposed and unexposed cells of the transcatheter valve, respectively. The proposed in silico modeling paradigm therefore simulated LV outflow hemodynamics in a time-efficient manner and demonstrated good agreement with clinical measurements. Future studies should investigate the ability of this paradigm to support clinical applications.

Fluid-Structure Interaction Simulation of an Intra-Atrial Fontan Connection

Simple Summary

A fluid-structure interaction (FSI) simulation of an intra-atrial Fontan connection was performed. Power loss and pressure drop results fluctuated less during the FSI simulation than during the simulation run with rigid walls, but there were no observable differences in time-averaged pressure drop, connection power loss or hepatic flow distribution. These results suggested that employing a rigid wall is a reasonable assumption when evaluating time-averaged hemodynamic quantities of the Fontan connection under resting breath-held flow conditions.

Abstract

Total cavopulmonary connection (TCPC) hemodynamics has been hypothesized to be associated with long-term complications in single ventricle heart defect patients. Rigid wall assumption has been commonly used when evaluating TCPC hemodynamics using computational fluid dynamics (CFD) simulation. Previous study has evaluated impact of wall compliance on extra-cardiac TCPC hemodynamics using fluid-structure interaction (FSI) simulation. However, the impact of ignoring wall compliance on the presumably more compliant intra-atrial TCPC hemodynamics is not fully understood. To narrow this knowledge gap, this study aims to investigate impact of wall compliance on an intra-atrial TCPC hemodynamics. A patient-specific model of an intra-atrial TCPC is simulated with an FSI model. Patient-specific 3D TCPC anatomies were reconstructed from transverse cardiovascular magnetic resonance images. Patient-specific vessel flow rate from phase-contrast magnetic resonance imaging (MRI) at the Fontan pathway and the superior vena cava under resting condition were prescribed at the inlets. From the FSI simulation, the degree of wall deformation was compared with in vivo wall deformation from phase-contrast MRI data as validation of the FSI model. Then, TCPC flow structure, power loss and hepatic flow distribution (HFD) were compared between rigid wall and FSI simulation. There were differences in instantaneous pressure drop, power loss and HFD between rigid wall and FSI simulations, but no difference in the time-averaged quantities. The findings of this study support the use of a rigid wall assumption on evaluation of time-averaged intra-atrial TCPC hemodynamic metric under resting breath-held condition.