PART 1 VISCOUS Flow in Pipes 113837
PART 1 VISCOUS Flow in Pipes 113837
PART 1 VISCOUS Flow in Pipes 113837
Turbomachines
Introduction
Fluid flow is classified as external or internal, depending on whether the fluid is forced to
flow over a surface or in a conduit.
The transport of a fluid (liquid or gas) in a closed conduit (commonly called a pipe if it is
of round cross section or a duct if it is not round) is extremely important in our daily
operations. Such applications range from the large, pipeline that carries crude oil e.g.
(from Chad to Cameroon), water pipes in our homes and the distribution system that
delivers the water from the city well to the house, and irrigation pipes that bring water
in farms.
The study of the motion of fluids that can be approximated as incompressible (such as
liquids, especially water, and gases at low speeds) is usually referred to as
hydrodynamics. Hydraulics is a subcategory of hydrodynamics that deals with liquid
flows in pipes and open channels. In this line, Applied Hydraulics is concerned with the use
of fluids to perform mechanical tasks. Example of such tasks are Hydraulic machines
consisting of machinery and tools that use liquid fluid power to do simple work, e.g. pump and turbines.
In this chapter we consider internal flow where the conduit is completely filled with the
fluid, and the flow is driven primarily by a pressure difference. In this context, we will
focus on the application of basic principles of fluid mechanics to the transportation of
incompressible flow of viscous fluids in pipes and ducts.
Page 1 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
fluids such as air at high velocities is typically turbulent. A flow that alternates
between being laminar and turbulent is called transitional. Osborne Reynolds (1842–
1912), a British scientist and mathematician, was the first to distinguish the difference
between these two classifications of flow.
After exhaustive experiments in the 1880s, Osborne Reynolds discovered that the
flow regime depends mainly on the ratio of inertial forces to viscous forces in the
fluid. This ratio is called the Reynolds number and is expressed for internal flow in a
circular pipe as
𝜌𝑉𝐷
𝑅𝑒 =
𝜇
where 𝑉 is the average velocity of the flow in the pipe, 𝐷 is a characteristic length of the
geometry (e.g. the diameter of a pipe), 𝜌 the fluid density, and 𝜇 the kinematic fluid
viscosity.
For flow through noncircular pipes, the Reynolds number is based on the hydraulic
diameter 𝐷ℎ defined as
4𝐴𝑐
𝐷ℎ =
𝑝
Where 𝐴𝑐 is the cross sectional area of the pipe, and 𝑝 is its wetted perimeter.
- For a cylindrical pipe of diameter 𝐷, the hydraulic diameter is
4(𝜋𝐷2 /4)
𝐷ℎ = =𝐷
𝜋𝐷
- For a rectangular duct of dimensions 𝑎 and 𝑏, the hydraulic diameter is
4(𝑎𝑏) 2𝑎𝑏
𝐷ℎ = =
2(𝑎 + 𝑏) (𝑎 + 𝑏)
which reduces to 𝑎 for a square duct.
The transition from laminar to turbulent flow does not depend only on the
geometry and flow conditions, it also depend on the degree of disturbance of the
flow by surface roughness, pipe vibrations, and fluctuations in the upstream flow.
For general engineering purposes, the following values are appropriate for flow in
pipes under most practical conditions:
- The flow is laminar if 𝑅𝑒 ≤ 2100
- The flow is turbulent if 𝑅𝑒 > 4000
- For Reynolds numbers between these two limits, the flow may switch
between laminar and turbulent conditions in a random fashion (transitional
flow)
Page 2 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
Figure: Entrance region, developing flow, and fully developed flow in a pipe system.
The region of the flow in which the effects of the viscous shearing forces caused by
fluid viscosity are felt is called the velocity boundary layer or just the boundary layer.
A virtual boundary surface divides the flow in a pipe into two regions: the boundary
layer region, in which the viscous effects and the velocity changes are significant, and
the irrotational (core) flow region, in which the frictional effects are negligible and
the velocity remains essentially constant in the radial direction. Flow in the entrance
region is called hydrodynamically developing flow since this is the region where the
velocity profile develops. The region beyond the entrance region in which the velocity
profile is fully developed and remains unchanged is called the hydrodynamically fully
developed region. The flow is said to be fully developed when the normalized
temperature profile remains unchanged as well. Hydrodynamically fully developed
flow is equivalent to fully developed flow when the fluid in the pipe is not heated or
Page 3 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
cooled since the fluid temperature in this case remains constant throughout. The
necessary condition for observing fully developed flow in a pipe is that the pipe
should be sufficiently long (relative to the entry length) so that the entry effect are
negligible.
To obtain a relationship giving the pressure gradient along the pipe in terms of the
shear stress at the walls of the pipe (wall shear stress), we write the force balance
on a volume element of radius 𝑅 and thickness 𝑧. Here, as the fluid velocity tends to
decrease as we approach the wall in any direction, the wall shear stress 𝜏𝑤 is
opposite to the fluid velocity. The equilibrium condition of the fluid element gives
𝜋𝑅2 𝑃(𝑧) − 𝜋𝑅2 𝑃(𝑧 + 𝑑𝑧) − 2𝜋𝑅𝑑𝑧𝜏(𝑅) = 0
Dividing by 𝜋𝑅𝑑𝑧 and rearranging, gives
𝑃(𝑧 + 𝑑𝑧) − 𝑃(𝑧)
𝑅 + 2𝜏(𝑅) = 0
𝑑𝑧
Taking the limits as 𝑑𝑧 → 0, we have
𝑑𝑃 𝜏(𝑅) −2𝜏𝑤
= −2 =
𝑑𝑧 𝑅 𝑅
As 𝑑𝑃/𝑑𝑧 is constant, the wall shear stress 𝜏𝑤 is also constant because the flow is
fully developed.
To obtain expression for the fluid velocity, we integrate twice the relation
𝜇 𝑑 𝑑𝑢 𝑑𝑃
(𝑟 ) =
𝑟 𝑑𝑟 𝑑𝑟 𝑑𝑥
hence
1 𝜕𝑝 2
𝑢(𝑟) = ( ) 𝑟 + 𝑐1 ln 𝑟 + 𝑐2
4𝜇 𝜕𝑧
Because the velocity must be finite at the center of the cylinder, 𝑐1 must the zero.
Applying the no-slip condition at the wall (𝑟 = 𝑅), 𝑢 = 0. Hence
𝟏 𝝏𝒑
𝑢(𝑟) = ( ) (𝒓𝟐 − 𝑹𝟐 )
𝟒𝝁 𝝏𝒛
Thus, the cross section velocity is parabolic. To obtain a relationship between the
volume flow rate 𝑄 passing through the tube and the pressure gradient, we consider a
flow through a differential, washer-shape ring. Since 𝑢 is constant on this ring, the
volume rate of flow through the differential area 𝑑𝐴 = (2𝜋𝑟)𝑑𝑟 is
𝑑𝑄 = 𝑢𝑑𝐴 = 𝑢(2𝜋𝑟)𝑑𝑟
and therefore,
𝑅
𝜋𝑅4 𝜕𝑝
𝑄 = ∫ 𝑢(2𝜋𝑟)𝑑𝑟 = − ( )
8𝜇 𝜕𝑧
0
Page 5 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
More practically, this relation can be expressed in terms of the pressure drop, Δ𝑝, which
occurs over a length 𝐿 along the tube. Since the flow is directed toward negative
pressure gradients, we set
𝜕𝑝 ∆𝑝
=− ;
𝜕𝑧 𝐿
hence
𝝅𝑫𝟒 ∆𝒑
𝑸=
𝟏𝟐𝟖𝝁𝑳
The above results confirm the following properties of laminar pipe flow. For a horizontal
pipe the flowrate is (a) directly proportional to the pressure drop, (b) inversely
proportional to the viscosity, (c) inversely proportional to the pipe length, and (d)
proportional to the pipe diameter to the fourth power.
This flow, the properties of which were first established experimentally by two independent
workers, G. Hagen 11797–18842 in 1839 and J. Poiseuille (1799–1869) in 1840, is termed
Hagen–Poiseuille flow. The above equation is commonly referred to as Poiseuille’s law. This
poiseille’s law specifies that for a given flow rate, the pressure drops and thus the required
pumping power is proportional to the length of the pipe and the viscosity of the fluid, but it
is inversely proportional to the fourth power of the radius(or diameter) of the pipe .
Therefore, the pumping power requirement for a laminar-flow in piping system can be
reduced by a factor of 16 by doubling the pipe diameter. Of course the benefits of the
reduction in the energy costs must be weighed against the increased cost of construction
due to using a larger-diameter pipe.
In terms of the average velocity vavg = 𝑄/𝜋𝑅2 , the Poiseuille relation becomes
𝑅2 ∆𝑝
vavg =
8𝜇𝐿
The maximum velocity occurs at the center of the tube and is given by
𝑅2 𝜕𝑝 𝑅 2 ∆𝑝
v𝑚𝑎𝑥 = − ( ) = = 2vavg
4𝜇 𝜕𝑧 4𝜇𝐿
Note that the velocity distribution can be written in terms of the maximum velocity as
𝑟2 𝑟2
𝑢(𝑟) = v𝑚𝑎𝑥 (1 − 2 ) = 2vavg (1 − 2 )
𝑅 𝑅
Page 6 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
Page 7 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
8𝜇𝐿vavg 32𝜇𝐿vavg
∆𝑃 = =
𝑅2 𝐷2
∆𝑃 is the pressure drop along the flow direction due to viscous effects. It represents an
irreversible pressure loss, and it is sometimes called pressure loss. In practice, it is
convenient to express the pressure loss for all types of fully developed internal flows
(laminar or turbulent flows, circular or noncircular pipes, smooth or rough surfaces,
horizontal or inclined pipes) as
2
𝐿 𝜌Vavg
∆𝑃 = 𝑓
𝐷 2
2
Where 𝜌Vavg /2 is the dynamic pressure and 𝑓 is the Darcy-Weisbach friction factor,
named after the French engineer Henry Darcy (1803–1858) and the German Julius
Weisbach (1806–1871), the two engineers who provided the greatest contribution to
its development. The frictions factor is a coefficient that takes into account the effects
of the geometry of the pipe and its roughness.
From the above relation, the friction factor can be expressed as
𝐷
∆𝑃 ( )
𝑓= 𝐿
2
𝜌Vavg
2
And since the pressure drop is related to the wall shear stress by the relation
2𝜏𝑤 𝐿 𝐿
∆𝑃 = = 4 𝜏𝑤 ,
𝑅 𝐷
the Darcy friction factor in terms of the wall shear stress becomes
𝐷
𝐿 ( ) 8𝜏 𝑤
𝑓 = 4 𝜏𝑤 𝐿2 =
𝐷 𝜌Vavg 𝜌V2avg
2
For laminar flows, the Darcy friction factor can be written terms of the Reynolds number
by equating
2
32𝜇𝐿vavg 𝐿 𝜌Vavg
∆𝑃 = =𝑓
𝐷2 𝐷 2
Solving for 𝑓 gives
64𝜇 64
𝑓= =
ρDvavg 𝑅𝑒
This equation shows that in laminar flow, the friction factor is a function of the Reynolds number
only and is independent of the roughness of the pipe surface (assuming, of course, that the roughness
is not extreme).
Page 8 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
In the analysis of piping systems, pressure losses are commonly expressed in terms of the
equivalent fluid column height, called the head loss ℎ𝐿 . The pressure head loss in the pipe
is the pipe head loss and is defined as
2
∆𝑃 𝐿 Vavg
ℎ𝐿 = =𝑓
𝛾 𝐷 2𝑔
The head loss ℎ𝐿 represents the additional height that the fluid needs to be raised by a pump in
order to overcome the frictional losses in the pipe. The head loss is caused by viscosity, and
it is directly related to the wall shear stress. Once the pressure loss (or head loss) is
known, the required pumping power to overcome the pressure loss is determined from
Ẇpump,loss = 𝑄∆𝑃 = 𝑄𝛾ℎ𝐿 = 𝑚̇𝑔ℎ𝐿 .
Page 9 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
ℎ𝐿 𝑚𝑎𝑗𝑜𝑟 , and the head loss in the various pipe components, termed the minor loss and
denoted ℎ𝐿 𝑚𝑖𝑛𝑜𝑟 . That is
ℎ𝐿 = ℎ𝐿 𝑚𝑎𝑗𝑜𝑟 + ℎ𝐿 𝑚𝑖𝑛𝑜𝑟
The head loss designations of “major” and “minor” do not necessarily reflect the relative
importance of each type of loss. For a pipe system that contains many components and a
relatively short length of pipe, the minor loss may actually be larger than the major loss.
1.4.1 Major Losses and the Effects of rough walls
It was not known until experiments in 1800 by Coulomb that surface roughness has an effect
on friction resistance. It turns out that the effect is negligible for laminar
pipe flow, and all the laminar formulas derived in this section are valid for rough walls also.
However, the friction factor in fully developed turbulent pipe flow depends on the
Reynolds number and the relative roughness 𝜀/𝐷, which is the ratio of the mean height of
roughness of the pipe to the pipe diameter.
A fundamental difference between laminar and turbulent flow is that the shear stress for
turbulent flow is a function of the density of the fluid, 𝜌. For laminar flow, the shear stress
is independent of the density, leaving the viscosity, 𝜇 as the only important fluid property.
Although the pressure drop for laminar pipe flow is found to be independent of the
roughness of the pipe, it is necessary to include this parameter when considering turbulent
flow. For turbulent flow there is a relatively thin viscous sublayer formed in the fluid near
the pipe wall and pressure drop is a function of the wall roughness. If we introduce the
relative roughness of the pipe, 𝜀/𝐷 ≪ 1, the pressure drop along a horizontal pipe can be
written as
ℓ 𝜌𝑉 2
∆𝑝 = 𝑓
𝐷 2
where 𝑓 is the friction factor, with
𝜀
𝑓 = 𝜙 ( , 𝑅𝑒 )
𝐷
1.4.2 The Moody Chart and the Colebrook Equation
For laminar fully developed flow, the value of the friction factor 𝑓 is simply 𝑓 = 64/𝑅𝑒 ,
independent of 𝜀/𝐷. For turbulent flow, the functional dependence of the friction factor on
the Reynolds number and the relative roughness, 𝑓 = 𝜙(𝜀/𝐷, 𝑅𝑒 ) is a rather complex one
that cannot, as yet, be obtained from a theoretical analysis. Available informations are so far
are obtained from an exhaustive set of experiments and usually presented in tabular,
graphical, and functional forms obtained by curve-fitting experimental data. Most such
experiments were conducted by J. Nikuradse in 1933, followed by the works of others where
friction factor was calculated from measurements of the flow rate and the pressure drop.
Page 10 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
In 1939, Cyril F. Colebrook (1910–1997) combined the available data for transition and
turbulent flow in smooth as well as rough pipes into the
following implicit relation known as the Colebrook
equation:
1 𝜀/𝐷 2.51
= −2.0 log ( + )
√𝑓 3.7 𝑅𝑒 √𝑓
This relation is valid for all turbulent flows and is
consistent with two previously derived relations dealing
with the friction factor in some limiting cases:
(i) The Prandtl equation which was derived assuming
a perfectly smooth pipe (𝜀/𝐷 → 0):
1
= 2.0 log(𝑅𝑒 √𝑓) − 0.8
√𝑓
(ii) The von Kármán equation that gives the friction
factor for fully rough turbulent flows or simply fully
turbulent flows (𝑅𝑒 → ∞):
1 𝜀/𝐷
= −2.0 log ( )
√𝑓 3.7
In 1942, the American engineer Hunter Rouse (1906–1996) verified Colebrook’s equation
and produced a graphical plot of 𝑓 as a function of 𝑅𝑒 and the product 𝑅𝑒 √𝑓. He also
presented the laminar flow relation and a table of commercial pipe roughness. Two years
later, Lewis F. Moody (1880–1953) redrew Rouse’s diagram into the form commonly used
today as the Moody chart.
The Moody chart presents the Darcy friction factor for pipe flow as a function of Reynolds
number and 𝜀/𝐷 over a wide range. Although this chart was developed for circular pipes, it
can also be used for noncircular pipes by replacing the diameter with the hydraulic
diameter.
It should also be noted that with advances in technology, newly manufactured pipes have
very small relative roughness as shown. However, after considerable use, most pipes
(because of a buildup of corrosion or scale) may have a relative roughness that is
considerably larger (perhaps by an order of magnitude) than that given.
A difficulty with its use is that it is implicit in the dependence of 𝑓; hence the determination
of the friction factor requires iteration. For the iteration process to be successful, an
acceptable approximation of the friction factor is needed. An approximate explicit relation
for 𝑓 was given by S. E. Haaland in 1983 as
1 𝜀/𝐷 1.11 6.9
= −1.8 log [( ) + ]
√𝑓 3.7 𝑅𝑒
Page 11 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
This explicit relation gives results within 2% of those obtained from the Colebrook relation;
hence it can serve as a good first guest for the iteration procedure when more accuracy is
needed.
Types of Fluid Flow Problems
In the design and analysis of piping systems that involve the use of the Moody chart
(or the Colebrook equation), we usually encounter three types of problems (the
fluid and the roughness of the pipe are assumed to be specified in all cases):
1. Determining the pressure drop (or head loss) when the pipe length and diameter
are given for a specified flow rate (or velocity)
2. Determining the flow rate when the pipe length and diameter are given for a
specified pressure drop (or head loss)
3. Determining the pipe diameter when the pipe length and flow rate are given for a
specified pressure drop (or head loss)
Problems of the second type and third type are commonly encountered in
engineering design (in the selection of pipe diameter, for example, that minimizes the
sum of the construction and pumping costs), but the use of the Moody chart with such
problems requires an iterative approach.
Problems of the first type are straightforward and can be solved directly by using the
Moody chart.
In problems of the second type, the diameter is given but the flow rate is unknown. A
good guess for the friction factor in that case is obtained from the completely
turbulent flow region for the given roughness. This is true for large Reynolds
numbers, which is often the case in practice. Once the flow rate is obtained, the
friction factor is corrected using the Moody chart or the Colebrook equation, and the
process is repeated until the solution converges. (Typically only a few iterations are
required for convergence to three or four digits of precision.)
In problems of the third type, the diameter is not known and thus the Reynolds
number and the relative roughness cannot be calculated. Therefore, we start
calculations by assuming a pipe diameter. The pressure drop calculated for the
assumed diameter is then compared to the specified pressure drop, and calculations
are repeated with another pipe diameter in an iterative fashion until convergence.
To avoid tedious iterations in head loss, flow rate, and diameter calculations,
Swamee and Jain (1976) proposed the following explicit relations that are accurate
to within 2 percent of the Moody chart:
−2
𝑄2 𝐿 𝜀/𝐷 𝜈𝐷 0.9 10−6 < 𝜀/𝐷 < 10−2
ℎ𝐿 = 1.07 {ln [ + 4.62 ( ) ]} ,
𝑔𝐷 5 3.7 𝑄 3 × 103 < 𝑅𝑒 < 3 × 108
Page 12 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
0.5 𝜀 0.5
𝑔𝐷 5 ℎ𝐿 𝐷 3.17𝜈 2 𝐿
𝑄 = −0.965 ( ) ln [ +( ) ], 𝑅𝑒 > 2000
𝐿 3.7 𝑔𝐷 3 ℎ𝐿
4.75 0.04
1.25
𝐿𝑄2 9.4
𝐿 5.2 10−6 < 𝜀/𝐷 < 10−2
𝐷 = 1.07 [𝜀 ( ) + 𝜈𝑄 ( ) ] ,
𝑔ℎ𝐿 𝑔ℎ𝐿 5 × 103 < 𝑅𝑒 < 3 × 108
𝐾𝐿 = 𝜙(geometry, 𝑅𝑒 )
That is Although 𝐾𝐿 is dimensionless, it often is not correlated in the literature with the
Reynolds number and roughness ratio but rather simply with the raw size of the pipe in, say,
inches. Almost all data are reported for turbulent flow conditions. A single pipe system may
have many minor losses. In this case, the net minor loss head is obtained as the sum of
individual losses
𝑉2
ℎ𝐿𝑚𝑖𝑛𝑜𝑟 = ∑ 𝐾𝐿
2𝑔
Table of some loss coefficients is given in the annex.
Once all the loss coefficients are available, the total head loss in a piping system is
determined from
V2avg 2
𝐿 Vavg
ℎ𝐿,𝑡𝑜𝑡𝑎𝑙 = ℎ𝐿, 𝑀𝑎𝑗𝑜𝑟 + ℎ𝐿,𝑚𝑖𝑛𝑜𝑟 = ∑ 𝐾𝐿 +∑𝑓
2𝑔 𝐷 2𝑔
Or in terms of volume discharges as
Q2 𝐿 Q
2
ℎ𝐿,𝑡𝑜𝑡𝑎𝑙 = ∑ 𝐾𝐿 2 + ∑ 𝑓
2A 𝑔 𝐷 2A2 𝑔
Piping systems often involve sudden or gradual expansion or contraction sections to
accommodate changes in flow rates or properties such as density and velocity. The
losses are usually much greater in the case of sudden expansion and contraction (or
wide-angle expansion) because of flow separation. By combining the equations of
mass, momentum, and energy balance, the loss coefficient for the case of a sudden
expansion is approximated as
2
Asmall
𝐾𝐿 = 𝛼 (1 − )
Alarge
where 𝛼 is the kinetic energy correction factor at the exit of the pipe, Asmall and Alarge are
the cross-sectional areas of the small and large pipes, respectively. Note that when a pipe
discharges into a reservoir, 𝐾𝐿 = 𝛼 as Alarge ≫ Asmall . No such relation exists for a sudden
contraction, and the 𝐾𝐿 values in that case must be read from a chart or table.
The losses due to expansions and contractions can be reduced significantly by
installing conical gradual area changers (nozzles and diffusers) between the small
and large pipes.
Piping systems also involve changes in direction without a change in diameter, and
such flow sections are called bends or elbows. The losses in these devices are due to
flow separation (just like a car being thrown off the road when it enters a turn too
fast) on the inner side and the swirling secondary flows that result. The losses during
changes of direction can be minimized by making the turn “easy” on the fluid by using
circular arcs (like 90° elbows) instead of sharp turns.
Page 14 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
Page 15 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
the layout of the system is prepared, the determination of the pipe diameters and the
pressures throughout the system, while remaining within the budget constraints, typically
requires solving the system repeatedly until the optimal solution is reached. Computer
modeling and analysis of such systems make this tedious task a simple chore.
Also, is that the total head loss through the system equals the sum of the head loss in each
pipe, including the minor losses:
ℎ𝐴→𝐵 = ℎ𝐿,1 + ℎ𝐿,2 + ℎ𝐿,3
In terms of the friction and minor losses in each pipe, we could rewrite this as
𝑉𝑗2
ℎ𝐴→𝐵,𝑚𝑖𝑛𝑜𝑟 = ℎ𝐿,1 + ℎ𝐿,2 + ℎ𝐿,3 = ∑ ℎ𝐿,𝑗 = ∑ 𝐾𝐿,𝑗
2𝑔
𝑗 𝑗
And the net major loss is
ℓ𝑗 𝑉𝑗2
ℎ𝐴→𝐵,𝑀𝑎𝑗𝑜𝑟 = ∑ 𝑓𝑗
𝐷𝑗 2𝑔
𝑗
If the flowrate is given, it is a straightforward calculation to determine the head loss or
pressure drop (Type I problem). If the pressure drop is given and the flowrate is to be
calculated (Type II problem), an iteration scheme is needed. In this situation none of the
friction factors, 𝑓𝑖 , are known, so the calculations may involve more trial-and-error attempts
than for corresponding single pipe systems. The same is true for problems in which the pipe
diameter (or diameters) is to be determined (Type III problems).
Page 16 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
Example:
Given a three-pipe system connected in series where the total pressure drop across the
system is Δ𝑃 = 150 𝑘𝑃𝑎 and the elevation drop is Δ𝑧 = 5.0 𝑚. For the pipes data listed in
the table below, and in case of a flow of water with density 𝜌 = 1000 𝑘𝑔/m3 and
kinematic viscosity 𝜈 = 1.02 × 10−6 m2 /s, determine (a) the total head loss across
the system (b) the relationships between the velocities 𝑉2 , 𝑉3 and 𝑉1 (c) the relation
between the friction factor 𝑓1 , 𝑓2 , 𝑓3 and the velocity 𝑉1 . (d) use an iterative technique
to find the values of the friction factors 𝑓1 ,
𝑓2 , 𝑓3 and the fluid velocities 𝑉1 , 𝑉2 , 𝑉3 . (e)
determine the Reynolds numbers and the
volume flow rate. Minor losses are
ssumed negligible.
1.5.2 Pipes in Parallel
Page 17 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
𝑄2 𝜋 2 𝑔𝐷𝑖5
ℎ𝐿 = 2, where 𝐶𝑖 =
(∑ √𝐶𝑖 /𝑓𝑖 ) 8ℓ𝑖
Since the 𝑓𝑖 vary with Reynolds number and roughness ratio, one begins by guessing values
of 𝑓𝑖 (fully rough values are recommended) and calculating a first estimate of ℎ𝐿 . Then each
pipe yields a flow-rate estimate 𝑄𝑖 = √𝐶𝑖 ℎ𝐿 /𝑓𝑖 and hence a new Reynolds number and a
better estimate of 𝑓𝑖 . Then the process is repeated until convergence.
The analysis of piping networks, no matter how complex they are, is based on two
simple principles:
1. Conservation of mass throughout the system must be satisfied. This is done by
requiring the total flow into a junction to be equal to the total flow out of the junction
for all junctions in the system. Also, the flow rate must remain constant in pipes
connected in series regardless of the changes in diameters.
2. Pressure drop (and thus head loss) between two junctions must be the same for all
paths between the two junctions. This is because pressure is a point function and it
cannot have two values at a specified point. In practice this rule is used by requiring
that the algebraic sum of head losses in a loop (for all loops) be equal to zero. (A head
loss is taken to be positive for flow in the clockwise direction and negative for flow in
the counterclockwise direction.)
Example:
Given a three-pipe system connected in parallel
where the total pressure drop across the system
is Δ𝑃 = 150 𝑘𝑃𝑎 and the elevation drop is Δ𝑧 =
5.0 𝑚. For the pipes data listed in the table
below, and in case of a flow of water with density
𝜌 = 1000 𝑘𝑔/m3 and kinematic viscosity
𝜈 = 1.02 × 10−6 m2 /s, determine (a) the total head loss across the system (b) the
relationships between the velocities 𝑉2 , 𝑉3 and 𝑉1 (c) the relation between the friction
factor 𝑓1 , 𝑓2 , 𝑓3 and the velocity 𝑉1 . (d) use an iterative technique to find the values of
the friction factors 𝑓1 , 𝑓2 , 𝑓3 and the fluid velocities 𝑉1 , 𝑉2 , 𝑉3 . (e) determine the
Reynolds numbers and the volume flow rate. Minor losses are assumed negligible.
Page 18 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
(mass flow and volume flow), and global properties (visualization of the entire flow
field).
The flow rate measurement techniques range from very crude to very elegant. The flow
rate of water through a garden hose, for example, can be measured simply by collecting
the water in a bucket of known volume and dividing the amount collected by the collection
time. A crude way of estimating the flow velocity of a river is to drop a float on the river
and measure the drift time between two specified locations. At the other extreme, some
flowmeters use the propagation of sound in flowing fluids while others use the
electromotive force generated when a fluid passes through a magnetic field
We shall concentrate in this section on velocity and volume flow measurements.
Page 19 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
flow to avoid significant errors that may be caused by misalignment. Because of the
slow response of the fluid-filled tubes leading to the pressure sensors, it is not useful
for unsteady flow measurements.
1.6.2 Obstruction Flowmeters: Orifice, Venturi, and Nozzle Meters
Consider incompressible steady flow of a fluid in a horizontal pipe of diameter 𝐷 that
is constricted to a flow area of diameter 𝑑, as shown. The mass balance equation
between a location before the constriction constriction (point 1) and the location
where constriction occurs (point 2) are written as
𝐴1 𝑉1 = 𝐴2 𝑉2 →
𝐴2
𝑉1 = ( ) 𝑉2
𝐴1
Or
𝑉1 = (𝑑2 /𝐷2 )𝑉2
and the Bernoulli equation gives
𝑃1 𝑉12 𝑃2 𝑉22
+ = +
𝜌 2 𝜌 2
Page 20 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
Of the numerous types of obstruction meters available, those most widely used are
orifice meters, flow nozzles, and Venturi meters. For standardized geometries, the
experimentally determined data for discharge coefficients are expressed as:
Page 21 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
2.1 8
91.71𝛽2.5
Orifice meter: 𝐶𝑑 = 0.5959 + 0.0312𝛽 − 0.184𝛽 + 0.75
𝑅𝑒𝐷
6.53𝛽0.5
Nozzle meter: 𝐶𝑑 = 0.9965 − + 0.5
𝑅𝑒𝐷
Where 𝑅𝑒𝐷 is the approach Reynolds number. These relations are valid for 0.25 ,
0.25 < 𝛽 < 0.75 , and 104 < 𝑅𝑒𝐷 < 107 . Precise values of 𝐶𝑑 depend on the
particular design of the obstruction, and
thus the manufacturer’s data should be
consulted when available. Also, the
Reynolds number depends on the flow
velocity, which is not known a priori.
Therefore, the solution is iterative in nature
when curve-fit correlations are used for 𝐶𝑑 .
The third and final type of obstruction meter is the venturi, named in honor of
Giovanni Venturi, an Italian physicist who first tested conical expansions and
contractions. The discharge coefficient for a standard venturi meter is given by
Venturi meter: 𝐶𝑑 = 0.9858 − 0.196𝛽4.5
The choice of meter depends on the loss and the cost. Other flow rate measurement
devices include Ultrasonic Flowmeters, Electromagnetic flowmeters, vortex
flowmeters and Laser Doppler Velocimetry.
Page 22 of 23
MECE 5204 Viscous flow in pipes
Turbomachines
Page 23 of 23