A relationship between Reynolds stresses and viscous dissipation: implications to red cell damage

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

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

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.

Deformation of Transcatheter Aortic Valve Prostheses: Implications for Hypoattenuating Leaflet Thickening and Clinical Outcomes

BACKGROUND

Although transcatheter aortic valve replacement (TAVR) therapy continues to grow, there have been concerns about the occurrence of hypoattenuating leaflet thickening (HALT), which may affect prosthesis function or durability. This study aimed to examine prosthesis frame factors and correlate their extent to the frequency of HALT and clinical outcomes.

METHODS

We prospectively examined 565 patients with cardiac computed tomography screening for HALT at 30 days after balloon-expandable SAPIEN3 and self-expanding EVOLUT TAVR. Deformation of the TAVR prostheses, asymmetric prosthesis leaflet expansion, prosthesis sinus volumes, and commissural alignment were analyzed on the postprocedural computed tomography. For descriptive purposes, an index of prosthesis deformation was calculated, with values >1.00 representing relative midsegment underexpansion. A time-to-event model was performed to evaluate the association of HALT with the clinical outcome.

RESULTS

Overall, HALT was present in 21% of SAPIEN3 patients and in 16% of EVOLUT patients at 30 days after TAVR. The occurrence of HALT was directly associated with greater prosthesis frame deformation (P<0.001), worse asymmetry of the leaflets (P<0.001), and smaller TAVR neosinus volumes (P<0.001). These relations were present in both prosthetic types and in all of their size ranges (all P<0.05). In multivariable analyses that include clinical variables previously associated with HALT (eg, anticoagulant therapy), variables of TAVR prosthesis deformation remained predictive of HALT. Although HALT was not associated with changes in prosthetic hemodynamics, its presence was associated with the risk of mortality at 1 year, with respect to greater incidences of all-cause mortality (hazard ratio, 2.98 [95% CI, 1.57-5.63]; P=0.001), cardiac death (hazard ratio, 4.58 [95% CI, 1.81-11.6]; P=0.001), and a composite outcome of all-cause mortality and heart failure hospitalization (hazard ratio, 1.94 [95% CI, 1.14-3.30]; P=0.02) with adjustment for age, sex, and comorbidities.

CONCLUSIONS

Nonuniform expansion of TAVR prostheses resulting in frame deformation, asymmetric leaflet, and smaller neosinus volume is related to occurrence of HALT in patients who undergo TAVR. These data may have implications for both prosthesis valve design and deployment techniques to improve clinical outcomes for these patients.

Comparison of performance of self-expanding and balloon-expandable transcatheter aortic valves

Objective

Methods

Results

Conclusions

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.

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.

Evaluation of Recirculation During Venovenous Extracorporeal Membrane Oxygenation Using Computational Fluid Dynamics Incorporating Fluid-Structure Interaction

Recirculation in venovenous extracorporeal membrane oxygenation (VV ECMO) leads to reduction in gas transfer efficiency. Studies of the factors contributing have been performed using in vivo studies and computational models. The fixed geometry of previous computational models limits the accuracy of results. We have developed a finite element computational fluid dynamics model incorporating fluid-structure interaction (FSI) that incorporates atrial deformation during atrial filling and emptying, with fluid flow solved using large eddy simulation. With this model, we have evaluated an extensive number of factors that could influence recirculation during two-site VV ECMO, and characterized their impact on recirculation, including cannula construction, insertion depth and orientation, VV ECMO configuration, circuit blood flow, and changes in volume, venous return, heart rate, and blood viscosity. Simulations revealed that extracorporeal blood flow relative to cardiac output, ratio of superior vena caval (SVC) to inferior vena caval (IVC) blood flow, position of the SVC cannula relative to the cavo-atrial junction, and orientation of the return cannula relative to the tricuspid valve had major influences (>20%) on recirculation fraction. Factors with a moderate influence on recirculation fraction (5%-20%) include heart rate, return cannula diameter, and direction of extracorporeal flow. Minimal influence on recirculation (<5%) was associated with atrial volume, position of the IVC cannula relative to the cavo-atrial junction, the number of side holes in the return cannula, and blood viscosity.

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.

Equivalent Scalar Stress Formulation Taking into Account Non-Resolved Turbulent Scales

PURPOSE

Cardiovascular engineering includes flows with fluid-dynamical stresses as a parameter of interest. Mechanical stresses are high-risk factors for blood damage and can be assessed by computational fluid dynamics. By now, it is not described how to calculate an adequate scalar stress out of turbulent flow regimes when the whole share of turbulence is not resolved by the simulation method and how this impacts the stress calculation.

METHODS

We conducted direct numerical simulations (DNS) of test cases (a turbulent channel flow and the FDA nozzle) in order to access all scales of flow movement. After validation of both DNS with literature und experimental data using magnetic resonance imaging, the mechanical stress is calculated as a baseline. Afterwards, same flows are calculated using state-of-the-art turbulence models. The stresses are computed for every result using our definition of an equivalent scalar stress, which includes the influence from respective turbulence model, by using the parameter dissipation. Afterwards, the results are compared with the baseline data.

RESULTS

The results show a good agreement regarding the computed stress. Even when no turbulence is resolved by the simulation method, the results agree well with DNS data. When the influence of non-resolved motion is neglected in the stress calculation, it is underpredicted in all cases.

CONCLUSION

With the used scalar stress formulation, it is possible to include information about the turbulence of the flow into the mechanical stress calculation even when the used simulation method does not resolve any turbulence.

Flow Dynamics and Mixing in Extracorporeal Support: A Study of the Return Cannula

Cannulation strategies in medical treatment such as in extracorporeal life support along with the associated cannula position, orientation and design, affects the mixing and the mechanical shear stress appearing in the flow field. This in turn influences platelet activation state and blood cell destruction. In this study, a co-flowing confined jet similar to a return cannula flow configuration found in extracorporeal membrane oxygenation was investigated experimentally. Cannula diameters, flow rate ratios between the jet and the co-flow and cannula position were studied using Particle Image Velocimetry and Planar Laser Induced Fluorescence. The jet was turbulent for all but two cases, in which a transitional regime was observed. The mixing, governed by flow entrainment, shear layer induced vortices and a backflow along the vessel wall, was found to require 9–12 cannula diameters to reach a fully homogeneous mixture. This can be compared to the 22–30 cannula diameters needed to obtain a fully developed flow. Although not significantly affecting mixing characteristics, cannula position altered the development of the flow structures, and hence the shear stress characteristics.

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.

Effect of membrane insertion for tricuspid regurgitation using immersed-boundary lattice Boltzmann method

BACKGROUND

Tricuspid regurgitation is treated by valve repair or replacement. However, these methods have limitations, and alternative treatment methods are therefore required.

OBJECTIVES

In this study, a new method of tricuspid valve treatment using artificial membrane insertion is analyzed. We performed tricuspid valve simulations using an artificial membrane inserted into the right ventricle (RV) or right atrium (RA).

METHODS

We use the lattice Boltzmann method with the immersed boundary condition to model the structural motion of the valve leaflet. The effect of membrane insertion is analyzed in terms of the stress, force, and impulse on the valve leaflet, along with the velocity, pressure, jet volume, and Reynolds stress in the flow field.

RESULTS

While the use of either membrane (RA or RV) leads to improved valve closure relative to the use of no membrane, the RV membrane is more effective than the RA membrane in achieving improved valve closure. In addition, a larger membrane area with a shorter distance between the leaflet and membrane increases membrane efficacy.

CONCLUSION

Our results suggest that membrane insertion can form an effective new method for the treatment of tricuspid regurgitation.

Fetal Transcatheter Trileaflet Heart Valve Hemodynamics: Implications of Scaling on Valve Mechanics and Turbulence

The scarcity of data available on the best approach for pulmonary fetal valve replacement or implantation necessitate an investigation on whether practices using adult transcatheter valves could be translated to fetal applications. The objective of this study is to evaluate the hemodynamic characteristics and the turbulent properties of a fetal sized trileaflet transcatheter pulmonary valve in comparison with an adult balloon-expandable valve in order to assess the possibility of designing valves for fetal applications using dynamic similarity. A 6 mm fetal trileaflet valve and a 26 mm SAPIEN 3 valve were assessed in a pulse duplicator. Particle image velocimetry was performed. Pressure gradient (ΔP), effective orifice area (EOA), regurgitant fractions (RF), pinwheeling indices (PI) and turbulent stresses were evaluated. ΔP was 8.56 ± 0.139 and 7.76 ± 0.083 mmHg with fetal valve and SAPIEN respectively (p < 0.0001); EOA was 0.10 ± 0.0007 and 2.1 ± 0.025 cm2 with fetal valve and SAPIEN respectively (p < 0.0001); RF with the fetal valve was 2.35 ± 1.99% and with SAPIEN 10.92 ± 0.11% (p < 0.0001); PI with fetal valve was 0.404 ± 0.01 and with SAPIEN 0.37 ± 0.07; The flow regime with the fetal valve was turbulent and Reynolds numbers reached about 7000 while those with the SAPIEN reached about 20,000 at peak velocity. Turbulent stresses were significantly higher with fetal valve compared with SAPIEN. Instantaneous viscous shear stresses with fetal valve were 5.8 times higher than those obtained with SAPIEN and Reynolds shear stresses were 2.5 times higher during peak systole. The fetal valve implantation leads to a turbulent flow (specific to this particular type and design of valve) regime unlike what is expected of a small valve with different flow properties compared to adult valves.

Large-Eddy Simulations of Flow in the FDA Benchmark Nozzle Geometry to Predict Hemolysis

PURPOSE

Modeling of hemolysis due to fluid stresses faces significant methodological challenges, particularly in geometries with turbulence or complex flow patterns. It is currently unclear how existing phenomenological blood-damage models based on laminar viscous stresses can be implemented into turbulent computational fluid dynamics simulations. The aim of this work is to generalize the existing laminar models to turbulent flows based on first principles, and validate this generalization with existing experimental data.

METHODS

A novel analytical and numerical framework for the simulation of flow-induced hemolysis based on the intermittency-corrected turbulent viscous shear stress (ICTVSS) is introduced. The proposed large-eddy simulation framework is able to seamlessly transition from laminar to turbulent conditions in a single flow domain by linking laminar shear stresses to dissipation of mechanical energy, accounting for intermittency in turbulent dissipation, and relying on existing power-law hemolysis models. Simulations are run to reproduce previously published hemolysis data with bovine blood in a benchmark geometry. Two sets of experimental data are relied upon to tune power-law parameters and justify that tuning. The first presents hemolysis measurements in a simple laminar flow, and the second is hemolysis in turbulent flow through the FDA benchmark nozzle. Validation is performed by simulation of blood injected into a turbulent jet of phosphate-buffered saline, with modifications made to account for the local concentration of blood.

RESULTS

Hemolysis predictions are found to be very sensitive to power-law parameters in the turbulent case, though a set of parameters is presented that both matches the turbulent data and is well-justified by the laminar data. The model is shown to be able to predict the general behavior of hemolysis in a second turbulent case. Results suggest that wall shear may play a dominant role in most cases.

CONCLUSION

The ICTVSS framework of generalizing laminar power-law models to turbulent flows shows promise, but would benefit from further numerical validation and carefully designed experiments.

Mixed hemodiafiltration reduces erythropoiesis stimulating agents requirement in dialysis patients: a prospective randomized study

BACKGROUND

Improved responsiveness to erythropoiesis stimulating agents (ESAs) in patients on on-line post-dilution hemodiafiltration (Post-HDF) compared with conventional hemodialysis (HD) was reported by some authors but challenged by others. This prospective, cross-over randomized study tested the hypothesis that an alternative infusion modality of HDF, mixed-dilution HDF (Mixed HDF), could further reduce ESAs requirement in dialysis patients compared to the traditional Post-HDF.

METHODS

One-hundred-twenty prevalent patients from 6 Dialysis Centers were randomly assigned to two six-months treatment sequences: A-B and B-A (A, Mixed HDF; B, Post-HDF). Primary outcome was comparative evaluation of ESA (darbepoetin alfa) requirement and ESA resistance. Treatments efficiency, iron and vitamins status, inflammation and nutrition parameters were monitored.

RESULTS

In sequence A, darbepoetin requirement decreased during Mixed HDF from 29.5 to 23.7 µg/month and increased significantly during Post-HDF (32.3 µg/month at 6th month) while, in sequence B, it increased during Post-HDF from 38.2 to 43.7 µg/month and decreased during Mixed HDF (23.9 µg/month at 6th month). Overall, EPO doses at 6 months on Mixed and Post-HDF were 23.8 and 38.4 µg/month, respectively, P < 0.01. A multiple linear model confirmed that Mixed HDF vs Post-HDF reduced significantly ESA requirement and ESA resistance (P < 0.0001), by a mean of 29% (CI 23-35%) in the last three months of the observation periods.

CONCLUSIONS

Mixed HDF decreased darbepoetin-alfa requirement in dialysis patients. This might help preventing the untoward side effects of high ESA doses, besides having a remarkable economic impact. Additional evidence is needed to confirm this potential benefit of Mixed-HDF.

The Implementation of an Adjustable Afterload Module for Ex Situ Heart Perfusion

PURPOSE

Windkessel impedance analysis has proven to be an effective technique for instituting artificial afterload on ex situ hearts. Traditional fixed parameter afterload modules, however, are unable to handle the changing contractile conditions associated with prolonged ex situ heart perfusion. In this paper, an adjustable afterload module is described comprising of three fully adjustable sub-components: a systemic resistor, a proximal resistor and a compliance chamber.

METHODS

Using a centrifugal pump, the systemic resistor and compliance chamber were subjected to testing across their operating ranges, whereby the predictability of resistance and compliance values was evaluated. The components were then assembled, and the full module tested on three separate porcine hearts perfused for 6 h with success defined by the ability to maintain physiological systolic and diastolic aortic pressures across flow rate variability.

RESULTS

For both the systemic resistor and compliance chamber, experimental measurements agreed with their theoretical equivalents, with coefficients of determination of 0.99 and 0.97 for the systemic resistor and compliance chamber, respectively. During ex situ perfusion, overall 95% confidence intervals demonstrate that physiological systolic (95-96.21 mmHg) and diastolic (26.8-28.8 mmHg) pressures were successfully maintained, despite large variability in aortic flow. Left ventricular contractile parameters, were found to be in line with those in previous studies, suggesting the afterload module has no detrimental impact on functional preservation.

CONCLUSIONS

We conclude that due to the demonstrable control of our afterload module, we can maintain physiological aortic pressures in a passive afterload working mode across prolonged perfusion periods, enabling effective perfusion regardless of contractile performance.

Influence of turbulent shear stresses on the numerical blood damage prediction in a ventricular assist device

The blood damage prediction in rotary blood pumps is an important procedure to evaluate the hemocompatibility of such systems. Blood damage is caused by shear stresses to the blood cells and their exposure times. The total impact of an equivalent shear stress can only be taken into account when turbulent stresses are included in the blood damage prediction. The aim of this article was to analyze the influence of the turbulent stresses on the damage prediction in a rotary blood pump's flow. Therefore, the flow in a research blood pump was computed using large eddy simulations. A highly turbulence-resolving setup was used in order to directly resolve most of the computed stresses. The simulations were performed at the design point and an operation point with lower flow rate. Blood damage was predicted using three damage models (volumetric analysis of exceeded stress thresholds, hemolysis transport equation, and hemolysis approximation via volume integral) and two shear stress definitions (with and without turbulent stresses). For both simulations, turbulent stresses are the dominant stresses away from the walls. Here, they act in a range between 9 and 50 Pa. Nonetheless, the mean stresses in the proximity of the walls reach levels, which are one order of magnitude higher. Due to this, the turbulent stresses have a small impact on the results of the hemolysis prediction. Yet, turbulent stresses should be included in the damage prediction, since they belong to the total equivalent stress definition and could impact the damage on proteins or platelets.

In vitro hemodynamic assessment of a novel polymeric transcatheter aortic valve

Transcatheter aortic valve replacement (TAVR) is a life-saving alternative to surgical intervention. However, the identification of features associated with poor outcomes, including residual paravalvular leakage (PVL), leaflet calcification, and subclinical leaflet thrombosis, are cause to be concerned about valve durablilty (Mylotte and Piazza, 2015a, 2015b; Dasi et al., 2017; Makkar et al., 2015; Kheradvar et al., 2015a). The aim of this study is to optimize the potential of a hyaluronan (HA) enhanced polymeric transcatheter aortic valve (HA-TAV) that has promised to reduce blood damage causing-turbulent flow while maintaining durability. HA-enhanced linear low-density polyethylene (LLDPE) leaflets were sutured to novel cobalt chromium stents, size 26 mm balloon expandable stents. Hemodynamic performance was assessed in a left heart simulator under physiological pressure and flow conditions and compared to a 26 mm Medtronic Evolut and 26 mm Edwards SAPIEN 3. High-speed imaging and particle image velocimetry (PIV) were performed. The HA-TAV demonstrated an effective orifice area (EOA) within one standard deviation of the leading valve, SAPIEN 3.The regurgitant fraction (RF) of the HA-TAV (11.23 ± 0.55%) is decreased in comparison the Evolut (15.74 ± 0.73%) and slightly higher than the SAPIEN 3 (10.92 ± 0.11%), which is considered trace regurgitation according to valve standards. A decreased number of higher principal Reynolds shear stresses were shown for the HA-TAV at each cardiac phase. The HA-TAV is directly comparable and in some cases superior to the leading commercially available prosthetic heart valves in in-vitro hemodynamic testing.

Protect, repair, destroy or sacrifice: a role of oxidative stress biology in inter-donor variability of blood storage?

Red blood cells (RBCs) have been historically regarded as a critical model to investigate cellular and oxidant stress biology. First of all, they are constantly exposed to oxidant stress, as their main function is to transport and deliver oxygen to tissues. Second, they are devoid of de novo protein synthesis capacity, which prevents RBCs from replacing irreversibly oxidised proteins with newly synthesised ones. As such, RBCs have evolved to (i) protect themselves from oxidant stress, in order to prevent oxidant damage from reactive species; (ii) repair oxidatively damaged proteins, through mechanisms that involve glutathione and one-carbon metabolism; (iii) destroy irreversibly oxidised proteins through proteasomal or protease-dependent degradation; and (iv) sacrifice membrane portions through mechanism of vesiculation. In this brief review we will summarize these processes and their relevance to RBC redox biology (within the context of blood storage), with a focus on how polymorphisms in RBC antioxidant responses could contribute to explaining the heterogeneity in the progression and severity of the RBC storage lesion that can be observed across the healthy donor population.

Assessing clinical implications and perspectives of the pathophysiological effects of erythrocytes and plasma free hemoglobin in autologous biologics for use in musculoskeletal regenerative medicine therapies. A review

Autologous biologics, defined as platelet-rich plasma (PRP) and bone marrow aspirate concentrate (BMC), are cell-based therapy treatment options in regenerative medicine practices, and have been increasingly used in orthopedics, sports medicine, and spinal disorders. These biological products are produced at point-of-care; thereby, avoiding expensive and cumbersome culturing and expansion techniques.

Numerous commercial PRP and BMC systems are available but reports and knowledge of bio-cellular formulations produced by these systems are limited. This limited information hinders evaluating clinical and research outcomes and thus making conclusions about their biological effectiveness. Some of their important cellular and protein properties have not been characterized, which is critical for understanding the mechanisms of actions involved in tissue regenerative processes. The presence and role of red blood cells (RBCs) in any biologic has not been addressed extensively. Furthermore, some of the pathophysiological effects and phenomena related to RBCs have not been studied. A lack of a complete understanding of all of the biological components and their functional consequences hampers the development of clinical standards for any biological preparation.

This paper aims to review the clinical implications and pathophysiological effects of RBCs in PRP and BMC; emphasizes hemolysis, eryptosis, and the release of macrophage inhibitory factor; and explains several effects on the microenvironment, such as inflammation, oxidative stress, vasoconstriction, and impaired cell metabolism.

•

Different biological formulations optimize disease specific regenerative treatment protocols.

•

Disintegrated RBC's release harmful components to regenerative therapy treatment vials.

•

The effectiveness of MSC injection depends on the quality of the bone marrow aspiration procedure.

•

PRP and BMC should contain minimal to no erythrocytes.

Use of Computational Fluid Dynamics to Analyze Blood Flow, Hemolysis and Sublethal Damage to Red Blood Cells in a Bileaflet Artificial Heart Valve

Artificial heart valves may expose blood to flow conditions that lead to unnaturally high stress and damage to blood cells as well as issues with thrombosis. The purpose of this research was to predict the trauma caused to red blood cells (RBCs), including hemolysis, from the stresses applied to them and their exposure time as determined by analysis of simulation results for blood flow through both a functioning and malfunctioning bileaflet artificial heart valve. The calculations provided the spatial distribution of the Kolmogorov length scales that were used to estimate the spatial and size distributions of the smallest turbulent flow eddies in the flow field. The number and surface area of these eddies in the blood were utilized to predict the amount of hemolysis experienced by RBCs. Results indicated that hemolysis levels are low while suggesting stresses at the leading edge of the leaflet may contribute to subhemolytic damage characterized by shortened circulatory lifetimes and reduced RBC deformability.

Regulation of erythrocyte Na+/K+/2Cl- cotransport by an oxygen-switched kinase cascade

Many erythrocyte processes and pathways, including glycolysis, the pentose phosphate pathway (PPP), KCl cotransport, ATP release, Na⁺/K⁺-ATPase activity, ankyrin-band 3 interactions, and nitric oxide (NO) release, are regulated by changes in O₂ pressure that occur as a red blood cell (RBC) transits between the lungs and tissues. The O₂ dependence of glycolysis, PPP, and ankyrin-band 3 interactions (affecting RBC rheology) are controlled by O₂-dependent competition between deoxyhemoglobin (deoxyHb), but not oxyhemoglobin (oxyHb), and other proteins for band 3. We undertook the present study to determine whether the O₂ dependence of Na⁺/K⁺/2Cl^- cotransport (catalyzed by Na⁺/K⁺/2Cl^- cotransporter 1 [NKCC1]) might similarly originate from competition between deoxyHb and a protein involved in NKCC1 regulation for a common binding site on band 3. Using three transgenic mouse strains having mutated deoxyhemoglobin-binding sites on band 3, we found that docking of deoxyhemoglobin at the N terminus of band 3 displaces the protein with no lysine kinase 1 (WNK1) from its overlapping binding site on band 3. This displacement enabled WNK1 to phosphorylate oxidative stress-responsive kinase 1 (OSR1), which, in turn, phosphorylated and activated NKCC1. Under normal solution conditions, the NKCC1 activation increased RBC volume and thereby induced changes in RBC rheology. Because the deoxyhemoglobin-mediated WNK1 displacement from band 3 in this O₂ regulation pathway may also occur in the regulation of other O₂-regulated ion transporters, we hypothesize that the NKCC1-mediated regulatory mechanism may represent a general pattern of O₂ modulation of ion transporters in erythrocytes.

Sinus Hemodynamics in Representative Stenotic Native Bicuspid and Tricuspid Aortic Valves: An In-Vitro Study

(1) The study’s objective is to assess sinus hemodynamics differences between stenotic native bicuspid aortic valve (BAV) and native tricuspid aortic valve (TrAV) sinuses in order to assess sinus flow shear and vorticity dynamics in these common pathological states of the aortic valve. (2) Representative patient-specific aortic roots with BAV and TrAV were selected, segmented, and 3D printed. The flow dynamics within the sinus were assessed in-vitro using particle image velocimetry in a left heart simulator at physiological pressure and flow conditions. Hemodynamic data calculations, vortex tracking, shear stress probability density functions and sinus washout calculations based on Lagrangian particle tracking were performed. (3) (a) At peak systole, velocity and vorticity in BAV reach 0.67 ± 0.02 m/s and 374 ± 5 s−1 versus 0.49 ± 0.03 m/s and 293 ± 3 s−1 in TrAV; (b) Aortic sinus vortex is slower to form but conserved in BAV sinus; (c) BAV shear stresses exceed those of TrAV (1.05 Pa versus 0.8 Pa); (d) Complete TrAV washout was achieved after 1.5 cycles while it was not for BAV. 4) In conclusion, sinus hemodynamics dependence on the different native aortic valve types and sinus morphologies was clearly highlighted in this study.

Determination of Reynolds Shear Stress Level for Hemolysis

Reynolds shear stress (RSS) has served as a metric for the effect of turbulence on hemolysis. Forstrom (1969) and Sallam and Hwang (1984) determined the RSS threshold for hemolysis to be 50,000 and 4,000 dyne/cm, respectively, using a turbulent jet. Despite the order of magnitude discrepancy, the threshold by Sallam and Hwang has been frequently cited for hemolytic potential in blood pumps. We recreated a Sallam apparatus (SA) to resolve this discrepancy and provide additional data to be used in developing a more accurate hemolysis model. Hemolysis was measured over a large range of Reynolds numbers (Re) (Re = 1,000-80,000). Washed bovine red blood cells (RBCs) were injected into the free jet of phosphate buffered saline, and hemolysis was quantified using a percent hemolysis, Hp = h (100 - hematocrit [HCT])/Hb, where h (mg/dl) is free hemoglobin and Hb (mg/dl) is total hemoglobin. Reynolds shear stress was calculated using two-dimensional laser Doppler velocimetry. Reynolds shear stress of ≥30,000 dyne/cm corresponding to Re of ≥60,000 appeared to cause hemolysis (p < 0.05). This RSS is an order of magnitude greater than the RSS threshold that Sallam and Hwang suggested, and it is similar to Forstrom's RSS threshold. This study resolved a long-standing uncertainty regarding the critical values of RSS for hemolysis and may provide a foundation for a more accurate hemolysis model.

Stented valve dynamic behavior induced by polyester fiber leaflet material in transcatheter aortic valve devices

OBJECTIVE

This study aims at assessing the global dynamic behavior, elastic deformability, closing energy and turbulence of rigid versus deformable stented (RS vs DS) valve systems with deformable and rigid textile materials (DT vs RT) through studying the stent-valve interaction compared to a bioprosthetic material in transcatheter aortic valves (TAV).

METHODS

Three 19 mm stented textile TAV designs (RS-DT, RS-RT and DS-RT) with different stent and leaflet properties were tested and compared with a control bioprosthetic TAV (RS-DB) in a left heart simulator flow loop under physiological pressure and flow. Particle Image Velocimetry and high speed imaging were performed. Pressure gradients (PG), leakage fractions (LF), Pinwheeling indices (PI), closing energy (E) and Reynolds shear stresses (RSS) were calculated.

RESULTS

(a) PGs and LFs were 11.86 ± 0.51 mmHg, 11.70 ± 0.34%; 8.84 ± 0.40 mmHg, 29.80 ± 0.76%; 11.59 ± 0.12 mmHg, 14.23 ± 1.64%; and 7.05 ± 0.09 mmHg, 12.08 ± 0.45% % for RS-DB, RS-DT, RS-RT and DS-RT respectively. (b) PIs were 15.79 ± 2.34%, 4.36 ± 0.84%, 2.47 ± 0.51% and 2.03 ± 0.33% for RS-DB, RS-DT, RS-RT and DS-RT respectively. (c) E is lowest for DS-RT (0.0010 ± 0.0002 J) followed by RS-RT (0.0017 ± 0.0002 J), RS-DB (0.0023 ± 0.0004 J) and highest with RS-DT (0.0036 ± 0.0007 J). (d) At peak systole lowest RSS was obtained with RS-DT (87.82 ± 0.58 Pa) and highest with DS-RT (122.98 ± 1.87 Pa).

CONCLUSION

PGs, LFs, PIs and E were improved with DS-RT compared to other textile TAVs and RS-DB. Despite achieving more RSS than the rest of TAVs, DS-RT still falls within the same range of RSS produced by the other 2 valves and control exceeding the threshold for platelet activation.

Large eddy simulation in a rotary blood pump: Viscous shear stress computation and comparison with unsteady Reynolds-averaged Navier-Stokes simulation

PURPOSE:

Numerical flow analysis (computational fluid dynamics) in combination with the prediction of blood damage is an important procedure to investigate the hemocompatibility of a blood pump, since blood trauma due to shear stresses remains a problem in these devices. Today, the numerical damage prediction is conducted using unsteady Reynolds-averaged Navier-Stokes simulations. Investigations with large eddy simulations are rarely being performed for blood pumps. Hence, the aim of the study is to examine the viscous shear stresses of a large eddy simulation in a blood pump and compare the results with an unsteady Reynolds-averaged Navier-Stokes simulation.

METHODS:

The simulations were carried out at two operation points of a blood pump. The flow was simulated on a 100M element mesh for the large eddy simulation and a 20M element mesh for the unsteady Reynolds-averaged Navier-Stokes simulation. As a first step, the large eddy simulation was verified by analyzing internal dissipative losses within the pump. Then, the pump characteristics and mean and turbulent viscous shear stresses were compared between the two simulation methods.

RESULTS:

The verification showed that the large eddy simulation is able to reproduce the significant portion of dissipative losses, which is a global indication that the equivalent viscous shear stresses are adequately resolved. The comparison with the unsteady Reynolds-averaged Navier-Stokes simulation revealed that the hydraulic parameters were in agreement, but differences for the shear stresses were found.

CONCLUSION:

The results show the potential of the large eddy simulation as a high-quality comparative case to check the suitability of a chosen Reynolds-averaged Navier-Stokes setup and turbulence model. Furthermore, the results lead to suggest that large eddy simulations are superior to unsteady Reynolds-averaged Navier-Stokes simulations when instantaneous stresses are applied for the blood damage prediction.

An in vitro evaluation of turbulence after transcatheter aortic valve implantation

BACKGROUND

This study aimed at assessment of post-transcatheter aortic valve (TAV) replacement hemodynamics and turbulence when a same-size SAPIEN 3 (Edwards Lifesciences Corp, Irvine, Calif) and Medtronic Evolut (Minneapolis, Minn) were implanted in a rigid aortic root with physiological dimensions and in a representative root with calcific leaflets obtained from patient computed tomography scans.

METHODS

TAV hemodynamics were studied by placing a SAPIEN 3 26-mm and an Evolut 26-mm in rigid aortic roots and representative root with calcific leaflets under physiological conditions. Hemodynamics were assessed using high-fidelity particle image velocimetry and high-speed imaging. Transvalvular pressure gradients (PGs), pinwheeling indices, and Reynolds shear stress (RSS) were calculated.

RESULTS

(1) PGs obtained with the Evolut and the SAPIEN 3 were comparable among the different models (10.5 ± 0.15 mm Hg vs 7.76 ± 0.083 mm Hg in the rigid model along with 13.9 ± 0.19 mm Hg vs 5.0 ± 0.09 mm Hg in representative root with calcific leaflets obtained from patient computed tomography scans respectively); (2) more pinwheeling was found in the SAPIEN 3 than the Evolut (0.231 ± 0.057 vs 0.201 ± 0.05 in the representative root with calcific leaflets and 0.366 ± 0.067 vs 0.122 ± 0.045 in the rigid model); (3) higher rates of RSS were found in the Evolut (161.27 ± 3.45 vs 122.84 ± 1.76 Pa in representative root with calcific leaflets and 337.22 ± 7.05 vs 157.91 ± 1.80 Pa in rigid models). More lateral fluctuations were found in representative root with calcific leaflets.

CONCLUSIONS

(1) Comparable PGs were found among the TAVs in different models; (2) pinwheeling indices were found to be different between both TAVs; (3) turbulence patterns among both TAVs translated according to RSS were different. Rigid aortic models yield more conservative estimates of turbulence; (4) both TAVs exhibit peak maximal RSS that exceeds platelet activation 100 Pa threshold limit.

Numerical Simulation and Experimental Validation of Swirling Flow in Spiral Vortex Ventricular Assist Device

Spiral Vortex Ventricular Assist Device (SV-VAD) supports cardiac patients with refractory heart failure. Unfortunately, thrombus formation and risk of stroke due to flow complications may lead to aggravated conditions. The hemodynamics of a continuous flow in the ventricular chamber of a SV-VAD can be analyzed using numerical simulation. Particle image velocimetry and laser Doppler anemometry are utilized for validating the simulated spiral flow in a transparent acrylic SV-VAD replica based on its cross-sectional averaged axial and tangential velocities. After validation, the relationship between swirling flow and blood cell damage is established by evaluating flow effect on thrombosis due to high shear stress. Based on our analysis, stagnancy of the flow within the SV-VAD is insignificant and its low shear stress minimizes hemolysis.

Simulated Transcatheter Aortic Valve Flow: Implications of Elliptical Deployment and Under-Expansion at the Aortic Annulus

Clinical use of transcatheter aortic valves (TAVs) has been associated with abnormal deployment, including oval deployment and under-expansion when placed into calcified aortic annuli. In this study, we performed an integrated computational and experimental investigation to quantify the impact of abnormal deployment at the aortic annulus on TAV hemodynamics. A size 23 mm generic TAV computational model, developed and published previously, was subjected to elliptical deployment at the annulus with eccentricity levels up to 0.68 and to under-expansion of the TAV at the annulus by up to 25%. The hemodynamic performance was quantified for each TAV deployment configuration. TAV opening geometries were fabricated using stereolithography and then subjected to steady forward flow testing in accordance with ISO-5840. Centerline pressure profiles were compared to validate the computational model. Our findings show that slight ellipticity of the TAV may not lead to degeneration of hydrodynamic performance. However, under large ellipticity, increases in transvalvular pressure gradients were observed. Under-expanded deployment has a much greater negative effect on the TAV hemodynamics compared with elliptical deployment. The maximum turbulent viscous shear stress (TVSS) values were found to be significantly larger in under-expanded TAVs. Although the maximum value of TVSS was not large enough to cause hemolysis in all cases, it may cause platelets activation, especially for under-expanded deployments.

Characterization of erythrocyte membrane tension for hemolysis prediction in complex flows

Hemolysis is a persistent issue with blood-contacting devices. Many experimental and theoretical contributions over the last few decades have increased insight into the mechanisms of hemolysis in both laminar and turbulent flows, with the ultimate goal of developing a comprehensive, mechanistic hemolysis model. Many models assume that hemolysis scales with a resultant, scalar stress representing all components of the fluid stress tensor. This study critically evaluates this scalar stress hypothesis by calculating the response of the red blood cell membrane to different types of fluid stress (laminar shear and extension, and three turbulent shear and extension cases), each with the same scalar stress. It was found that even though the scalar stress is the same for all cases, membrane tension varied by up to three orders of magnitude. In addition, extensional flow causes constant tension, while tank-treading in shear flow causes periodic tension, with tank-treading frequency varying by three orders of magnitude among the cases. For turbulent flow, tension also depends on eddy size. It is concluded, therefore, that scalar stress alone is inadequate for scaling hemolysis. Fundamental investigations are needed to establish a new index of the fluid stress tensor that provides reliable hemolysis prediction across the wide range of complex flows that occur in cardiovascular devices.

Characterization of Transition to Turbulence for Blood in a Straight Pipe Under Steady Flow Conditions

Blood is a complex fluid that, among other things, has been established to behave as a shear thinning, non-Newtonian fluid when exposed to low shear rates (SR). Many hemodynamic investigations use a Newtonian fluid to represent blood when the flow field of study has relatively high SR (>200 s-1). Shear thinning fluids have been shown to exhibit differences in transition to turbulence (TT) compared to that of Newtonian fluids. Incorrect prediction of the transition point in a simulation could result in erroneous hemodynamic force predictions. The goal of the present study was to compare velocity profiles near TT of whole blood and Newtonian blood analogs in a straight rigid pipe with a diameter 6.35 mm under steady flow conditions. Rheology was measured for six samples of whole porcine blood and three samples of a Newtonian fluid, and the results show blood acts as a shear thinning non-Newtonian fluid. Measurements also revealed that blood viscosity at SR = 200 s-1 is significantly larger than at SR = 1000 s-1 (13.8%, p < 0.001). Doppler ultrasound (DUS) was used to measure velocity profiles for blood and Newtonian samples at different flow rates to produce Reynolds numbers (Re) ranging from 1000 to 3300 (based on viscosity at SR = 1000 s-1). Two mathematically defined methods, based on the velocity profile shape change and turbulent kinetic energy (TKE), were used to detect TT. Results show similar parabolic velocity profiles for both blood and the Newtonian fluid for Re < 2200. However, differences were observed between blood and Newtonian fluid velocity profiles for larger Re. The Newtonian fluid had blunt-like velocity profiles starting at Re = 2403 ± 8 which indicated transition. In contrast, blood did not show this velocity profile change until Re = 2871 ± 104. The Newtonian fluid had large velocity fluctuations (root mean square (RMS) > 20%) with a maximum TKE near the pipe center at Re = 2316 ± 34 which indicated transition. In contrast, blood results showed the maximum TKE at Re = 2806 ± 109. Overall, the critical Re was delayed by ∼20% (p < 0.001) for blood compared to the Newtonian fluid. Thus, a Newtonian assumption for blood at flow conditions near transition could lead to large errors in velocity prediction for steady flow in a straight pipe. However, these results are specific to this pipe diameter and not generalizable since SR is highly dependent on pipe diameter. Further research is necessary to understand this relation in different pipe sizes, more complex geometries, and under pulsatile flow conditions.

Hydrodynamic Simulation of an Orbital Shaking Test for the Degradation Assessment of Blood-Contact Biomedical Coatings

Biomedical coatings are used to promote the wear resistance and the biocompatibility of a mechanical heart valve. An orbital shaking test was proposed to assess the durability of the coatings by the amount material eroded by the surrounding fluid. However, there is still a lack of understanding with regards to the shaker’s rotating conditions and the corresponding physiological condition. This study implemented numerical simulations by establishing a fluid dynamic model to evaluate the intensity of the shear stress under various rotating speeds and diameters of the shaker. The results are valuable to conduct in vitro tests for estimating the performance of biomedical coatings under real hemodynamic conditions and can be applied to other fluid-contact implants.

Turbulence measurements in an axial rotary blood pump with laser Doppler velocimetry

BACKGROUND

The implantation of rotary blood pumps as ventricular assist devices (VADs) has become a viable therapy for quite a number of patients with end-stage heart failure. However, these rotary blood pumps cause adverse events that are related to blood trauma. It is currently believed that turbulence in the pump flow plays a significant role. But turbulence has not been measured to date because there is no optical access to the flow space in rotary blood pumps because of their opaque casings.

METHODS

This difficulty is overcome with a scaled-up model of the HeartMate II (HM II) rotary blood pump with a transparent acrylic housing. A 2-component laser Doppler velocimetry (LDV) system was used for the measurement of time resolved velocity profiles and velocity spectra upstream and downstream of the rotor blades. Observing similarity laws, the speed and pump head were adjusted to correspond closely to the design point of the original pump - 10,600 rpm speed and 80 mmHg pressure head. A model fluid consisting of a water-glycerol mixture was used.

RESULTS

The measured velocity spectra were scalable by the Kolmogorov length and the Kolmogorov length was estimated to be between 14 and 24 µm at original scale, thus being about 1.5 to 3 times the size of a red blood cell.

CONCLUSIONS

It can be concluded that turbulence is indeed present in the investigated blood pump and that it can be described by Kolmogorov's theory of turbulence. The size of the smallest vortices compares well to the turbulence length scales as found in prosthetic heart valves, for example.

Investigating the impact of non-Newtonian blood models within a heart pump

A detailed computational fluid dynamics (CFD) study of transient, turbulent blood flow through a positive displacement left ventricular assist device is performed. Two common models for non-Newtonian blood flow are compared to the Newtonian model to investigate their impact on predicted levels of shear rate and wall shear stress. Given that both parameters are directly relevant to the evaluation of risk from thrombus and haemolysis, there is a need to assess the sensitivity to modelling non-Newtonian flow effects within a pulsatile turbulent flow, in order to identify levels of uncertainly in CFD. To capture the effects of turbulence, the elliptic blending Reynolds stress model is used in the present study, on account of superior performance of second moment closure schemes previously identified by the present authors. The CFD configuration includes two cyclically rotating valves and a moving pusher plate to periodically vary the chamber volume. An overset mesh algorithm is used for each instance of mesh motion, and a zero gap technique was employed to ensure full valve closure. The left ventricular assist device was operated at a pumping rate of 86 BPM (beats per minute) and a systolic duration of 40% of the pumping cycle, in line with existing experimental data to which comparisons are made. The sensitivity of the variable viscosity models is investigated in terms of mean flow field, levels of turbulence and global shear rate, and a non-dimensional index is used to directly evaluate the impact of non-Newtonian effects. The clinical relevance of the results is reported along with a discussion of modelling uncertainties, observing that the turbulent kinetic energy is generally predicted to be higher in non-Newtonian flow than that observed in Newtonian flow. Copyright © 2016 John Wiley & Sons, Ltd.

Assessment of turbulent viscous stress using ICOSA 4D Flow MRI for prediction of hemodynamic blood damage

Flow-induced blood damage plays an important role in determining the hemodynamic impact of abnormal blood flow, but quantifying of these effects, which are dominated by shear stresses in highly fluctuating turbulent flow, has not been feasible. This study evaluated the novel application of turbulence tensor measurements using simulated 4D Flow MRI data with six-directional velocity encoding for assessing hemodynamic stresses and corresponding blood damage index (BDI) in stenotic turbulent blood flow. The results showed that 4D Flow MRI underestimates the maximum principal shear stress of laminar viscous stress (PLVS), and overestimates the maximum principal shear stress of Reynolds stress (PRSS) with increasing voxel size. PLVS and PRSS were also overestimated by about 1.2 and 4.6 times at medium signal to noise ratio (SNR) = 20. In contrast, the square sum of the turbulent viscous shear stress (TVSS), which is used for blood damage index (BDI) estimation, was not severely affected by SNR and voxel size. The square sum of TVSS and the BDI at SNR >20 were underestimated by less than 1% and 10%, respectively. In conclusion, this study demonstrated the feasibility of 4D Flow MRI based quantification of TVSS and BDI which are closely linked to blood damage.

An Approach for Assessing Turbulent Flow Damage to Blood in Medical Devices

In this work, contributing factors for red blood cell (RBC) damage in turbulence are investigated by simulating jet flow experiments. Results show that dissipative eddies comparable or smaller in size to the red blood cells cause hemolysis and that hemolysis corresponds to the number and, more importantly, the surface area of eddies that are associated with Kolmogorov length scale (KLS) smaller than about 10 μm. The size distribution of Kolmogorov scale eddies is used to define a turbulent flow extensive property with eddies serving as a means to assess the turbulence effectiveness in damaging cells, and a new hemolysis model is proposed. This empirical model is in agreement with hemolysis results for well-defined systems that exhibit different exposure times and flow conditions, in Couette flow viscometer, capillary tube, and jet flow experiments.

Reversible binding of hemoglobin to band 3 constitutes the molecular switch that mediates O2 regulation of erythrocyte properties

Functional studies have shown that the oxygenation state of the erythrocyte regulates many important pathways, including glucose metabolism, membrane mechanical stability, and cellular adenosine triphosphate (ATP) release. Deoxyhemoglobin (deoxyHb), but not oxyhemoglobin, binds avidly and reversibly to band 3, the major erythrocyte membrane protein. Because band 3 associates with multiple metabolic, solute transport, signal transduction, and structural proteins, the hypothesis naturally arises that the O₂-dependent regulation of erythrocyte properties might be mediated by the reversible association of deoxyHb with band 3. To explore whether the band 3-deoxyHb interaction constitutes a "molecular switch" for regulating erythrocyte biology, we have generated transgenic mice with mutations in the deoxyHb-binding domain of band 3. One strain of mouse contains a "humanized" band 3 in which the N-terminal 45 residues of mouse band 3 are replaced by the homologous sequence from human band 3, including the normal human band 3 deoxyHb-binding site. The second mouse contains the same substitution as the first, except the deoxyHb site on band 3 (residues 12-23) has been deleted. Comparison of these animals with wild-type mice demonstrates that the following erythrocyte properties are controlled by the O₂-dependent association of hemoglobin with band 3: (1) assembly of a glycolytic enzyme complex on the erythrocyte membrane which is associated with a shift in glucose metabolism between the pentose phosphate pathway and glycolysis, (2) interaction of ankyrin with band 3 and the concomitant regulation of erythrocyte membrane stability, and (3) release of ATP from the red cell which has been linked to vasodilation.

Flow-Induced Damage to Blood Cells in Aortic Valve Stenosis

Valvular hemolysis and thrombosis are common complications associated with stenotic heart valves. This study aims to determine the extent to which hemodynamics induce such traumatic events. The viscous shear stress downstream of a severely calcified bioprosthetic valve was evaluated via in vitro 2D particle image velocimetry measurements. The blood cell membrane response to the measured stresses was then quantified using 3D immersed-boundary computational simulations. The shear stress level at the boundary layer of the jet flow formed downstream of the valve orifice was observed to reach a maximum of 1000-1700 dyn/cm(2), which was beyond the threshold values reported for platelet activation (100-1000 dyn/cm(2)) and within the range of thresholds reported for red blood cell (RBC) damage (1000-2000 dyn/cm(2)). Computational simulations demonstrated that the resultant tensions at the RBC membrane surface were unlikely to cause instant rupture, but likely to lead to membrane plastic failure. The resultant tensions at the platelet surface were also calculated and the potential damage was discussed. It was concluded that although shear-induced thrombotic trauma is very likely in stenotic heart valves, instant hemolysis is unlikely and the shear-induced damage to RBCs is mostly subhemolytic.

Lagrangian methods for blood damage estimation in cardiovascular devices--How numerical implementation affects the results

This paper evaluated the influence of various numerical implementation assumptions on predicting blood damage in cardiovascular devices using Lagrangian methods with Eulerian computational fluid dynamics. The implementation assumptions that were tested included various seeding patterns, stochastic walk model, and simplified trajectory calculations with pathlines. Post processing implementation options that were evaluated included single passage and repeated passages stress accumulation and time averaging. This study demonstrated that the implementation assumptions can significantly affect the resulting stress accumulation, i.e., the blood damage model predictions. Careful considerations should be taken in the use of Lagrangian models. Ultimately, the appropriate assumptions should be considered based the physics of the specific case and sensitivity analysis, similar to the ones presented here, should be employed.