Ventricular Flow Field Visualization During Mechanical Circulatory Support in the Assisted Isolated Beating Heart

Left ventricular flow dynamics by cardiac imaging techniques in heart failure patients: state of the art

BACKGROUND

The evaluation of left ventricular (LV) flow dynamics is a novel approach to assessing LV function that goes beyond traditional metrics. This approach has been applied to patients with heart failure (HF), providing valuable insights that are discussed in this review, with the aim of enhancing our understanding of LV function in the context of the HF syndrome.

METHODS

The analysis of LV flow dynamics is typically conducted using ultrasound and magnetic resonance imaging (MRI) techniques, primarily including particle image velocimetry echocardiography, Vector Flow Imaging, HyperDoppler, and four-dimensional flow MRI. A variety of parameters can be obtained that describe the geometry of the LV vortex, vorticity, kinetic energy, energy dispersion, as well as the amplitude and direction of the hemodynamic forces within the LV cavity.

RESULTS

In normal subjects, vortex formation plays a crucial role in optimizing LV filling, diastolic-systolic coupling, and energy transfer during systolic ejection. In patients with HF, alterations in vortex structure and dynamics have been associated with both systolic and diastolic LV dysfunction, demonstrating the potential to diagnose early LV dysfunction. Furthermore, these alterations have been linked to LV remodeling and thrombus formation. Several studies have also explored intracardiac flow metrics as biomarkers for guiding HF treatments, including pharmacological interventions, cardiac resynchronization therapy, and LV assist devices.

CONCLUSIONS

Currently available data suggest that the evaluation of LV flow dynamics can have diagnostic and prognostic utility in HF. However, large-scale, multicenter, and prospective studies are needed, particularly to validate therapeutic implications.

Influence of LVAD Cannula Outflow Graft Flow Rate and Location on Fluid-Particle Interactions and Thrombi Distribution: A Primary Numerical Study

A left ventricular assist device (LVAD) supports hemodynamics in heart failure patients. To deepen the understanding of hemodynamic changes and the movement of thrombi in the aorta, we examined three distinct LVAD blood flow rates across two implantation sites using the theory of computational fluid dynamics. Our findings revealed the complex dynamics of blood flow during cardiac systole under various scenarios. We also analyzed thrombi residence time and flow probabilities into aortic branches. Simulation results indicate that thrombi distribution in the aorta is significantly influenced by the location of the LVAD outflow graft and the flow rate. When the LVAD outflow graft is implanted into the ascending aorta, higher flow rates may reduce the risk of cerebral thrombosis. However, lower flow rates may reduce the risk of cerebral thrombosis while it is implanted into the descending aorta. The study may offer valuable insights into the LVAD implantation about the risk of cerebrovascular embolism.

A Computational Investigation of the Effects of Temporal Synchronization of Left Ventricular Assist Device Speed Modulation with the Cardiac Cycle on Intraventricular Hemodynamics

Patients with advanced heart failure are implanted with a left ventricular assist device (LVAD) as a bridge-to-transplantation or destination therapy. Despite advances in pump design, the risk of stroke remains high. LVAD implantation significantly alters intraventricular hemodynamics, where regions of stagnation or elevated shear stresses promote thrombus formation. Third generation pumps incorporate a pulsatility mode that modulates rotational speed of the pump to enhance in-pump washout. We investigated how the timing of the pulsatility mode with the cardiac cycle affects intraventricular hemodynamic factors linked to thrombus formation. Computational fluid dynamics simulations with Lagrangian particle tracking to model platelet behavior in a patient-specific left ventricle captured altered intraventricular hemodynamics due to LVAD implantation. HeartMate 3 incorporates a pulsatility mode that modulates the speed of the pump every two seconds. Four different timings of this pulsatility mode with respect to the cardiac cycle were investigated. A strong jet formed between the mitral valve and inflow cannula. Blood stagnated in the left ventricular outflow tract beneath a closed aortic valve, in the near-wall regions off-axis of the jet, and in a large counterrotating vortex near the anterior wall. Computational results showed good agreement with particle image velocimetry results. Synchronization of the pulsatility mode with peak systole decreased stasis, reflected in the intraventricular washout of virtual contrast and Lagrangian particles over time. Temporal synchronization of HeartMate 3 pulsatility with the cardiac cycle reduces intraventricular stasis and could be beneficial for decreasing thrombogenicity.

Benchtop Models of Patient-Specific Intraventricular Flow During Heart Failure and LVAD Support

The characterization of intraventricular flow is critical to evaluate the efficiency of fluid transport and potential thromboembolic risk but challenging to measure directly in advanced heart failure (HF) patients with left ventricular assist device (LVAD) support. The study aims to validate an in-house mock loop (ML) by simulating specific conditions of HF patients with normal and prosthetic mitral valves (MV) and LVAD patients with small and dilated left ventricle volumes, then comparing the flow-related indices result of vortex parameters, residence time (RT), and shear-activation potential (SAP). Patient-specific inputs for the ML studies included heart rate, end-diastolic and end-systolic volumes, ejection fraction, aortic pressure, E/A ratio, and LVAD speed. The ML effectively replicated vortex development and circulation patterns, as well as RT, particularly for HF patient cases. The LVAD velocity fields reflected altered flow paths, in which all or most incoming blood formed a dominant stream directing flow straight from the mitral valve to the apex. RT estimation of patient and ML compared well for all conditions, but SAP was substantially higher in the LVAD cases of the ML. The benchtop system generated comparable and reproducible hemodynamics and fluid dynamics for patient-specific conditions, validating its reliability and clinical relevance. This study demonstrated that ML is a suitable platform to investigate the fluid dynamics of HF and LVAD patients and can be utilized to investigate heart-implant interactions.

The role of innovative modeling and imaging techniques in improving outcomes in patients with LVAD

Heart failure remains a significant cause of mortality in the United States and around the world. While organ transplantation is acknowledged as the gold standard treatment for end stage heart failure, supply is limited, and many patients are treated with left ventricular assist devices (LVADs). LVADs extend and improve patients' lives, but they are not without their own complications, particularly the hemocompatibility related adverse events (HRAE) including stroke, bleeding and pump thrombosis. Mainstream imaging techniques currently in use to assess appropriate device function and troubleshoot complications, such as echocardiography and cardiac computed tomography, provide some insight but do not provide a holistic understanding of pump induced flow alterations that leads to HRAEs. In contrast, there are technologies restricted to the benchtop-such as computational fluid dynamics and mock circulatory loops paired with methods like particle image velocimetry-that can assess flow metrics but have not been optimized for clinical care. In this review, we outline the potential role and current limitations of converging available technologies to produce novel imaging techniques, and the potential utility in evaluating hemodynamic flow to determine whether LVAD patients may be at higher risk of HRAEs. This addition to diagnostic and monitoring capabilities could improve prevention and treatment of LVAD-induced complications in heart failure patients.

Impact of dobutamine stress on diastolic energetic efficiency of healthy left ventricle: an in vivo kinetic energy analysis

The total kinetic energy (KE) of blood can be decomposed into mean KE (MKE) and turbulent KE (TKE), which are associated with the phase-averaged fluid velocity field and the instantaneous velocity fluctuations, respectively. The aim of this study was to explore the effects of pharmacologically induced stress on MKE and TKE in the left ventricle (LV) in a cohort of healthy volunteers. 4D Flow MRI data were acquired in eleven subjects at rest and after dobutamine infusion, at a heart rate that was ∼60% higher than the one in rest conditions. MKE and TKE were computed as volume integrals over the whole LV and as data mapped to functional LV flow components, i.e., direct flow, retained inflow, delayed ejection flow and residual volume. Diastolic MKE and TKE increased under stress, in particular at peak early filling and peak atrial contraction. Augmented LV inotropy and cardiac frequency also caused an increase in direct flow and retained inflow MKE and TKE. However, the TKE/KE ratio remained comparable between rest and stress conditions, suggesting that LV intracavitary fluid dynamics can adapt to stress conditions without altering the TKE to KE balance of the normal left ventricle at rest.

Comparing left ventricular assist device inflow cannula angle between median sternotomy and thoracotomy using 3D reconstructions

BACKGROUND

Left ventricular assist device (LVAD) implantation via thoracotomy has many potential advantages compared to conventional sternotomy, including improved inflow cannula (IFC) positioning. We compared the difference in IFC angles, postoperative, and long-term outcomes for patients with LVADs implanted via thoracotomy and sternotomy.

METHODS

A single-center, retrospective analysis of 14 patients who underwent thoracotomy implantation was performed and matched with 28 patients who underwent sternotomy LVAD implantations for a total of 42 patients. Inclusion required a minimum LVAD support duration of 30 days and excluded concomitant procedures. A postoperative CT-chest was used to measure the angle the between the IFC and mitral valve in two-dimensions and results were compared with three-dimensional reconstruction using the same CT chest. Outcome data were extracted from medical records.

RESULTS

There was no significant difference in gender, INTERMACS score, BMI, or age between the two groups. Median cardiopulmonary bypass time was longer in the thoracotomy group compared to the sternotomy group, 107 min (86-122) versus 76 min (56-93), p < 0.01. 3D reconstructions revealed less deviation of the IFC away from the mitral valve in devices implanted via thoracotomy compared to sternotomy, median (IQR) angle 16.3° (13.9°-21.0°) versus 23.2° (17.9°-26.4°), p < 0.01. Rates of pump thrombosis, stroke, and gastrointestinal bleeding were not significantly different.

CONCLUSIONS

Devices implanted via thoracotomy demonstrated less deviation away from mitral valve. However, there was no difference in morbidity between the two approaches. 3D reconstruction of the heart is an innovative technique to measure angulation and is clinically advantageous when compared to 2D imaging.

Inflow cannula position as risk factor for stroke in patients with HeartMate 3 left ventricular assist devices

Background

A relation between the left ventricular assist device inflow cannula (IC) malposition and pump thrombus has been reported. This study aimed to investigate if the pump position, derived from chest X‐rays in HeartMate 3 (HM3) patients, correlates with neurological dysfunction (ND), ischemic stroke (IS), hemorrhagic stroke (HS) and survival.

Methods

This analysis was performed on routinely acquired X‐rays of 42 patients implanted with a HM3 between 2014 and 2017. Device position was quantified in patients with and without ND from frontal and lateral X‐rays characterizing the IC and pump in relation to spine, diaphragm or horizontal line. The primary end‐point was freedom from stroke and survival one‐year after HM3 implantation stratified by pump position.

Results

The analysis of X‐rays, 33.5 (41.0) days postoperative, revealed a significant smaller IC angle of HM3 patients with ND versus no ND (0.1° ± 14.0° vs. 12.9° ± 10.1°, p = 0.005). Additionally, the IC angle in the frontal view, IS: 4.1 (20.9)° versus no IS: 13.8 (7.5)°, p = 0.004 was significantly smaller for HM3 patients with IS. Using receiver operating characteristics derived cut‐off, IC angle <10° provided 75% sensitivity and 100% specificity (C‐statistic = 0.85) for predicting IS. Stratified by IC angle, freedom from IS at 12 months was 100% (>10°) and 60% (<10°) respectively (p = 0.002). No significant differences were found in any end‐point between patients with and without HS. One‐year survival was significantly higher in patients with IC angle >10° versus <10° (100% vs. 71.8%, p = 0.012).

Conclusions

IC malposition derived from standard chest X‐rays serves as a risk factor for ND, IS and worse survival in HM3 patients.

Effects of the atrium on intraventricular flow patterns during mechanical circulatory support

Simulations of the ventricular flow patterns during left ventricular assist device (LVAD) support are mainly performed with idealized cylindrical inflow, neglecting the influence of the atrial vortex. In this study, the influence of the left atrium (LA) on the intra-ventricular flow was investigated via Computational Fluid Dynamics (CFD) simulations. Ventricular flow was simulated by a combined Eulerian (carrier flow)/Lagrangian (particles) approach taking into account either the LA or a cylindrical inflow section to mimic a fully support condition. The flow deviation at the mitral valve, the blood low-velocity volume as well as the residence time and shear stress history of the particles were calculated. Inclusion of the LA deflects the flow at the mitral valve by 25°, resulting in an asymmetric flow jet entering the left ventricle. This reduced the ventricular low-velocity volume by 40% (from 6.4 to 3.9 cm³), increased (40%) the shear stress experienced by particles and correspondingly increased (27%) their residence time. Under the studied conditions, the atrial geometry plays a major role in the development of intraventricular flow patterns. A reliable prediction of blood flow dynamics and consequently thrombosis risk analysis within the ventricle requires the consideration of the LA in computational simulations.

Effect of Timings of the Lavare Cycle on the Ventricular Washout in an In Vitro Flow Visualization Setup

Left ventricular assist devices inherently alter the intraventricular flow field and create areas of blood stasis with potential thrombus formation. The Lavare cycle of the Medtronic HeartWare HVAD was designed to improve ventricular washout. This study aims to evaluate its effects on ventricular washout in a pulsatile in vitro setting with a focus on the timing of pump speed changes. Ventricular flow fields were obtained via particle image velocimetry in two modes: With constant left ventricular assist devices speed and with the Lavare cycle applied. The start of the Lavare cycle was shifted over an entire cardiac cycle, and ventricular washout was evaluated based on velocity fields, kinetic energy, and normalized pulsatility of flow fields. The ventricular flow fields showed dependence on the timing of the Lavare cycle and interaction between speed changes and the cardiac phase. Higher apical velocity was observed for speed decreases at the late E wave and for increases at mid systole by 29% (P = 0.002) and 61% (P < 0.001), respectively. Mean apical kinetic energy for these phases also increased by 21% (P = 0.0013) and 46% (P < 0.001). The Lavare cycle generally promotes higher apical washout and can specifically generate further improved washout if speed steps are applied at the correct timing on the cardiac cycle.

Pump position and thrombosis in ventricular assist devices: Correlation of radiographs and CT data

Malpositioning of left ventricular assist devices (LVAD) is a risk factor for thrombosis, but its identification from clinical imaging remains challenging. X-rays and CT scans were analyzed and parameters identified that correlated to pump thrombosis. Retrospective imaging data of patients (n = 115) with HeartmateII (HMII) or HVAD were analyzed in two groups (pump-thrombosis PT, n = 15 vs matched control group NT, n = 15) using routine X-rays and CT scans. In CT, directional deviations of the inflow cannula in three-chamber and two-chamber view (α and β angles) were identified. In HVAD PT frontal radiographs showed reduced pump body area and smaller minor axis (PT 41.3 ± 4.8 mm vs NT 34.9 ± 6.0 mm, p = 0.026), and in the lateral radiographs the visibility of the inflow cannula served as a predictive parameter for PT. In HMII patients, no parameters were associated with PT. The angle α differed significantly (NT −1.2 ± 7.5°, PT −22.0 ± 4.7°, p = 0.006) in HVAD patients. Further, correlations of x-ray parameters with CT angles α and β showed that radiographs can be used to identify malpositioned pumps. Well-aligned inflow cannula positions are essential. HVAD patients with a posterior rotation of the inflow cannula have a higher risk of pump thrombosis. This risk can reliably be identified from routine radiographs.

Intraventricular flow features and cardiac mechano-energetics after mitral valve interventions – feasibility of an isolated heart model

The aim of this work was the development of an isolated heart setup to delineate the interactions between intraventricular flow features, hemodynamic parameters and mechano-energetics after certain mitral valve therapies. Five porcine hearts were explanted and prepared for (i) edge-to-edge mitral valve repair, (ii) implantation of a rotatable biscupid mechanical valve prosthesis. Flow structures were visualized using echocardiography while hemodynamics was recorded in terms of pressures, flow rates and ventricular volume. Hemodynamic and cardiac mechano-energetics implied a marginal effect (<5%) of alternating leaflet orientation on ventricular pre-load and stroke work. After edge-to-edge repair, substantial variations in flow structures were observed. Beside promoting profound insights into fundamental physiologic mechanisms, the setup may be used for validation of computer aided therapy planning tools.

The influence of left ventricular assist device inflow cannula position on thrombosis risk

The use of left ventricular assist devices (LVADs) as a treatment method for heart failure patients has been steadily increasing; however, pathological studies showed presence of thrombi around the HeartWare ventricular assist device inflow cannula (IC) in more than 95% of patients after device explantation. Flow fields around the IC might trigger thrombus formation and require further investigation. In this study flow dynamics parameters were evaluated for different patient geometries using computational fluid dynamics (CFD) simulations. Left ventricular (LV) models of two LVAD patients were obtained from CT scans. The LV volumes of Patient 1 (P1) and Patient 2 (P2) were 264 and 114 cm3 with an IC angle of 20° and 9° from the mitral-IC tip axis at the coronal plane. The IC insertion site at the apex was central for P1, whereas it was lateral for P2. Transient CFD simulations were performed over 9 cardiac cycles. The wedge area was defined from the cannula tip to the wall of the LV apex. Mean velocity magnitude and blood stagnation region (volume with mean velocity <5 mm/s) as well as the wall shear stress (WSS) at the IC surface were calculated. Cardiac support resulted in a flow mainly crossing the ventricle from the mitral valve to the LVAD cannula for P2, while the main inflow jet deviated toward the septal wall in P1. Lower WSS at the IC surface and consequently larger stagnation volumes were observed for P2 (P1: 0.17, P2: 0.77 cm3 ). Flow fields around an LVAD cannula can be influenced by many parameters such as LV size, IC angle, and implantation site. Careful consideration of influencing parameters is essential to get reliable evaluations of the apical flow field and its connection to apical thrombus formation. Higher blood washout and lower stagnation were observed for a central implantation of the IC at the apex.

Validation of numerically simulated ventricular flow patterns during left ventricular assist device support

Intraventricular flow patterns during left ventricular assist device support have been investigated via computational fluid dynamics by several groups. Based on such simulations, specific parameters for thrombus formation risk analysis have been developed. However, computational fluid dynamic simulations of complex flow configurations require proper validation by experiments. To meet this need, a ventricular model with a well-defined inflow section was analyzed by particle image velocimetry and replicated by transient computational fluid dynamic simulations. To cover the laminar, transitional, and turbulent flow regime, four numerical methods including the laminar, standard k-omega, shear-stress transport, and renormalized group k-epsilon were applied and compared to the particle image velocimetry results in 46 different planes in the whole left ventricle. The simulated flow patterns for all methods, except renormalized group k-epsilon, were comparable to the flow patterns measured using particle image velocimetry (absolute error over whole left ventricle: laminar: 10.5, standard k-omega: 11.3, shear–stress transport: 11.3, and renormalized group k-epsilon: 17.8 mm/s). Intraventricular flow fields were simulated using four numerical methods and validated with experimental particle image velocimetry results. In the given setting and for the chosen boundary conditions, the laminar, standard K-omega, and shear–stress transport methods showed acceptable similarity to experimental particle image velocimetry data, with the laminar model showing the best transient behavior.