CE 540 Unsteady Flows in Rivers and Pipe Networks

Analysis of Pipe Networks

Applications: Water supply system, Pressure sewerage system, Water treatment plant
Theory:

A1V1

1. Continuity equation
2. Energy equation
V12
2g

A2V2

Energy line

hp
Hydraulic grade line

p1

V22
2g

p2

V1

V2

P
z1

datum

z2

V12
z1 

 hp
J 2g
p1

hL

V22
z2 

 hL
J 2g
L V2
f
D 2g
p2

Friction factor f is a function of Reynolds number (RN) and relative roughness

( k s / D ). It can be found from Moodys diagram or Colebrook equation:

1
f

k /D
2.51
2 log s

3.7
RN f

k / D 5.74
2 log s
 0.9
RN
3.7
Table 1. Equivalent sand roughness, ks in new pipes, (Moody, 1944)
Material
Equivalent sand roughness, ks (mm)
Riveted steel
0.9-9.0
Concrete
0.3-3.0
Ductile and cast iron
0.26
Galvanized iron
0.15
Asphalt-dipped ductile/cast iron
0.12
Commercial steel or wrought iron
0.046
Copper or brass tubing
0.0015
Glass, plastic (PVC)
|0

Jain's approximation :

1
f

Example: Water from a treatment plant is pumped into a distribution system at a rate
of 4.38 m3/s, a pressure of 480kPa, and a temperature of 20oC. The diameter of the
pipe is 750mm and is made of ductile iron.
(a) Calculate the friction factor by Colebrook equation
(b) Estimate the pressure 200m downstream of the treatment plant if the pipeline
remains horizontal.
Solution:

0.016

p2

270 kPa

Designing a pipe for a given flowrate and head loss

1. Assume f
2. Calculate
diameter
from
rearranged
Darcy-Weisbach

8 LQ 2

f
2
h L gS

equation:

1/ 5

3. Calculate Reynolds number, R N

VD

4. Calculate relative roughness, k s / D

5. Use RN and ks/D to calculate f from Colebrook equation.
6. Using new value of f repeat the procedure until new f agrees with the old one.
Example: A galvanized iron service pipe from a water main is required to deliver
200 L/s during a fire. If the length of the service pipe is 35m and the head loss in the
pipe is not to exceed 50m, calculate the minimum pipe diameter that can be used.
Solution: D= 136 mm
Three-reservoir problem
A
B
V1 , L1 , D1
V 2 , L2 , D 2

zA

zB
V3 , L3 , D3
datum

C
zC

Unknowns: V1 , V 2 , V3
Equations: One continuity and two energy equations

Other friction formulas:

Hazen-Williams formula (SI units)

0.85C H R 0.63 S 0.54

Mannings formula (SI units)

1 2 / 3 1/ 2
R S
n

Table 2. Pipe roughness coefficients

Pipe material
CH
Ductile and cast iron:
New, unlined
130
Old, unlined
80
Cement lined and seal coated
120

0.013
0.013

Steel:
Welded and seamless
Riveted
Concrete
Vitrified clay
Polyvinyl chloride (PVC)

0.015
0.015
0.013
0.009

120
110
110
-

The head loss can be expressed as:

hL
Mannings formula:

Darcy-Weisbach equation:

KQ x

Hazen-Williams formula:

x 1.75

(smooth pipe) to

x 1.85
x

(rough pipe)

Design of water distribution system

A municipal water distribution system is used to deliver water to the consumer.
Water is withdrawn from along the pipes in a pipe network system, for computational
purposes all demands on the system are assumed to occur at the junction nodes.
Pressure is the main concern in a water distribution system. At no time should the
water pressure in the system be so low that contaminated groundwater could enter the
system at points of leakage.
The total water demand at each node is estimated from residential, industrial and
commercial water demands at that node. The fire flow is added to account for
emergency water demand.

Water Demand
Total water demand at a node is equal to the average water demand per capita
multiplied by the population served by that node.
The population may be estimated by any of the available models of population
growth. These models use the previously available population data for future
projections.
For example, according to Arithmetic model:

P(t )

P0  kt

P(t) is the population at time t and P0 is the reference population.

Similarly a higher order polynomial function can be obtained as follows:

P(t )

P0  at  bt 2  ct 3  

and the population can be estimated at any time t.

Example: You are in the process of designing a water-supply system for a town, and
the design life of your system is to end in the year 2020. The population in the town
has been measured every 10 years since 1920 and is given below. Estimate the
population in the town using graphical extension and arithmetic growth projection.
Year
1920
1930
1940
1950
1960
1970
1980
1990

Population
125,000
150,000
150,000
185,000
185,000
210,000
280,000
320,000

Solution:
500000
400000
300000
200000
100000
0
1920

1940

1960

1980

2000

2020

From the graph, if the line is extended to year 2020, we get, the population as:

440,000

60000  5500t

By regression on the data for 1970 to 1990, we get the following expression:

P (t )

where, t is the time in years starting from year 1920. So, for year 2020, t= 100, and
we get the population,

490,000

Variations in demand
Maximum daily demand, Maximum hourly demand
Table 3. Typical demand factors
Condition
Range of demand factors
Daily average in maximum month
1.1-1.5
Daily average in maximum week
1.2-1.6
Maximum daily demand
1.5-3.0
Maximum hourly demand
2.0-4.0
Minimum hourly demand
0.2-0.6

Typical value
1.2
1.4
1.8
3.25
0.3

Fire demand

C i Oi ( X  P ) i

Insurance Services Office, Inc. (ISO, 1980) formula:

NFFi

where, NFFi is needed fire flow at location i, Ci is the construction factor based on
size and type of construction of the building, Oi is the occupancy factor reflecting the

kinds of material stored in the building (value range from 0.75 to 1.25), and ( X  P) i
is the sum of the exposure factor and communication factor that reflects the
proximity and exposure of other buildings (value range from 1.0 to 1.75).

C i ( L / min)

220 F Ai

Ai is the effective floor area in square meters, typically equal to the area of the largest
floor in the building plus 50% of the area of all other floors, F is a coefficient based
on the class of construction. The maximum value of Ci and typical F values are
given below:
Class of Construction
Description
Max. Ci(L/min)
F
1
Frame
1.5
30,000
2
Joisted masonry
1.0
30,000
3
Noncombustible
0.8
23,000
4
Masonry, noncombustible
0.8
23,000
5
Modified fire resistive
0.6
23,000
5

6
Fire resistive
0.6
Occupancy factors:
Combustibility class
Examples
C-1 Noncombustible
Steel or concrete products storage
C-2 Limited combustible
Apartments, mosques, offices
C-3 Combustible
Department stores, supermarkets
C-4 Free-burning
Auditoriums, warehouses
C-5 Rapid burning
Paint shops, upholstering shops

23,000
Oi
0.75
0.85
1.0
1.15
1.25

Average value of ( X  P) i is 1.4. The NFF should be rounded to the nearest

1000L/min if less than 9000 L/min and to the nearest 2000L/min if greater than
9000L/min.
Required fire flow durations:
Required fire flow (L/min)
Duration (h)
<9000
2
11000-13000
3
15000-17000
4
19000-21000
5
23000-26000
6
26000-30000
7
30000-34000
8
34000-38000
9
38000-45000
10
Example: Estimate the flowrate and volume of water required to provide adequate
fire protection to a 10-story noncombustible building with an effective floor area of
8000m2.
Solution: NFF=17000L/min
Volume= 4080m3.
Design periods and capacities in water supply systems:
Component
Design period(year)
Design capacity
Sources of supply:
River
Indefinite
Maximum daily demand
Wellfield
10-25
Maximum daily demand
Reservoir
25-50
Average annual demand
Pumps*:
Low-lift
10
Maximum daily demand
High-lift
10
Maximum hourly demand
Water treatment plant
10-15
Maximum daily demand
Service reservoir
20-25
Working storage + fire
demand + emergency
storage
Pipes
25-50
Greater of maximum daily
demand + fire flow and
maximum hourly demand
6

*: One reserve unit should be kept for the pumps.

Example: A metropolitan area has a population of 130000 people with an average
daily demand of 600 L/d/person. If the needed fire flow is 20000 L/min, estimate, (a)
the design capacities for the wellfield and the water-treatment plant, (b) the duration
that the fire flow must be sustained and the volume of water that must be kept in the
service reservoir in case of a fire and (c) the design capacity of the main supply
pipeline to the distribution system.
Solution: (a) 1.62 m3/s (b) 5940 m3 (c) 2.92 m3/s
Operating criteria for water supply systems
Primary functions of a water-supply system
1. Meet the demand while maintaining acceptable pressures
2. Supply water for fire without affecting the water supply to the rest of the
system
3. Provide sufficient level of redundancy to serve during emergency conditions
Minimum acceptable pressure:
Under normal conditions = 240 to 410 kPa
During fire or emergency = >140 kPa

Maximum pressure (not strict) : 650kPa.

Storage facilities
1. 20% to 25% of the maximum daily demand volume
2. Fire demand
3. Emergency storage (minimum storage equal to average daily system denmand)
4. Minimum height of water in the elevated storage tank based on minimum
piezometric head in the service area
5. Normal operating range 4.5 to 6.0 m.
Access
manhole

Roof vent

Normal operating range

Fire/emergency range
Overflow
pipe

Example: A service reservoir is to be designed for a water-supply system serving

250,000 people with an average demand of 600L/d/capita and a needed fire flow of
37000 L/min. Estimate the required volume of service storage.
Solution: 237600m3.
Example: A water-supply system is to be designed in an area where the minimum
allowable pressure in the distribution system is 300kPa. A hydraulic analysis of the
distribution network under average daily demand conditions indicates that the head
loss between the low-pressure service location, which has a pipeline elevation of
5.4m and the location of the elevated storage tank is 10m. Under maximum hourly
demand conditions, the head loss between the low-pressure service location and the
elevated storage tank is 12m. Determine the normal operating range for the water
stored in the elevated tank.
Solution: From 48m to 46 m.
Pipe Network Systems
T-2
T-1
P1

P4

C2
J-1 P-2

J-2 P-3

J-3
C3

P5

C1

J-4

P7

P6
P-8

C4

N eq

J-5

P-9

C5

J-6
C6

N j  N p  N f 1

Neq= Number of equations required

Nj= number of junction nodes
Nl= number of loops
Nf= number of Fixed Grade Nodes (e.g. elevated reservoirs)
In the above figure

N eq

6  2  2 1 9
8

Conditions to satisfy:
1. The algebraic sum of the pressure drops around any closed loop must be zero
2. The flow entering a junction must equal the flow leaving it.

Qa  '

Hardy Cross Method

Q
Q
Qa
'

Actual flow rate

Assumed flow rate
Correction

K (Q

 ') x

KQax  xK'Qax 1 
a

x 1
xK '2Qax  2   0
2

x
x 1
KQ
KxQ

'
a a

For small values of correction

'

KQ


xKQ

x
a

x 1
a

Example: Determine the discharge in each line and the pressure at each junction node
for the pipe network shown below using Hardy Cross method. The pipe and junction
data is given in the table below.

A
B
C2

P1

J-2

P-2
J-1

P-5

J-4

PP-3 4

C1

P7

P-6
J-3

C4

C3
Line
1
2
3
4
5
6
7

Nodes
A-1
1-2
1-3
2-3
2-4
3-4
4-B

Length(m)
300
250
200
220
180
250
270
Node
A
B
1
2
3
4

Diameter(cm)
30
20
20
15
20
20
25
Elevation(m)
126
123
96
99
93
90

From Darcy-Weisbach formula,

hL

f
0.015
0.019
0.019
0.02
0.019
0.019
0.017

K(s2/m5)
153.0
1226.5
981.2
4787.6
883.1
1226.5
388.4

Demand(L/s)
56
112
28
85

KQ 2 , where
8 fL
S 2 gD 5

HGL elevation at the end of pipe = HGL elevation at the beginning-Head loss

p J ( HGL El.  G.El.)

10

A
B
P1

J-2
P-2
J-1

P-5

P7
J-4

PP-3 4

P-6
J-3

Pipe
1
2
3
4
5
6
7
Node
A
1
2
3
4
B

Discharge(L/s)
167.9
58.4
53.5
16.4
37.2
9.1
113.1

Elevation(m)
126
96
99
93
90
123

Head loss (m)

5.61
4.78
3.22
1.55
1.42
0.13
5.97

HGL (m)

p(kPa)

120.39
115.61
117.17
117.03

239.3
163.0
237.1
265.2

11

Linear Method
Hardy Cross method requires an initial estimate of the flow that should be close to the
final solution.
Linear method is a very stable method in which all the equations are solved
simultaneously.
Each nonlinear head loss term is linearized using first two terms in Taylor series:

f (Q )

f (q) 

wf
(Q  q )
wq

q= flow rate from previous iteration

Q= unknown flow rate.

Kq n  nKq n1 (Q  q)

So,

KQ n
where

DQ  D'

D' (1  n) Kq n

nKq n1

For loops with pumps, nonlinear pump characteristic equation (LHS) is replaced with
linear equation (RHS):

AQ 2  BQ  H

Aq 2  Bq  H  (2 Aq  B)(Q  q)

where

2 Aq  B

E ' H  Aq 2

EQ  E '

Example: Write down the equations to solve the following network by linear method.

A
B
P1

C2
P-2
J-1

C1

J-2
P-5

P7
J-4

PP-3 4

P-6
J-3

C4

C3
Solution: First of all we assume the direction of flow in each pipe.

12

A
B
P-1

C2

P-7

J-2

P-2

P-5
J-4

P-4

J-1

P-6

P-3
C1

C4
J-3
C3

Continuity equations:

Q1  Q2  Q3  C1

Q2  Q4  Q5  C2

 Q5  Q6  Q7  C4

Q3  Q4  Q6  C3

0
0

Loop equations (clockwise positive):

K 2Q22  K 3Q32  K 4Q42

K 4Q42  K 5Q52  K 6Q62

0
0

 K1Q12  K 2Q22  K 5Q52  K 7Q72

D2Q2  D3Q3  D4Q4

zB  z A

 D2c  D3c  D4c

 D4c  D5c  D6c

The continuity equations are already linear, so linearizing the loop equations we get:

D4Q4  D5Q5  D6Q6

 D1Q1  D2Q2  D5Q5  D7Q7

C5
C6

z B  z A  D1c  D2c  D5c  D7c

C7

>A@>Q@ >C @
Where > A@ is coefficient matrix and >C @ is a column matrix of constants.

In matrix form the above equations can be written as follows:

13

 1
0

0
0

0
 D
1

1

0
0

1
0

1

1
0

D2

 D3

 D4

0
 D2

0
0

D4
0

Solving for the flow rates:

0
1
0

 D5
D5

0 Q1
0
0 Q2
1
0 Q3

 1 1 Q4
0
0 Q5

D6 0 Q6
0 D7 Q7
0

C1
C
2
C 3

C 4
C 5

C 6
C
7

>Q@ >A@1>C @

For the initial iteration a velocity of 1m/s may be assumed in each pipe. The
computed flow rates (Q) become estimated flow rates (q) for the next iteration. The
iterative process is continued until the equations converge to a solution. In other
words when the value of

Q q
N

Q
i

i 1

i 1

is within a specified limit.

14