Optimization of inflow waveform phase-difference for minimized total cavopulmonary power loss

The effect of respiration-driven flow waveforms on hemodynamic metrics used in Fontan surgical planning

OBJECTIVE

Poor total cavopulmonary connection (TCPC) hemodynamics have been hypothesized to be associated with long-term complications in Fontan patients. Image-based Fontan surgical planning has shown great potential as a clinical tool because it can pre-operatively evaluate patient-specific hemodynamics. Current surgical planning paradigms commonly utilize cardiac-gated phase contrast magnetic resonance (MR) imaging to acquire vessel flows. These acquisitions are often taken under breath-held (BH) conditions and ignore the effect of respiration on blood flow waveforms. This study investigates the effect of respiration-driven flow waveforms on patient-specific hemodynamics using real-time MR acquisitions.

METHODS

Patient-specific TCPCs were reconstructed from cardiovascular MR images. Real-time phase contrast MR images were acquired under both free-breathing (FB) and breath-held conditions for 9 patients. Numerical simulations were employed to assess flow structures and hemodynamics used in Fontan surgical planning including hepatic flow distribution (HFD) and indexed power loss (iPL), which were then compared between FB and BH conditions.

RESULTS

Differences in TCPC flow structures between FB and BH conditions were observed throughout the respiratory cycle. However, the average differences (BH - FB values for each patient, which are then averaged) in iPL and HFD between these conditions were 0.002 ± 0.011 (p = 0.40) and 1 ± 3% (p = 0.28), respectively, indicating no significant difference in clinically important hemodynamic metrics.

CONCLUSIONS

Respiration affects blood flow waveforms and flow structures, but might not significantly influence the values of iPL or HFD. Therefore, breath-held MR acquisition can be adequate for Fontan surgical planning when focusing on iPL and HFD.

In vitro validation of a self-driving aortic-turbine venous-assist device for Fontan patients

Background

Palliative repair of single ventricle defects involve a series of open-heart surgeries where a single-ventricle (Fontan) circulation is established. As the patient ages, this paradoxical circulation gradually fails, because of its high venous pressure levels. Reversal of the Fontan paradox requires an extra subpulmonic energy that can be provided through mechanical assist devices. The objective of this study was to evaluate the hemodynamic performance of a totally implantable integrated aortic-turbine venous-assist (iATVA) system, which does not need an external drive power and maintains low venous pressure chronically, for the Fontan circulation.

Methods

Blade designs of the co-rotating turbine and pump impellers were developed and 3 prototypes were manufactured. After verifying the single-ventricle physiology at a pulsatile in vitro circuit, the hemodynamic performance of the iATVA system was measured for pediatric and adult physiology, varying the aortic steal percentage and circuit configurations. The iATVA system was also tested at clinical off-design scenarios.

Results

The prototype iATVA devices operate at approximately 800 revolutions per minute and extract up to 10% systemic blood from the aorta to use this hydrodynamic energy to drive a blood turbine, which in turn drives a mixed-flow venous pump passively. By transferring part of the available energy from the single-ventricle outlet to the venous side, the iATVA system is able to generate up to approximately 5 mm Hg venous recovery while supplying the entire caval flow.

Conclusions

Our experiments show that a totally implantable iATVA system is feasible, which will eliminate the need for external power for Fontan mechanical venous assist and combat gradual postoperative venous remodeling and Fontan failure.

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.

Respiratory Effects on Fontan Circulation During Rest and Exercise Using Real-Time Cardiac Magnetic Resonance Imaging

BACKGROUND

It is known that respiration modulates cavopulmonary flows, but little data compare mean flows under breath-holding and free-breathing conditions to isolate the respiratory effects and effects of exercise on the respiratory modulation.

METHODS

Real-time phase-contrast magnetic resonance combined with a novel method to track respiration on the same image acquisition was used to investigate respiratory effects on Fontan caval and aortic flows under breath-holding, free-breathing, and exercise conditions. Respiratory phasicity indices that were based on beat-averaged flow were used to quantify the respiratory effect.

RESULTS

Flow during inspiration was substantially higher than expiration under the free-breathing and exercise conditions for both inferior vena cava (inspiration/expiration: 1.6 ± 0.5 and 1.8 ± 0.5, respectively) and superior vena cava (inspiration/expiration: 1.9 ± 0.6 and 2.6 ± 2.0, respectively). Changes from rest to exercise in the respiratory phasicity index for these vessels further showed the impact of respiration. Total systemic venous flow showed no significant statistical difference between the breath-holding and free-breathing conditions. In addition, no substantial difference was found between the descending aorta and inferior vena cava mean flows under either resting or exercise conditions.

CONCLUSIONS

This study demonstrated that inferior vena cava and superior vena cava flow time variance is dominated by respiratory effects, which can be detected by the respiratory phasicity index. However, the minimal respiration influence on net flow validates the routine use of breath-holding techniques to measure mean flows in Fontan patients. Moreover, the mean flows in the inferior vena cava and descending aorta are interchangeable.

Computational Analysis of the Mechanical Behaviors of Hemiarch and Total Arch Replacements

Surgery for aortic dissections or aneurysms can be extended into the aortic arch by hemiarch replacement (HAR) or total arch replacement (TAR). Although cardiovascular surgeons have been performing HAR and TAR for decades, the mechanical properties of HAR and TAR are not well understood. This study investigates the mechanical behaviors and stress distributions in HAR and TAR using a hybrid fluid-structure interaction analysis that combines computational fluid dynamics and structural static analysis. Geometrical information on the aortas of 11 subjects was extracted from contrast-enhanced computed tomography (CT) scan data. The CT images were imported into medical image processing software to reconstruct 3D models of the aortas. A 3D finite element model was employed to simulate aortas that receive HAR or TAR. The deformation of the great vessels and the stress distributions at both the vessels and the aortic grafts were calculated. The numerical results revealed that the aortas following TAR exhibited a lower level of stress than those following HAR. Higher stresses may cause arterial wall injury and increase the risk of rupture. Finite element analysis of the aortas and the aortic grafts provides useful information that helps physicians better understand the potential problems that may arise after various surgical procedures.

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.

Dynamic Mode Decomposition of Fontan Hemodynamics in an Idealized Total Cavopulmonary Connection

Univentricular heart disease is the leading cause of death from any birth defect in the first year of life. Typically, patients have to undergo three open heart surgical procedures within the first few years of their lives to eventually directly connect the superior and inferior vena cavae to the left and right pulmonary arteries forming the Total Cavopulmonary Connection or TCPC. The end result is a weak circulation where the single working ventricle pumps oxygenated blood to the body and de-oxygenated blood flows passively through the TCPC into the lungs. The fluid dynamics of the TCPC junction involve confined impinging jets resulting in a highly unstable flow, significant mechanical energy dissipation, and undesirable pressure loss. Understanding and predicting such flows is important for improving the surgical procedure and for the design of mechanical cavopulmonary assist devices. In this study, Dynamic Mode Decomposition (DMD) is used to analyze previously obtained Stereoscopic Particle Imaging Velocimetry (SPIV) data and Large Eddy Simulation (LES) results for an idealized TCPC. Analysis of the DMD modes from the SPIV and LES serve to both highlight the unsteady vortical dynamics and the qualitative agreement between measurements and simulations.

Regional variation of wall shear stress in ascending thoracic aortic aneurysms

The development of an ascending thoracic aortic aneurysm is likely caused by excessive hemodynamic loads exerted on the aneurysmal wall. Computational fluid-dynamic analyses were performed on patient-specific ascending thoracic aortic aneurysms obtained from patients with either bicuspid aortic valve or tricuspid aortic valve to evaluate hemodynamic and wall shear parameters, imparting aneurysm enlargement. Results showed an accelerated flow along the outer aortic wall with helical flow in the aneurysm center for bicuspid aortic valve ascending thoracic aortic aneurysms. In a different way, tricuspid aortic valve ascending thoracic aortic aneurysms exhibited normal systolic flow without substantial secondary pattern. Analysis of wall shear parameters evinced a high and locally varying wall shear stress on the outer aortic wall and high temporal oscillations in wall shear stress (oscillatory shear index) on either left or right side of aneurysmal aorta. These findings may explain the asymmetric dilatation typically observed in ascending thoracic aortic aneurysms. Simulations of a hypertensive scenario revealed an increase in wall shear stress upon 44% compared to normal systemic pressure models. Computational fluid-dynamics–based analysis may allow identification of wall shear parameters portending aneurysm dilatation and hence guide preventative intervention.

Computational fluid dynamics simulation to evaluate aortic coarctation gradient with contrast-enhanced CT

Coarctation of aorta (CoA) is a narrowing of the aorta leading to a pressure gradient (ΔP) across the coarctation, increased afterload and reduced peripheral perfusion pressures. Indication to invasive treatment is based on values of maximal (systolic) trans-coarctation ΔP. A computational fluid dynamic (CFD) approach is herein presented for the non-invasive haemodynamic assessment of ΔP across CoA. Patient-specific CFD simulations were created from contrast-enhanced computed tomography (CT) and appropriate flow boundary conditions. Computed ΔP was validated with invasive intravascular trans-CoA pressure measurements. Haemodynamic indices, including pressure loss coefficient (PLc), time-averaged wall shear stress (TAWSS) and oscillatory shear index (OSI), were also quantified. CFD-estimated ΔP values were comparable to the invasive ones. Moreover, the aorta proximal to CoA was exposed to altered TAWSS and OSI suggesting hypertension. PLc was found as a further geometric marker of CoA severity. Finally, CFD-estimated ΔP confirmed a significant reduction after percutaneous balloon dilatation and stenting of the CoA in one patient (e.g. from ΔP∼52 mmHg to ΔP∼3 mmHg). The validation of the ΔP computations with catheterisation measurements suggests that CFD simulation, based on CT-derived anatomical data, is a useful tool to readily quantify CoA severity.

Mock circulatory system of the Fontan circulation to study respiration effects on venous flow behavior

We describe an in vitro model of the Fontan circulation with respiration to study subdiaphragmatic venous flow behavior. The venous and arterial connections of a total cavopulmonary connection (TCPC) test section were coupled with a physical lumped parameter (LP) model of the circulation. Intrathoracic and subdiaphragmatic pressure changes associated with normal breathing were applied. This system was tuned for two patients (5 years, 0.67 m2; 10 years, 1.2 m2) to physiological values. System function was verified by comparison to the analytical model on which it was based and by consistency with published clinical measurements. Overall, subdiaphragmatic venous flow was influenced by respiration. Flow within the arteries and veins increased during inspiration but decreased during expiration, with retrograde flow in the inferior venous territories. System pressures and flows showed close agreement with the analytical LP model (p < 0.05). The ratio of the flow rates occurring during inspiration to expiration were within the clinical range of values reported elsewhere. The approach used to set up and control the model was effective and provided reasonable comparisons with clinical data.

Fontan conversion templates: patient-specific hemodynamic performance of the lateral tunnel versus the intraatrial conduit with fenestration

Intraatrial-conduit Fontan is considered a modification of both extracardiac and lateral-tunnel Fontan. In this study, the patient-specific hemodynamic performance of intraatrial-conduit and lateral-tunnel Fontan with fenestration, considered as conversion templates, was investigated based on the authors' patient cohort. Pulsatile computational fluid dynamics simulations were performed using patient-specific models of intraatrial-conduit and lateral-tunnel Fontan patients. Real-time "simultaneous" inferior and superior vena cava, pulmonary artery, and fenestration flow waveforms were acquired from ultrasound. Multiple hemodynamic performance indices were investigated, with particular focus on evaluation of the pulsatile flow performance. Power loss inside the lateral-tunnel Fontan appeared to be significantly higher than with the intraatrial-conduit Fontan for patient-specific cardiac output and normalized connection size. Inclusion of the 4-mm fenestration at a 0.24 L/min mean flow resulted in a lower cavopulmonary pressure gradient and less time-averaged power loss for both Fontan connections. Flow structures within the intraatrial conduit were notability more uniform than within the lateral tunnel. Hepatic flow majorly favored the left lung in both surgical connections: conversion from lateral-tunnel to intraatrial-conduit Fontan resulted in better hemodynamics with less power loss, a lower pressure gradient, and fewer stagnant flow zones along the conduit. This patient-specific computational case study demonstrated superior hemodynamics of intraatrial-conduit Fontan over those of lateral-tunnel Fontan with or without fenestration and improved performance after conversion of the lateral tunnel to the intraatrial conduit. The geometry-specific effect of the nonuniform hepatic flow distribution may motivate new rationales for the surgical design.

Numerical and experimental investigation of pulsatile hemodynamics in the total cavopulmonary connection

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.

Effect of flow pulsatility on modeling the hemodynamics in the total cavopulmonary connection

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.

Pulsatile venous waveform quality affects the conduit performance in functional and "failing" Fontan circulations

OBJECTIVE

To investigate the effect of pulsatility of venous flow waveform in the inferior and superior caval vessels on the performance of functional and "failing" Fontan patients based on two primary performance measures - the conduit power loss and the distribution of inferior caval flow (hepatic factors) to the lungs.

METHODS

Doppler angiography flows were acquired from two typical extra-cardiac conduit "failing" Fontan patients, aged 13 and 25 years, with ventricle dysfunction. Using computational fluid dynamics, haemodynamic efficiencies of "failing", functional, and in vitro-generated mechanically assisted venous flow waveforms were evaluated inside an idealised total cavopulmonary connection with a caval offset. To investigate the effect of venous pulsatility alone, cardiac output was normalised to 3 litres per minute in all cases. To quantify the pulsatile behaviour of venous flows, two new performance indices were suggested.

RESULTS

Variations in the pulsatile content of venous waveforms altered the conduit efficiency notably. High-frequency and low-amplitude oscillations lowered the pulsatile component of the power losses in "failing" Fontan flow waveforms. Owing to the offset geometry, hepatic flow distribution depended strongly on the ratio of time-dependent caval flows and the pulsatility content rather than mixing at the junction. "Failing" Fontan flow waveforms exhibited less balanced hepatic flow distribution to lungs.

CONCLUSIONS

The haemodynamic efficiency of single-ventricle circulation depends strongly on the pulsatility of venous flow waveforms. The proposed performance indices can be calculated easily in the clinical setting in efforts to better quantify the energy efficiency of Fontan venous waveforms in pulsatile settings.

Computer-Aided Patient-Specific Coronary Artery Graft Design Improvements Using CFD Coupled Shape Optimizer

This study aims to (i) demonstrate the efficacy of a new surgical planning framework for complex cardiovascular reconstructions, (ii) develop a computational fluid dynamics (CFD) coupled multi-dimensional shape optimization method to aid patient-specific coronary artery by-pass graft (CABG) design and, (iii) compare the hemodynamic efficiency of the sequential CABG, i.e., raising a daughter parallel branch from the parent CABG in patient-specific 3D settings. Hemodynamic efficiency of patient-specific complete revascularization scenarios for right coronary artery (RCA), left anterior descending artery (LAD), and left circumflex artery (LCX) bypasses were investigated in comparison to the stenosis condition. Multivariate 2D constraint optimization was applied on the left internal mammary artery (LIMA) graft, which was parameterized based on actual surgical settings extracted from 2D CT slices. The objective function was set to minimize the local variation of wall shear stress (WSS) and other hemodynamic indices (energy dissipation, flow deviation angle, average WSS, and vorticity) that correlate with performance of the graft and risk of re-stenosis at the anastomosis zone. Once the optimized 2D graft shape was obtained, it was translated to 3D using an in-house "sketch-based" interactive anatomical editing tool. The final graft design was evaluated using an experimentally validated second-order non-Newtonian CFD solver incorporating resistance based outlet boundary conditions. 3D patient-specific simulations for the healthy coronary anatomy produced realistic coronary flows. All revascularization techniques restored coronary perfusions to the healthy baseline. Multi-scale evaluation of the optimized LIMA graft enabled significant wall shear stress gradient (WSSG) relief (~34%). In comparison to original LIMA graft, sequential graft also lowered the WSSG by 15% proximal to LAD and diagonal bifurcation. The proposed sketch-based surgical planning paradigm evaluated the selected coronary bypass surgery procedures based on acute hemodynamic readjustments of aorta-CA flow. This methodology may provide a rational to aid surgical decision making in time-critical, patient-specific CA bypass operations before in vivo execution.