Chapter 4: Flow
In vivo, transportation of materials is controlled by the flow of a fluid such as liquid
or gas.

When the control of this flow does not work, life activity stagnates.

In thischapter, we learn fundamental matters of flow dynamics, such as resistance of bloodvessels.

Technological analysis of the flow phenomena of the biological fluid in vivo

is necessary for the design of the artificial heart and artificial blood vessel cooperatingwith a living body.

4.1 Fluid and solid
4.1.1 Fluid and pressure
In the stationary fluid, the force is uniformly transmitted to every direction.
vertical force at unit plane is called the “pressure”.

The

In each point of the stationary

fluid, the force at unit plane is constant regardless of the direction of the plane.

The“stress” in the solid (see 3.1.4), on the other hand, depends on the direction of the plane(Fig. 4.1).

Fig. 4.1: Pressure and stress.

A connecting tube, which enables the movement of the liquid between two positions© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

is called “communicating tube” (Fig. 4.2).

The pressure transmission through thestationary liquid in the communicating tube is applied to the remote measurement of thepressure.

The tube called a “catheter” is inserted into a blood vessel, filled with asaline solution (see 5.3.2), and applied to telemetry of the intravascular pressure.

Fig. 4.2: Telemetry of pressure by the communicating tube.

The mass per unit volume is called “density”.
uniform density in equilibrium.
the pressure.

The unit is kg m－3.

The fluid has

The volume of the “compressible fluid” is reduced by

The volume of the “incompressible fluid” is not reduced by thepressure.
The compressibility of the liquid is very small compared with that of the gas.

Thedensity of the incompressible fluid is constant regardless of the pressure.
When the density is constant, the law of conservation of mass leads to the law ofconservation of volume.

The law is described as the continuity equation 4.1 (Fig.

4.3).

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.3: Continuity.
A1 v1＝A2 v2(4.1)

In the equation (4.1), A is the flow passage cross-sectional area [m2], v is the flow
velocity [m s－1].

The equation shows that the flow velocity (v) increases, when theflow path cross-sectional area (A) decreases in the pipe.
The equation of Bernoulli (4.2) shows the law of conservation of mechanicalenergy in the form per unit volume.
(1/2) ρ v 2＋p ＝constant(4.2)

In equation (4.2), ρ is the density [kg m－3], v is the flow rate [m s－1], and p is thepressure [Pa]. The pressure (p) decreases at the fast flow velocity (v), when ρ isconstant in Eq. 4.2 (Fig. 4.4).

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.4: Expression of Bernoulli.

Fig. 4.5: Head drop.

In the gravitational field, fluid pressure p [Pa] is generated as the product of thedensity ρ [kg m－3], the gravitational acceleration g [m s－2] and the height h [m] (Fig.
4.5).
p＝ρ g hThe pressure by the head is proportional to the density.

(4.3)

The pressure by the head is© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japandifferent between the fluids of different density.

The head control valve detects thefluid pressure difference between the fluids [3, 20].

The device can be applied to theflow control in the shunt (the connecting tube between the intracranial and the atrium
or the peritoneal cavity) for hydrocephalus (Fig. 4.6).

Fig. 4.6: Shunt.

The blood is drained according to the head with the communicating tube.

The

flow can pass over the higher position through the communicating tube by “theprinciple of a siphon” (Fig. 4.7).

Fig. 4.7: Principle of siphon.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, JapanAt the room temperature (298 K) and at 1 atm (101.325 kPa), consider the
movement of water in accordance with the principles of a siphon (see the equation(4.3)).

In the gravitational field, the water of 10 m height generates the followingpressure:
1.0×103 kg m－3×9.8 m s－2×10 m＝98 kPa(4.4)

The water pressure at the point elevated by 10 m under 1 atm (101.325 kPa) is

101.325 kPa－98 kPa＝3.325 kPa(4.5)

The vapor pressure of water is 3.140 kPa at 298 K.
At the point higher than 10 m at 298 K, the pressure of the water is below the vapor
pressure, in which the water is vaporized.

The water cannot be drained over the heightof 10 m with the communicating tube.
The venous blood vessel wall has large compliance.

When the amount of blood in

the local vein decreases, the lumen becomes smaller.

Successively, the walls of thevessel stick together, and the lumen is embolized.
(a)).

This is called “collapse” (Fig. 4.8

The resistance of the flow path increases with the collapse.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.8: Collapse.

The collapse should be avoided at the drainage of blood from the venous system.
Several venous cannula tips are designed to avoid the collapse (Fig. 4.8 (b)) [21].

4.1.2 Elasticity and viscosity
When you hang a weight on a coil spring, the coil spring stretches.

When youremove the weight on the coil spring, the coil spring returns to its original length.

The“Hook elastic body” is the deformation model for the object similar to the coil spring(Fig. 4.9).

The solid is approximated to the Hook elastic body.

In the Hook elastic

body, the "strain" is proportional to the “stress” (see 3.1.5).

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.9: Hook elastic body.

The water does not return to its original shape and position after removal of theforce on the water.

Although a small force is enough to make the slow flow, a large

force is necessary to make the rapid flow.
Although the fluid such as liquid or gas is sheared by a force, it does not return tothe original shape and position after removal of the force.

In the fluid, the large force

increases the "deformation rate".
The change of the strain per the unit of time is called the “strain rate”.

The unitsof strain rate is s－1.
The deformation rate is quantitatively expressed by the rate at the shear between thesurfaces.

The shear rate (γ) is the quotient (Eq. 4.6) of the speed difference between

the surfaces (Δv) divided by the distance (y) between the surfaces (Fig. 4.10).

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.10: Shear rate.
γ＝Δv／y(4.6)

Unit of the shear rate is s－1.

For example, the speed difference of 360 km per hour

between two plates with the interval of 1 m makes the shear rate of 100 s－1
The shear force per the unit area is called the shear stress.
stress is Pa.

The unit of the shear

The quotient of the shear stress τ divided by the shear rate γ is called thecoefficient of the viscosity η (Eq. 4.7).
η＝τ／γThe unit of the coefficient of the viscosity is Pa s.

(4.7)

The unit of poise [P] issometimes used for the viscosity (Eq. 4.8).
1 P＝1 dyn cm－2 s＝0.1 N m－2 s＝0.1 Pa s

(4.8)

In the “Newtonian fluid”, the "shear stress" is proportional to the "shear rate" (Fig.
4.11). The Newtonian fluid is used as the simple model for the flow of the liquid or thegas.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.11: Newtonian fluid.

The coefficient of the viscosity represents the flow resistance of the fluid, anddepends on the temperature.

In the liquid, the coefficient of the viscosity decreases

with the increase of the temperature (Fig. 4.12).

The coefficients of the viscosity of

the water at 293 K and at 313 K are 1.002 × 10－3 Pa s and 0.653 × 10－3 Pa s,
respectively.

In the gas, conversely, the coefficient of the viscosity increases with theincrease of the temperature. The coefficient of the viscosity of the gas is much smaller
than the coefficient of the viscosity of the liquid.

In the air, the coefficients of theviscosity at 298 K and at 323 K are 18.2 × 10－6 Pa s and 19.3 × 10－6 Pa s, respectively[22].

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.12: Viscosity with temperature.

The viscosity can be explained by the interaction among the particles, which moveas units in the fluid.

The vibration of each particle at higher temperature induces shear

flow between particles, which decreases the viscosity.
For example, the oil becomes smoother with the higher temperature, and thicker
with the lower temperature.

In the extracorporeal circulation, the viscosity of theblood increases at the lower temperature.

The increase in viscosity is reduced bydiluting the blood with the plasma substitute.
The coefficient of the viscosity of the gas, on the other hand, increases with theincrease of the temperature. The velocity of the molecule of the gas increases with theincrease of the temperature, which increases the interaction such as collision among themolecules.

The interaction increases the coefficient of the viscosity.

In therespiration, the lower temperature reduces the coefficient of the viscosity of the air,which reduces the flow resistance in the respiratory tract.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

The blood has the higher coefficient of the viscosity at the lower shear rate.
shear stress is not proportional to the shear rate in the blood.

The

The fluid, in which thecoefficient of the viscosity varies with the shear rate, is called non-Newtonian fluid(Fig. 4.13).

Fig. 4.13: Viscosity with shear rate.

The blood contains a huge number of red blood cells.
blood cells is called hematocrit Ht.

The volume ratio of red

Ht is approximately equals to 40% in human.

The value is a large number, which is the same level with the volume ratio of the hard
spheres in the face-centered cubic lattice (see Question 3.3).
When the flow is slow, the interaction between the red blood cells increases the flow
resistance. The interaction among the red blood cells decreases at the fast flow, which
reduces the resistance of the flow.

In the fast flow, the deformation of each erythrocyte

is also contributes to the reduction of the flow resistance of the blood.
In the low shear region, red blood cells accumulate like a stack of coins, which iscalled rouleau (Fig. 4.14). The coefficient of the viscosity increases at the higher
hematocrit.

The viscosity of the blood changes with the thrombus formation. Thechange can be detected by a vibrating electrode (see 2.2.2), (Fig. 4.15) [23].

By the© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japanvibrating electrode, the differences can be detected not only on the impedance, but also
on the viscosity of the yolk and the albumen (Fig. 4.16) [7].

Fig. 4.14: Rouleau formation.

Fig. 4.15: Viscosity tracings with vibrating electrode.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.16: Measurement of local viscosity with vibrating electrode.

4.1.3 Viscoelasticity
The most of objects cannot deform instantaneously, and can restore the force.
They have the multi-properties of viscosity and elasticity.
viscoelasticity.
body.

This multi-property is called

The object, which shows viscoelastic property, is called viscoelastic

The deformation behaviors of the polymeric materials and the biological tissues

can be explained by viscoelasticity.
In Maxwell model, a viscous element is connected to an elastic element in series(Fig. 4.17). When a strain is applied to the model with a step function, a stress isgenerated in response to the strain at the elastic element.

The stress governs thedeformation speed at the viscous element.

The deformation at the viscous element

reduces the deformation at the elastic element.

The reduction of the deformation at theelastic element decreases the stress with time.

This phenomenon is called stressrelaxation.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.17: Maxwell model.

In Kelvin-Voigt model, a viscous element is connected with an elastic element inparallel (Fig. 4.18). When a stress is applied with a step function, the deformation
starts at the viscous element with a speed corresponding to the stress.

The gradualincrease of the strain at the elastic element increases the stress at the elastic element.
The amount of increase of the stress at the elastic element reduces the amount of thestress at the viscous element.

The reduction of the stress at the viscous element

reduces the deformation speed at the viscous element with time. When all of the stressis supported by the elastic element, the deformation of the elastic element is saturated.
The increase of deformation with time is called creep deformation.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.18: Kelvin-Voigt model.

A polymer solution has elasticity as well as viscosity.

Therefore, the solution

spreads when it is released from the narrow flow path (Pallas effect).

The surface ofthe solution near the rotating bar is lifted by the stress generated at the direction
different from that of the shear stress (Weissenberg effect).

The vortex and theturbulent flow resistance decrease (Toms effect) in the fluid.

4.2 Resistance of flow and distribution of velocity
4.2.1 Resistance of flow
Consider the resistance when the blood flows through a blood vessel.
volume per unit time is called “flow rate”.

The flowing

The flow rate Q [m3 s-1] increases inproportion to the pressure difference ΔP [Pa] between the upstream and the downstream.
In this case, the proportionality constant is the flow resistance Rf [Pa m-3 s].

Rf＝ΔP／Q(4.9)

This relationship corresponds to "an electrical resistance R when an electric current© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japanflows through a conductive wire" calculated as a value of the "the potential difference
ΔE between both ends of the conductive wire" divided by “the electric current I”.
R＝ΔE／I(4.10)“Systemic circulation resistance (total peripheral resistance) Rs” is a quotient of thedifference between “aortic pressure Pa” and “right atrial pressure Pr” divided by"cardiac output Qc" (Fig. 4.19).

Fig. 4.19: Circulation resistance.

Rs＝(Pa－Pr)／Qc(4.11)"Pulmonary circulation resistance (Rp)" is a quotient of the difference between"pulmonary arterial pressure (Pp)" and "left atrial pressure (Pl)" divided by "cardiac
output (Qc)".

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Rp＝(Pp－Pl)／Qc(4.12)

Since the pulmonary circulation has only one organ “the lung”, the resistance of thepulmonary circulation is lower than that of the systemic circulation: about one-fifth.
Since "cardiac output" is the common at both left and right ventricles, the pressure
difference at the pulmonary circulation is about one-fifth of that at the systemic
circulation.

When you count the bronchial circulation (from the aorta to thepulmonary vein), the output of the left ventricle is slightly bigger than that of the rightventricle.
Since the pressure difference [Pa] is divided by the flow rate [m3 s-1], the unit offlow resistance is [Pa m-3 s]. When [Pa] is replaced by the SI base unit [kg m-1 s-2], theunit of flow resistance becomes [kg m-4 s-1].
The aortic pressure of 13 kPa, the right atrial pressure of 1 kPa, and the cardiac
output of 10-4 m3 s-1 (= 6 l min-1) make the systemic circulation resistance of 12 × 107
Pa m-3 s.

With the similar equation, the pulmonary artery pressure of 3 kPa, the left

atrial pressure of 1 kPa, and the cardiac output of 10-4 m3 s-1 (= 6 l min-1) make thepulmonary circulation resistance of 2 × 107 Pa m-3 s.

In the pulsatile flow (see 4.3.1),

both the pressure and the flow rate vary periodically according to the pulsation cycle.
During discharge, the high pressure ejects the fluid.

During suction at theatmospheric pressure, on the other hand, 1 atm (101.325 kPa) is the maximum value forthe pressure difference.
the suction.

For this reason, it is not easy to increase the flow rate during

In the pumping of the heart, diastole is longer than systole (see question

4.6).
Electrical resistance R [Ω] of the metal wire is proportional to the length l [m], andis inversely proportional to the cross-sectional area A [m2].

R＝ρl／A(4.13)© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

The resistivity ρ [Ω m] is used to compare the magnitude of the resistance of metal.
The resistance of the flow increases with the length of the tube.

The resistance ofthe flow decreases with the increase of the cross-sectional area of the tube.

Is theresistance of the flow inversely proportional to the cross-sectional area (see 4.2.2)?

4.2.2 Hagen-Poiseuille flow
A flow, which has same velocity over the entire cross section of the circle pipe, iscalled “Plug flow”.

In the flow, the flow velocity in the center of the pipe is same asthat of the vicinity of the pipe wall (Fig. 4.20 (a)).

Fig. 4.20: Velocity distribution in pipe.

The flow velocity has the distribution at the cross section of the pipe, where the flow
velocity of the central position is faster than that of the vicinity of the wall.

A flow,

which has the parabolic velocity distribution at the cross section of the circle pipe, iscalled “Hagen-Poiseuille flow”.

In Hagen-Poiseuille flow, the velocity vector is asfollows: the velocity is the maximum at the center axis, and zero at the surface of thevessel wall.

The velocity vectors are parallel and symmetric with respect to the axis of

the pipe.

The envelope of the velocity vectors makes parabola, where the tip© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japancorresponds to the maximum value (Fig. 4.20 (b)).
In the following equations, the flow velocity distribution in the Hagen-Poiseuilleflow is introduced from a mechanical balance.
Let set concentric thin cylinders (length l, velocity v, radius r, and thickness dr) in
the circular pipe (radius a, and length l) (Fig. 4.21).
faster than that of the outer cylinder (Fig. 4.22).
slip between cylinders.

The speed of the inner cylinder is

The difference of the speed makes the

In the viscous flow, the slip makes the frictional force: withthe same direction to the flow by the inner cylinder, and with the opposite direction to
the flow by the outer cylinder, respectively.

When the cylinder makes a uniform linear

motion at a velocity v, the total force applied at the cylinder is zero.

Therefore, “theforce due to the pressure difference (ΔP) between the upstream side and the downstreamside” balances with “the frictional force F between the cylinders”.

F is the product of

the lateral area of the cylinder (2 π r l) and the shear stress τ.

Fig. 4.21: Force balance in cylinder in flow.
F＝2 π r l τ(4.14)

τ＝－η dv／dr(4.15)© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, JapanIn Eq. 4.15, η is the coefficient of viscosity, and dv / dr is the shear rate.
of r is the center of the circular pipe.

The origin

The velocity v has the maximum value at thecenter (r = 0). The velocity v decreases with the increase of r (Fig. 4.22).
reason, Eq. 4.15 has a negative sign at the right side.

For thisEq. 4.16 is introduced bysubstituting Eq. 4.15 into Eq. 4.14.

Fig. 4.22: Cylinders of fluid in flow through pipe.
F＝－2 π l η r dv／dr(4.16)

The frictional force dF, which is the difference between the frictional force withouter cylinder and that with the inner cylinder, resists the flow.
dF＝－2 π l η(1／dr)(r dv／dr) dr(4.17)

The force dF balances with the force derived from the pressure difference ΔP
between the upstream and the downstream, which are applied on the end surface area (2π r dr) of the thin cylinder.
dF＝2 π r dr ΔP(4.18)

Eq. (4.19) is derived from Eq. (4.17) and Eq. (4.18).

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

－2 π l η(1／dr)(r dv／dr) dr＝2 π r dr ΔP(4.19)

Eq. (4.19) is rewritten to Eq. (4.20).
－(1／dr)(r dv／dr) dr＝(ΔP／(l η))r dr(4.20)

Both sides of Eq. (4.20) are integrated in r.
－∫(1／dr)(r dv／dr) dr＝(ΔP／(l η)) ∫r dr(4.21)－r dv／dr＝(ΔP／(l η))(r2／2)＋C(4.22)In Eq. (4.22), C is an integration constant.
at r = 0 (center of the pipe).

C should be zero to apply the Eq. (4.22)

Thus, the Eq. (4.22) becomes Eq. (4.23).

r dv／dr＝－(ΔP／(l η))(r2／2)(4.23)The both sides of Eq. (4.23) are divided by r.
dv／dr＝－(ΔP／(l η))(r／2)(4.24)The left side of Eq. (4.24) is the tangent inclination of the envelope of the velocity
vectors.

dv／dr represents the shear rate γ.

γ is zero at the central axis, and the

maximum at the wall of the pipe.
Shear rate at the wall (γw) is calculated at r = a, where a is the inner radius of thepipe.
γw＝－(ΔP／(l η))(a／2)(4.25)© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

The both sides of the Eq. (4.24) are integrated by r.
∫dv＝－(ΔP／(l η))∫(r／2)dr(4.26)

v＝－(ΔP／(4 l η))r2＋D(4.27)

Equation (4.27) represents the flow velocity distribution in the cross section of thepipe.

D is an integration constant.

Since the flow velocity is zero at the wall, v = 0 at

r = a.
D＝(ΔP／(4 l η))a2(4.28)

Eq. (4.27) combined with Eq. (4.28) makes Eq. (4.29).
v ＝(ΔP／(4 l η))(a2－r2)(4.29)

Eq. (4.29) represents a parabolic flow velocity distribution.

The flow rate is

calculated by integrating the flow velocity in the cross-section of the channel.
𝑎

𝑄 = � 2𝜋𝜋 𝑑𝑑0𝛥𝛥 2(𝑎 − 𝑟 2 )4𝑙𝜂

𝜋 𝛥𝛥 𝑎 2𝑄=� (𝑎 𝑟 − 𝑟 3 )𝑑𝑑2𝑙𝜂 0𝑎

𝜋 𝛥𝛥 𝑎2 𝑟 2 𝑟 4
𝑄=
�
− �
2𝑙𝜂
2
4 0

Q ＝πa4ΔP／(8 l η)

(4.30)(4.31)(4.32)(4.33)

Eq. (4.33) is rewritten to Eq. (4.34).
ΔP ＝8 l η Q／(πa4)(4.34)© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, JapanEq. (4.34) is substituted into Eq. (4.25).
Γw ＝－4 Q／(πa3)(4.35)

Eq. (4.35) shows that the wall shear stress (Γw) is inversely proportional to the cubeof the radius at the constant flow rate.

By substituting Eq. (4.34) into Eq. (4.9),
Rf＝(8 l η)／(πa4)(4.36)In this distribution of flow velocities, resistance Rf [kg m-4 s-1] of flow in the straight
circular pipe is inversely proportional to the fourth power of the radius a [m].

When

the radius is reduced by 16%, the resistance of the flow is doubled.
Can the flow rate is intuitively understandable with the Illustration of the flow
velocity distribution?
Fig. 4.23.

The inner radius of the circular pipe is doubled from (a) to (b) in

At the same flow velocity distribution, the four times bigger sectional area

makes the four times bigger flow rate in “the plug flow” (Fig. 4.23 (a)).

In “theHagen-Poiseuille flow”, on the other hand, the increased flow velocity at the center
makes 16 times increase of the flow rate (Fig. 4.23 (b)).

In Fig. 4.23 (b), the whole

velocity distribution is extrapolated from the velocity distribution in the vicinity of thewall.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.23: Distribution of velocity.

When the flow velocity distribution of the Hagen-Poiseuille flow is applied to thehuman blood flow in the vessel with the inside diameter in each section (from 7 ×10－6
m to 0.025 m), the shear rate at the vascular wall is estimated by Eq. (4.35) as 60 ~ 800
s－1 [24].

The shear flow is one of the methods for continuous application of the mechanicalstimulus to the cells. Vascular endothelial cells cover the inner surface of the vessel
wall.

Under stimulation of the shear flow at the wall, orientation has been observed sothat the long axis of each cell makes orientation along the stream line.

On the other

hand, orientation of myoblasts has been observed perpendicular to the flow in theprocess of forming the myotubes by differentiation [16].

4.2.3 Requirement for Hagen-Poiseuille Flow (boundary conditions)“Hagen-Poiseuille flow” can be applied to a steady laminar flow in a sufficiently
long straight circular pipe without any branches (see 4.3.2).

The velocity distribution© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japanin Fig. 4.20(b) is realized, when viscous forces between the layers are equilibrium in theentire flow in the circular pipe.
Vessels have bends and branches.

The blood flow is pulsatile in the artery.

Sinceblood is a non-Newtonian fluid, the viscosity varies with the shear rate.

Thus, the flow

velocity distribution deviates from that of the Hagen-Poiseuille flow.

Even if above

conditions were considered, the flow resistance drastically increases with the reductionof the inner diameter of the vessel.
The flow velocity distribution in the cross-section shifts from Fig. 4.20 (b) in thevicinity of the inlet of the pipe, because the flow velocity distribution of the upstream istaken over.

At the downstream, the velocity distribution approaches to that of Fig.

4.20 (b).

The section is called as “inlet region”, and the length is called as “inlet

length” (Fig. 4.24).

The “inlet length” becomes shorter, when the radius r of thecircular pipe or the Reynolds number Re (see 4.3.2) is small.

Fig. 4.24: Inlet region.

In the viscous flow, the flow velocity profile becomes parabolic in the vicinity of thewall surface. The area is called as the “boundary layer” (Fig. 4.25).
The viscosity coefficient of the fluid can be calculated from the flow resistance atthe Hagen-Poiseuille flow through the capillary.

Measurement of the viscosity by the© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japancapillary, however, is difficult on the blood of non-Newtonian fluid, because the flow
velocity distribution deviates from the Hagen-Poiseuille flow.

The viscosity of

non-Newtonian fluid should be measured in the uniform shear field (see 4.2.4).

Fig. 4.25: Boundary layer.

4.2.4 Couette flow
Consider the flow of the fluid sandwiched between two walls with the distance of d.
One of the walls is moving at speed of v, and the other wall is stationary.
the wall moves at the same velocity of that of the wall.
stops.

The fluid at

The fluid at the stationary wall

The fluid sandwiched between the walls flows in the intermediate speed.

When the velocity at each position increases proportionally to the distance from thestationary wall, the flow is called “Couette flow” (Fig. 4.26).

In this case, theenvelope of the velocity vector of the flow makes a straight line.

Fig. 4.26: Couette flow.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, JapanIn the Couette flow, the shear rate γ is constant regardless of the distance from the
wall.
γ＝v／d

(4.37)In the fluid, which is sandwiched between the rotating cone and the stationary plate,the flow velocity distribution of the Couette flow type occurs (Fig. 4.27).

The distanced between the cone and plate increases in proportion to the distance r from the axis ofrotation.

When θ is very small,Fig. 4.27: Flow between rotating cone and stationary plate.

d＝r tanθ＝rθ(4.38)

θ is the angle [rad] between cone and plate.

The speed v [m s－1] of the conicalsurface is proportional to the distance r [m] from the axis of rotation.
v＝r ω
ω is the angular velocity [rad s－1].

(4.39)The shear rate γ is constant regardless of the© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japandistance r from the axis of rotation.
γ＝v／d＝r ω／(r θ)＝ω／θ

(4.40)A uniform shear rate, which is suitable for the measurement of non-Newtonian
fluids, is applied to the entire fluid (see 4.1.2).

In the cone-and-plate viscometer, thefluid is sandwiched between the rotating cone and the stationary plate (Fig. 4.28).

In

the viscometer, moment (torque) to keep the rotation of the cone at a constant speed is
proportional to the viscosity coefficient of the fluid (see Q. 4.9).

Fig. 4.28: Cone-plate viscometer.

The cone and plate viscometer can be applied to evaluate the effect of shear rate on
clot formation, and to measure increase of the flow resistance with clot formation (Fig.
4.29 (a)). The clotting time can be measured by the time t0 before the rise of the torque.
When the torque is increased from T0 to T1, the thrombus formation ability can beevaluated by the increase rate Rc (thrombus ratio) calculated by Eq. (4.41) (Fig. 4.29(b)) [25].

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.29: Clotting between rotating cone and stationary plate.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.29: Clot formation between rotating cone and stationary plate.
Rc＝(T1－T0)／T1(4.41)

When the torque increases with clot growth, Rc approaches to unity (Fig. 4.29 (c)).
When the shear rate is higher than 500 s － 1, Rc becomes lower than 0.5, whichcorresponds to the inhibition of the clot growth.

When the shear rate is lower than 100－1

s , on the other hand, Rc becomes higher than 0.7, which corresponds to the promotionof the clot growth (Fig. 4.29 (d)).
The stationary plane can be realized in the velocity distribution of the Couette typeflow between the clockwise disk and the counterclockwise disk. The principle isapplied to the “counter rotating rheoscope”, in which the floating object in the Couette
type flow is observed at the stationary plane (Fig. 4.30 (a, b)).

The deformation of thefloating object can be observed in the shear field of the fluid in the device [4].

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.30(a): Counter rotating rheoscope.

Fig. 4.30(b): Counter rotating rheoscope.

The deformation of erythrocytes can be observed during suspension in the shear
field at “counter rotating rheoscope”.

The deformability of erythrocytes changes fromgeneration in the bone marrow, while they circulate through the blood vessels.
Deformability of erythrocytes varies with the contents.

Deformability of eacherythrocyte can be measured, after sorting according to their density by centrifugation(see 6.1.1) [26].

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

4.2.5 Flow between parallel walls
The flow velocity distribution occurs between parallel plate walls: the flow velocity
is small near the wall, and maximum at the center.

The wall shear rate γ in the flow

velocity distribution, as shown in Fig. 4.31 (a), is able to be calculated by the similar
equations as in 4.2.2.

In Fig. 4.31 (a), the parallel velocity vectors are symmetrical

relative to the center plane between the parallel walls (with distance d).

In the flow

direction plane vertical to the wall, the envelope of the velocity vectors makes theparabola: the velocity is the maximum at the center of apex, and zero at the wall.

Fig. 4.31(a): Velocity distribution in flow between parallel walls.

Fig. 4.31(b): Force balance in flow between parallel walls.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

The origin is defined at the center between the parallel walls.
vertically towards the wall.

The y-axis is defined

A virtual thin flat plate sandwiched between y and y + dy

has the velocity of v, the width of b, and the length of l (Fig. 4.31 (b)).

When the plate

has uniform linear motion, “the friction force ΔF between the plates” balances with “theforce by the pressure difference ΔP between the upstream and the downstream”. The
friction force F is the product of the friction area of the plate (b l) and the shear stress τ.
F＝b l τ(4.42)

τ＝－η dv／dy(4.43)In Eq. (4. 43), η is the viscosity coefficient of the fluid.
maximum.

At the center (y ＝ 0), v is

The right side of the Eq. (4.43) has the negative sign, because the velocity

v decreases with y.
By substituting the Eq. (4.43) into Eq. (4.42),
F＝－b l η dv／dy(4.44)The difference of the frictional force dF between the outside of the thin flat plate
and the inside of the thin flat plate is the force resisting the flow.
dF＝－b l η(1／dy)(dv／dy)dy(4.45)

dF balances with the force by the pressure difference ΔP between the upstream andthe downstream loading at the end face of the thin flat plate (area of b dy).
dF＝b dy ΔP(4.46)

Combining Eq. (4.45) and Eq. (4.46),© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan－b l η(1／dy)(dv／dy)dy＝b dy ΔP(4.47)－(1／dy)(dv／dy)dy＝(ΔP／(l η))dy(4.48)

Both sides of Eq. (4.48) are integrated by y,－∫(1／dy)(dv／dy)dy＝(ΔP／(l η))∫dy(4.49)－dv／dy＝(ΔP／(l η))y＋C(4.50)In Eq. (4.50), C is an integration constant.
the center between the walls).

C = 0, because dv／dy = 0 at y = 0 (atTherefore, Eq. (4.50) becomes Eq. (4.51).

dv／dy＝－(ΔP／(l η))y(4.51)The left side of Eq. (4.51) is the tangent slope of the envelope of the velocity vector,
and represents the shear rate γ (=dv／dy).

γ is zero at the center, and the maximum at

the wall.
The shear rate at the wall (γw) (y ＝ d／2) is,
γw＝－(ΔP／(l η))(d／2)(4.52)The both sides of the Eq. (4.51) are integrated by y.
dv ＝－(ΔP／(l η))y dy(4.53)∫dv＝－(ΔP／(l η))∫y dy(4.54)v ＝－(ΔP／(2 l η))y 2＋D(4.55)© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, JapanEq. (4.55) represents the flow velocity distribution in the cross section. D is a
integration constant.

The flow velocity is zero at the wall (v＝0 at y＝d／2).

D＝(ΔP／(8 l η))d 2(4.56)

Eq. (4.56) is substituted into Eq. (4.55),
v＝(ΔP／(8 l η))(d 2－4y 2)(4.57)

Eq. (4.57) represents a parabolic flow velocity distribution.

By integrating theflow velocity in the channel cross-section, the flow rate is calculated.

𝑑2

𝑄 = � 2𝑏0

𝛥𝛥 2(𝑑 − 4𝑦 2 )𝑑𝑑
8𝑙𝜂
𝑑

𝑦3 2
𝑏 𝛥𝛥 2
�𝑑 𝑦 − 4 �
𝑄=
3 0
4𝑙𝜂

Q＝b d 3ΔP／(12 l η)

(4.58)(4.59)(4.60)

By modifying Eq. (4.60),
ΔP＝12 l η Q/(b d 3)

(4.61)

Eq. (4.61) is substituted into Eq. (4.52),
γw＝6 Q ／(b d 2)

(4.62)

Under the flow between parallel walls, the adhesive strength between the cells and© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

the wall of the scaffold can be estimated by the wall shear stress at the separation of thecells from the plates (Fig. 4.32) [14].

Microscopic observation of cells in the channel

between parallel walls (Figs. 4.33, 4.34) enables quantitative evaluation of the behaviorof cells: the effect of shear stress on migration, deformation, proliferation, orientationand differentiation of cells (Fig. 4.35, Fig. 4.36) [27].

Fig. 4.32: Deformation and exfoliation of cell in flow.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.33: Flow channel between parallel walls.

Fig. 4.34: Flow channel system with parallel wall for microscopic observation.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.35: Extension of cell.

Fig. 4.36: Movement, deformation, proliferation, orientation, and differentiation of cell.

4.2.6 Secondary flow
The viscosity of a fluid around a ball can be measured by the velocity of the ball

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japanfalling by the gravity (Fig. 4.37).

Fig. 4.37: Falling sphere.
η∝d 2(ρ1－ρ2 )／v(4.63)In Eq. (4.63), η is viscosity, d is diameter of the ball, v is falling velocity of the ball,and ρ1－ρ2 is the difference between the density of the ball and the fluid.
The influence of the tube wall is related to the diameter ratio d／D between the balland the circular tube.

Eq. (4.63) is established, when turbulence (see 4.3.2) does notoccur around the sphere.
Consider the flow around the moving sphere with rotation. Fluid in the vicinity ofthe surface of the rotating sphere is dragged to the rotational direction.

In the region

where the direction of dragging is same as that of moving, the density of the stream line
increases (see 4.3.2).

In the region of the higher density of the stream line, thepressure reduces (see 4.1.1, see Fig. 4.4).

By the reduction of the pressure, the ballreceives the force vertical to the moving direction (Fig. 4.38).

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.38: Magnus effect.

The cylinder rotating in the flow receives the similar effect.

The relative

movement of the rotating ball in the flow is same as the rotating ball moving in the fluid.
The effect is called “Magnus effect”.

The effect explains the curving direction of the

movement of the rotating ball.
The suspended particles rotate and shift to the center axis in the velocity distribution
in the pipe, where the velocity is high at the center and low near the wall.
called “concentration to axis”.

This isBehavior of suspended particles depends on thedensity of particle relative to that of the surrounding fluid.

By the effect, erythrocytes

flow near the central axis, and they do not flow in the vicinity of the wall (Fig. 4.39).
In this case, the destruction of erythrocytes decreases in the low shear field near thecenter.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.39: Axis concentration.

In the bending pipe, the vertical component at the main flow direction is generated
by the inertia (Fig. 4.40).

The flow component is called “the secondary flow”.

Fig. 4.40: Secondary flow in bend tube.

Consider the flow in the fluid sheared between “the side surface of the rotating innercylinder” and “the side surface of the stationary outer cylinder”, which have thecommon axis of symmetry.

The Couette type flow generates between the inner© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japancylinder and the outer cylinder.

The fluid at the vicinity of the inner cylinder isdirected to the outer cylinder by the centrifugal force.

The flow to the outside from theinside generates the continuous spiral secondary flow (Fig. 4.41). The spiral flow iscalled “Taylor vortex”.

Also in the flow between the rotating cone and the stationary

plate (see 4.2.4), the secondary flow of Taylor type vortex is generated by thecentrifugal force (Fig. 4.42).

Fig. 4.41: Secondary flow between cylinder (Taylor vortex).

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.42: Secondary flow between rotating cone and stationary plate.

The secondary flow with the centrifugal action can be decreased, on the other hand,by the outside rotation rather than the inside rotation: sandwich the fluid between thestationary inner cylinder and the rotating outside cylinder (Fig. 4.43), or between therotating concave cone and the stationary convex cone (Fig. 4.44).

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.43: Flow between rotating outer cylinder and stationary inner cylinder.

Fig. 4.44: Flow between stationary convex cone and rotating concave cone.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

The convex conical bottom surface of the stationary inner cylinder induces theuniform shear field even at the bottom (Fig. 4.43).
secondary flow.

The modification decreases the

The uniform shear rate throughout the fluid is effective for fatigue

failure testing of erythrocytes [15].

4.3 Steady flow and non-steady flow
4.3.1 Pulsatile flow
A flow, of which the velocity does not change with time, is called “steady flow”.
A flow, of which the velocity changes with time, is called “non-steady flow”.

Thevelocity of blood in the artery changes periodically correspond with the beat of the heart.
This kind of flow is called “pulsatile flow”, or “pulsating flow”.

Both pulsatile flow

and pulsating flow are included in the non-steady flow.
Do you remember the equation of motion?

The flow should be accelerated, when

the forces applied on the fluid are not balanced (see 4.2.2).

In a steady flow, the forceof the pressure difference is balanced with that of resistance (Fig. 4.21).

When theflow resistance is zero, the fluid flow continues without a pressure difference.
In the pulsatile flow, on the other hand, acceleration and deceleration arealternatively repeated.

During acceleration, the force due to the pressure difference

exceeds the force of resistance.

During deceleration, the force resistance exceeds theforce due to the pressure difference.
For example, consider the flow resistance at the aortic valve.

In the first half of thesystole (accelerated phase), the left ventricular pressure exceeds the aortic pressure.
The force due to the pressure difference is greater than that of the resistance in the aortic
valve.

The flow from the left ventricle into the aorta is accelerated (Fig. 4.45).

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.45: Pressure in pulsatile flow.

Subsequently, in the second half of the systole (deceleration phase), the left
ventricular pressure often falls below the aortic pressure (after the arrow in Fig. 4.40).
In other words, the pressure is higher at downstream than upstream.
difference decreases flow velocity.
the flow.

The pressure

The aortic valve closes, before the back motion of

The resistance of the aortic valve is very small, because of the expansion ofannulus of the aortic valve during opening (systolic phase of the heart).
When the flow path resistance at the aortic valve is large, the left ventricularpressure exceeds the aortic pressure even in the deceleration phase.
path resistance decelerates the flow.

The large flow

The greater flow resistance at the aortic valve canbe estimated by the timing, when the left ventricular pressure exceeds the aorticpressure in the deceleration period.

It is applied to the estimation of aortic valve

stenosis.
The intravascular pressure fluctuates in the pulsatile flow.

The fluctuation of theinternal pressure makes variation of the internal diameter of the vessel according to thecompliance of the vessel wall.

The local compliance governs the local wall

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japandeformation, which might make a step at the wall in the flow path (Fig. 4.46).

The

step at the connection of the artificial blood vessel with the native blood vessel causes
thrombus formation. To avoid clot formation, the compliance of the wall should be
adjusted to that of the native.

Fig. 4.46: Compliance of tube wall.

In the pulsatile flow, the periodical low shear rate inhibits the destruction of red
blood cells [15].

Since the stirring effect of the pulsatile flow washes away thestagnation region, the clot growth can be inhibited [28].

The periodical high wall

shear rate inhibits thrombus growth [29] (Fig. 4.47).

Effect of the flow on theadhesion of endothelial cells to the vascular wall is different between in the pulsatile
flow and in the steady flow [30].

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.47: Clot formation and hemolysis with shear rate.

4.3.2 Laminar flow and turbulent flow
When the inertial force is dominant than the viscous force, each layer in flow moves
independently and mixed with each other.
called “laminar flow”.

The flow with no mixture between layers is

The flow with mixture between layers, on the other hand, iscalled “turbulent flow” (Fig. 4.48).

Fig. 4.48: Tracing.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

The ratio of “inertial force” per “viscous force” is called “Reynolds number(Re)”.
Re＝ρv x／η

(4.64)In the equation (4.64), ρ is the density, v is the representative velocity, x is therepresentative length, and η is the viscosity coefficient.
Reynolds number is dimensionless (see 2.1.1), and has no units (see Q. 4.8), since it
is a ratio of the "force" to "force".

Reynolds number is an indicator for predicting

whether the flow is laminar or turbulent.

The Reynolds number, at which the flowtransitions from laminar to turbulent, is called as the “Critical Reynolds number”.
As the dimensionless number, Reynolds number is applied to the result of theexperiment at miniature model to the prediction of the flow at the actual size product(similarity law).

The simulation is planning at the same the Reynolds number.

Theturbulent flow occurs at the higher Reynolds number: the larger diameter of the bloodvessels, the higher velocity of the flow, and the lower viscosity.

With the average flow

velocity and the inside diameter at the various blood vessels, the Reynolds number inthe human blood circulation system is calculated between 0.0005 and 2000 [24].
At the steady flow in the sufficiently long straight pipe, v of the average flow
velocity at the cross-section, and the x of the inner diameter of the pipe, the critical
Reynolds number is in the range between 2000 and 4000. Since the blood flow in theartery is "pulsatile", the blood flow in the artery tends to generate turbulence.
A flow is able to be described by two ways of methods: tracing the velocity at a
fixed position, or tracing the velocity of a flowing particle.

The state of flow can bevisualized by the movement of the dye or the micro particles.
The line connecting the direction of flow is called as a streamline (Fig. 4.49).

Inthe laminar flow of the steady flow, streamline is consistent with the flow of the© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japantrajectory. Streamlines do not intersect.

When they intersect, the direction of flow

cannot be defined at the intersection.

Fig. 4.49: Streamline.

Both in the turbulent and in the unsteady flow, the streamline varies with thevariation of the flow velocity vector.

In these flows, the streamline at each moment

does not coincide with the tracings of the flow.
When the Reynolds number increases, the flow at each position in the fluid becomes
unsteady.

This is because non-steady peeling-vortex is generated by the inertia of thefluid (Fig. 4.50).

In the downstream of the obstruction, the peeling vortex is formed

alternately, which is called Karman vortex.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.50: Vortex.

A smooth wall of the flow path prevents the flow separation.

By smoothing theinner surface of the artificial ventricles, the disturbance of the blood flow can beprevented (Fig. 4.51).

The directions of the inflow and the outflow of the artificial

ventricle change the flow pattern in the artificial ventricle.

The flow pattern has aneffect on the thrombus formation (Fig. 4.52) [31].

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

Fig. 4.51: Artificial ventricle.

Fig. 4.52: Clot in artificial ventricle.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan

The mixing movement between layers increases the flow resistance in the turbulent
flow.

The stirring effect in the turbulent flow, on the other hand, increases thetransport efficiency of the material.

The agitating action of the fluid in the vicinity ofthe membrane increases the efficiency of the gas exchange in the oxygenator (theartificial lung) (see 5.3.1).

Questions
Q 4.1: Why you cannot apply the simple Kelvin-Voigt model to analysis of stresstracings, when you deform material by a stepwise strain?

Q 4.2: The “Newtonian Fluid” is flowing through a sufficiently long straight cylindricalpipe. How much decrease the resistance of the laminar flow, when the inner diameter
doubles?

Q 4.3: List cases in which the vascular resistance increases.

Q 4.4: Calculate the resistance of the laminar flow of the water through the capillary.
The capillary has the Inner radius a = 1 × 10-4 m and the length l＝0.1 m. Let thecoefficient of viscosity of the water η = 1 × 10-3 Pa s.
Q 4.5: Calculate the resistance of the pulmonary circulation, when the mean pulmonary
arterial pressure is 25 mmHg, the left atrial pressure is 5 mmHg, and the cardiac output
is 5.5 l min-1.

Q 4.6: Why the term of diastole is longer than that of systole in ventricular pumping
action?

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, JapanQ 4.7: In which case does the pulsatile-flow flow against the pressure gradient?

Q 4.8: Show that the Reynolds number is dimensionless from the definition equation.
Q 4.9: Express the torque “T” to rotate the cone at a rotational angular speed “ω”, when
the fluid of the viscosity coefficient “η” is sandwiched between the cone and the plate,using the radius “R” of the cone and the clearance angle “θ” between the cone and theplate in Figure 4.27. Calculate the torque “T” in the case of the following condition: η =0.005 Pa s, ω = 6 rad s-1, R = 0.02 m, and θ = 0.02 rad.

© Shigehiro Hashimoto 2013, Published by Corona Publishing Co., Ltd. Tokyo, Japan