Kuliah 4 5 Dan 6 Internal Incompressible Viscous Flow
Kuliah 4 5 Dan 6 Internal Incompressible Viscous Flow
VISCOUS FLOW
Nazaruddin Sinaga
Flow Measurements
2
3
4
Osborne Reynolds Experiment
5
6
Laminar and Turbulent Flow
7
Entrance Length
8
Developing Flow
9
10
11
12
13
14
15
16
2
1 5
17
18
19
20
21
22
23
24
25
26
27
Types of Flow
• The physical nature of fluid flow can be categorized into three types, i.e.
laminar, transition and turbulent flow. Reynolds Number (Re) can be
used to characterize these flow.
VD VD
Re (6.3)
where = density
= dynamic viscosity
= kinematic viscosity ( = /)
V = mean velocity
D = pipe diameter
du
dy
29
Pressure Loss During A Laminar Flow In A Pipe
ur
R 2 r 2 P
(6.7)
4 L
R 2 P
u max
4 L
It can be shown that the mean velocity is half the maximum velocity,
i.e. V=umax/2
Shear stress and velocity distribution in pipe for laminar flow
35
Fully Developed Laminar Flow
in a Pipe
• Volume Flow Rate
36
Fully Developed Laminar Flow
in a Pipe
• Average Velocity
Maximum Velocity
37
38
Turbulent Flow
This is the most commonly occurring flow in engineering
practice in which fluid particles move erratically causing
instantaneous fluctuations in the velocity components.
These fluctuations cause additional shear stresses. In this
type of flow both viscous and turbulent shear stresses exists.
Thus, the shear stress in turbulent flow is a combination of
laminar and turbulent shear stresses, and can be written as:
dU
lam tur
dy
where = dynamic viscosity
= eddy viscosity which is not a fluid property but depends upon
turbulence condition of flow.
39
40
Velocity Profile for Fully Developed
Pipe Flow
41
The velocity at any point in the cross-section will be proportional
to the one-seventh power of the distance from the wall, which
can be expressed as:
Uy 1/ 7
y
(6.10)
U CL R
U yU *
5.75 log10 5.5 (6.11a)
U*
y is the distance form the pipe wall, k is the surface roughness and
is the kinematic viscosity of the fluid.
44
Turbulent Velocity Profiles in Fully
Developed Pipe Flow
45
46
47
Head Losses
48
Head Losses
The momentum balance in the flow direction is thus given by
49
50
51
• Head Loss
53
54
Friction Factor
55
56
57
Nikuradse’s Experiments
60
61
Calculation of Minor Head Loss
• Minor Losses
– Examples: Inlets and Exits; Enlargements and Contractions;
Pipe Bends; Valves and Fittings
62
Pipe Entrance
• Developing flow
– Includes boundary layer
and core,
– viscous effects grow inward
from the wall
• Fully developed flow
– Shape of velocity profile is Entrance length Le Fully developed
Pressure flow region
same at all points along
pipe Entrance
pressure drop
Region of linear
pressure drop
• Where,
V2 L
hL h f hm f
D K m
m 2g m
66
67
68
69
70
71
72
73
74
Example 1
Water at 10C is flowing at a rate of 0.03 m3/s through a pipe. The pipe
has 150-mm diameter, 500 m long, and the surface roughness is estimated
at 0.06 mm. Find the head loss and the pressure drop throughout the
length of the pipe.
Solution:
From Table 1.3 (for water): = 1000 kg/m3 and =1.30x10-3 N.s/m2
V = Q/A and A=R2
A = (0.15/2)2 = 0.01767 m2
V = Q/A =0.03/.0.01767 =1.7 m/s
Re = (1000x1.7x0.15)/(1.30x10-3) = 1.96x105 > 2000 turbulent flow
To find , use Moody Diagram with Re and relative roughness (k/D).
k/D = 0.06x10-3/0.15 = 4x10-4
From Moody diagram, 0.018
The head loss may be computed using the Darcy-Weisbach equation.
L V2 500 x 1.7 2
hf 0.018 x 8.84m.
D 2g 0.15 x 2 x 9.81
The pressure drop along the pipe can be calculated using the relationship:
ΔP=ghf = 1000 x 9.81 x 8.84
ΔP = 8.67 x 104 Pa
75
Example 2
Determine the energy loss that will occur as 0.06 m3/s water flows
from a 40-mm pipe diameter into a 100-mm pipe diameter through
a sudden expansion.
Solution:
The head loss through a sudden enlargement is given by;
2
V
hm K a
2g
Q 0.06
Va 3.58 m / s
Aa (0.04 / 2) 2
76
Example 3
Calculate the head added by the pump when
the water system shown below carries a
discharge of 0.27 m3/s. If the efficiency of
the pump is 80%, calculate the power input
required by the pump to maintain the flow.
77
Solution:
Applying Bernoulli equation between section 1 and 2
P1 V12 P2 V2 2 (1)
z1 Hp z2 H L12
g 2g g 2g
V2 1000
H L1 2 0.015x 0 .5 0 .4 1
2g 0.4
V2
39.4
2g
V2
H p 230 200 39.4
From (2): 2 x 9.81
The velocity can be calculated using the
continuity equation:
Q 0.27
V 2.15 m / s
A 0.4 / 22
gQH p
p
Pin
gQH p 1000 x 9.81x 0.27 x39.3
Pin
p 0.8
80
Pipes in Series
81
Referring to Figure 6.11, the Bernoulli equation can be
written between points 1 and 2 as follows;
P1 V12 P2 V2 2
z1 z2 H L12 (6.18)
g 2g g 2g
Q1 = Q2 83
Pipes in Parallel
• A combination of two
or more pipes
connected between
two points so that the
discharge divides at
the first junction and
rejoins at the next is
known as pipes in Figure 6.12 Pipelines in
parallel
parallel. Here the
head loss between the
two junctions is the
same for all pipes.
84
Applying the continuity equation to the system;
Q1 = Qa + Qb = Q2 (6.19)
86
Pipe Network
The earliest systematic method of network analysis
(Hardy-Cross Method) is known as the head balance or
closed loop method. This method is applicable to
system in which pipes form closed loops. The outflows
from the system are generally assumed to occur at the
nodes junction.
For a given pipe system with known outflows, the
Hardy-Cross method is an iterative procedure based on
initially iterated flows in the pipes. At each junction
these flows must satisfy the continuity criterion, i.e.
the algebraic sum of the flow rates in the pipe meeting
at a junction, together with any external flows is zero.
87
• Assigning clockwise flows and their associated
head losses are positive, the procedure is as
follows:
Assume values of Q to satisfy Q = 0.
Calculate HL from Q using HL = K1Q2 .
If HL = 0, then the solution is correct.
88
If HL 0, then apply a correction factor,
Q, to all Q and repeat from step (2).
For practical purposes, the calculation is
usually terminated when HL < 0.01 m or Q
< 1 L/s.
A reasonably efficient value of Q for rapid
convergence is given by;
Q
H L
H
(6.21)
2 L
Q
89
Example 4
• A pipe 6-cm in diameter, 1000m long and with = 0.018 is
connected in parallel between two points M and N with
another pipe 8-cm in diameter, 800-m long and having =
0.020. A total discharge of 20 L/s enters the parallel pipe
through division at A and rejoins at B. Estimate the
discharge in each of the pipe.
90
Solution:
Continuity: Q = Q1 + Q2
(1)0.02 (0.06) 2 V12 (0.08) 2 V22
4 4
V1 1.778V2 7.074
0.018x1000 2 0.020x800 2
V1 V2
0.06 0.08
V1 0.8165V2 ( 2)
91
Q 2 A 2 V2 (0.08) 2 x 2.73
4
Q2 = 0.0137 m3/s
From (2):
V1 = 0.8165 V2 = 0.8165x2.73 = 2.23 m/s
Q1 = 0.0063 m3/s
92
Example 6.6
• For the square loop shown, find the discharge in all the
pipes. All pipes are 1 km long and 300 mm in diameter, with
a friction factor of 0.0163. Assume that minor losses can be
neglected.
93
Solution:
Assume values of Q to satisfy continuity equations all at nodes.
The head loss is calculated using; HL = K1Q2
HL = hf + hLm
But minor losses can be neglected: hLm = 0
Thus HL = hf
Head loss can be calculated using the Darcy-Weisbach equation
L V2
hf
D 2g
94
L V2
HL hf
D 2g
1000 V2
H L 0.0163 x x
0.3 2 x 9.81
Q2 Q2
H L 2.77 2.77 x
A2 2
2
x 0.3
4
H L 554Q 2
H L K' Q 2
K ' 554
First trial
AB 60 2.0 0.033
BC 40 0.886 0.0222
CD 0 0 0
AD -40 -0.886 0.0222
2.00 0.0774
Since H L
> 0.01 m, then correction has to be applied.
95
HL 2
Q 12.92 L / s
2 H L 2 x 0. 0774
Q
Second trial
96
Solution of Pipe Flow Problems
• Single Path
– Find p for a given L, D, and Q
Use energy equation directly
Terima kasih
104
Shear Stress in Pipes
• Steady, uniform flow in a pipe: momentum
flux is zero and pressure distribution across
pipe is hydrostatic, equilibrium exists
between pressure, gravity and shear forces
dp
Fs 0 pA ( p s ) A W sin 0 (D ) s
ds
dp dz
0 sA As 0 (D )s
ds ds
D d p
0 [ ( z )]
4 ds
D dh
0
4 ds
4 L 0
h1 h2 h f
D
0 F ( , V , , D, e)
4L
4 F ( 1 , 2 ) hf 0
D
Repeating variables : , V , D
4L e
e V 2 F (Re, )
1 Re; 2 ; 3 0 2 D D
D V
L V2 e
0 e 8F (Re, D )
F (Re, ) D 2g
V 2 D
L V2 e
e hf f f 8F (Re, )
0 V 2 F (Re, ) D 2g D
D
dV r dh
dy 2 ds
dV dV
dy dr
dV r dh
dr 2 ds
r dh
dV dr
2 ds
r 2 dh r 2 dh
V C C 0
4 ds 4 ds
2
r 2 dh r
V 0 1
• Velocity distribution in a pipe (laminar
4 ds r0
flow) is parabolic with maximum at
r 2 center.
V Vmax 1
r0
Discharge in Laminar Flow
dh 2
V ( r0 r 2 )
4 ds
dh 2
Q VdA 0r0 ( r0 r 2 )( 2rdr )
4 ds
r0
dh ( r 2 r02 ) 2
4 ds 2
0
r04 dh
Q
8 ds
D 4 dh
128 ds
Q
V
A
D 2 dh
V
32 ds
Head Loss in Laminar Flow
D 2 dh 32 LV
V hf
32 ds
D 2
dh 32
V 32 LV V 2 / 2
ds D 2
D 2 V 2 / 2
32
dh V ds L
D 2 64( )( ) V 2 / 2
V D D
32
h2 h1 V ( s2 s1 ) 64 L
D 2 ( ) V 2 / 2
Re D
h1 h2 h f
L V 2 64
hf f f
D 2 Re
32 LV
hf
D 2
EGL & HGL for Losses in a Pipe
• Entrances, bends, and other flow transitions
cause the EGL to drop an amount equal to
the head loss produced by the transition.
• EGL is steeper at entrance than it is
downstream of there where the slope is
equal the frictional head loss in the pipe.
• The HGL also drops sharply downstream of
an entrance
Ex. (10.8)
Given: Kerosene (S=0.94, =0.048 N-s/m2).
Horizontal 5-cm pipe. Q=2x10-3 m3/s.
Find: Pressure drop per 10 m of pipe.
Solution:
V2 p V2 p
α1 1 1 z1 hL α2 2 2 z2
2g γ 2g γ
32 μLV
hL
γD 2
32 μLV V22
0 0 0.5 α2 00
γD 2 2g
α2 2 32 μL
V2 V 0 .5 0
2g γD 2
2 2 32 * 4 * 10 5 * 10
V2 V 0 .5 0
2g 0.8 * 62.4 * (1 / 32) 2
V22 8.45V 16.1 0
V 1.60 ft / s
0.8 * 1.94 * 1.6 * (0.25 / 12)
Re 1293 (laminar)
4 * 10 5
Q V * A 1.6 * * (0.25/12)2 / 4 1.23 * 10 3 cfs
Ex. (10.34)
Given: Glycerin@ 20oC flows commercial
steel pipe.
Find: h
Solution: 12,300 N / m, 0.62 Ns / m 2
V2 p V2 p
α1 1 1 z1 hL α2 2 2 z 2
2g γ 2g γ
p1 p
z1 hL 2 z 2
γ γ
p p
h 1 z1 ( 2 z 2 ) hL
γ γ
VD VD 0.6 * 0.02
Re 23.5 (laminar)
5.1 * 10 4
32 μLV 32(0.62)(1)(0.6)
h h L 2.42 m
γD 2 12,300 * (0.02) 2
Ex. (10.43)
Given: Figure
Find: Estimate the elevation required in the upper
reservoir to produce a water discharge of 10 cfs
in the system. What is the minimum pressure in
the pipeline and what is the pressure there?
Solution:
V2 p V2 p V2 p V2 p
α1 1 1 z1 hL α2 2 2 z 2 α1 1 1 z1 hL αb b b zb
2g γ 2g γ 2g γ 2g γ
0 0 z1 hL 0 0 z 2 V2 p
0 0 z1 hL 1 * b b zb
L V 2 2g γ
hL K e 2 K b K E f
D 2g pb V2 L V 2
L 430 z1 zb b K e K b f
K e 0.5; K b 0.4 (assumed); K E 1.0; f 0.025 * 10.75 γ 2g D 2g
D 1
300 12.732
Q 10 133 110.7 1.0 0.5 0.4 0.025
V 12.73 ft / s 1 2 * 32.2
A / 4 * 12
1.35 ft
12.732
z1 100 0.5 2 * 0.4 1.0 10.75 133 ft pb 62.4 * ( 1.53) 0.59 psig
2 * 32.2
VD 12.73 * 1
Re 9 * 105
5
1.14 * 10
Ex. (10.68)
4
Given: Commercial steel pipe to carry 300 cfs of Relative roughness: k s 1.5 x10 0.00002
water at 60oF with a head loss of 1 ft per 1000 ft D 8.06
of pipe. Assume pipe sizes are available in
even sizes when the diameters are expressed in
inches (i.e., 10 in, 12 in, etc.). Get better estimate of f
Find: Diameter. VD
Re
Solution: 1.22 x10 5 ft 2 / s; k 1.5 x10 4 ft
s
Q
D
Assume f = 0.015 ( / 4) D 2 Q
( / 4) D
L V2
hf f 300
D 2g Re 3.9 x106
5
( / 4)(8.06)1.22 x10
1000 (Q /( / 4) D 2 ) 2
1 0.015 *
D 2g f=0.010
33,984 22,656
1 1
D5 D5
D 8.06 ft D 7.43 ft 89 in.
Use a 90 in pipe
Ex. (10.81)
Given: The pressure at a water main is 300 kPa Assume f = 0.020
gage. What size pipe is needed to carry water
1/ 5
from the main at a rate of 0.025 m3/s to a fL Q 2
factory that is 140 m from the main? Assume D 8
h f 2g
galvanized-steel pipe is to be used and that the
pressure required at the factory is 60 kPa gage 1/ 5
0.02 140 (0.025) 2
at a point 10 m above the main connection. 8 0.100 m
14.45 2
9.81
Find: Size of pipe.
Solution: k s 0.15
Relative roughness: 0.0015
D 100
L V2 L (Q /( / 4) D 2 ) 2
hf f f
D 2g D 2g Friction factor: f 0.022
1/ 5
fL Q 2 1/ 5
D 8 0.022
h f 2g D 0.100 0.102 m
0.020
V2 p V2 p
α1 1 1 z1 hL α2 2 2 z2
2g γ 2g γ Use 12 cm pipe
300,000 60,000
hf 10
9810 9810
h f 14.45 m
Ex. (10.83)
Given: The 10-cm galvanized-steel pipe is 1000 m
long and discharges water into the atmosphere.
The pipeline has an open globe valve and 4
threaded elbows; h1=3 m and h2 = 15 m.
Find: What is the discharge, and what is the
pressure at A, the midpoint of the line?
Solution:
V12 p1 V22 p2
α1 z1 hL α2 z2
2g γ 2g γ
L V2
0 0 12 (1 K e K v 4 K b f ) 00 V2 p V2 p
D 2g α A A A z A h L α2 2 2 z 2
2g γ 2g γ
D = 10-cm and assume f = 0.025
pA L V2
1000 2 15 ( 2 K b f )
24 g (1 0.5 10 4 * 0.9 0.025 )V γ D 2g
0.1
pA 500 (0.942) 2
24 g ( 2 * 0.9 0.025 ) 15 9.6 m
V2 γ 0 .1 2g
265.1
V 0.942 m / s p A 9810 * ( 9.26) 90.8 kPa
So f = 0.025
Ex. (10.95)
Given: If the deluge through the system shown is 2
cfs, what horsepower is the pump supplying to
the water? The 4 bends have a radius of 12 in
and the 6-in pipe is smooth.
Find: Horsepower
Solution:
V2 p V2 p
α1 1 1 z1 h p α2 2 2 z2 hL
2g γ 2g γ So f = 0.0135
V2 L 1700
0 0 30 h p 0 60 2 (1 0.5 4 K b f ) h p 60 30 1.611(1 0.5 4 * 0.19 0.0135 )
2g D (1 / 2)
Q 2 107.6 ft
V 10.18 ft / s
A ( / 4)(1 / 2) 2 Qh p
p 24.4 hp
V22 550
1.611 ft
2g
VD 10.18 * (1 / 2)
Re 4.17 x105
5
1.22 x10