Module 2, 3,4 and 5
Module 2, 3,4 and 5
Module 2, 3,4 and 5
PA - ρ Ahg – P0A = 0
P = P0 + ρ gh
Hydraulics
A = area
v = velocity
A1v1 = A2v2
“What comes in comes out.”
Av= V/s (volume flow rate) = constant
Bernoulli’s Equation
Consider the change in total energy of the fluid as it moves from the inlet to the
outlet.
Δ Etotal = Wdone on fluid - Wdone by fluid
Δ Etotal = (1/2mv22 + mgh1) – (1/2mv12 + mgh2)
Wdone on fluid - Wdone by fluid = (1/2mv22 + mgh1) – (1/2mv12 + mgh2)
P2V2 - P1V1 = (1/2mv22 + mgh1) – (1/2mv12 + mgh2)
P2 – P1 = (1/2ρ v12 + ρ gh1) – (1/2ρ v12 + ρ gh1)
A2 < A1 ; V2 > V1
According to Bernoulli’s Law, pressure at A2 is lower.
Choked flow: Because pressure cannot be negative,
total flow rate will be limited. This is useful in
controlling fluid velocity.
P2 + 1/2ρ v12 = P1 + 1/2ρ v12 ; ΔP = ρ/2*(v22 – v12)
Atomizer
(Demonstration)
∴ V2 = √(2g(h1-h2))
Q = Av2 = A √(2g(h1-h2))
Pitot – Static Tube
Discharge
Substituting value of v2 in above equation
.
co-efficient of discharge is high
Can be installed vertically, Horizontally, Inclined.
They are more accurate and can be used for a wide range of flows.
Around 90% of pressure drop can be recovered.
Assumptions
Assumptions made for deriving the expression of discharge through
the orifice meter is as mentioned here.
1. Fluid is ideal, i.e. in viscid and incompressible.
2. Fluid flow is steady and continuous
3. Fluid flow is irrotational
4. Frictionless inner surface
Expression for the rate of flow through Orifice
,
applying Bernoulli’s equation at section 1 and section 2.
Let A0 is the area of the orifice
Co-efficient of contraction, CC = A2/A0
Let us recall the continuity equation and we will have following equation
Thus we will use the value of CC in above equation of discharge Q and we will
have following result for rate of flow or discharge through orifice meter.
VENTURI METER
2
BASIC TERMS
• Velocity of fluid which passes through a
Flow Rate(v) given area per second (m/s).
•Converging part
•Throat
•Diverging part
5
6
WORKING PRINCIPLE
Venturi meter works under the principle of Bernoulli's
equation and Continuity equation.
7
Where, P Pressure
V Velocity
ρ1 Density of converging fluid
ρ2 Density of throat fluid
A1 Pipe area
A2 Throat area
V1 Velocity of converging fluid
V2 Velocity of throat
8
TYPES OF VENTURIMETER
•Horizontal venturimeter
9
Incompressible fluids
•B(m2,ρ2,p2,v2)
•A(m1,ρ1,p1,v1)
Z2
Z1
Datum
At pt:-A At pt:-B
½[(A2/A1)*V22-V22]=∆P/ρ
V22=(A12/A22-A12) * 2∆P/ρ
13
Q α √2∆P
Theoretical-- Q2=M.√(2∆Pg)/ω
Q2=Cd.E.M.A2√2g {hm(sg-1)-(Zx-Zy)}
15
CONSTRUCTION
The entry of the venture is cylindrical in shape to match the size of
the pipe through which fluid flows. This enables the venture to be
fitted to the pipe.
After the entry, there is a converging conical section with an
included angle of 19’ to 23’.
Following the converging section, there is a cylindrical section
with minimum area called as the throat.
After the throat, there is a diverging conical section with an
included angle of 5’ to 15’.
Openings are provided at the entry and throat of the venturi meter
for attaching a differential pressure sensor.
17
MANOMETER
The differential pressure sensor used here is
Manometer.
Manometer is a device to measure pressure.
A common simple manometer consists of a U shaped
tube of glass filled with some liquid.
Manometers measure a pressure difference by
balancing the weight of a fluid column between the two
pressures of interest. Large pressure differences are
measured with heavy fluids, such as mercury (high
density).
Small pressure differences, such as those experienced in
experimental wind tunnels or venturi flowmeters are
measured by lighter fluids such as water .
18
19
OPERATION
• P1 • P ses
• p2
• High • P1-p2
• Less
angle angle
22
SAMPLE PROBLEM
M=A1/√A12-A22
A1=π*d2/4=π*152/4=176.71 , A2=π*d2/4=π*7.52/4=44.178
M= 1.03
Q2=0.97
* 1* 1.03* 44.178√2*9.8*17.5(13.6-1)
=0.02901 m3/sec
23
PRESSURE IN PIPELINE
24
APPLICATIONS
ADVANTAGES
DISADVANTAGES
• Highly expensive
1
WORKING PRINCIPLE
●
CONSTRUCTION
●
PRINCIPLE
●
NUMERICALS
●
●It is a device used to measure the speed and flow rate or discharge
of fluid through a pipe.” Venturimeter is work on Bernoulli’s
equation.
● PRINCIPLE:
The principle of venturimeter is that when a fluid flows through
the venturimeter,
it accelerates in the convergent section and decelerates in the diverg
resulting in a drop in the static pressure followed by a pressure reco
in the flow direction.
USE OF VENTURIMETER
of Bernoulli’s equation.
●Bernoulli’s Statement: It states that in a steady, ideal
flow of an incompressible fluid, the total energy at
any point of the fluid is constant. The total energy
consists of pressure energy, kinetic energy and
potential energy or datum energy
Mathematically
PRESSURE ENERGY+KINETIC ENERGY+POTIENTAL
●
ENERGY=CONSTANT
●Here all the energies are taken per unit weight of the
fluid.
●The Bernoulli’s equation for the fluid passing through
the section 1 and 2 are given by
●Bernoulli’s equation for the fluid passing through the
section 1 and 2 are given by
p/(rho)*g+v2/2*g+z=contant
●
p1/(rho)*g+v12/2*g+z1=p2/(rho)*g+v22/2*g+z2
●
CONSTRUCTION
●A venturimeter is a device used for measuring the rate of flow of a
fluid flowing through a pipe.
● The main parts of a venturimeter are:
● Main parts of Venturimeter:-
● 1. Converging part
● 2. Throat
● 3. Diverging Part
CONVERGENT PART
● p1/(rho)*g+v12/2*g+z1=p2/(rho)*g+v22/2*g+z2
● As the pipe is horizontal, so z1 = z2
● p1/(rho)*g+v12/2*g=p2/(rho)*g+v22/2*g
● p1-p2/(rho)*g=v22-v12/2*g ...............................1
(p1 – p2)/ρg is the difference of pressure heads at section 1 and 2
and it is equal to h. So
h=(p1-p2)/(rho)*g ..................................................2
● Substituting this value of h (2) in equation (1), we get
● h=v22-v12/2*g ...........................3
● Applying Continuity Equation:
● a1*v1=a2*v2
● v1=a2*v2/a1 ......................................4
● Substituting this value of v1 in equation (2) and solving, we get
● v2=a1/sqrt(a12-a22)*sqrt(2*g*h)
● Since Q=a2*v2
● Q={(a1*a2)√(a12-a22)}*√(2*g*h)
Q is the theoretical discharge under ideal conditions. Actual
●
h=x[Sh/So-1]
●
●Case II :If differential Manometer contains liquid lighter than liquid
flowing through pipe:
● h=x[1-Sl/So]
● Where
● Sl=Specific Gravity of lighter liquid
● So=Specific gravity of liquid flowing through pipe
APPLICATION: Carburettor
●They are large in size and, therefore, where space is limited, they
can not be used.
● Initial costs, installation and expensive maintenance.
● Requires a long placement length. ...
● It can not be used in pipes of less than 7.5 cm in diameter.
● Maintenance is not easy
Problems:
Laminar Flow: In this type of flow, fluid particles moves along smooth straight parallel
paths in layers or laminas, with one layer gliding smoothly over an adjacent layer, the paths
transverse direction to the main flow. This irregular, fluctuating motion can be regarded as
Transitional
Turbulent
Types of flow depend on the Reynold number , ρVd
Re = --------
µ
Re < 2000 – flow is laminar
= 2
t
2 dt
=
1
t 0
Intensity of turbulence :
C lV 2
f
- So, Intensity of Turbulence = h f
2 gD
From the experimental measurement on turbulent flow through pipes, it has observed
h 1 2
OR h
f
g f
2 gD
Moody Diagram :
Developed to provide the friction factor for turbulent flow for various values of Relative roughness and
Reynold’s number!
From experimentation, in turbulent flow, the friction factor (or head loss) depends upon velocity of fluid
V, dia. of pipe D, density of fluid ρ, viscosity of fluid µ, wall roughness height ε.
So, f = f1 (V,D, ρ, µ, ε)
By the dimensional analysis,
f f1 VD
, D 0 , Where
D called relative roughness.
Key points about the Moody Diagram –
2. 2. f = 64/Nr.
4. Beyond 4000, for a given Nr, as the relative roughness term D/ε increases (less rough), friction
factor decreases
5. For given relative roughness, friction factor decreases with increasing Reynolds number till the
7. As relative roughness increases (less rough) – the boundary of the zone of complete turbulence
shifts (increases)
Co-efficient of friction in terms of shear
stress :
2 0 dF dAvu
(P1 – P2) = f OR t uv
V
2
dA dA
( co-efficient of frictionin terms of shear stress)
Shear stress in turbulent flow
In turbulent flow, fluid particles moves randomly, therefore it is impossible to trace the
Paths of the moving particles and represents it mathematically
u b
= mean velocity of particles moving along layer A v t
u = mean velocity of particles moving along layer B du du
b
dy dy
du du
u = u b - u a
dy dy
Shear stress in turbulent flow du du
dy dy
It is the shear stress exerted by layer A on b and t known as
k 2 y 2
Reynold’s stress.
Prandtl’s mixing length theory :
Prandtl’s assumed that distance between two layers in the transverse direction
(called mixing length l) such that the lumps of fluid particles from one layer could reach the other
Layer and the particles are mixed with the other layer in such a way that the momentum of the
du 2 2 du
k y
dy dy
du du k 2 y 2
where , (Viscosity)
dy dy
n = 0 for laminar flow.
du du
For highly turbulent flow,
dy dy .
Hydrodynamically Smooth and
Rough Pipe Boundaries
2
dv
k y
dy
2 2
0
y
v V Log *
C
K e
R
y
v vm ax 2.5V * Log
e
R
v v R
m ax
2.5V * Log
V *
e
y
Velocity Distribution for turbulent
flow
Velocity Distribution v y
8.5 2.5 log
e
in a hydrodynamically V*
Rough Pipes
Velocity Distribution for turbulent
flow in terms of average Velocity (V)
Velocity Distribution V V * R
1.75 2.5 log
in a hydrodynamically V* e
smooth pipe
4 flV 2
h f
2 gD
V 2 l
P Re, ,
2 D D
V
2
l
P Re,
2 D D
f
Re,
D
P D
- But f friction factor
lV 2
2
f Re,
D
0.079
f 5
,For 4000<Re<10 laminar sublayer in smooth pipe.
Re
1
4
1 R
2 log 1.74
4f 10
Introduction
External flows past objects encompass an extremely wide
variety of fluid mechanics phenomena. Clearly the character
of the flow field is a function of the shape of the body.
For a given shaped object, the characteristics of the flow
depend very strongly on various parameters such as size,
orientation, speed, and fluid properties.
According to dimensional analysis arguments, the character
of the flow should depend on the various dimensionless
parameters involved.
For typical external flows the most important of these
parameters are the Reynolds number, Re =UL/ν , where L– is
characteristic dimension of the body.
2
Introduction
For many high-Reynolds-number flows the flow field may be
divided into two region
i. A viscous boundary layer adjacent to the surface
ii. The essentially inviscid flow outside the boundary layer
W know that fluids adhere the solid walls and they take the
solid wall velocity. When the wall does not move also the
velocity of fluid on the wall is zero.
In region near the wall the velocity of fluid particles
increases from a value of zero at the wall to the value that
corresponds to the external ”frictionless” flow outside the
boundary layer
3
Introduction
Figure 6.1: Visualization of the flow around the car. It is visible the thin
layer along the body cause by viscosity of the fluid. The flow outside
the narrow region near the solid boundary can be considered as ideal
(inviscid).
4
Introduction
The concept of boundary layer was first introduced by a
German engineer, Prandtl in 1904.
According to Prandtl theory, when a real fluid flows past a
stationary solid boundary at large values of the Reynolds
number, the flow will be divided into two regions.
i. A thin layer adjoining the solid boundary, called the
boundary layer, where the viscous effects and rotation
cannot be neglected.
ii. An outer region away from the surface of the object
where the viscous effects are very small and can be
neglected. The flow behavior is similar to the upstream
flow. In this case a potential flow can be assumed.
5
Introduction
Since the fluid at the boundaries has zero velocity, there is a
steep velocity gradient from the boundary into the flow. This
velocity gradient in a real fluid sets up shear forces near the
boundary that reduce the flow speed to that of the boundary.
That fluid layer which has had its velocity affected by the
boundary shear is called the boundary layer.
For smooth upstream boundaries the boundary layer starts out as
a laminar boundary layer in which the fluid particles move in
smooth layers.
As the laminar boundary layer increases in thickness, it becomes
unstable and finally transforms into a turbulent boundary layer
in which the fluid particles move in haphazard paths.
When the boundary layer has become turbulent, there is still a
very thon layer next to the boundary layer that has laminar
motion. It is called the laminar sublayer.
6
Introduction
Fig. 6.2 The development of the boundary layer for flow over a
flat plate, and the different flow regimes. The vertical scale has
been greatly exaggerated and horizontal scale has been shortened.
7
Introduction
The turbulent boundary layer can be considered to consist
of four regions, characterized by the distance from the wall.
The very thin layer next to the wall where viscous effects
are dominant is the viscous sublayer. The velocity profile
in this layer is very nearly linear, and the flow is nearly
parallel.
Next to the viscous sublayer is the buffer layer, in which
turbulent effects are becoming significant, but the flow is
still dominated by viscous effects.
Above the buffer layer is the overlap layer, in which the
turbulent effects are much more significant, but still not
dominant.
Above that is the turbulent (or outer) layer in which
turbulent effects dominate over viscous effects.
8
Boundary layer thickness, δ
The boundary layer thickness is defined as the vertical
distance from a flat plate to a point where the flow velocity
reaches 99 per cent of the velocity of the free stream.
Another definition of boundary layer are the
➢ Boundary layer displacement thickness, δ*
➢ Boundary layer momentum thickness, θ
Boundary layer displacement thickness, δ*
Consider two types of fluid flow past a stationary horizontal
plate with velocity U as shown in Fig. 6.3. Since there is no
viscosity for the case of ideal fluid (Fig. 6.3a), a uniform
velocity profile is developed above the solid wall.
However, the velocity gradient is developed in the boundary
layer region for the case of real fluid with the presence of
viscosity and no-slip at the wall (Fig. 6.3b).
9
Boundary layer displacement thickness, δ*
Figure 6.3 Flow over a horizontal solid surface for the case of (a)
Ideal fluid (b) Real fluid
The velocity deficits through the element strip of cross
section b-b is U - u . Then the reduction of mass flow rate is
obtained as where b is the plate width.
The total mass reduction due to the presence of viscosity
compared to the case of ideal fluid
(6.1)
10
Boundary layer displacement thickness, δ*
However, if we displace the plate upward by a distance at
section a-a to give mass reduction of , then the deficit
of flow rates for the both cases will be identical if
(6.2)
11
Boundary layer displacement thickness, δ*
12
Boundary layer displacement thickness, δ*
The displacement thickness represents the vertical distance
that the solid boundary must be displaced upward so that the
ideal fluid has the same mass flow rate as the real fluid.
13
Boundary layer momentum thickness, θ
Another definition of boundary layer thickness, the boundary
layer momentum thickness θ, is often used to predict the drag
force on the object surface.
By referring to Fig. 6.3, again the velocity deficit through the
element strip of cross section b-b contributes to deficit in
momentum flux as
(6.3)
Thus, the total momentum reductions
(6.4)
Fig. 6.6.
Comparison of flow
past a sharp flat
plate at low and
high Reynolds
numbers: (a)
laminar, low-Re
flow; (b) high-Re
flow.
17
Reynolds Number and Geometry Effects
A high-Reynolds-number flow (Fig. 6.6b) is much more
amenable to boundary layer patching, as first pointed out by
Prandtl in 1904.
The viscous layers, either laminar or turbulent, are very thin,
thinner even than the drawing shows.
We define the boundary layer thickness δ as the locus of points
where the velocity u parallel to the plate reaches 99 percent of
the external velocity U.
The accepted formulas for flat-plate flow, and their approximate
ranges, are
(6.5)
18 .
Reynolds Number and Geometry Effects
where Rex = Ux/ν is called the local Reynolds number of the flow
along the plate surface. The turbulent flow formula applies for Rex
greater than approximately 106 .
Some computed values are shown below
20
Example 1
21
Boundary Layer: Momentum Integral Estimates
A shear layer of unknown thickness grows along the sharp flat plate
in Fig. 6.7. The no-slip wall condition retards the flow, making it
into a rounded profile u(x,y), which merges into the external
velocity U = constant at a “thickness” y = δ(x).
22
Boundary Layer: Momentum Integral Estimates
The drag force on the plate is given by the following momentum
integral across the exit plane:
(6.6)
where b is the plate width into the paper and the integration
is carried out along a vertical plane x = constant.
Equation (6.6) was derived in 1921 by Kármán, who wrote it
in the convenient form of the momentum thickness as:
(6.7)
23
Boundary Layer: Momentum Integral Estimates
(6.8)
(6.9)
(6.10)
(6.11)
(6.12)
25
Boundary Layer: Momentum Integral Estimates
(6.13)
(6.14)
26
Boundary Layer: Momentum Integral Estimates
This is only 10 percent higher than the known exact laminar-
plate-flow solution cf = 0.664/Rex1/2
The dimensionless quantity cf, called the skin friction coefficient,
is analogous to the friction factor f in ducts.
A boundary layer can be judged as “thin” if, say, the ratio δ/x is
less than about 0.1. This occurs at δ/x = 0.1 = 5.0/Rex1/2 or at Rex
= 2500.
For Rex less than 2500 we can estimate that boundary layer
theory fails because the thick layer has a significant effect on the
outer inviscid flow.
The upper limit on Rex for laminar flow is about 3 x106, where
measurements on a smooth flat plate show that the flow
undergoes transition to a turbulent boundary layer.
From 3 x106 upward the turbulent Reynolds number may be
arbitrarily large, and a practical limit at present is 5 x 1010 for oil
27
supertankers
Boundary Layer: Momentum Integral Estimates
For parallel flow over a flat plate, the pressure drag is zero,
and thus the drag coefficient is equal to the friction drag
coefficient, or simply the friction coefficient).
Once the average friction coefficient Cf is available, the
drag (or friction) force over the surface is determined from
28
Example 2
Are low-speed, small-scale air and water boundary layers really
thin? Consider flow at U =1 ft/s past a flat plate 1 ft long.
Compute the boundary layer thickness at the trailing edge for (a)
air and (b) water at 68F.
Solution
From Table νair = 1.61 E-4 ft2/s. The trailing-edge Reynolds
number thus is
Since this is less than 106, the flow is presumed laminar, and
since it is greater than 2500, the boundary layer is reasonably
thin. The predicted laminar thickness is
29
Example 2
30
Part II
Laminar and Turbulent Pipe Flow
31
Introduction
Fluid flow in circular and noncircular pipes is commonly
encountered in practice.
The hot and cold water that we use in our homes is pumped
through pipes. Water in a city is distributed by extensive
piping networks. Oil and natural gas are transported
hundreds of miles by large pipelines. Blood is carried
throughout our bodies by arteries and veins. The cooling
water in an engine is transported by hoses to the pipes in
the radiator where it is cooled as it flows.
The fluid in such applications is usually forced to flow by a
fan or pump through a flow section.
We pay particular attention to friction, which is directly
related to the pressure drop and head loss during flow
through pipes and ducts.
32
Introduction
The pressure drop is then used to determine the pumping
power requirement.
A typical piping system involves pipes of different diameters
connected to each other by various fittings or elbows to
route the fluid, valves to control the flow rate, and pumps to
pressurize the fluid.
The terms pipe, duct, and conduit are usually used
interchangeably for flow sections.
In general, flow sections of circular cross section are
referred to as pipes (especially when the fluid is a liquid),
and flow sections of noncircular cross section as ducts
(especially when the fluid is a gas) Small diameter pipes are
usually referred to as tubes.
33
Introduction
Most fluids, especially liquids, are
transported in circular pipes. This is
because pipes with a circular cross
section can withstand large pressure
differences between the inside and the
outside without undergoing significant
distortion.
Noncircular pipes are usually used in
applications such as the heating and
cooling systems of buildings where the
pressure difference is relatively small,
the manufacturing and installation
costs are lower, and the available space
is limited for ductwork.
34
Introduction
The fluid velocity in a pipe
changes from zero at the surface
because of the no-slip condition
to a maximum at the pipe center.
In fluid flow, it is convenient to
work with an average velocity
Vavg, which remains constant in
incompressible flow when the
cross-sectional area of the pipe is
constant.
The change in average velocity
due to change in density and
temperature and due to friction is
usually small and is thus
35 disregarded in calculations.
Introduction
36
LAMINAR AND TURBULENT FLOWS
Fluid flow in a pipe is streamlined at low
velocities but turns chaotic as the velocity is
increased above a critical value.
A laminar flow is characterized by smooth
streamlines and highly ordered motion, and
turbulent flow is characterized by velocity
fluctuations and highly disordered motion.
The transition from laminar to turbulent flow
does not occur suddenly; rather, it occurs over
some region in which the flow fluctuates between
laminar and turbulent flows before it becomes
fully turbulent.
Most flows encountered in practice are turbulent.
Laminar flow is encountered when highly viscous
fluids such as oils flow in small pipes or narrow
37 passages.
Reynolds Number
The transition from laminar to turbulent flow depends on the
geometry, surface roughness, flow velocity, surface temperature,
and type of fluid, among other things.
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
38
Reynolds Number
where Vavg = average flow velocity (m/s), D = characteristic
length of the geometry (diameter in this case, in m), and ν =
μ/ρ = kinematic viscosity of the fluid (m2/s).
Note that the Reynolds number is a dimensionless quantity.
At large Reynolds numbers, the inertial forces, which are
proportional to the fluid density and the square of the fluid
velocity, are large relative to the viscous forces, and thus
the viscous forces cannot prevent the random and rapid
fluctuations of the fluid.
At small or moderate Reynolds numbers, however, the
viscous forces are large enough to suppress these
fluctuations and to keep the fluid “in line.”
Thus the flow is turbulent in the first case and laminar in
the second.
39
Reynolds Number
The Reynolds number at which the flow becomes turbulent
is called the critical Reynolds number, Recr .
The value of the critical Reynolds number is different for
different geometries and flow conditions. For internal flow
in a circular pipe, the generally accepted value of the
critical Reynolds number is Recr = 2300.
For flow through noncircular pipes, the Reynolds number is
based on the hydraulic diameter Dh defined as
41
The Entrance Region
Consider a fluid entering a circular pipe at a uniform velocity.
Because of the no-slip condition, the fluid particles in the
layer in contact with the surface of the pipe come to a
complete stop.
This layer also causes the fluid particles in the adjacent layers
to slow down gradually as a result of friction.
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.
The hypothetical 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
42 essentially constant in the radial direction.
The Entrance Region
The thickness of this boundary layer increases in the flow
direction until the boundary layer reaches the pipe center and thus
fills the entire pipe.
The region from the pipe inlet to the point at which the boundary
layer merges at the centerline is called the hydrodynamic
entrance region, and the length of this region is called the
hydrodynamic entry length Lh.
43
The Entrance Region
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.
44
The Entrance Region
The velocity profile in the fully developed region is
parabolic in laminar flow and somewhat flatter (or fuller) in
turbulent flow due to eddy motion and more vigorous
mixing in the radial direction.
Entry Lengths
The hydrodynamic entry length is usually taken to be the
distance from the pipe entrance to where the wall shear
stress (and thus the friction factor) reaches within about 2
percent of the fully developed value.
In laminar flow, the hydrodynamic entry length is given
approximately as
45
Fig. The variation of wall shear stress in the flow direction for
flow in a pipe from the entrance region into the fully developed
46 region.
The Entrance Region
In turbulent flow, the intense mixing during random
fluctuations usually overshadows the effects of molecular
diffusion.
The hydrodynamic entry length for turbulent flow can be
approximated as [see Bhatti and Shah (1987) and Zhi-qing
(1982)]
48
Laminar Flow in Pipes
Flow in pipes is laminar for Re ≤ 2300, and that the flow is
fully developed if the pipe is sufficiently long (relative to
the entry length) so that the entrance effects are negligible.
In this section we consider the steady laminar flow of an
incompressible fluid with constant properties in the fully
developed region of a straight circular pipe.
In fully developed laminar flow, each fluid particle moves
at a constant axial velocity along a streamline and the
velocity profile u(r) remains unchanged in the flow
direction. There is no motion in the radial direction, and
thus the velocity component in the direction normal to flow
is everywhere zero. There is no acceleration since the flow
is steady and fully developed.
49
Laminar Flow in Pipes
Consider a ring-shaped
differential volume element of
radius r, thickness dr, and length
dx oriented coaxially with the
pipe, as shown in the Fig.
The volume element involves
only pressure and viscous effects
and thus the pressure and shear
forces must balance each other.
The pressure force acting on a
submerged plane surface is the
product of the pressure at the
centroid of the surface and the
surface area.
50
Laminar Flow in Pipes
A force balance on the volume element in the flow
direction gives
51
Laminar Flow in Pipes
52
Laminar Flow in Pipes
53
Laminar Flow in Pipes
The velocity profile u(r) is obtained by applying the
boundary conditions at r = 0 (because of
symmetry about the centerline) and u = 0 at r = R (the
no-slip condition at the pipe surface). We get
55
Pressure Drop and Head Loss
A quantity of interest in the analysis of pipe flow is the
pressure drop P since it is directly related to the power
requirements of the fan or pump to maintain flow.
We note that dP/dx = constant, and integrating from x = x1
where the pressure is P1 to x = x1 + L where the pressure is
P2 gives
56
Pressure Drop and Head Loss
Pressure drop due to viscous effects represents an
irreversible pressure loss, and it is called pressure loss PL
to emphasize that it is a loss.
Pressure drop is proportional to the viscosity μ of the fluid,
and P would be zero if there were no friction. Therefore,
the drop of pressure from P1 to P2 in this case is due
entirely to viscous effects.
In practice, it is found 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
57
Pressure Drop and Head Loss
where ρV2avg/2 is the dynamic pressure and f is the Darcy
friction factor also called the Darcy–Weisbach friction
factor
59
Pressure Drop and Head Loss
60
Pressure Drop and Head Loss
61
Inclined Pipes
Relations for inclined
pipes can be obtained in
a similar manner from a
force balance in the
direction of flow. The
only additional force in
this case is the
component of the fluid
weight in the flow
direction, whose
magnitude is
63
Inclined Pipes
It can also be shown that the average velocity and the volume
flow rate relations for laminar flow through inclined pipes
are, respectively,
64
Inclined Pipes
67
68
69
EXAMPLE 2. Pressure Drop and Head Loss in a Pipe
Water at 40°F (ρ = 62.42 lbm/ft3 and μ = 1.038 x 10-3 lbm/ft · s)
is flowing through a 0.12-in (= 0.010 ft) diameter 30-ft-long
horizontal pipe steadily at an average velocity of 3.0 ft/s.
Determine (a) the head loss, (b) the pressure drop, and (c) the
pumping power requirement to overcome this pressure drop.
70
71
72
Turbulent Flow in Pipes
Most flows encountered in engineering practice are turbulent,
and thus it is important to understand how turbulence affects
wall shear stress.
However, turbulent flow is a complex mechanism dominated
by fluctuations, and despite tremendous amounts of work
done in this area by researchers, the theory of turbulent flow
remains largely undeveloped.
Therefore, we must rely on experiments and the empirical or
semi-empirical correlations developed for various situations.
Turbulent flow is characterized by random and rapid
fluctuations of swirling regions of fluid, called eddies,
throughout the flow. These fluctuations provide an additional
mechanism for momentum and energy transfer.
73
Turbulent Flow in Pipes
In laminar flow, fluid particles flow in an orderly manner
along pathlines, and momentum and energy are transferred
across streamlines by molecular diffusion.
In turbulent flow, the swirling eddies transport mass,
momentum, and energy to other regions of flow much more
rapidly than molecular diffusion, greatly enhancing mass,
momentum, and heat transfer.
As a result, turbulent flow is associated with much higher
values of friction, heat transfer, and mass transfer
coefficients.
The eddy motion in turbulent flow causes significant
fluctuations in the values of velocity, temperature, pressure,
and even density (in compressible flow).
74
Turbulent Flow in Pipes
Fluctuations of the velocity component u with time at a
specified location in turbulent flow shown in fig below.
75
Turbulent Flow in Pipes
76
Turbulent Flow in Pipes
The total shear stress can be expressed conveniently as
77
Turbulent Flow in Pipes
The velocity gradients at the
wall, and thus the wall shear
stress, are much larger for
turbulent flow than they are
for laminar flow.
78
Turbulent Flow in Pipes
Turbulent Velocity Profile
The velocity profile is parabolic
in laminar flow but is much
fuller in turbulent flow, with a
sharp drop near the pipe wall.
Turbulent flow along a wall can
be considered to consist of four
regions, characterized by the
distance from the wall. The very
thin layer next to the wall where
viscous effects are dominant is
the viscous (or laminar or
linear or wall) sublayer.
79
Turbulent Flow in Pipes
The velocity profile in this layer is very nearly linear, and the
flow is streamlined.
Next to the viscous sublayer is the buffer layer, in which
turbulent effects are becoming significant, but the flow is still
dominated by viscous effects.
Above the buffer layer is the overlap (or transition) layer,
also called the inertial sublayer, in which the turbulent effects
are much more significant, but still not dominant.
Above that is the outer (or turbulent) layer in the remaining
part of the flow in which turbulent effects dominate over
molecular diffusion (viscous) effects.
Flow characteristics are quite different in different regions, and
thus it is difficult to come up with an analytic relation for
the velocity profile for the entire flow as we did for laminar
80
flow.
Turbulent Flow in Pipes
Numerous empirical velocity profiles exist for turbulent
pipe flow. Among those, the simplest and the best known is
the power-law velocity profile expressed as
81
Turbulent Flow in Pipes
The turbulent velocity
profile is fuller than the
laminar one, and it becomes
more flat as n (and thus the
Reynolds number)
increases.
The power-law profile
cannot be used to calculate
wall shear stress since it
gives a velocity gradient of
infinity there, and it fails to
give zero slope at the
centerline.
82
The Moody Chart
The friction factor in fully developed turbulent pipe flow
depends on the Reynolds number and the relative
roughness ε/D, which is the ratio of the mean height of
roughness of the pipe to the pipe diameter.
Colebrook equation
85
The Moody Chart
86
Turbulent Flow in Pipes
The Colebrook equation is implicit in f, and thus the
determination of the friction factor requires some iteration
unless an equation solver is used.
An approximate explicit relation for f was given by S. E.
Haaland in 1983 as
87
Turbulent Flow in Pipes
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
88 head loss)
Turbulent Flow in Pipes
Problems of the first type are straightforward and can be
solved directly by using the Moody chart.
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
unless an equation solver is used.
To avoid tedious iterations in head loss, flow rate, and
diameter calculations, Swamee and Jain proposed the
following explicit relations in 1976 that are accurate to
within 2 percent of the Moody chart:
89
Turbulent Flow in Pipes
90
EXAMPLE 3. Determining the Head Loss in a Water Pipe
91
which is greater than 4000. Therefore, the flow is turbulent. The
92
relative roughness of the pipe is calculated using the Table
93
The friction factor could also be determined easily from the
explicit Haaland relation. It would give f = 0.0172, which
is sufficiently close to 0.0174.
94
Minor Losses
The fluid in a typical piping system passes through various
fittings, valves, bends, elbows, tees, inlets, exits,
enlargements, and contractions in addition to the pipes.
These components interrupt the smooth flow of the fluid
and cause additional losses because of the flow separation
and mixing they induce.
In a typical system with long pipes, these losses are minor
compared to the total head loss in the pipes (the major
losses) and are called minor losses.
Although this is generally true, in some cases the minor
losses may be greater than the major losses. This is the
case, for example, in systems with several turns and valves
in a short distance.
95
Minor Losses
The head loss introduced by a completely open valve, for
example, may be negligible. But a partially closed valve may
cause the largest head loss in the system, as evidenced by the
drop in the flow rate.
Flow through valves and fittings is very complex, and a
theoretical analysis is generally not plausible.
Therefore, minor losses are determined experimentally, usually
by the manufacturers of the components.
Minor losses are usually expressed in terms of the loss
coefficient KL (also called the resistance coefficient), defined
as
101
→ α1 and α2 are kinetic energy correction factors
102
103
Example 5. Determine Head Loss
As shown in Fig. below, crude oil at 140 °F with γ =53.7 lb/ft3
and μ = 8 x 105 lb . s ft2 (about four times the viscosity of
water) is pumped across Alaska through the Alaskan pipeline, a
799-mile-long, 4-ft-diameter steel pipe, at a maximum rate of
Q = 2.4 million barrels day = 117 ft3 /s. Determine the
horsepower needed for the pumps that drive this large system.
104
Solution
From the energy equation we obtain
105
106
Example 6. Minor losses
Water at 10°C flows from a large reservoir to a smaller
one through a 5-cm diameter cast iron piping system,
as shown in Fig. below. Determine the elevation z1 for
a flow rate of 6 L/s.
107
108
109
Example 7
A horizontal pipe has an abrupt expansion from D1= 8 cm
to D2 = 16 cm. The water velocity in the smaller section is
10 m/s and the flow is turbulent. The pressure in the
smaller section is P1 = 300 kPa. Taking the kinetic energy
correction factor to be 1.06 at both the inlet and the outlet,
determine the downstream pressure P2, and estimate the
error that would have occurred if Bernoulli’s equation had
been used.
110
111
112
PIPING NETWORKS
Most piping systems encountered in practice such as the
water distribution systems in cities or commercial or
residential establishments involve numerous parallel and
series connections.
Pipes in Series
When the pipes are connected in series, the flow rate
through the entire system remains constant regardless of the
diameters of the individual pipes in the system. This is a
natural consequence of the conservation of mass principle
for steady incompressible flow.
The total head loss in this case is equal to the sum of the
head losses in individual pipes in the system, including the
minor losses.
113
PIPING NETWORKS
For a pipe that branches out into two (or more) parallel pipes
and then rejoins at a junction downstream, the total flow rate
is the sum of the flow rates in the individual pipes. The
pressure drop (or head loss) in each individual pipe connected
in parallel must be the same since and the
junction pressures PA and PB are the same for all the individual
pipes.
114
PIPING NETWORKS
115
PIPING NETWORKS
116
PIPING NETWORKS
Another type of multiple pipe system called a loop is shown in
Fig. In this case the flowrate through pipe (1) equals the sum
of the flowrates through pipes (2) and (3), or
for fluid that travels through pipes (1) and (3). These can be
combined to give This is statement of the fact that fluid
particles that travel through pipe (2) and particles that travel
through pipe (3) all originate from common conditions at the
junction (or node, N) of the pipes and all end up at the same
final conditions.
118
DIMENSIONAL ANALYSIS
(OPEN CHANNEL FLOW AND HYDRAULIC
MACHINERY)
UNIT – II
In this regard, dimensional analysis is an important tool that helps in correlating
analytical results with experimental data and to predict the prototype behavior from
the measurements on the model.
Dimensions and Units
In dimensional analysis we are only concerned with the nature of the
dimension i.e. its quality not its quantity.
maximum
The values of a, bthree
and c areor fourby variables
obtained only.
comparing the powers of the fundamental
dimension on both sides.
Rayeligh’s Method
Methodology:
Problem: Find the expression for Discharge Q in a open channel flow when Q is
To
depends on Arearelationship
define A and Velocity V. among variables
Let X is a function of X1 ,X2, X3 and mathematically it can be written as
X = f(X1, X2, X3)
Solution:
This method is used for determining the
Q = K.Aa.Vb 1
This can
where K is
beaalso
Non-dimensional
written as constant
expression for a variable which depends upon
X = K (X1a the
Substitute , X2bdimensions
, X3c ) whereon K
both
is constant
sides of equation
and a, b and
1 c are arbitrarily powers
M0 L3 T-1 = K. (L2)a.(LT-1)b
maximum three or four variables only.
The values
Equating powers
of a, bofand
M, cL,are
T on
obtained
both sides,
by comparing the powers of the fundamental
of T,on both sides. -1 = -b b=1
dimension
Power
Power of L, 3= 2a+b 2a = 2-b = 2-1 = 1
Substituting values of a, b, and c in Equation 1m
Q = K. A1. V1 = V.A
Rayeligh’s Method
Methodology:
To define relationship among variables
Let X is a function of X1 ,X2, X3 and mathematically it can be written as
X = f(X1, X2, X3)
This method is used for determining the
This can be also written as
expression
X = K (X , X , X for
1
a
2
b a variable
) where
3
c which
K is constant and a, b anddepends upon
c are arbitrarily powers
maximum
The values of a, bthree
and c areor fourby variables
obtained only.
comparing the powers of the fundamental
dimension on both sides.
Rayeligh’s Method
Methodology:
To define relationship among variables
Let X is a function of X1 ,X2, X3 and mathematically it can be written as
X = f(X1, X2, X3)
This method is used for determining the
This can be also written as
expression
X = K (X , X , X for
1
a
2
b a variable
) where
3
c which
K is constant and a, b anddepends upon
c are arbitrarily powers
maximum
The values of a, bthree
and c areor fourby variables
obtained only.
comparing the powers of the fundamental
dimension on both sides.
Buckingham’s π-Theorem
This method of analysis is used when number of variables are more.
Theorem:
If there are n variables in a physical phenomenon and those n variables contain m dimensions,
then variables can be arranged into (n-m) dimensionless groups called Φ terms.
Explanation:
If f (X1, X2, X3, ……… Xn) = 0 and variables can be expressed using m dimensions then
f (π1, π2, π3, ……… πn - m) = 0 where, π1, π2, π3, … are dimensionless groups.
Each π term contains (m + 1) variables out of which m are of repeating type and one is of non-
repeating type.
Each π term being dimensionless, the dimensional homogeneity can be used to get each π term.
Types of Similarities:
1. Geometric Similarity Length, Breadth, Depth, Diameter, Area,
Volume etc.,
2. Kinematic Similarity Velocity, Acceleration etc.,
3. Dynamic Similarity Time, Discharge, Force, Pressure Intensity,
Torque, Power
Geometric Similarity
The geometric similarity is said to be exist between the model and
prototype if the ratio of all corresponding linear dimensions in the model
and prototype are equal.
L P
B
P
D
P
Lr A P
Lr
2 V P
Lr
3
L m B
m D
m A m V m
V P
Vr a P
ar
V m a m
Also the directions of the velocities in the model and prototype should be same
Dynamic Similarity
The dynamic similarity is said exist between model and prototype if the
ratios of corresponding forces acting at the corresponding points are
equal
FP
Fr
Fm
Re Re
m p tr
Time Scale Ratio L r
Velocity, V = Length/Time T = L/V
V r
m p
t r
Fe Fe
m p
V m
V p
V r
VelocityScale Ratio L
r
g L
m m gL
p P
T Scale Ratio for Accele ration 1 Scale Ratio for Pressure Intensity Lr
r
F r
2.5
Qr Scale Ratio for Discharge Lr P r
Scale Ratio for Power Lr
3.5
Problems
1. In 1 in 40 model of a spillway, the velocity and discharge are 2 m/s and 2.5
m3/s. Find corresponding velocity and discharge in the prototype