The use of computational fluid dynamics in the development of ventricular assist devices

Analysis of non-physiological shear stress-induced red blood cell trauma across different clinical support conditions of the blood pump

Systematic research into device-induced red blood cell (RBC) damage beyond hemolysis, including correlations between hemolysis and RBC-derived extracellular vesicles, remains limited. This study investigated non-physiological shear stress-induced RBC damage and changes in related biochemical indicators under two blood pump clinical support conditions. Pressure heads of 100 and 350 mmHg, numerical simulation methods, and two in vitro loops were utilized to analyze the shear stress and changes in RBC morphology, hemolysis, biochemistry, metabolism, and oxidative stress. The blood pump created higher shear stress in the 350-mmHg condition than in the 100-mmHg condition. With prolonged blood pump operation, plasma-free hemoglobin and cholesterol increased, whereas plasma glucose and nitric oxide decreased in both loops. Notably, plasma iron and triglyceride concentrations increased only in the 350-mmHg condition. The RBC count and morphology, plasma lactic dehydrogenase, and oxidative stress across loops did not differ significantly. Plasma extracellular vesicles, including RBC-derived microparticles, increased significantly at 600 min in both loops. Hemolysis correlated with plasma triglyceride, cholesterol, glucose, and nitric oxide levels. Shear stress, but not oxidative stress, was the main cause of RBC damage. Hemolysis alone inadequately reflects overall blood pump-induced RBC damage, suggesting the need for additional biomarkers for comprehensive assessments.

Numerical methods for hemolysis and thrombus evaluation in the percutaneous ventricular assist device

BACKGROUND

A percutaneous ventricular assist device (pVAD) is an effective method to treat heart failure, but its complications, mainly hemolysis and thrombus formation, cannot be ignored. Accurate evaluation of hemolysis and thrombus formation in pVAD is essential to guide the development of pVAD and reduce the incidence of complications.

METHODS

This study optimized the numerical model to predict hemolysis and thrombus formation in pVAD. The hemolysis model is based on the power law function, and the multi-component thrombus prediction model is improved by introducing the von Willebrand factor.

RESULTS

The error between the numerical simulation and the hydraulic performance experiment is within 5%. The numerical results of hemolysis are in good agreement with those of in vitro experiments. Meanwhile, the thrombus location predicted by the numerical model is the same as that found in the in vivo experiment.

CONCLUSION

The numerical model suggested in this study may therefore accurately assess the possible hemolytic and thrombotic dangers in pVAD, making it an effective tool to support the development of pVAD.

Enhancing the implantation of mechanical circulatory support devices using computational simulations

Introduction: Patients with end-stage heart failure (HF) may need mechanical circulatory support such as a left ventricular assist device (LVAD). However, there are a range of complications associated with LVAD including aortic regurgitation (AR) and thrombus formation. This study assesses whether the risk of developing aortic conditions can be minimised by optimising LVAD implantation technique. Methods: In this work, we evaluate the aortic flow patterns produced under different geometrical parameters for the anastomosis of the outflow graft (OG) to the aorta using computational fluid dynamics (CFD). A three-dimensional aortic model is created and the HeartMate III OG positioning is simulated by modifying (i) the distance from the anatomic ventriculo-arterial junction (AVJ) to the OG, (ii) the cardinal position around the aorta, and (iii) the angle between the aorta and the OG. The continuous LVAD flow and the remnant native cardiac cycle are used as inlet boundaries and the three-element Windkessel model is applied at the pressure outlets. Results: The analysis quantifies the impact of OG positioning on different haemodynamic parameters, including velocity, wall shear stress (WSS), pressure, vorticity and turbulent kinetic energy (TKE). We find that WSS on the aortic root (AoR) is around two times lower when the OG is attached to the coronal side of the aorta using an angle of 45° ± 10° at a distance of 55 mm. Discussion: The results show that the OG placement may significantly influence the haemodynamic patterns, demonstrating the potential application of CFD for optimising OG positioning to minimise the risk of cardiovascular complications after LVAD implantation.

A homogenized two-phase computational framework for meso- and macroscale blood flow simulations

BACKGROUND AND OBJECTIVE

Owing to the complexity of physics linked with blood flow and its associated phenomena, appropriate modeling of the multi-constituent rheology of blood is of primary importance. To this effect, various kinds of computational fluid dynamic models have been developed, each with merits and limitations. However, when additional physics like thrombosis and embolization is included within the framework of these models, computationally efficient scalable translation becomes very difficult. Therefore, this paper presents a homogenized two-phase blood flow framework with similar characteristics to a single fluid model but retains the flow resolution of a classical two-fluid model. The presented framework is validated against four different sets of experiments.

METHODS

The two-phase model of blood presented here is based on the classical diffusion-flux framework. Diffusion flux models are known to be less computationally expensive than two-fluid multiphase models since the numerical implementation resembles single-phase flow models. Diffusion flux models typically use empirical slip velocity correlations to resolve the motion between phases. However, such correlations do not exist for blood. Therefore, a modified slip velocity equation is proposed, derived rigorously from the two-fluid governing equations. An additional drag law for red blood cells (RBCs) as a function of volume fraction is evaluated using a previously published cell-resolved solver. A new hematocrit-dependent expression for lift force on RBCs is proposed. The final governing equations are discretized and solved using the open-source software OpenFOAM.

RESULTS

The framework is validated against four sets of experiments: (i) flow through a rectangular microchannel to validate RBC velocity profiles against experimental measurements and compare computed hematocrit distributions against previously reported simulation results (ii) flow through a sudden expansion microchannel for comparing experimentally obtained contours of hematocrit distributions and normalized cell-free region length obtained at different flowrates and inlet hematocrits, (iii) flow through two hyperbolic channels to evaluate model predictions of cell-free layer thickness, and (iv) flow through a microchannel that mimics crevices of a left ventricular assist device to predict hematocrit distributions observed experimentally. The simulation results exhibit good agreement with the results of all four experiments.

CONCLUSION

The computational framework presented in this paper has the advantage of resolving the multiscale physics of blood flow while still leveraging numerical techniques used for solving single-phase flows. Therefore, it becomes an excellent candidate for addressing more complicated problems related to blood flow, such as modeling mechanical entrapment of RBCs within blood clots, predicting thrombus composition, and visualizing clot embolization.

Multi-Method Investigation of Blood Damage Induced By Blood Pumps in Different Clinical Support Modes

To investigate the effects of blood pumps operated in different modes on nonphysiologic flow patterns, cell and protein function, and the risk of bleeding, thrombosis, and hemolysis, an extracorporeal blood pump (CentriMag) was operated in three clinical modalities including heart failure (HF), venous-venous (V-V) extracorporeal membrane oxygenation (ECMO), and venous-arterial (V-A) ECMO. Computational fluid dynamics (CFD) methods and coupled hemolysis models as well as recently developed bleeding and thrombosis models associated with changes in platelet and von Willebrand factor (vWF) function were used to predict hydraulic performance and hemocompatibility. The V-A ECMO mode had the highest flow losses and shear stress levels, the V-V ECMO mode was intermediate, and the HF mode was the lowest. Different nonphysiologic flow patterns altered cell/protein morphology and function. The V-A ECMO mode resulted in the highest levels of platelet activation, receptor shedding, vWF unfolding, and high molecular weight multimers vWF (HMWM-vWF) degradation, leading to the lowest platelet adhesion and the highest vWF binding capacity, intermediate in the V-V ECMO mode, and opposite in the HF mode. The V-A ECMO mode resulted in the highest risk of bleeding, thrombosis, and hemolysis, with the V-V ECMO mode intermediate and the HF mode lowest. These findings are supported by published experimental or clinical statistics. Further studies found that secondary blood flow passages resulted in the highest risk of blood damage. Nonphysiologic blood flow patterns were strongly associated with cell and protein function changing, blood damage, and complications.

Optimization of a centrifugal blood pump designed using an industrial method through experimental and numerical study

With improved treatment of coronary artery disease, more patients are surviving until heart failure occurs. This leads to an increase in patients needing devices for struggling with heart failure. Ventricular assist devices are known as the mainstay of these devices. This study aimed to design a centrifugal pump as a ventricular assist device. In order to design the pump, firstly, the geometrical parameters of the pump, including the gap distance, blade height, and position of the outlet relative to the blade, were investigated. Finally, the selected configuration, which had all the appropriate characteristics, both hydraulically and physiologically, was used for the rest of the study. The study of the blade, as the main component in energy transfer to the blood, in a centrifugal pump, has been considered in the present study. In this regard, the point-to-point design method, which is used in industrial applications, was implemented. The designer chooses the relationship between the blade angles at each radius in the point-to-point method. The present study selected logarithmic and second-order relations for designing the blade's profile. In total, 58 blades were examined in this study, which differed regarding blade inlet and outlet angles and the relationship between angle and radial position. ANSYS CFX 17.0 software was utilized to simulate blades' performances, and a benchmark pump provided by the US Food and Drug Administration (FDA) was used to validate the numerical simulations. Then, the selected impeller from the numerical investigation was manufactured, and its performance was compared experimentally with the FDA benchmark pump. A hydraulic test rig was also developed for experimental studies. The results showed that among the blades designed in this study, the blade with an input angle of 45° and an output angle of 55°, which is designed to implement a logarithmic relationship, has the best performance. The selected impeller configuration can increase the total head (at least by 20%) at different flow rates compared to the FDA pump.

The Effect of Mechanical Circulatory Support on Blood Flow in the Ascending Aorta: A Combined Experimental and Computational Study

Percutaneous mechanical circulatory support (MCS) devices are designed for short-term treatment in cases of acute decompensated heart failure as a bridge to transplant or recovery. Some of the known complications of MCS treatments are related to their hemodynamics in the aorta. The current study investigates the effect of MCS on the aortic flow. The study uses combined experimental and numerical methods to delineate complex flow structures. Particle image velocimetry (PIV) is used to capture the vortical and turbulent flow characteristics in a glass model of the human aorta. Computational fluid dynamics (CFD) analyses are used to complete the 3D flow in the aorta. Three specific MCS configurations are examined: a suction pump with a counterclockwise (CCW) rotating impeller, a suction pump with a clockwise (CW) rotating impeller, and a discharge pump with a straight jet. These models were examined under varying flow rates (1-2.5 L/min). The results show that the pump configuration strongly influences the flow in the thoracic aorta. The rotating impeller of the suction pump induces a dominant swirling flow in the aorta. The swirling flow distributes the incoming jet and reduces the turbulent intensity near the aortic valve and in the aorta. In addition, at high flow rates, the local vortices formed near the pump are washed downstream toward the aortic arch. Specifically, an MCS device with a CCW rotating impeller induces a non-physiological CCW helical flow in the descending aorta (which is opposite to the natural helical flow), while CW swirl combines better with the natural helical flow.

Hemocompatibility and biophysical interface of left ventricular assist devices and total artificial hearts

ABSTRACT

Over the past 2 decades, there has been a significant increase in the utilization of long-term mechanical circulatory support (MCS) for the treatment of cardiac failure. Left ventricular assist devices (LVADs) and total artificial hearts (TAHs) have been developed in parallel to serve as bridge-to-transplant and destination therapy solutions. Despite the distinct hemodynamic characteristics introduced by LVADs and TAHs, a comparative evaluation of these devices regarding potential complications in supported patients, has not been undertaken. Such a study could provide valuable insights into the complications associated with these devices. Although MCS has shown substantial clinical benefits, significant complications related to hemocompatibility persist, including thrombosis, recurrent bleeding, and cerebrovascular accidents. This review focuses on the current understanding of hemostasis, specifically thrombotic and bleeding complications, and explores the influence of different shear stress regimens in long-term MCS. Furthermore, the role of endothelial cells in protecting against hemocompatibility-related complications of MCS is discussed. We also compared the diverse mechanisms contributing to the occurrence of hemocompatibility-related complications in currently used LVADs and TAHs. By applying the existing knowledge, we present, for the first time, a comprehensive comparison between long-term MCS options.

Investigating the Hydrodynamics of Intravenous Drug Infusions

One of the most common reported adverse events for intravenous (IV) infusions are infusion site reactions, ranging from redness and pain at the site of infusion to thrombophlebitis.  The connection between drug infusion and what drives these adverse events is not well understood. To aid in understanding these phenomena, it is crucial to accurately characterize the evolving hemodynamic environment of the infusion site when developing new intravenous formulations, as too rapid dilution may cause precipitation in the vein, while too little dilution might contribute to phlebitis. In this study, a Large-Eddy Simulation (LES) turbulence modeling inside a Computational Fluid Dynamics (CFD) framework has been used to simulate the flow and mixing characteristics of an infusion entering the bloodstream. This work represents the first such study reporting transient flow fields for intravenous infusions using LES CFD simulations with a realistic non-Newtonian blood model. The output of the CFD model closely resembled the flow and mixing patterns generated in benchtop tests for infusions into a blood analogue and water as the venous fluid across a wide range of flow rates.  These models were then investigated further to compare how changes to the fluid rheology model, needle orientation and needle position within the vein resulted in altered mixing regimes at different flow rates.

Design and optimisation of an Intra-Aortic Shrouded rotor axial pump

Undesirable side effects in patients with a LVAD (Left Ventricular Assist Device) pump fitted include blood damage, thrombosis, blood traumatisation, and End-Organ Disfunctions. These side effects have generally been attributed to the high wall shear stresses and the induced turbulent flow. In this study, we introduce a novel design to address these effects by lowering the rotational speed and providing an optimum flow path design to minimise blood damage. We present an initial scheme for a new Intra-Aortic Shrouded Rotary Axial Pump and develop a sequence of pump geometries, for which the Taguchi Design Optimisation Method has been applied. We apply CFD tools to simulate the pressure rise, pump performance, hydraulic efficiency, wall shear stress, exposure time and mass flow rate. A prototype pump has been tested in a mock cardiovascular circuit using a water-glycerol solution. The optimum design delivered the desired pressure/mass flow rate characteristics at a significantly low rpm (2900 rpm). As a result, the estimated blood damage index is low, matching the design requirements. The theoretical performance was matched by experimental results.

Effect of Hematocrit and Elevated Beat Rate on the 12cc Penn State Pediatric Ventricular Assist Device

Congenital heart disease affects approximately 40,000 infants annually in the United States with 25% requiring invasive treatment. Due to limited number of donor hearts and treatment options available for children, pediatric ventricular assist devices (PVADs) are used as a bridge to transplant. The 12cc pneumatic Penn State PVAD is optimized to prevent platelet adhesion and thrombus formation at patient nominal conditions; however, children demonstrate variable blood hematocrit and elevated heart rates. Therefore, with pediatric patients exhibiting greater variability, particle image velocimetry is used to evaluate the PVAD with three non-Newtonian hematocrit blood analogs (20%, 40%, and 60%) and at two beat rates (75 and 120 bpm) to understand the device's performance. The flow fields demonstrate a strong inlet jet that transitions to a solid body rotation during diastole. During systole, the rotation dissipates and reorganizes into an outlet jet. This flow field is consistent across all hematocrits and beat rates but at a higher velocity magnitude during 120 bpm. There are also minor differences in flow field timing and surface washing due to hematocrit. Therefore, despite patient differences in hematocrit or required pumping output, thorough surface washing can be achieved in the PVAD by altering operating conditions, thus reducing platelet adhesion potential.

Resistance valves in circulatory loops have a significant impact on in vitro evaluation of blood damage caused by blood pumps: a computational study

Background: Hemolysis and its complications are major concerns during the clinical application of blood pumps. In-vitro circulatory testing loops have been employed as the key procedure to evaluate the hemolytic and thrombogenic performance of blood pumps during the development phase and before preclinical in-vivo animal studies. Except for the blood damage induced by the pump under test, blood damage induced by loop components such as the resistance valve may affect the accuracy, reproducibility, and intercomparability of test results. Methods: This study quantitatively investigated the impact of the resistance valve on in vitro evaluation of blood damage caused by blood pumps under different operating points. A series of idealized tubing models under the resistance valve with different openings were created. Three pumps - the FDA benchmark pump, the HeartMate 3 LVAD, and the CH-VAD - were involved in hypothetical tests. Eight operating points were chosen to cover a relatively wide spectrum of testing scenarios. Computational fluid dynamics (CFD) simulations of the tubing and pump models were conducted at the same operating points. Results and Conclusion: Overall, hemolysis and platelet activation induced by a typical resistance valve are equivalent to 17%-45% and 14%-60%, respectively, of those induced by the pump itself. Both ratios varied greatly with flow rate, valve opening and pump models. Differences in blood damage levels between different blood pumps or working conditions can be attenuated by up to 45%. Thus, hemolysis and platelet activation induced by the resistance valve significantly affect the accuracy of in-vitro hemocompatibility evaluations of blood pumps. A more accurate and credible method for hemocompatibility evaluations of blood pumps will benefit from these findings.

The Relation Between Viscous Energy Dissipation and Pulsation for Aortic Hemodynamics Driven by a Left Ventricular Assist Device

Left ventricular assist device (LVAD) provides mechanical circulatory support for patients with advanced heart failure. Treatment using LVAD is commonly associated with complications such as stroke and gastro-intestinal bleeding. These complications are intimately related to the state of hemodynamics in the aorta, driven by a jet flow from the LVAD outflow graft that impinges into the aorta wall. Here we conduct a systematic analyses of hemodynamics driven by an LVAD with a specific focus on viscous energy transport and dissipation. We conduct a complementary set of analysis using idealized cylindrical tubes with diameter equivalent to common carotid artery and aorta, and a patient-specific model of 27 different LVAD configurations. Results from our analysis demonstrate how energy dissipation is governed by key parameters such as frequency and pulsation, wall elasticity, and LVAD outflow graft surgical anastomosis. We find that frequency, pulsation, and surgical angles have a dominant effect, while wall elasticity has a weaker effect, in determining the state of energy dissipation. For the patient-specific scenario, we also find that energy dissipation is higher in the aortic arch and lower in the abdominal aorta, when compared to the baseline flow without an LVAD. This further illustrates the key hemodynamic role played by the LVAD outflow jet impingement, and subsequent aortic hemodynamics during LVAD operation.

Tension controlled hollow-fiber winding machine for blood oxygenator prototypes

Blood oxygenators involve a complex network of hollow fibers for efficient gas exchange with blood. The optimal microstructural arrangement of these fibers is an ongoing research interest. While the fiber systems of commercial oxygenators are manufactured to address mass production, the research oxygenator prototypes demand more flexibility so that different design parameters can be tested. Here a hollow-fiber assembly system is designed and built for winding research grade extracorporeal blood oxygenator mandrels at different layout dimensions so that these different configurations can be evaluated for mass transfer capacity and blood damage. The hardware design and manufacturing details of this system presented together with its impact on the prototype oxygenator device assembly process. This in-house built system can wind thin fibers, having outer diameters ranging from 100 μm to 1 mm, at any specified winding angle continuously. A control system for fiber stress is also incorporated to eliminate fiber damage. Our system consists of three main units: (1) unwinding, (2) accumulator, and (3) winding systems, integrated together via the control software. The unwinding unit has a PID controller to maintain the position of the accumulator motor on the reference point by tuning the velocity of feeding fibers to the accumulator unit. Another PID controller preserves the desired tension value of the fibers by adjusting the position of the accumulator motor. Desired tension value is defined by the user and typically obtained through uniaxial testing of fibers. The control unit employs a "cascaded" PID controller since the PID controller in the accumulator unit maintains the tension and the PID controller in the unwinding unit controls the position of the accumulator motor. Finally, the winding unit utilizes two motors to wind the fibers over the outer diameter of a mandrel at the desired winding angle. The first motor drives the translational movement, and the second one provides mandrel rotation. The desired angles are achieved by tuning the synchronous movement of the winding motors. While the system is designed to produce assembled blood oxygenator mandrel prototypes, this concept is also applicable for producing cylindrical fiber-reinforced composite materials with specified fiber angles and stents winded on jigs.

Hemodynamic investigation of a novel rotary displacement blood pump for extracorporeal membrane oxygenation

Extracorporeal membrane oxygenation (ECMO) is a life support system used in the treatment of severe respiratory and circulatory failure. High shear stress caused by the high rotational speed of centrifugal blood pumps can cause hemolysis and platelet activation, which are among the major factors leading to the complications of the ECMO system. In this study, a novel blood pump named rotary displacement blood pump (RDBP), which can considerably reduce rotational speed and shear stress while ensuring the normal pressure flow relationship, was proposed. We employed computational fluid dynamics (CFD) analysis to investigate the performance of RDBP under adult ECMO support operating conditions (5 L/min with 350 mmHg). The efficiency and H-Q curves of the RDBP were calculated to evaluate its hydraulic performance, and pressure, flow patterns, and shear stress distribution were analyzed to estimate the hemodynamic characteristics in the pump. In addition, the modified index of hemolysis (MIH) was calculated for the RDBP based on a Eulerian approach. The hydraulic efficiency of the RDBP was 47.28%. The velocity distribution of flow field in the pump was relatively uniform. Most of the liquid (more than 75%) in the pump was exposed to low scale shear stress (<1 Pa), which was close to normal physiological conditions. The gap area was the main distribution location of high scale shear stress. The high wall shear stress (>9 Pa) volume fraction of the RDBP was small and located in the boundary areas between the rotor's edge and the housing. The MIH value of the RDBP was 9.87 ± 0.93 (mean ± SD). The RDBP can achieve better hydraulic efficiency and hemodynamic performance at lower rotational speed. The design of this novel pump is expected to provide a new direction for developing a blood pump for ECMO.

Kinetic and Dynamic Effects on Degradation of von Willebrand Factor

The loss of high molecular weight multimers (HMWM) of von Willebrand factor (vWF) in aortic stenosis (AS) and continuous-flow left ventricular assist devices (cf-LVADs) is believed to be associated with high turbulent blood shear. The objective of this study is to understand the degradation mechanism of HMWM in terms of exposure time (kinetic) and flow regime (dynamics) within clinically relevant pathophysiologic conditions. A custom high-shear rotary device capable of creating fully controlled exposure times and flows was used. The system was set so that human platelet-poor plasma flowed through at 1.75 ml/sec, 0.76 ml/sec, or 0.38 ml/sec resulting in the exposure time ( texp ) of 22, 50, or 100 ms, respectively. The flow was characterized by the Reynolds number (Re). The device was run under laminar (Re = 1,500), transitional (Re = 3,000; Re = 3,500), and turbulent (Re = 4,500) conditions at a given texp followed by multimer analysis. No degradation was observed at laminar flow at all given texp . Degradation of HMWM at a given texp increases with the Re. Re ( p < 0.0001) and texp ( p = 0.0034) are significant factors in the degradation of HMWM. Interaction between Re and texp , however, is not always significant ( p = 0.73).

Hemolysis performance analysis and a novel estimation model of roller pump system

OBJECTIVE

Hemolysis performance is a crucial criterion for roller pumps utilized in life supporting system. In this study, the factor of hemolysis for roller pumps was selected as the target, and an estimation formulation was built to evaluate its hemolysis.

METHODS

Several models were proposed and then simulated with the assistant of Computational fluid dynamics (CFD) framework. The hemolysis performance was calculated using the power law model based on CFD and the estimation model in accordance with geometry parameters proposed in this study. The results of the in vitro experiments were compared with the simulation results. Power law model with the lowest error was utilized in following analysis.

RESULTS

As indicated by the simulation result, the rotary speed most significantly affected the hemolysis performance of roller blood pumps, followed by roller number and diameter of tube. The index of hemolysis (IH) for roller blood pumps at a rotary speed of 20-100 rpm ranged from 8.73E-7 to 8.07E-5. The relative error of the estimation model (4.93%) was lower than of the power law model (6.78%).

CONCLUSION

The IH led by pumps shows a significant, nonlinear relationship with the rotary speed. The design of multiple rollers design is harmful for hemolysis performance and larger diameter of tube exhibits decreased hemolysis at constant flow rate. An estimation formula was proposed with lower relative error for roller pump with the same shell set, which exhibited reduced computation and elevated convenience. And it can be utilized in hemolysis estimation of roller pumps potentially.

Numerical analysis of different cardiac functions and support modes on blood damage potential in an axial pump

BACKGROUND

Cardiac functions and support modes of left ventricular assist device (LVAD) will influence the pump inner flow field and blood damage potential.

METHODS

Computational fluid dynamics (CFD) method and lumped-parameter-model (LPM) were applied to investigate the impacts of cardiac functions under full (9000 rpm) and partial (8000 rpm) support modes in an axial pump.

RESULTS

The constitution of hemolysis index (HI) in different components of the pump was investigated. HI was found to be more sensitive to positive incidence angles (i) compared with negative incidence angles in rotors. Negative incidence angles had little impact on HI both in rotors and the outlet guide vanes. The improved cardiac function made only a minor difference in HIave (estimated average HI in one cardiac cycle) by 9.88%, as the flow rate expanded mainly to higher flow range. Switching to partial support mode, however, would induce a periodic experience of severe flow separation and recirculation at low flow range. This irregular flow field increased HIave by 47.97%, remarkably increasing the blood damage potential.

CONCLUSION

This study revealed the relationship between the blade incidence angle i and HI, and recommended negative-incidence-angle blade designs as it yielded lower HI. Moreover, to avoid flow range below 50% of the design point, careful evaluations should be made before switching support modes as weaning procedures in clinical applications.

Dynamic characteristic modeling of left ventricular assist devices based on hysteresis effects

OBJECTIVE

The purpose of this study is to develop a new model for the dynamic characteristics of left ventricular assist devices (LVADs) interacting with the cardiovascular system under constant-speed modes.

METHODS

A new hysteresis model is established on the basis of the hysteresis effect and turbomachinery principles. The simulation results from the hysteresis model were compared with the inertia model. The in-vitro experiment results of a centrifugal pump (from literature) and the unsteady computational fluid dynamics (CFD) simulation results of an axial pump were used as the benchmarks.

RESULTS

Compared with the inertia model, at the partial support mode, the relative estimation error of the time to the maximum and minimum pump flow (Q) in the hysteresis model decreased at least 16.3% cardiac cycle (T_c) in the centrifugal pump and at least 1.9% T_c in the axial pump, indicating its ability to simulate more realistic Q fluctuations. Moreover, the hysteresis model could predict an accurate time distribution of different Q.

CONCLUSION

The hysteresis model provides a general calculation method for simulating the dynamic characteristics of constant-speed LVADs under interaction with the cardiovascular system. It is more accurate than the inertia model.

SIGNIFICANCE

The hysteresis model is helpful for the rapid estimation of unsteady dynamic characteristics in absence of a physical pump prototype at the preliminary design stage.

Quantitative Assessment of Aortic Hemodynamics for Varying Left Ventricular Assist Device Outflow Graft Angles and Flow Pulsation

Left ventricular assist devices (LVADs) comprise a primary treatment choice for advanced heart failure patients. Treatment with LVAD is commonly associated with complications like stroke and gastro-intestinal (GI) bleeding, which adversely impacts treatment outcomes, and causes fatalities. The etiology and mechanisms of these complications can be linked to the fact that LVAD outflow jet leads to an altered state of hemodynamics in the aorta as compared to baseline flow driven by aortic jet during ventricular systole. Here, we present a framework for quantitative assessment of aortic hemodynamics in LVAD flows realistic human vasculature, with a focus on quantifying the differences between flow driven by LVAD jet and the physiological aortic jet when no LVAD is present. We model hemodynamics in the aortic arch proximal to the LVAD outflow graft, as well as in the abdominal aorta away from the LVAD region. We characterize hemodynamics using quantitative descriptors of flow velocity, stasis, helicity, vorticity and mixing, and wall shear stress. These are used on a set of 27 LVAD scenarios obtained by parametrically varying LVAD outflow graft anastomosis angles, and LVAD flow pulse modulation. Computed descriptors for each of these scenarios are compared against the baseline flow, and a detailed quantitative characterization of the altered state of hemodynamics due to LVAD operation (when compared to baseline aortic flow) is compiled. These are interpreted using a conceptual model for LVAD flow that distinguishes between flow originating from the LVAD outflow jet (and its impingement on the aorta wall), and flow originating from aortic jet during aortic valve opening in normal physiological state.

Hemocompatibility and hemodynamic comparison of two centrifugal LVADs: HVAD and HeartMate3

Mechanical circulatory support using ventricular assist devices is a common technique for treating patients suffering from advanced heart failure. The latest generation of devices is characterized by centrifugal turbopumps which employ magnetic levitation bearings to ensure a gap clearance between moving and static parts. Despite the increasing use of these devices as a destination therapy, several long-term complications still exist regarding their hemocompatibility. The blood damage associated with different pump designs has been investigated profoundly in the literature, while the hemodynamic performance has been hardly considered. This work presents a novel comparison between the two main devices of the latest generation-HVAD and HM3-from both perspectives, hemodynamic performance and blood damage. Computational fluid dynamics simulations are performed to model the considered LVADs, and computational results are compared to experimental measurements of pressure head to validate the model. Enhanced performance and hemocompatibility are detected for HM3 owing to its design incorporating more conventional blades and larger gap clearances.

Numerical assessment of hemodynamic perspectives of a left ventricular assist device and subsequent proposal for improvisation

Due to the unavailability of donors, the use of left ventricular assist devices has emerged to be a reliable line of alternative treatment for heart failure. However, ventricular assist devices (VAD) have been associated with several postoperative complications such as thrombosis, hemolysis, etc. Despite considerable improvements in technology, blood trauma due to high shear stress generation has been a major concern that is largely related to the geometrical feature of the VAD. This study aims to establish the design process of a centrifugal pump by considering several variations in the geometrical feature of a base design using the commercial solver ANSYS-CFX. To capture the uncertain behavior of blood as fluid, Newtonian, as well as non-Newtonian (Bird-Carreau model), models are used for flow field prediction. To assess the possibility of blood damage maximum wall shear stress and hemolysis index have been estimated for each operating point. The results of the simulations yield an optimized design of the pump based on parameters like pressure head generation, maximum shear stress, hydraulic efficiency, and hemolysis index. Further, the design methodology and the steps of development discussed in the paper can serve as a guideline for developing small centrifugal pumps handling blood.

Multi-indicator analysis of mechanical blood damage with five clinical ventricular assist devices

PURPOSE

Device-induced blood damage contributes the hemolysis, thrombosis and bleeding complications in patients supported with ventricular assist device (VAD). This study aims to use a multi-indicator method to understand how devices causes blood damage and identify the "hot spots" of blood trauma within VADs.

METHODS

Computational fluid dynamics (CFD) methods were chosen to investigate the hemodynamic features of five clinical VADs (Impella 5.0, UltraMag, CHVAD, HVAD, and HeartMate II) under the same clinical support condition (flow rate of 4.5L/min, pressure head around 75 mmHg). A comprehensive multi-indicator evaluation method including hemodynamic parameters, hemolysis model, thrombotic potential model and bleeding probability model was used to analyze blood damage and assess the hemodynamic performance and hemocompatibility of these VADs.

RESULTS

Simulation results show that shear stress from 50 Pa to 100 Pa plays a major role in blood damage in Impella 5.0, UltraMag and CHVAD, while blood damage in HVAD and HeartMate II is mainly caused by shear stress greater than 100 Pa. Residence time was not the main factor for blood damage in Impella 5.0, and also makes a limited contribution to blood trauma in UltraMag and CHVAD, while it takes a critical role in elevating thrombotic potential in HVAD and HeartMate II. The distribution of regions of high hemolysis risk and high bleeding probability was similar for all these VADs and partially overlapped for high thrombotic potential regions. For Impella 5.0, regions with high hemolysis and bleeding risk were found mainly in the blade tip clearance and diffuser domains, high thrombotic potential regions were almost absent. For UltraMag, regions with high hemolysis, bleeding and thrombosis potential were found in two corners of the inlet pipe, the secondary flow passage, and the impeller eye. For CHVAD, the high-risk regions for hemolysis, bleeding and thrombosis are mainly in the inner side of the secondary flow passage and the middle region of the impeller passage. The narrow hydrodynamic clearance and impeller passage had a high risk of hemolysis and bleeding, and the clearance between the rotor and guide cone and the hydrodynamic clearance had high thrombotic potential. For HeartMate II, regions of high hemolysis risk and bleeding probability were found in the near-wall region of the straightener, the blade tip clearance and the diffuser domain. The corners of the inlet and outlet pipe and the straightener and diffuser regions had high thrombotic potential.

CONCLUSION

The risk of hemolysis, bleeding and thrombosis for these five VADs, in increasing order, was Impella 5.0, UltraMag, CHVAD, HVAD, and HeartMate II. Flow losses caused by the rotor mechanical movement, chaotic flow and narrow clearances increase the blood damage for all these VADs. The multi-indicator analysis can comprehensively evaluate the VAD performance with improved assessment accuracy of CFD.

Computations of the shear stresses distribution experienced by passive particles as they circulate in turbulent flow: A case study for vWF protein molecules

The stress distribution along the trajectories of passive particles released in turbulent flow were computed with the use of Lagrangian methods and direct numerical simulations. The flow fields selected were transitional Poiseuille-Couette flow situations found in ventricular assist devices and turbulent flows at conditions found in blood pumps. The passive particle properties were selected to represent molecules of the von Willebrand factor (vWF) protein. Damage to the vWF molecule can cause disease, most often related to hemostasis. The hydrodynamic shear stresses along the trajectories of the particles were calculated and the changes in the distribution of stresses were determined for proteins released in different locations in the flow field and as a function of exposure time. The stress distributions indicated that even when the average applied stress was within a safe operating regime, the proteins spent part of their trajectories in flow areas of damaging stress. Further examination showed that the history of the distribution of stresses applied on the vWF molecules, rather than the average, should be used to evaluate hydrodynamically-induced damage.

Development of the Centrifugal Blood Pump for a Hybrid Continuous Flow Pediatric Total Artificial Heart: Model, Make, Measure

Clinically-available blood pumps and total artificial hearts for pediatric patients continue to lag well behind those developed for adults. We are developing a hybrid, continuous-flow, magnetically levitated, pediatric total artificial heart (TAH). The hybrid TAH design integrates both an axial and centrifugal blood pump within a single, compact housing. The centrifugal pump rotates around the separate axial pump domain, and both impellers rotate around a common central axis. Here, we concentrate our development effort on the centrifugal blood pump by performing computational fluid dynamics (CFD) analysis of the blood flow through the pump. We also conducted transient CFD analyses (quasi-steady and transient rotational sliding interfaces) to assess the pump's dynamic performance conditions. Through modeling, we estimated the pressure generation, scalar stress levels, and fluid forces exerted on the magnetically levitated impellers. To further the development of the centrifugal pump, we also built magnetically-supported prototypes and tested these in an in vitro hydraulic flow loop and via 4-h blood bag hemolytic studies (n = 6) using bovine blood. The magnetically levitated centrifugal prototype delivered 0–6.75 L/min at 0–182 mmHg for 2,750–4,250 RPM. Computations predicted lower pressure-flow performance results than measured by testing; axial and radial fluid forces were found to be <3 N, and mechanical power usage was predicted to be <5 Watts. Blood damage indices (power law weighted exposure time and scalar stress) were <2%. All data trends followed expectations for the centrifugal pump design. Six peaks in the pressure rise were observed in the quasi-steady and transient simulations, correlating to the blade passage frequency of the 6-bladed impeller. The average N.I.H value (n = 6) was determined to be 0.09 ± 0.02 g/100 L, which is higher than desired and must be addressed through design improvement. These data serve as a strong foundation to build upon in the next development phase, whereby we will integrate the axial flow pump component.

Left ventricular assist device and pump thrombosis: the importance of the inflow cannula position

Pump thrombosis is a devastating complication after left ventricular assist device implantation. This study aims to elucidate the relation between left ventricular assist device implantation angle and risk of pump thrombosis. Between November 2010 and March 2020, 53 left ventricular assist device-patients underwent a computed tomography scan. Using a 3-dimensional multiplanar reformation the left ventricular axis was reconstructed to measure the implantation angle of the inflow cannula. All patients were retrospectively analyzed for the occurrence of pump thrombosis. In 10 (91%) patients with a pump thrombosis, the implantation angle was towards the lateral wall of the left ventricle. In only 20 patients (49%) of the patients without a pump thrombosis the inflow cannula pointed towards the lateral wall of the left ventricle. The mean angle in patients with a pump thrombosis was 10.1 ± 11.9 degrees towards the lateral wall of the left ventricle compared to 4.1 ± 19.9 degrees towards the septum in non-pump thrombosis patients (P = 0.005). There was a trend towards a significant difference in time to first pump thrombosis between patients with a lateral or septal deviated left ventricular assist device (hazard ratio of 0.15, P = 0.07). This study demonstrates that left ventricular assist device implantation angle is associated with pump thrombosis. Almost all patients in whom a pump thrombosis occurred during follow-up had a left ventricular assist device implanted with the inflow-cannula pointing towards the lateral wall of the left ventricle.

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.

In vitro study of red blood cell and VWF damage in mechanical circulatory support devices based on blood-shearing platform

Blood damage induced by mechanical circulatory support devices (MCSDs) remains a significant challenge to optimal clinical care. Although researchers have been conducting in vitro studies, the major determinant of blood damage is still unclear. An optimized capillary tube blood-shearing platform with custom designed parts was constructed to investigate the influence of two flow-dependent parameters (shear stress and exposure time) on the shear-induced damage of red blood cells and von Willebrand factor (VWF). Blood samples under different high shear stress and instantaneous exposure time were obtained by changing the flow rate and the length of capillary tube. Plasma free hemoglobin assay and immunoblotting of VWF were then performed on the sheared blood samples. The quantitative correlation between the hemolysis index and the two flow-dependent parameters was found following the power law mathematical model under the flow condition with high shear stress and instantaneous exposure time. The degradation of high molecular weight VWF was not obvious under high shear stress factor. However, the degradation of high molecular weight VWF was found as the result of the accumulation over exposure time under non-physiological shear stress, which was consistent with the different mechanism of VWF damage comparing to red blood cell damage. Compared to peak shear stress, exposure time has a greater effect on both red blood cell and VWF damage. To improve the hemocompatibility of MCSDs, it is more important to avoid regions of slow blood flow with non-physiological shear stress under laminar flow conditions.

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.

CFD analysis of the HVAD's hemodynamic performance and blood damage with insight into gap clearance

Mechanical circulatory support using ventricular assist devices has become commonplace in the treatment of patients suffering from advanced stages of heart failure. While blood damage generated by these devices has been evaluated in depth, their hemodynamic performance has been investigated much less. This work presents the analysis of the complete operating map of a left ventricular assist device, in terms of pressure head, power and efficiency. Further investigation into its hemocompatibility is included as well. To achieve these objectives, computational fluid dynamics simulations of a centrifugal blood pump with a wide-blade impeller were performed. Several conditions were considered by varying the rotational speed and volumetric flow rate. Regarding the device's hemocompatibility, blood damage was evaluated by means of the hemolysis index. By relating the hemocompatibility of the device to its hemodynamic performance, the results have demonstrated that the highest hemolysis occurs at low flow rates, corresponding to operating conditions of low efficiency. Both performance and hemocompatibility are affected by the gap clearance. An innovative investigation into the influence of this design parameter has yielded decreased efficiencies and increased hemolysis as the gap clearance is reduced. As a further novelty, pump operating maps were non-dimensionalized to highlight the influence of Reynolds number, which allows their application to any working condition. The pump's operating range places it in the transitional regime between laminar and turbulent, leading to enhanced efficiency for the highest Reynolds number.

Hemolytic Footprint of Rotodynamic Blood Pumps

OBJECTIVE

In preclinical examinations, rotodynamic blood pumps (RBPs) are predominantly evaluated at design-point conditions. In clinical practice, however, they run at diversified modes of operation. This study aimed at extending current preclinical evaluation of hemolytic profiles in RBPs toward broader, clinically relevant ranges of operation.

METHODS

Two implantable RBPs the HeartMate 3 (HM3) and the HeartWare Ventricular Assist Device (HVAD) were analyzed at three pump speeds (HM3: 4300, 5600, 7000rpm; HVAD: 1800, 2760, 3600rpm) with three flow rates (1-9L/min) per speed setting. Hemolysis measurements were performed in heparinized bovine blood. The delta free hemoglobin (dfHb) and the normalized index of hemolysis (NIH) served as hemolytic measures. Statistical analysis was performed by multiple comparison of the 9 operating conditions. Moreover, computational fluid dynamics (CFD) was applied to provide mechanistic insights into the interrelation between hydraulics and hemolysis by correlating numerically computed hydraulic losses with in-vitro hemolytic measures.

RESULTS

In both devices, dfHb increased toward increasing speeds, particularly during low but also during high flow condition. By contrast, in both RBPs magnitudes of NIH were significantly elevated during low flow operation compared to high flow conditions (p<0.0036). Maps of hemolytic metrics revealed morphologically similar trends to in-silico hydraulic losses (r>0.793).

CONCLUSIONS

While off-design operation is associated with increased hemolytic profiles, the setting of different operating conditions render a preclinical prediction of clinical impact with current hemolysis metrics difficult.

SIGNIFICANCE

The identified increase in hemolytic measures during episodes of off-design operation is highlighting the need to consider worst-case operation during preclinical examinations.

Hemolytic Performance in Two Generations of the Sputnik Left Ventricular Assist Device: A Combined Numerical and Experimental Study

Background: Currently, left ventricular assist devices (LVADs) are a successful surgical treatment for patients with end-stage heart failure on the waiting list or with contraindicated heart transplantation. In Russia, Sputnik 1 LVAD was also successfully introduced into clinical practice as a bridge-to-transplant and a destination therapy device. Development of Sputnik 2 LVAD was aimed at miniaturization to reduce invasiveness, optimize hemocompatibility, and improve versatility for patients of various sizes. Methods: We compared hemolysis level in flow path of the Sputnik LVADs and investigated design aspects influencing other types of blood damage, using predictions of computational fluid dynamics (CFD) and experimental assessment. The investigated operating point was a flow rate of 5 L/min and a pressure head of 100 mm Hg at an impeller rotational speed of 9100 min⁻¹. Results: Mean hemolysis indices predicted with CFD were 0.0090% in the Sputnik 1 and 0.0023% in the Sputnik 2. Averaged values of normalized index of hemolysis obtained experimentally for the Sputnik 1 and the Sputnik 2 were 0.011 ± 0.003 g/100 L and 0.004 ± 0.002 g/100 L, respectively. Conclusions: Obtained results indicate obvious improvements in hemocompatibility and sufficiently satisfy the determined miniaturization aim for the Sputnik 2 LVAD development.

Development of Inspired Therapeutics Pediatric VAD: Computational Analysis and Characterization of VAD V3

PURPOSE

Pediatric heart failure patients remain in critical need of a dedicated mechanical circulatory support (MCS) solution as development efforts for specific pediatric devices continue to fall behind those for the adult population. The Inspired Pediatric VAD is being developed as a pediatric specific MCS solution to provide up to 30-days of circulatory or respiratory support in a compact modular package that could allow for patient ambulation during treatment.

METHODS

Hydrodynamic performance (flows, pressures), impeller/rotor mechanical properties (torques, forces), and flow shear stress and residence time distributions of the latest design version, Inspired Pediatric VAD V3, were numerically predicted and investigated using computational fluid dynamics (CFD) software (SolidWorks Flow Simulator).

RESULTS

Hydrodynamic performance was numerically predicted, indicating no change in flow and pressure head compared to the previous device design (V2), while displaying increased impeller/rotor torques and translation forces enabled by improved geometry. Shear stress and flow residence time volumetric distributions are presented over a range of pump rotational speeds and flow rates. At the lowest pump operating point (3000 RPM, 0.50 L/min, 75 mmHg), 79% of the pump volume was in the shear stress range of 0-10 Pa with < 1% of the volume in the critical range of 150-1000 Pa for blood damage. At higher speed and flow (5000 RPM, 3.50 L/min, 176 mmHg), 65% of the volume resided in the 0-10 Pa range compared to 2.3% at 150-1000 Pa.

CONCLUSIONS

The initial computational characterization of the Inspired Pediatric VAD V3 is encouraging and future work will include device prototype testing in a mock circulatory loop and acute large animal model.

Distribution and history of extensional stresses on vWF surrogate molecules in turbulent flow

The configuration of proteins is critical for their biochemical behavior. Mechanical stresses that act on them can affect their behavior leading to the development of decease. The von Willebrand factor (vWF) protein circulating with the blood loses its efficacy when it undergoes non-physiological hemodynamic stresses. While often overlooked, extensional stresses can affect the structure of vWF at much lower stress levels than shear stresses. The statistical distribution of extensional stress as it applies on models of the vWF molecule within turbulent flow was examined here. The stress on the molecules of the protein was calculated with computations that utilized a Lagrangian approach for the determination of the molecule trajectories in the flow filed. The history of the stresses on the proteins was also calculated. Two different flow fields were considered as models of typical flows in cardiovascular mechanical devises, one was a Poiseuille flow and the other was a Poiseuille–Couette flow field. The data showed that the distribution of stresses is important for the design of blood flow devices because the average stress can be below the critical value for protein damage, but tails of the distribution can be outside the critical stress regime.

Hemolysis Testing In Vitro: A Review of Challenges and Potential Improvements

Many medical devices such as cardiopulmonary bypass systems, mechanical heart valves, or ventricular assist devices are intended to come into contact with blood flow during use. In vitro hemolysis testing can provide valuable information about the hemocompatibility of prototypes and thus help reduce the number of animal experiments required. Such tests play an important role as research and development tools for objective comparisons of prototypes and devices as well as for the extrapolation of their results to clinical outcomes. Therefore, it is important to explore and provide new ways to improve current practices. In this article, the main challenges of hemolysis testing are described, namely the difficult blood sourcing, the high experimental workload, and the low reproducibility of test results. Several approaches to address the challenges identified are proposed and the respective literature is reviewed. These include the replacement of blood as the "shear-sensitive fluid" by alternative test fluids, the replacement of sparse, manual sampling and blood damage assessment by a continuous and automated monitoring, as well as an analysis of categories and causes of variability in hemolysis test results that may serve as a structural template for future studies.

Blood Pumps for Extracorporeal Membrane Oxygenation: Platelet Activation During Different Operating Conditions

Extracorporeal membrane oxygenation (ECMO) is a therapy used in severe cardiopulmonary failure. Blood is pumped through an artificial circuit exposing it to nonphysiologic conditions, which promote platelet activation and coagulation. Centrifugal pumps used at lower flow rates than their design point may lose pump efficiency and increase the risk of hemolysis. In this study, thrombogenic properties of two ECMO pumps designed for adult and neonatal use were evaluated using simulations in different flow scenarios. Three scenarios, adult pump in adult mode (4 L/min), adult pump in baby mode (300 ml/min), and neonatal pump used in its design point (300 ml/min), were simulated using computational fluid dynamics. The flow was numerically seeded with platelets, whose activation state was computed considering the stress history that acted along their respective path lines. Statistical distributions of activation state and residence time were drawn. The results showed that using the adult pump in baby mode increased the fraction of platelets with higher activation state confirming that low-pump flow rate impacts thrombogenicity. The neonatal pump showed a backflow at the inlet, which carried platelets in a retrograde motion contributing to an increased thrombogenic potential compared with the adult mode scenario.

Physiology of Continuous-Flow Left Ventricular Assist Device Therapy

The expanding use of continuous-flow left ventricular assist devices (CF-LVADs) for end-stage heart failure warrants familiarity with the physiologic interaction of the device with the native circulation. Contemporary devices utilize predominantly centrifugal flow and, to a lesser extent, axial flow rotors that vary with respect to their intrinsic flow characteristics. Flow can be manipulated with adjustments to preload and afterload as in the native heart, and ascertainment of the predicted effects is provided by differential pressure-flow (H-Q) curves or loops. Valvular heart disease, especially aortic regurgitation, may significantly affect adequacy of mechanical support. In contrast, atrioventricular and ventriculoventricular timing is of less certain significance. Although beneficial effects of device therapy are typically seen due to enhanced distal perfusion, unloading of the left ventricle and atrium, and amelioration of secondary pulmonary hypertension, negative effects of CF-LVAD therapy on right ventricular filling and function, through right-sided loading and septal interaction, can make optimization challenging. Additionally, a lack of pulsatile energy provided by CF-LVAD therapy has physiologic consequences for end-organ function and may be responsible for a series of adverse effects. Rheological effects of intravascular pumps, especially shear stress exposure, result in platelet activation and hemolysis, which may result in both thrombotic and hemorrhagic consequences. Development of novel solutions for untoward device-circulatory interactions will facilitate hemodynamic support while mitigating adverse events. © 2021 American Physiological Society. Compr Physiol 12:1-37, 2021.

Individual Variability in von Willebrand Factor Fragility in Response to Shear Stress: A Possible Clue for Predicting Bleeding Risk

Acquired von Willebrand syndrome (AVWS), characterized by reduced von Willebrand factor (VWF) large multimers, has recently been implicated as the principal mechanism underlying bleeding in patients implanted with left ventricular assist devices (LVADs). Hematological severity of AVWS varies among patients, even if an identical device is implanted. We investigated whether this diversity in hematological severity is due to individual variability in VWF fragility, according to responses to incremental shear stress. Whole-blood samples were sheared at 20,000-40,000 s-1 shear rate, an index of shear stress, using a custom-made shear stressor that could generate shear stress compatible with that produced by an LVAD. The degree of VWF large multimers degradation was evaluated using the VWF large multimer index. A significant inverse correlation was observed between the VWF large multimer index and LVAD-compatible magnitudes of shear stress: the VWF large multimer indices were 68.5 ± 18.3, 48.0 ± 13.9, 33.9 ± 12.1, 23.7 ± 7.9, and 18.7% ± 8.7% at 20,000, 25,000, 30,000, 35,000, and 40,000 s-1 of shear rates, respectively (P < 0.0001). Furthermore, experimental VWF large multimer index values were compatible with those derived from patients with implanted LVADs (median; 28.9%). Finally, reduction in the VWF large multimer index corresponding to shear stress showed individual variation. We demonstrated that the combined use of a novel high shear stress loading device and quantitative evaluation of VWF large multimers may predict risk of bleeding before LVAD implantation.

Experimental Hemodynamics within the Penn State Fontan Circulatory Assist Device

For children born with a single functional ventricle, the Fontan operation bypasses the right ventricle by forming a four-way total cavopulmonary connection adapting the existing ventricle for the systemic circulation. However, upon adulthood, many Fontan patients exhibit low cardiac output and elevated venous pressure, eventually requiring a heart transplantation. Despite efforts to develop a Fontan pump or use an existing ventricular assist device for failing Fontan support, there is still no device designed or tested for subpulmonary support. Penn State University is developing a hydrodynamically levitated Fontan circulatory assist device (FCAD) for bridge-to-transplant or destination therapy. The FCAD hemodynamics, at both steady and pulsatile conditions for three pump operating conditions, were quantified using particle image velocimetry to determine the velocity magnitudes and Reynolds normal and shear stresses. Data were acquired at three planes (0 mm and ±25% of the radius) for the inferior and superior vena cavae inlets and the pulmonary artery outlet. The inlets had a blunt velocity profile that became skewed towards the collecting volute as fluid approached the rotor. At the outlet, regardless of the flow condition, a high-velocity jet exited the volute and moved downstream in a helical pattern. Turbulent stresses observed at the volute exit were influenced by the rotor's rotation. Regardless of inlet conditions, the pump demonstrated advantageous behavior for clinical use with a predictable flow field and a low risk of platelet adhesion and hemolysis based on calculated wall shear rates and turbulent stresses, respectively.

Hemocompatibility Assessment Platform Drive System Design: Trade-off between Motor Performance and Hemolysis

Left ventricular assist devices (LVADs) have long been used to treat adults with heart failure, but LVAD options for pediatric patients with heart failure are lacking. Despite the urgent need for long-term, implantable pediatric LVADs, design challenges such as hemolysis, pump thrombosis, and bleeding persist. We have developed a Hemocompatibility Assessment Platform (HAP) to identify blood trauma from individual LVAD components. A HAP would aid in refining pump components before in vivo testing, thereby preventing unnecessary animal sacrifice and reducing development time and cost. So that the HAP does not confound hemolysis data, the HAP drive system consists of an enlarged air-gap motor coupled to a magnetic levitation system. Although it is known that an enlarged air gap motor will have diminished performance, while the larger gap in the motor will cause less blood damage, the trade-offs are not fully characterized. Therefore, in this study we evaluated these trade-offs to determine an optimal rotor diameter for the HAP drive motor. The motor performance was characterized with an experimental method by determining the torque constant for the HAP drive motor with varied rotor diameters. The torque threshold was set as 10 mNm to achieve a nominal current of 3.5A. Hemolysis in the HAP drive motor gap was estimated by calculating scalar shear stress generated in the HAP motor gap analytically and numerically. A design criterion of 30 Pa was selected for scalar shear stress to achieve minimal hemolysis and platelet activation in the HAP drive system.Clinical Relevance- We evaluated a Hemocompatibility Assessment Platform for developing LVAD prototypes that can best balance motor performance and hemocompatibility. This design method can assist with optimizing the drive system during the research stage and illustrates how motor geometry can be tuned to reduce blood trauma.

Mock circulatory loops used for testing cardiac assist devices: A review of computational and experimental models

Heart failure is a major health risk, and with limited availability of donor organs, there is an increasing need for developing cardiac assist devices (CADs). Mock circulatory loops (MCL) are an important in-vitro test platform for CAD's performance assessment and optimisation. The MCL is a lumped parameter model constructed out of hydraulic and mechanical components aiming to simulate the native cardiovascular system (CVS) as closely as possible. Further development merged MCLs and numerical circulatory models to improve flexibility and accuracy of the system; commonly known as hybrid MCLs. A total of 128 MCLs were identified in a literature research until 25 September 2020. It was found that the complexity of the MCLs rose over the years, recent MCLs are not only capable of mimicking the healthy and pathological conditions, but also implemented cerebral, renal and coronary circulations and autoregulatory responses. Moreover, the development of anatomical models made flow visualisation studies possible. Mechanical MCLs showed excellent controllability and repeatability, however, often the CVS was overly simplified or lacked autoregulatory responses. In numerical MCLs the CVS is represented with a higher order of lumped parameters compared to mechanical test rigs, however, complex physiological aspects are often simplified. In hybrid MCLs complex physiological aspects are implemented in the hydraulic part of the system, whilst the numerical model represents parts of the CVS that are too difficult to represent by mechanical components per se. This review aims to describe the advances, limitations and future directions of the three types of MCLs.

Large eddy simulation as a fast and accurate engineering approach for the simulation of rotary blood pumps

An accurate representation of the flow field in blood pumps is important for the design and optimization of blood pumps. The primary turbulence modeling methods applied to blood pumps have been the Reynolds-averaged Navier-Stokes (RANS) or URANS (unsteady RANS) method. Large eddy simulation (LES) method has been introduced to simulate blood pumps. Nonetheless, LES has not been widely used to assist in the design and optimization of blood pumps to date due to its formidable computational cost. The purpose of this study is to explore the potential of the LES technique as a fast and accurate engineering approach for the simulation of rotary blood pumps. The performance of "Light LES" (using the same time and spatial resolutions as the URANS) and LES in two rotary blood pumps was evaluated by comparing the results with the URANS and extensive experimental results. This study showed that the results of both "Light LES" and LES are superior to URANS, in terms of both performance curves and key flow features. URANS could not predict the flow separation and recirculation in diffusers for both pumps. In contrast, LES is superior to URANS in capturing these flows, performing well for both design and off-design conditions. The differences between the "Light LES" and LES results were relatively small. This study shows that with less computational cost than URANS, "Light LES" can be considered as a cost-effective engineering approach to assist in the design and optimization of rotary blood pumps.

Design and Computational Evaluation of a Pediatric MagLev Rotary Blood Pump

Pediatric heart failure (HF) patients have been a historically underserved population for mechanical circulatory support (MCS) therapy. To address this clinical need, we are developing a low cost, universal magnetically levitated extracorporeal system with interchangeable pump heads for pediatric support. Two impeller and pump designs (pump V1 and V2) for the pediatric pump were developed using dimensional analysis techniques and classic pump theory based on defined performance criteria (generated flow, pressure, and impeller diameter). The designs were virtually constructed using computer-aided design (CAD) software and 3D flow and pressure features were analyzed using computational fluid dynamics (CFD) analysis. Simulated pump designs (V1, V2) were operated at higher rotational speeds (~5,000 revolutions per minute [RPM]) than initially estimated (4,255 RPM) to achieve the desired operational point (3.5 L/min flow at 150 mm Hg). Pump V2 outperformed V1 by generating approximately 30% higher pressures at all simulated rotational speeds and at 5% lower priming volume. Simulated hydrodynamic performance (achieved flow and pressure, hydraulic efficiency) of our pediatric pump design, featuring reduced impeller size and priming volume, compares favorably to current commercially available MCS devices.

Video-based valve motion combined with computational fluid dynamics gives stable and accurate simulations of blood flow in the Realheart total artificial heart

BACKGROUND

Patients with end-stage, biventricular heart failure, and for whom heart transplantation is not an option, may be given a Total Artificial Heart (TAH). The Realheart® is a novel TAH which pumps blood by mimicking the native heart with translation of an atrioventricular plane. The aim of this work was to create a strategy for using Computational Fluid Dynamics (CFD) to simulate haemodynamics in the Realheart®, including motion of the atrioventricular plane and valves.

METHODS

The accuracies of four different computational methods for simulating fluid-structure interaction of the prosthetic valves were assessed by comparison of chamber pressures and flow rates with experimental measurements. The four strategies were: prescribed motion of valves opening and closing at the atrioventricular plane extrema; simulation of fluid-structure interaction of both valves; prescribed motion of the mitral valve with simulation of fluid-structure interaction of the aortic valve; motion of both valves prescribed from video analysis of experiments.

RESULTS

The most accurate strategy (error in ventricular pressure of 6%, error in flow rate of 5%) used video-prescribed motion. With the Realheart operating at 80 bpm, the power consumption was 1.03 W, maximum shear stress was 15 Pa, and washout of the ventricle chamber after 4 cycles was 87%.

CONCLUSIONS

This study, the first CFD analysis of this novel TAH, demonstrates that good agreement between computational and experimental data can be achieved. This method will therefore enable future optimisation of the geometry and motion of the Realheart®.

The sheep as a pre-clinical model for testing intra-aortic percutaneous mechanical circulatory support devices

The save deployment of intra-aortic percutaneous mechanical circulatory support devices is highly dependent on the inner aortic diameter. Finding the anatomically and ethically most suitable animal model for performance testing of new pMCS devices remains challenging. For this study, an ovine model using adult ewes of a large framed breed (Swiss White Alpine Sheep) was developed to test safety, reliability, and biocompatibility of catheter-mounted mechanical support devices placed in the descending thoracic aorta. Following the drawback of fluctuating aortic diameter and device malfunction in the first four animals, the model was improved by stenting the following animals with an aortic stent. Stenting the animals with an intra-aortic over the balloon stent was found to standardize the experimental set-up and to avoid early termination of the experiment due to non-device related issues.

Stochastic simulation of the FDA centrifugal blood pump benchmark

In the present study, the effect of physical and operational uncertainties on the hydrodynamic and hemocompatibility characteristics of a centrifugal blood pump designed by the U.S. food and drug administration is investigated. Physical uncertainties include the randomness in the blood density and viscosity, while the operational uncertainties are composed of the pump rotational speed, mass flow rate, and turbulence intensity. The non-intrusive polynomial chaos expansion has been employed to conduct the uncertainty quantification analysis. Additionally, to assess each stochastic parameter's influence on the quantities of interest, the sensitivity analysis is utilized through the Sobol' indices. For numerical simulation of the pump's blood flow, the SST [Formula: see text] turbulence model and a power-law model of hemolysis were employed. The pump's velocity field is profoundly affected by the rotational speed in the bladed regions and the mass flow rate in other zones. Furthermore, the hemolysis index is dominantly sensitive to blood viscosity. According to the results, pump hydraulic characteristics (i.e., head and efficiency) show a more robust behavior than the hemocompatibility characteristics (i.e., hemolysis index) regarding the operational and physical uncertainties. Finally, it was found that the probability distribution function of the hemolysis index covers the experimental measurements.