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

Effect of impeller rotational phase on the FDA blood pump velocity fields

BACKGROUND

The Food and Drug Administration (FDA) blood pump is an open-source benchmark cardiovascular device introduced for validating computational and experimental performance analysis tools. The time-resolved velocity field for the whole impeller has not been established, as is undertaken in this particle image velocimetry (PIV) study. The level of instantaneous velocity fluctuations is important, to assess the flow-induced rotor vibrations which may contribute to the total blood damage.

METHODS

To document these factors, time-resolved two-dimensional PIV experiments were performed that were precisely phase-locked with the impeller rotation angle. The velocity fields in the impeller and in the volute conformed with the previous single blade passage experiments of literature.

RESULTS

Depending on the impeller orientation, present experiments showed that volute outlet nozzle flow can fluctuate up to 34% during impeller rotation, with a maximum standard experimental uncertainty of 2.2%. Likewise, the flow fields in each impeller passage also altered in average 33.5%. Considerably different vortex patterns were observed for different blade passages, with the largest vortical structures reaching an average core radii of 7 mm. The constant volute area employed in the FDA pump design contributes to the observed velocity imbalance, as illustrated in our velocity measurements.

CONCLUSIONS

By introducing the impeller orientation parameter for the nozzle flow, this study considers the possible uncertainties influencing pump flow. Expanding the available literature data, analysis of inter-blade relative velocity fields is provided here for the first-time to the best of our knowledge. Consequently, our research fills a critical knowledge gap in the understanding of the flow dynamics of an important benchmark cardiovascular device. This study prompts the need for improved hydrodynamic designs and optimized devices to be used as benchmark test devices, to build more confidence and safety in future ventricular assist device performance assessment studies.

Global sensitivity analysis with multifidelity Monte Carlo and polynomial chaos expansion for vascular haemodynamics

Computational models of the cardiovascular system are increasingly used for the diagnosis, treatment, and prevention of cardiovascular disease. Before being used for translational applications, the predictive abilities of these models need to be thoroughly demonstrated through verification, validation, and uncertainty quantification. When results depend on multiple uncertain inputs, sensitivity analysis is typically the first step required to separate relevant from unimportant inputs, and is key to determine an initial reduction on the problem dimensionality that will significantly affect the cost of all downstream analysis tasks. For computationally expensive models with numerous uncertain inputs, sample-based sensitivity analysis may become impractical due to the substantial number of model evaluations it typically necessitates. To overcome this limitation, we consider recently proposed Multifidelity Monte Carlo estimators for Sobol' sensitivity indices, and demonstrate their applicability to an idealized model of the common carotid artery. Variance reduction is achieved combining a small number of three-dimensional fluid-structure interaction simulations with affordable one- and zero-dimensional reduced-order models. These multifidelity Monte Carlo estimators are compared with traditional Monte Carlo and polynomial chaos expansion estimates. Specifically, we show consistent sensitivity ranks for both bi- (1D/0D) and tri-fidelity (3D/1D/0D) estimators, and superior variance reduction compared to traditional single-fidelity Monte Carlo estimators for the same computational budget. As the computational burden of Monte Carlo estimators for Sobol' indices is significantly affected by the problem dimensionality, polynomial chaos expansion is found to have lower computational cost for idealized models with smooth stochastic response.

The impact of small movements with dual lumen cannulae during venovenous extracorporeal membrane oxygenation: A computational fluid dynamics analysis

BACKGROUND AND OBJECTIVES

Venovenous Extracorporeal Membrane Oxygenation (VV ECMO) provides respiratory support to patients with severe lung disease failing conventional medical therapy. An essential component of the ECMO circuit are the cannulas, which drain and return blood into the body. Despite being anchored to the patient to prevent accidental removal, minor cannula movements are common during ECMO. The clinical and haemodynamic consequences of these small movements are currently unclear. This study investigated the risk of thrombosis and recirculation caused by small movements of a dual lumen cannula (DLC) in an adult using computational fluid dynamics.

METHODS

The 3D model of an AVALON Elite DLC (27 Fr) and a patient-specific vena cava and right atrium were generated for an adult patient on ECMO. The baseline cannula position was generated where the return jet enters the tricuspid valve. Alternative cannula positions were obtained by shifting the cannula 5 and 15 mm towards inferior (IVC) and superior (SVC) vena cava, respectively. ECMO settings of 4 L/min blood flow and pulsatile flow at SVC and IVC were applied. Recirculation was defined as a scalar value indicating the infused oxygenated blood inside the drainage lumen, while thrombosis risk was evaluated by shear stress, stagnation volume, washout, and turbulent kinetic energy.

RESULTS

Recirculation for all models was less than 3.1 %. DLC movements between -5 to 15 mm increased shear stress and turbulence kinetic energy up to 24.7 % and 11.8 %, respectively, compared to the baseline cannula position leading to a higher predicted thrombosis risk. All models obtained a complete washout after nine seconds except for when the cannula migrated 15 mm into the SVC, indicating persisting stasis and circulating zones.

CONCLUSION

In conclusion, small DLC movements were not associated with an increased risk of recirculation. However, they may increase the risk of thrombosis due to increased shear rate, turbulence, and slower washout of blood. Developing effective cannula securement devices may reduce this risk.

Turbulent flow field in maglev centrifugal blood pumps of CH-VAD and HeartMate III: secondary flow and its effects on pump performance

Secondary flow path is one of the crucial aspects during the design of centrifugal blood pumps. Small clearance size increases stress level and blood damage, while large clearance size can improve blood washout and reduce stress level. Nonetheless, large clearance also leads to strong secondary flows, causing further blood damage. Maglev blood pumps rely on magnetic force to achieve rotor suspension and allow more design freedom of clearance size. This study aims to characterize turbulent flow field and secondary flow as well as its effects on the primary flow and pump performance, in two representative commercial maglev blood pumps of CH-VAD and HeartMate III, which feature distinct designs of secondary flow path. The narrow and long secondary flow path of CH-VAD resulted in low secondary flow rates and low disturbance to the primary flow. The flow loss and blood damage potential of the CH-VAD mainly occurred at the secondary flow path, as well as the blade clearances. By contrast, the wide clearances in HeartMate III induced significant disturbance to the primary flow, resulting in large incidence angle, strong secondary flows and high flow loss. At higher flow rates, the incidence angle was even larger, causing larger separation, leading to a significant decrease of efficiency and steeper performance curve compared with CH-VAD. This study shows that maglev bearings do not guarantee good blood compatibility, and more attention should be paid to the influence of secondary flows on pump performance when designing centrifugal blood pumps.

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.

A multi-constituent model for assessing the effect of impeller shroud on the thrombosis potential of a centrifugal blood pump

Centrifugal blood pumps can be used for treating heart failure patients. However, pump thrombosis has remained one of the complications that trouble clinical treatment. This study analyzed the effect of impeller shroud on the thrombosis risk of the blood pump, and predicted areas prone to thrombosis. Multi-constituent transport equations were presented, considering mechanical activation and biochemical activation. It was found that activated platelets concentration can increase with shear stress and adenosine diphosphate(ADP) concentration increasing, and the highest risk of thrombosis inside the blood pump was under extracorporeal membrane oxygenation (ECMO) mode. Under the same condition, ADP concentration and thrombosis index of semi-shroud impeller can increase by 7.3% and 7.2% compared to the closed-shroud impeller. The main reason for the increase in thrombosis risk was owing to elevated scalar shear stress and more coagulation promoting factor-ADP released. The regions with higher thrombosis potential were in the center hole, top and bottom clearance. As a novelty, the findings revealed that impeller shroud can influence mechanical and biochemical activation factors. It is useful for identifying potential risk regions of thrombus formation based on relative comparisons.

In silico parametric analysis of femoro-jugular venovenous ECMO and return cannula dynamics: In silico analysis of femoro-jugular VV ECMO

BACKGROUND

Increasingly, computational fluid dynamics (CFD) is helping explore the impact of variables like: cannula design/size/position/flow rate and patient physiology on venovenous (VV) extracorporeal membrane oxygenation (ECMO). Here we use a CFD model to determine what role cardiac output (CO) plays and to analyse return cannula dynamics.

METHODS

Using a patient-averaged model of the right atrium and venae cava, we virtually inserted a 19Fr return cannula and a 25Fr drainage cannula. Running large eddy simulations, we assessed cardiac output at: 3.5-6.5 L/min and ECMO flow rate at: 2-6 L/min. We analysed recirculation fraction (Rf), time-averaged wall shear stress (TAWSS), pressure, velocity, and turbulent kinetic energy (TKE) and extracorporeal flow fraction (EFF = ECMO flow rate/CO).

RESULTS

Increased ECMO flow rate and decreased CO (high EFF) led to increased Rf (R = 0.98, log fit). Negative pressures developed in the venae cavae at low CO and high ECMO flow (high CR). Mean return cannula TAWSS was >10 Pa for all ECMO flow rates, with majority of the flow exiting the tip (94.0-95.8 %).

CONCLUSIONS

Our results underpin the strong impact of CO on VV ECMO. A simple metric like EFF, once supported by clinical data, might help predict Rf for a patient at a given ECMO flow rate. The return cannula imparts high shear stresses on the blood, largely a result of the internal diameter.

Evaluation of a new magnetically suspended centrifugal neonatal pump in healthy animals using a veno-venous extracorporeal membrane oxygenation configuration

BACKGROUND

The objective of this animal study was to evaluate the hemodynamic performance of a new centrifugal pump for extra-corporeal membrane oxygenation (ECMO) support in neonates.

METHODS

Six healthy swines were supported with veno-venous ECMO with the New Born ECMOLife centrifugal pump (Eurosets, Medolla, Italy) at different flow rates: 0.25, 0.5, 0.6, and 0.8 L/min; three animals were evaluated at low-flows (0.25 and 0.5 L/min) and three at high-flows (0.6 and 0.8 L/min). Each flow was maintained for 4 hours. Blood samples were collected at different time-points. Hematological and biochemical parameters and ECMO parameters [flow, revolutions per minute (RPM), drainage pressure, and the oxygenator pressure drop] were evaluated.

RESULTS

The increase of the pump flow from 0.25 to 0.5 L/min or from 0.6 to 0.8 L/min required significantly higher RPM and produced significantly higher pump pressures [from 0.25 to 0.5 L/min: 1470 (1253-1569) versus 2652 (2589-2750) RPM and 40 (26-57) versus 125 (113-139) mmHg, respectively; p < .0001 for both - from 0.60 to 0.8 L/min: 1950 (1901-2271) versus 2428 (2400-2518) RPM and 66 (62-86) versus 106 (101-113) mmHg, respectively; p < .0001 for both]. Median drainage pressure significantly decreased from -18 (-22; -16) mmHg to -55 (-63; -48) mmHg when the pump flow was increased from 0.25 to 0.5 L/min (p < .0001). When pump flow increased from 0.6 to 0.8 L/min, drainage pressure decreased from -32 (-39; -24) mmHg to -50 (-52; -43) mmHg, (p < .0001). Compared to pre-ECMO values, the median levels of lactate dehydrogenase, d-dimer, hematocrit, and platelet count decreased after ECMO start at all flow rates, probably due to hemodilution. Plasma-free hemoglobin, instead, showed a modest increase compared to pre-ECMO values during all experiments at different pump flow rates. However, these changes were not clinically relevant.

CONCLUSIONS

In this animal study, the "New Born ECMOLife" centrifugal pump showed good hemodynamic performance. Long-term studies are needed to evaluate biocompatibility of this new ECMO pump.

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.

Hemostatic Complications During Neonatal Extracorporeal Membrane Oxygenation: Roller Pump and Centrifugal Pump Driven Circuits

Recently three different neonatal extracorporeal membrane oxygenation (ECMO) circuits have been employed in our clinic. These circuits were compared for clotting and bleeding complications. Initially, we used an ECMO circuit containing a roller pump and venous bladder without severe complications. Manufacturing of circuit components was discontinued, necessitating the replacement of this circuit by a circuit with a centrifugal pump with 3/8 inch inlet and outlet. Acute increase of oxygenator resistance requiring emergency changeout became unexpectedly a regularly occurring complication. The increase in resistance was suspected to be caused by oxygenator clotting, although oxygenator function was preserved. To prevent this complication, we changed to a levitating centrifugal pump with 1/4 inch inlet and outlet, after which no oxygenator malfunction has been observed. Macroscopic and electron microscopic analysis demonstrates that small clots are formed within the circuit, presumably in or near the centrifugal pump, which are transported to the oxygenator and clog up the hollow fiber layer at the inlet side, barely penetrating the oxygenator beyond this first layer. Our results suggest that low blood velocities accompanied with recirculation of blood within or near the centrifugal pump and/or heat generation within the pump could contribute to the formation of these clots.

Hemodynamic and recirculation performance of dual lumen cannulas for venovenous extracorporeal membrane oxygenation

Venovenous extracorporeal membrane oxygenation (ECMO) can be performed with two single lumen cannulas (SLCs) or one dual-lumen cannula (DLC) where low recirculation fraction ([Formula: see text]) is a key performance criterion. DLCs are widely believed to have lower [Formula: see text], though these have not been directly compared. Similarly, correct positioning is considered critical although its impact is unclear. We aimed to compare two common bi-caval DLC designs and quantify [Formula: see text] in several positions. Two different commercially available DLCs were sectioned, measured, reconstructed, scaled to 27Fr and simulated in our previously published patient-averaged computational model of the right atrium (RA) and venae cavae at 2-6 L/min. One DLC was then used to simulate ± 30° and ± 60° rotation and ± 4 cm insertion depth. Both designs had low [Formula: see text] (< 7%) and similar SVC/IVC drainage fractions and pressure drops. Both cannula reinfusion ports created a high-velocity jet and high shear stresses in the cannula (> 413 Pa) and RA (> 52 Pa) even at low flow rates. Caval pressures were abnormally high (16.2-23.9 mmHg) at low flow rates. Rotation did not significantly impact [Formula: see text]. Short insertion depth increased [Formula: see text] (> 31%) for all flow rates whilst long insertion only increased [Formula: see text] at 6 L/min (24%). Our results show that DLCs have lower [Formula: see text] compared to SLCs at moderate-high flow rates (> 4 L/min), but high shear stresses. Obstruction from DLCs increases caval pressures at low flow rates, a potential reason for increased intracranial hemorrhages. Cannula rotation does not impact [Formula: see text] though correct insertion depth is critical.

The Disconnect Between Extracorporeal Circulation and the Microcirculation: A Review

Extracorporeal circulation (ECC) procedures, such as cardiopulmonary bypass (CPB) and extracorporeal membrane oxygenation (ECMO), take over the function of one or more organs, providing clinicians time to treat underlying pathophysiological conditions. ECMO and CPB carry significant mortality rates for patients, despite prior decades of research focused on the resulting failure of critical organs. Since the focus of these procedures is to support blood flow and provide oxygen-rich blood to tissues, a shift in research toward the effects of ECMO and CPB on the microcirculation is warranted. Along with provoking systemic responses, both procedures disrupt the integrity of red blood cells, causing release of hemoglobin (Hb) from excessive foreign surface contact and mechanical stresses. The effects of hemolysis are especially pronounced in the microcirculation, where plasma Hb leads to nitric oxide scavenging, oxidization, formation of reactive oxygen species, and inflammatory responses. A limited number of studies have investigated the implications of ECMO in the microcirculation, but more work is needed to minimize ECMO-induced reduction of microcirculatory perfusion and consequently oxygenation. The following review presents existing information on the implications of ECMO and CPB on microvascular function and proposes future studies to understand and leverage key mechanisms to improve patient outcomes.

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.

Investigation of the influence of blade configuration on the hemodynamic performance and blood damage of the centrifugal blood pump

PURPOSE

The design and optimization of centrifugal blood pumps is crucial for improved extracorporeal membrane oxygenation system performances. Secondary flow passages are common in centrifugal blood pumps, allowing for a high volume of unstable flow. Traditional design theory offers minimal guidance on the design and optimization of centrifugal blood pumps, so it's critical to understand how design parameter variables affect hydraulic performances and hemocompatibility.

METHODS

Computational fluid dynamics (CFD) was employed to investigate the effects of blade number, blade wrap angle, blade thickness, and splitters on pressure head, hemolysis, and platelet activation state. Eulerian and Lagrangian features were used to analyze the flow fields and hemocompatibility metrics such as scalar shear stress, velocity distribution, and their correlation.

RESULTS

The equalization of frictional and flow losses allow impellers with more blades and smaller wrap angles to have higher pressure heads, whereas the trade-off between the volume of high scalar shear stress and exposure time allows impellers with fewer blades and larger blade wrap angles to have a lower HI; there are configurations that increase the possibility of platelet activation for both number of blades and wrap angles. The hydraulic performance and hemocompatibility of centrifugal blood pumps are not affected by blade thickness. Compared to the main blades, a splitters can improve the blood compatibility of a centrifugal blood pump with a small reduction in pressure head, but there is a trade-off between the length and location of the splitter that suppresses flow losses while reducing the velocity gradient. According to correlation analysis, pressure head, HI, and the volume of high shear stress were all substantially connected, and exposure time had a significant impact on HI. The platelet activation state was influenced by the average scalar shear stress and the volume of low velocity.

CONCLUSION

The findings reveal the impact of design variables on the performance of centrifugal blood pumps with secondary flow passages, as well as the relationship between hemocompatibility, hydraulic performance, and flow characteristics, and are useful for the design and optimization of this type of blood pump, as well as the prediction of clinical complications.