Theory of non-Newtonian viscosity of red blood cell suspension: effect of red cell deformation

The effects of the deformation of red blood cells on non-Newtonian viscosity of a concentrated red cell suspension are investigated theoretically. To simplify the problem an elastic spherical shell filled with an incompressible Newtonian fluid is considered as a model of a normal red cell. The equation of the surface of the shell suspended in a steady simple shear flow is calculated on the assumption that the deformation from a spherical shape is very small. The relative viscosity of a concentrated suspension of such particles is obtained based on the "free surface cell" method proposed by Happel. It is shown that the relative viscosity decreases as the shear rate increases.     

Cyclic loading on the red cell membrane in a shear flow: a possible cause of hemolysis

The fluid force acting on single human red cells in a high shear flow was analyzed. A two-dimensional elliptical microcapsule as a model of the deformed red cells was adopted to numerically calculate the distributions of the shear forces on both sides of the cell membrane. It is theoretically shown that the cell membrane undergoes an unsteady cyclic loading under the rotational motion around the interior. The mechanism leading to blood cell trauma is examined by repeatedly loading the continuously moving cell membrane.     

Analysis of red blood cell motion through cylindrical micropores: effects of cell properties.

Filtration through micropores is frequently used to assess red blood cell deformability, but the dependence of pore transit time on cell properties is not well understood. A theoretical model is used to simulate red cell motion through cylindrical micropores with diameters of 3.6, 5, and 6.3 microns, and 11-microns length, at driving pressures of 100-1000 dyn/cm2. Cells are assumed to have axial symmetry and to conserve surface area during deformation. Effects of membrane shear viscosity and elasticity are included, but bending resistance is neglected. A time-dependent lubrication equation describing the motion of the suspending fluid is solved, together with the equations for membrane equilibrium, using a finite difference method. Predicted transit times are consistent with previous experimental observations. Time taken for cells to enter pores represents more than one-half of the transit time. Predicted transit time increases with increasing membrane viscosity and with increasing cell volume. It is relatively insensitive to changes in internal viscosity and to changes in membrane elasticity except in the narrowest pores at low driving pressures. Elevating suspending medium viscosity does not increase sensitivity of transit time to membrane properties. Thus filterability of red cells is sensitively dependent on their resistance to transient deformations, which may be a key determinant of resistance to blood flow in the microcirculation.     

Motion and deformation of a red blood cell in a shear flow: a two-dimensional simulation of the wall effect

The motion and deformation of a single red blood cell in a simple shear flow between two parallel walls is studied theoretically. A two-dimensional deformable microcapsule is adopted as a model for the cell, which has a thin moving membrane, like a tank-tread, around the interior and is deformed into an elliptical shape with a constant area. Applying the finite element method to the Stokes equations, the tank-tread motion and deformation is determined in a stationary motion, under fluid dynamic interaction between the cell and the walls. It is shown that the motion and deformation of the microcapsule crucially depends on the channel width between the two walls. As the width decreases, the microcapsule is more elongated and the frequency of tank-tread motion decreases at a constant shear rate. In addition, the angle of inclination decreases at the low range of the viscosity ratio of internal to external fluids and increases at the high range. The results obtained are compared with experimental observations and applied to the behavior of cells under mutual interaction.     

Aggregation and disaggregation of red blood cells

The aggregation of red blood cells may be analyzed as an interaction of an adhesive surface energy and the elastic stored energy which results from deformation of the cell. The adhesive surface energy is the work required to separate a unit adhered area and is the resultant of adhesive forces due to the bridging molecules that bind the cells together and the electrostatic repulsion due to surface charge. The elastic strain energy in the case of the red blood is associated with the membrane elasticity only since the interior of the cell is liquid. The membrane elasticity is due both to bending stiffness and shear. The area expansion is small and may be neglected. These assumptions allow realistic computation of red cell shapes in rouleaux. The disaggregation of rouleaux requires an external force which must overcome the adhesive energy and also supply additional elastic energy of deformation. Depending on the geometry, the initial effect of elastic energy may tend to aid disaggregation. In a shear flow, the stresses on a suspended rouleau alternately tend to compress and to disaggregate the cells if they are free to rotate. This introduces a time dependence so that viscous effects due to the viscosity of the cell membrane, the cell cytoplasm and the external fluid may play a role in determining whether disaggregation proceeds to completion or not.     

Tank-tread frequency of the red cell membrane: dependence on the viscosity of the suspending medium

Single human red cells were suspended in media with viscosities ranging from 12.9 to 109 mPa s and subjected to shear flow ranging from 1/s to 290/s in a rheoscope. This is a transparent cone-plate chamber adapted to a microscope. The motion of the membrane around red cells oriented in a steady-state fashion in the shear field (tank-tread motion) was videotaped. The projected length and width of the cells as well as the frequency of tank-tread motion were measured. One-thousand eight-hundred seventy-three cells of three blood donors were evaluated. The frequency increased with the mean shear rate in an almost linear fashion. The slope of this dependence increased weakly with the viscosity of the suspending medium. No correlation was found between the frequency and four morphological red cell parameters: the projected length and width of the cells as well as the ratio and the square root of the product of these quantities. The energy dissipation within the red cell membrane was estimated based on the measured parameters and compared to the energy dissipation in the undisturbed shear flow. At constant mean shear rate the rise of the energy dissipation with viscosity is slower whereas at constant viscosity the rise with the shear rate is steeper than in the undisturbed shear flow. A fit of the data collected in this work to a theoretical red cell model might allow one to determine intrinsic mechanical constants in the low deformation regime.     

Effects of aggregation on the flow properties of red blood cell suspensions in narrow vertical tubes

The flow properties of aggregating red cell suspensions flowing at low rates through vertical tubes with diameters from 30 microns to 150 microns are analyzed using a theoretical model. Unidirectional flow is assumed, and the distributions of velocity and red cell concentration are assumed to be axisymmetric. A three-layer approximation is used for the distribution of red cells, with a cylindrical central core of aggregated red cells moving with uniform velocity, a cell-free marginal layer near the tube wall, and an annular region located between the core and the marginal layer containing suspended non-aggregating red cells. This suspension is assumed to behave approximately as a Newtonian fluid whose viscosity increases exponentially with red cell concentration. Physical arguments concerning the mechanics of red cell attachment to, and detachment from the aggregated core lead to a kinetic equation for core formation. From this kinetic equation and the equation for conservation of red cell volume flux, a relationship between core radius and pressure gradient is obtained. Then the relative viscosity is calculated as a function of pseudo-shear rate. At low flow rates, it is shown that the relative viscosity decreases with decreasing flow and that the dependence of relative viscosity on shear rates is more pronounced in larger tubes. It is also found that the relative viscosity decreases with increasing aggregation tendency of suspension. These theoretical predictions are in good qualitative and quantitative agreement with experimental results.     

Determination of red blood cell membrane viscosity from rheoscopic observations of tank-treading motion.

Measurements of the dimensions and membrane rotational frequency of individual erythrocytes steadily tank-treading in a rheoscope are used to deduce the surface shear viscosity of the membrane. The method is based on an integral energy principle which says that the power supplied to the tank-treading cell by the suspending fluid is equal to the rate at which energy is dissipated by viscous action in the membrane and cytoplasm. The integrals involved are formulated with the aid of an idealized mathematical model of the tank-treading red blood cell (RBC) (Keller and Skalak, 1982, J. Fluid Mech., 120:24-27) and evaluated numerically. The outcome is a surface-averaged value of membrane viscosity which is representative of a finite interval of membrane shear rate. The numerical values computed show a clear shear-thinning characteristic as well as a significant augmentation of viscosity with cell age and tend toward agreement with those determined for the rapid phase of shape recovery in micropipettes (Chien, S., K.-L. P. Sung, R. Skalak, S. Usami, and A. Tozeren, 1978, Biophys. J., 24:463-487). The computations also indicate that the rate of energy dissipation in the membrane is always substantially greater than that in the cytoplasm.     

Detection of red cell aggregation by low shear rate viscometry in whole blood with elevated plasma viscosity

The viscosity of whole blood measured at low shear rates is determined partly by shear resistance of the red cell aggregates present, stronger aggregation increasing the viscosity in the absence of other changes. Effects of cell deformability can confound interpretation and comparison in terms of aggregation, however, particularly when the plasma viscosity is high. We illustrate the problem with a comparison of hematocrit-adjusted blood from type 1 diabetes patients and controls in which it is found the apparent and relative viscosities at a true shear rate of 0.20 s-1 are lower in the patient samples than age matched controls, in spite of reports that aggregation is increased in such populations. Because the plasma viscosities of the patients were higher on average than controls, we performed a series of experiments to examine the effect of plasma protein concentration and viscosity on normal blood viscosity. Dilution or concentration by ultrafiltration of autologous plasma and viscosity measurements at low shear on constant hematocrit red cell suspensions showed (a) suspension viscosity at 0.25 and 3 s-1 increased monotonically with plasma protein concentration and viscosity but (b) the relative viscosity increased, in concert with the microscopic aggregation grade, up to a viscosity of approximately 1.25 mPa-s but above this the value the relative viscosity no longer increased as the degree of aggregation increased in concentrated plasmas. It is suggested that in order to reduce cell deformation effects in hyperviscous pathological plasmas, patient and control plasmas should be systematically diluted before hematocrit is adjusted and rheological measurements are made. True shear rates should be calculated. Comparison of relative viscosities at low true shear rates appears to allow the effects of red cell aggregation to be distinguished by variable shear rate viscometry in clinical blood samples.     

Static equilibrium configurations of a model red blood cell

Summary The membrane of the red blood cell is modeled as a fluid shell which resists bending and changes in area. The differential equations governing the mechanical equilibrium of such a membrane are derived and axisymmetric solutions are obtained numerically.     

Shape changes of the human red cell studied by aqueous two-phase partition

Aqueous two-phase systems have been used to study the human red cell during metabolically induced shape changes. When the discoid character of the cells was lost in favour of echinocytic forms, the partition increased both in charge-sensitive and in charge-insensitive two-phase systems. Reversal of the shape transformation by ATP repletion not only led to shape recovery but also restored the initial partition. Therefore it is apparent that red cells exhibit a shape-dependent partition behaviour. As the partition is dependent on surface properties (such as charge and hydrophobicity) of the partitioned material, the results show that the shape changes caused rearrangement of the membrane and thereby exposure of or greater accessibility of binding groups on the cell surface. The similar partition behaviour in the charge-sensitive and charge-insensitive phase systems show that the increased partition was caused mainly by increased hydrophobic interactions between the cells and the upper phase. The observed partition behaviour therefore suggests that the echinocytic cells acquire a higher affinity for the upper phase by repacking the lipid bilayer or at least the outer leaflet into a less efficient packed and thus more fluid structure.     

Temperature dependence of the viscoelastic recovery of red cell membrane.

The time-dependent recovery of an elongated red cell is studied as a function of temperature. Before release, the elongated cell is in static equilibrium where external forces are balanced by surface elastic force resultants. Upon release, the cell recovers its initial shape with a time-dependent exponential behavior characteristic of a viscoelastic solid material undergoing large ("finite") deformation. The recovery process is characterized by a time constant, tc, that decreases from approximately 0.27 s at 6 degrees C to 0.06 s at 37 degrees C. From this measurement of the time constant and an independent measurement of the shear modulus of surface elasticity for red cell membrane, the value for the membrane surface viscosity as a function of temperature can be calculated.     

Elastic thickness compressibilty of the red cell membrane

We have used an ultrasensitive force probe and optical interferometry to examine the thickness compressibility of the red cell membrane in situ. Pushed into the centers of washed-white red cell ghosts lying on a coverglass, the height of the microsphere-probe tip relative to its closest approach on the adjacent glass surface revealed the apparent material thickness, which began at approximately 90 nm per membrane upon detection of contact (force approximately 1-2 pN). With further impingement, the apparent thickness per membrane diminished over a soft compliant regime that spanned approximately 40 nm and stiffened on approach to approximately 50 nm under forces of approximately 100 pN. The same force-thickness response was obtained on recompression after retraction of the probe, which demonstrated elastic recoverability. Scaled by circumferences of the microspheres, the forces yielded energies of compression per area which exhibited an inverse distance dependence resembling that expected for flexible polymers. Attributed to the spectrin component of the membrane cytoskeleton, the energy density only reached one thermal energy unit (k(B)T) per spectrin tetramer near maximum compression. Hence, we hypothesized that the soft compliant regime probed in the experiments represented the compressibility of the outer region of spectrin loops and that the stiff regime < 50 nm was the response of a compact mesh of spectrin backed by a hardcore structure. To evaluate this hypothesis, we used a random flight theory for the entropic elasticity of polymer loops to model the spectrin network. We also examined the possibility that additional steric repulsion and apparent thickening could arise from membrane thermal-bending excitations. Fixing the energy scale to k(B)T/spectrin tetramer, the combined elastic response of a network of ideal polymer loops plus the membrane steric interaction correlated well with the measured dependence of energy density on distance for a statistical segment length of approximately 5 nm for spectrin (i.e., free chain end-to-end length of approximately 29 nm) and a hardcore limit of approximately 30 nm for underlying structure.     

Influence of preparative procedures on the membrane viscoelasticity of human red cell ghosts

The effects of systematic variations in the preparative procedures on the membrane viscoelastic properties of resealed human red blood cell ghosts have been investigated. Ghosts, prepared by hypotonic lysis at 0 degrees C and resealing at 37 degrees C, were subjected to: measurement of the time constant for extensional recovery (tc); measurement of the membrane shear elastic modulus (mu) via three separate techniques; determination of the membrane viscosity (eta m) via a cone-plate Rheoscope. Membrane viscosity was also determined as eta m = mu X tc. Compared to intact cells, ghosts had shorter tc, regardless of their residual hemoglobin concentration (up to 21.6 g/dl). However, prolonged exposure to hypotonic media did increase their recovery time toward the intact cell value. The shear elastic modulus, as judged by micropipette aspiration of membrane tongues (mu p), was similar for all ghosts and intact cells. This result, taken with the tc data, indicates that ghosts have reduced membrane viscosity. Rheoscopic analysis also showed that eta m was reduced for ghosts, with the degree of reduction (approx. 50%) agreeing well with that estimated by the product mu p X tc. However, flow channel and pipette elongation estimates indicated that the ghost membrane elastic modulus was somewhat elevated compared to intact cells. We conclude that: ghosts have reduced membrane viscosity; ghosts have membrane rigidities close to intact cells, except possibly when the membrane is subjected to very large strains; the reduction in eta m is not directly related to the loss of hemoglobin; prolonged exposure of ghosts to low-ionic strength media increases the membrane viscosity toward its initial cellular level. These data indicate that the mechanical characteristics of ghost membranes can be varied by changing the methods of preparation and thus have potential application to further studies of the structural determinants of red cell membrane viscoelasticity.     

Multiphase hemodynamic simulation of pulsatile flow in a coronary artery

A multiphase transient non-Newtonian three-dimensional (3-D) computational fluid dynamics (CFD) simulation has been performed for pulsatile hemodynamics in an idealized curved section of a human coronary artery. We present the first prediction, to the authors' knowledge, of particulate buildup on the inside curvature using the multiphase theory of dense suspension hemodynamics. In this study, the particulates are red blood cells (RBCs). The location of RBC buildup on the inside curvature correlates with lower wall shear stress (WSS) relative to the outside curvature. These predictions provide insight into how blood-borne particulates interact with artery walls and hence, have relevance for understanding atherogenesis since clinical observations show that atherosclerotic plaques generally form on the inside curvatures of arteries. The buildup of RBCs on the inside curvature is driven by the secondary flow and higher residence times. The higher viscosity in the central portion of the curved vessel tends to block their flow, causing them to migrate preferentially through the boundary layer. The reason for this is the nearly neutrally buoyant nature of the dense two-phase hemodynamic flow. The two-phase non-Newtonian viscosity model predicts greater shear thinning than the single-phase non-Newtonian model. Consequently, the secondary flow induced in the curvature is weaker. The waveforms for computed hemodynamic parameters, such as hematocrit, WSS, and viscosity, follow the prescribed inlet velocity waveforms. The lower oscillatory WSS produced on the inside curvature has implications for understanding thickening of the intimal layer.     

Theoretical model and experimental study of red blood cell (RBC) deformation in microchannels

The motion and deformation of red blood cells (RBCs) flowing in a microchannel were studied using a theoretical model and a novel automated rheoscope. The theoretical model was developed to predict the cells deformation under shear as a function of the cells geometry and mechanical properties. Fluid dynamics and membrane mechanics are incorporated, calculating the traction and deformation in an iterative manner. The model was utilized to evaluate the effect of different biophysical parameters, such as: inner cell viscosity, membrane shear modulus and surface to volume ratio on deformation measurements. The experimental system enables the measurement of individual RBCs velocity and their deformation at defined planes within the microchannel. Good agreement was observed between the simulation results, the rheoscope measurements and published ektacytometry results. The theoretical model results imply that such deformability measuring techniques are weakly influenced by changes in the inner viscosity of the cell or the ambient fluid viscosity. However, these measurements are highly sensitive to RBC shear modulus. The shear modulus, estimated by the model and the rheoscope measurements, falls between the values obtained by micropipette aspiration and laser trapping. The study demonstrates the integration of a theoretical model with a microfabricated device in order to achieve a better understanding of RBC mechanics and their measurement using microfluidic shear assays. The system and the model have the potential of serving as quantitative clinical tools for diagnosing deformability disorders in RBCs.     

Haemoglobin and the red cell membrane.

The results presented here indicate that haemoglobin is an integral part of the red cell membrane. The haemoglobin content of the membrane is highly dependent on the Ca++ content of the membrane in health and disease. Changes in the red cell interior alter the whole organization of the membrane and are even reflected in the binding of immunoglobulins to the red cell surface. The preferential binding of Hb-s A2 and S to the membrane has been confirmed. This phenomenon cannot be explained by differences in the charge between these haemoglobins and Hb A.     

Alteration of red cell membrane viscoelasticity by heat treatment: effect on cell deformability and suspension viscosity

The membrane shear elastic modulus (mu) and the time constant for extensional shape recovery (tc) were measured for normal, control human red blood cells (RBC) and for RBC heat treated (HT) at 48 degrees C. Three separate methods for the measurement of mu were compared (two used a micropipette and one employed a flow channel), and the membrane viscosity (n) was calculated from the relation n = mu. tc. The deformability of HT and control cells was evaluated using micropipette techniques, and the bulk viscosity of RBC suspensions at 40% hematocrit was measured. The shear elastic modulus, or "membrane rigidity", was more than doubled by heat treatment, although both the absolute value for mu and the estimate of the increase induced by heat treatment varied depending on the method of measurement. Heat treatment caused smaller increases in membrane viscosity and in membrane bending resistance, and only minimal changes in cell geometry. The deformability of HT cells was reduced: 1) the pressure required for cell entry (Pe) into 3 micrometers pipettes was increased, on average, by 170%; 2) at an aspiration pressure (Pa) exceeding Pe, longer times were required for cell entry into the same pipettes. However, when Pa was scaled relative to the mean entry pressure for a given sample (i.e, Pa/Pe), entry times were similar for control and HT cells. Bulk viscosity of HT RBC suspensions was elevated by approximately 12% on average (shear rates 75 to 1500 inverse seconds). These findings suggest that alteration of RBC membrane mechanical properties, similar to those induced by heat treatment, would most affect the in vivo circulation in regions where vessel dimensions are smaller than cellular diameters.     

Effects of sedimentation of small red blood cell aggregates on blood flow in narrow horizontal tubes

The flow properties of aggregating red cell suspensions flowing at low flow rates through horizontal tubes are analyzed using a theoretical model. The effects of sedimentation of small aggregates, which will be formed at comparatively high flow rates, on the relative apparent viscosity are considered. In the case in which a large number of small aggregates are formed in a suspension flowing through a horizontal tube, it seems that red cells are transported as a concentrated suspension through the bottom part of the tube because of sedimentation of aggregates. A two-layer flow model is used for the distribution of red cells. It consists of plasma in the upper part and a concentrated red cell suspension in the bottom part of the tube divided by a smooth and horizontal interface. It is assumed that the suspension is a Newtonian fluid whose viscosity increases exponentially with hematocrit. The velocity distribution, the relative apparent viscosity and the flux of red cells are calculated as functions of width of plasma layer for a different discharge hematocrit. The theoretical results are compared with the results obtained from experimental data. The relative apparent viscosity increases rapidly with an increasing degree of sedimentation over a wide range of plasma layer widths.     

Ektacytometric analysis of factors regulating red cell deformability

Photometric analysis of laser diffraction patterns has been used to obtain quantitative measurements of deformability of specifically modified normal red cells. Variation of deformability with suspending medium osmolality and with applied shear stress was used to distinguish between changes in internal viscosity, surface area-to-volume ratio, and viscoelastic properties of the membrane in their influence on whole cell deformability.