PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
J Biomech Eng. Author manuscript; available in PMC 2010 May 11.
Published in final edited form as:
PMCID: PMC2867477
NIHMSID: NIHMS197376

Blood Flow in Small Curved Tubes

Abstract

Blood flow in small curved tubes is modeled by the two-fluid model where a relatively cell-free fluid layer envelops a fluid core of higher viscosity. The parameters in the model are successfully curve fitted to experimental data for straight tubes. The curved tube equations are then solved by perturbation theory. It was found that curvature in general lowers the tube resistance, but increases the shear stress near the inside wall.

1 Introduction

Blood can be considered as homogeneous in large blood vessels, but in small tubes the rheology shows considerable non-homogeneous behavior [1,2]. Due to the particulate nature of red blood cells, there is an almost cell-free layer of plasma near the vessel wall. The major effect of this layer is to decrease the apparent viscosity or resistance (the Fahraeus-Lindqvist effect), especially in medium to moderately small blood vessels. For vessels of 50-1000 microns, the phenomenon can be adequately described by the two-fluid model. First proposed by Vand [3] for general suspensions, the model assumes a homogeneous fluid core enclosed by a fluid annulus of lower viscosity. The two-fluid model has been applied to blood flow in straight tubes, notably by Hayes [4] and Sharan and Popel [5]. The theoretical derivation of the resistance formula of this model and our simple curve fit with blood flow data in small straight tubes are given in the Appendix.

In the microvasculature there are numerous instances where the blood vessels are curved or tortuous [6], especially in the diseased state. Also in biomechanical instrumentation blood may be transported through small curved tubes. It is the aim of this study to investigate the effect of curvature on the wall shear and the resistance of the flow through such a tube.

We shall use a coordinate system first proposed by Dean [7,8] to treat the flow in a circular tube of small constant curvature. The inertial effects are neglected since the Reynolds number in the microcirculation is typically very small (of order 10−3). The Dean number, proportional to (curvature)(Reynolds number)2, is even smaller. The secondary recirculation, which is prevalent for large curved tubes, would be entirely absent.

Formulation

Let the tube be of radius a, with a constant centerline curvature K. An orthogonal system can then be constructed from the elemental distance squared (Fig. 1a):

Fig. 1
(a) The curved tube coordinates (r,θ,s). The direction s is along the center line. (b) Axial cross section showing two fluid regions and the blunted velocity profile.
|dx|2=(dr)2+r2(dθ)2+L2(ds)2
(1)

where the center line scale factor is

L1Krcosθ
(2)

(r,θ) are local polar coordinates and s is the coordinate along the center line. Since the inertial terms are absent, due to symmetry one can show only the velocity component w(r,θ) in the s direction exists. Using the scale factors from Eq. (1) (e.g. Batchelor [9]) the Stokes equation reduces to:

μr{[rL(Lw)r]r+[1rL(Lw)θ]θ}=GL
(3)

Here G is the constant center-line pressure gradient (assumed negative along s), and μ is the viscosity. Fig. 1b shows the two regions. Region I is the core containing red blood cells with velocity wI and viscosity μI and Region II is the cell-free layer of thickness t from the vessel wall and with velocity wII and viscosity μII. Eq. (3), with the same pressure gradient, is solved separately for the two regions. The boundary conditions are that on the axis (at r = 0) wI is bounded, on the wall (at r = a) there is no slip, i.e. wII(a,θ) = 0 and velocities and shear stresses match on the interface at r = a-t:

wI=wII
(4)
μIL(wIL)r=μIIL(wIIL)r
(5)

Solution

The solution for the two-fluid flow in a straight circular tube is given in the Appendix. We shall perturb from this solution. Let the ratio of the tube radius to the radius of curvature of the center line be small, i.e.:

aK=ε
(6)

where ε[double less-than sign]1. (Actual values of ε is between 0 and 0.3 obtained from Ref. [6]). The other variable are normalized to order unity as follows:

ηr/a
(7)
wI=U(w0+εw1+ε2w2+)
(8)
wII=αU(u0+εu1+ε2u2+)
(9)

Here U = |G|a2/4μI is a normalized pressure gradient, α [equivalent] μIII, and w0,u0 is the normalized straight tube solution (Eqs. (A2-A4)):

w0=Cη2,u0=1η2
(10)

where C is defined in Eq. (A4). Upon substitution of Eqs. (8,9) into Eqs. (3) the first-order equation is found to be

4ηcosθ+1η[η(w1ηw0cosθ)η+η2cosθw0η]η+1η2(w1ηcosθw0)θθ=0
(11)

The first-order equation for wII is similar to Eq. (11), only with w replaced by u. Together with the perturbed boundary conditions Eqs. (4,5) the solution is:

w1=cosθ[C1η34η3]
(12)
u1=cosθ[34(ηη3)+C2(η1η)]
(13)

where δ[equivalent]t/a and

C1={6+(α1)[18δ(16α+21)δ2+(4α+3)(4δ3δ4)]}4[2+2(α1)δ(α1)δ2]
(14)
C2=7(α1)δ(2δ)(1δ)24[2+2(α1)δ(α1)δ2]
(15)

Notice the outer region velocity u1 is zero on the tube wall but may include terms which are singular at the center. Since w1,u1 are periodic in θ, they do not contribute to the net flow rate, which is of higher order. However, for the second-order flow rate correction only the non-periodic part is needed.

Let an over bar denote the average with respect to the angle θ. The averaging is then applied to the second-order terms of Eq. (3). The result is:

2η3[η(12ηw0+cosθw1¯w2η¯)]η=0
(16)

with a similar equation for u2. After some work, the solution is:

w2¯=C3+C+C14η21132η4
(17)
u2¯=7+4C216(1η2)+1132(1η4)+C4lnη
(18)

where

C3=132[2+2(α1)δ(α1)δ2]×{2+(α1)[10δ+(8α21)δ28(2α3)δ3+2(7α+8)δ46(α1)δ5+(α1)δ628αδ(2δ)(1δ)2ln(1δ)]}
(19)
C4=21(α1)(1δ)2δ(2δ)8[2+2(α1)δ(α1)δ2]
(20)

The net flow rate is then integrated

F=2πa2U[q0+ε2q2+O(ε4)]
(21)

where

q0=01δw0ηdη+α1δ1u0ηdη=α4[1(1δ)4(11α)]
(22)
q2=01δw2¯ηdη+α1δ1u2¯ηdη=1192[2+2(α1)δ(α1)δ2]×{2+(α1)[14δ(576α+43)δ2+4(432α+19)δ3(1856α+85)δ4+(832α+62)δ5(120α+29)δ68(α1)δ7+(α1)δ8]}
(23)

The above computations are facilitated by a computer program with symbolic capabilities

Results and Discussions

Figure 2 shows some typical velocity distributions. The zeroth-order is a blunted parabola form which has been observed in experiments using straight tubes [2]. The first-order correction is due to the curving of the vessel. Contrary to high Reynolds number flows, the velocity near the inner surface is increased while that near the outer surface is decreased. Consequently the shear stress is higher on the inside wall near point A in Fig. 1a. The shear stress on the vessel wall is:

Fig. 2
Typical velocity profiles (δ=0.1,α=3.2,θ=0). The first-order correction w1 due to curvature is anti-symmetric.
τ=μIIwIIr=|G|a4(u0η+εu1η+)|η=1=|G|a4[2+εcosθ(32+7(α1)δ(2δ)(1δ)22[2+2(α1)δ(α1)δ2])]
(24)

Whether higher local shear stress is the cause of vessel tortuosity is still uncertain at present. Fig. 3a shows a typical zeroth-order net flow rate q0 versus the thickness-radius ratio. For smaller vessels (larger δ) the flow rate can be as much as 50% more due to the lowered viscosity of the cell-free layer. Fig. 3b shows a typical first-order flow rate correction q2 due to curvature effects. Notice q2 becomes negative when δ is larger than 0.026. This means the resistance of a larger curved vessel is lower than that of a straight vessel, but the resistance of a smaller curved vessel is higher than a straight vessel of the same size. These properties would persist even when the curvature is not small.

Fig. 3
The flow rate as a function of δ. (α=3.317). q0 is for the straight tube, q2 is the correction due to curvature.

We can also solve Eq. (3) numerically. This is not recommended since the parameter ε is indeed small and the perturbation method is adequate. Furthermore perturbation solutions, being exact, do show more clearly the dependence of the various parameters.

Acknowledgments

This research is partially supported by NIH Grant RR 01243.

Appendix-Two-Fluid Flow in a Straight Tube

For a straight tube the curvature K is zero and thus the scale factor L = 1. Eq. (3) becomes the Poisson equation:

μr[(rwr)r+1rwθθ]=G
(A1)

Let δ[equivalent]t/a and α[equivalent]μIII. The axisymmetric solution that satisfies the conditions of boundedness, no slip, and Eqs. (4,5) is:

wI=|G|a24μI[C(r/a)2]
(A2)
wII=|G|a24μII[1(r/a)2]
(A3)

where

C=(1δ)2α[(1δ)21]
(A4)

The total flow can be integrated

Q=0atwI2πrdr+atawII2πrdr=π|G|a48μII[1(1δ)4(11α)]
(A5)

In order to compare with experiments, we define an apparent viscosity

μapp=π|G|D4128Q
(A6)

where D = 2a is the inside tube diameter, and a relative apparent viscosity

μrelμappμII
(A7)

Using Eq. (A5) we find

μrel=11(1δ)4(11/α)
(A8)

This is essentially the result of Vand [3]. The apparent viscosity is that of the core when the thickness of the cell-free layer is zero, and is the plasma viscosity when the thickness equal the radius of the tube. Using in vitro experiments (in glass tubes) of 19 sources compiled by Pries et al. [1] and Secomb [10] we curve-fitted Eq. (A8) for each hematocrit as shown in Fig. 4. The fitted values of α and t = D δ/2 are shown in Table 1.

Fig. 4
Relative apparent viscosity as a function of diameter ——glass tube experiments [10], — — — Eq. (A8)
Table 1
Values of α and t fitted to experimental data in Fig. 4

The error of the curve fit is within 3% for diameter D from 40 to 1000 microns. We keep in mind that any shear dependence from rouleaux formation, cell deformation, rotation, wall effects, nonparticle free plasma layer, etc are included in the apparent curve fit of the two parameters. Also included is the relation between the tube hematicrit and the feed hematcrit, upon which the data was based. Thus it is not necessary to consider such complicated effects for the two-fluid model. For in vivo experiments with blood vessels instead of glass tube, the experiments show higher apparent viscosity, probably due to the uneven glycocalyx wall. The curves however can be fitted similarly.

Contributor Information

C. Y. Wang, Departments of Mathematics and Physiology, Michigan State University, East Lansing, MI 48824.

J. B. Bassingthwaighte, Center for Bioengineering, University of Washington, Seattle, WA 98195.

References

1. Pries AR, Secomb TW, Gaehtgens P. Biophysical Aspects of Blood Flow in the Microvasculature. Cardiovasc Res. 1966;32:654–667. [PubMed]
2. Bishop JJ, Popel AS, Intaglietta M, Johnson PC. Rheological Effects of Red Blood Cell Aggregation in the Venous Network: A Review of Recent Studies. Biorheology. 2001;38:263–274. [PubMed]
3. Vand V. Viscosity of Solutions and Suspensions I Theory. J Phys Chem-US. 1948;52:277–299. [PubMed]
4. Hayes RH. Physical Basis of the Dependence of Blood Viscosity on Tube Radius. Am J Physiol. 1960;198:1193–1200. [PubMed]
5. Sharan M, Popel AS. A Two-Phase Model for Flow of Blood in Narrow Tubes With Increased Effective Viscosity Near the Wall. Biorheology. 2001;38:415–428. [PubMed]
6. Nicola PD, Maggi G, Tassi G. Microcirculation-An Atlas. Schattauer; Stuttgart: 1983.
7. Dean WR. Note on the Motion of a Fluid in a Curved Pipe. Phil Mag Ser. 1927;74:208–233.
8. Berger SA, Talbot L, Yao LS. Flow in Curved Pipes. Annu Rev Fluid Mech. 1983;15:461–512.
9. Batchelor GK. An Introduction to Fluid Dynamics. Cambridge Univ Press; 1967. Appendix 2.
10. Secomb TW. Mechanics of Blood Flow in the Microcirculation. In: Ellington CP, Pedley TJ, editors. Biological Fluid Dynamics. Soc. Exp. Biol.; Cambridge, MA: 1995. pp. 305–321.