7/9/2013

TOPIC 2: FLOW IN PIPES AND

CHANNELS

OBJECTIVES

1. Calculate the friction factor for a pipe using the

Colebrook-White equation.

2. Undertake head loss, discharge and sizing

calculations for single pipelines.

3. Use head-loss vs discharge relationships to

calculate flow in pipe networks.

4. Relate normal depth to discharge for uniform

flow in open channels.

FLOW REGIMES
laminar

VD
Re 
turbulent ν

V = average (or bulk) velocity

D = diameter

For a pipe, Recrit ≈ 2300

1
 7/9/2013

DEVELOPMENT LENGTH

Ldevel 0.06 Re (laminar )

 1/6
D 4.4 Re (turbulent)

PIPE FLOW: BALANCE OF FORCES

l
 pressure
p 
r direct
ion of
 gravity
flow
 friction
z
p+p
mg


p( πr 2 )  ( p  Δp)(πr 2 )  mg
 sin θ  τ
(2πr Δl )  0
    
pressure force weight friction

m  ρπr Δl 2

  Δp(πr 2 )  ρπr 2 g Δz  τ(2πr Δl )  0
Δz  
sin θ   
Δl   Δ( p  ρgz )(πr 2 )  τ(2πr Δl )  0
Δ( p  ρgz ) τ
 2 0
Δl r
1 Δp*
p*  p  ρgz τ r
2 Δl

(Downstream) pressure Gradient: G  

Δp *

dp* τ   12 Gr
Δl dl

LAMINAR PIPE FLOW

Balance of forces: τ   12 Gr (stress  pressure gradient)

du
Viscous stress: τμ (stress  velocity gradient)
dr

du

1G
r R
dr 2μ r
Gr 2
u  constant u  0 on r  R
4μ

G 2
u (R  r 2 )
4μ

r2 GR2
u  u0 (1  ) u0 
R2 4μ

2
 7/9/2013

EXAMPLE, PAGE 4
G 2
u (R  r 2 )
4μ
Δp * ρgh f
Pressure gradient: G 
L L

Find, from the velocity distribution given above:

(a) the centreline velocity, u0;

(b) the average velocity V;

(c) the flow rate, Q, in terms of head loss and pipe diameter;
L V2
(d) the friction factor λ, defined by h f  λ ( )
D 2g
as a function of Reynolds number, Re.

QUESTIONS
Which forces are in balance in steady pipe flow?

Pressure, gravity, friction

How can one combine the effects of pressure and weight?

Via piezometric pressure p* = p + ρgz

How do we convert between pressure and head?

p
p = ρgh h
ρg
How do we define (a) dynamic pressure; (b) dynamic head?
V2
2 ρV
1 2

2g

How do we define the skin-friction coefficient?

τw
cf  1
2 ρV 2

PIPE FLOW: BALANCE OF FORCES

L
 pressure
w
p
 gravity
directi
on of  friction
flow
z
 p+p
πD 2
 Δp    sin θ  τ
mg w  πDL  0
mg  4  
 weight wall friction
net pressure force

πD 2 
mρ L
4   πD πD
2 2

  Δp   ρgΔz   τ w  πDL
Δz  4 4
sin θ  
L 
πD 2
 Δ( p  ρgz )   τ w  πDL
4
L
p*  p  ρgz  Δp*  4 τw
D

L 1
Definition of skin-friction coefficient: τ w  c f ( 12 ρV 2 ) Δp *  4c f ( ρV 2 )
D 2

3
 7/9/2013

DARCY-WEISBACH EQUATION
L 1
Δp *  λ ( ρV 2 )
D 2

L
Pressure loss due to friction = λ  dynamic pressure
D

L V2
hf  λ ( )
D 2g
L
Head loss due to friction = λ  dynamic head
D

λ  4c f

hf = frictional head loss

λ = friction factor
L = length of pipe
D = diameter
V = average velocity ( Q)

EXAMPLE SHEET, Q2

A 0.75 m diameter pipe carries 0.6 cumec.

At point A, elevation 40 m, a Bourdon gauge fitted

to the pipe records 1.75 bar, while at point B,
elevation 34 m and 1.5 km along the pipe a
similar gauge reads 2.1 bar.

Determine the flow direction and calculate the

friction factor.

CALCULATING THE FRICTION FACTOR

L V2
Defined by: hf  λ ( )
D 2g

64
Laminar flow: λ
Re

Turbulent flow – two limits:

1 Re λ
Smooth:  2.0 log10
λ 2.51
1 3.7 D
Rough:  2.0 log10
λ ks

1  k 2.51 
Colebrook-White Equation:  2.0 log10  s  
λ  3.7 D Re λ 

4
 7/9/2013

TYPICAL ROUGHNESS FOR COMMERCIAL

PIPES

Material ks (mm)
Riveted steel 0.9 – 9.0
Concrete 0.3 – 3.0
Wood stave 0.18 – 0.9
Cast iron 0.26
Galvanised iron 0.15
Asphalted cast iron 0.12
Commercial steel or wrought iron 0.046
Drawn tubing 0.0015
Glass 0 (smooth)

MOODY CHART
0.10
0.09 ks/D
0.08
0.05
0.07
Transition

0.06 0.03
0.05 0.02
Laminar
0.04 l = 64/Re
0.01
l

0.03 0.005

0.002
0.02 0.001

0.0004
0.0002
0.0001
smooth-walled limit
0.00001 0.00005
0.01
1.0E02 1.0E03 1.0E04 1.0E05 1.0E06 1.0E07
Re = VD/n

OTHER LOSSES

Loss coefficient K
head loss  K  (dynamic head)
V2
hK
2g

Commercial pipe fittings Entry/exit losses

Fitting K Configuration K
Globe valve 10 Bell-mouthed entry 0
Gate valve – wide open 0.2 Abrupt entry 0.5
Gate valve – ½ open 5.6 Protruding entry 1.0
90 elbow 0.9 Bell-mouthed exit 0.2
Side outlet of T-junction 1.8 Abrupt enlargement 0.5

5
 7/9/2013

PIPELINE CALCULATIONS

D h
Q

Main Design Parameters:

Head loss: h
Quantity of flow: Q
Diameter: D

Other Parameters:
Length: L
Roughness: ks
Kinematic viscosity: ν
Minor loss coefficient: K

Method: available head = sum of head losses along the pipe

CALCULATION FORMULAE

1. Head Losses

L V2
h  (λ  K )( )
D 2g

2. Loss coefficients
e.g. friction factor (Colebrook-White):

1  k 2.51 
 2.0 log10  s  
λ  3.7 D Re λ 

HEADS AT THE ENDS OF PIPES

z1
Smooth exit to a downstream reservoir:
H1 = z1
H2 = z2
z2

No residual dynamic head at exit.

Free jet to atmosphere (or abrupt exit to a tank):

H1 = z1
z1
H2 = z2 + V22/2g

Dynamic head must be included at exit.

z2
V2

6
 7/9/2013

TYPICAL PIPELINE CALCULATIONS

Type 1 – flow
Know: diameter D, head h
Find: discharge Q

Easy!

Type 2 – head
Know: diameter D, discharge Q
Find: head h

Solve Colebrook-White equation (iteratively)

Type 3 – size
Know: discharge Q, head h
Find: diameter D

Solve Colebrook-White and head-loss equations simultaneously and iteratively

EXAMPLE SHEET, Q8
Crude oil (specific gravity 0.86, kinematic viscosity 9.010-5 m 2 s–1) is to be
pumped from a barge to a large storage tank. The pipeline is horizontal and of
diameter 250 mm, length 400 m and roughness 0.1 mm. It enters the tank 8 m
below the level of oil in the tank. When the control valve is fully open the static
pressure at pump delivery is 3105 Pa gauge. Ignore minor losses due to pipe
fittings, entrance/exit losses etc.
Control
valve 8m
Pump

Storage tank
Barge

Find:
(a) (using hydrostatics) the gauge pressure where the pipe enters the tank;
(b) (from the pressures at the two ends) the head loss along the pipeline;
(c) the volumetric flow rate in the pipeline.

If the pump delivery pressure remains the same but a valve reduces the flow by
half, find:
(d) the head loss at the valve;
(e) the power loss at the valve.

EXAMPLE SHEET, Q5
(a) A pipeline is to be constructed to bring water from an upland storage
reservoir to a town 30 km away, at an elevation 150 m below the water
level of the reservoir. In summer the pipeline must be able to convey up to
5000 cubic metres per day. If the pipe is fabricated from material of
roughness 0.3 mm, find the required diameter.

(b) During the winter, water requirements fall to only 3000 cubic metres per
day and the excess head available can be used to drive a small turbine. If
the turbine has an efficiency of 75% find the maximum power output.

The Colebrook-White equation is

1  k 2.51 
 2.0 log10  s  
λ  3.7 D Re λ 
where λ = friction factor, ks = roughness, D = pipe diameter, V = average velocity,
Re  VD/ν = Reynolds number.

For water, take density ρ = 1000 kg m –3 and kinematic viscosity ν = 1.010–6 m 2 s–1.

7
 7/9/2013

EXAMPLE, PAGE 15

A reservoir is to be used to supply water to a factory 5 km

away. The water level in the reservoir is 60 m above the
factory. The pipe lining has roughness 0.5 mm. Minor
losses due to valves and pipe fittings can be
accommodated by a loss coefficient K = 80. Calculate the
minimum diameter of pipe required to convey a discharge
of 0.3 m3 s–1.

GRAPHICAL REPRESENTATION OF HEAD

p V2
Energy Grade Line (EGL) z Total head
ρg 2g

p
Hydraulic Grade Line (HGL) z Piezometric head
ρg

GRAPHICAL REPRESENTATION OF HEAD

energy
grade lin
e
hydrauli
c grade
line
reservoir V2/2g
Pipe friction only
pipeline
p/g reservoir

entry loss
Pipe friction with minor EGL
HGL
losses (exaggerated), exit loss

including change in pipe pipeline

diameter.

EGL
HGL
Pumped system
pipeline

pump

8
 7/9/2013

PIPE NETWORKS: EXAMPLE

A C

Which way does the flow go in pipe BD?

PIPE NETWORKS: ELECTRICAL ANALOGUE

90  10 

10 V A C 0V

10  90 

What are the voltages at B and D?

Which way does the current go in BD?

PIPE NETWORKS: BASIC RULES

1. Continuity: at any junction,

Qin   Qout
total flow in = total flow out

2. Each point has a unique head, H

3. Each pipe has a head-loss vs discharge (resistance)

relation:
h = αQ2

9
 7/9/2013

ELECTRICAL ANALOGY

Continuity; unique head ↔ Kirchoff’s Laws

head H ↔ potential V
discharge Q ↔ current I

Resistance law:
head loss H  Q2 ↔ potential difference V  I

What are the hydraulic analogues of:

 a resistor?
 a capacitor?
 an inductor?
 a transistor?

PIPES IN SERIES AND PARALLEL

1 2
Pipes in series

Q = Q1 = Q2 same flow
H = H1 + H2 add head changes
α  α1  α 2
R  R1  R2

1
Pipes in parallel
2

H = H1 = H2 same head change

Q = Q1 + Q2 add flows 1 1 1
 
α α1 α2
1 1 1
 
R R1 R2

JUNCTION PROBLEMS: METHOD

A

B
J ?
C
Method: Adjust HJ until net flow out of J = 0

(0) Establish the head vs discharge relations for all pipes

H  H  αQ2 etc.
J A JA

(1) Guess an initial value of head at the junction, HJ.

(2) Calculate flow rates in all pipes, QJA etc.

(3) Calculate net flow out of junction, QJA  QJB  QJC

(4) Adjust the head at the junction, HJ, until net flow out of junction = 0
QJA  QJB  QJC  0

10
 7/9/2013

EXAMPLE SHEET, Q17

In a water-storage scheme three reservoirs A, B and C are connected by a single
junction J as shown. The water levels in A, B and C are 300 m, 200 m and 140 m
respectively. The pipeline properties are given below. Friction factors may be
assumed constant and minor losses may be neglected.
300 m
A
Pipeline JA JB JC
Length L (m) 5000 3000 4000
Diameter D (m) 0.4 0.25 0.3 200 m
Friction factor λ 0.015 0.03 0.02 B
140 m
C
J

Calculate the total flow in each pipe and the direction of flow in pipe JB if:

(a) there is a valve-regulated flow of 50 L s–1 to reservoir C but water flows freely
under gravity in the other pipes;

(b) water flows freely under gravity in all pipes.

FLOW IN PIPES AND OPEN CHANNELS

PIPE FLOW OPEN-CHANNEL FLOW

Fluid: LIQUIDS or GASES LIQUIDS (free surface)

Driven by: PRESSURE, GRAVITY or BOTH GRAVITY (down slope)

Size: DIAMETER HYDRAULIC RADIUS

Volume: FILLS pipe Depends on DEPTH

Equations: DARCY-WEISBACH (head loss) MANNING’S FORMULA

COLEBROOK-WHITE (friction factor)

NORMAL FLOW

h Q

 Normal flow = steady, uniform flow

(constant-depth flow under gravity)

 At best an approximation for rivers / natural channels

 For any given Q there is a particular normal depth

11
 7/9/2013

NORMAL FLOW
EGL
2
V /2g
hf
HGL (fr
ee surfa
ce): p =
0
h

In normal flow:
 Equal hydrostatic pressure forces at any cross section
 Downslope component of weight balances bed friction
 Channel bed, free surface (= HGL) and EGL are parallel;
i.e. loss of fluid head equals drop in height
 Usual to assume small slopes

PART 1: BALANCE OF FORCES

A = area of fluid cross-section

P = wetted perimeter A
P

L downslope component of weight = friction on sides

mg sin θ  τb  wetted surface area

ρALg sin θ  τb PL
mg

b A
 ρg sin θ  τb
P

cross  sectional area A

Hydraulic radius (*** depends on depth ***): Rh  
wetted perimeter P

Normal-flow relationship: τb  ρgRh S

PART 2: EXPRESSION FOR FRICTION

τb  ρgRh S Rh is the hydraulic radius

c f ( 12 ρV 2 )  ρgRh S definition of the skin-friction coefficient

2g
V2  Rh S
cf

Chézy’s Formula: V  C Rh S

Robert Manning (compilation of experimental data):

C  Rh1/ 6  function of roughness
Rh1/ 6

n

1 2 / 3 1/ 2
Manning’s Formula: V Rh S
n

12
 7/9/2013

MANNING’S ROUGHNESS COEFFICIENT

Channel type Surface n (m–1/3 s)
Glass 0.01
Brass 0.011
Steel, smooth 0.012
painted 0.014
riveted 0.015
Cast iron 0.013
Concrete, finished 0.012
Artificial lined channels
unfinished 0.014
Planed wood 0.012
Clay tile 0.014
Brickwork 0.015
Asphalt 0.016
Corrugated metal 0.022
Rubble masonry 0.025
Clean 0.022
Gravelly 0.025
Excavated earth channels
Weedy 0.03
Stony, cobbles 0.035
Clean and straight 0.03
Natural channels Sluggish, deep pools 0.04
Major rivers 0.035
Pasture, farmland 0.035
Light brush 0.05
Floodplains
Heavy brush 0.075
Trees 0.15

CALCULATION FORMULAE (SUMMARY)

1 2 / 3 1/ 2
Manning’s Formula: V Rh S
n

V = average velocity
n = Manning’s roughness parameter
S = slope (gradient)
A cross- sectionalarea
Rh = hydraulic radius  
P wetted perimeter
Method
For a given channel:
1. Write area A and perimeter P as functions of a parameter (often depth, h)
2. Calculate hydraulic radius
3. Calculate average velocity
4. Calculate quantity of flow

Two Main Types of Problem

 Given h find Q
 Given Q find h

EXAMPLE SHEET, Q20

A V-shaped channel with sides sloping at 30º to the horizontal

has a gradient of 1 m in 100 m and an estimated Manning’s n
of 0.012 m–1/3 s.

Calculate:

(a) the discharge for a depth of 0.5 m;

(b) the depth when the discharge is 2 m3 s–1.

13
 7/9/2013

EXAMPLE SHEET, Q21

A concrete pipe 750 mm in diameter is laid to a gradient of 1

in 200. The estimated value of Manning’s n is 0.012 m–1/3 s.

Calculate the discharge when:

(a) the pipe is full;

(b) the depth is 90% of maximum.

Explain why the answer in (b) exceeds that in (a).

EXAMPLE SHEET, Q25

A culvert used to divert run-off has a rectangular cross

section with base width 0.4 m and side heights of 0.3 m.
Manning’s coefficient may be taken as n = 0.012 m–1/3 s.

(a) Find the minimum slope S necessary to carry a

discharge Q = 0.3 m3 s–1.

(b) If the slope from part (a) is doubled for the same
discharge, calculate depth of flow.

CONVEYANCE
1 2 / 3 1/ 2 A
(a) Manning’s formula: V  Rh S Rh 
n P
(b) Discharge: Q  VA
2/3
1 A
Q   S1/ 2  A
n P

Q  KS 1/ 2 1 A5 / 3
K
n P2/3
conveyance

For compound channels (e.g. river plus flood plain) simply add the conveyances:
1 2 3

flood plain river flood plain

Keff  K1  K 2  K3

14
 7/9/2013

COMMON SHAPES OF CHANNEL

rectangle trapezoid circle

h h
  R
b b h

h2
area A bh bh 
tan α R 2 (θ  12 sin 2θ)

wetted perimeter P 2h 2Rθ

b  2h b
sin α

15