Modeling and prediction of flow-induced hemolysis: a review

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

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

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.

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

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

Analysis of the Suitability of an Effective Viscosity to Represent Interactions Between Red Blood Cells in Shear Flow

Many methods to computationally predict red blood cell damage have been introduced, and among these are Lagrangian methods that track the cells along their pathlines. Such methods typically do not explicitly include cell-cell interactions. Due to the high volume fraction of red blood cells (RBCs) in blood, these interactions could impact cell mechanics and thus the amount of damage caused by the flow. To investigate this question, cell-resolved simulations of red blood cells in shear flow were performed for multiple interacting cells, as well as for single cells in unbounded flow at an effective viscosity. Simulations run without adjusting the bulk viscosity produced larger errors unilaterally and were not considered further for comparison. We show that a periodic box containing at least 8 cells and a spherical harmonic of degree larger than 10 are necessary to produce converged higher-order statistics. The maximum difference between the single-cell and multiple-cell cases in terms of peak strain was 3.7%. To achieve this, one must use the whole blood viscosity and average over multiple cell orientations when adopting a single-cell simulation approach. The differences between the models in terms of average strain were slightly larger (maximum difference of 6.9%). However, given the accuracy of the single-cell approach in predicting the maximum strain, which is useful in hemolysis prediction, and its computational cost that is orders of magnitude less than the multiple-cell approach, one may use it as an affordable cell-resolved approach for hemolysis prediction.

Incorporating Unresolved Stresses in Blood Damage Modeling: Energy Dissipation More Accurate Than Reynolds Stress Formulation

OBJECTIVE

Reynolds Averaged Navier Stokes (RANS) models are often used as the basis for modeling blood damage in turbulent flows. To predict blood damage by turbulence stresses that are not resolved in RANS, a stress formulation that represents the corresponding scales is required. Here, we compare two commonly employed stress formulations: a scalar stress representation that uses Reynolds stresses as a surrogate for unresolved fluid stresses, and an effective stress formulation based on energy dissipation.

METHODS

We conducted unsteady RANS simulations of the CentriMag blood pump with three different closure models and a Large Eddy Simulation (LES) for reference. We implemented both stress representations in all models and compared the resulting total stress distributions in Eulerian and Lagrangian frameworks.

RESULTS

The Reynolds-stress-based approach overestimated the contribution of unresolved stresses in RANS, with differences between closure models of up to several orders of magnitude. With the dissipation-based approach, the total stresses predicted with RANS deviated by about 50% from the LES reference, which was more accurate than only considering resolved stresses.

CONCLUSION

The Reynolds-stress-based formulation proved unreliable for estimating scalar stresses in our RANS simulations, while the dissipation-based approach provided an accuracy improvement over simply neglecting unresolved stresses.

SIGNIFICANCE

Our results suggest that dissipation-based inclusion of unresolved stresses should be the preferred choice for blood damage modeling in RANS.

An in silico analysis of unsteady flow structures in a microaxial blood pump under a pulsating rotation speed

BACKGROUND AND OBJECTIVE

Ventricular assist devices (VADs) are generally designed to perform continuous flow. However, it has been proven that continuous flow, which is not a physiological hemodynamic state, may cause severe complications such as gastrointestinal bleeding, pulmonary hypertension, and ventricular suction. For these reasons, many pulsating blood pump control strategies have been proposed and have the potential for application in percutaneous ventricular assist devices (pVADs) or microaxial blood pumps. A few cases report extra hemolysis when introducing pulsating speed, while none involve blood pumps. This research's primary purpose is to evaluate the potential hemolysis of pVAD under pulsating flow conditions.

METHODS

First, the pulsating flow state is deduced using a heart failure model and varying speed. The heart model is established according to the pathology state collected from a clinical check. The rotation speed and boundary physical state are set to fit the heart failure model. The computational fluid dynamics (CFD) method with the hemolysis prediction model is performed. Furthermore, we used proper orthogonal decomposition (POD) analysis to reconstruct the flow field and obtain more details about shearing and transporting effects.

RESULTS

(1) As a variable rotational speed was introduced, no significant gain in hemolysis accumulation appeared in pVAD. This is quite different from long-term implantable VADs. (2) Pulsation affects hemolysis mainly through pressure (or normal stress). Variable rotational speed affects hemolysis mainly through flow instability. (3) Variable rotational speed will increase the instability and influence hemolysis by transporting and shearing effects, while the transporting effect is more significant.

CONCLUSIONS

The unsteady flow state will affect the spatial distribution of hemolysis, which should be taken into account during control strategy and impeller shape design.

Effect of flow rate ratio and positioning on a lighthouse tip ECMO return cannula

Extracorporeal membrane oxygenation is a life-saving support therapy in the case of cardiopulmonary refractory failure. Its use is associated to complications due to the presence of artificial surfaces and supraphysiological stress conditions. Thus, knowledge of the fluid structures associated to each component can give insight into sources of blood damage. In this study, an experimentally validated numerical study of a conventional lighthouse tip cannula in return configuration was carried out to characterize the flow structures using water or a Newtonian blood analog with different flow rate ratios and cannula positioning and their influence on hemolysis. The results showed that strong shear layers developed where the jets from the side holes met the co-flow. Stationary backflow regions at the vessel wall were also present downstream of the cannula. In the tilted case, the recirculation was much more pronounced on the wide side and almost absent on the narrow side. Small vortical backflow structures developed at the side holes which behaved like obstacles to the co-flow, creating pairs of counter-rotating vortices, which induced locally higher risk of hemolysis. However, global hemolysis index did not show significant deviations. Across the examined flow rate ratios, the holes on the narrow side consistently reinfused a larger fraction of fluid. A radial force developed in the tilted case in a direction so as to recenter the cannula in the vessel.

Development of an efficient, ready to use, blood platelet-release device based on two new flow regime parameters: The periodic hydrodynamic loading and the shear stress accumulation

In vitro production of blood platelets for transfusion purposes is gaining interest. While platelet production is now possible on a laboratory scale, the challenge is to move towards industrial production. Attaining this goal calls for the development of platelet release devices capable of producing large quantities of platelets. To this end, we have developed a continuous-flow platelet release device composed of five spherical chambers each containing two calibrated cones placed in a staggered configuration. Following perfusion of proplatelet-bearing cultured megakaryocytes, the device achieves a high yield of about 100 bona-fide platelets/megakaryocyte, at a flow rate of ∼80 mL/min. Performances and operating conditions comply with the requirements of large-scale platelet production. Moreover, this device enabled an in-depth analysis of the flow regimes through Computational Fluid Dynamics (CFD). This revealed two new universal parameters to be taken into account for an optimal platelet release: i.e. a periodic hydrodynamic load and a sufficient accumulation of shear stress. An efficient 16 Pa.s shear stress accumulation is obtained in our system at a flow rate of 80 mL/min.

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.

Turbulence in surgical suction heads as detected by MRI

BACKGROUND

Blood loss is common during surgical procedures, especially in open cardiac surgery. Allogenic blood transfusion is associated with increased morbidity and mortality. Blood conservation programs in cardiac surgery recommend re-transfusion of shed blood directly or after processing, as this decreases transfusion rates of allogenic blood. But aspiration of blood from the wound area is often associated with increased hemolysis, due to flow induced forces, mainly through development of turbulence.

METHODS

We evaluated magnetic resonance imaging (MRI) as a qualitative tool for detection of turbulence. MRI is sensitive to flow; this study uses velocity-compensated T1-weighted 3D MRI for turbulence detection in four geometrically different cardiotomy suction heads under comparable flow conditions (0-1250 mL/min).

RESULTS

Our standard control suction head Model A showed pronounced signs of turbulence at all flow rates measured, while turbulence was only detectable in our modified Models 1-3 at higher flow rates (Models 1 and 3) or not at all (Model 2).

CONCLUSIONS

The comparison of flow performance of surgical suction heads with different geometries via acceleration-sensitized 3D MRI revealed significant differences in turbulence development between our standard control Model A and the modified alternatives (Models 1-3). As flow conditions during measurement have been comparable, the specific geometry of the respective suction heads must have been the main factor responsible. The underlying mechanisms and causative factors can only be speculated about, but as other investigations have shown, hemolytic activity is positively associated with degree of turbulence. The turbulence data measured in this study correlate with data from other investigations about hemolysis induced by surgical suction heads. The experimental MRI technique used showed added value for further elucidating the underlying physical phenomena causing blood damage due to non-physiological flow.

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.

A microfluidic method to investigate platelet mechanotransduction under extensional strain

Background

Blood platelets have evolved a complex mechanotransduction machinery to rapidly respond to hemodynamic conditions. A variety of microfluidic flow-based approaches have been developed to explore platelet mechanotransduction; however, these experimental models primarily focus on the effects of increased wall shear stress on platelet adhesion events and do not consider the critical effects of extensional strain on platelet activation in free flow.

Objectives

We report the development and application of a hyperbolic microfluidic assay that allows for investigation of platelet mechanotransduction under quasi-homogenous extensional strain rates in the absence of surface adhesions.

Methods

Using a combined computational fluid dynamic and experimental microfluidic approach, we explore 5 extensional strain regimes (geometries) and their effect on platelet calcium signal transduction.

Results

We demonstrate that in the absence of canonical adhesion, receptor engagement platelets are highly sensitive to both initial increase and subsequent decrease in extensional strain rates within the range of 747 to 3319/s. Furthermore, we demonstrate that platelets rapidly respond to the rate of change in extensional strain and define a threshold of ≥7.33 × 10⁶/s/m, with an optimal range of 9.21 × 10⁷ to 1.32 × 10⁸/s/m. In addition, we demonstrate a key role of both the actin-based cytoskeleton and annular microtubules in the modulation of extensional strain-mediated platelet mechanotransduction.

Conclusion

This method opens a window onto a novel platelet signal transduction mechanism and may have potential diagnostic utility in the identification of patients who are prone to thromboembolic complications associated with high-grade arterial stenosis or are on mechanical circulatory support systems, for which the extensional strain rate is a predominant hemodynamic driver.

Addressing inadequate blood flow during normothermic regional perfusion for in-situ donation after circulatory death grafts preservation

Donation after circulatory death (DCD) has emerged as attainable strategy to tackle the issue of organ shortage, expanding the donor pool. The DCD concept has been applied to the multiple declinations of circulatory arrest, as per the Modified Maastricht Classification. Notwithstanding, whichever the scenario, DCD donors experience a variable warm ischemia time whose correlation with graft dysfunction is ascertained. This applies to both "controlled" (cDCD) donors (i.e., the timespan from the withdrawal of life-sustaining therapies to the onset of in-situ perfusion), and "uncontrolled" DCD (uDCD) (i.e., the low-flow period during cardiopulmonary resuscitation - CPR). This sums up to the no-flow time from cardiac arrest to the start of CPR for uDCD donors, and to the no-touch period for both uDCDs and cDCDs. Static and hypothermic storage may not be appropriate for DCD grafts. In order to overcome this ischemic insult, extracorporeal membrane oxygenation devices are adopted to guarantee the in-situ grafts preservation by means of techniques such as the normothermic regional perfusion (NRP) which consists in a selective abdominal perfusion obtained via the endovascular or surgical occlusion of the thoracic aorta. The maintenance of an adequate pump flood throughout NRP is therefore a sine qua non to accomplish the DCD donation. The issue of insufficient pump flow during NRP is prevalent and clinically significant but its management remains technically challenging and not standardized. Hereby we propose a systematic algorithmic approach to address this relevant occurrence.

A novel model for hemolysis estimation in rotating impeller blood pumps considering red blood cell aging

For blood pumps with a rotating vane-structure, hemolysis values are estimated using a stress-based power-law model. It has been reported that this method does not consider the red blood cell (RBC) membrane's shear resistance, leading to inaccurate estimation of the hemolysis value. The focus of this study was to propose a novel hemolysis model which can more accurately predict the hemolysis value when designing the axial flow blood pump. The movement behavior of a single RBC in the shear flow field was simulated at the mesoscale. The critical value of shear stress for physiological injury of RBCs was determined. According to the critical value, the equivalent treatment of RBC aging was studied. A novel hemolysis model was established considering the RBC's aging and the hemolysis' initial value. The model's validity was verified under the experimental conditions of shear stress loading and the conditions of the shear flow field of the blood pump. The results showed that compared with other hemolysis models for estimating the hemolysis value of blood pumps, the novel hemolysis model proposed in this paper could effectively reduce the estimation error of the hemolysis value and provide a reference for the optimal design of rotary vane blood pumps.

Hemolytic-Pred: A machine learning-based predictor for hemolytic proteins using position and composition-based features

Objective

The objective of this study is to propose a novel in-silico method called Hemolytic-Pred for identifying hemolytic proteins based on their sequences, using statistical moment-based features, along with position-relative and frequency-relative information.

Methods

Primary sequences were transformed into feature vectors using statistical and position-relative moment-based features. Varying machine learning algorithms were employed for classification. Computational models were rigorously evaluated using four different validation. The Hemolytic-Pred webserver is available for further analysis at http://ec2-54-160-229-10.compute-1.amazonaws.com/.

Results

XGBoost outperformed the other six classifiers with an accuracy value of 0.99, 0.98, 0.97, and 0.98 for self-consistency test, 10-fold cross-validation, Jackknife test, and independent set test, respectively. The proposed method with the XGBoost classifier is a workable and robust solution for predicting hemolytic proteins efficiently and accurately.

Conclusions

The proposed method of Hemolytic-Pred with XGBoost classifier is a reliable tool for the timely identification of hemolytic cells and diagnosis of various related severe disorders. The application of Hemolytic-Pred can yield profound benefits in the medical field.

Pilot scale production and characterization of next generation high molecular weight and tense quaternary state polymerized human hemoglobin

Polymerized human hemoglobin (PolyhHb) is being studied as a possible red blood cell (RBC) substitute for use in scenarios where blood is not available. While the oxygen (O2 ) carrying capacity of PolyhHb makes it appealing as an O2 therapeutic, the commercial PolyhHb PolyHeme® (Northfield Laboratories Inc.) was never approved for clinical use due to the presence of large quantities of low molecular weight (LMW) polymeric hemoglobin (Hb) species (<500 kDa), which have been shown to elicit vasoconstriction, systemic hypertension, and oxidative tissue injury in vivo. Previous bench-top scale studies in our lab demonstrated the ability to synthesize and purify PolyhHb using a two-stage tangential flow filtration purification process to remove almost all undesirable Hb species (>0.2 µm and <500 kDa) in the material, to create a product that should be safer for transfusion. Therefore, to enable future large animal studies and eventual human clinical trials, PolyhHb synthesis and purification processes need to be scaled up to the pilot scale. Hence in this study, we describe the pilot scale synthesis and purification of PolyhHb. Characterization of pilot scale PolyhHb showed that PolyhHb could be successfully produced to yield biophysical properties conducive for its use as an RBC substitute. Size exclusion high performance liquid chromatography showed that pilot scale PolyhHb yielded a high molecular weight Hb polymer containing a small percentage of LMW Hb species (<500 kDa). Additionally, the auto-oxidation rate of pilot scale PolyhHb was even lower than that of previous generations of PolyhHb. Taken together, these results demonstrate that PolyhHb has the ability to be seamlessly manufactured at the pilot scale to enable future large animal studies and clinical trials.

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.

Contraction of the rigor actomyosin complex drives bulk hemoglobin expulsion from hemolyzing erythrocytes

Erythrocyte ghost formation via hemolysis is a key event in the physiological clearance of senescent red blood cells (RBCs) in the spleen. The turnover rate of millions of RBCs per second necessitates a rapid efflux of hemoglobin (Hb) from RBCs by a not yet identified mechanism. Using high-speed video-microscopy of isolated RBCs, we show that electroporation-induced efflux of cytosolic ATP and other small solutes leads to transient cell shrinkage and echinocytosis, followed by osmotic swelling to the critical hemolytic volume. The onset of hemolysis coincided with a sudden self-propelled cell motion, accompanied by cell contraction and Hb-jet ejection. Our biomechanical model, which relates the Hb-jet-driven cell motion to the cytosolic pressure generation via elastic contraction of the RBC membrane, showed that the contributions of the bilayer and the bilayer-anchored spectrin cytoskeleton to the hemolytic cell motion are negligible. Consistent with the biomechanical analysis, our biochemical experiments, involving extracellular ATP and the myosin inhibitor blebbistatin, identify the low abundant non-muscle myosin 2A (NM2A) as the key contributor to the Hb-jet emission and fast hemolytic cell motion. Thus, our data reveal a rapid myosin-based mechanism of hemolysis, as opposed to a much slower diffusive Hb efflux.

Towards an Innovative Sensor in Smart Capsule for Aerial Drones for Blood and Blood Component Delivery

Aerial drone technology is currently being investigated worldwide for the delivery of blood components. Although it has been demonstrated to be safe, the delivered medical substances still need to be analyzed at the end of the flight mission to assess the level of haemolysis and pH prior to the use in a patient. This process can last up to 30 min and prevent the time saved using drone delivery. Our study aims to integrating an innovative sensor for the haemolysis and pH detection into the Smart Capsule, an already demonstrated technology capable of managing transfusion transport through drones. In the proposed scenario, the haemolysis is evaluated optically by a minilysis device using LED-photodetector combination. The preliminary validation has been demonstrated for both the thermal stability of the Smart Capsule and the haemolysis detection of the minilysis device prototype. Firstly, the onboard temperature test has shown that the delivery system is capable of maintaining proper temperature, even though the samples have been manipulated to reach a higher temperature before inserting into the Smart Capsule. Then, in the laboratory haemolysis test, the trend of linear regression between the outputs from the spectrophotometer and the minilysis prototype confirmed the concept design of the minilysis device.

A cell-resolved, Lagrangian solver for modeling red blood cell dynamics in macroscale flows

When red blood cells (RBCs) experience non-physiologically high stresses, e.g., in medical devices, they can rupture in a process called hemolysis. Directly simulating this process is computationally unaffordable given that the length scales of a medical device are several orders of magnitude larger than that of a RBC. To overcome this separation of scales, the present work introduces an affordable computational framework that accurately resolves the stress and deformation of a RBC in a spatially and temporally varying macroscale flow field such as those found in a typical medical device. The underlying idea of the present framework is to treat RBCs as one-way coupled tracers in the macroscale flow by capturing the effect of the flow on their dynamics but neglecting their effect on the flow at the macroscale. As a result, the RBC dynamics are simulated after those of the flow in a postprocessing step by receiving the fluid velocity gradient tensor measured along the RBC trajectory as the input. To resolve the fluid velocity in the immediate vicinity of the RBC as well as the motion of the membrane, we employ the boundary integral method coupled to a structural solver. The governing equations are discretized in space using spherical harmonics, yielding spectral integration accuracy. The predictions produced by this formulation are in good agreement with those obtained from simulations of spherical capsules in shear flows and optical tweezers experiments. The accuracy of the present method is evaluated using unbounded shear flow as a benchmark. Its computational cost grows proportional to p ⁵, where p is the degree of the spherical harmonic. It also exhibits a fast convergence rate that is approximately O ( p 6 ) for p ⪅ 20.

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

BACKGROUND AND OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Effect of Conical Spiral Flow Channel and Impeller Parameters on Flow Field and Hemolysis Performance of an Axial Magnetic Blood Pump

For a blood pump, the blood flow channel and impeller parameters directly affect the performance of the pump and the resulting blood circulation. The flow channel in particular has a great impact on the hydraulic performance of the pump (e.g., flow and pressure), which directly determines the overall performance of the blood pump. Traditional bearing-supported blood pumps can cause mechanical damage to blood cells, leading to hemolysis and thrombosis. In this study, therefore, we designed a conical spiral axial blood pump with magnetic levitation. The blood pump was supported by electrodynamic bearings in the radial direction and electromagnetic bearings in the axial direction. The impeller and the front and rear hubs were integrated to minimize blood stagnation and reduce the formation of thrombosis. The hub had a conical spiral flow channel design, which not only reduced the size of the impeller but also increased blood flow and pressure while meeting the design requirements. Computational fluid dynamics (CFD) analysis was used to analyze the flow field of the axial blood pump, a power function model was used to establish a hemolysis prediction model, and the particle tracking method was used to obtain the flow trajectories of individual blood cells, thereby predicting hemolysis-related performance of the blood pump. The simulation results showed that the main high shear stress area in the blood pump was located in the impeller inlet and the clearance between the top of the impeller and the inner chamber of the blood pump. When the hub taper angle of the blood pump was 0.72° and the clearance was 0.3 mm, the average hemolysis prediction value was 0.00216. This prediction value was smaller than that of traditional axial blood pumps. These findings can provide an important reference for the structural design of axial blood pumps and for reducing the hemolysis prediction value.

Impact of Altered Vena Cava Flow Rates on Right Atrium Flow Characteristics

The right atrium (RA) combines the superior (SVC) and inferior vena cava (IVC) flows. Treatments like extracorporeal membrane oxygenation (ECMO) and hemodialysis by catheter alter IVC/SVC flows. Here we assess how altered IVC/SVC flow contributions impact RA flow. Four healthy volunteers were imaged with CT, reconstructed and combined into a patient-averaged model. Large Eddy Simulations (LES) were performed for a range of IVC/SVC flow contributions (30-70% each, increments of 5%) and common flow metrics were recorded. Model sensitivity to reconstruction domain extent, constant/pulsatile inlets and hematocrit was also assessed. Consistent with literature, a single vortex occupied the central RA across all flowrates with a smaller counter-rotating vortex, not previously reported, in the auricle. Vena cava flow was highly helical. RA turbulent kinetic energy (TKE) (P=0.027) and time-averaged wall shear stress (WSS) (P<0.001) increased with SVC flow. WSS was lower in the auricle (2 Pa, P<0.001). WSS in the vena cava were equal at IVC/SVC =65/35%. The model was highly sensitive to the reconstruction domain with cropped geometries lacking helicity in the vena cavae, altering RA flow. RA flow was not significantly affected by constant inlets or hematocrit. The rotational flow conventionally described in the RA is confirmed however a new, smaller vortex was also recorded in the auricle. When IVC flow dominates, as is normal, TKE in the RA is reduced and WSS in the vena cavae equalize. Significant helicity exists in the vena cava, a result of distal geometry and this geometry appears crucial to accurately simulating RA flow.

Rational Approach to Chatter in Veno-venous Extracorporeal Membrane Oxygenation to Limit Fluid Administration: An Algorithmic Description

Krishnamurthy V, Karanth S. Rational Approach to Chatter in Venovenous Extracorporeal Membrane Oxygenation to Limit Fluid Administration: An Algorithmic Description. Indian J Crit Care Med 2022;26(2):244–245.

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.

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

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

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

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

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.

Erythrocyte morphological symmetry analysis to detect sublethal trauma in shear flow

The viscoelastic properties of red blood cells (RBC) facilitate flexible shape change in response to extrinsic forces. Their viscoelasticity is intrinsically linked to physical properties of the cytosol, cytoskeleton, and membrane-all of which are highly sensitive to supraphysiological shear exposure. Given the need to minimise blood trauma within artificial organs, we observed RBC in supraphysiological shear through direct visualisation to gain understanding of processes leading to blood damage. Using a custom-built counter-rotating shear generator fit to a microscope, healthy red blood cells (RBC) were directly visualised during exposure to different levels of shear (10-60 Pa). To investigate RBC morphology in shear flow, we developed an image analysis method to quantify (a)symmetry of deforming ellipsoidal cells-following RBC identification and centroid detection, cell radius was determined for each angle around the circumference of the cell, and the resultant bimodal distribution (and thus RBC) was symmetrically compared. While traditional indices of RBC deformability (elongation index) remained unaltered in all shear conditions, following ~100 s of exposure to 60 Pa, the frequency of asymmetrical ellipses and RBC fragments/extracellular vesicles significantly increased. These findings indicate RBC structure is sensitive to shear history, where asymmetrical morphology may indicate sublethal blood damage in real-time shear flow.

Computational Fluid-Structure Interaction Study of a New Wave Membrane Blood Pump

PURPOSE

Wave membrane blood pumps (WMBP) are novel pump designs in which blood is propelled by means of wave propagation by an undulating membrane. In this paper, we computationally studied the performance of a new WMBP design (J-shaped) for different working conditions, in view of potential applications in human patients.

METHODS

Fluid-structure interaction (FSI) simulations were conducted in 3D pump geometries and numerically discretized by means of the extended finite element method (XFEM). A contact model was introduced to capture membrane-wall collisions in the pump head. Mean flow rate and membrane envelope were determined to evaluate hydraulic performance. A preliminary hemocompatibility analysis was performed via calculation of fluid shear stress.

RESULTS

Numerical results, validated against in vitro experimental data, showed that the hydraulic output increases when either the frequency or the amplitude of membrane oscillations were higher, with limited increase in the fluid stresses, suggesting good hemocompatibility properties. Also, we showed better performance in terms of hydraulic power with respect to a previous design of the pump. We finally studied an operating point which achieves physiologic flow rate target at diastolic head pressure of 80 mmHg.

CONCLUSION

A new design of WMBP was computationally studied. The proposed FSI model with contact was employed to predict the new pump hydraulic performance and it could help to properly select an operating point for the upcoming first-in-human trials.

Elongational Stresses and Cells

Fluid forces and their effects on cells have been researched for quite some time, especially in the realm of biology and medicine. Shear forces have been the primary emphasis, often attributed as being the main source of cell deformation/damage in devices like prosthetic heart valves and artificial organs. Less well understood and studied are extensional stresses which are often found in such devices, in bioreactors, and in normal blood circulation. Several microfluidic channels utilizing hyperbolic, abrupt, or tapered constrictions and cross-flow geometries, have been used to isolate the effects of extensional flow. Under such flow cell deformations, erythrocytes, leukocytes, and a variety of other cell types have been examined. Results suggest that extensional stresses cause larger deformation than shear stresses of the same magnitude. This has further implications in assessing cell injury from mechanical forces in artificial organs and bioreactors. The cells' greater sensitivity to extensional stress has found utility in mechanophenotyping devices, which have been successfully used to identify pathologies that affect cell deformability. Further application outside of biology includes disrupting cells for increased food product stability and harvesting macromolecules for biofuel. The effects of extensional stresses on cells remains an area meriting further study.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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.

Heuristic optimization of impeller sidewall gaps-based on the bees algorithm for a centrifugal blood pump by CFD

Optimization studies on blood pumps that require complex designs are gradually increasing in number. The essential design criteria of centrifugal blood pump are minimum shear stress with maximal efficiency. The geometry design of impeller sidewall gaps (blade tip clearance, axial gap, radial gap) is highly effective with regard to these two criteria. Therefore, unlike methods such as trial and error, the optimal dimensions of these gaps should be adjusted via a heuristic method, giving more effective results. In this study, the optimal gaps that can ensure these two design criteria with The Bees Algorithm (BA), which is a population-based heuristic method, are investigated. Firstly, a Computational Fluid Dynamics (CFD) analysis of sample pump models, which are selected according to the orthogonal array and pre-designed with different gaps, are performed. The dimensions of the gaps are optimized through this mathematical model. The simulation results for the improved pump model are nearly identical to those predicted by the BA. The improved pump model, as designed with the optimal gap dimensions so obtained, is able to meet the design criteria better than all existing sample pumps. Thanks to the optimal gap dimensions, it has been observed that compared to average values, it has provided a 42% reduction in aWSS and a 20% increase in efficiency. Moreover, original an approach to the design of impeller sidewall gaps was developed. The results show that computational costs have been significantly reduced by using the BA in blood pump geometry design.

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.

Platelet adhesion potential estimation in a normal and diseased coronary artery model: effects of shear stress magnitude versus shear stress history

Platelets play a salient role in the pathogenesis of coronary diseases; primarily through their adhesion to other platelets, endothelial cells and plasma proteins. It is necessary for platelets to activate in order for them to adhere to these different substrates. One of the key regulatory mechanical factors in platelet activation is shear stress, which has been shown to alter multiple platelet functions through the activation of mechanoreceptors. Our goal was to investigate how different numerical shear stress tracking techniques affect platelet adhesion estimates within physiologically relevant computational models. Previously, we developed a physiological coronary artery computational fluid dynamics model. Shear stress waveforms, obtained from these models, were used to monitor in vitro platelet and endothelial cell adhesion marker expression. In this work, the adhesion marker expression data was regressed to obtain numerical functions for receptor expression predictions. These functions were input into a customized adhesion model utilizing different shear stress tracking techniques. For the normal vascular conditions and minimal pathological disease models, shear stress tracking did not significantly affect the adhesion estimates. However, for the severe pathological model, the two shear stress tracking methods had vastly different estimates. Therefore, shear stress tracking methods must be chosen accurately to predict platelet adhesion potentials for accurate modeling techniques.

Model Verification and Error Sensitivity of Turbulence-Related Tensor Characteristics in Pulsatile Blood Flow Simulations

Model verification, validation, and uncertainty quantification are essential procedures to estimate errors within cardiovascular flow modeling, where acceptable confidence levels are needed for clinical reliability. While more turbulent-like studies are frequently observed within the biofluid community, practical modeling guidelines are scarce. Verification procedures determine the agreement between the conceptual model and its numerical solution by comparing for example, discretization and phase-averaging-related errors of specific output parameters. This computational fluid dynamics (CFD) study presents a comprehensive and practical verification approach for pulsatile turbulent-like blood flow predictions by considering the amplitude and shape of the turbulence-related tensor field using anisotropic invariant mapping. These procedures were demonstrated by investigating the Reynolds stress tensor characteristics in a patient-specific aortic coarctation model, focusing on modeling-related errors associated with the spatiotemporal resolution and phase-averaging sampling size. Findings in this work suggest that attention should also be put on reducing phase-averaging related errors, as these could easily outweigh the errors associated with the spatiotemporal resolution when including too few cardiac cycles. Also, substantially more cycles are likely needed than typically reported for these flow regimes to sufficiently converge the phase-instant tensor characteristics. Here, higher degrees of active fluctuating directions, especially of lower amplitudes, appeared to be the most sensitive turbulence characteristics.

Ex vivo assessment of erythrocyte tolerance to the HeartWare ventricular assist device operated in three discrete configurations

Despite technological advances in ventricular assist devices (VADs) to treat end-stage heart failure, hemocompatibility remains a constant concern, with supraphysiological shear stresses an unavoidable reality with clinical use. Given that impeller rotational speed is related to the instantaneous shear within the pump housing, it is plausible that the modulation of pump speed may regulate peak mechanical shear stresses and thus ameliorate blood damage. The present study investigated the hemocompatibility of the HeartWare HVAD in three configurations typical of clinical applications: standard systemic support left VAD (LVAD), pediatric support LVAD, and pulmonary support right VAD (RVAD) conditions. Two ex vivo mock circulation blood loops were constructed using explanted HVADs, in which pump speed and external loop resistance were manipulated to reflect the flow rates and differential pressures reported in configurations for standard adult LVAD (at 3150 rev⸱min^-1 ), pediatric LVAD (at 2400 rev⸱min^-1 ), and adult RVAD (at 1900 rev⸱min^-1 ). Using bovine blood, the mock circulation blood loops were tested at 37°C over a period of 6 hours (consistent with ASTM F1841-97) and compared with static control. Hemocompatibility assessments were conducted for each test condition, examining hematology, hemolysis (absolute and normalized index), osmotic fragility, and blood viscosity. Regardless of configuration, continuous exposure of blood to the VAD over the 6-hour period significantly altered hematological and rheological blood parameters, and induced increased hemolysis when compared with a static control sample. Comparison of the three operational VAD configurations identified that the adult LVAD condition-associated with the highest pump speed, flow rate, and differential pressure across the pump-resulted in increased normalized hemolysis index (NIH; 0.07) when compared with the lower pump speed "off-label" counterparts (NIH of 0.04 in pediatric LVAD and 0.01 in adult RVAD configurations). After normalizing blood residence times between configurations, pump speed was identified as the primary determinant of accumulated blood damage; plausibly, blood damage could be limited by restricting pump speed to the minimum required to support matched cardiac output, but not beyond.

Sublethal Supraphysiological Shear Stress Alters Erythrocyte Dynamics in Subsequent Low-Shear Flows

Blood is a non-Newtonian, shear-thinning fluid owing to the physical properties and behaviors of red blood cells (RBCs). Under increased shear flow, pre-existing clusters of cells disaggregate, orientate with flow, and deform. These essential processes enhance fluidity of blood, although accumulating evidence suggests that sublethal blood trauma-induced by supraphysiological shear exposure-paradoxically increases the deformability of RBCs when examined under low-shear conditions, despite obvious decrement of cellular deformation at moderate-to-higher shear stresses. Some propose that rather than actual enhancement of cell mechanics, these observations are "pseudoimprovements" and possibly reflect altered flow and/or cell orientation, leading to methodological artifacts, although direct evidence is lacking. This study thus sought to explore RBC mechanical responses in shear flow using purpose-built laser diffractometry in tandem with direct optical visualization to address this problem. Freshly collected RBCs were exposed to a mechanical stimulus known to drastically alter cell deformability (i.e., prior shear exposure (PSE) to 100 Pa × 300 s). Samples were subsequently transferred to a custom-built slit-flow chamber that combined laser diffractometry with direct cell visualization. Cell suspensions were sheared in a stepwise manner (between 0.3 and 5.0 Pa), with each step being maintained for 15 s. Deformability and cell orientation indices were recorded for small-scatter Fraunhofer diffraction patterns and also visualized RBCs. PSE RBCs had significantly decreased visualized and laser-derived deformability at any given shear stress ≥1 Pa. Novel, to our knowledge, observations demonstrated that PSE RBCs had increased heterogeneity of direct visualized orientation with flow vector at any shear, which may be due to greater vorticity and thus instability in 5-Pa flow compared with unsheared control. These findings indicate that shear exposure and stress-strain history can alter subsequent RBC behavior in physiologically relevant low-shear flows. These findings may yield insight into microvascular disorders in recipients of mechanical circulatory support and individuals with hematological diseases that alter physical properties of blood.

Characterization of anisotropic turbulence behavior in pulsatile blood flow

Turbulent-like hemodynamics with prominent cycle-to-cycle flow variations have received increased attention as a potential stimulus for cardiovascular diseases. These turbulent conditions are typically evaluated in a statistical sense from single scalars extracted from ensemble-averaged tensors (such as the Reynolds stress tensor), limiting the amount of information that can be used for physical interpretations and quality assessments of numerical models. In this study, barycentric anisotropy invariant mapping was used to demonstrate an efficient and comprehensive approach to characterize turbulence-related tensor fields in patient-specific cardiovascular flows, obtained from scale-resolving large eddy simulations. These techniques were also used to analyze some common modeling compromises as well as MRI turbulence measurements through an idealized constriction. The proposed method found explicit sites of elevated turbulence anisotropy, including a broad but time-varying spectrum of characteristics over the flow deceleration phase, which was different for both the steady inflow and Reynolds-averaged Navier–Stokes modeling assumptions. Qualitatively, the MRI results showed overall expected post-stenotic turbulence characteristics, however, also with apparent regions of unrealizable or conceivably physically unrealistic conditions, including the highest turbulence intensity ranges. These findings suggest that more detailed studies of MRI-measured turbulence fields are needed, which hopefully can be assisted by more comprehensive evaluation tools such as the once described herein.

CFD Assisted Evaluation of In Vitro Experiments on Bearingless Blood Pumps

OBJECTIVE

This paper presents a comparative study of computational fluid dynamics (CFD) simulations and in vitro hemolysis examinations of a bearingless centrifugal blood pump. The outcomes of the in vitro study are analyzed with the help of CFD hemolysis models.

METHODS

Several pump prototypes were manufactured and tested. Each model was implemented in a CFD framework and simulated with different Eulerian hemolysis models. The outcomes are compared to experimental data. The model achieving the highest correlation is used to explain the in vitro outcomes in detail.

RESULTS

It is shown that a double-stage model achieves the best correlation. The sensitivity of the simulation is considerably lower than that of in vitro tests. The CFD model reveals that most of the cell destruction is caused in the radial gap of the pump. Further critical regions are the bottom volume and the shroud clearance gap. Only 0.5% of the priming volume is subject to overcritical shear stress.

CONCLUSION

Cell compatibility can be improved by increasing the radial gap, lowering the shroud and hub clearance gaps, and increasing the fillet radius of the inlet nozzle. CFD models can be used to examine the cell damage effects and help to further improve the pump design.

SIGNIFICANCE

This paper compares different Eulerian CFD hemolysis models, parameter sets, and equivalent shear stresses to several in vitro hemolysis tests. The sensitivity of the models is compared to that of in vitro studies. It is shown that CFD simulations have their limitations but can help with interpreting the outcomes of in vitro studies.

On the Discretization of the Power-Law Hemolysis Model

Flow-induced hemolysis remains a concern for blood-contacting devices, and computer-based prediction of hemolysis could facilitate faster and more economical refinement of such devices. While evaluation of convergence of velocity fields obtained by computational fluid dynamics (CFD) simulations has become conventional, convergence of hemolysis calculations is also essential. In this paper, convergence of the power-law hemolysis model is compared for simple flows, including pathlines with exponentially increasing and decreasing stress, in gradually expanding and contracting Couette flows, in a sudden radial expansion and in the Food and Drug Administration (FDA) channel. In the exponential cases, convergence along a pathline required from one to tens of thousands of timesteps, depending on the exponent. Greater timesteps were required for rapidly increasing (large exponent) stress and for rapidly decreasing (small exponent) stress. Example pathlines in the Couette flows could be fit with exponential curves, and convergence behavior followed the trends identified from the exponential cases. More complex flows, such as in the radial expansion and the FDA channel, increase the likelihood of encountering problematic pathlines. For the exponential cases, comparison of converged hemolysis values with analytical solutions demonstrated that the error of the converged solution may exceed 10% for both rapidly decreasing and rapidly increasing stress.

Development and accuracy evaluation of a degree of occlusion visualization system for roller pumps used in cardiopulmonary bypass

In roller pumps used for cardiopulmonary bypass (CPB), the degree of blockage within the tube resulting from compression of the tube by the rollers, or the degree of occlusion, is closely related to hemolysis, with both tight occlusive and non-occlusive degrees promoting hemolysis. There are as yet no international standards regarding methods of adjusting occlusiveness, and the amount of mechanical stress exerted upon blood remains unknown. To prevent hemolysis during CPB using roller pumps, there is a need to clarify and quantitatively assess the mechanical stress of the occlusiveness of the roller pump. In this study, we have developed a degree of occlusion quantification system which constructs the flow channel shape within an occluded tube from red optical density images, and we have verified the validity of this system. Utilizing a linear actuator, an acrylic roller and raceway, a solution colored with simulated blood powder, and a 3/8-inch vinyl chloride tube, this system uses a camera to capture red optical density images within an occluded tube and constructs the tube flow channel shape using a formula manipulation system. To verify the accuracy of this system, we compared the thickness of a cross-section of the flow channel constructed with the degree of occlusion quantification system with the thickness of a cross-section of silicone cured under the same occlusion conditions. Our experiments indicated that for areas with a small tube gap, this system can construct highly accurate three-dimensional shapes and obtain quantitative indicators assessing the degree of occlusion.

Coarse-Grained Modeling of Pore Dynamics on the Red Blood Cell Membrane under Large Deformations

Transient pore formation on the membrane of red blood cells (RBCs) under high mechanical tensions is of great importance in many biomedical applications, such as RBC damage (hemolysis) and mechanoporation-based drug delivery. The dynamic process of pore formation, growth, and resealing is hard to visualize in experiments. We developed a mesoscale coarse-grained model to study the characteristics of transient pores on a patch of the lipid bilayer that is strengthened by an elastic meshwork representing the cytoskeleton. Unsteady molecular dynamics was used to study the pore formation and reseal at high strain rates close to the physiological ranges. The critical strain for pore formation, pore characteristics, and cytoskeleton effects were studied. Results show that the presence of the cytoskeleton increases the critical strain of pore formation and confines the pore growth. Moreover, the pore recovery process under negative strain rates (compression) is analyzed. Simulations show that pores can remain open for a long time during the high-speed tank-treading induced stretching and compression process that a patch of the RBC membrane usually experiences under high shear flow. Furthermore, complex loading conditions can affect the pore characteristics and result in denser pores. Finally, the effects of strain rate on pore formation are analyzed. Higher rate stretching of membrane patch can result in a significant increase in the critical areal strain and density of pores. Such a model reveals the dynamic molecular process of RBC damage in biomedical devices and mechanoporation that, to our knowledge, has not been reported before.

A novel device for elimination of cancer cells from blood specimens

Circulating tumor cells (CTCs) are derivatives of solid cancerous lesions that detach from the tumor mass and enter the blood circulation. CTCs are considered to be the precursors of metastasis in several cancer types. They are present in the blood of cancer patients as single cells or clusters, with the latter being associated with a higher metastatic potential. Methods to eliminate CTCs from the bloodstream are currently lacking. Here, we took advantage of the lower shear stress-resistance of cancer cells compared to blood cells, and developed a device that can eliminate cancer cells without blood damage. The device consists of an axial pump and a coupled rotating throttle, controllable to prevent local blood flow impairment, yet maintaining a constant shear performance. When processing cancer cells through our device, we observe cancer cell-cluster disruption and viability reduction of single cancer cells, without noticeable effects on human blood cells. When injecting cancer cell-containing samples into tumor-free recipient mice, processed samples fail to generate metastasis. Together, our data show that a selective disruption of cancer cells is possible while preserving blood cells, paving the way towards the development of novel, implantable tools for CTC disruption and metastasis prevention.

Multiscale modeling of hemolysis during microfiltration

In this paper, we propose a multiscale numerical algorithm to simulate the hemolytic release of hemoglobin (Hb) from red blood cells (RBCs) flowing through sieves containing micropores with mean diameters smaller than RBCs. Analyzing the RBC damage in microfiltration is important in the sense that it can quantify the sensitivity of human erythrocytes to mechanical hemolysis while they undergo high shear rate and high deformation. Here, the numerical simulations are carried out via lattice Boltzmann method and spring connected network (SN) coupled by an immersed boundary method. To predict the RBC sublytic damage, a sub-cellular damage model derived from molecular dynamic simulations is incorporated in the cellular solver. In the proposed algorithm, the local RBC strain distribution calculated by the cellular solver is used to obtain the pore radius on the RBC membrane. Index of hemolysis (IH) is calculated by resorting to the resulting pore radius and solving a diffusion equation considering the effects of steric hinderance and increased hydrodynamic drag due to the size of the hemoglobin molecule. It should be mentioned that current computational hemolysis models usually utilize empirical fitting of the released free hemoglobin (Hb) in plasma from damaged RBCs. These empirical correlations contain ad hoc parameters that depend on specific device and operating conditions, thus cannot be used to predict hemolysis under different conditions. In contrast to the available hemolysis model, the proposed algorithm does not have any empirical parameters. Therefore, it can predict the IH in microfilter with different sieve pore sizes and filtration pressures. Also, in contrast to empirical correlations in which the Hb release is related to shear stress and exposure time without considering the physical processes, the proposed model links flow-induced deformation of the RBC membrane to membrane permeabilization and hemoglobin release. In this paper, the cellular solver is validated by simulating optical tweezers experiment, shear flow experiment as well as an experiment to measure RBC deformability in a very narrow microchannel. Moreover, the shape of a single RBC at the rupture moment is compared with corresponding experimental data. Finally, to validate the damage model, the results obtained from our parametric study on the role of filtration pressure and sieve pore size in Hb release are compared with experimental data. Numerical results are in good agreement with experimental data. Similar to the corresponding experiment, the numerical results confirm that hemolysis increases with increasing the filtration pressure and reduction in pore size on the sieve. While in experiment, the RBC pore size cannot be measured, the numerical results can quantify the RBC pore size. The numerical results show that at the sieve pore size of 2.2 μm above 25 cm Hg, RBC pore size is above 75 nm and RBCs experience rupture. More importantly, the results demonstrate that the proposed approach is independent from the operating conditions and it can estimate the hemolysis in a wide range of filtration pressure and sieve pore size with reasonable accuracy.

Prediction of mechanical hemolysis in medical devices via a Lagrangian strain-based multiscale model

This work introduces a new Lagrangian strain-based model to predict the shear-induced hemolysis in biomedical devices. Current computational models for device-induced hemolysis usually utilize empirical fitting of the released free hemoglobin (Hb) in plasma from damaged red blood cells (RBCs). These empirical correlations contain parameters that depend on specific device and operating conditions, thus cannot be used to predict hemolysis in a general device. The proposed algorithm does not have any empirical parameters, thus can presumably be used for hemolysis prediction in various blood-wetting medical devices. In contrast to empirical correlations in which the Hb release is related to the shear stress and exposure time without considering the physical processes, the proposed model links flow-induced deformation of the RBC membrane to membrane permeabilization and Hb release. In this approach, once the steady-state numerical solution of blood flow in the device is obtained under a prescribed operating condition, sample path lines are traced from the inlet of the device to the outlet to calculate the history of the shear stress tensor. In solving the fluid flow, it is assumed that RBCs do not have any influence on the flow pattern. Along each path line, shear stress tensor will be input into a coarse-grained (CG) RBC model to calculate the RBC deformation. Then the correlations obtained from molecular dynamics (MD) simulations are applied to relate the local areal RBC deformation to the perforated area on the RBC membrane. Finally, Hb released out of transient pores is calculated over each path line via a diffusion equation considering the effects of the steric hindrance and increased hydrodynamic drag due to the size of the Hb molecule. The total index of hemolysis (IH) is calculated by integration of released Hb over all the path lines in the computational domain. Hemolysis generated in the Food and Drug Administration (FDA) nozzle and two blood pumps, that is, a CentriMag blood pump (a centrifugal pump) and HeartMate II (an axial pump), for different flow regimes including the laminar and turbulent flows are calculated via the proposed algorithm. In all the simulations, the numerical predicted IH is close to the range of experimental data. The results promisingly indicate that this multiscale approach can be used as a tool for predicting hemolysis and optimizing the hematologic design of other types of blood-wetting devices.

Hemolysis associated with Impella heart pump positioning: In vitro hemolysis testing and computational fluid dynamics modeling

INTRODUCTION

Short-term mechanical circulatory support devices provide temporary hemodynamic support in heart failure and are increasingly used to enable recovery or as a bridge to decision. Blood damage with mechanical circulatory support devices is influenced by many factors, including the magnitude and duration of shear stress and obstruction to blood flow. This study aimed to evaluate the effects of the Impella CP® heart pump positioning on hemolysis using in vitro hemolysis testing and computational fluid dynamics modeling.

METHODS

The in vitro hemolysis testing was conducted per the recommended Food and Drug Administration and American Society for Testing and Materials guidelines. The bench hemolysis testing and computational fluid dynamics simulation analysis were performed for both normal operating (outlet unobstructed) and outlet-obstructed condition of Impella CP (mimicking outlet on the aortic valve due to improper positioning).

RESULTS

The modified index of hemolysis was 2.78 ± 0.69 at normal operating conditions compared to 18.7 ± 7.8 when the Impella CP outlet was obstructed (p = 0.002). Computational fluid dynamics modeling showed about three times increase in exposure time to regions of high shear stress when the Impella CP outlet was obstructed compared to unobstructed condition, thus supporting the experimental observations.

CONCLUSION

Based on these results, it is recommended to ensure proper placement of Impella CP via regular monitoring using echocardiographic guidance or other methods to minimize the risk of hemolysis associated with an obstructed outflow.