Effects of pulmonary afterload on the hemodynamics after the hemi-Fontan procedure

Hemi-Fontan and bidirectional Glenn operations result in flow-mediated viscous energy loss at the time of stage II palliation

Background

Superior cavopulmonary connection (SCPC) for stage II palliation of hypoplastic left heart syndrome (HLHS) is achieved most frequently by either a bidirectional Glenn (BDG) or hemi-Fontan (HF) operation. The comparison of flow hemodynamic efficiency at the region of surgical reconstruction and in proximal pulmonary arteries has been evaluated primarily using computational modeling techniques with conflicting reports. The purpose of this descriptive study was to compare flow hemodynamics following stage II (BDG vs HF) using 4-dimensional flow magnetic resonance imaging (4D-Flow MRI) with particular focus on flow-mediated viscous energy loss (EL') under matched hemodynamic conditions.

Methods

Patients with hypoplastic left heart syndrome (HLHS) who underwent either HF or BDG as part of stage II palliation underwent pre-Fontan 4D-Flow MRI. Patients were matched by the pulmonary vascular resistance index, net superior vena cava (SVC) flow, right pulmonary artery (RPA) and left pulmonary artery (LPA) size, and age. Maximum EL' throughout the cardiac cycle was calculated along the SVC-RPA and SVC-LPA tracts.

Results

Eight patients who underwent HF as part of their stage II single ventricle palliation were matched with 8 patients who underwent BDG. There were no differences between the 2 groups in median volumetric indices, including end-diastolic volume (P = .278) and end-systolic volume (P = .213). Moreover, no differences were observed in ejection fraction (P = .091) and cardiac index (P = .324). There also were no differences in peak EL' measured along the SVC-RPA tract (median, 0.05 mW for HF vs 0.04 mW for BDG; P = .365) or along the SVC-LPA tract (median, 0.05 mW vs 0.04 mW; P = .741).

Conclusions

The second stage of surgical palliation of HLHS using either HF or BDG results in similar flow-mediated viscous energy loss throughout the SCPC junction. 4D-Flow MRI and computational methods should be applied together to investigate flow hemodynamic patterns throughout the Fontan palliation and overall efficiency of the Fontan circuit.

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.

NUMERICAL STUDY OF BIDIRECTIONAL GLENN WITH UNILATERAL PULMONARY ARTERY STENOSIS

Purpose: Hypoplastic left heart syndrome (HLHS) is a congenital heart disease and is usually associated with pulmonary artery stenosis. The superior vena cava-to-pulmonary artery (bidirectional Glenn) shunt is used primarily as a staging procedure to the total cava-to-pulmonary connection for single-ventricle complex. When HLHS coexists with pulmonary artery stenosis, the surgeons then face a multiple problem. This leads to high demand of optimized structure of Glenn surgery. The objective of this article is to investigate the influence of various anastomotic structures and the direction of superior vena cava (SVC) in Glenn on hemodynamics under pulse inflow conditions and try to find an optimal structure of SVC in Glenn surgery with unilateral pulmonary artery stenosis.

Method: First, 3D patient-specific models were constructed from medical images of a HLHS patient before any surgery by using the commercial software Mimics, and another software Free-form was used to deform the reconstructed models in the computer. Four 3D patient-specific Glenn models were constructed: model-1 (normal Glenn), model-2 (lean the SVC back to the stenotic pulmonary artery), model-3 (lean the SVC towards the stenotic pulmonary artery), model-4 (add patch at junction of the SVC toward stenosis at pulmonary artery). Second, a lumped parameter model (LPM) was established to predict boundary conditions for computational fluid dynamics (CFD). In addition, numerical simulations were conducted using CFD through the finite volume method. Finally, hemodynamic parameters were obtained and evaluated.

Results: It was showed that model-4 have relatively balanced vena cava blood perfusion into the left pulmonary artery (LPA) and right pulmonary artery (RPA), this may be due to less helical flow and the patch at junction of the SVC. Near stenosis of pulmonary artery, model-4 performed with the higher wall shear stress (WSS), which would benefit endothelial cell function and gene expression. In addition, results showed that model-4 performed with the lower oscillatory shear index (OSI) and wall shear stress gradient (WSSG), which would decrease the opportunity of vascular intimal hyperplasia.

Conclusion: It is benefited that surgeons adds patch at junction of the SVC towards stenosis at pulmonary artery. These results can impact the surgical design and planning of the Glenn surgery with unilateral pulmonary artery stenosis.

Three-dimensional simulations in Glenn patients: clinically based boundary conditions, hemodynamic results and sensitivity to input data

While many congenital heart defects can be treated without significant long term sequelae, some achieve successful palliation as their definitive endpoints. The single-ventricle defect is one such defect and leaves the child with only one operational ventricle, requiring the systemic and the pulmonary circulations to be placed in series through several operations performed during early childhood. Numerical simulations may be used to investigate these hemodynamic conditions and their relation to post-operative sequelae; however, they rely heavily on boundary condition prescription. In this study, we investigate the impact of hemodynamic input data uncertainties on simulation results. Imaged-based patient-specific models of the multi-branched pulmonary arteries and superior vena cava were built for five cavopulmonary connection (i.e. Glenn) patients. Magnetic resonance imaging and catheterization data were acquired for each patient prior to their Fontan surgery. Inflow and outflow boundary conditions were constructed to match available clinical data and resulted in the development of a framework to incorporate these types of clinical data into patient-specific simulations. Three-dimensional computational fluid dynamics simulations were run and hemodynamic indicators were computed. Power loss was low (and efficiency very high) and a linear correlation was found between power loss and cardiac index among the five patients. Other indicators such as low wall shear stress were considered to better characterize these patients. Flow was complex and oscillatory near the anastomosis, and laminar in the smaller branches. While common trends were seen among patients, results showed differences among patients, especially in the 3D maps, strengthening the importance of patient-specific simulations. A sensitivity analysis was performed to investigate the impact of input data (clinical and modeling) to construct boundary conditions on several indicators. Overall, the sensitivity of the output indicators to the input data was small but non-negligible. The sensitivity of commonly used hemodynamic indicators to compare patients is discussed in this context. Power efficiency was much more sensitive to pressure variation than power loss. To increase the precision of such indicators, mean flow split between right and left lungs needs to be measured with more accuracy with higher priority than refining the model of how the flow is distributed on average among the smaller branches. Although ± 10% flow split imprecision seemed reasonable in terms of patient comparison, this study suggests that the common practice of imposing a right pulmonary artery/left pulmonary artery flow split of 55%/45% when performing patient specific simulations should be avoided. This study constitutes a first step towards understanding the hemodynamic differences between pre- and post Fontan surgery, predicting these differences, and evaluating surgical outcomes based on preoperative data.

In vitro hemodynamic investigation of the embryonic aortic arch at late gestation

This study focuses on the dynamic flow through the fetal aortic arch driven by the concurrent action of right and left ventricles. We created a parametric pulsatile computational fluid dynamics (CFD) model of the fetal aortic junction with physiologic vessel geometries. To gain a better biophysical understanding, an in vitro experimental fetal flow loop for flow visualization was constructed for identical CFD conditions. CFD and in vitro experimental results were comparable. Swirling flow during the acceleration phase of the cardiac cycle and unidirectional flow following mid-deceleration phase were observed in pulmonary arteries (PA), head-neck vessels, and descending aorta. Right-to-left (oxygenated) blood flowed through the ductus arteriosus (DA) posterior relative to the antegrade left ventricular outflow tract (LVOT) stream and resembled jet flow. LVOT and right ventricular outflow tract flow mixing had not completed until approximately 3.5 descending aorta diameters downstream of the DA insertion into the aortic arch. Normal arch model flow patterns were then compared to flow patterns of four common congenital heart malformations that include aortic arch anomalies. Weak oscillatory reversing flow through the DA junction was observed only for the Tetralogy of Fallot configuration. PA and hypoplastic left heart syndrome configurations demonstrated complex, abnormal flow patterns in the PAs and head-neck vessels. Aortic coarctation resulted in large-scale recirculating flow in the aortic arch proximal to the DA. Intravascular flow patterns spatially correlated with abnormal vascular structures consistent with the paradigm that abnormal intravascular flow patterns associated with congenital heart disease influence vascular growth and function.

Modeling the Fontan circulation: where we are and where we need to go

The Fontan procedure and its subsequent modifications over the past 30 years can be described as a class of surgical procedures for patients born with complex congenital heart disease exhibiting a single-ventricle physiology. The long-term outcome for children currently undergoing a Fontan procedure remains worrisome because of multiple late morbidities observed. Despite significant modeling efforts spanning three decades, improvements to the Fontan procedure have occurred without comprehensive validation from these modeling studies. Careful examination shows that modeling studies to date offer only a "glimpse through a keyhole" into understanding and modeling a representative range of the variations in anatomy and physiology that exist in Fontan patients. Suggestions for future investigations are provided.

Toward optimal hemodynamics: computer modeling of the Fontan circuit

The construction of efficient designs with minimal energy losses is especially important for cavopulmonary connections. The science of computational fluid dynamics has been increasingly used to study the hemodynamic performance of surgical operations. Three-dimensional computer models can be accurately constructed of typical cavopulmonary connections used in clinical practice based on anatomic data derived from magnetic resonance scans, angiocardiograms, and echocardiograms. Using these methods, the hydraulic performance of the hemi-Fontan, bidirectional Glenn, and a variety of types of completion Fontan operations can be evaluated and compared. This methodology has resulted in improved understanding and design of these surgical operations.

Congenital Heart Disease Associated Pulmonary Arterial Hypertension

Pulmonary arterial hypertension associated with congenital heart disease, although common (15%-30%) in all-comers with congenital heart disease, is variable in terms of clinical manifestations, severity of associated pulmonary arterial hypertension, and response to therapy and outcomes (depending on lesion anatomy, pulmonary circulation flows and pressures, and presence and timings of surgeries). Evaluation includes imaging and catheterization. Surgical or another anatomic correction may be desirable after rigorous preinterventional assessment. Patients who are not surgical candidates or who fail to improve early or late after surgery may have the potential to respond to idiopathic pulmonary arterial hypertension therapies. Lung or heart/lung transplantation remains an option for selected recalcitrant patients.

Wall shear stress is the primary mechanism of energy loss in the Fontan connection

Long-term outcome following the Fontan operation may be affected by the amount of energy lost as blood flows through the anastomosis geometry. A method for detailed quantification of energy loss is applied to computational simulations of the flow in an atriopulmonary and a total cavopulmonary model. Five types of flow (near wall, slow recirculation, medium speed vortices, collision, and streamlined flow) are identified and their energy losses quantified. The presence of recirculation regions decreases the efficiency of the atriopulmonary model, and a region of increased energy loss is seen in the collision region in the total cavopulmonary model. However, the most significant energy loss is through wall shear stress, which is maximal in areas where there is rapid, near wall flow.

Efficiency differences in computational simulations of the total cavo-pulmonary circulation with and without compliant vessel walls

The Fontan operation is a palliative surgical procedure performed on children born with congenital defects of the heart that have yielded only a single functioning ventricle. The total cavo-pulmonary connection (TCPC) is the most popular variant of the Fontan procedure. The objective of the study was to quantify and compare the efficiency of numerical models of the TCPC with rigid versus elastic vessel wall models. The pressure drop and power loss through both type TCPC models was measured. Significant differences in efficiencies exist between rigid versus elastic numerical models. We have shown incorporating elasticity into numerical models of the total cavo-pulmonary connection is important when determining circuit efficiencies.

Computational simulations of the total cavo-pulmonary connection: insights in optimizing numerical solutions

The Fontan procedure is a palliative surgical technique that is used to treat patients with congenital heart defects that include complex lesions such as those with a hypoplastic ventricle. In vitro, in vivo, and computational models of a set of modifications to the Fontan procedure, called the total cavopulmonary connection (TCPC), have been developed. Using these modeling methods, attempts have been made at finding the most energy efficient TCPC circuit. Computational modeling has distinct advantages to other modeling methods. However, discrepancies have been found in validation studies of TCPC computational models. There is little in the literature available to help explain and correct for such discrepancies. Differences in computational results can occur when choosing between steady flow versus transient flow numerical solvers. In this study transient flow solver results were shown to be more consistent with results from previous TCPC in vitro experiments. Using a transient flow solver we found complex fluctuating flow patterns can exist with steady inflow boundary conditions in computational models of the TCPC. To date such findings have not been reported in the literature. Furthermore, our computational modeling results suggest fluctuating flow patterns as well as the magnitudes of these secondary flow structures diminish if the TCPC offset between vena cavae is increased or if flanged connections are added. An association was found between these modifications and improvements in TCPC circuit flow efficiencies. In summary, development of accurate computational simulations in the validation process is critical to efforts in finding the most efficient TCPC circuits, efforts aimed at potentially improving the long term outcome for Fontan patients.

Physics-Driven CFD Modeling of Complex Anatomical Cardiovascular Flows?A TCPC Case Study

Recent developments in medical image acquisition combined with the latest advancements in numerical methods for solving the Navier-Stokes equations have created unprecedented opportunities for developing simple and reliable computational fluid dynamics (CFD) tools for meeting patient-specific surgical planning objectives. However, for CFD to reach its full potential and gain the trust and confidence of medical practitioners, physics-driven numerical modeling is required. This study reports on the experience gained from an ongoing integrated CFD modeling effort aimed at developing an advanced numerical simulation tool capable of accurately predicting flow characteristics in an anatomically correct total cavopulmonary connection (TCPC). An anatomical intra-atrial TCPC model is reconstructed from a stack of magnetic resonance (MR) images acquired in vivo. An exact replica of the computational geometry was built using transparent rapid prototyping. Following the same approach as in earlier studies on idealized models, flow structures, pressure drops, and energy losses were assessed both numerically and experimentally, then compared. Numerical studies were performed with both a first-order accurate commercial software and a recently developed, second-order accurate, in-house flow solver. The commercial CFD model could, with reasonable accuracy, capture global flow quantities of interest such as control volume power losses and pressure drops and time-averaged flow patterns. However, for steady inflow conditions, both flow visualization experiments and particle image velocimetry (PIV) measurements revealed unsteady, complex, and highly 3D flow structures, which could not be captured by this numerical model with the available computational resources and additional modeling efforts that are described. Preliminary time-accurate computations with the in-house flow solver were shown to capture for the first time these complex flow features and yielded solutions in good agreement with the experimental observations. Flow fields obtained were similar for the studied total cardiac output range (1-3 1/min); however hydrodynamic power loss increased dramatically with increasing cardiac output, suggesting significant energy demand at exercise conditions. The simulation of cardiovascular flows poses a formidable challenge to even the most advanced CFD tools currently available. A successful prediction requires a two-pronged, physics-based approach, which integrates high-resolution CFD tools and high-resolution laboratory measurements.

Computational fluid dynamics simulations in realistic 3-D geometries of the total cavopulmonary anastomosis: the influence of the inferior caval anastomosis

Fluid dynamics of Total Cavo-Pulmonary Connection (TCPC) were studied in 3-D models based on real dimensions obtained by Magnetic Resonance (MR) images. Models differ in terms of shape (intra- or extra-cardiac conduit) and cross section (with or without patch enlargement) of the inferior caval (IVC) anastomosis connection. Realistic pulsatile flows were submitted to both the venae cavae, while porous portions were added at the end of the pulmonary arteries to reproduce the pulmonary afterload. The dissipated power and the flow distribution into the lungs were calculated at different values of pulmonary arteriolar resistances (PAR). The most important results are: i) power dissipation in different TCPC designs is influenced by the actual cross sectional area of the IVC anastomosis and ii) the inclusion of a patch minimizes the dissipated power (range 4-13 mW vs. 14-56 mW). Results also show that the perfusion of the right lung is between 15% and 30% of the whole IVC blood flow when the PAR are evenly distributed between the right and the left lung.

Computational fluid dynamics in the evaluation of hemodynamic performance of cavopulmonary connections after the Norwood procedure for hypoplastic left heart syndrome

OBJECTIVE

Computational fluid dynamics have been used to study the hemodynamic performance of surgical operations, resulting in improved design. Efficient designs with minimal energy losses are especially important for cavopulmonary connections. The purpose of this study was to compare hydraulic performance between the hemi-Fontan and bidirectional Glenn procedures, as well as the various types of completion Fontan operations.

METHODS

Three-dimensional models were constructed of typical hemi-Fontan and bidirectional Glenn operations according to anatomic data derived from magnetic resonance scans, angiocardiograms, and echocardiograms. Boundary conditions were imposed, and fluid dynamics were calculated from a mathematic code. Power losses, flow distribution to each lung, and pressures were measured at three predetermined levels of pulmonary arteriolar resistance. Models of the lateral tunnel, total cavopulmonary connection, and extracardiac conduit completion Fontan operations were constructed, and power losses, total flow distribution, vena caval and pulmonary arterial pressures, and flow distribution of inferior vena caval return were calculated.

RESULTS

The hemi-Fontan and bidirectional Glenn procedures performed nearly identically, with similar power losses and nearly equal flow distributions to each lung at all levels of pulmonary arteriolar resistance. However, the lateral tunnel Fontan procedure as performed after the hemi-Fontan operation had lower power losses (6.9 mW, pulmonary arteriolar resistance 3 units) than the total cavopulmonary connection (40.5 mW) or the extracardiac conduit (42.9 mW), although the inclusion of an enlargement patch toward the right in the total cavopulmonary connection was effective in reducing the difference (10.0 mW). Inferior vena caval flow to the right lung was 52% for the lateral tunnel, compared with 19%, 30%, 19%, and 15% for the total cavopulmonary connection, total cavopulmonary connection with right-sided enlargement patch, extracardiac conduit, and extracardiac conduit with a bevel to the left lung, respectively.

CONCLUSIONS

According to these methods, the hemi-Fontan and bidirectional Glenn procedures performed equally well, but important differences in energy losses and flow distribution were found after the completion Fontan procedures. The superior hydraulic performance of the lateral tunnel Fontan operation after the hemi-Fontan procedure relative to any other method may be due to closer to optimal caval offset achieved in the surgical reconstruction.