Background: The development of clinically applicable fluid-structure interaction (FSI) models of the left heart is inherently challenging when using in vivo cardiovascular magnetic resonance (CMR) data for validation, due to the lack of a well-controlled system where detailed measurements of the ventricular wall motion and flow field are available a priori. The purpose of this study was to (a) develop a clinically relevant, CMR-compatible left heart physical model; and (b) compare the left ventricular (LV) volume reconstructions and hemodynamic data obtained using CMR to laboratory-based experimental modalities. Methods: The LV was constructed from optically clear flexible silicone rubber. The geometry was based off a healthy patient's LV geometry during peak systole. The LV phantom was attached to a left heart simulator consisting of an aorta, atrium, and systemic resistance and compliance elements. Experiments were conducted for heart rate of 70 bpm. Wall motion measurements were obtained using high speed stereo-photogrammetry (SP) and cine-CMR, while flow field measurements were obtained using digital particle image velocimetry (DPIV) and phase-contrast magnetic resonance (PC-CMR). Results: The model reproduced physiologically accurate hemodynamics (aortic pressure = 120/80 mmHg; cardiac output = 3.5 L/min). DPIV and PC-CMR results of the center plane flow within the ventricle matched, both qualitatively and quantitatively, with flow from the atrium into the LV having a velocity of about 1.15 m/s for both modalities. The normalized LV volume through the cardiac cycle computed from CMR data matched closely to that from SP. The mean difference between CMR and SP was 5.5 ± 3.7 %. Conclusions: The model presented here can thus be used for the purposes of: (a) acquiring CMR data for validation of FSI simulations, (b) determining accuracy of cine-CMR reconstruction methods, and (c) conducting investigations of the effects of altering anatomical variables on LV function under normal and disease conditions.
Computational fluid dynamics (CFD) tools have been extensively applied to study the hemodynamics in the total cavopulmonary connection (TCPC) in patients with only a single functioning ventricle. Without the contraction of a sub-pulmonary ventricle, pulsatility of flow through this connection is low and variable across patients, which is usually neglected in most numerical modeling studies. Recent studies suggest that such pulsatility can be non-negligible and can be important in hemodynamic predictions. The goal of this work is to compare the results of an in-house numerical methodology for simulating pulsatile TCPC flow with experimental results. Digital particle image velocimetry (DPIV) was acquired on TCPC in vitro models to evaluate the capability of the CFD tool in predicting pulsatile TCPC flow fields. In vitro hemodynamic measurements were used to compare the numerical prediction of power loss across the connection. The results demonstrated the complexity of the pulsatile TCPC flow fields and the validity of the numerical approach in simulating pulsatile TCPC flow dynamics in both idealized and complex patient specific models.
The ultimate goal of Fontan surgical planning is to provide additional insights into the clinical decision-making process. In its current state, surgical planning offers an accurate hemodynamic assessment of the pre-operative condition, provides anatomical constraints for potential surgical options, and produces decent post-operative predictions if boundary conditions are similar enough between the pre-operative and post-operative states. Moving forward, validation with post-operative data is a necessary step in order to assess the accuracy of surgical planning and determine which methodological improvements are needed. Future efforts to automate the surgical planning process will reduce the individual expertise needed and encourage use in the clinic by clinicians. As post-operative physiologic predictions improve, Fontan surgical planning will become an more effective tool to accurately model patient-specific hemodynamics.
Computational fluid dynamic (CFD) simulations are widely utilized to assess Fontan hemodynamics that are related to long-term complications. No previous studies have systemically investigated the effects of using different inlet velocity profiles in Fontan simulations. This study implements real, patient-specific velocity profiles for numerical assessment of Fontan hemodynamics using CFD simulations. Four additional, artificial velocity profiles were used for comparison: (1) flat, (2) parabolic, (3) Womersley, and (4) parabolic with inlet extensions [to develop flow before entering the total cavopulmonary connection (TCPC)]. The differences arising from the five velocity profiles, as well as discrepancies between the real and each of the artificial velocity profiles, were quantified by examining clinically important metrics in TCPC hemodynamics: power loss (PL), viscous dissipation rate (VDR), hepatic flow distribution, and regions of low wall shear stress. Statistically significant differences were observed in PL and VDR between simulations using real and flat velocity profiles, but differences between those using real velocity profiles and the other three artificial profiles did not reach statistical significance. These conclusions suggest that the artificial velocity profiles (2)–(4) are acceptable surrogates for real velocity profiles in Fontan simulations, but parabolic profiles are recommended because of their low computational demands and prevalent applicability.
Cardiovascular simulations have great potential as a clinical tool for planning and evaluating patient-specific treatment strategies for those suffering from congenital heart diseases, specifically Fontan patients. However, several bottlenecks have delayed wider deployment of the simulations for clinical use; the main obstacle is simulation cost. Currently, time-averaged clinical flow measurements are utilized as numerical boundary conditions (BCs) in order to reduce the computational power and time needed to offer surgical planning within a clinical time frame. Nevertheless, pulsatile blood flow is observed in vivo, and its significant impact on numerical simulations has been demonstrated. Therefore, it is imperative to carry out a comprehensive study analyzing the sensitivity of using time-averaged BCs. In this study, sensitivity is evaluated based on the discrepancies between hemodynamic metrics calculated using time-averaged and pulsatile BCs; smaller discrepancies indicate less sensitivity. The current study incorporates a comparison between 3D patient-specific CFD simulations using both the time-averaged and pulsatile BCs for 101 Fontan patients. The sensitivity analysis involves two clinically important hemodynamic metrics: hepatic flow distribution (HFD) and indexed power loss (iPL). Paired demographic group comparisons revealed that HFD sensitivity is significantly different between single and bilateral superior vena cava cohorts but no other demographic discrepancies were observed for HFD or iPL. Multivariate regression analyses show that the best predictors for sensitivity involve flow pulsatilities, time-averaged flow rates, and geometric characteristics of the Fontan connection. These predictors provide patient-specific guidelines to determine the effectiveness of analyzing patient-specific surgical options with time-averaged BCs within a clinical time frame.
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.