Study of flow-induced hemolysis using novel Couette-type blood-shearing devices

The dynamics of red blood cells traversing slits of mechanical heart valves under high shear

Hemolysis, including subclinical hemolysis, is a potentially severe complications of mechanical heart valves (MHVs), which leads to shortened red blood cell (RBC) lifespan and hemolytic anemia. Serious hemolysis is usually associated with structural deterioration and regurgitation. However, the shear stress in MHVs' narrow leakage slits is much lower than the shear stress threshold causing hemolysis and the mechanisms in this context remain largely unclear. This study investigated the hemolysis mechanism of RBCs in cell-size slits under high shear rates by establishing in vitro microfluidic devices and a coarse-grained molecular dynamics (CGMD) model, considering both fluid and structural effects simultaneously. Microfluidic experiments and computational simulation revealed six distinct dynamic states of RBC traversal through MHVs' microscale slits under various shear rates and slit sizes. It elucidated that RBC dynamic states were influenced by not only by fluid forces but significantly by the compressive force of slit walls. The variation of the potential energy of the cell membrane indicated its stretching, deformation, and rupture during traversal, corresponding to the six dynamic states. The maximum forces exerted on membrane by water particles and slit walls directly determined membrane rupture, serving as a critical determinant. This analysis helps in understanding the contribution of the slit walls to membrane rupture and identifying the threshold force that leads to membrane rupture. The hemolysis mechanism of traversing microscale slits is revealed to effectively explain the occurrences of hemolysis and subclinical hemolysis.

Investigation of the Effect of High Shear Stress on Mesenchymal Stem Cells Using a Rotational Rheometer in a Small-Angle Cone-Plate Configuration

Within the healthy human body, cells reside within the physiological environment of a tissue compound. Here, they are subject to constant low levels of mechanical stress that can influence the growth and differentiation of the cells. The liposuction of adipose tissue and the subsequent isolation of mesenchymal stem/stromal cells (MSCs), for example, are procedures that induce a high level of mechanical shear stress. As MSCs play a central role in tissue regeneration by migrating into regenerating areas and driving regeneration through proliferation and tissue-specific differentiation, they are increasingly used in therapeutic applications. Consequently, there is a strong interest in investigating the effects of shear stress on MSCs. In this study, we present a set-up for applying high shear rates to cells based on a rotational rheometer with a small-angle cone-plate configuration. This set-up was used to investigate the effect of various shear stresses on human adipose-derived MSCs in suspension. The results of the study show that the viability of the cells remained unaffected up to 18.38 Pa for an exposure time of 5 min. However, it was observed that intense shear stress damaged the cells, with longer treatment durations increasing the percentage of cell debris.

Computational fluid dynamics simulating of the FDA benchmark blood pump with different coefficient sets and scaler shear stress models used in the power-law hemolysis model

Hemolysis is the most important issue to consider in the design and optimization of blood-contacting devices. Although the use of Computational Fluid Dynamics (CFD) in hemolysis prediction studies provides convenience and has promising potential, it is an extremely challenging process. Hemolysis predictions with CFD depend on the mesh, implementation method, coefficient set, and scalar-shear-stress model. To this end, an attempt was made to find the combination that would provide the most accurate result in hemolysis prediction with the commonly cited power-law based hemolysis model. In the hemolysis predictions conducted using CFD on the Food and Drug Administration (FDA) benchmark blood pump, 3 different scalar-shear-stress models, and 5 different coefficient sets with the power-law based hemolysis model were used. Also, a mesh independence test based on hemolysis and pressure head was performed. The pressure head results of CFD simulations were compared with published pressure head of the FDA benchmark blood pump and a good agreement was observed. In addition, results of CFD-hemolysis predictions which are conducted with scalar-shear-stress model and coefficient set combinations were compared with experimental hemolysis data at three operating conditions such as 6-7 L/min flow rates at 3500 rpm rotational speeds and 6 L/min at 2500 rpm. One of the combinations of the scalar-shear-stress model and the coefficient set was found to be within the error limits of the experimental measurements, while all other combinations overestimated hemolysis.

Assessment of haemolysis models for a positive-displacement total artificial heart

The assessment and reduction of haemolysis within mechanical circulatory support (MCS) remains a concern with regard to device safety and regulatory approval. Numerical methods for predicting haemolysis have typically been applied to rotary MCS devices and the extent to which these methods apply to positive-displacement MCS is unclear. The aim of this study was to evaluate the suitability of these methods for assessing haemolysis in positive-displacement blood pumps. Eulerian scalar-transport and Lagrangian particle-tracking approaches derived from the shear-based power-law relationship were used to calculate haemolysis in a computational fluid dynamics model of the Realheart total artificial heart. A range of power-law constants and their effect on simulated haemolysis were also investigated. Both Eulerian and Lagrangian methods identified the same key mechanism of haemolysis: leakage flow through the bileaflet valves. Whilst the magnitude of haemolysis varied with different power-law constants, the method of haemolysis generation remained consistent. The Eulerian method was more robust and reliable at identifying sites of haemolysis generation, as it was able to capture the persistent leakage flow throughout the entire pumping cycle. This study paves the way for different positive-displacement MCS devices to be compared across different operating conditions, enabling the optimisation of these pumps for improved patient outcomes.

Empirical and Computational Evaluation of Hemolysis in a Microfluidic Extracorporeal Membrane Oxygenator Prototype

Microfluidic devices promise to overcome the limitations of conventional hemodialysis and oxygenation technologies by incorporating novel membranes with ultra-high permeability into portable devices with low blood volume. However, the characteristically small dimensions of these devices contribute to both non-physiologic shear that could damage blood components and laminar flow that inhibits transport. While many studies have been performed to empirically and computationally study hemolysis in medical devices, such as valves and blood pumps, little is known about blood damage in microfluidic devices. In this study, four variants of a representative microfluidic membrane-based oxygenator and two controls (positive and negative) are introduced, and computational models are used to predict hemolysis. The simulations were performed in ANSYS Fluent for nine shear stress-based parameter sets for the power law hemolysis model. We found that three of the nine tested parameters overpredict (5 to 10×) hemolysis compared to empirical experiments. However, three parameter sets demonstrated higher predictive accuracy for hemolysis values in devices characterized by low shear conditions, while another three parameter sets exhibited better performance for devices operating under higher shear conditions. Empirical testing of the devices in a recirculating loop revealed levels of hemolysis significantly lower (<2 ppm) than the hemolysis ranges observed in conventional oxygenators (>10 ppm). Evaluating the model's ability to predict hemolysis across diverse shearing conditions, both through empirical experiments and computational validation, will provide valuable insights for future micro ECMO device development by directly relating geometric and shear stress with hemolysis levels. We propose that, with an informed selection of hemolysis parameters based on the shear ranges of the test device, computational modeling can complement empirical testing in the development of novel high-flow blood-contacting microfluidic devices, allowing for a more efficient iterative design process. Furthermore, the low device-induced hemolysis measured in our study at physiologically relevant flow rates is promising for the future development of microfluidic oxygenators and dialyzers.

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.

Experimental validation of the power law hemolysis model using a Couette shearing device

BACKGROUND

The study of blood trauma, such as hemolysis in blood-carrying devices, is crucial due to the high incidence of adverse events like alteration of blood function, bleeding, and multi-organ failure. The extent of flow-induced hemolysis, predominantly influenced by stress duration and intensity, is described by established model parameters based on the power law approach. In recent years, various parameters were determined using different Couette shearing devices and donor species. However, they have not been validated due to limited experimental data.

METHODS

This study provides hemolysis measurements in a Couette shearing device and evaluates the suitability of different power law parameters. The revised Couette shearing device generates well-defined dynamic stress loads that are repeatedly applied to blood samples at a defined temperature. Human blood samples with an adjusted hematocrit of 30%, were tested with varying repetitions (20 to 80 times). The half-sinusoidal stress loads had amplitudes of 73 to 140 Pa and exposure times of 24 msec per repetition. The parameters of five common power law hemolysis approaches were then compared with the experimental data.

RESULTS

The prediction with the power law model parameters C = 3.458 × 10-6, α = 0.2777 and β = 2.0639 showed a good agreement with the experimental results.

CONCLUSION

The effect of multiple short-time stresses on hemolysis was investigated to validate the power law hemolysis model with the Couette shearing device of this study.

Evaluation of erythrocyte membrane oxidation due to their exposure to shear flow generated by extracorporeal blood pump

Several similarities have been found between shear stress-induced erythrocyte damage and physiological aging of erythrocytes in terms of elevated mechanical fragility, increased erythrocyte aggregation, and decreased membrane surface charge. Accordingly, we hypothesized that blood pump circulation, which generates shear stress, would accelerate erythrocyte aging, manifesting as oxidation. Therefore, the purpose of this study was to investigate the effect of blood pump circulation on erythrocyte oxidation. Fresh porcine blood was acquired from a slaughterhouse and anticoagulated with sodium citrate. About 500 mL of anticoagulated whole blood was circulated for 180 min in an in vitro test circuit comprising a BP-80 blood pump with a pump speed and a pump pressure head of 100-120 mmHg. A blood sample was taken at the start of the circulation and 180 min afterward. The hemolysis level and oxidation amount of the erythrocyte membrane were analyzed and compared between samples. Hemolysis increased with the prolongation of shear exposure inside the pump circuit. After 180 min of blood pumping in circuit, the oxidation level of the erythrocyte membrane showed an increase of 0.1 nmol/mg protein. Moreover, the membrane oxidation levels of sheared erythrocytes were greater than those of control erythrocytes. These results suggest that blood pump circulation accelerates erythrocyte aging and give us a greater understanding of the effects of blood pump perfusion.

Fluid Flow Analysis of Neonatal Dual-Lumen Cannulas for Venovenous Extracorporeal Membrane Oxygenation

Hemolysis persists as a common and serious problem for neonatal patients on extracorporeal membrane oxygenation (ECMO). Since the cannula within the ECMO circuit is associated with hemolysis-inducing shear stresses, real-world internal fluid flow measurements are urgently needed to understand the mechanism and confirm computational estimates. This study appears to be the first experimental study of fluid flow inside commercial ECMO dual-lumen cannulas (DLCs) and first particle image velocimetry (PIV) visualization inside a complicated medical device. The internal geometries of four different opaque neonatal DLCs, both atrial and bicaval positioning geometries each sized 13 Fr and 16 Fr, were replicated by three-dimensional printing clear lumen scaled-up models, which were integrated in a circuit with appropriate ECMO flow parameters. PIV was then used to visualize two-dimensional fluid flow in a single cross section within the models. An empirical model accounting for shear stress and exposure time was used to compare the maximum expected level of hemolysis through each model. The maximum measured peak shear stress recorded was 16±2 Pa in the top arterial bicaval 13 Fr model. The atrial and 16 Fr cannula models never produced greater single-pass peak shear stress or hemolysis than the bicaval and 13 Fr models, respectively, and no difference was found in hemolysis at two different flow rates. After 5 days of flow, small DLC-induced hemolysis values for a single pass through each cannula were modeled to linearly accumulate and caused the most severe hemolysis in the bicaval 13 Fr DLC. Engineering and clinical solutions to improve cannula safety are proposed.

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.

Practical implications of the erroneous treatment of exposure time in the Eulerian hemolysis power law model

BACKGROUND

Eulerian and Lagrangian power-law formulations are both widely used for computational fluid dynamics (CFD) to predict flow-induced hemolysis in blood-contacting medical devices. Both are based on the same empirical power-law correlation between hemolysis and the shear stress and exposure time. In the Lagrangian approach, blood damage is predicted by tracking both the stress and exposure time along a finite number of pathlines in the domain. In the Eulerian approach, a scalar transport equation is solved for a time-linearized damage index within the entire domain. Previous analytical work has demonstrated that there is a fundamental problem with the treatment of exposure time in the Eulerian model formulation such that the only condition under which the model correctly represents the true exposure time is in a flow field with no streamwise velocity variation. However, the practical implications of this limitation have yet to be thoroughly investigated.

METHODS

In this study, we demonstrate the inaccuracy of Eulerian hemolysis power-law model predictions due to the erroneous treatment of exposure time by systematically considering four benchmark test cases with increasing degrees of flow acceleration: Poiseuille flow through a straight tube, inclined Couette flow, and flow through a converging tube with two different convergence ratios.

RESULTS

Compared with Lagrangian predictions, we show that, as flow acceleration becomes more pronounced, the resultant inaccuracy in the Eulerian predictions increases. For the inclined Couette flow case, there is a small degree of flow acceleration that yields a discrepancy in the range of 10% between Lagrangian and Eulerian predictions. For flows with a larger degree of acceleration, as occurs in the converging tube flow cases, the discrepancy is considerably larger (up to 257%).

CONCLUSION

The inaccuracy of hemolysis predictions due to the erroneous treatment of exposure time in the Eulerian power-law model can be significant when there is streamwise velocity variation in the flow field. These results may partially explain the extremely large variability in CFD hemolysis predictions reported in the literature between Lagrangian and Eulerian models.

Cell-scale hemolysis evaluation of intervenient ventricular assist device based on dissipative particle dynamics

Most of the existing hemolysis mechanism studies are carried out on the macro flow scale. They assume that the erythrocyte membranes with different loads will suffer the same damage, which obviously has limitations. Thus, exploring the hemolysis mechanism through the macroscopic flow field information is a tough challenge. In order to further understand the non-physiological shear hemolysis phenomenon at the cell scale, this study used the coarse-grained erythrocytes damage model at the mesoscopic scale based on the transport dissipative particle dynamics (tDPD) method. Combined with computational fluid dynamics the hemolysis of scalarized shear stress (τ) in the clearance of "Impella 5.0" was evaluated under the Lagrange perspective and Euler perspective. The results from the Lagrange perspective showed that the change rate of scaled shear stress (τ˙) was the most critical factor in damaging RBCs in the rotor region of "Impella 5.0"and other transvalvular micro-axial blood pumps. Then, we propose a dimensionless number Dk with time integration based on τ˙ to evaluate hemolysis. The Dissipative particle dynamics simulation results are consistent with the Dk evaluation results, so τ˙ may be an important factor in the hemolysis of VADs. Finally, we tested the hemolysis of 30% hematocrit whole blood in the "Impella 5.0" shroud clearance from the Euler perspective. Relevant results indicate that because of the wall effect, the RBCs near the impeller side are more prone to damage, and most of the cytoplasm is also gathered at the rotor side.

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.

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.

Hemolysis prediction in bio-microfluidic applications using resolved CFD-DEM simulations

BACKGROUND AND OBJECTIVE

Hemolysis, namely hemoglobin leakage from red blood cells (RBCs), is one of the major sources of incorrect results in clinical tests, especially when passive microfluidics is involved. This is due to small characteristic dimensions which could cause strong RBCs deformation. Prediction of hemolysis is essential in the design and optimization of lab-on-a-chip devices for cell sorting and plasma separation. The aim of this work is to provide a numerical simulation tool this purpose applicable to real-scale bio-microfluidic devices with affordable computational cost.

METHODS

Blood is modelled as a suspension of biological cells, mainly RBCs, in liquid plasma assumed as a Newtonian, incompressible carrier fluid. Therefore, the physics of cells and carrier fluid is coupled by means of an immersed boundary concept known as resolved CFD-DEM. In this approach, the Navier-Stokes equations are numerically solved through a finite volume method with an additional penalty term to account for the presence of RBCs. RBCs' positions and velocities are updated by solving Newton and Euler equations for conservation of linear and angular momentum. To model the RBCs deformation, a reduced-order model is employed, where each RBC is represented by a clump of overlapping rigid spheres connected by fictional numerical bonds, whose properties are tuned to reproduce the ones of RBCs viscoelastic membrane. This coupled approach allows access to cell-level information and facilitates the usage of strain-based hemolysis models.

RESULTS

Different micro-channel geometries and blood hematocrits are simulated, to explore the influence of these factors on RBCs damage. Statistical analysis is performed to extract relevant biophysical quantities from numerical simulations such as hemolysis index distribution at the channel exit. Finally, the effect of carrier fluid viscosity is studied in relation to cell-cell interactions.

CONCLUSIONS

Simulation results show that hemolysis occurrence is almost independent of the hematocrit values in the microchannel, implying the possibility to speed up calculation using low hematocrit values. Nevertheless, using whole blood viscosity for the carrier fluid overestimates the value of the hemolysis index by almost one order of magnitude.

Red blood cell sublethal damage: haemocompatibility is not the absence of haemolysis

Blood is a complex fluid owing to its two-phase suspension of formed cellular elements within a protein-rich plasma. Vital to its role in distributing nutrients throughout the circulatory system, the mechanical properties of blood - and particularly red blood cells (RBC)-primarily determine bulk flow characteristics and microcirculatory flux. Various factors impair the physical properties of RBC, including cellular senescence, many diseases, and exposure to mechanical forces. Indeed, the latter is increasingly relevant following the advent of modern life support, such as mechanical circulatory support (MCS), which induce unique interactions between blood and artificial environments that leave blood cells with the signature of aging, albeit accelerated, and crucially underlie various serious complications, including death. Accumulating evidence indicates that these complications appear to be associated with mechanical shear forces present within MCS that are not extreme enough to overtly rupture cells, yet may still induce "sublethal" injury and "fatigue" to vital blood constituents. Impaired RBC physical properties following elevated shear exposure-a hallmark of sublethal injury to blood-are notable and may explain, at least in part, systemic complications and premature mortality associated with MCS. Design of optimal next-generation MCS devices thus requires consideration of biocompatibility and blood-device interactions to minimize potential blood complications and promote clinical success. Presented herein is a contemporary understanding of "blood damage," with emphasis on shear exposures that alter microrheological function but do not overtly destroy cells (ie, sublethal damage). Identification of key cellular factors perturbed by supraphysiological shear exposure are examined, offering potential pathways to enhance design of MCS and blood-contacting medical devices.

Influence of Inlet Boundary Conditions on the Prediction of Flow Field and Hemolysis in Blood Pumps Using Large-Eddy Simulation

Inlet boundary conditions (BC) are one of the uncertainties which may influence the prediction of flow field and hemolysis in blood pumps. This study investigated the influence of inlet BC, including the length of inlet pipe, type of inlet BC (mass flow rate or experimental velocity profile) and turbulent intensity (no perturbation, 5%, 10%, 20%) on the prediction of flow field and hemolysis of a benchmark centrifugal blood pump (the FDA blood pump) and a commercial axial blood pump (Heartmate II), using large-eddy simulation. The results show that the influence of boundary conditions on integral pump performance metrics, including pressure head and hemolysis, is negligible. The influence on local flow structures, such as velocity distributions, mainly existed in the inlet. For the centrifugal FDA blood pump, the influence of type of inlet BC and inlet position on velocity distributions can also be observed at the diffuser. Overall, the effects of position of inlet and type of inlet BC need to be considered if local flow structures are the focus, while the influence of turbulent intensity is negligible and need not be accounted for during numerical simulations of blood pumps.

Intensification of Middle- and High-Molecular-Weight Toxins Removal in Dialysis Process

A novel extracorporeal circuit for hemodialysis is described, which enhances the removal of middle- and high-molecular-weight toxins. To enhance the removal of especially high-molecular toxins, a recirculation pathway for dialysate flow is added to conventional circuit. The influence of an increase in the ratio of recirculation to dialysate flow rate and the ultrafiltration rate (UFR) on the removal of toxins was investigated to evaluate the removal of different toxins in clinical conditions. Removal of toxins was also modeled by an analytical method and solved by the MATLAB software (The MathWorks, Inc., Natick, MA). A significant increase in removal of urea (up to 31%) and vitamin B12 (11%) was achieved when the UFR is low (≤50 ml/h) or zero. The model showed an excellent agreement with the experimental results, which indicates its applicability for the removal of different toxins in an extracorporeal circuit. Increase in recirculation flow, while adjusting the UFR near zero, improves the mass transfer coefficient and can lead to enhanced especially middle- and high-molecular-weight toxin removal.

Analysis of shear stress related hemolysis in a ventricular assist device

BACKGROUND

Implantable devices such as ventricular assist devices provide appropriate treatment for patients with advanced heart failure. Unfortunately, these devices still have many problems, particularly related to blood damage.

OBJECTIVE

The aim of this research is to examine two new ventricular assist devices in terms of induced shear stress, exposure time, and induced hemolysis.

METHOD

Reverse engineering was used on multiple axial flow ventricular assist devices to collect all the details related to the designs (diameters, lengths, blade angles…), which were used to build two prototypes: Model A and Model B.

RESULTS

The obtained results were close to a large extent, except for static pressure rise, where the difference was clear.

CONCLUSION

Compared with what has been published in other studies, the overall performance of both models was excellent.

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.

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.

Flow simulation-based particle swarm optimization for developing improved hemolysis models

The improvement and development of blood-contacting devices, such as mechanical circulatory support systems, is a life saving endeavor. These devices must be designed in such a way that they ensure the highest hemocompatibility. Therefore, in-silico trials (flow simulations) offer a quick and cost-effective way to analyze and optimize the hemocompatibility and performance of medical devices. In that regard, the prediction of blood trauma, such as hemolysis, is the key element to ensure the hemocompatibility of a device. But, despite decades of research related to numerical hemolysis models, their accuracy and reliability leaves much to be desired. This study proposes a novel optimization path, which is capable of improving existing models and aid in the development of future hemolysis models. First, flow simulations of three, turbulent blood flow test cases (capillary tube, FDA nozzle, FDA pump) were performed and hemolysis was numerically predicted by the widely-applied stress-based hemolysis models. Afterward, a multiple-objective particles swarm optimization (MOPSO) was performed to tie the physiological stresses of the simulated flow field to the measured hemolysis using an equivalent of over one million numerically determined hemolysis predictions. The results show that our optimization is capable of improving upon existing hemolysis models. However, it also unveils some deficiencies and limits of hemolysis prediction with stress-based models, which will need to be addressed in order to improve its reliability.

Examining the universality of the hemolysis power law model from simulations of the FDA nozzle using calibrated model coefficients

Computational fluid dynamics (CFD) is widely used to predict mechanical hemolysis in medical devices. The most popular hemolysis model is the stress-based power law model that is based on an empirical correlation between hemoglobin release from red blood cells (RBCs) and the magnitude of flow-induced stress and exposure time. Empirical coefficients are traditionally calibrated using data from experiments in simplified Couette-type blood-shearing devices with uniform-shear laminar flow and well-defined exposure times. Use of such idealized coefficients in simulations of real medical devices with complex hemodynamics is thought to be a primary reason for the historical inaccuracy of absolute hemolysis predictions using the power law model. Craven et al. (Biomech Model Mechanobiol 18:1005-1030, 2019) recently developed a CFD-based Kriging surrogate modeling approach for calibrating empirical coefficients in real devices that could potentially be used to more accurately predict absolute hemolysis. In this study, we use the FDA benchmark nozzle to investigate whether utilizing such calibrated coefficients improves the predictive accuracy of the standard Eulerian power law model. We first demonstrate the credibility of our CFD flow simulations by comparing with particle image velocimetry measurements. We then perform hemolysis simulations and compare the results with in vitro experiments. Importantly, the simulations use coefficients calibrated for the flow of a suspension of bovine RBCs through a small capillary tube, which is relatively comparable to the flow of bovine blood through the FDA nozzle. The results show that the CFD predictions of relative hemolysis in the FDA nozzle are reasonably accurate. The absolute predictions are, however, highly inaccurate with modified index of hemolysis values from CFD in error by roughly three orders of magnitude compared with the experiments, despite using calibrated model coefficients from a relatively similar geometry. We rigorously examine the reasons for the inaccuracy that include differences in the flow conditions in the hemolytic regions of each device and the lack of universality of the hemolysis power law model that is entirely empirical. Thus, while the capability to predict relative hemolysis is valuable for product development, further improvements are needed before the power law model can be relied upon to accurately predict the absolute hemolytic potential of a medical device.

Role of thrombin to non-physiological shear stress induced platelet activation and function alternation

OBJECTIVE

Non-physiological shear stress (NPSS) and thrombin have two distinct mechanisms for activating platelets. NPSS in mechanically assisted circulation (MAC) devices can cause platelet dysfunction, e.g., by shedding its key receptors. In addition, patients with heart failure have increased levels of thrombin generation, which may further affect the NPSS-induced platelet dysfunction, resulting in device-associated complications. This study aimed to assess the combined effect of NPSS and thrombin in platelet activation, expression of adhesion receptors on the platelet surface, and alterations of platelet aggregation.

METHODS

Fresh human blood from healthy donors was divided into two groups; one group was treated by adding 0.01 U/mL thrombin, and another group not treated with thrombin served as a control comparison. They were then pumped through a novel blood shearing device which produces similar shear stress conditions to those in the MAC devices. Three levels of NPSS (i.e., 75, 125, and 175 Pa) with a 1.0 s exposure time were selected for the shearing conditions. Expression of platelet activation markers (PAC-1, activated GPIIb/IIIa and CD62P, platelet surface P-selectin) were investigated along with the shedding of platelet receptors (GPIb, GPIIb/IIIa, and GPVI), generation of platelet microparticles, and Phosphatidylserine (PS)-positive platelets detected by flow cytometry. Platelet aggregation (induced by collagen/ristocetin) was measured by Lumi-aggregometry.

RESULTS

Platelet receptors were shed after exposure to NPSS showing a positive correlation with the level of shear stress. The generation of platelet microparticles and PS-positive platelets also increased with greater NPSS. Elevated NPSS decreased the platelet aggregation capacity. Platelet activation level increased with greater NPSS. Being treated by thrombin can further exacerbate these characteristics under same level of NPSS, except that platelet activation level drastically dropped after the exposure to 175 Pa NPSS in the thrombin-treated blood.

CONCLUSION

After being treated by thrombin, platelets became more susceptible to NPSS, resulting in more receptor shedding, platelet microparticles, and PS-positive platelets, thus limiting platelet aggregation capacity after exposure to NPSS. Platelet activation, in terms of PAC-1 and P-selectin, is an interim status competing between the expression and shedding of these makers/receptors. When platelets have reached a saturation level of activation, exposure to excessive NPSS can potentially impair activation.

Small-animal blood exchange is an emerging approach for systemic aging research

We describe a small-animal blood exchange approach developed for aging research as an alternative to heterochronic parabiosis or plasma injections. In parabiosis, animals are surgically coupled, which has several disadvantages, including difficulty controlling experimental procedure, the effects of shared organs, environmental enrichment from jointly exploring the housing enclosure, involuntary exercise and an imprecise onset of blood sharing. Likewise, in plasma injections, the added volumes need to be small, and there is little flexibility in changing the relative contributions of ectopic to endogenous blood components. These factors complicate the conclusions and interpretations, including the identification of key mechanisms and molecular or cellular determinants. Our approach, where blood is exchanged between animals without them being surgically coupled, is less invasive than parabiosis. The percentage of exchanged blood or other exchanged fluids is known and precise. The age of plasma and cells can be mixed and matched at all desired relative contributions to the endogenous systemic milieu, and the onset of the effects can be accurately delineated. In this protocol, we describe the preparatory and animal surgery steps required for small-animal blood exchange in mice and compare this process with parabiosis and plasma injections. We also provide the design, hardware and software for the blood exchange device and compare automated and manual exchange methods. Lastly, we report mathematical modeling of the dilution of blood factors. The fluid exchange takes ~30 min when performed by a well-trained biomedical scientist; the entire process takes ~2 h.

Targeting biophysical cues to address platelet storage lesions

Platelets play vital roles in vascular repair, especially in primary hemostasis, and have been widely used in transfusion to prevent bleeding or manage active bleeding. Recently, platelets have been used in tissue repair (e.g., bone, skin, and dental alveolar tissue) and cell engineering as drug delivery carriers. However, the biomedical applications of platelets have been associated with platelet storage lesions (PSLs), resulting in poor clinical outcomes with reduced recovery, survival, and hemostatic function after transfusion. Accumulating evidence has shown that biophysical cues play important roles in platelet lesions, such as granule secretion caused by shear stress, adhesion affected by substrate stiffness, and apoptosis caused by low temperature. This review summarizes four major biophysical cues (i.e., shear stress, substrate stiffness, hydrostatic pressure, and thermal microenvironment) involved in the platelet preparation and storage processes, and discusses how they may synergistically induce PSLs such as platelet shape change, activation, apoptosis and clearance. We also review emerging methods for studying these biophysical cues in vitro and existing strategies targeting biophysical cues for mitigating PSLs. We conclude with a perspective on the future direction of biophysics-based strategies for inhibiting PSLs. STATEMENT OF SIGNIFICANCE: Platelet storage lesions (PSLs) involve a series of structural and functional changes. It has long been accepted that PSLs are initiated by biochemical cues. Our manuscript is the first to propose four major biophysical cues (shear stress, substrate stiffness, hydrostatic pressure, and thermal microenvironment) that platelets experience in each operation step during platelet preparation and storage processes in vitro, which may synergistically contribute to PSLs. We first clarify these biophysical cues and how they induce PSLs. Strategies targeting each biophysical cue to improve PSLs are also summarized. Our review is designed to draw the attention from a broad range of audience, including clinical doctors, biologists, physical scientists, engineers and materials scientists, and immunologist, who study on platelets physiology and pathology.

Physical Properties of Blood and their Relationship to Clinical Conditions

It has been long known that blood health heavily influences optimal physiological function. Abnormalities affecting the physical properties of blood have been implicated in the pathogenesis of various disorders, although the exact mechanistic links between hemorheology and clinical disease manifestations remain poorly understood. Often overlooked in current medical practice, perhaps due to the promises offered in the molecular and genetic era, the physical properties of blood which remain a valuable and definitive indicator of circulatory health and disease. Bridging this gap, the current manuscript provides an introduction to hemorheology. It reviews the properties that dictate bulk and microcirculatory flow by systematically dissecting the biomechanics that determine the non-Newtonian behavior of blood. Specifically, the impact of hematocrit, the mechanical properties and tendency of red blood cells to aggregate, and various plasma factors on blood viscosity will be examined. Subsequently, the manner in which the physical properties of blood influence hemodynamics in health and disease is discussed. Special attention is given to disorders such as sickle cell disease, emphasizing the clinical impact of severely abnormal blood rheology. This review expands into concepts that are highly topical; the relation between mechanical stress and intracellular homeostasis is examined through a contemporary cell-signaling lens. Indeed, accumulating evidence demonstrates that nitric oxide is not only transported by erythrocytes, but is locally produced by mechanically-sensitive enzymes, which appears to have intracellular and potentially extracellular effects. Finally, given the importance of shear forces in the developing field of mechanical circulatory support, we review the role of blood rheology in temporary and durable mechanical circulatory support devices, an increasingly utilized method of life support. This review thus provides a comprehensive overview for interested trainees, scientists, and clinicians.

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.

The influence of non-conformal grid interfaces on the results of large eddy simulation of centrifugal blood pumps

BACKGROUND

Computational fluid dynamics has been widely used to assist the design and evaluation of blood pumps. Discretization errors associated with computational grid may influence the credibility of numerical simulations. Non-conformal grid interfaces commonly exist in rotary machines, including rotary blood pumps. Should grid size across the interface differ greatly, large errors may occur.

METHODS

This study explored the effects of non-conformal grid interface on the prediction of the flow field and hemolysis in blood pumps using large eddy simulation (LES). Two benchmarks, a nozzle model and a centrifugal blood pump were chosen as test cases.

RESULTS

This study found that non-conformal grid interfaces with considerable change of grid sizes led to discontinuities of flow variables and brought errors to metrics such as pressure head (7%) and hemolysis (up to 14%).

CONCLUSIONS

The results on the full unstructured grid are more accurate with negligible changes of flow variables across the non-conformal grid interface. A full unstructured grid should be employed for centrifugal blood pumps to minimize the influence of non-conformal grid interfaces for LES simulations.

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.

The Effect of Blood Viscosity on Shear-Induced Hemolysis using a Magnetically Levitated Shearing Device

INTRODUCTION

Blood contacting medical devices, including rotary blood pumps, can cause shear-induced blood damage that may lead to adverse effects in patients. Due in part to an inadequate understanding of how cell-scale fluid mechanics impact red blood cell membrane deformation and damage, there is currently not a uniformly accepted engineering model for predicting blood damage caused by complex flow fields within ventricular assist devices (VADs).

METHODS

We empirically investigated hemolysis in a magnetically levitated axial Couette flow device typical of a rotary VAD. The device is able to accurately control the shear rate and exposure time experienced by blood and to minimize the effects of other uncharacterized stresses. Using this device, we explored the effects of both hematocrit and plasma viscosity on shear-induced hemolysis to characterize blood damage based on the viscosity-independent shear rate, rather than on shear stress.

RESULTS

Over a shear rate range of 20,000-80,000 1/s, the Index of Hemolysis (IH) was found to be dependent upon and well-predicted by shear rate alone. IH was independent of hematocrit, bulk viscosity, or the suspension media viscosity, and less correlated to shear stress (MSE=0.46-0.75) than to shear rate (MSE=0.06-0.09).

CONCLUSION

This study recommends that future investigations of shear-induced blood damage report findings with respect to the viscosity-neutral term of shear rate, in addition to the bulk whole blood viscosity measured at an appropriate shear rate relevant to the flow conditions of the device.

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.

Discrete responses of erythrocytes, platelets, and von Willebrand factor to shear

Despite decades of technological advancements in blood-contacting medical devices, complications related to shear flow-induced blood trauma are still frequently observed in clinic. Blood trauma includes haemolysis, platelet activation, and degradation of High Molecular Weight von Willebrand Factor (HMW vWF) multimers, all of which are dependent on the exposure time and magnitude of shear stress. Specifically, accumulating evidence supports that when blood is exposed to shear stresses above a certain threshold, blood trauma ensues; however, it remains unclear how various constituents of blood are affected by discrete shears experimentally. The aim of this study was to expose blood to discrete shear stresses and evaluate blood trauma indices that reflect red cell, platelet, and vWF structure. Citrated human whole blood (n = 6) was collected and its haematocrit was adjusted to 30 ± 2% by adding either phosphate buffered saline (PBS) or polyvinylpyrrolidone (PVP). Viscosity of whole blood was adjusted to 3.0, 12.5, 22.5 and 37.5 mPa·s to yield stresses of 3, 6, 9, 12, 50, 90 and 150 Pa in a custom-developed shearing system. Blood samples were exposed to shear for 0, 300, 600 and 900 s. Haemolysis was measured using spectrophotometry, platelet activation using flow cytometry, and HMW vWF multimer degradation was quantified with gel electrophoresis and immunoblotting. For tolerance to 300, 600 and 900 s of exposure time, the critical threshold of haemolysis was reached after blood was exposed to 90 Pa for 600 s (P < 0.05), platelet activation and HMW vWF multimer degradation were 50 Pa for 600 s and 12 Pa for 300 s respectively (P < 0.05). Our experimental results provide simultaneous comparison of blood trauma indices and thus also the relation between shear duration and magnitude required to induce damage to red cells, platelets, and vWF. Our results also demonstrate that near-physiological shear stress (<12 Pa) is needed in order to completely avoid any form of blood trauma. Therefore, there is an urgent need to design low shear-flow medical devices in order to avoid blood trauma in this blood-contacting medical device field.

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.

The Design and Evaluation of a Portable Extracorporeal Centrifugal Blood Pump

In recent years, blood pumps have become the bridge to heart transplantation for patients with heart failure. Portability and wearability of blood pumps should be considered to ensure patient satisfaction in everyday life. To date, the focus has been on the development of portable and wearable peripheral components, little attention has been paid to the portable and wearable performance of the blood pump itself. This study reported a novel design of a wearable and portable extracorporeal centrifugal blood pump. Based on an in-house centrifugal maglev blood pump, the wearable and portable blood pump was designed with parallel inlet and outlet pipes to improve the wearable performance. A ring cavity was set at the inlet to convert the circumferential velocity of the inlet pipe to an axial velocity. The hydraulic and hemolytic performance of the baseline and portable blood pumps were analyzed and compared. Compared with the baseline pump, the hydrodynamic and hemolytic performance of the portable pump has been maintained without serious degradation. The results of this study will improve the life quality of patients with heart failure, and enhance the clinical benefits of artificial heart.

Potential Factors for Poor Reproducibility of In Vitro Hemolysis Testing

In vitro testing of hemolysis is essential for the validation and development of ventricular assist devices. However, as many factors influence hemolysis, such tests' inter- and intralaboratory reproducibility is poor. In this work, CentriMag blood pumps were used to conduct a hemolysis study according to ASTM F1841 with blood from 23 bovine donors. Complementary blood analysis, including cell count, plasma composition, and viscosity, was performed to identify factors relevant to the variability of hemolysis testing results. Three strategies were tested to improve reproducibility: albumin supplementation, maintaining glucose concentration, and replacement of plasma with plasma-like buffer. Differences in red blood cell stability among donors were responsible for the largest portion of the total variance. Hematocrit varied widely among donors, and its adjustment to a standard value led to the artificial introduction of between-donor differences, especially in viscosity. It seems likely, that a more careful selection of donors with similar characteristics or repeated blood collection from the same donor could improve reproducibility. However, no direct correlations were found between the susceptibility to hemolysis and individual donor or blood characteristics in this study. The addition of albumin and glucose had a negligible effect while washing blood samples with artificial plasma significantly reduced mean hemolysis, although not its variation. The findings contribute to the understanding of variability in hemolysis experiments and give reason to question the common standard operating procedures, such as hemodilution or maintaining glucose concentration. To confirm the factors identified here, additional studies isolating the effects of individual factors are necessary.

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.

On the Optimization of a Centrifugal Maglev Blood Pump Through Design Variations

Centrifugal blood pumps are usually designed with secondary flow paths to avoid flow dead zones and reduce the risk of thrombosis. Due to the secondary flow path, the intensity of secondary flows and turbulence in centrifugal blood pumps is generally very high. Conventional design theory is no longer applicable to centrifugal blood pumps with a secondary flow path. Empirical relationships between design variables and performance metrics generally do not exist for this type of blood pump. To date, little scientific study has been published concerning optimization and experimental validation of centrifugal blood pumps with secondary flow paths. Moreover, current hemolysis models are inadequate in an accurate prediction of hemolysis in turbulence. The purpose of this study is to optimize the hydraulic and hemolytic performance of an inhouse centrifugal maglev blood pump with a secondary flow path through variation of major design variables, with a focus on bringing down intensity of turbulence and secondary flows. Starting from a baseline design, through changing design variables such as blade angles, blade thickness, and position of splitter blades. Turbulent intensities have been greatly reduced, the hydraulic and hemolytic performance of the pump model was considerably improved. Computational fluid dynamics (CFD) combined with hemolysis models were mainly used for the evaluation of pump performance. A hydraulic test was conducted to validate the CFD regarding the hydraulic performance. Collectively, these results shed light on the impact of major design variables on the performance of modern centrifugal blood pumps with a secondary flow path.

Hemolysis estimation in turbulent flow for the FDA critical path initiative centrifugal blood pump

Hemolysis in medical devices and implants has been a primary concern over the past fifty years. Turbulent flow in particular can cause cell trauma and hemolysis in such devices. In this work, the effects of turbulence on red blood cell (RBC) damage are examined by simulating the flow field through a centrifugal blood pump that has been identified as a case study through the critical path initiative of the US Food and Drug Administration (FDA). In this study, a new model was employed to predict hemolysis in the turbulent flow environment in the pump selected for the FDA critical path initiative. The operating conditions for a centrifugal blood pump were specified by the FDA for rotational speeds of 2500 and 3500 rpm. The model is based on the analysis of the smaller eddies within the turbulent flow field, since it is assumed that turbulent flow eddies with sizes comparable to the dimensions of the RBCs lead to cell trauma. The Kolmogorov length scale of the velocity field is used to identify such small eddies. Using model parameters obtained in prior work through comparisons to capillary and jet flow, it is found that hemolysis for the 2500-rpm pump was predicted well, while hemolysis for the 3500-rpm pump was overpredicted. Results indicate refinement of the model and empirical constants with better experimental data is needed.

Neutrophil injury and function alterations induced by high mechanical shear stress with short exposure time

High mechanical shear stresses (HMSS) can cause damage to blood, which manifests as morphologic changes, shortened life span, biochemical alterations, and complete rupture of blood cells and proteins, leading to the alterations of normal blood function. The aim of this study is to determine the state of neutrophil activation and function alterations caused by HMSS with short exposure time relevant to ventricular assist devices. Blood from healthy donors was exposed to three levels of HMSS (75Pa, 125Pa, and 175Pa) for a short exposure time (0.5 s) using our Couette-type blood-shearing device. Neutrophil activation (Mac-1, platelet-neutrophil aggregates) and surface expression levels of two key functional receptors (CD62L and CD162) on neutrophils were evaluated by flow cytometry. Neutrophil phagocytosis and transmigration were also examined with functional assays. Results showed that the expression of Mac-1 on neutrophils and platelet-neutrophil aggregates increased significantly while the level of CD62L expression on neutrophils decreased significantly after the exposure to HMSS. The Mac-1 expression progressively increased while the CD62L expression progressively decreased with the increased level of HMSS. The level of CD162 expression on neutrophils slightly increased after the exposure to HMSS, but the increase was not significant. The phagocytosis assay data revealed that the ability of neutrophils to phagocytose latex beads coated with fluorescently labeled rabbit IgG increased significantly with the increased level of HMSS. The transmigration ability of neutrophils slightly increased after the exposure to HMSS, but did not reach a significant level. In summary, HMSS with a short exposure time of 0.5 seconds could induce neutrophil activation, platelet-neutrophil aggregation, shedding of CD62L receptor, and increased phagocytic ability. However, the exposure to the three levels of HMSS did not cause a significant change in neutrophil transmigration capacity and shedding of CD162 receptor on neutrophils.

Efficient encapsulation of functional proteins into erythrocytes by controlled shear-mediated membrane deformation

Red blood cells (RBCs) are attractive carriers of biomolecular payloads due to their biocompatibility and the ability to shelter their encapsulated cargo. Commonly employed strategies to encapsulate payloads into RBCs, such as hypotonic shock, membrane fusion or electroporation, often suffer from low throughput and unrecoverable membrane impairment. This work describes an investigation of a method to encapsulate protein payloads into RBCs by controlling membrane deformation either transiently or extendedly in a microfluidic channel. Under the optimized conditions, the loading efficiency of enhanced green fluorescent protein into mouse RBCs increased was about 2.5- and 4-fold compared to that with osmotic entrapment using transient and extended deformation, respectively. Significantly, mouse RBCs loaded with human arginase exhibit higher enzymatic activity and membrane integrity compared to their counterparts loaded by osmotic entrapment. These features together with the fact that this shear-mediated encapsulation strategy allows loading with physiological buffers highlight the key advantages of this approach compared to traditional osmotic entrapment.

Mathematical models for shear-induced blood damage based on vortex platform

Non-physiological shear stress in Ventricular Assist Device (VAD) is considered to be an important trigger of blood damage, which has become the biggest shackle for clinical application. The researches on blood damage in literature were limited to qualitative but did not make much quantitative analysis. The purpose of this study was to investigate the quantitative influence of two flow-dependent parameters: shear stress (rotational speed) and exposure time on the shear-induced damage of red blood cells and von Willebrand Factor (vWF). A vortex blood-shearing platform was constructed to conduct in vitro experiments. Free hemoglobin assay and vWF molecular weight analysis were then performed on the sheared blood samples. MATLAB was used for regression fitting of original experimental data. The quantitative correlations between the hemolysis index, the degradation of high molecular weight vWF and the two flow-dependent parameters were found both following the power law model. The mathematic models indicated that the sensitivity of blood damage on red blood cells and vWF to exposure time was both greater than that of shear stress. Besides, the damage of vWF was more serious than that of red blood cells at the same flow condition. The models could be used to predict blood damage in blood-contacting medical devices, especially for the slow even stagnant blood flow regions in VAD, thus may provide useful guidance for VAD development and improvement. It also indicated that the vortex platform can be used to study the law of blood damage for the simple structure and easy operation.

Production of erythrocyte microparticles in a sub-hemolytic environment

Microparticles are produced by various cells due to a number of different stimuli in the circulatory system. Shear stress has been shown to injure red blood cells resulting in hemolysis or non-reversible sub-hemolytic damage. We hypothesized that, in the sub-hemolytic shear range, there exist sufficient mechanical stimuli for red blood cells to respond with production of microparticles. Red blood cells isolated from blood of healthy volunteers were exposed to high shear stress in a microfluidic channel to mimic mechanical trauma similar to that occurring in ventricular assist devices. Utilizing flow cytometry techniques, both an increase of shear rate and exposure time showed higher concentrations of red blood cell microparticles. Controlled shear rate exposure shows that red blood cell microparticle concentration may be indicative of sub-hemolytic damage to red blood cells. In addition, properties of these red blood cell microparticles produced by shear suggest that mechanical trauma may underlie some complications for cardiovascular patients.

Enhanced discrete phase model for multiphase flow simulation of blood flow with high shear stress

Shear stress is often present in the blood flow within blood-contacting devices, which is the leading cause of hemolysis. However, the simulation method for blood flow with shear stress is still not perfect, especially the multiphase flow model and experimental verification. In this regard, this study proposes an enhanced discrete phase model for multiphase flow simulation of blood flow with shear stress. This simulation is based on the discrete phase model (DPM). According to the multiphase flow characteristics of blood, a virtual mass force model and a pressure gradient influence model are added to the calculation of cell particle motion. In the experimental verification, nozzle models were designed to simulate the flow with shear stress, varying the degree of shear stress through different nozzle sizes. The microscopic flow was measured by the Particle Image Velocimetry (PIV) experimental method. The comparison of the turbulence models and the verification of the simulation accuracy were carried out based on the experimental results. The result demonstrates that the simulation effect of the SST k-ω model is better than other standard turbulence models. Accuracy analysis proves that the simulation results are accurate and can capture the movement of cell-level particles in the flow with shear stress. The results of the research are conducive to obtaining accurate and comprehensive analysis results in the equipment development phase.