The Influence of Rotary Blood Pump Speed Modulation on the Risk of Intraventricular Thrombosis

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.

Study on the hemodynamic effects of different pulsatile working modes of a rotary blood pump using a microfluidic platform that realizes in vitro cell culture effectively

Rotary blood pumps (RBPs) operating at a constant speed generate non-physiologic blood pressure and flow rate, which can cause endothelial dysfunction, leading to adverse clinical events in peripheral blood vessels and other organs. Notably, pulsatile working modes of the RBP can increase vascular pulsatility to improve arterial endothelial function. However, the laws and related mechanisms of differentially regulating arterial endothelial function under different pulsatile working modes are still unclear. This knowledge gap hinders the optimal selection of the RBP working modes. To address these issues, this study developed a multi-element in vitro endothelial cell culture system (ECCS), which could realize in vitro cell culture effectively and accurately reproduce blood pressure, shear stress, and circumferential strain in the arterial endothelial microenvironment. Performance of this proposed ECCS was validated with numerical simulation and flow experiments. Subsequently, this study investigated the effects of four different pulsation frequency modes that change once every 1-4-fold cardiac cycles (80, 40, 80/3, and 20 cycles per min, respectively) of the RBP on the expression of nitric oxide (NO) and reactive oxygen species (ROS) in endothelial cells. Results indicated that the 2-fold and 3-fold cardiac cycles significantly increased the production of NO and prevented the excessive generation of ROS, potentially minimizing the occurrence of endothelial dysfunction and related adverse events during the RBP support, and were consistent with animal study findings. In general, this study may provide a scientific basis for the optimal selection of the RBP working modes and potential treatment options for heart failure.

Numerical simulation analysis of multi-scale computational fluid dynamics on hemodynamic parameters modulated by pulsatile working modes for the centrifugal and axial left ventricular assist devices

Continuous flow (CF) left ventricular assist devices (LVAD) operate at a constant speed mode, which could result in increased risk of adverse events due to reduced vascular pulsatility. Consequently, pump speed modulation algorithms have been proposed to augment vascular pulsatility. However, the quantitative local hemodynamic effects on the aorta when the pump is operating with speed modulation using different types of CF-LVADs are still under investigation. The computational fluid dynamics (CFD) study was conducted to quantitatively elucidate the hemodynamic effects on a clinical patient-specific aortic model under different speed patterns of CF-LVADs. Pressure distribution, wall shear stress (WSS), time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT), and velocity were calculated to compare their differences at constant and pulsatile speeds under centrifugal and axial LVAD support. Results showed that pulse pressure on the aorta was significantly larger under pulsatile speed mode than that under constant speed mode for both CF-LVADs, indicating enhanced aorta pulsatility, as well as the higher peak blood flow velocity on some representative slices of aorta. Pulsatile speed modulation enhanced peak WSS compared to constant speed; high TAWSS region appeared near the branch of left common carotid artery and distal aorta regardless of speed modes and CF-LVADs but these regions also had low OSI; RRT was almost the same for all the cases. This study may provide a basis for the scientific and reasonable selection of the pulsatile speed patterns of CF-LVADs for treating heart failure patients.

In Vitro Investigation of the Effect of the Timing of Left Ventricular Assist Device Speed Modulation on Intraventricular Flow Patterns

Thromboembolic events remain a common complication for left ventricular assist device (LVAD) patients. To prevent in-pump thrombosis, third-generation LVADs use speed modulation, which is not synchronized with the native left ventricle (LV) contractility. This study aims to investigate the effect of speed modulation on intraventricular flow patterns, and specifically, the impact of timing relative to pressure variations in the LV. Stereo-particle image velocimetry measurements were performed in a patient-derived LV implanted with an LVAD, for different timings of the speed modulation and speed. Speed modulation has a strong effect on instantaneous afterload and flowrate (-16% and +20%). The different timings of the speed modulation resulted in different flowrate waveforms, exhibiting different maxima (5.3-5.9 L/min, at constant average flowrate). Moreover, the timing of the speed modulation was found to strongly influence intraventricular flow patterns, specifically, stagnation areas within the LV. These experiments highlight, once more, the complex relationship between LVAD speed, hemodynamic resistance, and intraventricular pressure. Overall, this study demonstrates the importance of considering native LV contractility in future LVAD controls, to improve hemocompatibility and reduce the risk of thromboembolic complications.

A Mathematical Model of Artificial Pulse Synchronization for the HeartMate3 Left Ventricular Assist Device

Constant speed control of rotary LVADs attenuates vascular pulsatility, which has been linked to clinical complications such as thrombus formation, bleeding, and valvular dysfunction. Speed modulation can improve pulsatility and washout, but optimization requires coordination with the native heartbeat. A simple mathematical model of the left ventricle-left ventricular assist device (LV-LVAD) flow interaction was developed that sums the individual contributions of the native LV and the HeartMate3 artificial pulse (AP) to predict the total systemic flow. The model flow and pulsatility predictions results were in good agreement with experimental data from a mock circulatory loop measured for full bypass support conditions. The model was used to evaluate three schemes for optimizing the synchronization of the AP with the native heart. The optimized interaction occurred when the AP speed increase occurred during contraction, resulting in a doubling of flow pulsatility, and corresponded to an increase in the area enclosed by the dynamic pressure-flow relation. The model provides a simple tool for exploring the optimization of LVAD speed modulation that can reduce the time and expense of mock loop studies during the development process.

Left atrial appendage occlusion in ventricular assist device patients to decrease thromboembolic events: A computer simulation study

Atrial fibrillation (AF) is a common comorbidity in left ventricular assist device (LVAD) patients and has been identified as a risk factor for thromboembolic stroke. Blood stagnation within the left atrial appendage (LAA) is considered a possible major source of thrombosis and clinical studies have shown reduced thromboembolic risk after LAA occlusion (LAAO). Therefore, this study aims to investigate the effect of LAAO on thrombosis-related parameters using patient-specific simulations. Left ventricular and left atrial geometries of an LVAD patient were obtained from computed tomography and combined with hemodynamic data with either sinus rhythm (SR) or AF generated by a lumped parameter model. In four simulations applying contractile walls, stagnation volume and blood residence times were evaluated with or without AF and with or without LAAO. Reduced atrial contraction in AF resulted in unfavorable flow dynamics within the left atrium. The average atrial velocity was lower for the AF simulation when compared to SR, resulting in a 55% increase in the atrial stagnation volume (from 4.2 to 6.5 cm³). Moreover, blood remained in the LAA for more than 8 cardiac cycles. After LAAO the atrial stagnation decreased from 4.2 to 1.4 cm³ for SR and from 6.5 to 2.3 cm³ for the AF simulation. A significant stagnation volume was found in the LAA for both SR and AF, with larger values occurring with AF. These regions are known as potential sources for thrombus formation and can be diminished by LAAO. This significantly improved the thrombus-related flow parameters and may also lower the risk of thromboembolic events from the appendage.

The hemodynamics and blood trauma in axial blood pump under different operating model

Speed modulation of blood pump has been proved to help restore vascular pulsatility and implemented clinically during treatment for cardiac failure. However, its effect on blood trauma has not been studied thoroughly. In this paper, we study the flow field of an axial pump FW-X under the modes of co-pulse, counter pulse and constant speed to evaluate the blood trauma. Based on the coupling model of cardiovascular system and axial blood pump, aortic pressure and the pump flow were obtained and applied as the boundary conditions at the pump outlet and inlet. The level of shear stress and hemolysis index were derived from computational fluid dynamics (CFD) simulation. Results showed the constant speed mode had the lowest shear stress level and hemolytic index at the expense of diminished pulsatility. Compared with the constant speed mode, the hemolysis index of co-pulse and counter pulse mode was higher, but it was helpful to restore vascular pulsatility. This method can be easily incorporated in the in vitro testing phase to analyze and decrease a pump's trauma before animal experimentation, thereby reducing the cost of blood pump development.

Fuzzy-based modeling and speed optimization of a centrifugal blood pump using a modified and constrained Bees algorithm

BACKGROUND AND OBJECTIVE

Side effects that may occur when using blood pumps for treatment of patients are the main limitations on pump rotational speed determination. Efforts are being made to reduce side effects in both design and usage procedures. In determining the pump speed for treatment, decreasing the pressure on the main artery and preserving the valve functions are taken into consideration. In addition to these, the parameters considered for design which include pump efficiency and mechanical effects on blood cells, should also be taken into consideration. In this study, the aim is to obtain the optimum pump speed for the maximum hydraulic efficiency and minimum wall shear stresses that occur inside the pump.

METHODS

Blood pump modeling based on fuzzy logic is created on the hydraulic performance data of a centrifugal blood pump, whose design, CFD analysis, manufacture and experimental testing have been performed previously. Using this fuzzy logic model, the optimum pump speeds were determined using the Bees Algorithm, an intuitive optimization algorithm, in the operating range 1-7 L/min fluid flow rate. In the optimization process, the aim is to achieve minimum shear stress with maximal efficiency. Intravascular pressure limits (90-160 mm-Hg) were set as pressure constraints.

RESULTS

The optimum operating point is obtained as a 3350 rpm pump speed and a 4.35 L/min flow rate. At this operating point, CFD simulation is performed, and maximum wall shear stress was found to be 1458 Pa and its efficiency as 34.2%.

CONCLUSIONS

In addition to the parameters commonly used in the pump speed optimization of blood pumps, the use of wall shear stresses and pump efficiency can provide certain improvements.

A New Mathematical Numerical Model to Evaluate the Risk of Thrombosis in Three Clinical Ventricular Assist Devices

(1) Background: Thrombosis is the main complication in patients supported with ventricular assist devices (VAD). Models that accurately predict the risk of thrombus formation in VADs are still lacking. When VADs are clinically assisted, their complex geometric configuration and high rotating speed inevitably generate complex flow fields and high shear stress. These non-physiological factors can damage blood cells and proteins, release coagulant factors and trigger thrombosis. In this study, a more accurate model for thrombus assessment was constructed by integrating parameters such as shear stress, residence time and coagulant factors, so as to accurately assess the probability of thrombosis in three clinical VADs. (2) Methods: A mathematical model was constructed to assess platelet activation and thrombosis within VADs. By solving the transport equation, the influence of various factors such as shear stress, residence time and coagulation factors on platelet activation was considered. The diffusion equation was applied to determine the role of activated platelets and substance deposition on thrombus formation. The momentum equation was introduced to describe the obstruction to blood flow when thrombus is formed, and finally a more comprehensive and accurate model for thrombus assessment in patients with VAD was obtained. Numerical simulations of three clinically VADs (CH-VAD, HVAD and HMII) were performed using this model. The simulation results were compared with experimental data on platelet activation caused by the three VADs. The simulated thrombogenic potential in different regions of MHII was compared with the frequency of thrombosis occurring in the regions in clinic. The regions of high thrombotic risk for HVAD and HMII observed in experiments were compared with the regions predicted by simulation. (3) Results: It was found that the percentage of activated platelets within the VAD obtained by solving the thrombosis model developed in this study was in high agreement with the experimental data (r² = 0.984), the likelihood of thrombosis in the regions of the simulation showed excellent correlation with the clinical statistics (r² = 0.994), and the regions of high thrombotic risk predicted by the simulation were consistent with the experimental results. Further study revealed that the three clinical VADs (CH-VAD, HVAD and HMII) were prone to thrombus formation in the inner side of the secondary flow passage, the clearance between cone and impeller, and the corner region of the inlet pipe, respectively. The risk of platelet activation and thrombus formation for the three VADs was low to high for CH-VAD, HVAD, and HM II, respectively. (4) Conclusions: In this study, a more comprehensive and accurate thrombosis model was constructed by combining parameters such as shear stress, residence time, and coagulation factors. Simulation results of thrombotic risk received with this model showed excellent correlation with experimental and clinical data. It is important for determining the degree of platelet activation in VAD and identifying regions prone to thrombus formation, as well as guiding the optimal design of VAD and clinical treatment.

An Accelerated Thrombosis Model for Computational Fluid Dynamics Simulations in Rotary Blood Pumps

Purpose

Thrombosis ranks among the major complications in blood-carrying medical devices and a better understanding to influence the design related contribution to thrombosis is desirable. Over the past years many computational models of thrombosis have been developed. However, numerically cheap models able to predict localized thrombus risk in complex geometries are still lacking. The aim of the study was to develop and test a computationally efficient model for thrombus risk prediction in rotary blood pumps.

Methods

We used a two-stage approach to calculate thrombus risk. The first stage involves the computation of velocity and pressure fields by computational fluid dynamic simulations. At the second stage, platelet activation by mechanical and chemical stimuli was determined through species transport with an Eulerian approach. The model was compared with existing clinical data on thrombus deposition within the HeartMate II. Furthermore, an operating point and model parameter sensitivity analysis was performed.

Results

Our model shows good correlation (R² > 0.93) with clinical data and identifies the bearing and outlet stator region of the HeartMate II as the location most prone to thrombus formation. The calculation of thrombus risk requires an additional 10–20 core hours of computation time.

Conclusion

The concentration of activated platelets can be used as a surrogate and computationally low-cost marker to determine potential risk regions of thrombus deposition in a blood pump. Relative comparisons of thrombus risk are possible even considering the intrinsic uncertainty in model parameters and operating conditions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13239-021-00606-y.

Mathematical and Computational Modeling of Device-Induced Thrombosis

Given the extensive and routine use of cardiovascular devices, a major limiting factor to their success is the thrombotic rate that occurs. This both poses direct risk to the patient and requires counterbalancing with anticoagulation and other treatment strategies, contributing additional risks. Developing a better understanding of the mechanisms of device-induced thrombosis to aid in device design and medical management of patients is critical to advance the ubiquitous use and durability. Thus, mathematical and computational modelling of device-induced thrombosis has received significant attention recently, but challenges remain. Additional areas that need to be explored include microscopic/macroscopic approaches, reconciling physical and numerical timescales, immune/inflammatory responses, experimental validation, and incorporating pathologies and blood conditions. Addressing these areas will provide engineers and clinicians the tools to provide safe and effective cardiovascular devices.

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.

Understanding the influence of left ventricular assist device inflow cannula alignment and the risk of intraventricular thrombosis

BACKGROUND

Adverse neurological events associated with left ventricular assist devices (LVADs) have been suspected to be related to thrombosis. This study aimed to understand the risks of thrombosis with variations in the implanted device orientation. A severely dilated pulsatile patient-specific left ventricle, modelled with computational fluid dynamics, was utilised to identify the risk of thrombosis for five cannulation angles. With respect to the inflow cannula axis directed towards the mitral valve, the other angles were 25° and 20° towards the septum and 20° and 30° towards the free wall.

RESULTS

Inflow cannula angulation towards the free wall resulted in longer blood residence time within the ventricle, slower ventricular washout and reduced pulsatility indices along the septal wall. Based on the model, the ideal inflow cannula alignment to reduce the risk of thrombosis was angulation towards the mitral valve and up to parallel to the septum, avoiding the premature clearance of incoming blood.

CONCLUSIONS

This study indicates the potential effects of inflow cannulation angles and may guide optimised implantation configurations; however, the ideal approach will be influenced by other patient factors and is suspected to change over the course of support.

Tipless transseptal cannula concept combines improved hemodynamic properties and risk-reduced placement: An in silico proof-of-concept

As a leading cause of death worldwide, heart failure is a serious medical condition in which many critically ill patients require temporary mechanical circulatory support (MCS) as a bridge-to-recovery or bridge-to-decision. In many cases, the TandemHeart system is used to unload the left heart by draining blood from the left atrium (LA) to the femoral artery via a transseptal multistage cannula. However, even though the correct positioning of the cannula is crucial for a safe treatment, the long cannula tip currently used in transseptal cannulas complicates positioning, making the cannula vulnerable to displacement during MCS. To overcome these limitations, we propose the development of a new tipless transseptal cannula with improved hemodynamic properties. We discuss the tipless cannula concept by comparing it to the common multistage cannula concept using computational fluid dynamics simulations and assess the flow field in the LA, the wall shear stresses (WSS), and the pressure loss. Across the two distinct time points of end-systole and end-diastole and two drainage flow rates of 3.5 and 5.0 L/min, we find a more homogeneous inlet flow pattern for the tipless cannula concept, accompanied by a remarkably reduced area of platelet-activating WSS (up to 10-times smaller area compared to the multistage cannula). Moreover, pressure loss is up to 14.5% lower in the tipless cannula concept, confirming overall improved hemodynamic properties of the tipless cannula concept. Finally, a diameter-dependent study reveals that lower WSS and pressure losses can be further reduced by large-lumen designs for any simulation setting. Overall, our results suggest that a tipless cannula concept remedies the crucial disadvantages of a long-tip multistage cannula by reducing the risk of misplacement, and it furthermore promotes optimized hemodynamics. With this successful proof-of-concept, we underscore the potential for and encourage the realization of further experimental investigations regarding the development of a tipless transseptal cannula for MCS.

Experimental investigation of the influence of the hydraulic performance of an axial blood pump on intraventricular blood flow

Blood flow inside the left ventricle (LV) is a concern for blood pump use and contributes to ventricle suction and thromboembolic events. However, few studies have examined blood flow inside the LV after a blood pump was implanted. In this study, in vitro experiments were conducted to emulate the intraventricular blood flow, such as blood flow velocity, the distribution of streamlines, vorticity and the standard deviation of velocity inside the LV during axial blood pump support. A silicone LV reconstructed from computerized tomography (CT) data of a heart failure patient was incorporated into a mock circulatory loop (MCL) to simulate human systemic circulation. Then, the blood flow inside the ventricle was examined by particle image velocimetry (PIV) equipment. The results showed that the operating conditions of the axial blood pump influenced flow patterns within the LV and areas of potential blood stasis, and the intraventricular swirling flow was altered with blood pump support. The presence of vorticity in the LV from the thoracic aorta to the heart apex can provide thorough washing of the LV cavity. The gradually extending stasis region in the central LV with increasing blood pump support is necessary to reduce the thrombosis potential in the LV.

Evaluation of cardiac beat synchronization control for a rotary blood pump on valvular regurgitation with a mathematical model

We have studied the cardiac beat synchronization (CBS) control for a rotary blood pump (RBP) and revealed that it can promote pulsatility and reduce cardiac load. Besides, patients with LVAD support sometimes suffer from aortic and mitral regurgitation (AR and MR). A control method for the RBP should be validated in wider range of conditions to clarify its benefits and pitfalls prior to clinical application. In this study, we evaluated pulsatility and cardiac load reduction obtained with the CBS control on valvular failure conditions with a mathematical model. Diastolic assist could reduce cardiac load on the left ventricle by decreasing external work of the ventricle even in MR cases while it was not so effective in AR cases. Systolic assist can still promote pulsatility in AR and MR cases; however, aortic valve function should be carefully confirmed since pulse pressure can be wider not due to systolic assist but to AR.

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

Investigations of ventricular flow patterns during mechanical circulatory support are limited to in vitro flow models or in silico simulations, which cannot fully replicate the complex anatomy and contraction of the heart. Therefore, the feasibility of using echocardiographic particle image velocimetry (Echo-PIV) was evaluated in an isolated working heart setup. Porcine hearts were connected to an isolated, working heart setup and a left ventricular assist device (LVAD) was implanted. During different levels of LVAD support (unsupported, partial support, full support), microbubbles were injected and echocardiographic images were acquired. Iterative PIV algorithms were applied to calculate flow fields. The isolated heart setup allowed different hemodynamic situations. In the unsupported heart, diastolic intra-ventricular blood flow was redirected at the heart’s apex towards the left ventricular outflow tract (LVOT). With increasing pump speed, large vortex formation was suppressed, and blood flow from the mitral valve directly entered the pump cannula. The maximum velocities in the LVOT were significantly reduced with increasing support. For the first time, cardiac blood flow patterns during LVAD support were visualized and quantified in an ex vivo model using Echo-PIV. The results reveal potential regions of stagnation in the LVOT and, in future the methods might be also used in clinical routine to evaluate intraventricular flow fields during LVAD support.

Artificial Organs 2018: A Year in Review

In this Editor's Review, articles published in 2018 are organized by category and summarized. We provide a brief reflection of the research and progress in artificial organs intended to advance and better human life while providing insight for continued application of these technologies and methods. Artificial Organs continues in the original mission of its founders "to foster communications in the field of artificial organs on an international level." Artificial Organs continues to publish developments and clinical applications of artificial organ technologies in this broad and expanding field of organ Replacement, Recovery, and Regeneration from all over the world. Peer-reviewed special issues this year included contributions from the 13th International Conference on Pediatric Mechanical Circulatory Support Systems and Pediatric Cardiopulmonary Perfusion edited by Dr. Akif Undar, and the 25th Congress of the International Society for Mechanical Circulatory Support edited by Dr. Marvin Slepian. Additionally, many editorials highlighted the worldwide survival differences in hemodialysis and perspectives on mechanical circulatory support and stem cell therapies for cardiac support. We take this time also to express our gratitude to our authors for offering their work to this journal. We offer our very special thanks to our reviewers who give so generously of time and expertise to review, critique, and especially provide meaningful suggestions to the author's work whether eventually accepted or rejected. Without these excellent and dedicated reviewers the quality expected from such a journal could not be possible. We also express our special thanks to our Publisher, John Wiley & Sons for their expert attention and support in the production and marketing of Artificial Organs. We look forward to reporting further advances in the coming years.