A newly designed oscillating viscometer for blood viscosity measurements

A newly designed type of oscillating viscometer is described. The viscometer consists of either a tube or a rod oscillating at a resonance frequency with amplitudes in the micro- and nanometer range. A fluid flowing through the tube or surrounding the rod damps the torsional oscillations. The increase in the damping depends on the viscosity of the fluid and is used to determine viscosity. It was found that viscosity measurements are feasible during blood flow. This new type of viscometer may be useful to the study of biophysical properties of blood at the wall surface during flow and give new insights into blood flow. The device allows direct viscosity measurement on blood directly as it is drawn from the vein through the tube without any anticoagulant.     

Effects of flow geometry on blood viscoelasticity

The viscoelastic properties of blood are dominated by microstructures formed by red cells. The microstructures are of several types such as irregular aggregates, rouleaux, and layers of aligned cells. The dynamic deformability of the red cells, aggregation tendency, cell concentration, size of confining vessel and rate of flow are determining factors in the microstructure. Viscoelastic properties, viscosity and elasticity, relate to energy loss and storage in flowing blood while relaxation time and Weissenberg number play a role in assessing the importance of the elasticity relative to the viscosity. These effects are shown herein for flow in a large straight cylindrical tube, a small tube, and a porous medium. These cases approximate the geometries of the arterial system: large vessels, small vessels and vessels with many branches and bifurcations. In each case the viscosity, elasticity, relaxation time and Weissenberg number for normal human blood as well as blood with enhanced cell aggregation tendency and diminished cell deformability are given. In the smaller spaces of the microtubes and porous media, the diminished viscosity shows the possible influence of the F?hraeus-Lindqvist effect and at high shear rates, the viscoelasticity of blood shows dilatancy. This is true for normal, aggregation enhanced and hardened cells.     

Blood flow: Theory,effective viscosity and effects of particle distribution

A study is made of blood flow by assuming that the blood constitutes a suspension of cells in plasma instead of a simple homogeneous fluid. A macroscopic theory governing the motion of plasma in a plasma-cell system is derived from the local volume averaging method for a system without mass transfer between the phases, and its characteristic length is much larger than the size of the cells. The equations governing the motion of the local averaged fluid quantities include one additional term in the equation of motion and two additional terms in the energy equation. These terms represent, respectively, the force exerted upon the fluid by the particles, and the rate of heat transfer and work kone upon the fluid by the particles. The theory is applied to obtain the effective viscosity as the explicit function of the volume concentration of the cells by assuming that the cells behave like rigid spherical particles with slip-collision, and the plasma is an compressible Newtonian fluid. Comparison with existing experimental results shows a good agreement. The theory is also used to obtain the effects of cell distribution upon the overall effective viscosity in a circular tube. The quantitative result shows that there is a decrease in overall effective viscosity as the concentration of cells increases toward the center of the tube, and the overall effective viscosity is smaller than the flow with evenly distributed cells.     

Numerical Simulation of Blood Flow Through Microvascular Capillary Networks

A numerical method is implemented for computing blood flow through a branching microvascular capillary network. The simulations follow the motion of individual red blood cells as they enter the network from an arterial entrance point with a specified tube hematocrit, while simultaneously updating the nodal capillary pressures. Poiseuille’s law is used to describe flow in the capillary segments with an effective viscosity that depends on the number of cells residing inside each segment. The relative apparent viscosity is available from previous computational studies of individual red blood cell motion. Simulations are performed for a tree-like capillary network consisting of bifurcating segments. The results reveal that the probability of directional cell motion at a bifurcation (phase separation) may have an important effect on the statistical measures of the cell residence time and scattering of the tube hematocrit across the network. Blood cells act as regulators of the flow rate through the network branches by increasing the effective viscosity when the flow rate is high and decreasing the effective viscosity when the flow rate is low. Comparison with simulations based on conventional models of blood flow regarded as a continuum indicates that the latter underestimates the variance of the hematocrit across the vascular tree.     

Microcontinuum model for pulsatile blood flow through a stenosed tube

The effects of polar nature of blood and pulsatility on flow through a stenosed tube have been analysed by assuming blood as a micropolar fluid. Linearized solutions of basic equations are obtained through consecutive applications of finite Hankel and Laplace transforms. The analytical expressions for axial and particle angular velocities, wall shear stress, resistance to flow and apparent viscosity have been obtained. The axial velocity profiles for Newtonian and micropolar fluids have been compared. The interesting observation of this analysis is velocity, in certain parts of cycle, for micropolar fluid is higher than Newtonain fluid. Variation of apparent viscosity eta a with tube radius shows both inverse Fahraeus-Lindqvist and Fahraeus-Lindqvist effects. Finally, the resistance to flow and wall shear stress for normal and diseased blood have been computed and compared.     

Particulate nature of blood determines macroscopic rheology: a 2-D lattice Boltzmann analysis

Historically, predicting macroscopic blood flow characteristics such as viscosity has been an empirical process due to the difficulty in rigorously including the particulate nature of blood in a mathematical representation of blood rheology. Using a two-dimensional lattice Boltzmann approach, we have simulated the flow of red blood cells in a blood vessel to estimate flow resistance at various hematocrits and vessel diameters. By including white blood cells (WBCs) in the flow, we also calculate the increase in resistance due to white cell rolling and adhesion. The model considers the blood as a suspension of particles in plasma, accounting for cell-cell and cell-wall interactions to predict macroscopic blood rheology. The model is able to reproduce the Fahraeus-Lindqvist effect, i.e., the increase in relative apparent viscosity as tube size increases, and the Fahraeus effect, i.e., tube hematocrit is lower than discharge hematocrit. In addition, the model allows direct assessment of the effect of WBCs on blood flow in the microvasculature, reproducing the dramatic increases in flow resistance as WBCs enter short capillary segments. This powerful and flexible model can be used to predict blood flow properties in any vessel geometry and with any blood composition.     

Effects of velocity jumps on blood flow in large arteries

When a human being experiences a sudden velocity change, the blood flow is disturbed. A theoretical analysis to predict the effects of sudden velocity changes on blood flow in large arteries is presented. The situations is modelled as a one-dimensional flow problem in a viscoelastic tube where the fluid viscosity convective term in the equation of motion and nonlinearity in the elastic modulus of the tube wall are neglected. The governing equations of the model are solved by Laplace transformation. The computed results show that relatively high blood pressures, capable of harming circulation, are produced even by relatively moderate velocity jumps.     

Linear and nonlinear analyses of pulsatile blood flow in a cylindrical tube

A non-Newtonian shear-thinning constitutive relation is proposed to study pulsatile flow of whole blood in a cylindrical tube. The constitutive relation, which satisfies the principle of material frame indifference, is derived from viscometric data obtained from whole blood over a range of hematocrits. Assuming axisymmetric flow in a rigid cylindrical tube of constant diameter, a second-order, nonlinear partial differential equation governing the axial velocity component is obtained. Imposing a periodic pressure gradient, the governing equation was solved numerically using finite difference methods over a range of Stokes values and hematocrits. For a forcing frequency of 1 Hz, results are presented over tube diameters ranging between 0.1 and 2 cm and over hematocrits ranging between 10 and 80%. For a given hematocrit, velocity profiles predicted for the non-Newtonian model under sinusoidal forcing reveal attenuated volume flow rate and enhanced vorticity transport over the tube cross-section relative to a Newtonian fluid having a viscosity corresponding to the high shear-rate limit. For moderate to high Stokes numbers, consistent with flow in large arteries, our results revealed a viscosity distribution that was nearly time invariant. An analytic solution was obtained for a fluid having arbitrarily prescribed radially varying, temporally invariant viscosity and density distributions under arbitrary periodic pressure forcing. Close agreement was observed between our numerical and analytical results when the imposed viscosity distribution was chosen to approximate the time-averaged viscosity distribution predicted by the shear-thinning non-Newtonian model. For St > or approximately= 100, the disparity between our results and those of a Newtonian fluid of constant viscosity grows with a decreasing ratio of the DC to AC components of the pressure-gradient amplitude below 50%. In particular, for any purely oscillatory pressure-gradient (vanishing DC component), the Womersley solution is a particularly poor predictor of the amplitude and phase of wall shear rate for over half of the flow cycle. Under such circumstances, the analytical models presented here provide a simple and accurate means of estimating instantaneous wall shear rate, knowing only the pressure gradient and hematocrit.     

The effects of erythrocythemia on blood viscosity,maximal systemic oxygen transport capacity and maximal rates of oxygen consumption in an amphibian

Graded erythrocythemia was induced by isovolemic loading of packed red blood cells in the toad, Bufo marinus. Blood viscosity, hematocrit, hemoglobin concentration, maximal aortic blood flow rate and maximal rates of oxygen consumption were determined after each load. Blood viscosity was related to hematocrit in the expected exponential manner; ln eta = 0.43 + 0.035 Hct. Maximal blood flow rates in the dorsal aorta were inversely proportional to blood viscosity and fit predictions of the Poiseuille-Hagen flow formula. The effect of increased blood viscosity was to reduce aortic pulse volume, but not maximal heart rate. Maximal systemic oxygen transport capacity (aortic blood flow rate X hemoglobin concentration X O2 binding capacity of hemoglobin) was linearly correlated with the maximal rate of oxygen consumption. These date indicate that optimal hematocrit theory is applicable for maximal blood flow rates in vivo, and that systemic oxygen transport is the primary limitation to aerial VO2 max in amphibians.     

Effect of plasma exchange on blood viscosity and cerebral blood flow.

The effects of plasma exchange using a low viscosity plasma substitute on blood viscosity and cerebral blood flow were investigated in eight subjects with normal cerebral vasculature. Plasma exchange resulted in significant reductions in plasma viscosity, whole blood viscosity, globulin and fibrinogen concentration without affecting packed cell volume. The reduction in whole blood viscosity was more pronounced at low shear rates suggesting an additional effect on red cell aggregation. Despite the fall in viscosity there was no significant change in cerebral blood flow. The results support the metabolic theory of autoregulation. Although changes in blood viscosity appear not to alter the level of cerebral blood flow under these circumstances, plasma exchange could still be of benefit in the management of acute cerebrovascular disease.     

Dynamics of blood flow: modeling of Fåhraeus and Fåhraeus–Lindqvist effects using a shear-induced red blood cell migration model

Blood flow in micro capillaries of diameter approximately 15–500 μm is accompanied with a lower tube hematocrit level and lower apparent viscosity as the diameter decreases. These effects are termed the Fåhraeus and Fåhraeus–Lindqvist effects, respectively. Both effects are linked to axial accumulation of red blood cells. In the present investigation, we extend previous works using a shear-induced model for the migration of red blood cells and adopt a model for blood viscosity that accounts for the suspending medium viscosity and local hematocrit level. For fully developed hematocrit profiles (i.e., independent of axial location), the diffusion fluxes due to particle collision frequency and viscosity gradients are of equal magnitude and opposite directions. The ratio of the diffusion coefficients for the two fluxes affects both the Fåhraeus and Fåhraeus–Lindqvist effects and is found related to the capillary diameter and discharge hematocrit using a well-known data-fit correlation for apparent blood viscosity. The velocity and hematocrit profiles were determined numerically as functions of radial coordinate, tube diameter, and discharge hematocrit. The velocity profile determined numerically is consistent with the derived analytical expression and the results are in good agreement with published numerical results and experimental data for hematocrit ratio and hematocrit and velocity profiles.     

Asymmetric flows of spherical particles in a cylindrical tube

To study the rheological behavior of blood cells in various flow patterns through narrow vessels, we analyzed numerically the motion of blood cells arranged in one row or two rows in tube flow, at low Reynolds numbers. The particles are assumed to be identical rigid spheres placed periodically along the vessel axis at off-axis positions with equal spacings. The flow field of the suspending fluid in a circular cylindrical tube is analyzed by a finite element method applied to the Stokes equations, and the motion of each particle is simultaneously determined by a force-free and torque-free condition. In both cases of single- and two-file arrangements of the particles, their longitudinal and angular velocities are largely affected by the radial position and the axial spacing between neighboring particles. The apparent viscosity of the asymmetric flows is higher than that of the symmetric flow where particles are located on the tube centerline, and this is more pronounced when particles are placed farther from the tube centerline and when the axial distance between neighboring particles is reduced.     

Plug Effect of Erythrocytes in Capillary Blood Vessels

As an idealized problem of the motion of blood in small capillary blood vessels, the low Reynolds number flow of plasma (a newtonian fluid) in a circular cylindrical tube involving a series of circular disks is studied. It is assumed in this study that the suspended disks are equally spaced along the axis of the tube, and that their centers remain on the axis of the tube and that their faces are perpendicular to the tube axis. The inertial force of the fluid due to the convective acceleration is neglected on the basis of the smallness of the Reynolds number. The solution of the problem is derived for a quasi-steady flow involving infinitesimally thin disks. The numerical calculation is carried out for a set of different combinations of the interdisk distance and the ratio of the disk radius to the tube radius. The ratio of the velocity of the disk to the average velocity of the fluid is calculated. The different rates of transport of red blood cells and of plasma in capillary blood vessels are discussed. The average pressure gradient along the axis of the tube is computed, and the dependence of the effective viscosity of the blood on the hematocrit and the diameter of the capillary vessel is discussed.     

A Mathematical Study on Three Layered Oscillatory Blood Flow Through Stenosed Arteries

A mathematical model is constructed to examine the characteristics of three layered blood flow through the oscillatory cylindrical tube (stenosed arteries).The proposed model basically consists three layers of blood (viscous fluids with different viscosities) named as core layer (red blood cells),intermediate layer (platelets/white blood cells) and peripheral layer (plasma).The analysis was restricted to propagation of small-amplitude harmonic waves,generated due to blood flow whose wave length is larger compared to the radius of the arterial segment.The impacts of viscosity of fluid in peripheral layer and intermediate layer on the interfaces,average flow rate,mechanical efficiency,trapping and reflux are discussed with the help of numerical and computational results.This model is the generalized form of the preceding models.On the basis of present discussion,it is found that the size of intermediate and peripheral layers reduces in expanded region and enhances in contracted region with the increasing viscosity of fluid in peripheral layer,whereas,opposite effect is observed for viscosity of fluid in intermediate layer.Final conclusion is that the average flow rate and mechanical efficiency increase with the increasing viscosity of fluid in both layers,however,the effects of the viscosity of fluid in both layers on trapping and reflux are opposite to each other.     

A couple stress model of blood flow in the microcirculation

A simple mathematical model depicting blood flow in the capillary is developed with an emphasis on the permeability property of the blood vessel based on Starling's hypothesis. In this study the effect of inertia has been neglected in comparison with the viscosity on the basis of the smallness of the Reynolds number of the flow in the capillary. The capillary blood vessel is approximated by a circular cylindrical tube with a permeable wall. The blood is represented by a couple stress fluid. With such an ideal model the velocity and pressure fields are determined. It is shown that an increase in the couple stress parameter increases the resistance to the flow and thereby decreases the volume rate flow. A comparison of the results with those of the Newtonian case has also been made.     

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.     

A problem in the mathematical biophysics of blood circulation: I

The standard approximation method used in mathematical biophysics is applied to the problem of flow of an incompressible viscous fluid in an elastic distensible tube. It is found that the wall of the tube may perform damped transversal harmonic oscillations due to that flow. The phenomenon is independent of the viscosity, the latter contributing only a damping factor. While, due to rather rough approximations, the practical applicability of the equation derived is rather limited, it is suggested that they may give a clue to the understanding of vibrations of the walls of blood vessels which long ago have been suggested as the possible source of some hemic murmurs.     

The influence of suspending phase viscosity on the passage of red blood cells through capillary-size micropores

Much attention has been paid to the study of blood flow in long, narrow tubes. While the influence of tube diameter and driving pressure have been examined in detail, the influence of suspending phase viscosity has generally been assumed only to affect the blood viscosity in a linearly proportional manner, hence the practice of normalizing apparent blood viscosity values by the suspending phase viscosity to give a relative viscosity (e.g., Pries et al., 1992). While this assumption is probably valid for long tubes, it apparently does not hold for blood flow in short tubes (and by extension also for flow in short or branching capillary segments in vivo) in which RBC deformation plays a more significant role. In this paper we present a series of experiments using the Cell Transit Analyzer (CTA) in which the influence of driving pressure and suspending phase viscosity on RBC passage through short, narrow tubes has been systematically evaluated. Over the range studied (1 to 10 cm water), the influence of driving pressure was found to be unremarkable, in that RBC velocity scaled directly and linearly with pressure. This finding is consistent with previous studies. However, a distinct intercept was observed in the linear relationship between RBC pore transit time and suspending phase viscosity, which presumably arises as a consequence of RBC deformation either at the pore entrance or within the pore. Two simple mathematical models for the suspending phase-viscosity/transit-time relationship were considered. The results show that making CTA measurements over a range of suspending medium viscosities is a simple and practical way to obtain additional information about RBC mechanical properties.     

A two-phase model for flow of blood in narrow tubes with increased effective viscosity near the wall

A two-phase model for the flow of blood in narrow tubes is described. The model consists of a central core of suspended erythrocytes and a cell-free layer surrounding the core. It is assumed that the viscosity in the cell-free layer differs from that of plasma as a result of additional dissipation of energy near the wall caused by the red blood cell motion near the cell-free layer. A consistent system of nonlinear equations is solved numerically to estimate: (i) the effective dimensionless viscosity in the cell-free layer (beta), (ii) thickness of the cell-free layer (1-lambda) and (iii) core hematocrit (H(c)). We have taken the variation of apparent viscosity (mu(app)) and tube hematocrit with the tube diameter (D) and the discharge hematocrit (H(D)) from in vitro experimental studies [16]. The thickness of the cell-free layer computed from the model is found to be in agreement with the observations [3,21]. Sensitivity analysis has been carried out to study the behavior of the parameters 1-lambda, beta, H(c), B (bluntness of the velocity profile) and mu(app) with the variation of D and H(D).