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

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.

The Impact of Rotor Axial Displacement Variation on Simulation Accuracy of Fully Magnetic Levitation Centrifugal Blood Pump

The rotor axial displacement of the full magnetic levitation blood pump varies with the operating conditions. The effect of rotor axial displacement on simulation results is unclear. This study aimed to evaluate the effect of rotor axial displacement on the predicted blood pump flow field, hydraulic performance, and hemocompatibility through simulation. This study used the CentriMag blood pump as a model, and conducted computational fluid dynamics simulations to assess the impact of rotor displacement. Considering rotor axial displacement leads to opposite results regarding predicted residence time and thrombotic risk compared with not considering rotor axial displacement. Not considering rotor axial displacement leads to deviations in the predicted values, where the effects on the flow field within the blood pump, ratio of secondary flow, and amount of shear stress >150 Pa are significant. The variation in the back clearance of the blood pump caused by the ideal and actual rotor displacements is the main cause of the above phenomena. Given that the rotor axial displacement significantly impacts the simulation accuracy, the effect of rotor axial displacement must be considered in the simulation.

Optimization design of semi-open impeller based on thrombogenicity in a blood pump

BACKGROUND

Blood clots are composed of aggregated fibrin and platelets, and thrombosis is the body's natural response to repairing injured blood vessels or stopping bleeding. However, when this process is activated abnormally, such as in a mechanical blood pump, it can lead to excessive thrombus formation. Therefore, how to avoid or reduce the probability of thrombus formation is an important indicator of the stable operation of a blood pump.

METHODS

In this paper, Lagrangian particle tracking trajectories are simulated to study platelet transport in a blood pump. The design of the thrombus blood pump was optimized using an orthogonal design method based on three factors: inlet angle, outlet angle, and blade number. The effect of blood pump pressure, rotational speed, impeller outlet angle, inlet angle, and number of blades on thrombus formation was analysed using Fluent software. The thrombogenic potential was derived by analyzing the trajectory and flow parameters of platelet particles in the blood pump, as well as the statistical parameters of residence time and stress accumulation thrombus in the platelet pump.

RESULTS

When the impeller inlet angle is 30°, the outlet angle is 20°, and the number of blades is 6, the probability of thrombus formation is minimized in the orthogonal design method, aligning with the requirements for blood pump performance.

CONCLUSIONS

These design parameters serve as a numerical guideline for optimizing the geometry of the semi-open impeller in blood pumps and provide a theoretical foundation for subsequent in vitro experiments.

[Progress in the analysis of hemolysis and coagulation models for interventional micro-axial flow blood pumps]

Interventional micro-axial flow blood pump is widely used as an effective treatment for patients with cardiogenic shock. Hemolysis and coagulation are vital concerns in the clinical application of interventional micro-axial flow pumps. This paper reviewed hemolysis and coagulation models for micro-axial flow blood pumps. Firstly, the structural characteristics of commercial interventional micro-axial flow blood pumps and issues related to clinical applications were introduced. Then the basic mechanisms of hemolysis and coagulation were used to study the factors affecting erythrocyte damage and platelet activation in interventional micro-axial flow blood pumps, focusing on the current models of hemolysis and coagulation on different scales (macroscopic, mesoscopic, and microscopic). Since models at different scales have different perspectives on the study of hemolysis and coagulation, a comprehensive analysis combined with multi-scale models is required to fully consider the influence of complex factors of interventional pumps on hemolysis and coagulation.

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.

Optimizing the design of axial flow pump blades based on fluid characteristics

Non-physiological blood flow conditions in axial blood pumps lead to some complications, including hemolysis, platelet activation, thrombosis, and embolism. The high speed of the axial blood pump destroys large amounts of erythrocytes, thereby causing hemolysis and thrombosis. Thus, this study aims to reduce the vortices and reflux in the flow field by optimizing the axial blood pump. The axial blood pump and arterial flow field were modeled by the finite element method. The blood was assumed to be incompressible, turbulent, and Newtonian. The SST k-ω turbulence model was used. The frozen rotor method was also used to calculate the snapshot of motion. Many vortices and reflux exist in the flow field of the blood pump without optimization. The improved flow field had almost no vortex and reflux, thereby reducing the exposure time of blood. The optimized blood pump had little influence on the pressure field and shear stress field. The optimized blood pump mainly reduced the vortex, reflux, and the risk of thrombosis in the flow field. The flow field characteristics of an axial blood pump were studied, and the results showed the risk of thrombosis and hemolysis in the blood pump. In accordance with the relationship between the blade shape and the flow field, the blade of the blood pump was optimized, reducing the vortex and reflux of the flow field, as well as the risk of thrombosis.

Clot embolization studies and computational framework for embolization in a canonical tube model

Despite recent advances in the development of computational methods of modeling thrombosis, relatively little effort has been made in developing methods of modeling blood clot embolization. Such a model would provide substantially greater understanding of the mechanics of embolization, as in-vitro and in-vivo characterization of embolization is difficult. Here, a method of computationally simulating embolization is developed. Experiments are performed of blood clots formed in a polycarbonate tube, where phosphate-buffered saline is run through the tube at increasing flow rates until the clot embolizes. The experiments revealed embolization can be initiated by leading edge and trailing edge detachment or by non-uniform detachment. Stress-relaxation experiments are also performed to establish values of constitutive parameters for subsequent simulations. The embolization in the tube is reproduced in silico using a multiphase volume-of-fluid approach, where the clot is modeled as viscoelastic. By varying the constitutive parameters at the wall, embolization can be reproduced in-silico at varying flow rates, and a range of constitutive parameters fitting the experiments is reported. Here, the leading edge embolization is simulated at flow rates consistent with the experiments demonstrating excellent agreement in this specific behavior.

A mathematical model for assessing shear induced bleeding risk

PURPOSE

The objective of this study is to develop a bleeding risk model for assessing device-induced bleeding risk in patients supported with blood contact medical devices (BCMDs).

METHODS

The mathematical model for evaluating bleeding risk considers the effects of shear stress on von Willebrand factor (vWF) unfolding, high molecular weight multimers-vWF (HMWM-vWF) degradation, platelet activation and receptor shedding and platelet-vWF binding ability. Functions of the effect of shear stress on the above factors are fitted/employed and solved by the Eulerian transport equation. An axial flow-through Couette device and two clinical VADs which are HeartWare Ventricular Assist Device (HVAD) and HeartMate II (HM II) blood pump were employed to perform the simulation to evaluate platelet receptor shedding (GPIbα and GPIIb/IIIa), loss of HWMW-vWF, platelet-vWF binding ability and bleeding risk for validating the accuracy of our model.

RESULTS

The platelet-vWF binding ability after being subjected to high shear region in the axial flow-through Couette device predicted by our bleeding model was highly consistent with reported experimental data. As indicated by our CFD simulation results in the axial flow-through Couette device, it can find that an increase in shear stress led to a decrease in the adhesion ability of platelets on vWF, while the binding ability of vWF with platelets first increase and then decrease as shear stress elevates gradually beyond a threshold. The factor of exposure time can enhance the effect of shear stress. Additionally, the shear-induced bleeding risk predicted by our model increases with increasing shear stress and exposure time in an axial flow-through Couette device. As indicated by our numerical model, the bleeding risk in HVAD was higher than HMII, which is highly consistent with the meta-analysis based on clinical statistics. Our simulation investigations in these two clinical VADs also found that HVAD caused a higher rate of platelet receptor shedding and lower damage to HWMW-vWF than HeartMate II. The high shear stress generated in the narrow and turbulent regions of both VADs was the underlying cause of device-induced bleeding.

CONCLUSION

In this study, the shear-induced bleeding risk predicted by our bleeding model in axial flow-through Couette device and two clinical VADs is consistent or highly correlated with experimental and clinical findings, which proves the accuracy of our bleeding model. Our bleeding model can be used to aid the development of new BCMDs with improved functional characteristics and biocompatibility, and help to reduce risk of device-induced adverse events in patients.

The impact of rotor configurations on hemodynamic features, hemocompatibility and dynamic balance of the centrifugal blood pump: A numerical study

To investigate the effect of rotor design configuration on hemodynamic features, hemocompatibility and dynamic balance of blood pumps. Computational fluid dynamics was employed to investigate the effects of rotor type (closed impeller, semi-open impeller), clearance height and back vanes on blood pump performance. In particular, the Eulerian hemolysis model based on a power-law function and the Lagrangian thrombus model with integrated stress accumulation and residence time were applied to evaluate the hemocompatibility of the blood pump. This study shows that compared to the closed impeller, the semi-open impeller can improve hemolysis at a slight sacrifice in head pressure, but increase the risk of thrombogenic potential and disrupt rotor dynamic balance. For the semi-open impeller, the pressure head, hemolysis, and axial thrust of the blood pump decrease with increasing front clearance, and the risk of thrombosis increases first and then decreases with increasing front clearance. Variations in back clearance have little effect on pressure head, but larger on back clearance, worsens hemolysis, thrombogenic potential and rotor dynamic balance. The employment of back vanes has little effect on the pressure head. All back vanes configurations have an increased risk of hemolysis in the blood pump but are beneficial for the improvement of the rotor dynamic balance of the blood pump. Reasonable back vanes configuration (higher height, wider width, longer length and more number) decreases the flow separation, increases the velocity of blood in the back clearance, and reduces the risk of blood pooling and thrombosis. It was also found that hemolysis index (HI) was highly negatively correlated with pressure difference between the top and back clearances (r = -.87), and thrombogenic potential was positively correlated with pressure difference between the top and back clearances (r = .71). This study found that rotor type, clearance height, and back vanes significantly affect the hydraulic performance, hemocompatibility and rotor dynamic balance of centrifugal blood pumps through secondary flow. These parameters should be carefully selected when designing and optimizing centrifugal blood pumps for improving the blood pump clinical outcomes.

Impact of volute design features on hemodynamic performance and hemocompatibility of centrifugal blood pumps used in ECMO

BACKGROUND

The centrifugal blood pump volute has a significant impact on its hemodynamic performance hemocompatibility. Previous studies about the effect of volute design features on the performance of blood pumps are relatively few.

METHODS

In the present study, the computational fluid dynamics (CFD) method was utilized to evaluate the impact of volute design factors, including spiral start position, volute tongue radius, inlet height, size, shape and diffuser pipe angle on the hemolysis index and thrombogenic potential of the centrifugal blood pump.

RESULTS

Correlation analysis shows that flow losses affect the hemocompatibility of the blood pump by influencing shear stress and residence time. The closer the spiral start position of the volute, the better the hydraulic performance and hemocompatibility of the blood pump. Too large or too small volute inlet heights can worsen hydraulic performance and hemolysis, and higher volute inlet height can increase the thrombogenic potential. Small volute sizes exacerbate hemolysis and large volute sizes increase the thrombogenic risk, but volute size does not affect hydraulic performance. When the diffuser pipe is tangent to the base circle of the volute, the best hydraulic performance and hemolysis performance of the blood pump is achieved, but the thrombogenic potential is increased. The trapezoid volute has poor hydraulic performance and hemocompatibility. The round volute has the best hydraulic and hemolysis performance, but the thrombogenic potential is higher than that of the rectangle volute.

CONCLUSION

This study found that the hemolysis index shows a significant correlation with spiral start position, volute size, and diffuser pipe angle. Thrombogenic potential exhibits a good correlation with all the studied volute design features. The flow losses affect the hemocompatibility of the blood pump by influencing shear stress and residence time. The finding of this study can be used to guide the optimization of blood pump for improving the hemodynamic performance and hemocompatibility.

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.