Influence of bypass angles on extracardiac Fontan connections: a numerical study

Patient-specific closed-loop model of the fontan circulation: Calibration and validation

The Fontan circulation, designed for managing patients with a single functional ventricle, presents challenges in long-term outcomes. Computational methods offer potential solutions, yet their application in cardiology practice remains largely unexplored. Our aim was to assess the ability of a patient-specific, closed-loop, reduced-order blood flow model to simulate pulsatile blood flow in the Fontan circulation. Using one-dimensional models, we simulated the aorta, superior and inferior venae cavae, and right and left pulmonary arteries, while lumping heart chambers and remaining vessels into zero-dimensional models. The model was calibrated with patient-specific haemodynamic data from combined cardiac catheterisation and magnetic resonance exams, using a novel physics-based stepwise methodology involving simpler open-loop models. Testing on a 10-year-old, anesthetised patient, demonstrated the model's capability to replicate pulsatile pressure and flow in the larger vessels and ventricular pressure. Average relative errors in mean pressure and flow were 2.9 % and 3.6 %, with average relative point-to-point errors (RPPE) in pressure and flow at 5.2 % and 16.0 %. Comparing simulation results to measurements, mean aortic pressure and flow values were 50.7 vs. 50.4 mmHg and 41.6 vs. 41.9 ml/s, respectively, while ventricular pressure values were 28.7 vs. 27.4 mmHg. The model accurately described time-varying ventricular volume with a RPPE of 2.9 %, with mean, minimum, and maximum ventricular volume values for simulation results vs. measurements at 59.2 vs. 58.2 ml, 38.0 vs. 37.6 ml, and 76.0 vs. 74.4 ml, respectively. It provided physiologically realistic predictions of haemodynamic changes from pulmonary vasodilation and atrial fenestration opening. The new model and calibration methodology are freely available, offering a platform to virtually investigate the Fontan circulation's response to clinical interventions and explore potential mechanisms of Fontan failure. Future efforts will concentrate on broadening the model's applicability to a wider range of patient populations and clinical scenarios, as well as testing its effectiveness.

Numerical Simulation of Hemodynamics in Two Models for Total Anomalous Pulmonary Venous Connection Surgery

OBJECTIVE

To numerically compare the prospective hemodynamic outcomes between a new window surgery and a traditional surgery in the treatment of supracardiac total anomalous pulmonary venous connection (S-TAPVC).

METHODS

A 3D geometry model, composed of pulmonary vein (PV) and left atrium (LA), was reconstructed based on summarized data with S-TAPVC. Two surgery models were established based on this model. One is the traditional surgery model, where an elliptical anastomosis was created by incising and stitching the LA and the common vein (CV) along the axis of the CV. The other is the new window surgery model, where the CV was incised with an H-shaped orifice, and LA was incised with a transposed H-shaped orifice, and then the orifice edges were stitched like a window. Two models with a relative cross sectional area (RCSA) of 300 mm2/m2 and 500 mm2/m2 were established, which correspond to traditional surgery and window surgery. Numerical simulation of hemodynamics was carried out. The velocity, left atrium and pulmonary vein pressure, the pressure difference of anastomosis and the energy conversion efficiency were analyzed to evaluate the prospective hemodynamic outcomes of these two operations.

RESULTS

Window surgery presented a lower blood flow velocity, pressure difference, and the WSS at the anastomosis, compared to traditional surgery. In terms of energy loss, the power conversion efficiency of window surgery was significantly higher than that of traditional surgery, with 66.8% and 53.5%, respectively.

CONCLUSION

The new window surgery demonstrates a lower pressure difference of anastomosis and higher energy conversion efficiency, which may be a better choice compared with the traditional surgery for S-TAPVC patient.

Physiological Fontan Procedure

Objective: The conventional Fontan circulation deviates the superior vena cava (SVC = 1/3 of the systemic venous return) toward the right lung (3/5 of total lung volume) and the inferior vena cava (IVC = 2/3 of the systemic venous return) toward the left lung (2/5 of total lung volume). A "physiological" Fontan deviating the SVC toward the left lung and the IVC toward the right lung was compared with the conventional setting by computational fluid dynamics, studying whether this setting achieves a more favorable hemodynamics than the conventional Fontan circulation. Materials and Methods: An in-silico 3D parametric model of the Fontan procedure was developed using idealized vascular geometries with invariant sizes of SVC, IVC, right pulmonary artery (RPA), and left pulmonary artery (LPA), steady inflow velocities at IVC and SVC, and constant equal outflow pressures at RPA and LPA. These parameters were set to perform finite-volume incompressible steady flow simulations, assuming a single-phase, Newtonian, isothermal, laminar blood flow. Numerically converged finite-volume mass and momentum flow balances determined the inlet pressures and the outflow rates. Numerical closed-path integration of energy fluxes across domain boundaries determined the flow energy loss rate through the Fontan circulation. The comparison evaluated: (1) mean IVC pressure; (2) energy loss rate; (3) kinetic energy maximum value throughout the domain volume. Results: The comparison of the physiological vs. conventional Fontan provided these results: (1) mean IVC pressure 13.9 vs. 14.1 mmHg (= 0.2 mmHg reduction); (2) energy loss rate 5.55 vs. 6.61 mW (= 16% reduction); (3) maximum kinetic energy 283 vs. 396 J/m3 (= 29% reduction). Conclusions: A more physiological flow distribution is accompanied by a reduction of mean IVC pressure and by substantial reductions of energy loss rate and of peak kinetic energy. The potential clinical impact of these hemodynamic changes in reducing the incidence and severity of the adverse long-term effects of the Fontan circulation, in particular liver failure and protein-losing enteropathy, still remains to be assessed and will be the subject of future work.

Evaluation of hemodynamic performance of total cavopulmonary connection (TCPC) with porous inserts

Infants born with univentricular heart disease undergo Fontan surgery to establish separate systemic and pulmonary circulations. This surgery results in better blood circulation across a single ventricle that supplies oxygenated blood to the body and passively returns venous blood to the lungs through the total cavopulmonary connection (TCPC). Reducing the pressure drop across the TCPC during Fontan circulation helps in reducing the work load of univentricular heart, and various designs have been proposed for this purpose. The goal of this work is to analyze the effect of placing a porous insert at an appropriate position in the pulmonary artery, on the pressure drop across the TCPC. A 3D computational model of a total TCPC connection provided with a porous insert is developed and solved by finite volume method, under assumptions of steady, laminar, and Newtonian flows. The effects of the porous medium properties-porosity and permeability-across the connection, are analyzed. Compared to the no-porous medium case, TCPC with the porous medium insert exhibits a maximum reduction of 27% in energy loss for the flow range studied. The porous medium used in TCPC connection lowers the energy dissipation by curtailing the flow recirculation zones across the connection. The influences of the diameter of the blood vessel, total cardiac output, and the thickness, permeability, and position of porous media on energy loss are analyzed. The criteria to select the porous medium properties and position for a given Fontan geometry are also determined.

Patient-specific in vitro models for hemodynamic analysis of congenital heart disease - Additive manufacturing approach

Non-invasive hemodynamic assessment of total cavopulmonary connection (TCPC) is challenging due to the complex anatomy. Additive manufacturing (AM) is a suitable alternative for creating patient-specific in vitro models for flow measurements using four-dimensional (4D) Flow MRI. These in vitro systems have the potential to serve as validation for computational fluid dynamics (CFD), simulating different physiological conditions. This study investigated three different AM technologies, stereolithography (SLA), selective laser sintering (SLS) and fused deposition modeling (FDM), to determine differences in hemodynamics when measuring flow using 4D Flow MRI. The models were created using patient-specific MRI data from an extracardiac TCPC. These models were connected to a perfusion pump circulating water at three different flow rates. Data was processed for visualization and quantification of velocity, flow distribution, vorticity and kinetic energy. These results were compared between each model. In addition, the flow distribution obtained in vitro was compared to in vivo. The results showed significant difference in velocities measured at the outlets of the models that required internal support material when printing. Furthermore, an ultrasound flow sensor was used to validate flow measurements at the inlets and outlets of the in vitro models. These results were highly correlated to those measured with 4D Flow MRI. This study showed that commercially available AM technologies can be used to create patient-specific vascular models for in vitro hemodynamic studies at reasonable costs. However, technologies that do not require internal supports during manufacturing allow smoother internal surfaces, which makes them better suited for flow analyses.

Hemodynamic study of TCPC using in vivo and in vitro 4D Flow MRI and numerical simulation

INTRODUCTION

Altered total cavopulmonary connection (TCPC) hemodynamics can cause long-term complications. Patient-specific anatomy hinders generalized solutions. 4D Flow MRI allows in vivo assessment, but not predictions under varying conditions and surgical approaches. Computational fluid dynamics (CFD) improves understanding and explores varying physiological conditions. This study investigated a combination of 4D Flow MRI and CFD to assess TCPC hemodynamics, accompanied with in vitro measurements as CFD validation. 4D Flow MRI was performed in extracardiac and atriopulmonary TCPC subjects. Data was processed for visualization and quantification of velocity and flow. Three-dimensional (3D) geometries were generated from angiography scans and used for CFD and a physical model construction through additive manufacturing. These models were connected to a perfusion system, circulating water through the vena cavae and exiting through the pulmonary arteries at two flow rates. Models underwent 4D Flow MRI and image processing. CFD simulated the in vitro system, applying two different inlet conditions from in vitro 4D Flow MRI measurements; no-slip was implemented at rigid walls. Velocity and flow were obtained and analyzed. The three approaches showed similar velocities, increasing proportionally with high inflow. Atriopulmonary TCPC presented higher vorticity compared to extracardiac at both inflow rates. Increased inflow balanced flow distribution in both TCPC cases. Atriopulmonary IVC flow participated in atrium recirculation, contributing to RPA outflow; at baseline, IVC flow preferentially traveled through the LPA. The combination of patient-specific in vitro and CFD allows hemodynamic parameter control, impossible in vivo. Physical models serve as CFD verification and fine-tuning tools.

HEMODYNAMIC EFFECTS OF THE ANASTOMOSES IN THE MODIFIED BLALOCK–TAUSSIG SHUNT: A NUMERICAL STUDY USING A 0D/3D COUPLING METHOD

The modified Blalock–Taussig (BT) shunt is a palliative surgery which can help the tetralogy of Fallot (TOF) patient increase the blood oxygen saturation by interposing a systemic-to-pulmonary artery shunt. Two typical anastomotic shapes are frequently used in clinical practice: the end-to-side (ETS) and the side-to-side (STS) anastomosis. This paper examines the hemodynamic influence of the anastomotic shape in the modified BT shunt. Three models with different anastomotic shapes were reconstructed. The ETS anastomoses were applied in the first model. For the innominate artery (IA) and the pulmonary artery (PA) in the second model, the ETS and the STS anastomosis were applied, respectively. Finally, the STS anastomoses were applied in the third model. The 0D/3D coupling method was used to perform a numerical simulation by coupling the three-dimensional (3D) artery model with a zero-dimensional (0D) lumped parameter model for the cardiovascular system. The simulation results showed that the perfusion into the left and right PA in Model 1 was unbalanced. Swirling flow appeared in the shunt in Model 3, but the shunt flow rate in Model 3 was lower. The ETS anastomosis at the PA may cause unbalanced blood perfusion into the left and right PA. Conversely, the STS anastomosis can make the blood perfusion more balanced. Otherwise, the STS anastomosis at the IA could generate a swirling flow in the shunt which may provide a better hemodynamic environment while decreasing the pulmonary perfusion.

Patient-specific assessment of cardiovascular function by combination of clinical data and computational model with applications to patients undergoing Fontan operation

The assessment of cardiovascular function is becoming increasingly important for the care of patients with single-ventricle defects. However, most measurement methods available in the clinical setting cannot provide a separate measure of cardiac function and loading conditions. In the present study, a numerical method has been proposed to compensate for the limitations of clinical measurements. The main idea was to estimate the parameters of a cardiovascular model by fitting model simulations to patient-specific clinical data via parameter optimization. Several strategies have been taken to establish a well-posed parameter optimization problem, including clinical data-matched model development, parameter selection based on an extensive sensitivity analysis, and proper choice of parameter optimization algorithm. The numerical experiments confirmed the ability of the proposed parameter optimization method to uniquely determine the model parameters given an arbitrary set of clinical data. The method was further tested in four patients undergoing the Fontan operation. Obtained results revealed a prevalence of ventricular abnormalities in the patient cohort and at the same time demonstrated the presence of marked inter-patient differences and preoperative to postoperative changes in cardiovascular function. Because the method allows a quick assessment and makes use of clinical data available in clinical practice, its clinical application is promising.

Hemodynamic performance of the Fontan circulation compared with a normal biventricular circulation: a computational model study

The physiological limitations of the Fontan circulation have been extensively addressed in the literature. Many studies emphasized the importance of pulmonary vascular resistance in determining cardiac output (CO) but gave little attention to other cardiovascular properties that may play considerable roles as well. The present study was aimed to systemically investigate the effects of various cardiovascular properties on clinically relevant hemodynamic variables (e.g., CO and central venous pressure). To this aim, a computational modeling method was employed. The constructed models provided a useful tool for quantifying the hemodynamic effects of any cardiovascular property of interest by varying the corresponding model parameters in model-based simulations. Herein, the Fontan circulation was studied compared with a normal biventricular circulation so as to highlight the unique characteristics of the Fontan circulation. Based on a series of numerical experiments, it was found that 1) pulmonary vascular resistance, ventricular diastolic function, and systemic vascular compliance play a major role, while heart rate, ventricular contractility, and systemic vascular resistance play a secondary role in the regulation of CO in the Fontan circulation; 2) CO is nonlinearly related to any single cardiovascular property, with their relationship being simultaneously influenced by other cardiovascular properties; and 3) the stability of central venous pressure is significantly reduced in the Fontan circulation. The findings suggest that the hemodynamic performance of the Fontan circulation is codetermined by various cardiovascular properties and hence a full understanding of patient-specific cardiovascular conditions is necessary to optimize the treatment of Fontan patients.

Transient hemodynamic changes upon changing a BCPA into a TCPC in staged Fontan operation: a computational model study

The clinical benefits of the Fontan operation in treating single-ventricle defects have been well documented. However, perioperative mortality or morbidity remains a critical problem. The purpose of the present study was to identify the cardiovascular factors that dominate the transient hemodynamic changes upon the change of a bidirectional cavopulmonary (Glenn) anastomosis (BCPA) into a total cavopulmonary connection (TCPC). For this purpose, two computational models were constructed to represent, respectively, a single-ventricle circulation with a BCPA and that with a TCPC. A series of model-based simulations were carried out to quantify the perioperative hemodynamic changes under various cardiovascular conditions. Obtained results indicated that the presence of a low pulmonary vascular resistance and/or a low lower-body vascular resistance is beneficial to the increase in transpulmonary flow upon the BCPA to TCPC change. Moreover, it was found that ventricular diastolic dysfunction and mitral valve regurgitation, despite being well-known risk factors for poor postoperative outcomes, do not cause a considerable perioperative reduction in transpulmonary flow. The findings may help physicians to assess the perioperative risk of the TCPC surgery based on preoperative measurement of cardiovascular function.