Analysis of Flow in Water Distribution Networks PDF
Analysis of Flow in Water Distribution Networks PDF
FLOW in WATER
DISTRIBUTION
NETWORKS
TECHNOMIC
^PUBLISHING CO., INC.,
T ANCASTER • BASET ,
Analysis of Flow in Water Distribution Networks
aTECHNOMICpublication
Preface xv
t INTRODUCTION
1.1 Introduction 1
1.2 Historical Development 2
1.3 Fluid Properties 3
1.3.1 Density, Specific Weight, Specific Gravity
1.3.2 Compressibility 5
1.3.3 Viscosity 6
1.3.4 Other Properties 8
1.3.5 Properties of Water 8
1.4 Continuity Principle 8
1.5 Energy Principle 14
1.5.1 Forms of Energy 15
1.5.2 Bernoulli's Theorem 16
1.5.3 Hydraulic and Energy Grade Lines 20
1.6 Momentum Principle 20
4 EQUIVALENT PIPES 69
4.1 Introduction 69
4.2 Pipes in Series 69
4.2.1 Discharge Known 70
4.2.2 Total Head Loss Known 73
4.3 Pipes in Parallel 75
4.3.1 Head Loss Known 11
4.3.2 Total Discharge Known 78
Table of Contents ix
CALIBRATION 385
14.1 Introduction 385
14.1.1 Definition 385
14.1.2 Reasonable Agreement 386
14.2 Data Collection and Preparation 387
14.2.1 Skeletonization of Network 387
14.2.2 Heads at Source Nodes 389
14.2.3 Supply Rates at Source Nodes 389
14.2.4 Consumption at Demand Nodes 389
14.2.5 Pipe Details 389
14.3 Calibration Methods 391
14.3.1 Methods Adjusting Pipe Coefficients 392
14.3.2 Methods Simultaneously Adjusting Pipe Coefficients and
Nodal Demands 395
14.4 Practical Considerations 405
14.4.1 Fine-Tuning the Model 405
14.4.2 Sources of Error in Field Data Collection 406
14.4.3 Field Equipment 408
Nomenclature 425
Answers 435
References 441
Index 453
T is taken for granted that water will always be available in the desired
I quantity when we open faucets, without our bothering to know how it is
brought to us. It is only when we do not get a sufficient quantity of water
that we start grumbling about the water supply engineers. This is the
natural human tendency and it is observed in all walks of life. It is there-
fore necessary for water supply engineers to have a good knowledge of
water supply engineering in general and water distribution systems in par-
ticular. This will help in the proper planning, analysis and design of water
distribution systems. The availability of water at different localities under
different demand patterns and at different water levels in the reservoirs, the
dependability of the system under fire flow requirements or pipe breaks,
and the expansion of the existing water distribution systems to accom-
modate new residential localities or industries are some examples where
the proper understanding of network behaviour would be necessary. It was
with this intention that a two-semester course on water distribution sys-
tems was introduced at Nagpur University—an elective subject in the final
year of the undergraduate civil engineering curriculum, and a compulsory
one for the post graduate courses in hydraulic engineering and environ-
mental engineering disciplines. The present text is the outgrowth of the
lecture notes prepared for them over the past several years.
The emphasis throughout the text has been on explaining the fundamen-
tal concepts. Numerous examples have been provided in the text. The ex-
amples are simple enough so that they can be tackled by calculators but are
complex enough to explain the finer points. Problems are provided, wher-
ever necessary, for the readers to test their knowledge. Answers are listed
towards the end of the book (pp. 435-440). No reference to computers or
computer programs is made in the text (except in the end) because the
reader, having understood the fundamental concepts, can easily formulate
programs in BASIC or FORTRAN languages. It is however emphasized
that the distribution networks encountered in practice are not simple
xv
xvi Preface
PRAMOD R. BHAVE
n
Introduction
1.1 INTRODUCTION
Density
The mass per unit volume of a substance is termed its density or specific
mass, and is generally denoted by Q (rho). Its dimensions are ML"3 and
unit is g/cm3 or kg/m3 (1 g/cm3 = 1,000 kg/m3).
Specific Weight
The weight per unit volume of a substance is termed its specific weight,
usually denoted by y (gamma). As it is a force per unit volume, its dimen-
4 INTRODUCTION
1 = Qg (1-1)
Specific Gravity
Example 1.1
Solution
19 040
.'. Mass of the liquid = ' = 1.9409 kg
9ol
= 93 °
Fluid Properties
19 040
Specific weight = 2Q87 ' X 1Q_3 = 9,123 N/m 3
19 040
Specific gravity = 2Q'473 = 0-93
1.3.2 COMPRESSIBILITY
Example 1.2
Find the specific weight of water at a depth of 100 m from the water sur-
face. Take K = 2.07 x 109 Pa.
Solution
V = Vo + dV
-*('•?)
6 INTRODUCTION
981
• v= vi\ -*P\ =-*—(} '°°° \
°\ Kj 9,810 \ 2.07 x 109j
=
m X °" 9 526 ^
.'. Specific weight of water at a 100-m depth
9,810
= 9,814.7 N/m 3 .
0.999 526
(Even at a 100-m depth from the water surface, the increase in the
specific weight of water is hardly 0.05 %. Therefore, considering water as
incompressible is valid in practice.)
1.3.3 VISCOSITY
dv
T = „ - (1.3)
v = 0.5y - If
in which v = velocity in m/s at a point y meters away from the inner sur-
face of the pipe wall. Find the shear stress at the pipe wall and at 0.05 m
from the pipe wall. Also determine the total resistance in a 100-m length
of the pipe.
Solution
dv
From Newton's law of viscosity, Equation (1.3), T = /*—
As v = 0.5y - f, ^ = 0.5 - 2y
•• ( ? ) = 0.5 s-
\dy/(^o)
At pipe wall, T0 = 0.04 X 0.5 = 0.02 Pa.
Since water is the primary fluid that is considered in this text, some
properties of water at different temperatures are given in Table 1.1.
Consider a control volume bounded by a stream tube and two cross sec-
tions 1-1 and 2-2 as shown in Figure 1.1. Let the fluid enter through section
1-1 and leave through section 2-2. If the flow is steady, the rate of accu-
mulation of the fluid mass between sections 1-1 and 2-2 is zero. Thus, the
fluid mass rates passing through sections 1-1 and 2-2 are equal. Summing
up the mass flow rate over a collection of stream tubes, the continuity
equation for mass flow rate for steady flow in a pipe can be written as
M= QQ = Q I vdA =
= Q
Q \\ vdA (1.6)
«Mi J
Q = I vdA
vdA == I\ vdA (1.7)
J At J/A2
If F, and V2 are the average velocities at sections 1-1 and 2-2 with cross-
sectional areas At and A 2 , respectively, the continuity equations for mass
10 INTRODUCTION
flow rate and volumetric flow rate in incompressible flow become, respec-
tively
EM = LQQ =0 (1.10)
LQ = LAV = 0 (1.11)
Example 1.4
.Q2
T - * • 2 m/s
pipe.1
= 0-020
Continuity Principle 11
flow in pipes 1 and 3 are 2 m/s towards J and 3 m/s away from J, whereas
the discharge in pipe 4 is 0.020 m3/s away from J. Determine the velocity
and discharge in pipe 2.
Solution
For Pipe 1,
Qx = AM = J ( 0 . 3 ) 2 X 2 = 0.1414 nrVs
Assuming that the flow in pipe 2, i.e., Q2, is towards the junction and
taking the flow towards the junction as positive and that away from the
junction as negative, the continuity equation for the volumetric flow rate
gives
Qi + ft - Ga - QA = 0
or Q2 = - 0 . 0 2 7 2 m3/s
Solution
and the velocity of flow through it is v = 1.2 (y/r)u7. Therefore, the dis-
charge 5Q passing through it is given by
2 4x
6Q = 2TT (r - y) 6y X 1.2 (ylr)"7 = —jr(ryU7 - / " ) 8y
Q
" Jo '" 50
49
Q = — ,r(0.1)2 = 0.03079 nvVs
and
O 49
the average velocity V = — = — = 0.98 m/s
y4 50
Example 1.6
Reservoir,1
Reservoir, 2
60 TT
9
20 30 20 20
r
20 30 10
and the supplies at the junctions by inward arrows. All the discharges are
in liters per minute, L/min. Determine the supply rates from the two reser-
voirs and the discharges in the pipes.
Solution
Example 1.7
Determine the supply rate at junction 6 and the discharges in pipes of the
distribution network shown in Figure: Example 1.7.
0-1 -
14 INTRODUCTION
Solution
By adding the continuity principle to the entire network, the supply rate
at node 6 = 0.2 + 0.3 + 0.1 + 0.3 - 0.5 = 0.4.
Assuming the flow directions in the pipes as shown in Figure: Exam-
ple 1.7, we get, 2,-4 = 0.5 - G,-2; Q4-s = Q^ - 0.1 = 0.4 - 0,- 2 ;
ft-s = Qi-2 - 0.2 - 0 2 . 3 ; 06-5 = 0.3 - 0 4 - 5 - Q2-s = 0.3 -
(0.4 - e,_2) - (e,_ 2 - 0.2 - Q2-3) = 0.1 + Q2.3; Q6.3 = 0.3 -
G2-3.
At junction 6, the supply rate = Q6-s + Q6-3 = (0.1 + Q2-3) +
(0.3 — Q2-3) = 0.4, as obtained from the consideration of the entire net-
work. In this example:
(1) The supply rate at junction 6 can be found either by considering the
network as a whole or by determining the discharge in pipes 6-5 and
6-3 and then considering junction 6.
(2) The absolute values for the pipe discharges cannot be determined by
simply applying the continuity principle. They can be expressed
only in terms of discharges in some pipes. In the solution herein, the
pipe discharges are expressed in terms of the discharge in pipes 1-2
or 2-3 or both.
Potential Energy
Flow Energy
Reservoir
PA-
d L
Figure 1.2 Flow energy.
16 INTRODUCTION
Kinetic Energy
The energy possessed by a fluid mass due to its velocity is termed kinetic
energy and is equal to mass times V2/2. Thus, the kinetic energy possessed
by a unit mass is V2/2.
It can be easily seen that each of the above forms of the available energy
per unit mass has the dimensions of FL/M or L2T~2 and units of N • m/kg
or J/kg. These forms of energy are interchangeable and their sum consti-
tutes the total energy E per unit mass of a fluid. Thus,
p V2
£ per unit mass = gZ + — + -y- (1-12)
Instead of considering the energy per unit mass, many times the energy
is considered per unit weight. As weight is g times mass, by dividing each
term of Equation (1.15) by g, Equation (1.15) can be written as
7 = Qg; and hL - —
Each term in Equation (1.16) represents energy per unit weight and
therefore has the dimensions of FL/F, i.e., L and unit of N-m/N, i.e.,
m. Each term represents a height of liquid column and is therefore called,
in practice, a head. Thus, Z represents the potential head, position head,
or elevation head; ply represents the pressure head; and V2/2g the velocity
head. As the liquid would rise to a height of Z 4- {ply) in a piezometer
when it is introduced at right angles to the direction of flow in a flowing
liquid, the sum Z + {ply) is termedpiezometric head. It is also called the
hydraulic head. The sum Z + {ply) + {V2/2g) is the total head.
Between sections 1-1 and 2-2 of a pipeline, a pump may be fitted to sup-
ply energy or a turbine may be fitted to extract energy. Representing such
externally supplied or extracted energy per unit mass by Ee{Ee = Ep for
pumps and — ET for turbines) and head by he{he = hp for pumps and — hT
for turbines), the generalized forms of the Bernoulli's equation for unit
mass and unit weight are, respectively,
Example 1.8
Solution
(0 2) X
Q \ ' ^
Velocity at section 2, V2 = ^ = = 7.2 m/s
80 X 1,000 1.82
= + +
° 9,810 Y^9Jl =8J2° m
° f Water
7 2g
/ 24.3 X 1,000 7 22
= 6 + I - " ' " ~ 1' I + ~ — 1 r r ^ 7 = 6.165 m of water
\ ~ 9,810
Energy Principle 19
Since the total head at section 1 is greater than that at section 2, the
direction of flow is from section 1 to section 2.
Loss of head = /?, - h2 = 8.320 - 6.165 = 2.155 m
Example 1.9
A pump lifts 0.15 m3 of water per second. The diameter of the pump in-
let is 300 mm and the vacuum pressure at the inlet is 150 mm of mercury.
The outlet diameter of the pump is 200 mm and the pressure at the outlet
is 137.4 kN/m2. If the exit of the pump is 1 m higher than the inlet and the
power of the motor driving the pump is 35 kW, find the efficiency of the
pump.
Solution
:
Similarly, at outlet, V2 = = 4.775 m/s
I X 0.2-
The pressure head at the inlet is 150 mm of mercury, vacuum.
13 6
.'. — = - T^L X - = -2.04 m of water
7 1,000
Applying Bernoulli's equation to the inlet (section 1) and the outlet (sec-
tion 2) and taking the datum through the inlet,
y 2g 7 2g
2
2.122 137.4 4.7752
•"• ° + ( - 2 0 4 ) +
19^ +
*' = ' + W +
19^
.-. hp = 17.979 m
™ ,. . . . 0.15 X 1000 X 9.81 X 17.979
The energy supplied by the pump = Trwi
= 26.456 kW
The momentum principle for a fluid mass can be derived from Newton's
second law of motion which states that the net force acting on a body in
a fixed direction is equal to the rate of increase of the momentum of the
body in that direction.
Consider the control volume bounded by a stream tube and two cross
sections 1-1 and 2-2 as shown in Figure 1.3. Let the fluid enter through sec-
tion 1-1 and leave through section 2-2, with average velocities K, and V2,
respectively. For steady flow, the mass of fluid entering at section 1-1 is
equal to that leaving at section 2-2, and this, for unit time is QQ. Thus,
from Newton's second law of motion, the net force in the ^-direction, i .e.,
Fx is given by
Example 1.10
respectively, and the diameter at the exit end (section 2-2) is 150 mm. De-
termine the force required to keep the bend stable. Neglect the loss of head
in the bend.
Solution
0.35
at the entrance, V, = = 4.951 m/s and
0.35
at the exit, V2 = = 19.806 m/s
7T
(0.15) 2
P2 = ^Q(V\ - VI)
= 1.101 x 105 Pa
Momentum Principle
IK \
:. 294 x 103\-r x 0.3 2 \ - Fx - 1.101 x 105
•K
x j x 0.15 | cos 60° = 1000
In the j-direction
2.1 INTRODUCTION
HEN a real fluid flows through a pipe, part of the total energy of
W the fluid is spent in maintaining the flow. This used-up energy is con-
verted to thermal energy due to internal friction and turbulence. Such
a conversion, which is the loss of energy as far as its utility is concerned,
is usually expressed in the form of head of liquids and therefore termed
head loss.
The head loss in a pipe is classified into two categories:
The frictional head loss in a pipe is due to the viscosity of the fluid and
the turbulence of the flow and is present throughout the length of the pipe.
As the frictional head loss in a long pipe is relatively larger than the other
head losses, the frictional head loss is also termed major head loss.
When there is a sudden or gradual change in the boundaries of the fluid
or when there is a local obstruction to flow, the flow pattern changes. This
results in a change in the magnitude, direction or distribution of the veloc-
ity of flow. Such a change introduces additional head loss which is of a
local nature and is usually much less than the frictional head loss in a long
pipe. It is therefore termed minor head loss.
It must, however, be remembered that the classification into major and
minor head losses is rather relative, and therefore for a pipeline of small
length and having many minor appurtenances, the total minor head loss
can be more than the frictional head loss. In water distribution networks,
however, the pipelines are of considerable length and therefore the terms
major head loss and minor head loss can be used.
Frictional head loss is considered in this chapter and the head loss due
to other causes, i.e., minor head loss, is considered in the next chapter.
25
26 FRICTIONAL HEAD LOSS IN PIPES
The term VDQ//JL is the Reynolds number Re. The dimensionless param-
eter elD is known as the relative roughness.
Introducing a coefficient of friction / such that
j , M (2.3)
Equation (2.2) can be written as
J. Weisbach [10] was the first to write a resistance equation for pipes in
the form of Equation (2.4). As H. P. G. Darcy carried out considerable
work on pipe flow, his name is also associated with Equation (2.4). There-
fore, Equation (2.4) is now commonly known as the Darcy-Weisbach for-
mula.
For noncircular pressure conduits, D in Equation (2.4) is replaced by AR
in which R is the hydraulic radius defined as the ratio of the cross-sectional
area of flow A, and the wetted perimeter P, i.e., R = A/P.
Nikuradse's Experiments on Artificially Roughened Pipes 27
16 fLQ2 SfLQ2
h( = (2.5)
fLQ2
hf = (2.6)
12.1 D 5
o-io
0O8
0-06
0-04
0-02
- Laminar Transitional
—flow »i flow
0-01
of spacing and shape aspects \e' ID and m in Equation (2.1)] were not stud-
ied in Nikuradse's experiments, the experiments did prove the validity of
the relative roughness concept.
The following conclusions can be drawn from Nikuradse's experiments
(Figure 2.1):
(1) When Re < 2,000, i.e., when the flow is laminar,/is independent
of elD and depends only on Re.
(2) When 2,000 < Re < 4,000, i.e., when the flow is in a state of tran-
sition from laminar to turbulent flow, / increases as Re increases.
(3) When Re > 4,000, turbulence sets in and the effect of elD may be
felt. However, as long as the pipe wall roughness is within the lami-
nar sublayer, the pipe behaves as a smooth pipe. The flow is termed
smooth turbulent flow. The value of/depends on Re only and is
given by the smooth pipe curve.
(4) Depending upon the value of elD, the curve branches out from the
smooth pipe curve. The larger the value of elD, the earlier the
branching takes place. As the value of Re increases, the effect of elD
also increases. This continues until the boundary roughness has its
complete effect on the flow. This flow is termed transitional turbulent
flow. In transitional turbulent flow, the value of/depends on both Re
and elD.
(5) For very large values of Re, the effect of viscosity is negligible, and
/becomes independent of Re and remains constant for a particular
value of elD. This flow is termed rough turbulent flow and /depends
only on the elD value.
Stanton, in collaboration with Pannell, conducted experiments on sev-
eral pipes of different diameters and materials using several fluids. There-
fore, the diagram as in Figure 2.1 is sometimes known as the Stanton dia-
gram or the Stanton-Pannell diagram.
roughness is the roughness which provides the value of/that is the same
as the limiting value for the rough turbulent flow in Nikuradse's artificially
roughened pipes. The values of the equivalent surface roughness e for
pipes of some common materials are given in Table 2.1 [9,13-15].
Johnson [16] provided a diagram for commercial pipes using several
nondimensional groups. Rouse [17] presented a graph of \l\ff versus Kar-
man number K = Re \ff = D3n \l2gSlv upon which the Colebrook
transition curves [12] for various values of the relative roughness were
superimposed. Moody suggested that Rouse should convert the diagram to
the more customary / versus Re form. In the belief that the Prandtl-
Karman parameters were the more significant ones and therefore the form
used by him was superior to that suggested by Moody, Rouse refused to
convert the diagram [18]. Whereupon Moody [19] published his proposed
form in the widely read, "Transactions of the American Society of
Mechanical Engineers." This diagram, shown in Figure 2.2 is now univer-
sally known as the Moody diagram and gives the values of/ for different
values of Re and elD.
Several other investigators such as Powell [20], Ackers [21], Thiruveng-
dam [22], Ramalho, Tiller and Berry [23], Bewtra [24], Barr and Smith
[25], Asthana [26], Li [27) and Barr [28] also presented charts to tackle
different flow problems. However, the Moody diagram is more popular
and widely used in practice.
It is worthwhile to remember that the concept of equivalent surface
roughness on which the Moody diagram is based has some limitations. For
example, the concept of equivalent surface roughness implies that only the
height of surface irregularities significantly affects the flow. However, the
Equivalent Surface
Pipe Material Roughness e, mm
<u|Q
0-000 01
.8
"6 8103 2 4 6 8B 4 2 3 4 6 4 6 4 6 8»7 4 6 810
Reynolds number Re = ^ j
This value of/is the same as that obtained theoretically from the Hagen-
Poiseuille equation for laminar flow.
For transitional flow (2,000 < Re < 4,000) the value of/is uncertain.
However, from experiments it is found to lie between 0.03 and 0.08 for
commercial pipes.
f- ^ ^ (2
(£)
Equation (2.9) is known as Karman-Prandtl equation for smooth turbu-
lent flow. This equation is implicit and the value of/must be obtained by
trial and error.
Several investigators such as Schimemi, Barnes, Blair, Lamont [30] and
Ger and Holley [31] have suggested relationships similar to that of Blasius
for smooth turbulent flow.
For the transitional turbulent flow,/is a function of Re and elD. For this
flow Colebrook [12] suggested a relationship, known commonly as the
Colebrook-White equation:
(2IOa)
or
(2.11a)
Explicit Relationships for Friction Coefficient 33
or
which has an accuracy of ± 5 percent for 4,000 < Re < 107, and
e/D < 0.01. Nahavandi and Catanzaro [35] used this relationship in pipe
network analysis.
Wood [36] suggested a relationship
f = a + b Rec (2.15)
(2.16)
(2.18)
which has an average accuracy of ±0.58 percent for 4000 < Re < 108
and 0 < e/D < 5 X KT2 [34].
Jain [40] proposed
(p)
Swamee and Jain [41] proposed
/ = 7 / ^ ^ C^M, (2-2D
5 .72 \ V
e°>}\
Re°
e 5.02 e 5.02 e
^ log ^ ^
Re log (3 JD
•^T = -2 log ^ ^ ^ ^ ^
3 JD Re 3 JD
(2.23)
which is rather complicated but extremely accurate for 4,000 < Re <
108, andO < (e/D) < 5 X 10"2, and thus gives an excellent explicit form
of the Colebrook-White equation.
Haaland [44] proposed
6 9
• ' - " " (2.24)
which has an average accuracy of ±0.33 percent for 4,000 < Re < 108,
and 0 < (e/D) < 5 X 10"2.
Chen [34] proposed
which has an average accuracy of ±0.30 percent for 4,000 < Re < 108
and 0 < (e/D) < 5 X 10"2.
36 FRICTIONAL HEAD LOSS IN PIPES
(2.26)
(A + BY
where
1
A = 2.457 In 0 9
7\ - 0.27 e
,Re) ) D
and
B= /37,530y
e/D == 0.01
4
10 0.04313 0.04232 0.04235 0.04393 0.04406 0.04398 0.04404 0.04313 0.04304 0.04398 0.04404
106 0.03796 ' 0.03772 0.03872 0.03801 0.03801 0.03796 0.03801 0.03796 0.03804 0.03801 0.03799
108 0.03790 0.03766 0.03865 0.03791 0.03791 0.03785 0.03791 0.03790 0.03798 0.03791 0.03788
e/D = 0.001
104 0.03238 0.03263 0.03179 0.03249 0.03270 0.03263 0.03267 0.03237 0.03217 0.03257 0.03269
106 0.01994 0.02067 0.02099 0.02004 0.02003 0.02001 0.02003 0.01994 0.01994 0.02004 0.02002
10s 0.01964 0.02043 0.02043 0.01964 0.01964 0.01962 0.01964 0.01964 0.01968 0.01964 0.01963
e/D = 0.0001
4
10 0.03104 0.03120 0.03176 0.03096 0.03118 0.03111 0.03115 0.03102 0.03099 0.03105 0.03118
106 0.01344 0.01343 0.01415 0.01354 0.01351 0.01350 0.01351 0.01344 0.01333 0.01353 0.01351
103 0.01200 0.01244 0.01215 0.01201 0.01201 0.01200 0.01201 0.01200 0.01202 0.01201 0.01201
FRICTIONAL HEAD LOSS IN PIPES
3o
h
10.68 LQl
f =
in which L and D are in meters and Q in cubic meters per second. In prac-
tice, the pipe diameter D may also be given in millimeters or centimeters
and the discharge Q in liters per minute, million liters per day and several
other units. Instead of converting the pipe diameter to meters and pipe dis-
charge to cubic meters per second for the entire network, it may be ad-
visable to replace 10.68 in Equation (2.28) by an appropriate value. Such
values are given in Table 2.3.
The Hazen-Williams formula is derived for R = 0.3 m and 5 = 1/1000.
Therefore, the HW coefficient CHW is dependent on the R and 5 values
and also the flow conditions. The CMW values for common pipe materials
as recommended by Lamont [15,30] after examining 372 records are given
in Table 2.4. These values are for a velocity of flow of about 0.9 m/s.
These values should be corrected as shown in Table 2.5 when the veloc-
ities are considerably above or below 0.9 m/s. Lamont emphasizes that the
HW formula is not suitable when the CHW values are appreciably be-
low 100.
The HW formula is more suitable for smooth pipes and therefore for
medium to large diameter new pipes [30]. Although it is not quite suitable
for old pipes, the formula is used in practice by reducing the CHW values,
as described later in Section 2.J2.
Scobey [46] assembled results of more than 1000 determinations made
by himself and others on 147 pipelines ranging in diameters from 25 mm
to 5470 mm and suggested a relationship similar to the HW formula but
with the exponent of S as 0.53 [47]. White [48] on the other hand, from
several tests on the glazed and unglazed pipelines, found an average value
of 0.54 for the exponent of 5, thereby confirming the exponent of S in the
HW formula [48].
TABLE 2.3. Values of Constant in Hazen-Williams Formula
Discharge, in
Pipe diameter,
D, in m3/s m 3 /min m3/h m3/d Us L/min ML/d
3 6 9 5
m 10.68 5.438 x 10" 2.769 x 10" 7.694 x 10" 2.969 x 10" 1.512 x 10-s 2.768 x 10-3
cm 5 .869 x 101° 2.988 x 107 1.522 x 104 42.28 1.631 x 10s- 83.07 1.521 x 107
6 6
mm 4 .351 x 1015 2.215 x 1012 1.128 x 109 3.134 x 10 1.209 x 1Qio 6 .158 x 10 1.128 x 1012
TABLE 2.4. Values of Hazen-Williams Coefficient for Clean New Pipes of Different Materials.
V = 143.5 C R fl 0 6 5 7 5 S 0 5 5 2 5 (2.29)
/ /}1.8099
mm and for velocities ranging from 0.3 m/s to 6.0 m/s. They have pre-
sented a chart for the effect of water temperature on CR values and have
also presented a nomograph for direct estimation of the loss of head for
flow of water at about 20°C O = 1 X 10"6 m2/s) for smooth pipes for
different diameters and discharges [49]. As the use of the modified HW
formula requires determination of CR values which in turn depend upon
the Reynolds number and relative roughness of the pipe, it is preferable to
use the rational, dimensionally homogeneous and simple Darcy-Weisbach
formula for which the coefficient of friction can be explicitly determined
as described earlier.
The usual range of the values of N for common pipe materials is given
in Table 2.6 [13]. Manning formula is more suitable for rough pipes with
N > 0.015 as stated by Lamont [30] and therefore it is more suitable for
open channel flow rather than pipe flow.
Powell [50] and Williams [51] have"shown that the Manning formula, as
is presently used, was not proposed by Manning. Gauckler in 1867, Hagen
PVC 0.008-0.011
Brass, copper, glass, lead, prestressed concrete 0.009-0.012
Concrete 0.010-0.017
Wood stave 0.011 -0.013
Welded steel 0.012-0.013
Coated cast iron • 0.012-0.014
Uncoated cast iron 0.013-0.015
Galvanized iron, riveted steel 0.015-0.017
Comparison of Head Loss Formulas 43
_ fLQ2 _ 10-68L6 1 8 5 2
1
~ 12.ID ~ Cligf2*)4-87
5
from which
013
=
129.2P
129.2P
' /""1.8521 /TO. 148 ^ ' )
Q°
^ ~ C\^2Re° 148
£>° °18 (2
' -*
_
10.29N ;
' ~ 12.ID5 ~ D1
44 FRICTIONAL HEAD LOSS IN PIPES
from which
/ = mDuT ( 2 - 36 )
From Equation (2.35) it is clear that/depends upon Re to some extent
and D to a lesser extent in addition to the CHW value. Therefore, the prac-
tice of using the same value of CHW for different diameters and discharges
is not proper. However, if the CHW value is selected using Table 2.4 and
modified for velocity using Table 2.5, better results can be obtained. Simi-
larly, the modified HW formula would also give better results.
Equation (2.36) shows that the Manning's coefficient N is a function of
the pipe diameter and is independent of Re. Therefore, it would be proper
to use the Manning formula for rough turbulent flow only.
The frictional head loss in a pipe can be expressed by a general head loss
formula,
h, = RQ" (2.37)
Head Loss
Formula Ra n
Darcy-Weisbach — or
12.1 D 5
10 68L 1
Hazen-Williams • 18S2n487 -852
.. . 10.29A/2/.
Manning — 2
a
R value is for L and D in meters and Q in cube meters per second.
Simple Pipe Flow Problems 45
When the Hazen-Williams or Manning head loss formulas are used, the
solution is quite simple because the unknown parameter hf, Q, or D can
be determined directly by using Equation (2.28) or Equation (2.33). How-
ever, when the Darcy-Weisbach head loss formula is used, the solution may
become complicated as the coefficient of friction / depends on the dis-
charge Q and diameter D. Therefore, even though the determination of hf
is simple as Q and D and therefore/can be easily determined, the deter-
mination of Q or D is rather complicated as / is unknown, and therefore
a trial and error procedure may be necessary.
Here Q, L, D, the fluid and therefore the kinematic viscosity v, and the
pipe material and therefore e, are known, and hf is unknown. Therefore,
using the known parameters, the Reynolds number Re and relative
roughness elD are determined. Using one of the explicit formulas, the
coefficient of friction/is determined. Later, using Equation (2.5) or Equa-
tion (2.6), the head loss h, is determined.
Example 2.1
Solution
v = Q= 5/3600
A
9 .302 x 1 0 - nWs
Q 860
„ 64 64
and
Example 2.2
Water flows through a 300-mm diameter, l-km long new cast iron pipe.
If the discharge is 0.1 m3/s, determine the head loss in the pipe by all three
head loss formulas.
Solution
e 0.26
= 8.667 x 10"4
D 300
_ 4.351 x 10' 5 LQ 1 8 5 2
f /""1.852 r»4.87
= 10.69 m
It is seen that, for this example, the DW and HW formulas give nearly
the same head loss while the Manning formula gives a higher value.
48 FRICTIONAL HEAD LOSS IN PIPES
For this case, hf, L, D, v and e are known and the discharge Q is un-
known. From the known e and D values, the relative roughness elD is de-
termined. Rearranging the Darcy-Weisbach formula, Equation (2.5),
(2.38)
Re =
•KDV
or
(2.39)
(2.40)
3.7 D
which is now explicit and gives the value of/directly. This value of/ is
substituted in Equation (2.38) to determine the discharge Q.
Example 2.3
Solution
From Table 2.1, for PVC pipe, e = 0.0021 mm = 2.1 X 10"6 m. Fur-
ther, v = 1.12 X 10"6 m2/s, L = 2000 m, D = 0.15 m, and h, = 30 m.
Substituting these values in Equation (2.40),
Simple Pipe Flow Problems 49
x
* " " ~" = 0.02994, say 0.03 m3/s
* 4
For this case, hf, Q, L, v and e are known and the diameter D is un-
known. Since the diameter D is unknown, both the relative roughness elD
and the Reynolds number Re are unknown. Therefore, the coefficient of
friction/is also unknown.
From a rearrangement of Equation (2.5)
<241)
to = X I - § (2.42,
Example 2.4
Find the diameter of a new 1500-m long cast iron pipe to carry water at
a rate of 3000 L/min, if the permissible head loss is 20 m.
Solution
3000
Q = 3000 L/min = ^ ^ = 0.05 m'/s
D5 _ 8 X 1500 X (0.05)2 .
7T2 X 9.81 X 20 *
or
D* = 0.01549/ (1)
_ 40 1 ' _ 4 x 0.05 \_
e
~ i r i > D ~ - i r X \ x \0'b D
or
f
= ' I 0.00026025 5.74 V? = ° ° 2 1 8 9
8 ++
° \ 33.7
7 xx 0.1987
01987 320,380O9/J
320,
Head Loss Due to Uniformly Decreasing Discharge 51
l6/ r
e,-<e,-a>f
Integrating and simplifying, the head loss in the entire pipe of length L
is given by
L 16/L
\ - Qi\ G. + Qi
52 FRICTIONAL HEAD LOSS IN PIPES
Laterals.
(a)
1,
Qi-K -rQ
f
(b)
Figure 2.3 Head loss in pipe with uniformly decreasing discharge: (a) pipe with service connec-
tions; (b) pipe with uniformly decreasing discharge.
the average discharge. The head loss for the pipe in this case would be
Example 2.5
Solution
/
= r / o.ooi3025 5.74^p = 002277
[ °8\ 3.7 +
139,993°9 j j
Using Equation (2.43),
_ 16 X 0.02277 X 1000
' ~ 7T2 X 19.62 X (0.2) 5
= 2.846 m
Herein the difference between the two head loss values is insignificant.
The water carrying capacity of pipes decreases with age. Cast iron (CI)
pipes are subject to corrosion and even smooth, cement-lined and bitumen-
lined pipes are affected by the formation of slime. In general, the loss in
carrying capacity is caused by (1) a decrease in the area of the cross section
due to the accumulation of deposits on the interior of the pipes, (2) an in-
crease in the roughness, or (3) by both. Biological growth, silting, in-
crustation, and tuberculation are the most common forms of such deposits,
which may vary in thickness from 1 to more than 10 mm.
For the annual growth rate of roughness in old CI pipes, Colebrook and
White [53] suggested a formula in terms of the pH of water. Their formula
in SI units is
tween the annual growth rate and the Langlier saturation index, which in
turn depends upon the pH, alkalinity, and the calcium content of water.
The growth rate of roughness as observed by Lamont [15,30] in coated cast
iron pipes is given in Table 2.8.
It must, however, be remembered that it is not possible to precisely pre-
dict the annual growth rate. Figure 2.4 shows the effect of age on the HW
coefficients. The solid curves show the anticipated HW coefficients as
given by the HW tables of Williams and Hazen [45]. Data from the ob-
servations of the Committee, New England Water Works Association
[54], the values on mains in the Los Angeles water system [55], and the
average curve [56] are shown in the figure. Some observations taken on
water mains in India as reported by Ghatpande and Dhekane [57], Raman
et al. [58], and Dhabadgaonkar and Thakre [59] are also shown in the
figure. The HW coefficients as calculated by Lamont [15,30] are given in
Table 2.9.
Hudson [60] studied several pipes in seven cities of U.S.A. and con-
cluded . . .
(4) Softening of water results in deposits and loss in carrying capacity.
(5) Alignment of a pipe is a major factor determining carrying capacity.
140
100mm
10 20 - 30 40 50 60 70
Pipe Age in Years
30 years old
1 slight 100 106 112 117 120
2 moderate 83 90 97 102 107
3 appreciable 59 70 78 83 89
4 severe 41 50 58 66 73
60 years old
1 slight 90 97 102 107 112
2 moderate 69 79 85 92 96
3 appreciable 49 58 66 72 78
4 severe 30 39 48 56 62
100 years;old
1 slight 81 89 95 100 104
2 moderate 61 70 78 83 89
3 appreciable 40 49 57 64 71
4 severe 21 30 39 48 54
(6) The CHW value depends upon so many unknown factors that it is im-
possible to estimate it accurately.
(7) A head-loss test must be conducted on each section of main to deter-
mine the working coefficient of carrying capacity in order to have an
accurate design coefficient to use in solving hydraulic problems.
It is thus evident that the correct estimation of the CHW value (and natu-
rally the/and N values) for old pipes is extremely difficult and would re-
quire a large number of field observations to arrive at good results.
Example 2.6
A 200-mm diameter, 1-km long coated cast-iron pipe connects two res-
ervoirs that have a 20 m difference in their water levels. Estimate the dis-
charge through the pipe when the pipe is new, 30 years old, 60 years old,
and 100 years old. Use DW and HW formulas. Assume v for
water = 1 x 10"6 m2/s and the degree of corrosive attack of water as
moderate.
Solution: Darcy-Weisbach Formula
• When the pipe is new, from Table 2.1, e = 0.13 mm and therefore
Q = 0.06436 m3/s.
• When the pipe is 30 years old, from Table 2.8, e = 2.41 mm,
and therefore Q = 0.04368 nrVs
• When the pipe is 60 years old, from Table 2.8, e = 4.70 mm, and
therefore Q = 0.03862 m3/s
• When the pipe is 100 years old, from Table 2.8, e = 7.75 mm,
and therefore Q = 0.03482 m3/s
Solution: Hazen-Williams Formula
The values of Q given by the DW and HW formulas are nearly the same.
In the DW formula, the relationship between Q and e is not simple, but as
e increases, / increases, and therefore Q decreases. In the HW formula,
however, Q oc CHW as seen from Equation 2.28.
2.13 PROBLEMS
3.1 INTRODUCTION
By making several assumptions Borda observed that the minor head loss
due to sudden enlargement (Figure 3.1), i.e., hse, may be theoretically ex-
pressed as
(K, - v2y
K = z (3.2)
2g
in which K,, V2 = velocity in smaller and larger pipes, respectively [4].
Several investigators including Brightmore as stated by King and Brater
[13] and Archer [61] investigated this head loss experimentally. They
observed that Borda's formula gives lower head loss values than those de-
termined experimentally for lower velocities and smaller differences in di-
ameters and higher values for the opposite conditions. Therefore, Archer
[61] suggested another formula, from which King and Brater [13] estimated
values of the coefficient for sudden enlargement Kse for different pipe di-
ameter ratios D2IDX for use in the formula
K = Kse p - (3.3)
Vl2 — V22
K =^ 2g (3.4)
—V, ~~ V2
Diameter, D,
1. Sudden Enlargement = K —
/ \l\ - V\
Gradual Enlargement hqe = Kg
\ " 2g )
Exit
\ ••" ~ " 2g
3ell-mouthed 0.1
Sharp-cornered 1
Pipe into still water or air 1
(free discharge) (V2 - 0)
(continued)
61
TABLE 3.1. (continued).
5. Gradual Contraction = K :
2g
Ordinary 0.25
Bell-mouthed 0.10
Streamlined 0.04
Bell-mouthed 0.04
Slightly rounded 0.23
Sharp edged 0.5
Re-entrant (pipe projecting into tank) 1.0
Strainer and foot valve 2.5
a. Elbows—45°
Flanged—regular 0.20-0.30
Flanged—long radius 0.18-0.20
Screwed—regular 0.30-0.42
b. Elbows-90°
Flanged—regular 0.21-0.30
Flanged—long radius 0.18-0.20
Intersection of two cylinders
(welded pipe- not rounded) 1.25-1.8
Screwed—short radius 0.9
Screwed—medium radius 0.75
Screwed—long radius • 0.6
c. Miter Bends
Angle of bend, 9
5° 0.016- 0.024
10° 0.034- 0.044
15° 0.042- 0.062
22.5° 0.066- •0.154
• 30° 0.130- •0.165
45° 0.236- 0.320
60° 0.471- 0.684
90° 1.129- 1.265
d. Return Bend (2 Nos. 90°)
Flanged—regular 0.38
Flanged—long radius 0.25
Screwed 2.2
62
TABLE 3.1. (continued).
8. Tees
/ V2 \
9. Obstructions h0 - Ko —
I 2gj
Pipe-area to flow-area-at-
obstruction ratio
1.1 0.21
1.4 1.15
1.6 2.40
2.0 5.55
3.0 15.0
4.0 27.3
5.0 42.0
6.0 57.0
7.0 72.5
10.0 121.0
0. Flow Meters
a. Venturimeters
Throat to Inlet Diameter Ratio Long Tube Type Short Tube Type
63
TABLE 3.1. (continued).
1 V 2 \
11. Valves KV
T)
a. Gate Valves
fully open 0.19
three-fourths open 1.15
one-half open 5.6
one-quarter open 24.0
b. Butterfly Valves
Closure angle, 9
0° (fully open) 0.3
10° 0.46
20° 1.38
30° 3.6
40° 10
50° 31
60° 94
c. Diaphragm Valves
fully open 2.3
three-fourths open 2.6
one-half open 4.3
one-quarter open 21.0
d. Plug Globe or Stop Valve
fully open 4.0
three-fourths open 4.6
one-half open 6.4
one-quarter open 780.0
e. Check (Reflux) Valves
Swing check (fully open) 2.5
Ball type (fully open) 2.5-3.5
Horizontal lift type 8-12
f. Foot Valve with Strainer 2.5
g: Pressure Reducing Valve 10.0
64
Sudden Contraction 65
Based on the experimental study of Gibson [62] and Parker [63], King
and Brater [13] estimated the values of Kge for different values of the angle
of enlargement 0 . Some typical values are given in Table 3.1.
3.4 EXIT
V1 — V2
hex = Kex^~ (3.5)
in which K,x = head loss coefficient for exit. The value of Kex depends on
the type of exit, as shown in Table 3.1. When a pipe discharges into still
water or air, the kinetic energy of the outgoing water is lost, and therefore,
K.x = 1 and V2 = 0 in Equation (3.5), giving
(3.6)
—
' --.
• ^
rVena contracta,
—v, —V2
A2
m—"
Area, A,
Figure 3.3 Sudden contraction.
coefficient for sudden contraction. The value of C and therefore of Kst de-
pends mostly on the diameters of the small and large pipes.
Using the experimental results of Brightmore. and suggestion of Parker
[63], King and Brater [13] estimated the values of K,c for different values
of diameter ratios. Some typical values are given in Table 3.1.
V\
K = *„, ~ (3-8)
2H
in which Kgc = head loss coefficient for gradual contraction, and V2 =
velocity at the small end. The values of KKU for different types of reducers
are given in Table 3.1.
3.7 ENTRANCE
V2
/»,„ = K,,,, . (3.9)
2
#
in which Ken = head loss coefficient for the pipe entrance, and V = aver-
age velocity in the pipe. The head loss at the entrance is due to separation
and therefore depends on the shape of the entrance. Head loss coefficients
for different entrance shapes are given in Table 3.1.
Tees 67
Even though the magnitude of the average velocity remains the same in
bends and elbows, the direction of velocity changes. This change requires
a radial inward force to provide radial inward acceleration. This causes
distortion in the velocity distribution which in turn produces additional
turbulence as well as double spiral or secondary motion that may persist
for a downstream distance up to 100 pipe diameters until dissipation due
to viscous forces takes place.
Minor head loss in bends and elbows. hh. is given by
hb = K,X (3.10)
in which Kb = head loss coefficient for bends and elbows, and V = aver-
age velocity in the pipe. The value of K,, depends upon the angle of the
bend 9 , the ratio of the radius of the bend R to the diameter of the pipe
D (i.e., on relative radius RID), and the relative roughness of the pipe elD.
The value of Kb increases with the angle of the bend and the relative rough-
ness, but decreases with the RID ratio. The Reynolds number has prac-
tically no influence on the coefficient.
Several investigators such as Fuller. Beij. Anderson, and Straub as
stated by King and Brater [13] have carried out experiments and provided
the values of Kb. Ito [64,65] provided information for flow and head loss
in bends. Jeppson [14] provided curves showing the effect of the relative
radius RID and the angle of bend 9 on the value of K,,. Values of Kh for
some typical bends and elbows as compiled by Amirtharajah [66] are
given in Table 3.1.
3.9 TEES
The head loss in tees depends on the division of flow and the extent of
the change in direction of flow. The head loss in tees, /?,, can be expressed
as
V1
h, = K,— (3.11)
2#
in which K, = head loss coefficient for lees. Values of A", for some typical
tees are given in Table 3.1 [66].
68 MINOR HEAD LOSS IN PIPES
3.10 OBSTRUCTIONS
in which Kfm = head loss coefficient for the flow meters. Values of Kfm for
some typical venturimeters and orifice meters are given in Table 3.1 [66].
3.12 VALVES
Valves are provided in pipelines to control the rate of flow, the direction
of flow, or the pressure. The head loss through valves, hv, can be expressed
as
V1
hv = K.— (3.14)
2
S
in which Kv = head loss coefficient for valves. Values of Kv for gate
valves, butterfly valves, diaphragm valves, plug globe or stop valves,
check valves, foot valves, and pressure reducing valves are given in Table
3.1 [66-68].
Equivalent Pipes
4.1 INTRODUCTION
Q. = a , x = 1, 2, . . . (4.1)
(2) The head loss in the equivalent pipe is the sum of the head losses in
the individual pipes. Therefore, neglecting minor head losses,
K = K + K + ••• = £ K (4-2)
X
69
70 EQUIVALENT PIPES
Li L2 U
U U U
ht, hf, K
Oi = O2 = O3
feLe y fxLx
-QS - L -fir (4-3)
i-'X /A A-.
' Wk _ y mk
Two types of problems arise in series pipe systems: (1) determining the
total head loss and the equivalent pipe when the discharge is known and
(2) determining the discharge through the system and the equivalent pipe
when the total head loss is known.
Since the discharge is known, the head loss in individual pipes and
flierrfore the total head loss can be easily obtained. Similarly, the length
Pipes in Series 71
Example 4.1
* L • i f •0M0867
422J20
"SL
= x x 0.15 x iSL x 10- =
Using the Zigrang-Sylvester relationship, Equation (2.23), / 300 =
0.02043, and/,so = 0.02303.
(1) Head loss in the system
h, = hfl + hfl
(2) For the 200-mm diameter equivalent new cast iron pipe,
e 0.26
= 0.0013
D 200
Le = 471.6 m
giving
D, = 0.3930/° 2 (1)
4 x 0.05 63408
•KD. X 1.004 X 10"' De
/ 200 100 \ _
- n ?487
((0.3) 4 8 7 +
(0.15) 4 8 7 / ~
When the total head loss in a series system is known, for the DW head
loss formula, a trial and error procedure is necessary to determine the
head loss in individual pipes, the discharge in the system, and the length
or the diameter of the equivalent pipe.
Example 4.2
In the series system of Example 4.1, the total head loss is 10.00 m. Find
(1) the discharge through the system, (2) the length of an equivalent new
cast iron pipe of 200 mm diameter, and (3) the diameter of a 300-m long
equivalent new cast iron pipe. Use the Zigrang-Sylvester explicit relation-
ship for estimating/in the Darcy-Weisbach head loss formula. Also find all
the answers using theHazen-Williams head loss formula. Assume CHW =
130 for all pipes.
h, = hh + hfl
• , Q _ /• X 200 x g 2 / 2 x 100 x Q1
12.1 X (0.3) 5 12.1 X (0.15) 5
(2) The equivalent pipe carries the same discharge as the series combi-
nation.
.-. Qe = 0.06162 nrVs
= 0.0013
D 200
and
/ = 0.02158
But
0.02158 x L , x (0.06162)2
= 10
12.1 x (0.2)5
which gives
• Le = 472.5 m
/ x 300 x (0.06162)2
1 _
' i2.i x m
which gives
De = 0.3933/? 2 (1)
and
4 x 0.06162 78144
_ 6 {
e
~ I x 1.004 x 10" De '
which gives
Q = 0.06864m3/s
as in Example 4.1
(2) The discharge in the equivalent pipe is the sum of the discharges in
the individual pipes. Thus,
a (4-8)
». L 2 , f 2 , Q 2 , h f ,
C H w,.£>, 21
£0.54 -E L?-
(4.10)
D8X'3
NM n =E
~ u
NxLxn
(4.11)
-m, (4.12)
loss is known and (2) determining individual pipe discharges, the head
loss, and dimensions of the equivalent pipe when the total discharge is
known.
Using the known head loss and other parameters, the discharge in each
pipe is calculated. A trial and error procedure is necessary if the DW head
loss formula is used. These discharges, when added together, give the total
discharge. From this, the equivalent pipe dimensions can be obtained.
Example 4.3
15
% - 1 = °<*>
Taking/, = 0.02 for the first trial,
0.02 X 300 X Q\
12.1 X (0.2) 5
i. - i =ommi
Q = 0.4451 nrVs; Re = 1.889.000;/= 0.01806: and/.. = 82.18 m
Similarly,
(2) For the equivalent pipe, the diameter = 300 mm, and C,m = 130.
This gives Le = 94.13 m.
4.3.2 TOTAL DISCHARGE KNOWN
(4) Determine Re, elD and then/for each pipe and then go to step 2
(5) Terminate the procedure when/values stabilise
(6) Determine Qu Q2, . . . ; and h,
Example 4.4
In the 3-pipe parallel system of Example 4.3, the total discharge through
the system is 18 ML/d. Find the distribution of flow and the head loss us-
ing the DW head loss formula. Take v for water = 1 mnWs.
Pipes in Parallel 79
Find the answers using the HW formula with HW coefficients for the
three pipes as given in lixample 4.3. Also find the HW coefficient of a
250-m long equivalent pipe of 300 mm diameter.
Solution
300 25 40
Q± _ ° Qj - ° @'
(0.2f ~ (0.3r ~ (0.25) s
giving
r5.3 = °- 0 3 7 6 5 mJ;s
giving
Example 4.5
Find the diameter of a 900-m long equivalent pipe (CHW = 100) to re-
place the series-parallel system shown in Figure: Example 4.5. The length,
diameter and CHW coefficient values for the pipes are as follows: Pipe
1—300 m, 250 mm, 120; Pipe 2 ^ 0 0 m, 300 mm, 130; Pipe 3—200 m,
200 mm, 100; Pipe 4—500 m, 400 mm, 130; and Pipe 5—300 m, 250
mm, 80.
Minor Loss Elements 81
Solution
900 300
(120) 1M2
,852/n "><;\4.87
(0.25)
400 200
(130)1 8 " (0.3) 487 (100)' 852 (0.2) 4 8
' s * /** y\ v< oci /A ^\ \ A on •
9 0 0 5 0 0 300
(100)1-852 Dt*7 ~ (130) 1852 (0.4) 4 87 ++
(80) 1852 (0.25) 4 87
giving De = 0.2839 m
Now two pipes, each 900-m long, having CHW = 100 and of diameters
0.2495 m and 0.2938 m are connected in parallel. Therefore from Equa-
tion (4.10),
or
63
m = (0.2495) 2 " + (0.2839) 2 •"
giving De = 0.3483 m
When minor loss elements such as valves, meters, and bends are present
in a pipeline, the head loss in them can be considered by the methods
82 EQUIVALENT PIPES
KYL _ f^k YL
2g ~ D " 2g
which gives
AL = f (4.13)
or
AL =
SKQ2 _ IO.29N2ALQ2
ir2gir ~ D 16/3
giving
KD*/3
A L < 4 I 5 )
It is thus seen, that for Manning's head loss formula, the equivalent pipe
length is independent of the discharge.
AL = 35KD (4.16)
provides a good approximation for equivalent pipe lengths for minor loss
elements.
Example 4.6
A 300-m long, 200-mm diameter new cast iron pipe contains a globe
valve (K = 10) and a meter (K = 2). If the discharge in the pipe is 3000
L/min, determine the equivalent pipe length using DW (e = 0.26 mm),
TABLE 4.1. Equivalent Pipe Lengths.
Ratio of Equivalent
Pipe Length
Appurtenance to Diameter
1. Sudden Enlargement
D,ID2 = 3:4 7
1:2 22
1:4 31
2. Sudden Contraction
D,ID2 = 2:1 12
3. Entrance: Sharp-Edged 18
4. Bends and Elbows
a. Elbows—45°
Flanged—regular
Flanged—long radius
Screwed—regular
b. Elbows—90°
Flanged—regular 24
Flanged—long radius 18
Cast—regular 32
Cast—long radius 20
c. Welded-90°
bend radius
= 0.5 38
pipe diameter
1.0 17
1.5 12
2.0 10
3.0 9
d. Return bend—180°
regular 73
large radius 39
Tee
Line flow 60
Branch flow 300
6. Valve
a. Gate valve
Fully open 6.7
Three-fourths open 39
One-half open 190
One-quarter open 800
b. Globe valve 320
c. Swing check 80
Minor Loss Elements 85
Solution
3000
Q = 3000 L/min = ^ ^ = 0.05 mVs
and
Re = 4X0.05 =
7T x 0.2 x 1 x 10-6
A , 12 X 0.2 4 ' 3
AL = =u2Jm
124.5 x (omoy
so that the equivalent length = 412.7 m. When the discharge doubles, the
corresponding answers are 412.4 m, 433.8 m, and 412.7 m, respectively.
Using Equation (4.16), L = 35 X 12 X 0.2 = 84 m, so that the
86 EQUIVALENT PIPES
equivalent pipe length = 384 m. Using Table 4.1, and assuming that the
equivalent pipe lengths are proportional to the K values. /. = (320 +
64) X 0.2 = 76.8 m, so that the equivalent pipe length = 376.8 ni.
4.6 PROBLEMS
Problem 4.1: Three new CI pipes with the following lengths and di-
ameters: 200 m, 200 mm; 300 m. 150 mm; and 400 m, 250 mm, are con-
nected in series. If this pipeline carries water at a rate of 4,000 L/min.
find: (1) head loss in the system. (2) the length of an equivalent new CI
pipe 300 mm in diameter, and (3) the diameter of a 900-m long equivalent
new CI pipe. Solve the example using the DW, HW, and Manning formu-
las.
Problem 4.2: In the series system of Problem 4.1, the total head loss
is 20 m. Find: (1) the discharge through the system, (2) the length of
an equivalent new CI pipe 300 mm in diameter, and (3) the diameter of a
900-m long equivalent new CI pipe. Solve the example using the DW, HW,
and Manning formulas.
Problem 4.3: The three pipes of Problem 4.1 are connected in parallel.
If the head loss through the system is 15 m, find: (1) the total discharge
through the system, (2) the length of an equivalent new CI pipe 300 mm
in diameter, and (3) the diameter of a 400-m long equivalent new CI pipe.
Solve the example using the DW and HW formulas.
Problem 4.4: In the 3-pipe parallel system of Problem 4.3, the total dis-
charge is 0.18 m 3 /s. Determine the distribution of flow and the head loss
using the DW and HW formulas.
Problem4.5: The series-parallel system shown by Figure: Problem 4.5
consists of 6 pipes having the following length, diameter and HW coeffi-
cient: Pipe 7—300 m, 250 mm, 120: Pipe 2—400 m, 200 mm, 110; Pipe
3—600 m, 200 mm, 100: Pipe 4—200 m. 200 mm, 90; Pipe 5—300 m,
250 mm, 110; and Pipe 6—250 m. 200 mm, 120. Find: (1) the diameter of
an equivalent 600-m long pipe with an HW coefficient of 130 and (2) the
HW coefficient of a 600-m long, 300-mm diameter equivalent pipe.
Problem 4.6: Solve Problem 4.5 if the six pipes are connected as
shown in Figure: Problem 4.6.
Problems 87
Problem 4.7: In the series-pipe system of Problem 4.1, the first pipe
contains a globe valve (K = 10), the second pipe contains a meter (K =
2). and third pipe contains some minor loss elements (total A' = 25). If
the pipe line carries water at a rate of 4000 I./min. rind: (I) the head loss
in the system and (2) the diameter of a 900-m long equivalent new CI pipe.
Solve the example using the DW and HW formulas.
Reservoirs, Pumps, and Special Valves
5.1 RESERVOIRS
Example 5.1
JW.L. 120.00 m
_ W . L 100.00 m
and discharges into the atmosphere at point 3. The length, diameter and
CHW coefficient values for the pipes are as follows: Pipe 1—300 m, 300
mm, 100; Pipe 2—150 m, 200 mm, 130; and Pipe 3—200 m, 300 mm,
100. Determine the discharges in pipes 1 and 2 and the hydraulic gradient
level (HGL) values at points J and 3, when the outflow at point 3 is: (1) 0.2
m3/s and (2) 0.4 nrVs.
Reservoirs 91
Solution
For the pipe diameter in millimeters and discharge in cubic meters per
second, using Table 2.3, the resistance constant for pipe 1,
•«*•-*
Case 1
As q3 = 0.2 m3/s < 0.2720 m J /s, reservoir 2 receives water. Let the
discharge in pipe 1 be Q. From the flow continuity at 7, the discharge in
pipe 2, from 7 to 2, i.e., Q2 = Q - 0.2. Let the HGL at 7 be //,.
For pipe 1, 222.910 1852 = 120 - //,; and for pipe 2, 493.90 X
(Q - 0.2) 1852 = Hj - 100. Adding these two equations, 222.910 1852 +
493.90(0 - 0.2) 1852 = 20. (This equation can also be directly obtained
by considering path 1-7-2, so that the algebraic sum of the head losses
in pipes 17 and 72 is equal to the difference in water levels in reservoirs
1 and 2.)
Solving by trial and error, Q = 0.2548 m3/s.
Case 2
As q3 = 0.4 m3/s > 0.2720 mVs, reservoir 2 supplies water. Let the
discharge in pipe 1 be Q so that Q2 (from 2 to 7) = 0.4 — Q. Con-
sidering path 1-7-2, 120 - 222.91 G 1852 + 493.90(0.4 - Q) 1 8 5 2 = 100.
Solving by trial and error, Q — 0.3109 m3/s.
IT
Example 5.2
Solution
/120 92\1/1852
ft
- hasr) = a3262 mVs
and
= ° 2 5 7 0 m"/
As g , > <93, / 6 > ft and Q, = Q2 + Q3. Let Q3 = Q and Q2 = mQ
so that <2, = (1 + w)(2. Considering path 1-7-2,
852
222.911(1 + m)QY + 493.9O(m0' 8 5 2 = 28 (1)
(a) (b)
Figure 5.2 Four reservoir system: (a) with one junction point, (b) with two junction points.
Reservoirs 95
This procedure can be extended to systems having more than four reser-
voirs. However, such systems form a particular case of general pipe net-
works, and therefore network analysis procedures, described in later chap-
ters, can better be employed instead of the trial-and-error procedure.
Example 5.3
Solution
/ ioo - 90 \ 0 5
Q2 (from J, to 2) = — = 0 . 1 5 8 1 nWs
\ 400 J
5.2 PUMPS
5.2.1 INTRODUCTION
7 2#
+ E (friction losses)
in which hp — net head delivered by the pump. The velocity heads are
head ~i z ...
6'7 8
Figure 5.3 Typical pump installation.
Pumps 97
-o y
o
QJ /
x: / /
£ .^fTbtal
V)
______^_—-""hea Gloss
>*
CO
Static head
t
Discharge Qp-»-
Figure 5 4 Typical system head-discharge eur\c.
As long as the water level in the two reservoirs [Figure (5.3)] remains
unchanged, the static head remains constant. However, the total head loss
which comprises frictional head losses and minor head losses varies with
the discharge Qp through the pump. Thus, the system head /;, can be ex-
pressed as
A curve showing the relationship between the system head hx and the
discharge Qp is known as the system head-discharge curve. A typical
system head-discharge curve is shown in Figure 5.4.
Figure 5.5 in which the system head-discharge curve of Figure 5.4 is also
shown. The point of intersection of the system-head curve and the pump-
head curve is the actual operating point having discharge Qp that gives
K = hT.
Head-discharge curves for different types of pumps are shown in Figure
5.6. The pump head-discharge curve in Figure 5.6(a) is strictly mono-
tonically decreasing, i.e., the head decreases as the discharge increases.
Thus, for a given head there is only one value for the discharge (point a).
Such a curve is a stable head-discharge curve and is obtained for screw
pumps. The head-discharge curves in Figures 5.6(b) and 5.6(c) which are
obtained for centrifugal pumps and half-axial pumps, respectively, are not
monotonically decreasing. Thus, for the same head, two discharges
[points b and c in Figure 5.6(b)l or three discharges [points d, e, and f in
Figure 5.6(c)] may exist. Such curves having two or more discharges for
the same head are unstable head-discharge curves.
Two typical system curves are also shown in Figure 5.6. In Figure 5.6(a),
there is always one intersection point and thus only one discharge exists for
a particular head. For pump curves in Figure 5.6(b) and 5.6(c), even
though system curve 2 has one intersection point, system curve 1 has two
and three intersection points, respectively. Thus, the condition hp = /?,,
which gives the operating points, gives no insight as to which of these
operating points is more likely to occur. Therefore, a distribution network
having several pumps having unstable head-discharge curves would have
multiple operating points, and the analysis and operation of such networks
would be uncertain as stated by Collins et al. [75,76]. In practice, however,
improvement in the system design with provision of pumps having stable
head-discharge curves (except perhaps at very small percentage of dis-
charge) can ensure a single operating point and thus a unique steady-state
flow condition. Therefore, it is quite common in pipe network analysis to
assume that the pumps have stable head-discharge curves. Therefore, it is
assumed throughout in this text that the head-discharge curves for pumps
are stable giving unique operating points.
5.2.4.1 hp as a Function of Qp
B_
P = QP + — (5.5)
2A
hp = AGl + \ H o - ^ \ (5.6)
5.2.4.2 Qp as a Function of hp
QP = Go - ahbp (5.9)
However, this changes the nature of the curve, i.e., Qp becomes a convex
function of hp instead of being a concave function. Thus, fitting a curve
given by Equation (5.10) would be improper.
Using Equation (5.4) and its modified version, Equation (5.6),
or
Pumps 101
For Equation (5.11), the nature of the curve is maintained and therefore
curve fitting is accurate and convergence is also improved.
If Equation (5.7) is used, a much simpler relationship between Qp and
hp is obtained. Thus,
HP - h,. (5.12)
which maintains the nature of the curve, is fairly accurate, and also im-
proves convergence.
By applying polynomial regression to several pump curves, Fietz [80.
67] observed that a third degree polynomial (hp as a polynomial of Q,) is
usually the best fit. However, Jeppson and Davis |68J have observed that
a second degree polynomial is usually sufficient in practice.
Typical head-discharge data and the corresponding curve are shown in
Figure 5.7(a). The discharge-head curve for the same data is shown in
20
£ . 1 5 --
O
5 10 15 20
Pump head, m
(b)
Figure 5.7 Pump data and curves: (a) head-discharge curve, (b) discharge-head curve.
102 RESERVOIRS, PUMPS. AND SPECIAL VALVES
Figure 5.7(b). It is seen that both the curves are concave. Assuming the
working range of the pump to be 0.1-0.3 m3/s the analytical relationships
described earlier are given below:
Power, kW
T70
Efficiency -60
%
r100 -50
•90 -40
•80 -30
-70 -20
•60 -10
1
•-50
0-05 0-10 0-15 0-20
Discharge or Capacity, n
Figure 5.8 Typical pump characteristic curves.
a
d. PumpBX
•"~-<\ -Combination
t PumpA\^
\
• Capacity •Capacity
(a) (b)
Figure 5.9 Two pump combinations: (a) series, (b) parallel.
104 RESERVOIRS, PUMPS. AND SPECIAL VALVES
40 -i 1 r
Tank full, source at low level •
N 2&3
1&3
1&2
Pump/1
-Tank empty, source at high level
. CD
XI
without exerting excessive pressures in the pipes and pump casings of the
lower floors.
Parallel combinations are commonly used for lifting water that varies in
quantity (from minimum demand to maximum demand) and head (source
at high level and service reservoir empty to source at low level and service
reservoir nearly full). Figure 5.10 shows the head-discharge curves for a
three-unit parallel combination [81,82]. Pumps of different capacities are
installed so that they can be run in different combinations to obtain maxi-
mum efficiency.
Example 5.4
Assuming the working range for the pump to be 0.1-0.3 rrrYs, from Sec-
tion 5.2.4,
For the first trial, assuming the flow in the pipe to be rough turbulent,
from Equation (2.11b),
0.26
- , . \ / = 0.01897
3.7 X 300
0.01897X200
12.1 X (0.3)5
AQ_ 4 X 0.1944
Re = = 825,060
•nDv TT x 0.3 x 1 x 10-6
Delivery Pipe
Suction pipe * - (Pipe 2)
(Pipe
90° elbow 90° elbow
K = 0.9 K = 0.90
W.L. pump Gate Valve
80.00 m fully open
Foot Valve K = 0.19
K = 2.50
The suction pipe is 30 m long, 300 mm in diameter, and the delivery pipe
is 150 m long, 250 mm in diameter. Both pipes are new cast iron with
e = 0.26 mm. Several appurtenances with their K values are located on
the pipeline as shown in the figure. The pump head-discharge data are
given in Table 1: Example 5.5. Develop the system head-discharge curve
106 RESERVOIRS, PUMPS, AND SPECIAL VALVES
Discharge, Head,
m3/s m
0 30.8
0.02 30.9
0.04 30.8
0.06 29.0
0.08 27.0
0.10 23.0
0.12 17.3
0.14 9.0
and determine the steady state head and flow. Use the Zigrang-Sylvester
relationship for determining/for the DW formula. Take v for water = 1.3
mm 2 /s.
Solution
8
> < 3 0 , 8 X 150 1
2 _2 2 w 2
•K X 9.81 X {03f - " 7T n oi
X 9.81 v^
X /r\ IJ\SJ
(0.25) *^
8 8 1
2 r-v4 ( t - v a l v e + Kcibov,) H " — (Aj,a,c.valvc + 2ft e |[, o w + A^.xil) \Q
X 9.8? X
= 97.93g 2
.-. System head, hT = 10.00 + (1020/, + 12691/2 + 97.93)g 2
Pumps 107
30'
25
-§20
<v ating head 15-7 m
/
"Q-15
0- 10 r
5 Operating discharge,0-124rr^7
0 002 004 0-06 0-08 0-10 0-12 0-14
Discharge.mVs
gives the operating point for which head = 15.7 m; and discharge =
0.124 m3/s.
Example 5.6
WL 90.00 m
characteristic curve is as given in Figure 5.7. The pipe and other data are
the same as that in Example 5.1. Determine the different values as in Ex-
ample 5.1, if the outflow at 3 is (1) 0.1 m3/s and (2) 0.4 nrVs.
108 RESERVOIRS, PUMPS, AND SPECIAL VALVES
0, £/»„ K K
m^/s h m m m m m
0 0 0 0 10.00 10.00
0.02 0.02280 0.02286 0.13 0.04 0.17 10.00 10.17
0.04 0.02177 0.02152 0.47 0.16 0.63 10.00 10.63
0.06 0.02053 0.02100 1.04 0.35 1.39 10.00 11.39
0.08 0.02018 0.02073 1.82 0.63 2.45 10.00 12.45
0.10 0.01995 0.02056 2.81 0.98 3.79 10.00 13.79
0.12 0.01980 0.02044 4.03 1.41 5.44 10.00 15.44
0.13 0.01974 0.02040 4.72 1.66 6.38 10.00 16.38
0.14 0.01969 0.02036 5.46 1.92 7.38 10.00 17.38
Solution
hp = 17.19 - 49.20C; 8
Therefore,
Following the procedure outlined for case 1, the direction of flow in pipe
2 is found to be from 2 to J, and the required values obtained are as fol-
lows: Ql = 0.2530 m 3 /s; Q2 = 0.1470 m3/s; //, = 85.84 m; and
H3 = 58.61 m.
Example 5.7
WL 92.00 m
Sump
The pump characteristic curve is as given in Figure 5.7. The pipe and other
data are the same as that in Example 5.2. Determine the different values
as in Example 5.2, when the sump W.L. is (1) 80.00 m and (2) 98.00 m.
Solution
= 80.00, i.e.,
Valves and the head losses through them are considered in Chapter 3.
Though all valves cause loss of head when flow takes place through them,
some valves are simple as far as their inclusion in pipe network analysis
is concerned. For example, globe valves, gate valves, and such other
valves control the magnitude of flow through them. They can be con-
sidered in the analysis by using appropriate head loss coefficients (Table
3.1) depending upon how far open the valves are. It is assumed that the
coefficients remain the same regardless of the magnitude of the discharge
and the direction of flow through them. Some valves, on the other hand,
have multiple actions. They control the direction of flow through them and
reduce the water pressure to a preset value. Check valves and pressure re-
ducing valves are examples of this type and require special treatment in the
analysis. These valves are now considered.
A check valve allows flow in one direction only. It opens when the flow
is in the desired direction and closes when the flow is in the opposite direc-
tion, thus preventing flow. Check valves are mainly used on suction and
delivery sides of pumps. Afoot valve is a check valve fitted on the inlet of
Special Valves 111
w
(a)
Hi>Hj
(b)
Qx=O
(C)
Figure 5.11 Check valve: (a) fitted in pipe connecting nodes / and./, (b) //, > //,-, (c) //, < //,.
the suction pipe of a centrifugal pump. It holds water in the suction pipe
and thus obviates or reduces the priming effort.
Different situations arising due to the presence of a check valve are
shown in Figure 5.11. Figure 5.11(a) shows a check valve fitted in pipe x,
connecting nodes / and j , with HGL values //, and Hj, respectively. The
check valve in Figure 5.11(a) allows flow from node i to node j . When
Hi > Hj, flow takes place from / to j and Qx > 0 [Figure 5.11(b)J. When
Hi < Hj, the check valve closes and prevents flow from j to i so that
& = 0 [Figure 5. ll(c)J.
Example 5.8
1 WLiOO.OOm
2
Solution
r>
8A: 8 x 2.5
ir2gD4 7T2 X 9.81 x (0.2) 4
= 129.10
112 RESERVOIRS, PUMPS, AND SPECIAL VALVES
From Example 5.1, reservoir 2 starts supplying water when q3 > 0.2720
rrrVs.
Case 1
Case 2
" 2 5 8
Pipe,i 4 7 10
3 6 9
8 7 6
Figure 5.12 A network.
Hset
d
(Q)
i Hi>H S e t »H d >H j
(b)
9-
i Hi>H d >Hset>Hj j
(c)
o- CU=O c
i Hi>Hj>H Se t J
(d)
9-
Hi<Hj
(e)
Figure 5.13 Pressure reducing valve: (a) fitted in pipe connecting nodes 1 and ;, (b) H, >
««, a Hd > Hj, (c) //, > Hd> HM > H,, (d) H, > H, > H^, (e) H, < Ht.
114 RESERVOIRS, PUMPS, AND SPECIAL VALVES
The following four modes of operation may occur for the PRV:
(1) Hi > Hsa > Hd > Hj [Figure 5.13(b)]: As H > Hj, flow takes
place from / toy. Further, as Hd < //sel, the PRV is inoperative. If
A"PRv is the head-loss coefficient, the resistance constant of the PRV,
i.e.,
8ATp RV
**PRV —
'H — H \05
u
r» I •'set 1
"PRV = I ^ I >
(3) Hi > Hj > Hxx [Figure 5.13(d)]: It may happen that a PRV is
bypassed so that //, becomes more than //„,,. As //, > //,-, flow
should take place from / toy. However, as the intervening //se, is less
than Hj, the PRV prevents flow, i.e., Qx = 0. The pressures on the
upstream and downstream ends of the PRV increase so that Hu = //,
and Hd = Hj.
(4) Hi < Hj [Figure 5.13(e)]: As //, < //,, flow direction would have
been fromy to / in the absence of a PRV. However, the PRV shuts off
when Hi < Hj and prevents flow from j to /. The PRV behaves as a
check valve so that Qx = 0, //„ = #, and Hd = Hj.
Example 5.9
W.L. 120-00 m
JW.L. 100-00 m
Special Valves 115
data remain the same as in Example 5.1. Determine the discharges in pipes
1 and 2 and all the HGL values for the following conditions:
(1) q3 = 0.4 m3/s; Hxll = 110.00 m; and HM.2 = 96.00 m
(2) q3 = 0.4 nWs; //set_, = 80.00 m; and //,cl2 = 75.00 m
(3) q3 = 0.3 m3/s; HM.X = 110.00 m; and Hxt2 = 80.00 m
For these PRV settings, what is the limiting value of q3 at which
reservoir 2 starts supplying water?
(4) q3 = 0.1 m3/s; HmA = 110.00 m; and Hsa2 = 80.00 m
Solution
10201
*
and
*PRV"2 =
^ x 9.8^ x (0.2)4
Case 1
Case 2
Since the PRV settings are low compared to the respective reservoir
water levels, assume that both the PRVs are operative. Further, from Case
1, Gi + Qi = 0.4. .'. Hdl = 80.00 m, and Hd2 = 75.00 m. For path
dt-J-d2, 80.00 - 148.610! 8 " + 329.27(0.4 - (?,)' 852 = 75.00 which
gives £>, = 0.2669 m 3 /s; and Q2 = 0.1331 m3/s.
For these discharges, if the PRVs were inoperative, Hdt = 120.00 —
74.30 x (0.2669)' 852 - 102.01 X (0.2669)2 = 106.30 m > 80.00 m.
.'. PRV-1 is operative and reduces the pressure from 113.56 m ( = //„,) to
80.00 m. Similarly, Hd2 = 100 - 164.63 X (0.1331)1852 - 516.42 x
(0.1331)2 = 86.92 m > 75.00 m. .'. PRV-2 is also operative and reduces
the pressure from 96.07 m ( = HU2) to 75.00 m.
The assumption that both the PRVs are operative is correct, and there-
fore, £>, = 0.2669 nrVs; Q2 = 0.1331 m3/s. The different HGL values
are: HUI = 113.56 m; Hdl ( = // set ,) = 80.00 m; HU2 = 96.07 m; Hd2
( = // scl2 ) = 75.00 m; Hj = 67.13 m; and H3 = 39.90 m. In this case
both the PRVs are in a pressure reducing condition, i.e., in mode 2 of op-
eration.
Case 3
Case 4
5.4 PROBLEMS
Problem 5.1: The length, diameter, and CHW coefficient for the pipes
in the system of Figure: Example 5.1 are as follows: Pipe 1—400 m, 350
mm, 80; Pipe 2—200 m, 250 mm, 110; and Pipe 3—300 m, 300 mm, 90.
Determine: (1) Demand at point 3 so that the discharge in pipe 2 is zero
and (2) discharges in pipes 1 and 2 and the HGL values at J and 3 when
the demand at 3 is (a) twice the value in (1) and (b) one-half the value in
(1).
Problem 5.2: For the three-reservoir system shown in Figure 5.1,
//, = 125.00 m, H2 = 103.00 m, and H3 = 87.00 m. If the length, diam-
eter, and CHW coefficient for the pipes are as given in Problem 5.1, find
the discharge in each pipe and the HGL at J.
Problem 5.3: In the four-reservoir system of Figure 5.2(a), //, =
130.00 m, H2 = 120.00 m, H3 = 110.00 m, and H4 = 100.00 m. If the
length, diameter, and CHw values for the pipes are: Pipe 1—200 m, 300
mm, 120; Pipe 2—300 m, 250 mm, 100; Pipe 3—400 m, 250 mm, 90;
and Pipe 4—250 m, 200 mm, 110; find the discharge in each pipe and the
HGL value at J.
Problem 54: The reservoirs and pipes of Problem 5.3 are connected to
form the system shown by Figure: 5.2(b). If the water levels in the reser-
voirs are the same as those in Problem 5.3, and if pipe 7, J2 is 300 m long,
300 mm in diameter with CHW = 110, find the discharge in each pipe and
the HGL values at Jt and J2-
Problem 5.5: A sump (water level = 107.20 m) is connected to an
elevated service reservoir (W.L. = 121.20 m) through a 500-m long, 350-
mm diameter new CI pipe (e = 0.26 mm). The pump fitted in this pipe
has discharge-head characteristics as follows: 0 m 3 /s-50.2 m, 0.01
m 3 /s-50.3 m, 0.02 nr7s-50.2 m, 0.03 m 3 /s^9.5 m, 0.04 m3/s-48.1 m, 0.05
m 3 /s-46.4 m, 0.06 m 3 /s-43.5 m, 0.07 nWs^lO.l m, 0.08 m 3 /s-36.4 m, 0.09
118 RESERVOIRS, PUMPS, AND SPECIAL VALVES
6.1 DEFINITIONS
Source
Source
a
A looped network contains loops (Figure 6.4). It must have at least one
source as in a single-source network, but may have more sources as in a
multisource network. A looped network must have at least one sink.
(Why?) The node that serves as the sink depends upon the nodal demands.
For example, node F in Figure 6.4 is a sink. However, if demand at D in-
creases, the flow direction in links DE and DF may reverse and node D
would behave as a sink. Although the direction of flow in a link is fixed for
a particular demand pattern, it may change for another demand pattern.
Thus, for a node, a link which behaves as a supply link for one demand
pattern may behave as a distribution link for another demand pattern. Sim-
ilarly, when the direction of flow changes in a link, the upstream node
becomes the downstream node and vice versa.
Pipe networks seldom consist of loops only. In practice, they are of com-
posite type consisting of loops, branches, and serial parts. However, it is
usually possible to reduce a composite network to a looped network for
analysis, and therefore composite networks are referred to as looped net-
works.
In serial and branching networks, only one path is available for
transporting water from a source to a particular node. Thus, if a link of
such a network is closed for repairs or replacement, all demand nodes sit-
uated downstream of this link are completely cutoff from the source.
Therefore, supply through a serial or branching network is less reliable
when compared to that through -a looped one where alternate paths exist.
However, serial and branching networks are cheaper and usually easier to
analyze than looped networks. Serial and branching networks are used for
industrial water supply purposes, sprinkler irrigation projects, and for dis-
C D
Figure 6.4 Looped network.
Parameters 123
6.3 PARAMETERS
6.3.1 CONFIGURATION
The pipe lengths are obtained from the known geometrical layout of the
network. When nodes are connected by links consisting of pipes in series,
in parallel, and in series-parallel combinations, such links are usually re-
placed by equivalent pipes in network analysis.
The pipe diameters are either known or calculated for equivalent pipes.
However, some unknown pipe diameters can be obtained from the analysis
as explained later.
and N values do not depend on pipe discharges and therefore remain con-
stant during the analysis. The Darcy-Weisbach friction factor/, however, is
a function of the Reynolds number and therefore of the pipe discharge, and
thus it must be reevaluated when the pipe discharge changes.
The nodal demands fluctuate with the time of day, day of the week, and
the season and thus are variable [1,82]. However, it is a common practice
to assume that the demands remain constant in the analysis. Such analysis
is for steady state conditions and is termed static analysis as opposed to
dynamic analysis, in which fluctuations in demand are considered.
As the inflows at the source, nodes depend upon the nodal demands
which are steady in static analysis, the source supply pattern is also steady
in the static analysis.
The hydraulic gradient levels or simply the heads at the demand nodes
are mostly unknown and obtained from the analysis. Therefore, they are
generally treated as unknown parameters.
When service reservoirs are the source nodes, the network is gravity-fed
and the heads at the source nodes are known. Thus, the source heads are
known and remain fixed in the static analysis.
In case of pumped source nodes, the water level in the sumps is known
and is assumed to remain fixed. However, the head supplied by the pump
to the network is a function of the discharge and therefore is not known.
As stated in Section 2.9, the length L, diameter D and the roughness
Parameter Interrelationships 125
The main elements of a hydraulic network are the nodes, pipes, and
loops. Even though they can be labelled arbitrarily, a procedure for their
labelling is suggested herein.
Let a multiple source branching network have M source nodes and N de-
mand nodes. It is preferable to label the source nodes initially from 1 to M
and label the demand nodes later from M + I to M + N (= J). The
pipes are labelled in any convenient order.
Consider a looped network shown in Figure 6.5. Let the network have
M source nodes, labelled 1 . . . M, and N demand nodes, labelled
M + 1, . . . , M + N ( = / ) , so that the node labels are j , j = 1 . . . J
126 TYPES AND PARAMETERS
D ' D
X = J + C - 1 (6.1)
The head loss in a pipe can be expressed by a general head loss relation-
ship, Equation (2.37). Thus, for all pipes in the network,
h , = Ht - Hj = RXQ"X, x = 1 X (6.2)
in which hx = head loss in pipe x; //,-, //, = HGL values at the upstream
node / and downstream nodej of pipe x, respectively; Qx = discharge in
pipex; Rx = resistance constant of pipe x (Table 2.7); and n = an expo-
nent (Table 2.7). Since n * 1 in most cases, Equation (6.2) represents a
nonlinear relationship between hx and Qx.
The direction of flow in a pipe is assumed during the network analysis.
However, when the direction of flow changes, Qx becomes - Qx. There-
fore, to evaluate hx when the flow direction reverses and Qx becomes — Qx,
Equation (6.2) is written as
hx = H, - Hj = 'RX\ Qx (6.3)
Hi > Hj. When the flow direction reverses, Q,, and therefore hx, become
negative, and //, < //,.
Equation (6.3) is sometimes expressed as
K = RL Qx, x = 1 . . . X (6.4)
RL = R,\ G x | - \ J C = 1 . . . X (6.5)
Equation (6.4) is the linearized form of the general head loss relation-
ship, expressing hx as a linear function of Q,. In Equations (6.2) and (6.3)
the resistance constant Rx has a fixed value for a pipe when the Hazen-
Williams or the Manning's head loss formulas are used. However, the
modified resistance constant Rx depends upon Q and therefore it must be
reevaluated when the value of Qx changes.
Equation (6.2) can also be expressed as
To take care of the change in the flow direction, Equation (6.6) can be ex-
pressed as
,x = 1. . .X (6.9)
R\'
Q, + qj = 0,j=l...J (6.10)
J
' R ) + qj = O,j= I . . .J (6.11)
i connected
toy through x
Hi Hj
~ , , , +qj = 0J=l...J (6.12)
/connected Dl/nl LJ U \
J
t o ; through j: ' !
C (H - Hj) + qj = 0 , j = 1 . . . J (6.13)
/ connected
toy through x
For all the loops of a network, the algebraic sum of the head losses in
the pipes forming a loop must be zero. Thus,
2 K = 2 R*Qi = 0, c = I, II . . . C (6.14)
Rules for Solvability of Pipe Networks 129
R
E ^ = E 'Q' = 0, c = I, II . . . C (6.15)
Shamir and Howard were perhaps the first to consider the problem of
solvability of pipe networks when all three types of unknowns, nodal
heads, nodal flows, and pipe characteristics, are present. They initially
suggested certain rules [88] and later updated them [89-90]. Their rules
are as follows:
(1) The total number of the unknowns must be equal to the total number
of the nodes in the network.
(2) At least one nodal head must be known.
(3) A node having an unknown consumption must be connected to at
least one other node with a known consumption.
(4) A subsystem consisting of a pipe with unknown resistance and its
two end nodes, must have no more than one additional unknown—
one of the heads or consumptions at the two end nodes.
(5) Considering any node, at least one of the following must be un-
known: the consumption at the node, the head at the node itself or at
any adjacent node, or the resistance of a pipe connected at the node.
These rules were developed heuristically by considering different com-
binations of unknowns, and no rigorous mathematical proof was provided
for their validity.
DO TYPES AND PARAMETERS
Using graph theory, Gofman and Rodeh [91] also developed certain rules
for solvability of networks having pipes with unknown characteristics.
(They termed such pipes as head generators.) However, in developing their
rules they assumed that a pipe network "virtually always has a unique solu-
tion in the absence of head generators." Therefore, their rules are not gen-
eral as they do not deal with the number and distribution of the unknown
nodal heads and unknown nodal flows. Their rules are as follows:
(1) The number of nodes with known heads must exceed the number of
head generators by one.
(2) For a given network G, there exists a connected network G' such
that:
a. All head generators and nodes with known heads belong to G'
b. G' is loopless
c. A path in G' connecting two nodes of known nodal heads con-
tains at least one head generator
Gofman and Rodeh have rigorously proved the validity of these rules.
However, these rules are not general as stated earlier. If a network is
solvable in the absence of head generators, it would still remain solvable
when the head generators are introduced such that they satisfy these rules.
The unknown number rule is the same as rule 1 of Shamir and Howard.
The minimum-number rule is similar to rule 2 of Shamir and Howard but
is more comprehensive because it also avoids the possibility of the
specified nodal flows being inconsistent. The assignment rule is similar to
rule 5 of Shamir and Howard but is more comprehensive and in combina-
tion with the path rule guarantees that all the available equations are inde-
pendent. Rules 3 and 4 of Shamir and Howard are not essential as shown
in Example 6.1. The path rule has no similar rule in the set of rules pro-
posed by Shamir and Howard, but is similar to the rules proposed by Gof-
man and Rodeh.
132 TYPES AND PARAMETERS
Example 6.1
Node,1
Pipe, 2
5 8
10
13 9 l4
10
8
Solution
Case 2
1 Hi 1,2,4 2,4 4
2 Qi 1 1 — — —
3 H2 1,2,3,5 — 2,5 — 5
4 Ft? 2,3 — 2 2
5 H3 2,3,6,7 — 2,6. — 6
6 93 3 3 — — —
7 R7 4.5 — 4,5 — 4,5
8 9/ 7 7 — — —
9 H8 4,8,9 _. 4,8 4.8
10 99 9 9 — — -
Case 3
1 H, 1,2,4 1,4 1
2 a? 2 2 — — —
3 9s 5 5 — — —
4 Re 3,7 — 3,7 — 3,7
5 96 6 6 — — —
6 H, 3,7,10 — 3,7,10 — 3,7,10
7 fi7 4,5 — 4 4 —
8 H8 4,8,9 4.8 8
133
(continued!
TABLE: EXAMPLE 6.1. (continued).
9 9 9 _ _ _
10 7,10 — 7,10 — 7,10
Case 4
1 Qi 1 1 — — —
2 H2 1,2,3,5 — 2,5 — 5
3 R2 2,3 — 2 2 —
4 H3 2,3,6,7 — 2,6 — 6
' 5 Q3 3 3 — — —
6 fly 4,5 — 4,5 — 4,5
7 9/ 7 7 — — —
8 H8 4,8,9 — 4,8 — 4,8
9 Q9 9 9 — — —
10 Q10 10 10 — — —
TABLE: EXAMPLE 6.1. (continued).
! I I I I I I I
I I I I I I I I
Case 2
4 4
1
5 5
2
6 6
3
nil nil nil
7
Case 3
1 - - - 1
2
_ _ - - - 5
- 3,7 3 - - 3
- — — — — 6
- 3,7,10 7 — — 7
_ _ — — - 4
8 8
c - - - -
in (continued)
TABLE: EXAMPLE 6.1. (continued).
' 9
— 7,10 10 — — 10
Case 4
1
5 I 5
— — — — 2
— — —
—
6
— — —
— 6
3
• —
4 4 —
— 4
—
— — —
_
—
4,8 — 8
_8 7
8
— 9
— — — — — 10
Problems 137
6.7 PROBLEMS
Problem 6.2: In the network of Figure: Problem 6.2, the unknowns are
as follows:
(1) qu Hu R u R2, H2, q3, R12
(2) Hi, H2, R2, q4, Hs, H7, R9
(3) q3, H3, R3, q6, H6, q7, H7
(4) q2, q4, H4, q7, H7, R9, Ri2
(5) q2, H2, q4, H4, q7, H7, R n
Problem 6.3: In the network of Figure: Example 6.1, the unknowns are
as follows:
(1) <7i> <?2, H2, R4, Hs, R6, R7, qs, H9, q10, Ri3
(2) qu q2, H2, R3, R4, Hs, H6, q7, qs, H9
(3) qu Hu q3, qs, Hs, q6, H7, q9, q10, Hl0
(4) qu Hi, q3, qs, H5, q6, H7, Hs, q9, qi0
(5) q2, H2, H4, R4, R6, H7, qs, H9, q10, Rl2
(6) Hi, qu R2, q4, Hs, q6, Hs, R9, qw, R u
(7) Hi, qx, H4, q4, H7, qs, H9, H10, Ri3, Ri4
(8) Hi, q2, q4, R s , Hs, Hs, H9, qi0, Hw, Ri4
(9) q2, H2, H4, R s , q6, H6, q9, qi0, R u , Ri3
Formulation of Equations
7.1 INTRODUCTION
The seventh combination in which some R, some H and all (except one)
q values are known can be looked upon as a particular case of the general
combination 3 and therefore would not be considered separately.
Thus, combinations 1, 2, and 3 are important combinations in the analy-
sis of pipe networks. Even though combinations 1 and 2 are particular
cases of the general combination 3, they are considered separately because
of their common occurrence in practice. This also makes it easier to de-
velop the subject and understand the formulation of the equations and the
methods of their solution. The formulation of equations, therefore, is con-
sidered for combinations 1, 2, and 3, in that order, and is also illustrated
through typical networks. Initially branching networks and subsequently
looped networks are considered for each combination.
As the flows at all the demand nodes are known, starting from the end
nodes and successively applying the node-flow continuity relationship
[Equation (6.10)], the flows and head losses in all the pipes and the supply
at the source node are determined. As the head at the source node is
known, moving along the paths and successively subtracting the head
losses in pipes, the heads at all demand nodes are determined. Thus, it is
quite simple to analyze such a network.
Loop, I II
7.4.2.1 0 Equations
The unknown pipe flows, Qx, x = 1 . . . X are taken as the basic un-
known parameters in formulating the Q equations.
General Formulation. Applying the node-flow continuity relationship
[Equation (6.10)] for the demand nodes
= 2 . . .J (7.1)
x connected
lo j
2 Ro = 0, c = I, II, . . . C (7.2)
x6,
Demand,
HGL, 0-4
1 Pipe, 2
100-00
Q2 R=20 Q5 40 5
SI
a r
a
Loop, I II
2 3 3 6
Q3 1Q Q6 \
0-3 0-1
0-1
Figure 7.2 One source, five-demand-node looped network.
figure. The directions of the unknown pipe flows are assumed as shown in
the figure.
Considering the inflow to a node as positive and outflow from it as
negative, the node-flow continuity relationships for nodes 2 . . . 6 are, re-
spectively,
0 - 0 - Q6 + 0.1 =0 (7.3b)
Qs + Qi - 0.4 = 0 (7.3d)
Q6 - Qi - 0.1 = 0 (7.3e)
Moving in the clockwise direction along the two basic loops I and II, the
loop-head loss relationships yield, respectively
F7 = 4 0 0 - 5 0 0 - 2 0 0 + 3 0 0 = 0 (7.3g)
ing the linearized form. Equation (6.15), the loop-head loss relationships
are
7.4.2.2 H Equations
The unknown nodal heads, //,, j = 2 . . . J are taken as the basic un-
known parameters in formulating H equations.
General Formulation. By applying the node-flow continuity relation-
ship, Equation (6.11), to nodes 2 . . . J, the H equations are
in which //,, //, = heads at upstream node and downstream node, respec-
tively.
The linearized version of the H equation is
C = ,x = 1... X (7.4c)
( l
*»» H H " ^
Kox tli — tlj
"-0.3=0 (7.5a)
146 FORMULATION OF EQUATIONS
(H4 H, , . , „ . „„ , _ Q4 = 0 (? 5d)
40 / \ 50
Equations (7.5a) . . . (7.5e) provide five nonlinear equations for the five
basic unknowns H2 . . . H6.
The linearized version is
F6 = C,'(100 - H2) - Q(H2 - H3) - 0.3 = 0 (7.5f)
Equations (7.5f) . . . (7.5j) provide five linearized equations for the five
basic unknowns H2 . . . H6.
7.4.2.3 AO Equations
Note that as pipe 4 is common to both loops, the corrections AC?, and
AC?n are applied to the assumed discharge 0.05 m3/s in pipe 4.
Equations (7.7a) and (7.7b) provide two nonlinear equations for the two
basic AC? parameters, AC?, and AC?,i-
7.4.2.4 AH Equations
ox
innected toL
i connected
j through x \ I. o)
Illustration. For the network of Figure 7.2, Hoi = 100 m, a given value.
Let the assumed Mj values be: ^H?. — 97 m, ,// 3 = 95 m, {H4 = 94 m,
iHs = 90 m, and ,// 6 = 92 m. The basic unknown parameters are
148 FORMULATION OF EQUATIONS
-0.3=0 (7.9a)
[(94+AH,) -(90+AHs)]""
- 0.2 = 0 (7.9c)
40
[(94 + Atf4) - (90 + AH,)] "" [(92 + A//6) - (90 + A//5)] "-
40 50
L J L J
- 0.4 = 0 (7.9d)
[(95 + AH3) - (92 + A//6)l u" [(92 + A//.) - (90 + A///55; ) l " "
50 J
- 0.1 = 0 (7.9e)
. \\ ?
9 2 11
2 6
4 5
D J
1
3I
10
I
5
3 12
13
12
-I
15 +<
iou +< 1oi5
(a) (b)
Figure 7.3 Multisource branching network: (a) complete network, (b) reduced network.
sistance constants of all pipes are known and they are / ? , ! . . . Rol4 for
pipes 1 . . . 14, respectively. The heads at the source nodes are known and
are Hol . . . Ho3 for source nodes 1 . . . 3, respectively. The demands at
the demand nodes are known and are qo4 . . . qol5, for nodes 4 . . . 15, re-
spectively.
The discharges in pipes 2, 3, 5, 7, 8, 9, 12, 13, and 14 can be easily deter-
mined. These pipes of known discharges are removed from the network,
and the outflows at the intermediate nodes 2, 4, 7, and 12 are correspond-
ingly modified as shown in Figure 7.3b. As q2 and q3 are unknown, the dis-
charges in the pipes of the reduced network are unknown. However, these
discharges must be such that they satisfy the node-flow continuity relation-
ship and completely utilize the available head loss along the paths. Thus,
the head loss along the path 1-4-7-2 must equal Hol — HB2 and that along
path 1-4-7-12-3 must equal HoX — Ho3. Therefore, pseudo pipes connecting
nodes 1 and 2, and 1 and 3 are introduced as shown by dashed lines in
Figure 7.3(b). Thus, two pseudo loops, I and II, are formed as shown in
the figure.
As shown herein, a multisource branching network having M source
nodes can be converted to a looped network having M — \ pseudo loops.
Therefore, the formulation of equations would be similar to that of looped
networks.
C+1
M-1
M+1
C+M-1
M+N
N. 1
7.5.2.1 0 Equations
(<Z<>3 • • • (lob) and //values (Hot, Ho2) are as shown in the figure. The direc-
tions of the unknown pipe flows are assumed as shown in the figure.
From the node-flow continuity relationship for nodes 3 . . . 6, respec-
tively,
F, = 0 - 0 - 0.3 = 0 (7.12a)
Fi = 0 + 0 - 0 - 0.2 = 0 (7.12b)
F3 = 0 + 0 - 0.4 = 0 (7.12c)
FA = 0 - 0 - 0.1 = 0 (7.12d)
F6 = 40(2? - 5 0 0 - 2 0 0 + 3 0 0 = 0 (7.12f)
Demand,0-2 m3/s
q. Head,
\1OO-OO
Q3 ]Q/'i*o6 Q
in q i
(7.5a), and (7.5c) . . . (7.5g), respectively. As one nodal head is known in-
stead of one nodal flow, Equation (7.5b) is replaced by Equation (7.12g).
In general, therefore, for each introduction of a known nodal head in place
of a known nodal flow, one nonlinear loop-head loss equation replaces a
linear node-flow continuity equation.
Using the linearized form, Equations (7.12e) . . . (7.12g) become re-
spectively,
The unknown nodal heads are taken as the basic unknown parameters in
formulating H equations. Sometimes, the unknown nodal flows are also
taken as the basic unknown parameters to form H-q equations. The pseudo
pipes and the pseudo loops are not necessary in formulating the H or H-q
equations.
General Formulation. By applying the node-flow continuity relationship
to the demand nodes labelled M + 1 to M + N, the H equations are
C = ,x = 1... X (7.13c)
Multisource Networks with Known Pipe Resistances 153
(7.14a)
/ connected I
j through x
F. = \ - * - ^ \ + F ^ M -0.4 = 0 (7.15c)
-0.1=0 <7.15d)
(
fj QC\ I In /QC t / \ 1/n / QC _ ZJ \ II n
In formulating the AQ equations, the pseudo pipes and pseudo loops are
introduced as is done in formulating the Q equations. It is also necessary
to initialize the pipe flow, i.e., assume the pipe flows for the first iteration.
The basic unknown parameters are AQC, c = I, II, . . . C, C +
1 . . . C + M - 1.
General Formulation. From the loop-head loss relationship for all the
real and pseudo loops
S /? o ,(,a + £ IAGC)- = 0, c = I, I I , . . . C, C + 1, . . . C + M - 1
(7.17)
in which the summations are taken as explained earlier for Equation (7.6).
Equation (7.17) is the same as Equation (7.6) except that it is now applied
to pseudo loops also. Because a known q value is replaced by a known H
Multisource Networks with Known Pipe Resistances 155
7.5.2.4 AH Equations
i connected to
j through x
, x, . + Aft) - (90
^3 —
40
- 0.1 = 0 (7.20d)
Equations (7.20a) . . . (7.20d) provide four independent nonlinear
equations for the four basic unknowns Aft . . . AH6. Note that these
equations are the same as Equations (7.9a) and (7.9c) . . . (7.9e) except
the change in node labels. As in the case of H equations, for each addi-
tional known H value, the number of AH equations decreases by one.
1 100.00 Qi 1 12.30
2 95.32 0.15 2 R2
3 94.78 0.31 3 48.60
4 94.30 0.30 4 24.60
5 H5 0.13 5 18.60
6 He 0.07 6 Re
7 7 49.60
— — — 8 58.10
— — — 9 14.30
X =J + C- 1 (7.21)
Further, let 7UH denote the number of nodes with unknown head, JVq
the number of nodes with unknown flow (inflow or outflow) and XVR
denote the number of pipes with unknown resistance. According to the
unknown-number rule for network solvability,
Note that when both, the head and theflow,are unknown at a node, the
node is included in both, Jvli and JLq.
The formulation of the different equations is explained and illustrated
through the network of Figure 7.6. The unknown parameters are qu R2,
Hs, H6, R6, q7 and H7. The values of the known parameters are given in
Table 7.1.
Note that the network contains two unknown resistances: R2 and R6.
Further, the heads and flows are known at nodes 2, 3 and 4, while they are
unknown at node 7.
The network has seven nodes and seven unknowns, and therefore the
unknown-number rule is satisfied. //, . . . H4 are known and q, and q7 are
unknown, and therefore the minimum-number rule is also satisfied. The
assignment rule is satisfied as the assignment of unknowns to nodes can be
done as follows: qi to node 1, q7 to 7, H7 to 2, R6 to 3, H6 to 4, H5 to 5,
and R2 to 6. Three other ways of assignment are also possible:
158 FORMULATION OF EQUATIONS
(1) g, to 1, q7 to 7, H7 to 2, R6 to 3, H6 to 4, R2 to 5, and // 5 to 6
(2) gi to 1, q7 to 7, /?6 to 2, // 5 to 3, H7 to 4, // 6 to 5, and R2 to 6
(3) qx to 1, ^7 to 7, fl6 to 2, Hs to 3, W7 to 4, /?2 to 5, and Hb to 6
General Formulation. All the pipe flows Qx, x = 1 . . . X and the un-
known pipe resistances, Ru, u = 1 . . . XVR are treated as the basic un-
known parameters in formulating the Q-R equations. Thus, the total
number of the unknowns is X + XUR. Applying the node-flow continuity
relationship to nodes with known nodal flows, J — JUq linear equations
can be obtained. Because the network has Cbasic loops, Cnonlinear loop-
head loss equations can be obtained. Further, as the network has J — 7UH
nodes with known heads, J — Jul — 1 head-loss equations for the
pseudo loops can be obtained. Thus, the total number of the available
equations is J — JUq + C + J — Ji:H — 1 which, using Equations
(7.21) and (7.22) is equal to X + XUR, the number of basic unknowns.
Thus, an adequate number of Q-R equations can be formulated to evaluate
the basic Q-R unknowns.
As is shown in the illustration that follows, it is possible to reduce the
X + Xy,R number of Q-R equations to X number of Q equations.
Illustration. In the network of Figure 7.6, X — 9 and XUR = 2. There-
fore, eleven Q-R equations must be available.
Node.l Pipe,5 2
IV
(b)
Figure 7.7 Loops and pseudo loops for network of Figure 7.6: (a) basic loops, (b) pseudo loops,
(c) overlapping loop, (d) replaced overlapping loop.
F2 = 0 4 + 0 6 - 0 7 - 0.31 = 0 (7.23b)
^ = 03 + 07 - 0« - 0.30 = 0 (7.23c)
^5 = 02 - 03 - 0.07 = 0 (7.23e)
The three basic loops I . . . Ill shown in Figure 7.7(a) yield, respec-
tively
Because the network has four nodes with known heads, three pseudo
loops can be selected. It is preferable to select them such that they do not
contain pipes 2 and 6 with unknown resistances. Selecting pseudo loops
IV . . . VI as shown in Figure 7.7(b), the loop-head loss relationship
yields, respectively
- 2 4 . 6 0 0 - 12.300 = 0 (7.231)
(7.24a)
95 32 - H-,'
71
14.30 -0.15=0 (7.24b)
94.78 - 94.30 \ 1 / n
- 0.31 = 0 (7.24c)
49.60
94.30 - HAUn
- 0.30 = 0 (7.24d)
58.10
R
162 FORMULATION OF EQUATIONS
'Hs - H(
F
* = iT^ + ^"^ ~ q7 =0 (7.24g)
H6 - 94.30 \ 1/n
- 0.20 = 0 (7.24h)
48.60
which now does not contain the unknown resistance R2 and contains only
the basic //unknowns. Equation (7.24h) can be obtained by merging nodes
5 and 6 and then writing the node-flow continuity equation for this merged
node 5-6. This merged node 5-6 has pipe 1 with inflow and pipes 4 and 3
with outflows at node 5-6. The total external flow at this merged node is
0.13 + 0.07 = 0.20 m3/s. This technique of merging nodes at the ends of
a pipe of unknown resistance to obtain a merged node is similar to the
technique of merging basic loops to obtain an overlapping loop. Both these
techniques eliminate the pipe of unknown resistance from the correspond-
ing equations.
Similarly, by combining Equations (7.24b) and (7.24c), i.e., by merging
nodes 2 and 3, we get
assumed Q and R values be: ,<2i = 0.7 rrrVs; iQ2 = 0.3 m 3 /s;
, g 3 = 0.23 m 3 /s; ,Q4 = 0.27 nWs; ,QS = 0.6 m°7s; , g 6 = 0.3 mVs;
,(?7 = 0.26 m 3 /s; ,(?„ = 0.19 m 3 /s; ,(?, = 0.15 m 3 /s; ^ = 25; and
i/?6 = 30. The AQ-AR equations for loops I . . . VI will be, respectively
- 100.00 = 0 (7.25e)
Thus, six independent AQ-AR equations are obtained for six unknowns:
A/?2, A/?6, AQm . . . AGvi.
To obtain AQ equations, since there are four AQ unknowns, four equa-
tions containing only AQ terms and not AR terms must be formed. Equa-
tions (7.25d), (7.25e), and (7.250 do not contain AR terms. Equations
(7.25a) and (7.25b) contain AR terms but they are used to adjust the
assumed R6 and R2 values, respectively. As Equation (7.25c) contains a AR
term, it is replaced by an equation for a loop that does not contain either
pipe 2 or pipe 6. Therefore, consider the overlapping loop as shown in
Figure (7.7c), or simply the pseudo loop shown in Figure (7.7d). Thus,
Equations (7.25a), (7.25b), (7.25c), and (7.25f) would now modify to, re-
spectively,
» - M - £ r "»]"•-0.3.-.
100.00 - (98.00 + AH5)V/n [(98.00 + AH5) - 94.781""
_
12.30 | I 24.60 J
- 0.13 = 0 (7.26e)
- 0.46 = 0 (7.26i)
Consider the network shown in Figure 7.8 with a pump provided in pipe
1 to lift water from the sump and supply it to the network. The network
also contains an elevated reservoir at node 2, connected to the network
through pipe 7.
Consider for the time being that reservoir 2 and pipe 7 are absent so that
the network is a single-source pumped network. As all the nodal demands
are known, the discharge supplied by the pump is known, and using the
head-discharge relationship, the head supplied by the pump is known. This
is then similar to the single-source networks described earlier.
The network of Figure 7.8 has two source nodes—a pump and a reser-
voir—and therefore pseudo loop, loop III, is introduced.
Inclusion of Pumps 169
2
Elevated
Reservoir,
H02
7.7.1.1 0 Equations
The network has seven pipes and therefore seven independent Q equa-
tions are necessary. Node-flow continuity equations at nodes 3 . . . 6, and
loop-head loss equations for real loops I and II provide six equations. The
seventh equation is available from the pseudo loop as follows:
7.7.1.2 H Equations
hp — H3 = (7.28)
from which
(7.29)
170 FORMULATION OF EQUATIONS
or
1. + H - H3
(7.33)
K
op
/"
i<a - qo3 = 0 (7.34)
Thus, hp need not be taken as the fifth unknown parameter. The un-
known parameters are H3 . . . H6 for which the three node-flow continuity
equations at 4, 5, and 6 and Equation (7.34) provide the required four H
equations.
Inclusion of Pumps 171
7.7.1.3 AQ Equations
The pipe flows are initially assumed so that they satisfy the node-flow
continuity relationships at nodes 3 . . . 6. Thus, the AQ equations for
loops I . . . Ill are, respectively,
7.7.1.4 AH Equations
The four H unknowns, H3 . . . H6, and the pump head hp are initially
assumed, and corrections AH3 . . . AH6 and Ahp are applied so that the
corrected values are iH3 + AH3 . . . !# 6 + AH6, and ihp + Ahp. These
values are then used in node-flow continuity equations at nodes 4, 5, and
6 in Equations (7.30) and (7.31). If Equation (5.7) is used as a head-
discharge relationship for the pump, hp need not be assumed as explained
earlier. The node-flow continuity equations at nodes 4 . . . 6 and 3 [Equa-
tion (7.34)] provide four AH equations for the four basic AH unknowns,
AH3 . . . AH6.
Booster pumps are provided within the networks and therefore they lie
on pipes forming real or pseudo loops. The formulation of different equa-
tions, however, is similar to that for supply pumps and therefore it is not
described separately.
172 FORMULATION OF EQUATIONS
We have seen in Chapter 5 that when check valves are provided in pipes
they permit flow in one direction. Different situations regarding the provi-
sion of check valves are (1) check valve in an external pipe, (2) check valve
in an internal pipe, and (3) check valves in several pipes. Formulation of
different types of equations for these situations is explained herein.
7.8.1.1 0 Equations
Figure 7.9 Network with check valve: (a) check valve in external pipe, (b) check valve in internal
pipe.
Inclusion of Check Kilves 173
negative, Q7 is made zero. This case would be as if pipe 7 does not exist,
and the network has six pipes. Therefore, Q7 is made zero in the node-flow
continuity equations for node 5 and 6, i.e., Equations (7.3d) and (7.3e), re-
spectively. Further, as pipe 7 does not exist, loop II also does not exist and
therefore the head loss equation for this loop, i.e., Equation (7.3g), is
dropped. Thus, for the six unknowns Qx . . . Q6, the available six Q equa-
tions are Equations (7.3a) . . . (7.3c), modified Equations (7.3d) and
(7.3e), and Equation (7.3f).
7.8.1.2 H Equations
7.8.1.3 AO Equations
The AQ equations are the same as Equations (7.7a) and (7.7b). However,
during any stage of the solution of the AQ equations if the application of
the correction AQn to the then assumed discharge in pipe 7 reverses the
flow direction, i.e., for the fh iteration if ,Q7 + ,AQn < 0, the correc-
tion tAQu is restricted so that ,Q7 + tAQu is zero. Thus, for the next
iteration, t+1£?7 = 0. For example, for the first iteration, as iQ7 = 0.05
m3/s from node 6 to node 5, if iAQn > 0.05 nrVs, in the clockwise direc-
tion, then ,AQ,i is restricted to 0.05 m3/s in the clockwise direction.
7.8.1.4 AW Equations
7.8.2.1 Q Equations
If the procedure described earlier for the check valve in an external pipe
is followed, there is no problem for the node-flow continuity equations.
However, as pipe 4 is common to both loops, the procedure would elimi-
nate both the loops, and the number of available equations would be five
for the six Q unknowns. This can be avoided by externalizing pipe 4, i.e.,
keeping it in only one loop. This can be achieved by considering two
loops: (1) loop I, an overlapping loop comprising pipes 1, 2, 5, 7, 6, and
3 and (2) loop II comprising pipes 1, 2, 4 and 3 or pipes 4 , 5 , 7, and 6. The
Q equations are then formulated as explained earlier for a check valve in
an external pipe.
7.8.2.2 H Equations
7.8.2.3 AO Equations
7.8.2.4 AH Equations
If check valves exist in two pipes, say pipes 7 and 4, or in pipes 5 and
4, the check valves can be externalized by considering the overlapping
loop and the basic loop I, and the procedure described earlier can be fol-
lowed.
2 3 3
(d)
d,
(f)
Figure 7.10 Networks with pressure-reducing valves: (a) PRV in external pipe, (b) considered
loops, (c) PRV in internal pipe, (d) considered loops, (e) PRVs in several pipes without isolating
a portion, (f) considered loops, (g) PRVs in several pipes isolating a portion, (h) considered
loops.
From this it is clear that a PRV can be replaced by (1) a source node d
Inclusion of Pressure Reducing Valves Yll
with known HGL equal //se, and inflow Q7 and (2) a sink node u with
unknown head and outflow Q7.
7.9.1.1 0 Equations
The node-flow continuity equations for nodes 2 . . . 6 and the head loss
equation for loop I are given by Equations (7.3a) . . . (7.3f), respectively.
However, the loop-head loss equation, Equation (7.3g), is not applicable
because of condition 1. The PRV is replaced by source and sink nodes, and
a pseudo loop is formed as shown in Figure 7.10(b). Therefore, for this
pseudo loop, loop II
and
F, - ( « ^ * ) - ( « - = ; * ) - 0..-0 (7.37b,
The AQ equations can be written for basic loop I and pseudo loop II
shown in Figure 7.10(b). When the correction AQ, is obtained for loop I
178 FORMULATION OF EQUATIONS
Equations (7.7a) and (7.38) provide two AQ equations for two correc-
tions, AQi and AQn. As stated earlier AQ will be applied to pipes
1 . . . 4 of loop I and AQt to" pipes 4 . . . 7 of real loop II.
7.9.1.4 AH Equations
and
- 0.1 = 0 (7.39b)
7.9.2.1 0 Equations
F6 = 2 0 0 + 4 0 0 - 5 0 0 - 2 0 0 - 1 0 0 - 4 0 0 = 0 (7.40a)
and
7.9.2.2 H Equations
The H equations for the flow continuity at nodes 2, 5, and 6 are the same
as Equations (7.5a), (7.5d), and (7.5e), respectively. The H equations for
nodes 3 and 4 modify to, respectively,
1
~\ 10 / \ Ro.d4 ) [ 20 j
+ 0.1 = 0 (7.41a)
180 FORMULATION OF EQUATIONS
and
100 - HAV" + / HM - H4
3
~ ' 20 / • \ R^4 ) \ 40
- 0.2 = 0 (7.41b)
7.9.2.3 AQ Equations
The AQ equations for the overlapping and pseudo loops are, respec-
tively,
and
7.9.2.4. AH Equations
and
Problems 181
7.10 PROBLEMS
2 0-08m% - 0-06
as demand nodes with demands as shown in the figure. The length, diam-
eter, and HW coefficient of the pipes are as follows: Pipe 1—300 m, 250
mm, 120; pipe 2—250 m, 300 mm, 110; pipe 3—300 m, 200 mm, 120;
pipe 4—400 m, 100 mm, 100; pipe 5—200 m, 200 mm, 90; and pipe
6—300 m, 150 mm, 100; respectively.
Problem 7.2: The two-source network shown in Figure: Problem 7.2
has nodes 1 and 2 as source nodes with HGL values 110.00 m and 100.00
7.8 has nodes 1 and 2 as source nodes with HGL values of 105.00 m and
100.00 m, respectively. Nodes 3 . . . 6 are demand nodes with demands as
shown in the figure. The R values for the pipes are shown underlined in the
figure. The head loss in pipes is given by h = RQl8&2 Qi in meters, Q in
cubic meters per second). A PRV (K = 10) is fitted in pipe 6 (diame-
ter = 400 mm) just downstream of node 5 as shown in the figure, with a
set pressure head of 90.00 m.
Hardy Cross Method
8.1 INTRODUCTION
A iQi
i Roi
1
oc o
E en
q
e —,
k Rok
,Qk B
(1) Only one equation is considered at a time from the available set of
AQ equations.
(2) The effect of adjacent loops is neglected and therefore each AQ equa-
tion contains only one AQ term as unknown. For example, in Equa-
tion (7.7a) of Chapter 7, the correction AQ, is retained and AQ,, is
discarded. Similarly, in Equation (7.7b) the correction AQn is re-
tained and AQ, is discarded. Thus, except the correction for the
loop under consideration, corrections for other adjoining loops are
assumed to be zero.
(3) Each term of the modified AQ equation is expanded in a Taylor's
series with only the first-power AQ terms retained, and higher-power
AQ terms are neglected.
Now consider a single loop of a distribution network as shown in Figure
8.1 (assumption 1). Let it consist of four pipes—i,j, k, and m, with known
resistances Rol, Roj, Rok, and Rom, respectively. Let the assumed pipe dis-
charges satisfying the node-flow continuity relationships be ,£?., iQj, i2*>
and iQm in pipes i, j , k, and m, respectively, in the direction as shown in
the figure. Let IAQ be the clockwise correction along the loop as shown in
the figure. The AQ equation for this loop, after disregarding the effect of
adjacent loops (assumption 2) is
- iA0- = 0 (8.1)
Expanding in a Taylor's series and neglecting higher-power terms in
(assumption 3), we get
1 U
~ Rorn.Q;-1 + ^ - n ' , 0 - ' + i^-n-ifiT 1 + iC-n-ifiE" 1
(8.3)
Note that in Equation (8.3) the numerator gives the algebraic sum of the
head losses in pipes along the loop, the head loss being positive for pipes
/ and j and negative for pipes k and m. The denominator, however, is the
sum of all positive terms because when the direction of flow in a pipe is
negative, the correction is also negative. Therefore, Equation (8.3) can be
written as
E R- • i<2"
AQ = - ^ r 7 (8-4)
n • R^ • iQ"'1
• ,Ql
, c = I, H . . . C (8.5)
» • * « • 'Q"'1
As each term of the numerator represents the head loss in a pipe, Equa-
tion (8.5) can also be written in terms of the head losses in pipes as
EA
, c = I, II . . . C (8.6)
n • ,hx
but considers each loop separately for the calculation of the corrections as
in the Hardy Cross method. However, the simultaneous approach which is
more popular than the sequential approach continues to be known as the
Hardy Cross method and is also so explained in textbooks of fluid
mechanics, hydraulics, and water supply engineering, for example, refer-
ences [9,82]. Herein also the simultaneous approach is considered and
termed the Hardy Cross method. Liu has also shown that convergence is
faster, and therefore, fewer iterations are required in the simultaneous ap-
proach than in the sequential approach.
The approach based on AQ equations tries to make the algebraic sum of
the head losses along loops equal to zero. In other words, it tries to balance
the nodal heads obtained from different paths. For example, when the loop
of Figure 8.1 is balanced, the head obtained at node B from that at node
A by considering the path consisting of pipes / and j must be the same as
that obtained by considering the path consisting of pipes m and k. There-
fore, this approach based on AQ equations is known as the method of bal-
ancing heads.
The Hardy Cross method of balancing heads is based on the three
assumptions cited earlier. The third assumption that neglects the higher-
power AQ terms in Taylor's expansion is tested in Example 8.1. The net-
work consists of a single loop so that the effect of the first two assumptions
is absent.
Example 8.1
node, \0-4
Hol = 100.00 m, Rol = 8, Ro2 = 10, Ro3 = 20, Ro4 = 15, qo2 = 0.2
m3/s, qo3 = 0.4 m3/s, and qo4 = 0.3 nrVs. Determine the pipe discharges
and the heads at demand nodes if the head loss in a pipe is given by
h = RQ2.
Method of Balancing Heads 189
Solution
The loop corrections are initially calculated using Equation (8.4) for all
loops of a network and are then simultaneously applied to complete one
iteration. When a pipe is common to more than one loop, the corrections
of all these loops are applied to the discharge in the pipe to obtain revised
discharge for the next iteration. The iterative procedure is continued until
satisfactory convergence is reached.
Example 8.2
Pipe,
X ,0, ,AO
Iteration 1
Iteration 2
4.4575
2 10 0.55 3.0250 11.00
32.10
3 20 0.35 2.4500 14.00
= -0.1389
4 15 -0.05 a -0.0375 1.50
L 4.4575 32.10
Iteration 3
Iteration 4
q o2 =0-8nrr 3 /s
Hoi=100-00m
Solution
Let the assumed direction and the magnitude of discharges in the pipes
be as shown in the figure. The analysis is carried out as shown in Table:
Example 8.2 in which the first two iterations are given. Convergence is
reached in five iterations and the solution is as follows: Qt = 1.6745 m3/s;
Q2 = 1.3255 m 3 /s; Q3 = -0.2216 nrVs; Q4 = 0.6529 nrVs; Qs =
0.6471 nrVs; q, = 3.0000 m 3 /s; H2 = 87.010 m; H3 = 84.285 m; and
H4 = 86.518 m.
Note that as pipe 3 is common to both loops, loops I and II, both the cor-
rections AQi and AQU are applied to the discharge in it. Convergence is
slightly delayed because the effect of loop II on loop I and that of loop I
on loop II are neglected. In general, the number of iterations increases as
the size of the network increases because of the interaction among the
loops.
Because the network has multiple sources, pseudo loops, one less than
the number of sources, are introduced. The loop corrections are calculated
for the real loops using Equation (8.5). Applying the procedure used in
deriving Equation (8.5) to the pseudo loop consisting of starting node S,
end node E, and real and pseudo pipes as shown in Figure 8.2, the correc-
tion term for the fh iteration of this pseudo loop is given by
-Q" + HF - Hs
xeP
(8.7)
,Qr
TABLE: EXAMPLE 8.2.
(Clockwise Positive, Anticlockwise Negative.)
Loop I Loop II
Iteration 1
Iteration 2
Eo- 1
pseudo loop, p
in which HE, Hs = heads at end node and starting node, respectively, and
subscript p denotes a pseudo loop. The term HE — Hs can be looked upon
as the head loss in the pseudo pipe and therefore Equation (8.5) can be
used in place of Equation (8.7). The head loss in the pseudo pipe remains
fixed throughout the iterative procedure. Note that the denominator con-
tains the terms for the real pipes only and no term for the pseudo pipe.
Example 8.3
Solution
Assume the first-trial pipe discharges as ,(?, = 0.35 m 3 /s, tQ2 = 0.45
m /s, ,03 = 0.05 m 3 /s, ,Q 4 = 0.05 m 3 /s, XQS = 0.30 m 3 /s, tQ6 = 0.20
3
m 3 /s, and tQ7 = 0.10 nvVs along the directions shown in the figure. The
first two iterations are shown in Table: Example 8.3. Note that pipe 4 is
TABLE: EXAMPLE 8.3.
(Clockwise Positive, Anticlockwise Negative.)
Iteration 1
Iteration 2
common to loops I and II, while pipes 1 and 3 are common to loop I and
pseudo loop III. HE - Hs = 95.00 - 100.00 = -5.00 m. After several
iterations the solution is <2i = 0.3246 m3/s, Q2 = 0.4715 m3/s, Q3 =
0.0246 m3/s, Q4 = -0.0219 m3/s, Q5 = 0.2496 m'/s, Q6 = 0.2503 nVVs,
Q7 = 0.1504 m 3 /s, H3 = 95.010 m, HA = 95.025 m, H5 = 91.956 m,
H6 = 93.459 m, q, = 0.7961 m3/s, and q2 = 0.2039 m3/s.
This method of balancing heads for the analysis of networks with un-
known pipe resistances solves either the AQ-AR equations or the reduced
AQ equations. (The formulation of both these types of equations is ex-
plained in Chapter 7.) When the approach is based on AQ-AR equations,
the assumed unknown pipe resistances are also simultaneously adjusted
along with the assumed pipe discharges, and the method is implicit. When
the method is based on the reduced AQ equations, only the pipe discharges
are assumed and adjusted. After obtaining satisfactory convergence, the
unknown pipe resistances are calculated in the end. This approach is ex-
plicit.
r-,-4
H OA
'0,
RL
IQU
pseudo loop
AQ - - ^ x = i,j, k, m
L, I n- Rox- .QT 1 | + | n • ,ft, • .QT 1 |
(8.8)
+ RoMQrr, + ,AG)"
Equation (8.9) is used for adjusting the assumed value of pipe resis-
tance, i.e., iRu. Therefore correction ,A/?U is retained and correction ,A(?
is disregarded. Therefore, Equation (8.9) becomes
It* • iQl + Ron, • ,(£, + (,/?„ + tARH) • XQL + HoB - HoA = 0 (8.10)
in which the subscript x also includes the pseudo pipe, the head loss in
which is HoB — HoA and remains fixed throughout. Thus, the numerator
represents the algebraic sum of the head losses in all pipes (including the
Method of Balancing Heads 197
pseudo pipe) of the pseudo loop. Note that the numerator may be positive
or negative and so also the denominator depending upon the direction of
flow in pipe u.
Generalizing Equation (8.11) for all pseudo loops used to adjust un-
known pipe resistances, Ru, u = 1 . . . U, the correction for the t'h itera-
tion is given by
Example 8.4
The head loss in a pipe of the network shown in Figure: Example 8.4(a)
is given by h = RQ2. Analyze the network by head-balancing procedure
using (1) implicit and (2) explicit approaches.
Solution
The unknown parameters are qx, H3, H4, R4 and Hs. All the network
solvability rules are checked and found to be satisfied.
1. Implicit Approach. For this network, X = 6 and / UH = 3. Thus,
% ~ -AJH = 6 - 3, i.e., three independent loops are available. Out of
these three loops, X lR , i.e., one loop will be used to adjust the assumed
R4 value and the remaining two loops to adjust the pipe discharges.
Herein, the two basic loops, loops I and II in Figure: Example 8.4(b) are
198 HARDY CROSS METHOD
0-5 q
node,
100-00m'
Loop, I 3 3 II 2
(b)
Loop,I
used to adjust the pipe discharges and the third loop, loop III is used to ad-
just the unknown resistance R4. Let the assumed pipe discharges be ,Q, =
0.7 m3/s, ,Q 2 = 0.5 m3/s, , & = 0.2 nrVs, ,Q4 = 0.4 nrVs, tQ5 = 0.3
m3/s, and ,Q6 = 0.3 m3/s. Let the assumed pipe resistance be ,V?4 = 20.
The first two iterations of the iterative procedure are given in Table 1:
Example 8.4. Continuing the procedure for several iterations satisfactory
convergence is reached. The solution is: Qt = 0.6851 m 3 /s; Q2 = 0.5149
m 3 /s; Q3 = 0.2003 m 3 /s; Q4 = 0.3854 m 3 /s; Q5 = 0.3146 m 3 /s; Q6 =
0.2854 m 3 /s; R4 = 14.85; q, = 1.2 mVs; H3 = 97.349 m 3 /s; H4 =
96.245 m; and // 5 = 94.040 m.
2. Explicit Approach. In the explicit approach, two loops, loops I and II,
shown in Figure: Example 8.4(c) are considered for AQ corrections, and
the calculations are shown in Table 2: Example 8.4. As nodal flow at node
2 is known and fixed, when the correction AQ]U is applied to pipes 2 and
6 of loop II, the node-flow continuity equation at node 2 is violated.
Therefore, to maintain the node-flow continuity relationship at node 2 and
also at nodes 5, 4, and 1, the correction AQn is extended to pipes 5, 4,
and 1. Therefore, for iteration 2, 2Qi = 0.7 - 0.010 (correction for loop
TABLE 1: EXAMPLE 8.4.
(Implicit Approach.)
Loop ! Loop II Loop
2 2
Pipe 0 RO 2RQ Pipe RO 2RQ Pipe Q RQ2 Q2
Iteration 1
Loop I Loop II
Iteration 1
0.32 0.22
= -0.010 1A = -0.015
32.2 14.8
Iteration 2
- 0.324 -0.106
2 AQ. = - 0.010 2AQ,, = - = 0.007
32.85 15.06
I) — 0.015 (correction for loop II) = 0.675 m3/s. Thus, loop II compris-
ing real pipes 2 and 6 is used to calculate AQn but the correction is ap-
plied to the overlapping loop comprising pipes 2, 6, 5, 4 and 1. The itera-
tive procedure is continued until satisfactory convergence is reached. In
the end, loop III is used to calculate the value of ft, and the final solution
as given earlier is obtained.
Example 8.5
1. Implicit Approach. The selected three loops are shown in Figure: Ex-
ample 8.5(a). Loops I and II are the same as loops I and II of Figure: Ex-
ample 8.4(b), respectively. Further, as in Example 8.4, these two loops are
Method of Balancing Heads 201
Loop I
202 HARDY CROSS METHOD
According to this procedure the AQ equations for the given and assumed
values are
and
Equations (3) and (4) are solved iteratively, and the first two iterations
are given in Table: Example 8.5. Note that the correction AQ, is applied
to pipes 1, 4, 5, 6, and 2, while AQn is applied to pipes 1, 2, and 3. Fur-
ther, AQn is obtained from loop III in which the RQ2 values are taken for
pipes 6 and 2 while the IRQ value is taken for pipe 2 only as it is the pipe
common to loops II and in.
8.2.5 NETWORKS WITH PUMPS
The AQ equations for the real or pseudo loops contain hp and therefore
an expression of Qp as explained in Chapter 7. The balancing of heads for
such loops is carried out as shown in the following example.
Example 8.6
Analyze the network shown in Figure: Example 8.6 using the Darcy-
Weisbach formula for the head loss in pipes and the Swamee-Jain explicit
Method of Balancing Heads 203
Iteration 1
0.696 0.22
iAQ, = - = -0.015 iAON = -0.022
47.08 10
Iteration 2
-0.251 -0.234
2AQ,, = - 0.022
47.24 10.74
relationship for the friction coefficient. The pipe details are as given
in Table 1: Example 8.6. Three points on the head-discharge curve for
the pump have the following discharge to head characteristics: 0.1
nvVs - 11.89 m; 0.2 nWs - 11.01 m; and 0.3 m3/s - 9.37 m, respec-
tively.
Sump E.S.R.
W.L? 80-00 m W.LrlOO-OOm
204 HARDY CROSS METHOD
Solution
Using the pump data, the head-discharge equation for the pump is ob-
tained as
and
Gp = Qp - 0.0353 (3)
and
Loop
c Pipe x iAO
4 - 1.55 154.6
II 5 0.06 0.0215 835 -3.00 100.2 -0.0031
6 0.04 0.0241 3933 6.29 314.7
r 1.74 569.5
The formulation of the AQ equations for networks with check valves and
pressure reducing valves has been explained in Chapter 7. The loops are
selected such that these valves are externalized.
For check valves the correction is restricted so that the discharge in the
pipe containing a check valve does not change direction.
When a PRV is present in a network, it is initially presumed that it is in
operative condition. An appropriate pseudo loop is formed for calculating
the AQ correction and is applied to the loop such that flow continuity is
maintained. The procedure is explained in Example 8.7.
Example 8.7
Analyze the network shown in Figure: Example 8.7(a). The R values for
5 0-4
100-00 m,
a -1
Loop I
the pipes are as given in the figure. Pipe 7 having diameter of 350 mm con-
tains a PRV, just downstream of node 3 with a set pressure head of 80.00
m. The head loss in pipes is given by h = RQ2.
Solution
Loop I Loop I
Assume that the nodal demands are 0.2, 1.5, 0.2, 0.2, and 0.1 m3/s at
nodes 2, 3, 4, 5, and 6, respectively. Even though initially it is assumed
that the PRV is operative and the flow direction in pipe 7 is from node 3
to node 6, during the iterative procedure, the correction changes the flow
direction in pipe 7. Herein the PRV acts as a check valve and prevents flow
from node 6 to node 3. The discharge in pipe 7 is zero, and the subsequent
correction in loop II (or loop III) is clockwise, and therefore it is restricted
to zero. Thus, the correction is applicable to loop I only. The final solution
is Ql = 2.2 nrVs, Q2 = 0.7595 m3/s, Q3 = 1.2405 m3/s, Q4 = -0.2595
m3/s, Q5 = 0.3 m3/s, Q6 = 0.1 nrVs, and Q7 = 0. H, = 84.16 m and
H6 = 85.93 m.
Example 8.8
Assume that in the network of Figure: Example 8.7, a PRV is also pro-
vided in pipe 5 (500 mm diameter) just downstream of node 4, as shown
in Figure: Example 8.8(a). The pressure-head settings of PRV-1 and PRV-2
are 60.00 m and 65.00 m, respectively. Analyze the network.
100-OOm 1
r
Solution
It is presumed that both the PRVs are operative and are replaced by
source and sink nodes as shown in Figure: Example 8.8(b). The PRVs have
separated nodes 5 and 6 from the main network. This separated network,
loop II, is analyzed and the pipe discharges are Q<, = 0.7847 m3/s,
Method of Balancing Flows 209
The method of balancing flows |101| solves the AH equations on the fol-
lowing assumptions:
(1) Only one equation is considered at a time from the available set of
AH equations.
(2) The effect of the adjacent nodes is neglected and therefore each AH
equation contains only one AH term as unknown. Thus, except the
correction for the node under consideration, corrections for other
connected nodes are assumed to be zero.
(3) Each term of the modified AH equation is expanded in a Taylor's
series and only the first-power AH terms are retained and the higher-
power AH terms are neglected.
[ Roc,
UHj + AHj) -
_ _
ft* J I ft- J (8.13)
tjoj —
or
[ 1HJ-1Hk lAHJV" [ , ]
(8.14)
, ; = M + I,. . . ,M + N (8.16)
During the iterative procedure, when the head //, and that at a connected
node are equal, Qlnh becomes infinite for the pipe connecting them. In
this case, such a pipe is omitted for this iteration.
Method of Balancing Flows 211
As this approach tries to make the algebraic sum of the flows at a node
equal to zero, i.e., it tries to balance the total inflow and total outflow, it
is known as the method of balancing flows.
Example 8.9
Solution
Assumed
Head, h,, Q, Correction,
Of
Node m Pipe m m3/s n • hij m
"External outflow.
term corresponding to the unknown assigned to nodey, we can get, for the
general t'h iteration,
p.
,AHj = zz when unknown //, is assigned to node j
(8.17)
i connected
toy
TJ TJ
in which ,Fj = residue of inflow over outflow at node j , for the t'h itera-
tion.
Equation (8.17) is the same as Equation (8.16) except that it includes a
term for the flow in a pipe of unknown resistance, assumed as f/?,v for the
f" iteration. Equation (8.18) gives the additive correction to the assumed
Hi value, for the fh iteration, while Equation (8.19) directly gives the Rtj
value for the (t + 1)" iteration. Equations (8.17) . . . (8.19) are used in the
iterative procedure to determine the unknown H and R parameters. Care
should be taken, during the iterative procedure, to use proper signs for the
various terms.
Figure 8.5 Part of distribution network having pipe with unknown resistance.
214 HARDY CROSS METHOD
(7.26i) provide three independent AH equations for the three basic un-
knowns, AH6, A// 5 , and AH7, respectively. Therefore, AHe is obtained at
node 4 using Equation (8.18); A// 5 is obtained at merged node 5-6 using
Equation (8.17); and AH7 is obtained at merged node 2-3 using Equation
(8.18).
Example 8.10
Solution
1. Implicit Approach. The //and R unknowns in the network are H3, H4,
R4, and Hs. To satisfy the assignment rule, H3 is assigned to node 3, H4 to
node 4, R4 to node 5, and H5 to node 2. Therefore, AH3 will be obtained
at node 3, A// 4 at node 4, AR4 at node 5, and AHS at node 2. Considering
the assumed flow directions in pipes, shown in Figure: Example 8.4, let
the first-trial values be ,// 3 = 97.00 m, ,// 4 = 96.00 m, ,// 5 = 95.00 m,
and ,/?4 = 20. The first iteration of the iterative procedure is given in
Table 1: Example 8.10. The second-trial values obtained at the end of first
iteration are 2 // 3 = 97.30 m, 2H4 - 96.59 m, 2HS = 95.60 m, and
2R4 = 5.51. Continuing the procedure for several iterations, the solution
obtained in Example 8.4 is reached. The same solution is obtained if the
alternate assignment of H3 to node 2, H4 to node 3, R4 to node 4, and Hs
to node 5 is made.
2. Explicit Approach. Pipe 4 with unknown resistance R4 has nodes 4
and 5 as end nodes. Therefore these two nodes are merged and the total
outflow is 1.2 mVs. For the flow directions assumed in the figure, Pipes 1,
3, and 5 are the supply pipes to this merged node 4-5. The unknown
assignments are H3 to node 3, H4 to merged node 4-5, and Hs to node 2.
The first iteration of the iterative procedure is given in Table 2: Example
8.10. The first-trial values are the same as those assumed for the implicit
approach. The second-trial values obtained at the end of first iteration are
2H3 = 97.30 m, 2 // 4 = 95.88 m, and 2 // 5 = 95.60 m. After several itera-
tions, the solution given earlier is obtained.
Assumed
Value of
Assigned Assigned Connected O
Node Unknown Unknown Pipe m m3/s n-hH Parameter
CM CM
CO
2 H5 95.00 m 8 1.00 -0.3536 AHi - 0.602 m
3.00 -0.2739 0.0457
LO
40
0.6000e
E -0.0275
a
External outflow.
B
Using Equation (8.17).
c
Firsttrial assumed value.
d
Using Equat.on (8.19).
e
External inflow.
'Using Equation (8.18).
TABLE 2: EXAMPLE 8.10.
Assumed
Value of
Assigned Connected
Node Assigned Unknown, Pipe QH AW,,
i Unknown m ij m m3/s n-h,. m
(2) Carry out one iteration and obtain the pipe discharges and the AH
corrections
(3) Taking the discharge in the pipe in which the pump is situated as the
discharge delivered by the pump, find the head delivered by the
pump
(4) Update the unknown nodal heads
(5) Using the updated nodal heads (step 4) and the head delivered by the
pump (step 3), carry out the next iteration (step 2)
(6) Continue the iterative procedure until satisfactory convergence is
reached
Example 8.11
Solution
Let the first-trial value for the head delivered by the pump, ,/?p be
10.75 m, a value corresponding to XQP — §.Tl rrvVs, assumed in Example
8.6. Considering the assumed flow directions in pipes (Figure: Example
8.6), let ^ 3 = 85.00 m; ,// 4 = 83.00 m; ,// 5 = 80.00 m; and ,// 6 =
90.00 m. The coefficient of friction is taken as 0.02 for all pipes, and the
corresponding R values are obtained.
The first iteration is shown in Table: Example 8.11: hi-3 = tHi +
lhp - tH3 = 80.00 + 10.75 - 85.00 = 5.75 m. Since <2,-3 obtained
after the first iteration is 0.1679 rrvVs, 2hp = 12 + 2.7 X 0.1679 -
38.2 X 0.16792 = 11.38 m. Using the obtained pipe discharges, the/and
R values of the pipes are updated, the heads for nodes 3 . . . 6 are cor-
rected, and second iteration is carried out.
After several iterations, the solution given in Example 8.6, with
H3 = 86.49 m; HA = 80.71 m; H, = 76.73 m; and H6 = 89.20 m is ob-
tained.
As stated in Example 8.6, if the head-discharge relationship given by
Equation (5.7) is used for the pump, the pump can be made an integral part
of the pipe in which it is situated. The sump level is increased by Hp and
the pipe resistance constant by Rp.
In the case of a check valve, as long as the flow in the pipe is in the
permissible direction, additional resistance offered by the check valve is
considered for the pipe in which it is present, and the normal iterative pro-
218 HARDY CROSS METHOD
Assumed
Node Value of Pipe On AH,,
/ Hr m ij m m3/s n-hH m
3 85.00 1-3 204 5.75 0.1679 0.0146 1.32
3-4 2066 2.00 -0.0311 0.0078
3-5 8707 5.00 -0.0240 0.0024
-0.0800 a -
E 0.0328 0.0248
4 83.00 3-4 _ — 0.0311 0.0078 -1.23
4-5 3265 3.00 -0.0303 0.0051
6-4 775 7.00 0.0950 0.0068
-0.1200 3 —
-0.0242 0.0197
5 80.00 3-5 — 0.0240 0.0024 -1.01
4-5 _ _ 0.0303 0.0051
6-5 3265 10.00 0.0553 0.0028
-0.1200 3 —
s -0.0104 0.0103
6 90.00 2-6 157 10.00 0.2524 0.0126 -0.81
6-4 — — - 0.0950 0.0068
6-5 - _ - 0.0553 0.0028
-0.1200a —
-0.0179 0.0222
a
External nodal flows—inflow positive, outf ow negative.
cedure is followed. However, if the flow direction reverses, the pipe con-
taining the check valve is excluded from the next iteration.
In the case of a pressure reducing valve, it is initially assumed that it is
in the operative mode and is therefore replaced by source and sink nodes.
The discharge in the pipe downstream of the PRV is determined, and this
discharge is taken as additional outflow at the upstream node of the pipe.
The iterative procedure is carried out and after adjusting the nodal heads,
a check is made at the end of each iteration to verify that the PRV has re-
mained in the operative mode. If the PRV becomes inoperative, the origi-
nal pipe is introduced in the iterative procedure but its resistance is in-
creased to include the additional resistance offered by the PRV. On the
other hand, after adjusting the nodal heads if it is observed that the head
at the downstream end of the pipe is more than that at the upstream end,
the PRV functions as a check valve and therefore the pipe is excluded from
the next iteration.
Thus, for check valves and PRVs it is necessary to ascertain their mode
of operation before the next iteration is carried out.
Method of Balancing Flows 219
Example 8.12
Example 8.13
Solution
As explained in Example 8.8, the nodes 5 and 6 are separated from the
main network. The analysis of the separated network gives Hs = 59.46
m, H6 = 57.68 m, Q4_5 = 0.7847 mVs, and ft-6 = 0.2153 m3/s. The
main network is now analyzed with outflows of 0.2 + 0.2153, i.e., 0.4153
nvVs at node 3 and 0.5 + 0.7847, i.e., 1.2847 rrrVs at node 4. The solution
is H2 = 92.00 m, H3 = 89.93 m, and H4 = 86.23 m.
220 HARDY CROSS METHOD
The usual Hardy Cross method, described so far, is based on three as-
sumptions: (1) considering one loop (node) correction equation at a time,
(2) neglecting the effect of adjacent loops (nodes), and (3) neglecting the
higher-power correction terms in the Taylor's series. As seen in Example
8.1, the third assumption doesn't have much effect on convergence. The
method is now modified herein so that the second assumption of neglect-
ing the effect of adjacent loops (nodes) is not made, i.e., the method takes
into consideration the effect of the adjacent loops (nodes) in deriving the
correction terms.
E *~ • >Q"
x
n • Rox • tQ7
n - R o x - ,QTX
,AQaxm (8.21)
n • Rox • ,QT
Modified Hardy Cross Method 221
iAQ 0 i m
i iQi
E Roi 'o1
iAQ arnrT1 ,-j iAQ Q jm
E iQk
'k ' Rok
iAQ o km
The first term on the right-hand side of Equation (8.21) represents the
usual correction tAQc given by Equation (8.5). The second term is the
summation of the influences of the adjacent-loop corrections. The term
within parentheses is termed the adjacent loop influence factor and is
denoted by Fax. Thus, Fai is a factor representing the influence of an adja-
cent loop connected through pipe x. Equation (8.21) can, therefore, be
written as
To evaluate the value of ,AQcm, the values of ,AQaxm are necessary. How-
ever, these values are also unknown and depend upon other loops. There-
fore, it would be necessary to solve simultaneously all such loop-flow cor-
rection equations. However, because only one loop is considered at a time
(assumption 1), ,AQaxm is replaced by ,AQax, the corrections obtained by
the usual method, and only half the effect is considered to compensate for
the replacement of tAQaxm by tAQax. Thus, Equation (8.22) modifies to
(8.23)
Example 8.14
Solution
= 0.024 + ^ | T T g 9 (-0.051)
or
r.ff.-.ff,-
R,.,,
-+
,AHim = ,AHin
Q_
•on neci ed
t" / (8.26)
The first term on the right-hand side of Equation (8.26) represents the
usual correction ,AHj given by Equation (8.16). The second term is the
summation of the influences of the adjacent-node corrections. The term
within parentheses is termed the adjacent node influence factor and is
denoted by Fai. Thus, Fai is a factor representing the influence of the adja-
cent node /. Equation (8.26) can now be written as
Herein also, to evaluate the value of ,AHjm, the values of ,A//,m are nec-
essary. However, these values are also unknown and depend upon other
nodes. Therefore, it would be necessary to solve simultaneously all such
nodal-head correction equations. However, since only one node is con-
sidered at a time, ,AH,m is replaced by ,Aff,-, and only half the effect is con-
sidered. Equation (8.27) now modifies to
Equation (8.28) gives the modified AH correction term, and because the
effect of adjacent nodes is considered, it converges faster than the usual
method. However, the computational work per iteration has increased.
Example 8.15
Solution
= -1.54 m
= 1.42 m
1 [0.0536 .. „ . 0.0653
= -°-68 + T 0J189 ( 1 - 5 2 ) + aTT89 ( 0 - 4 7 )
= - 0.21 m
AH
,A// 6m
- n 47 , [ [° 0 7 7 9 (0)
- 0.47 + 2 [ 0 J 4 3 2
,m+, °- 0653
= 0.31 m
flows) and also the pipe resistances if they are unknown. The computa-
tional work involved depends upon the initial values of these parameters.
If the initial values are close to the final values, the computational work
and therefore the number of iterations is less. Therefore, a good selection
of the initial values is desirable.
For better initial guess of pipe flows, Conklin [105] suggested the use of
a distributional factor for each pipe as a function of its diameter and
length. Williams [ 106] used distributional factors for initializing pipe flows
in his studies on convergence problems in network analysis. Two other ap-
proaches based on (1) a minimum spanning tree and (2) a minimum path
length tree can also be used. Using pipe resistance constants, these ap-
proaches convert the looped network into a branching configuration, i.e.,
a tree. (If the DW head ioss formula is used, the coefficient of friction
equal to 0.02 can be used for all pipes, and pipe resistance constants are
evaluated.) Pipes not on the tree are assigned zero flow and therefore the
discharges in pipes on the tree are obtained as described in Section 7.4.1
of Chapter 7.
250 j
a—
10 120 ' 200
o o CD
00
C\J
40 6 160 7 250
200 ,
140 '
s CD
CD
80 160 ( 200
10 11 12
(a)
10 130 33U
o- 1 5* '
2 9
60 460
'6 7
5 3*
oo 6
7 JOj 1 8, . 1 \t —•
10 11 12 80 150 310 510
(b) (c)
Figure 8.7 Selection of initial pipe flows by conversion of looped network to Iree: (a) using mini-
mum spanning tree approach; (b) using minimum path length tree approach.
pipes show the order in which the pipes are added. The same minimum
spanning tree is obtained if any other node is selected as a starting node.
Considering the pipe resistances instead of the pipe lengths, the lengths
of the different paths connecting a demand node to a source node are ob-
tained, and the path with the minimum path length is selected. This con-
cept is used in obtaining a tree. The algorithm, explained with the help of
the network in Figure 8.7(a) is as follows:
Consider all nodes directly connected to the source node (nodes 2 and
5), obtain their path lengths (paths 1-2 and 1-5 with path lengths of 10, and
20, respectively), and select the node with the least path length (node 2
with path length 10). The two nodes and the connecting pipe form a partial
spanning tree (nodes 1 and 2 and pipe 1-2). Next consider all nodes di-
rectly connected to the nodes on the partial spanning tree (nodes 5, 6, and
3), obtain their path lengths (20, 70, and 130 for nodes 5, 6, and 3, respec-
tively), and select the node with the least path length (node 5 with path
length 20). Add it to the partial spanning tree to obtain a new spanning tree
and continue the procedure. When all nodes are connected to the partial
spanning tree, the resulting network is the required tree. The tree obtained
for the network of Figure 8.7(a) is shown in Figure 8.7(c). The numbers
near the nodes show path lengths, and those near the pipes show the order
Convergence Problems 227
in which the pipes are connected. The minimum path length tree approach
can also be used for multiple source networks [86] in which a separate tree
is obtained for each source node.
Both these methods have some merits and limitations [107,109] but in
general provide good starting values.
The initial value for the unknown resistance of a pipe can be taken as the
average value of the resistances of the pipes connected to the pipe under
consideration at its end nodes. For example, in the network of Figure
8.7(a), the initial value of the resistance constant of pipe 2-3, if unknown,
can be taken as 87.5, the average of 10, 60. 80, and 200, the resistance con-
stants of pipes 1-2, 6-2, 3-7, and 3-4, respectively.
Even though the Hardy Cross method is simple, easy to apply, and
therefore, suitable for hand calculation, its greatest drawback is the pres-
ence of the convergence problems. It is observed to converge very slowly
and therefore requires large number of iterations. The number of iterations
increases as the size of the network increases. In some situations the
method starts diverging and does not converge at all.
In addition to the size of the network, the number of iterations required
for convergence depends upon (1) the convergence criterion. (2) the pipe
resistance constants, and (3) the pipe discharges. Lower values of conver-
gence criterion, larger values of some pipe resistances (small diameters or
large lengths of pipes or both) in comparison with the R values of other
pipes, and lower values of discharges in some pipes tend to increase the
number of iterations. The convergence criterion depends upon the accu-
racy one desires, consistent with the time, effort, and money one is pre-
pared to spend. Techniques to improve convergence when the delay is due
to unusually large pipe resistances or unusually low pipe discharges have
been developed by several investigators [110-113]. Some of them are briefly
described herein.
228 HARDr' CROSS METHOD
Voyles and Wilke [110] have shown that the convergence is improved in
the method of balancing heads when the loops are selected such that the
sum of the common flow resistance factors is minimum. For example, in
Figure 7.5 of Chapter 7 if ^, 4 is large (300 instead of 30), it is preferable
to select loop 2 comprising pipes 1, 2, 5, 7, 6, and 3 instead of pipes 4 . . .
7. In other words, pipes with large R values are externalized.
For the method of balancing flows, Dillingham [111] treated a pipe with
a large R value as an elongated joint and applied the computed head cor-
rection to both the end nodes of the pipe.
When the discharge in a pipe is extremely low, the head loss through it
is negligible. This results in a very large value of Q/h for this pipe in the
denominator of the AH correction equation [Equation (8.16)]. This gives
an extremely small value for the AH correction, thus delaying conver-
gence. Dillingham suggested the use of the linear head loss relationship
h = RQ for such a pipe so that the corresponding Q/h term can never ex-
ceed the MR value of this pipe and thus convergence is improved.
Williams [106] exclusively studied the problem of convergence in the
method of balancing heads and suggested the use of a convergence acceler-
ation factor A, given by
in which pipes represent the number of pipes in the network, i.e., X; and
total pipes represent the total number of pipes used in the iteration. Thus,
a pipe common to two loops is considered twice. For example, for the net-
work of Figure 7.5, pipes would equal 7, while total pipes would equal 10,
since pipes 1, 3, and 4 are common to two loops. The correction obtained
at the end of each iteration, ,AQC [Equation (8.5)] is multiplied by A and
the resulting value A • ,AQC is used to modify the pipe flows.
Instead of using a common fector, different factors for individual loops
or nodes can also be used. The modified Hardy Cross method described
earlier modifies individual corrections. McCormick and Bellamy [112] and
Van den Berg as reported by Epp and Fowler [107] suggested special mea-
sures to maximize individual corrections and thereby improve conver-
gence.
For the method of balancing flows, Kootattep and Aya [113] have sug-
gested the use of buffered successive over relaxation, i.e., the use of an
accelerator and decelerator.
If the present correction ,AHj is in the same direction as that of the pre-
vious correction t -iA//,, i.e., if both the corrections are of the same sign,
the present correction is accelerated to enhance convergence. On the other
hand, if the present correction is in the reverse direction and therefore of
Problems 229
in which
The appropriate value of accelerator A i is between 1.3 and 1.4 and that
of decelerator A2 between 0.5 and 1.0. They have recommended At as 1.3
and A2 as 0.5. They have named the method as buffered successive over re-
laxation because the work of the decelerator looks like that of a buffer
against excessive over correction.
This principle can also be used for the method of balancing heads.
8.7 PROBLEMS
9.1 INTRODUCTION
b
F(a + b) = F(a) + bF'(a) + ^rF" y
2! ' ' (n - 1)!
231
232 NEWTON-RAPHSON METHOD
in which F', F'",..., F" denote the first, second, . . . , « " " derivative of
F, respectively, and the last term denotes the remainder after n terms. Us-
ing the remainder after two terms, the finite Taylor's series becomes
Now if tAx is relatively small in comparison with ,x, (tAx)2 will still be
smaller, and therefore we can neglect the remainder after two terms. The
resulting expression is
,F'(x)-,Ax = -,F(x) (9.5)
from which the expression for the correction is
[ F(x)
Basic Concepts 233
F1(xltx1) = 0 (9.9a)
and
F2(x,,x2) = 0 (9.9b)
Let the t'h trial values of*! and x2, be ,x, and ,x2 and the corrections be
,Ax, and ,Ax 2 , respectively. Following the procedure outlined earlier, and
assuming ,Ax, and ,Ax2 to be small, and therefore neglecting their
234 NEWTON-RAPHSON METHOD
(9.10a)
and
i i i i
3F,
A.x, F,
dXi dx2 ' ' • dxn
dx2 dxn
The first matrix is the coefficient matrix, also called the Jacobian of the
« functions F,, F2, . . ., Fn with respect to the n variables xt, x2 xn
for the t'h iteration. The second matrix is a column matrix representing the
corrections Ax,, Ax2, . . ., Axn for the /'* iteration, and the matrix on the
right-hand side is a column matrix giving the values or the residues of the
functions Fu F2, . . ., Fn for the t'h iteration.
The application of the Newton-Raphson method to the analysis of water
distribution networks was first proposed by Martin and Peters [116] and
since then it has been extensively used in research and in the field f88,91,
107,112,117-119]. The Newton-Raphson method can be applied to any type
of equations described in Chapter 7. However, it is commonly applied in
practice to the nodal-head equations, i.e., H equations, and to the loop-
flow correction equations, i.e., AQ equations.
Example 9.1
Solution
The network reproduced herein as Figure: Example 9.1, has two source
0-1
nodes 1 and 2, with known heads, Hol = 100.00 m and Ho2 = 95.00 m.
The heads H3 . . . H6 are unknown.
Because n = 1.85, (l/n) = 0.54. Therefore, the //equations at nodes
3 . . . 6 are, respectively,
F3 =
0.2 = 0 (2)
(3)
95 - H6 - // 5
F6 = -0.1 = 0 (4)
20 50
Head Equations 237
dH3 ~ (40)
°- 5 ^(ioo-ftr-«(-i)- (10) 05 - 95)-
dF4 0.54
^^(100-//4)-046(-l)
8H4 (20)
0.54 ( f t - Hs)-°46
(40)0.54
dF4
dH6
dH4 ~ (40)°54 (
~ (40)
4 5)
0.54
(H6 -
dH6
3F
dH3
-=O;4£ =O
m.
dF6 0.54
(H6 -
dH5 (50)°-s'
0.54 0.54
(95 - - Hsy
(5O)05'
Let the first iteration heads for the nodes be YH3 = 96.00 m, tH4 =
94.00 m, ,/Zs = 92.00 m, and ,/7 6 = 93.00 m as assumed in Example 8.9.
Substituting these first-iteration H values and evaluating the values of the
elements of the Jacobian and also the F3 . . . F6 values, Equation (5) be-
comes, for the first iteration
-0.1946 0 0 0
0 -0.1866 0.0536 0
0 0.0536 -0.1189 0.0653
0 0 0.0653 -0.1432
-0.3000
iA// 4 0.2827
X (6)
.A// 5 -0.0813
0.0673
(2) The elements of the coefficient matrix are symmetrical about the
main diagonal, i.e., (d/\/d//,) = (dF,/3//,), i =£ j . If node / is con-
nected to node y through pipe ij, the value of (3F,/d//,) is obtained
through the term denoting discharge in pipe ij. Because the same
pipe is considered for evaluating {dFjIdF,), these two elements are
equal. Thus, only the elements on the main diagonal and in the trian-
gle above it, i.e., upper triangle, need be evaluated. The elements in
the lower triangle can be written because of symmetry.
(3) The elements on the main diagonal, dFJdHj, i = j are numerically
the same as the corresponding Y.(Quln • hu) values in Table: Exam-
ple 8.9. For example, 3FJdHA and T,{QJn • hiA) are numerically
equal (0.1866). Similarly, the nondiagonal elements, i.e., dFJdHj,
i ± j , are the same as the corresponding Qijln • hu values in Table:
Example 8.9. For example, 8FJdH5 = Q4-s/n • h4-s = 0.0536.
Further, the values in Example 8.15 (modified Hardy Cross method)
are identical with the corresponding values in this example. How-
ever, the Newton-Raphson method solves the equations exactly,
while the modified Hardy Cross method solves them approximately.
Example 9.2
Solution
1. H-q-R Equations
(2)
100 -
10
(3)
inn
IUU H \o.5 Iu u \o.5
~ "4 1 , {^3 - HA - 0.5 = 0
I + I
27.5 I
(4)
98 - // 5 \ 0 5
- 0.7 = 0 (5)
40
For the five corrections, AH3, AH4, A// 5 , Aqu and A/?4, the five equa-
tions in matrix form are
dH3
0 0 A// 4
dF3
'\dH3 1 1
Head Equations 241
11
I dF4 dF4 6t4
dH3 dH4 3R4
dFs
AR4
m4 dHs dR4
Let the first-iteration values be the same as those in Example 8.10 so that
tH3 = 97.00 m, ,// 4 = 96.00 m, ,H5 = 95.00 m, and ,/? 4 = 20. Let
,<7i = 1.00 mVs. [Actually, applying the flow continuity equation to the en-
tire network, we get q, — 1.2 rrrVs. Therefore, qx could have been treated
as a known and the unknowns could be the remaining four parameters for
which Equations (2) . . . (5) are the required equations. However, tqt is
deliberately taken different from 1.2 m3/s to illustrate the method.] There-
fore, the five equations for the first iteration are
0.0913 0.0884 () 1 0
0.1768 0 0.()456 0 0
0.3634 0.0953 () 0 0
0.0953 -0.2955 0. 118 0 0.00559
0 0.1118 -0. 1574 0 -0.00559
1
AH3 ' -o.
2548 '
AH4 -0.0274
X AHS 0.1106 (7)
Aqi 0.1742
AR* . L -0.2025 .
The H-R equations are Equations (2) . (5) for which the correction
equations in matrix form are
dF2 BF2
AH3 F2
dH3
\dF3 dF^
! dH3 dH4 \\*H4\
242 NEWTON-RAPHSON METHOD
i i (8)
dF4 dF4 dF4 3F4
A// 5 F4
dH3 dH4 dHs dR4
27.5 40
(9)
Thus, Equations (2), (3), and (9) provide three equations for H3, H4,
and Hs. The corresponding correction equations are
0 AH3
dH3
dF3 dF3
dH4
AH4 F3 (10)
dH3
It has been shown in Chapter 7 that an additional equation for the head
supplied by a pump can be formulated [for example, Equation (7.31)].
However, as explained in Example 8.11, the head delivered by a pump, hp,
Head Equations 243
is initially assumed, and the head at the upstream node of the pipe in
which the pump is located is correspondingly increased. If the head-
discharge relationship is assumed for the pump such that it is consistent
with the head loss relationship for pipes, the pump can be made an integral
part of the pipe in which the pump is located.
Example 9.3
Solution
Sump,
W.L.=80-00 m : ESR
W.L=100-00m
Figure: Example 8.6 reproduced herein as Figure: Example 9.3 are as fol-
lows:
80.00 + hB- H3\05 (H3 - HA05 - H5\05
R3
- 0.08 = 0 (1)
<H3 - HMS
- 0.12 = 0 (2)
F t = !±=JlL
- 0.12 = 0 (3)
- 0.12 = 0 (4)
244 NEWTON-RAPHSON METHOD
Because the pump is situated in pipe 1, the head supplied by it, hp, is in-
cluded in pipe 1 as shown in Equation (1). As explained in Example 8.11,
the R values of the pipes change from iteration to iteration. They are ini-
tially obtained by taking / = 0.02 and are updated for each subsequent
iteration. Further, as hp is also assumed and updated, the nodal heads
H3 . . . H6 are the variables for which Equations (1) . . . (4) are available.
The equations for the AH corrections in matrix form are
For the first iteration, XH3 = 85.00 m, ,// 4 = 83.00 m, ,// 5 = 80.00 m,
and ,// 6 = 90.00 m. Further, ,hp = 10.75 m, ,R, = 204, ,R2 = 2,066,
>R3 = 8,707, ,/? 4 = 3,265, ,/? 5 = 755, ,/? 6 = 3,265, and ,/? 7 = 157, as
in Example 8.11. Using these first-iteration values, the elements of the co-
efficient matrix in Equation (1) are calculated and the equations are solved
to get the head corrections and subsequently the nodal heads for the sec-
ond iteration. The pipe discharges obtained in the first iteration are used
in the second iteration to determine the friction coefficients, the pipe resis-
tances, and the power supplied by the pump 2hp . The iterative procedure
is continued and the solution given earlier is finally obtained.
PRV and the assumed head of the downstream node of the pipe. This dis-
charge is then taken as the additional outflow at the upstream node of the
pipe. If the PRV becomes inoperative or prevents flow in the reverse direc-
tion, the procedure described earlier for the check valve is followed. A
check is made at the end of each iteration to verify that the PRV has re-
mained in the assumed mode.
Example 9.4
Solution
Using the known head for node 1 and the set-pressure head of the PRV,
0-3 rcfe
100-OOnn t
1
° R=2 2
0-6
o-2 y\ 5Q 6
100.00 - // 2 \ 0 5 H2 -
F2 =
- 0.3 = 0 (1)
H2 - H3\°-s (H3 - K
F3 = - (0.2 + ft.,,) = 0 (2)
40
- 0.5 = 0
40
(3)
#4 - / / 5 \ ° 5 / t t , - // 6 \ 0 5
- 0.4 = 0 (4)
9 / \ 12
in which Q3-u = discharge in the pipe on the upstream of the PRV and
Rd-6 = resistance of the pipe on the downstream of the PRV.
As the PRV is situated just downstream of node 3, Rd-t = R3~6 = 50.
Assuming the first-iteration H values as in Example 8.12, Q3-u = Qd-6 =
0.3162 nWs. Equations (1) . . . (5) are solved by the Newton-Raphson
method, the corrections are evaluated, and the nodal heads H2 . . . H6 are
corrected. Discharge in pipe 3-6 is determined, the pressure head at the
downstream end of the PRV is calculated, and the mode of the PRV is de-
termined. If the PRV is in the operative mode, Q3-u is determined and the
next iteration is carried out. If the PRV becomes inoperative, Q3-u in
Equation (2) and the term 1(80.00 - H6)/Rd-6]°* in Equation (5) are re-
placed by [(//, - //6)/(/?3-6 + /?PRV)J°-S. If the PRV behaves as a check
valve, these terms in Equations (2) and (5) are dropped. Thus, at the end
of each iteration, the operational mode of the PRV is determined and the
appropriate forms of Equations (2) and (5) are selected.
Example 9.5
Solution
3, Q 3
, i
.2 Q* R
3 J0_ /T95-00 6 22
q 0-1
III z i
. i
as the loop flow corrections in the clockwise direction for the three loops
I . . . Ill, respectively, the AQ equations are
F, = 20(0.45 + AQ,) 185 - 30(0.05 - AQ, + AQ,,)' 85
The basic unknowns are AQ, . . . AQ,,, which will be assumed and
successively corrected by applying corrections. Let these additive correc-
tions be z, . . . 2ni to the assumed AQ, . . . AQ,,, values, respectively.
The corrections z, . . . Zm will be successively evaluated until they
become less than some prespecified value, and the AQ, . . . AQ,,, values
stabilize. The correction equations are
dF,,, dF,,,
3AG, 3AGn
248 NEWTON-RAPHSON METHOD
ir = 20 X 1.85(0.45 + AC?,)085
dAQ,
dFm
= 40 X 1.85(0.35 - AC?, + AC?,n)08
SAG,,
For the first iteration taking ,AC?i = iAC?>> = iAC?m = 0, Equation (4)
becomes
iAC?,,, and iAQm, respectively. Because the XAQ values are zero, the 2AC?
values are 2A£?i = 0.0146 m3/s, 2AC?n = -0.0495 m 3 /s, and 2AC?,n =
— 0.0098 m 3 /s. Taking these AC? values, the elements of the coefficient
matrix and the F values are reevaluated to get correction equations, similar
to Equation (14), for the next iteration. The iterative procedure is con-
tinued and the results are shown in Table: Example 9.5. The answers given
earlier in Example 8.3 are finally obtained.
9.4.2 NETWORKS WITH UNKNOWN PIPE RESISTANCES
The application of the loop equations to networks with unknown pipe re-
sistances was first suggested by Gofman and Rodeh [91]. They used graph
theory in formulating the appropriate equations.
As seen in Chapter 7, we can formulate either the AQ-AR equations or
the AC? equations for networks with unknown pipe resistances. When AC?-
TABLE: EXAMPLE 9.5.
AR equations are formulated and solved, the unknown pipe resistances are
also assumed and corrected along with the pipe discharges and the method
is implicit. When the AQ equations are formulated and solved, the method
is explicit. The unknown pipe resistances are evaluated in the end after sat-
isfactory convergence is reached.
Example 9.6
Solution
Let the corrections to AQ, and AQ,, be z, and z,,, respectively, and the
correction to AR4 be r4. The correction equations in matrix form are
dF, dF,
u F,
dAQ, dAQ,,
(8)
For the first iteration, taking ,Ag, = AQu = 0, Equation (8) becomes
which yields ,2, = -0.00027 m3/s and ,z,, = -0.01468 nWs. This gives
2AQ, = -0.00027 m3/s and 2AQU = -0.01468 mVs.For these AQ
TABLE: EXAMPLE 9.6.
values, as the 2F, and 2F,, values are negligible, the final AQ values are
AQ, = -0.00027 mVs and Ag,, = -0.01468 nrVs. These corrections
are applied to loop I (pipes 1, 2, and 3) and loop II (pipes 6, 2, 1, 4, and
5) to obtain the final Ql . . . Q6 values that are as obtained earlier in Ex-
ample 8.4. Loop III comprising pipes 3 . . . 6, the overlapping loop com-
prising pipes 6, 2, 1,4, and 5 or the path comprising pipes 1, 4, and 5 can
then be used to evaluate the value of R4 that is as obtained earlier in Exam-
ple 8.4.
The three loops shown in Figure: Example 8.6 are considered and the
first-iteration R values are taken as shown in Table 2: Example 8.6. The
pipe discharges are also taken as shown therein. From Equations (2) and
(3) of the solution of Example 8.6, the head delivered by the pump is
given by
Using Equation (2), the AQ equations for the first iteration are
from which ,z, = 0.02397 m3/s, ,z,, = 0.01276 m3/s, and ,?,„ = 0.05289
m3/s. Therefore, 2A£>, = 0.02397 m3/s, 2A<2n = 0.01276 rrrVs,
and 2A<2,,, = 0.05289 m3/s. This gives 2Ql = 0.1671 nWs, 2Q2 = 0.05108
nrVs, 2 03 = 0.03603 m J /s, 2Q4 = 0.03121 nWs, 2QS = 0.1001 m3/s, 2Qb
= 0.05276 m3/s, and 2Q7 = 0.2729 m3/s. Using these discharges, the/and
R values of pipes for iteration 2 are obtained and Equations (3) . . . (5) are
rewritten using the revised R values. Equation (6) is then reformulated,
solved, and the z values are obtained. The iterative procedure is continued
and the answers given in Example 8.6 are finally obtained.
reverse direction, the procedure described earlier for check valves is fol-
lowed. A check must be made to ascertain the mode of the PRV.
Example 9.8
Solution
Note that, for maintaining continuity, the correction AQn is applied to the
pipes of loop III, i.e., pipes 4 . . . 7. Therefore, the correction AQn is ap-
plied to pipe 4 in Equation (1), but is not applied to pipes 1 and 2 in Equa-
tion (2).
The correction equations are
(3)
256 NEWTON-RAPHSON METHOD
For the first iteration taking ,A£?, = iAQu = 0, Equation (3) becomes
and
In the Hardy Cross method described in the last chapter, each loop or
node was considered separately, and the effect of adjacent loops or con-
necting nodes was neglected. Thus, each equation was an independent
equation containing only one unknown and therefore could be solved
separately. On the other hand, in the Newton-Raphson method described
in this chapter, all nodes or loops are considered simultaneously. Simi-
larly, in the linear theory method described in the next chapter, all pipes
Labelling Network Elements 257
Epp and Fowler applied the Newton-Raphson method for the solution of
AQ equations, and therefore proposed a labelling algorithm for the loops
of a network. Their algorithm is explained with the help of Figure 9.2 in
which the loops are initially labelled with letters for reference.
Their algorithm is as follows:
A list is made of all loops that contain a node of degree 2. (The degree
of a node is n if n is the number of pipes connected to it.) Usually, these
258 NEWTON-RAPHSON METHOD
L M N 0 P
IY VII X XUI XV
F G H J K
II V VIII XI XIV
A B
c D E
I III
YI IX XII
Figure 9.2 Labelling loops of a network using an approach suggested by Epp and Fowler.
loops occur at the corners of the network (Loops A, E, L, and P). All
pseudo loops are then added to this list. The first loop of this list is labelled
I (one). Its neighbors are labelled II, III, IV, . . . etc. (Loop A is labelled
I and adjacent loops F and B are labelled II and III, respectively.) Now
consider loop II. Its unlabelled neighbors are L and G. The labelled neigh-
bors of L consist of single loop II, whereas loop G has two labelled neigh-
bors, namely II and III. Thus, loop L which has the fewest labelled neigh-
boring loops is labelled IV and loop G is labelled V. All the neighbors of
loop II have now been labelled and therefore loop HI is considered next.
As loop III has only one unlabelled neighbor C, loop C is automatically
labelled VI. The process is continued until all loops are labelled as shown
in Figure 9.2.
In this simple illustrative network, the half bandwidth is three and is op-
timal. However, the above algorithm may not always produce a minimum
bandwidth. Therefore, the above procedure is repeated starting with each
loop having a node of degree 2. The numbering that produces the smallest
bandwidth is then selected as the best loop labelling.
Jeppson and Davis [83] have suggested an algorithm that can be used to
band the coefficient matrix of Q, AQ, or H equations. Initially, the pipes,
loops, or nodes are arbitrarily numbered and the algorithm suggests how
they should be renumbered for minimum bandwidth. The algorithm is ex-
plained for H equations for a one-source, nine-demand node network
shown in Figure 9.3. The source node is labelled 0 (zero) and the demand
nodes are labelled 1 . . . 9 as shown so that the H unknowns are //, . . .
H9. The nine //equations obtained for the node-flow continuity relation-
ship at nodes 1 . . . 9 are labelled (1) . . . (9), respectively, so that the
equation numbers are identical with the node numbers for which they are
written. (This is the reason why the source node is labelled 0 herein). The
unknowns in Equations (1) . . . (9) are shown by • in Figure 9.4(a). This
i
(2)
U) (1)
Figure 9.3 Network for illustrating the algorithm of Jeppson and Davis for banding a coefficient
matrix.
Unknown Unknown
1 2 3 4 5 6 7 8 9 A; AI 1 2 3 4 5 6 7 8 9
1 • 9 3 1 3 * •
it 2 • • 82 2-2 • • •
S.3 • • • 6 1 3-5.1 • • •
• • • • 10 4 i'5 4 • • • •
• • 10 5 sl5 • • •
• • • • 14 8 6 iG 7 • • • • •
§7 • • • • • 106 7i 8 • • •- •
,3"8 • • • 10 7 8^6 • • • •
9 • • • 15 9 9 9 • • •
9 7 5 7 12 14 10 12 16
(0) 0
1
B! 4 3 1 2 6 8 5
(bl
7 9
Position Position
On ^
1 2 3 4 5 6 7 1 2 3 4 5 6 7 8 9
Unknown 8 9 £
Unknown
3 4 2 1 7 5 8 3 4 2 17 5 8 6 9 A-
3 6 9 Ai
4 Ai * ~ 3 * • •
1 1 4
62 2 ^ 2 • • • • 6
&
•S.1
• 10 4 3 -%_ 4 0 • • • 7
• 7 3 4 b 1 • • • 10
-05 • 10 5 5 -§ 5 • • • 10
£7 10 6 6 2 7 • • • • • 10
£8 • • 13 7 7 ;§ 8 • • • • 13
iS"6 14 8 8
^ 6 • • • • 14
9 • • 16 9 9 9 • 16
Bi 4 7 7 9 10 11 13 13 16
(c) d)
Figure 9A Procedure for obtaining a banded coefiicient matrix: (a) original coefficient matrix;
(b) matrix after changing positions of rows (equations); (c) matrix after changing positions of col-
umns (unknowns); (d) finally obtained banded coefficient matrix.
260 NEWTON-RAPHSON METHOD
figure, therefore, also depicts the nonzero elements in the coefficient ma-
trix. The numbers across the top denote the positions of the nine un-
knowns, Hi . . . H9, while the numbers at the left denote the positions of
the nine equations, Equations (1) . . . (9).
The iterative procedure for reducing the bandwidth consists of two
steps: (1) changing the positions of the equations, and (2) changing the
positions of the unknowns.
1. Changing Positions of Equations. Each row (equation) is given a
weight equal to the sum of the position numbers of the first and last non-
zero elements in that equation, as shown in column At [Figure 9.4(a)]. The
weights in column At are next arranged in ascending order of magnitude
(column A /) and the equations are changed to positions given by this order
[Figure 9.4(b)].
When weights of two or more arrays are equal, they are assigned posi-
tions arbitrarily. It is presumed that the ascending order is maintained as
long as the weight of an array is not less than the weight in the earlier
array.
2. Changing Positions of Unknowns. Each column (unknown) is given a
weight equal to the sum of the position numbers of the top and bottom non-
zero elements in the column as shown by row fi, [Figure 9.4(b)]. The
weights in row 5, are next arranged in ascending order of magnitude (row
B-) and the positions of the unknowns are reassigned [Figure 9.4(c)].
These two steps of giving weights and reassigning the positions of the
equations or unknowns in the ascending order of the weights is repeated
iteratively until all the equations or unknowns are arranged in ascending
order and no change in their positions occurs, or until some maximum al-
lowed number of such iterations have been completed to avoid alternation
of the positions of two or more equations or unknowns between consecu-
tive iterations.
One more iteration of changing positions of the rows of the matrix in
Figure 9.4(c) results in the matrix shown in Figure 9.4(d). For this matrix
it is observed that all #, values are in the ascending order. Since all At val-
ues are already in the ascending order, as shown again by column At in
Figure 9.4(d), the iterative procedure is complete.
The original coefficient matrix of Figure 9.4(a) is symmetrical while that
in Figure 9.4(d) is not. For a symmetrical matrix, the relative positions of
each equation and corresponding unknown are identical. Therefore, to
maintain the identical positions, the ordering of the unknowns is the same
as that of the equations. Therefore, only At (or Bt) are obtained and the
positions of the equations and unknowns are identically changed. For the
network of Figure 9.3, the different symmetrical coefficient matrices at dif-
ferent ordering stages are shown in Figure 9.5. The node labels as obtained
from the final coefficient matrix for the symmetrical matrix with minimum
bandwidth [Figure 9.5(b)] are shown in parentheses in Figure 9.3.
Problems 261
Position Position
4 5 6 89 I
Unknown Unknown
8 6 9 AiA'i S. 3 2 4 1 7 5 8 6 9 A ;
3 • • • 5 1 1_ 3 • • • 4
£2 • • • • 7 2 2^2 • • • • 6
• • • 9 4 3-S.4 • • • 7
• • • • 7 3 4 b 1 • • • 9
• • • 12 6 5 "§7 • • • • • 10
5 • • • • • 105 6*2 5 • • • 11
28 • • • • 12 7 7z 8 • • • 13
• • • • 148 8uf 6 • • • • 14
• • • 169 9 9 • • • 16
Bi 5 7 9 7 12 10 12 14 16 Bi 4 6 7 9 10 11 13 14 16
(a) (b)
Figure 9.5 Procedure for obtaining symmetrical banded coefficient matrix: (a) after first order-
ing; (b) finally obtained symmetrical banded coefficient matrix.
9.6 PROBLEMS
10.1 INTRODUCTION
E have seen in Chapter 7 that at least some (if not all) equations in-
W volved in pipe network analysis are nonlinear. Because no direct
method is available for their solution, these equations are linearized and
then solved. Naturally, the solution is approximate, and therefore it is cor-
rected, and the iterative procedure is continued until satisfactory accuracy
is reached.
The Newton-Raphson method described in the previous chapter expands
the nonlinear terms in Taylor's series, neglects the residues after two
terms, and thereby considers only the linear terms. As such, the Newton-
Raphson method linearizes the nonlinear equations through partial differ-
entiation. This method is general and works even when the nonlinear
equations are transcendental containing exponential, trigonometric, hyper-
bolic, or logarithmic terms. The nonlinearity in the equations for pipe net-
work analysis, however, is algebraic, uniform and simple: the variables are
raised to the same, nonunity exponent. For example, the nonlinear Q equa-
tions obtained from the loop-head loss relationship contain terms, RXQ",
and the H equations obtained from the node-flow continuity relationship
contain terms, [(//, - //,• )//?<,•]"". In both these types of nonlinearity, the
value of n is the same for all terms—1.852 in the Hazen-Williams pipe-
head loss formula and 2.0 in the Darcy-Weisbach and Manning formulas.
This is a helpful feature and the nonlinear term can be easily linearized by
merging a part of the nonlinear term into the pipe resistance constant as in
Equation (6.5). This principle was first suggested and used by Mcllroy
[121], Marlow et al. [122], and Muir [123]. Later, Wood and Charles [124]
developed it, and it is now widely used in practice [124-129].
Even though the principle can be used to linearize all types of equations
including A£> equations [121,130], the method is applied in practice to
pipe-discharge equations, i.e., Q equations, and the nodal-head equations,
i.e., H equations.
263
264 LINEAR THEORY METHOD
The pipe discharges are taken as the basic unknowns in formulating the
Q equations. The application of the linear theory method for the solution
of these equations is described herein for different situations.
")G. = 0, c = 1 . . . C (10.1)
2 MO, = 0, c = 1 . . . C (10.2)
in which ,/?,' = modified resistance constant of pipe x for the t'h iteration
and is given by
2 R^Q: = 0, c = 1 . . . C (10.5)
Pipe Discharge Equations 265
,Gx = I , M , a X a , - . ) i 0 5 , t = 2, 3 . . . (10.6)
Q ( 2) gx
,6x = "" 2 """ , x = 1 . . . X (10.8)
Instead of taking the average of the obtained values in the previous two
iterations, taking the average of the assumed and obtained values in the
previous iteration gives rapid convergence as shown later and also by Ex-
ample 10.1. Therefore,
Example 10.1
Solution
Head, demand,O-2rry's
Q, - Q3 - 0.3 = 0 (1)
0 5 + 0 7 - 0.4 = 0 (3)
0 6 - 0 7 - 0.1 = 0 (4)
From the two basic loops I and II and the pseudo loop HI, we have, re-
spectively,
[Equations (1) . . . (7) are the same as Equations (7.12a). . . (7.12g) ex-
cept that n is replaced by 1.85.]
For the first iteration, XQX = . . . = xQi = 1 and therefore, Equations
(5) . . . (7) become, respectively,
4 0 0 , + 1003 - 5 = 0 (10)
Pipe Discharge Equations 267
Pipe Assumed Obtained Assumed Obtained Assumed Obtained Assumed Obtained Assumed Obtained
Pipe Assumed Obtained Assumed Obtained Assumed Obtained Assumed Obtained Assumed Obtained
XUR and the approach is implicit. When the Q equations are formulated
and solved, the number of the unknowns is X and the approach is explicit.
In the loop-head loss equations for the implicit approach, the terms indi-
cating the head loss through pipes of known resistance are linearized as
usual for pipe discharges and those indicating head loss through pipes of
unknown resistance are linearized for pipe resistances. Thus,
o-smys c O-7
t 4
R4
node.1^
a
evil
to m
WO-OOmpN. a'
to m
0-6
^ ^
Qi ~ Q3 + Qe - 0.6 = 0 (1)
Qi + Q3 - QA - 0.5 = 0 (2)
Q4 + Q, - 0.7 = 0 (31
or by
for the pseudo loop comprising pipes 1, 4, and 5. Thus, Equations (1) . . .
(4), (7), and (12) provide six 0 equations for the six unknowns 0 , . . . 0 6 .
These equations are then solved by the linear theory method to obtain the
final 0 values and subsequently the value of R3.
Example 10.3
Sump
W.L.=80-00
274 LINEAR THEORY METHOD
Solution
The network has seven pipes and therefore the basic 0 unknowns are
Qi . . . Q7. From the node-flow continuity relationship at nodes 3, 4, 5,
and 6, of Figure: Example 8.6, reproduced herein as Figure: Example 10.3,
we get, respectively,
G. - Gi - a - 0.08 = o (i)
03 + Q* + ft - 0.12 = 0 (3)
Taking/ = 0.02 for all pipes, thefirst-iterationR values of the pipes are
obtained as shown in Table: Example 10.3. Using these R values, the head
loss equations for loops I . . . Ill are, respectively,
Pipe iOx i ' . iRx ifli Q«,, if, 2RX 2R'X Ox(2) 3Ox »f. 3RX 3R'X 0x,3,
1 1 0.02 204 204 0.1467 0.1467 0.01859 189.7 27.82 0.1779 0.1623 0.01853 189.0 30.68 0.1617
2 1 0.02 2066 2066 0.0433 0.0433 0.02180 2252 97.51 0.0581 0.0507 0.02167 2238 113.48 0.0511
3 1 0.02 8707 8707 0.0234 0.0234 0.02351 10235 239.5 0.0398 0.0316 0.02324 10118 319.72 0.0306
4 1 0.02 3265 3265 0.0350 0.0350 0.02316 3782 132.4 0.0292 0.0321 0.02323 3792 121.7 0.0327
5 1 0.02 775 775 0.1117 0.1117 0.02121 821.7 91.8 0.0911 0.1014 0.02125 823.3 83.48 0.1016
6 1 0.02 3265 3265 0.0615 0.0615 0.02388 3899 239.8 0.0510 0.0562 0.02392 3905 219.46 0.0567
7 1 0.02 157 157 0.2933 0.2933 0.01774 139.5 40.9 0.2621 0.2777 0.01776 139.7 38.80 0.2783
276 LINEAR THEORY METHOD
Example 10.4
Solution
5Q 6 o-6
replaced by source and sink nodes and a pseudo loop (loop II in Figure:
Example 8.7) is formed. Thus 0 equations for nodes 2 . . . 6 and loops I
and II are, respectively,
G. - Qi - & - 0.3 = o (1)
Q6 + Q7 - 0.6 = 0 (5)
Equations (6) and (7) are linearized and solved simultaneously with
Equations (1) . . . (5). Equations (6) and (7) are modified and the iterative
procedure is continued. Finally, the answers given in Example 8.7 are ob-
tained. As initially assumed, the PRV is in the operative mode.
When the PRV becomes inoperative for another pressure setting given
in Example 8.7, loop II is replaced by loop III shown in Figure: Example
8.7. Therefore, Equation (7) is replaced by
Note that the resistance constant of pipe 7 now includes the resistance
constant of the PRV also.
278 LINEAR THEORY METHOD
Equations (1) . . . (6) and (8) are solved by linear theory and the an-
swers given in Example 8.7 are obtained.
For the changed nodal flows in Example 8.7, the PRV behaves as a check
valve. Here, pipe 7 containing the PRV is dropped. Therefore, the term Q7
is dropped in Equations (1) . . . (6). Further, as loop II containing pipe 7
is now absent, Equation (7) or (8) for loop II is also absent. Therefore, the
required Q equations for basic unknowns Qi . . . Q6 are provided by
Equations (1) . . . (6) in which Q7 is absent. These six equations are
solved by linear theory and the answers given in Example 8.7 are obtained.
The nodal heads are taken as the basic unknowns in the formulation of
//equations as explained in Chapter 7. The application of the linear theory
method for the solution of these equations for different situations is ex-
plained herein.
10.3.1 NETWORKS WITH KNOWN PIPE RESISTANCES
H — H\l/"
J
—^ -\ + qoj = 0, for a l l ; of known q3 (10.12)
/
connected to
j through x
in which //, = Hoi and //, = Hoj for source nodes. Therefore, for the tth
iteration, Equation (10.12) can be written as
C" — ,x = 1 . . .X (10.15)
1^ x —
or
,0, x = l (10.16)
in which ,QX, thx = discharge and head loss in pipe x for the fh iteration,
respectively.
To begin the iterative procedure, it is necessary to initialize the dis-
charge and head loss in a pipe, i.e., to select the values of XQX and thx,
x = 1 . . . X. Collins and Johnson [125] suggest that tQx should be
chosen to correspond to the Reynolds number of 200,000. As Re =
AQI-KDV, we get
= 5O,0OOTDXI»,JC = 1 (10.17)
Taking r = l x 10" 6 m 2 / s , w e g e t
lrD
r> ' 1 v (10.18)
i Q* = - j o " , x = l . . . X
in which the discharge and diameter are in cubic meters per second and
meters, respectively.
Issacs and Mills [126] suggest that the initial discharge may be set to the
same value for all pipes. They observed that an initial pipe discharge rang-
ing from 0.001 m3/s to 1 m3/s, when used in a test network of sixty pipes,
had no significant effect on the solution [126]. Alternatively, they also sug-
gest using an initial velocity of 1 m/s in all pipes to initialize the pipe flows
[126-127]. (The initial velocity in the approach of Collins and Johnson is
0.2/D,.) The value of the initial head loss xhx corresponding to the
assumed discharge i Qx is then evaluated from
I A, = Ra = 1 . • •X (10.19)
tain Hj values at the end of iteration 1, i.e., Hjm values. These Hj{X)
values are then used to determine C^(i) values from Equation (10.15) or
Equation (10.16). In the iterative procedure, if the values of ,C'X are taken
equal to ,C,' (t - 1) , overcorrection occurs as observed by Collins and John-
son. Therefore, they suggest the use of the average of the assumed and ob-
tained values. Thus, for the t'h iteration,
The iterative procedure is continued until the assumed and obtained C'x
values and therefore the //, values are sufficiently close. The unknown
nodal flows are obtained in the end.
If the unknown nodal flows are also considered as basic unknowns in ad-
dition to the unknown nodal heads as done by Collins and Johnson, we get
H-q equations, one for each node. These equations are linearized and
solved iteratively so that all the unknown nodal heads and unknown nodal
flows are obtained simultaneously at the end of each iteration.
As shown by Collins and Johnson [125], a similarity exists between a
structural network and a hydraulic network such as a water distribution
network. These similarities are as follows:
(1) The algebraic sum of the forces at a joint of a structural network
must be zero to maintain equilibrium. Similarly, the algebraic sum of
the flows at a node of a hydraulic network must be zero to maintain
flow continuity.
(2) The displacement of the ends of the members meeting at a joint of a
structural network must be the same. Similarly, the heads at the ends
of the pipes meeting at a node of a hydraulic network must be the
same.
(3) The force-displacement relationship specified by the geometric and
elastic properties of the members must be satisfied in a structural
network. Thus,
F = Ku (10.21)
Q = C •h (10.22)
Example 10.5
Solve Example 8.3 (Example 8.9) by the linear theory method based on
(1) H equations and (2) H-q equations.
Solution
1. H Equations
Head
Example 10.5, the linearized H equations written for the node-flow conti-
nuity relationship at nodes 3 . . . 6 are, respectively,
a +a o o o H3
o a +a +a -a o H4
o -CL c; + a -a Hs
o o -a a + c; H6
As the pipe diameters are not given, we assume the discharges in pipes
to be 0.1 m 3 /s. The corresponding hx and iC« values are obtained from
Equations (10.19) and (10.16), respectively, and are shown in Table: Exam-
ple 10.5. Therefore, for the first iteration Equation (9) becomes
Iteration 1 Iteration 2
2 c;
R<>
Pipe A ic; *«2, Q*,2,
C
x(2,
1 40 0.1 0.5650 0.177 4.34 0.301 0.0694 0.1232 4.52 0.308 0.0681
2 20 0.1 0.2825 0.354 3.02 0.360 0.1192 0.2366 3.37 0.382 0.1133
3 10 0.1 0.1413 0.708 0.66 0.230 0.3486 0.5283 0.48 0.194 0.4036
4 30 0.1 0.4238 0.236 -1.98 -0.230 0.1162 0.1761 -1.63 -0.207 0.1271
5 40 0.1 0.5650 0.177 2.27 0.212 0.0934 0.1352 2.31 0.214 0.0927
6 20 0.1 0.2825 0.354 0.28 0.100 0.3554 0.3547 0.53 0.141 0.2651
7 50 0.1 0.7063 0.142 0.01 0.010 1.0013 0.5717 0.15 0.043 0.2885
284 LINEAR THEORY METHOD
tained. The average of , C'x and C*'m is taken as 2CX' for the second itera-
tion. Therefore, Equation (9) now becomes
2. H-q Equations
Equations (12), (13), and (1) . . . (4) give the required linearized H-q
equations. These equations are expressed in the matrix form as
Ci + d 0 -Ci -Ci 0 0
0 a + c; + c6' -Ci -c; 0 -Ci
— Ci -Ci ci + a 0 0 0
-Ci C '
L-4
0 a +-Cia + a -Ci 0
0 0 0 C :.' + c; -Ci
0 — Ci 0 0
-a Ci + Ci
100"
95
H3 -0.3
(14)
H4 -0.2
H, -0.4
H6 -0.1
Nodal Head Equations 285
Using the XC'X values from Table: Example 10.5, Equation (14) becomes
100 ' •
95 q2
H3 = -0.3
X (15)
H4 -0.2
H5 -0.4
H6 -0.1
known at source nodes and unknown at demand nodes. The column matrix
on the right-hand side represents the known flows at demand nodes and
unknown flows at source nodes, the preceding sign for the nodal flows be-
ing negative for outflows and positive for inflows.
The procedure for writing the coefficient matrix as described above for
H-q equations [Equation (14)] can also be followed for writing the coeffi-
cient matrix for H equations [Equation (9)]. However, when the node
under reference, node^, is connected to another node of known head, node
i through pipe ij, the corresponding element of the column matrix on the
right-hand side is increased by the product of Hoi and C,y. For example, as
node 4 is connected to nodes 1 and 2 of known heads of 100 m and 95 m,
through pipes 1 and 4, respectively, the corresponding element in the
matrix on the right-hand side becomes lOOCV + 95C4' - 0.2 as shown in
Equation (9). The coefficient matrix for H equations is symmetrical but
the row sums and column sums are no longer zero.
Equation (6.2) giving the head loss in a pipe is linearized for Rx and
therefore can be expressed as
or
f R, (10.24)
or
Qx = | ,B, | • Rx (10.25)
Nodal Head Equations 287
in which
Example 10.6
Solve Example 8.4 (Example 8.10) by the linear theory method based on
(1) implicit and (2) explicit approaches.
Solution
1. Implicit Approach
0.6 - 98(C5' + i
98C6' - 0.6 (6)
-0.5
98C5' - 0.7
Initially, let the t Qx values be 0.1 m3/s. The corresponding xhx and i C'x
values are evaluated and shown in Table: Example 10.6. Therefore, for the
first iteration Equation (6) becomes
146.4
121.9 (7)
-0.5
23.8
2. Exploit Approach
Iteration 1 Iteration 2
1 8 0.1 0.08C 1.250 0.01 0.0354 3.536 2.3930 0.147 0.1356 0.9221
2 10 0.1 0.100 1.000 1.19 0.3450 0.290 0.6450 1.315 0.3626 0.2758
3 27.5 0.1 0.275 0.364 -1.18 -0.2071 0.1755 0.2697 -1.168 -0.2061 0.1764
5 40 0.1 0.40C 0.250 6.46 0.4019 0.0622 0.1561 7.449 0.4315 0.0579
6 8 0.1 0.08C 1.250 -0.81 -0.3182 0.3928 0.8214 -0.685 -0.2926 0.4272
Iteration 1 Iteration 2
-C«' 0 -C5' H3
c2' + a + c6' -a o H4
These equations are iteratively solved and the answers given in Example
8.6 are obtained.
The formulation of nodal head equations for networks with pumps has
been explained in Chapter 7. As stated therein, it is preferable to use the
head-discharge relationship that is consistent with the head loss relation-
ship for the pipes.
Example 10.7
Solve Example 8.6 (Example 8.11) by the linear theory method based on
head equations.
1
Sump ESR.
W.L=8O-0Om WL=i00-00m
0-12
Nodal Head Equations 291
Solution
c + a + c, -a -a o
-a a + a + CL -C; -a
-a -a a +a +a -a
o -c -a CL + a + a
The H equations in the matrix form as given by Equation (5) could have
been directly written as explained in Example 10.5. Herein, the sump head
of 80 m is increased by the pump head hp in the first term of the matrix on
the right-hand side.
Initially, the Reynolds number is taken as 200,000 in all pipes as pro-
posed by Collins and Johnson [125]. Therefore, discharges in pipes, i.e.,
the iQ, values, are obtained using Equation (10.18), and are shown in
Table: Example 10.7. Corresponding to these discharges, the lfx values are
obtained from the Swamee-Jain pipe-head loss explicit relationship. The
corresponding thx values are obtained from the Darcy-Weisbach head loss
relationship, and finally the XC'X values are obtained using Equation
(10.16).
The head developed by the pump is given by
Iteration 1 Iteration 2
Pipe A '*> C
x(1> 'x(2)
1 0.04712 0.01976 0.4476 0.10527 0.693 0.02009 0.05815 0.08390 0.09458 1.115 0.01902 0.07581 0.06799
2 0.03142 0.02213 2.2570 0.01392 -1.072 0.02245 -0.02150 0.02006 0.01699 0.353 0.02354 0.01205 0.03413
3 0.02356 0.02350 5.6785 0.00415 1.889 0.02400 0.01345 0.00712 0.00563 3.428 0.02359 0.01827 0.00533
4 0.02356 0.02350 2.1295 0.01106 2.961 0.02318 0.02798 0.00945 0.01025 3.075 0.02316 0.02852 0.00927
5 0.03142 0.02213 0.8463 0.03713 4.516 0.02129 0.07400 0.01639 0.02676 5.437 0.02124 0.08129 0.01495
6 0.02356 0.02449 2.2192 0.01062 7.477 0.02391 0.04377 0.00585 0.00823 8.512 0.02388 0.04673 0.00549
7 0.05498 0.01930 0.4590 0.11978 3.065 0.01804 0.14697 0.04795 0.08386 4.001 0.01794 0.16837 0.04208
*P
— — 12.04 — 12.03 — — — — 11.99 — — —
Nodal Head Equations 293
H3 9.609
H4 -0.120
(7)
H5 -0.120
H6 11.858
Solving these equations H3(1) = 91.347 m, # 4 ( , , = 92.419 m,
#5(i> = 89.458 m, and H6(l) = 96.935 m. These HjW values are used to
determine hx{1) values, fz(1) values with the help of Equation (2.40), and
then QxM values with the help of Equation (2.38). The Cx'(i> values are
obtained from Equation (10.16).
For the second iteration the 2CX values are taken as the average of the
XC'X and C7(i) values as given by Equation (10.20) and the iterative proce-
dure is continued. Table: Example 10.7 shows the first two iterations.
The iterative procedure is continued and the answers given in Example
8.6 are obtained.
Example 10.8
Solve Example 8.7 (Example 8.12) by the linear theory method based on
H equations.
294 LINEAR THEORY METHOD
Solution
5 0-4
100-00m
0-2"3 50 g 0-6
Using the known head of 100 m at node 1 and the set pressure head of 80
m at d, the downstream of the PRV, the linearized H equations for nodes
2 . . . 6 of the network for the t'h iteration are, respectively,
Ci — Ci 0 0
c +a + a -Ci -C'd.
-c: ci + a 0
-Ci a + a + ci -Ci 0
0
-a0 -Ci
0
-ci a + Ci
0 0 a +c
ft 100C - 0.3
ft • -80CJ-6 - 0.2
X ft ,= -0.5 (6)
ft -0.4
80C'-6 - 0.6 .
determined and solving Equation (6) we get the values of H2 <, ( . . . H6{1).
The discharge in pipe 7 is evaluated and the pressure head at the down-
stream end of the PRV is determined to check the mode of the PRV. For
the given nodal flows and the pressure head setting of the PRV, the PRV
remains in the operative mode throughout the iterative procedure and the
final solution as given in Example 8.7 (Example 8.12) is obtained.
At the pressure head setting of 89 m, the PRV becomes inoperative dur-
ing the iterative procedure. At this stage, CJ. 6 is replaced by C 7 ', corre-
sponding to the increased resistance from R7 to R7 + /?PRV, and the
known pressure head setting of 89 m is replaced by unknown pressure
head H3. Equation (6) therefore becomes
:» + ci •Ci -C2' 0 0
Ci a +-Cia + a -Ci 0 -Ci
Ci a +-Cia + a -Ci 0
0 0 ci + a -Ci
0 -Ci 0 -Ci ci + a
H2 ' 100C,' - 0.3
H3 -0.2
X Ht -0.5 (7)
H5 -0.4
u -0.6
or
= g". Q3
~3" (10.29)
from which
AGLT =
AQLT ,= n • A0 NR (10.33)
2<2=IQ+^AQLT (10.34)
or
( 1 0 .36)
z* - 2
,Q, = 2 (10.37)
Let the head Ht be known and fixed and the assumed head at nodey for
the first iteration be tHj with correction A//, so that
or
- lHJ) - AH^'"
(10.41)
(10.42)
Therefore, we have
:F M - it,.,l - U0.44,
or
These HJ{1)m values are used to calculate the modified values, i.e.,
hx{i)m, Qx{nm and Cx'(1)m values, so that
lC; + C;
2 C = 2 """ (10.50)
_ [ / / • ( , / / , + 2AHj)]
*-x(l)m — DO.5 (IV.
or
R'o
300 LINEAR THEORY METHOD
Substituting the values of ,C I ' and CV(1)m from Equations (10.51) and
(10.52), respectively, in Equation (10.20), we get
s
(H, - tHj)-°- + [(//, - 1HJ) -
C - - (10.54)
Expanding the term within brackets in the Taylor's series and consider-
ing the first two terms, we get
_ 1
2
* 2
(10.55)
or
(10.56)
Now, the term on the right-hand side of Equation (10.56) represents the
first two terms of expansion in Taylor's series of
Thus, for the /'* iteration, instead of using the average of the assumed
and modified obtained values in the previous iteration, the obtained value
in the previous iteration can be used. Therefore, for the fh iteration,
rC; according to
H3 H4 R* «s H3 HA RA H5
Iteration m m m m m m
10.5 PROBLEMS
11.1 INTRODUCTION
4 0
(a) . (b)
Figure 11.1 Networks: (a) hydraulic network; (b) analogous electrical network.
V = RI (11.1)
Camp and Hazen [143] were perhaps the first to introduce electrical ana-
lyzers for the analysis of flow through hydraulic networks. The analyzers
used linear resistors (V = RI) to represent the pipes of the hydraulic net-
works. To obtain the analogy of the behaviour of hydraulic networks in
these analyzers, the resistances of the linear resistors are successively ad-
justed by trial and error so that ultimately the resistances of these resistors
correspond to RQ"'1 for all pipes. Thus, the network analysis through ana-
lyzers with linear resistors is similar to network analysis through linear
theory. The linear theory linearizes and solves the nonlinear equations
mathematically, while an electrical analyzer with linear resistors linear-
izes and solves the nonlinear equations electrically. Because of the neces-
sity of successive adjustments, these analyzers give the result indirectly.
Therefore electrical analyzers with linear resistors are also termed as in-
direct-reading electrical analyzers.
Although satisfactory results can be obtained with these analyzers for
networks having a large number of nodes and pipes, these analyzers suffer
from the following limitations:
(1) The analyzers involve a trial-and-error procedure.
(2) There is no guarantee of balancing and therefore of convergence.
(3) Because each pipe is represented by different resistors at various
stages of balancing, the resistances of the resistors keep on changing.
No resistor can be assigned a definite resistance. Thus, the storage
and proper selection of the resistors become difficult.
306 OTHER ANALYSIS METHODS
(4) The first settings of the linear resistors give preliminary values of
currents and voltages. These values must be recorded and adjusted
for all the resistors of the network, even though the discharge and
head loss in only a few pipes may be required.
The analyzer of Camp and Hazen [143] used dc supply but later Camp
[135] alternatively suggested the use of ac supply. Reid and Wolfenson
[144] have described the general features of these analyzers. Perry et al.
[145], and Suryaprakasam et al. [146] have offered valuable techniques for
speedy convergence. Application of these analyzers to gas networks has
been described by Clennon [147].
V = Rl" (11.2)
The resistivity of the tungsten filament depends on its temperature, and the
distribution of the temperature along the filament is affected by radiation
of heat to the surroundings and conduction of heat to the leads. As the cur-
rent through a filament increases, the temperature rises which in turn in-
creases the resistance. This causes the voltage to drop as the desired-
nonlinear function of the current. Fluistors resemble incandescent lamps
but are primarily designed to give nonlinear resistance characteristics
rather than to serve as an efficient source of light.
Tungsten filaments provide nonlinear characteristic up to their melting
point reached at about 12 volts. However, the normal working range is kept
between 0 to 5 volts for longer life. Test results have shown that the ob-
tained nonlinear characteristic of the fluistor matches excellently with the
required nonlinear characteristic except for voltages lower than one. How-
ever, these lower voltages would cause an overall error of about 0.5 percent
for the entire network.
The average voltage drop per fluistor is kept at about 2.5 V per fluistor
in the normal operation of the analyzers. Therefore, a hydraulic network
having 40 pipelines in series as in a 20 X 20 grid network can be
simulated by an electrical network having approximately 100 V across its
most widely separated terminals. Networks representing 120 pipes in
series can operate on 300 V without resorting to ac current and isolating
transformers.
Fluistors having approximately 200 different coefficient values ranging
from 0.058 to 500 (a ratio of roughly 1:9,000) would serve as a normal
stock to simulate hydraulic networks. Fluistors are designed so that resis-
tances of successive fluistors are spaced at approximately 5 percent, i.e.,
the resistance of each successive fluistor is approximately 1.05 times that
of the preceding lower one. This practice enables the selection of a fluistor
within an accuracy of about 2.5 percent in its resistance. As the resistance
of a fluistor resembles the resistance constant of a pipe which in turn de-
pends upon its length, diameter, and roughness characteristic, a single
fluistor can represent pipes of different lengths, diameters, and roughness
characteristics, but having the same value for their resistance constants.
For a given voltage drop across the filament, the temperature at the mid-
point of the filament is approximately the same for all fluistors irrespective
of the value of the actual resistance of the fluistor. Therefore, the bright-
ness of the filament can serve as an indication of the voltage drop across
the filament and thereby the head loss in the corresponding pipe. Thus, the
operator, simply by observing the filaments, can know the portions of the
network having large head losses. He may measure only these voltage
drops, estimate corresponding head losses, and suggest suitable changes,
if necessary, in the network.
V. Raman [148] has successfully used small commercial lamp bulbs and
308 OTHER ANALYSIS METHODS
auto lamp bulbs as nonlinear resistors for small networks and provided a
description of the fabrication and performance of such analyzers.
The analyzers with nonlinear resistors are analogous to hydraulic net-
works in their nonlinear behavior and can directly give the required results.
They are therefore called direct-reading electrical analyzers. Because they
were first suggested by McDroy, they are also known as Mcllroy's electrical
analyzers.
Mcllroy has described in detail several aspects such as the preparation
of the network map, selection of voltage-head loss ratios, current-flow rate
scales, operation of the analyzers, and interpretation and accuracy of the
analyzer results.
Network analysis through electrical analyzers does not involve any com-
putational work, provides immediate visual indication of the sections of a
network where head losses are excessive, and also provides ease in consid-
ering alternate designs and different operating conditions. The Hardy
Cross method was the only systematic method prevalent for network anal-
ysis when electrical analyzers were developed in the 1940s, and therefore
electrical analyzers became quite popular. According to McPheron [133],
eighteen Mcllroy electrical analyzers were in use in 1961 in the U.S.A.
These included those at Cornell, Washington, and other universities, used
for research and consultancy work; at Baltimore and Philadelphia, used
for water distribution networks; at southern California and other gas com-
panies, used for gas utility networks; and at the U.S. Bureau of Mines,
used for mine ventilation studies. Several investigators have described the
performances of these electrical analyzers [149-152].
The electrical analyzers require a special setup and thus have high initial
cost. They also require the services of trained personnel. Their accuracy
is within ± 3 percent of the total head loss in a network. Therefore, when
digital computers became readily available, because they provided faster,
cheaper, and more accurate results and could also consider elements other
than pipes in network analysis, they replaced electrical analyzers. Because
of the further advancement in computer technology and also the introduc-
tion of Newton-Raphson and linear theory methods as additional network-
analysis tools, electrical analyzers are almost obsolete as far as the analy-
sis of flow in water distribution networks is concerned.
Fox and Keech [153] considered the steady state of flow in a distribution
network as an ultimate stage of the unsteady flow developed in it after the
start-up of the network. Therefore, they suggested an approach based on
Analysis through Unsteady Behaviour during Start-Up 309
unsteady flow equations for the steady state analysis of distribution net-
works. Usually, the methods of solving unsteady flow problems are more
complex and take more time as compared to steady state problems. How-
ever, as only the ultimate steady state behaviour is required and not the in-
termediate transient states, the approach can be simplified.
The basic unsteady flow equations for flow of water in a pipeline of neg-
ligible slope are [153,154]
dH aV dV fV V|
_L
T"
1/
V + = 0 (11.3)
dx ~~dx~ ~dl ID
and
dh r cl dv
dx + dt
-o (11.4)
g ox
in which V = velocity of flow; H = piezometric head; g — gravitational
acceleration; D = pipe diameter; and c = wave speed. Since we are not
interested in the transient behaviour at intermediate stages, conversion of
Equations (11.3) and (11.4) into finite difference form yields two character-
istic equations for the ultimate steady state flow in a pipe. These two equa-
tions are [153]
c(Vd - Vu) fLVu | V,
Hd = //„ — (11.5)
g 2gD
u = Hd Q* I"'1 01.8)
In the characteristic equations, the wave speed c depends upon the fluid
properties, pipe material, pipe thickness, and pipe anchorage conditions.
The actual value of c would be necessary in an unsteady flow analysis.
However, for the ultimate steady state condition, any value of c would give
the same result, and therefore c can be taken as a round value of 1,000 m/s
for water.
When a reservoir is located at the upstream end or the downstream end
of a pipe, the appropriate characteristic must be used. Consider, for exam-
ple, the three-reservoir network shown in Figure 11.3 in which three reser-
voirs 1, 2, and 3, having water levels //„,, Ho2, and Ho3 are connected to a
junction J through pipes, 1, 2, and 3, respectively. For pipe 1, the forward
characteristic gives
Ac
Hj = Hol — (11.9)
Ac
j = Ho2 + 2» - Q2d) + Q2d I""1 (11.10)
Ac
j = Ho3 + 3» - Q3d) 6 M I"'1 (11.11)
Example 11.1
Solution
120 - Hj HJ - 92 Hj - 80
1442.1 ^" ~ 3244.7 •ana ^ - ~ 1442.1
Substituting these values in the continuity equation for node J, i.e., Equa-
tion (11.12), and solving it we get Hj = 98.55 m. Therefore, Qld = 0.01488
m 3 /s, ftu = 0.00202 m3/s, and Q3u = 0.01286 m3/s. These values are sub-
stituted in Equations (1) . . . (3), the equations are resolved and new val-
ues of Hj, Q,d, ftu, and ftu are obtained. The iterative procedure is con-
tinued until the discharges in all pipes at their upstream and downstream
ends are nearly equal and the steady state condition is reached. After sev-
eral iterations, the answers given in Example 5.2 can be obtained.
312 OTHER ANALYSIS METHODS
11.4.1 INTRODUCTION
Ex = KQxhx (11.13)
SE
= Ti I KQ'h' I (U-14)
The discharges in pipes, i.e., the Qx values, should be such that they
satisfy the node-flow continuity relationship at all demand nodes. There-
fore, the optimization model for a single-source network can now be ex-
pressed as
Q* = q*-Q* (11.19b)
+ Ro3\Q3\"+i (11.20)
= 0 = (n + l)[Rol | qol + Q
[It can be shown that d2SE/dQl > 0, Equation (11.20) is strictly convex
and therefore, Equation (11.21) gives the value of Q3 for the minimum value
of SE.]
Reintroducing Qt and Q2, from Equations (11.19a) and (11.19b) in Equa-
tion (11.21) and simplifying, we get
which is precisely the loop-head loss relationship giving the algebraic sum
of the head losses in the three pipes forming the loop as zero. Thus, when
the expenditure of the system energy for a loop is minimum, the pipe dis-
charges are adjusted so that the loop-head loss relationship is satisfied. As
the node-flow continuity relationships are already satisified and as the
loop-head loss relationship is also now satisfied, the solution of the op-
timization model yields the hydraulic network analysis.
The proof can be extended to single-source, multiple-looped networks.
Selecting one pipe discharge in each basic loop as a basic decision varia-
ble, the optimization model can be recast to an unconstrained optimization
model with the number of basic decision variables equal to the number of
basic loops. By partially differentiating with respect to each basic decision
316 OTHER ANALYSIS METHODS
variable, and simplifying, we get the loop-head loss equations for the
loops. Thus, the solution of the optimization model having pipe discharges
as the decision variables yields the hydraulic network analysis.
The proof can also be extended to multiple-source networks. Herein also
one pipe discharge is selected as a basic decision variable from each real
loop, and one source-node flow as a basic decision variable from each
pseudo loop. The optimization model is then recast as an unconstrained
optimization model and is partially differentiated and simplified to get the
loop-head loss equations as a necessary condition for the optimal solution.
Example 11.2
= O-8rn/s
0-9
Solution
+ 6 | Q4 (1)
04 + 05 - 1.3 = 0 (3)
Q2 - Qs - Qs - 0.9 = 0 (4)
0 , = 2.1 - Q3 - Qs (5)
0 2 = 0.9 + 0 3 + Qs (6)
Minimize: SE = 5 | 2.1 - Q3 - Q5 | 2 8 5 2
+ 8 I 0.9 + 03 + 05
+ 8 I 03 + 6 I 1.3 - 05
+ 5 I 05 (8)
Figure 11.6 Network with ground node and pseudo pipes for the content model.
Analysis through Optimization 319
n+\ 1 «+l
Minimize: Content C = Ro G. R02 G2
1G3 n+l 1
R«\ G4
Gs "+1 + (n +
/ TJ TJ
1)^2 {noa — no 2 )
- G. - G2 = 0 (ii .23b)
Herein the first five terms of the objective function, Equation (11.23a)
represent the energy loss (in discharge-head loss units) in real pipes 1 . . .
5 of the network, respectively, and the last two terms represent (n + 1)
times the energy loss in the pseudo pipes. Note that the node-flow conti-
nuity constraints are written for all nodes of the network.
To prove that the solution of this optimization model gives the network
analysis, consider Qu Qs, and q2 as the basic decision variables. Using the
node-flow continuity constraints, let us express G2, Qi, Q*, and #, in terms
°f Qu Qs, and q2 so that
G4 = ? 2 - G s (11.24c)
Substituting the values of Q2, Q3, Q4, and q, from Equations (11.24a)
320 OTHER ANALYSIS METHODS
dC
= (n + \)[-Ro3 - Qs -
Gl Ro
Equations (11.27a), (11.27b) and (11.27d) are precisely the loop head loss
equations for the two real loops 1-3-4-1, 3-2-4-3, and the pseudo loop 1-4-
2-G-l, respectively. Since all the node-flow continuity and loop-head loss
equations are satisfied, the solution of the optimization model gives the
network analysis. For simplicity, H^ can be taken as zero, so that the gen-
eral content optimization model can be expressed as
- (n + 1) E<ijH
j
OJ (j e source nodes)
(11.28a)
For the network of Figure 11.5, Rol =2, Rol = 1.2, Ro3 = 3, Ro4 =
1.6, Ro5 = 2.4, qo3 = 6 L/s, q^ = 4 L/s, //„, = 100.00 m, and Hol =
95.00 m. The head loss in a pipe is given by h = RQl 85 in which h is in
meters, and Q in liters per second. Analyze the network using the content
model approach.
Solution
0 . - Q3 + 04 - 6 = 0 (4)
and 0 2 + Q3 + 0 5 - 4 = 0 (5)
The solution of this nonlinear optimization model would give the mini-
mum value of the content C. Alternatively, taking 0 , , 0 5 , and q2 as the
basic decision variables, and expressing Q2, Q3, 0 4 , and qx in terms of Qu
Qs, and q2, we get the unconstrained optimization model given by
Minimize: Content C = 2 \ 0.
Q, 2.8S _|_
1.2 10 - 0. - 92
3 a + < ? 2 " G. - 6
2 ,
1.6 <?2 - 0s
2.8S
+ 2. 4 c ) 1 2.85
03 = - 0 . 8 6 3 L/s, H3 = 87.00 m,
05 = 1.598 L/s,
The content model approach described in the earlier section was based
on equations having pipe discharges and nodal flows as basic unknowns,
i.e., on Q-q equations. Colins et al. [75] have also proposed a model,
termed the co-content model, that is based on equations having the un-
known nodal heads as the basic unknowns, i.e., based on //equations. The
Analysis through Optimization 323
unknown pipe flows are expressed in terms of the nodal heads and the
known pipe resistances, so that the energy loss in pipe x is given by
Ex = Q , - h x = 1V
(11.29)
ox
Minimize:
H.t - H3 Hol - H4
Co-Content (CC) -
Rir
H3 - H4 Ho2 - H3
U _
°2
The first five terms of the objective function [Equation (11.30)] represent
Figure 11.7 Network with ground node and pseudo pipes for the co-content model.
324 OTHER ANALYSIS METHODS
the energy loss in real pipes 1 . . . 5 of the network, respectively, and the
last two terms represent [(1/n) + 1J times the energy loss in the pseudo
pipes. Note that there are no constraints and therefore we get directly an
unconstrained optimization model in two basic decision variables H2 and
H3. Therefore, for minimization of the optimization model, partially dif-
ferentiating Equation (11.30) with respect to H2 and H3, we get
+ — H3 H3 - H4
t£-(i '"- oi
— H3
Ro + J =o (11.31a)
dCC
dK HHl- — HA H3 — H4
Ho2 — H4
= 0
(11.31b)
which are precisely the node-flow continuity equations for nodes 3 and 4,
respectively. Therefore, the solution of the co-content optimization model
gives the values of the unknown nodal heads H3 and H4 such that the node-
flow continuity relationships are satisfied. When the nodal heads are the
basic decision variables, the loop-head loss constraints are automatically
satisfied. Further, as the solution now also satisfies the node-flow conti-
nuity relationships, the solution of the co-content optimization model
gives the network analysis. For simplicity, Hoc can be taken as zero so that
the general co-content optimization model can be expressed as
Minimize:
Co-Content CC -
ft -
+ - +
(11.33)
Problems 325
in which //, = Hoi and Ht = Hoj for source nodes. Note that the second
summation is only for nodes with known flows.
Example 11.4
Solution
Taking // oC = 0, n = 1.85, and substituting the known values for the net-
work of Figure 11.7, the co-content model is expressed as
Minimize:
Co-Content CC =
100 - ft I ' 5 4 100 - fl, I
(1.2) 0 5 4
(2)0.54
95 - H3
LJ u I 1.54 (1.6)° 54
(3)054
1.54(67/3
54
95 - H4 | '
(24)0.54
The solution of the optimization H the minimum value of the
model gives
co-content as 1455.299 Ls"1- m for H3 = 87.00 m, and H4 = 89.29 m.
Other answers are as given in Example 11.3.
The number of the unknowns in the general content model (Example
11.3) is seven: five for the real pipes and two for the pseudo pipes. In addi-
tion, it has four constraints, one for each node. Thus, the number of un-
known decision parameters is maximum as in the formulation of Q equa-
tions in network analysis. When the content model is converted to an
unconstrained one, the number of decision variables reduces to three, one
for each loop, as in the formulation of AQ equations. The co-content
model is an unconstrained one and has only two decision variables for two
unknown nodal heads as in the formulation of H equations.
11.5 PROBLEMS
Problem 114: Solve Problem 5.2 using the system energy approach.
Problem 11.5: Solve Problem 5.2 using the content model approach
and find the minimum content of the network.
Problem 11.6: Using the solution of Example 8.2, find the minimum
content of the network, and, by slightly changing the solution, verify that
the content of the network is minimum for the obtained solution.
Problem 11.7: Solve Problem 5.2 using the co-content model approach
and find the minimum co-content of the model.
Problem 11.8: Using the solution of Example 8.3, find the minimum
co-content of the network, and, by slightly changing the solution, verify
that the co-content of the network is minimum for the obtained solution.
Dynamic Analysis
12.1 INTRODUCTION
Type of Demand:
Residential
Commercial
Office Buildings
Hospitals
00 —
0 2 6 8 10 12 2 L, 6 8 10 12
Midnight A.M. Noon RM. Midnight
objectives can be achieved by carrying out the analysis of the network over
a period of 24-48 hours under changing demand patterns, reservoir water
levels, and boundary conditions. Such an analysis of the network is called
an extended period simulation, extended period analysis, or simply dy-
namic analysis.
The necessity of dynamic analysis for water distribution networks was
recognised by Shamir and Howard [88], Tart [159], and Demoyer, Gilman,
and Goodman [160]. However, Rao and Bree [161] and Rao, Markel, and
Bree [158] were perhaps the first to suggest a general, systematic readily
available procedure for carrying out a dynamic analysis of water distribu-
tion networks. Their approach is described in the next section and is fol-
lowed, in a subsequent section, by an improvement suggested by Bhave
[162].
12.2.1 INTRODUCTION
12.2.2 PROCEDURE
t + At are used to calculate the volume change in the reservoirs and the
predictor iteration is repeated. The predictor-corrector iterative procedure
is continued until the predicted and corrected reservoir water levels agree
to the desired degree of accuracy. Thus, the static analysis of the network
from time t is linked to that for time t + At through an iterative, pre-
dictor-corrector integration step, thereby giving the dynamic analysis of
the network for time interval (t, t + At). The schedule of the pump and
valve settings and the effect of the boundary conditions, if any, are also
considered during this step. The reservoir flow rates and water levels for
time t + At obtained at the end of the iterative integration step become the
input for the next time interval.
The procedure for carrying out the dynamic analysis of a network by the
indirect method consists of the following steps:
Step 1. At time t, the following data are available:
a. Reservoir HGL values Hr(t) and reservoir water volumes Vr{t) in
which r 6 M, the set of reservoirs
b. Nodal demands <£,(?) and qj(t + At) in which y E N, the set of demand
nodes
c. If the network contains boundary elements such as pumps that bring
water into the network from external sources and booster pumps that
bring water into the network from neighboring pressure zones, then the
flows qb(t) or the HGL values Hb(t) and the relationship between qb and
Hb at all the boundary elements, b, b G B, the set of boundary ele-
ments. The flows and the HGL values at the boundary elements may be
a function of the network demand and of the water purchase contracts.
Step 2. Using the data in step 1, the static analysis of the network is ob-
tained for time t. This solution provides the HGL values Hj(t),j <E N for
the demand nodes and the outflow rates qr(i), r € Mfor the reservoirs.
Step 3. Assuming the outflow rate qr(t) for reservoir r to be constant
during the time interval (t, t + At), the volume decrease in reservoir water
is obtained as
AVr(t, t + At) = qr(t)At, r = 1 . . . M (12.1)
Step 4. The net total outflow from all the reservoirs is computed as
M
Similarly, the net outflow from all the boundary elements is computed as
B
AV
» = £ ?»«A* (12-3)
Iterative Method 331
As the total outflow from the reservoirs and boundary elements is the in-
flow to the distribution network, the total inflow to the network is LAVr +
LAVb.
Step 5. By plotting a curve of the demand rate versus time for the entire
region under consideration, the cumulative total demand D(t, t + At) is
obtained from the area under the demand curve for time interval (t,
t + At).
Step 6. The error in water volume balance for the network is obtained as
rEM
At
[qb(t) + qb(t + At)] Y + D(t, t + At)
bEB
(12.7)
332 DYNAMIC ANALYSIS
in which Ec = corrected error for the entire network. This error is reallo-
cated to reservoir r, r G M, in proportion to the average flow rate. Thus,
+ + At)
eM, t + At) = - ^ *<' X Ec(t, t + AO
L '(f + )
(12.8)
Using this corrected reservoir volume, the HGL for reservoir r for time
t + At, i.e., Hr(t + At), is computed for all r, r E M, from the respec-
tive capacity-elevation relationship Vr = fr(Hr), r G M.
Step 12. The difference between the predicted and corrected HGL values
is estimated for all r, r G M. If the difference is more than the permissible
limit, the predictor-corrector integration step is repeated for the same time
interval. If the difference is within the permissible limit for all reservoirs,
the dynamic analysis for time interval (t, t + At) is complete. Finally, the
obtained corrected values Hrc(t + At) and Vrc(t + At) are set to Hr(t +
At) and Vr(t + At) respectively, for all reservoirs. These values then serve
as the starting values for the next time interval.
Step 13. If there is a switch of a pump or valve in the time interval (t,
t + At), then the time t + A't (t < t + A't < t + At) at which the
switching occurs is determined. The time interval is reduced from At to
A't and the integration procedure is carried out for two time intervals
(f, t + A't) and (t + A't, t + At).
The entire procedure is cycled through the various time intervals to
cover the entire period of analysis. This completes the dynamic analysis of
the network.
Example 12.1
The network shown in Figure: Example 12.1 has nodes 1 and 2 as reser-
voirs and nodes 3 . . . 6 as demand nodes. The resistance constants in the
Hazen-Williams head loss formula, hf = RQ1852 for pipes 1 . . . 7, are
shown underlined in the figure. The elevation-capacity data for the reser-
voirs are given in Table 1: Example 12.1, while the nodal demand rates at
an interval of 4 hr are given in Table 2: Example 12.1. From the demand
Iterative Method 333
Node,1
I rj
curve for the entire network, total network demand volumes are obtained
and given in Table 3: Example 12.1. Reservoir 1 is filled from an external
source during time interval 4-12 hr at a constant rate of 14 m 3 /min, while
reservoir 2 is filled during time interval 12-20 hr at a constant rate of 10
m 3 /min. At t = 0 hr, the water levels in reservoirs 1 and 2 are 101.200 m,
and 102.000 m, respectively. Both reservoirs are floating on the network.
Carry out the dynamic analysis of the network from 0 hr to 4 hr taking
it as one time step. (Complete dynamic analysis for the total period of
0-24 hr is carried out in Example 12.4.)
Solution
0 100.00 98.00
1,000 100.60 98.80
2,000 101.20 99.60
3,000 101.80 100.40
4,000 102.40 101.20
5,000 103.00 102.00
6,000 103.60 102.80
334 DYNAMIC ANALYSIS
Total Network
Time Interval Demand Volume
in Hours in Cubic Meters
0-4 1,488
4-8 2,256
8-12 2,304
12-16 2.184
16-20 1.992
20-24 1,296
Iterative Method 335
1 8
.-. elp(0, 4) = ^ X 672 = 115.60 m3
and
and
i« no
Hip(4) = 101.20 + y ~ X (101.80 - 101.20) = 101.354 m
and
2
#2P(4) = 100.40 + . 5 ^ 2 X (101.20 - 100.40) = 100.605 m
1 ,UUU
240
X - y - + 1,488 = 96 m3
= 29.99 m3 (outflow)
and
-3.9849 - 3.9914
X
e2c(0, 4) - (Q 5 g 4 9 _ 3 9 g 4 9 ) + ( _ 4 2086 - 3.9914)
= 66.01 m3 (outflow)
Step 11. The reservoir volumes are corrected using the corrector
equation:
M4) = ^.(0) + fa,(0) + 9l(4)] X y + elc(0, 4)
= 1,535.17 m3
and
K2c(4) = V2(0) + fo2(0) + q2(4)] X y + e2c(0, 4)
= 3,976.83 m3
Direct Method 337
and
As the absolute values of the difference between the predicted and cor-
rected water levels for the reservoirs are more than 0.001 m, an assumed
permissible value, the corrected water levels are taken as the predicted
values and the next predictor-corrector integration iteration (steps 9 . . .
12) is carried out. The iterative procedure (steps 9 . . . 12) is continued
and the iteration results are shown in Table 4: Example 12.1. The solution
of the dynamic analysis for time interval 0-4 hr is <2i(4) = 0.9525
nrVmin, Q2(4) = 2.1981 nrVmin, Q3{4) = -1.0474 nrVmin, QA(4) =
1.3983 nrVmin, Qs(4) = 1.9962 m3/min, Q6(4) = 2.6036 m3/ min,
G,(4) = 1.0037 nWmin, q,{4) = -3.1506 nrVmin, q2(4) = -5.0494
nrVmin, K,(4) = 1671.7 m3, V2{4) = 3841.2 m \ //,(4) = 101.003 m,
#2(4) = 101.073 m, H3(4) = 100.637 m, H4(4) = 100.143 m, // 5 (4) =
98.704 m, and H6(4) = 99.308 m.
tor integration iteration (steps 3 . . .8). These values are then successively
corrected (steps 9 . . . 12) so that the predicted reservoir HGL values Hrp
(t + At) and the predicted reservoir flows qn,(t + At) reasonably agree
with the corresponding corrected values Hrc(t + At) and qrc(t + At), re-
spectively, for all reservoirs. However, as has been recently shown by
Bhave [162], it is possible to obtain a relationship between Hr(t + At) and
qr(t + At) and use it in the analysis so that the static analysis for time
t + At can be obtained directly from the static analysis for time t thereby
avoiding the predictor-corrector iterative procedure.
12.3.1 THEORY
AK ( A
Hr(t + At) = Hr(t) + f !:f, * ° (12.12)
JrlJirU)\
average rate of demand for the entire network. If we multiply these values
by time interval At, we get the total volumes, and therefore,
Equation (12.13) states that the total volume of water supplied by the res-
ervoirs and boundary elements during the time interval (t, t + At) as ob-
tained from the average rate of supply is equal to the total volume of out-
flow from the demand nodes of the network during the same time interval
(t, t + At) as obtained from the average rate of demand.
Now, from the network demand curve, the area under the curve for time
interval (t, t + At) gives the total volume of water consumed and therefore
the total volume of water actually supplied by the demand nodes. This ac-
tual volume is D(t, t + At) instead of
as estimated from the average rate basis. If C(t, t + At) is taken as a cor-
rection factor to convert the estimated total supply to the actual total sup-
ply by the demand nodes, we have
„,
(12.14)
At
2rf qAt) + (t + A?) T
j j
Distributing the difference between the actual and estimated total net-
work demand to all reservoirs and boundary elements in proportion to
their withdrawal rates, the supply of volume of water by reservoir r in time
interval (f, t + At) can be expressed as
As stated earlier, the term/ r '[// r (/)] represents the cross-sectional area
of reservoir r at HGL Hr(t). For reservoirs having constant cross-sectional
area, the term ft[Hr(t)] remains constant throughout the entire dynamic
analysis. For reservoirs having variable cross-sectional area, the value of
fr[Hr(t)] is computed for each starting time t of each time step, and it is
assumed that the value remains constant during the time step, i.e., the time
interval it, t + At). The correction factor C(t, t + At) is the same for all
reservoirs and boundary elements for a particular time interval (f, t + At),
but will be different for different time steps and therefore is computed for
each time step.
All terms except qr(t + At) and Hr(t + At) in Equations (12.16) and
(12.17) are either known or evaluated. Thus, these equations give relation-
ships between reservoir HGL values and reservoir outflow rates at time
t + At.
If a reservoir is filled at a constant rate of qfr during time interval (t,
t + At), Equations (12.16) and (12.17) would modify to, respectively,
and
it + At) - Hr{t)] X MHAt)]
A
'> = ^ Ar
C(f, f + AO y
(12.19)
Direct Method 341
If the supply rate from a boundary element remains constant during time
interval (t, t + At), i.e., if it is qb(f) throughout the time interval (t, t +
At), a correction must not be applied to it. Therefore, while evaluating the
correction factor given by Equation (12.14), this known and fixed volume
qb(t)At is subtracted from the numerator and denominator of the right-
hand side of Equation (12.14).
If the supply to the network is from a clearwell having constant water
level, this known and fixed water level is used in the analysis for time
t + At. [Herein, the question of establishing and using a relationship be-
tween Hr(t + At) and qr{t + At) does not arise.] The value of qr(t + At)
for the clearwell is obtained from the analysis. The supply of water from
the clearwell during the time interval (t, t + At) is obtained using Equa-
tion (12.15).
Equations (12.16). . . (12.24) are used in formulating the required equa-
tions.
12.3.2 PROCEDURE
The procedure for carrying out the dynamic analysis of a network by the
direct method consists of the following steps.
Step 1. Same as step 1 of the indirect method.
Step 2. Same as step 2 of the indirect method.
Step 3. Using the nodal demand rates at times t, and t + At, determine
342 DYNAMIC ANALYSIS
the total outflow from the network during the time interval (f, t + At), as
obtained from the average rate of demand.
Step 4. Same as step 5 of the indirect method.
Step 5. Compute the value of the correction factor C(t, t + At) using
Equation (12.14).
Step 6. Formulate the appropriate equations from Equations (12.16) . . .
(12.24).
Step 7. Using the equations obtained in step 6 and other equations, as
described in Chapter 7, formulate the equations required for carrying out
the static analysis for time t + At.
Step 8. Carry out the static analysis for time t + At.
Step 9. The system is checked for preset switch points for controlling
valves and pumps. If one of these is switched in the time interval (t, t +
At), then go to step 10; otherwise, the dynamic analysis for time step (t +
At) is complete and go to step 3 to continue the dynamic analysis for the
next time step.
Step 10. If there is a switch of a pump or a valve in the time interval (t,
t + At) then determine the time t + A't (t < t + A't < t + Af) at which
the switching occurs. Reduce the time interval from At to A't and carry
out the dynamic analysis (steps 3 . . . 9) for time interval (t, t + A't) and
then for time interval (t + A't, t + At).
The direct method for carrying out dynamic analysis is illustrated in Ex-
ample 12.2.
Example 12.2
Carry out the dynamic analysis of the network of Example 12.1 by the
direct method for time interval (0, 4 hr), taking it as one time step.
Solution
and
and
As seen from Figure 1: Example 12.1, for the assumed directions of flow
in pipes, the rate of flow (inflow positive, outflow negative) for the reser-
voirs is given by
and
Steps 7 and 8. Equations (1) . . . (6) are used in formulating the equa-
tions necessary for carrying out static analysis at t = 4 hr. Replacing
1.852, the index in the pipe-head loss relationship by n, for simplicity, we
can derive the necessary equations as follows:
HEquations. The //equations are formulated by writing node-flow con-
tinuity equations at all nodes of the network. Therefore, using Equation
(2), the H equation for node 1 is given by
0.4 | [ 0.2
= 0 (7a)
344 DYNAMIC ANALYSIS
= 0 (7b)
[Note that q, and q2 represent flows for reservoirs. Because they are out-
flows, as seen from Equations (5) and (6), respectively, they are preceded
here by a negative sign.]
For nodes 3 . . . 6, respectively,
» (7e,
-1.600-0 (7f,
From the loop-head loss equations for the basic loops I and II, we get,
respectively,
Similarly, from the loop-head loss relationship for the pseudo loop, i.e.,
loop III,
Similarly,
Substituting the values of //,(4) and H2(4) from Equations (8h) and (8i),
respectively, in Equation (8g) and simplifying, we get
only the Q equation for the pseudo loop has changed. The total number of
Q equations for the dynamic analysis has remained the same, i.e., seven,
as in the ordinary static analysis.
AQ Equations. The values of Qi(4) . . . Q7(4) are assumed so that they
satisfy the node-flow continuity equations, Equations (8a) . . . (8d). Let us
denote these values by suffix o as usual, and let us omit 4 and the enclosing
parentheses for simplicity. If AQ, . . . AQm are the loop-flow corrections
(clockwise positive) for loops I . . . I l l , respectively, Equations (8e) . . .
(8g) become, respectively,
- /f,(4) = 0 (9c)
Substituting the values of 7/^4) and H2(4) from Equations (8h) and (8i),
respectively, in Equation (9c), we get
Example 12.3
Obtain the dynamic analysis for the network of Figure: Example 8.6 (re-
produced herein as Figure: Example 12.3) for time interval 0-2 hr. The
nodal demands at t = 0 hr are as given in Example 8.6. The nodal de-
mands at t = 2 hr are q3(2) = 0.12 m 3 /s, qA{2) = 0.10 m3/s, qs(2) = 0.18
nvVs, and q6(2) = 0.10 m3/s. The reservoir at node 2 has a constant cross-
348 DYNAMIC ANALYSIS
Sump Reservoir
sectional area of 1,000 m2. The water level in the sump at node 1 remains
constant at 80 m. The total network-demand volume, as obtained from the
demand curve, during the time interval 0-2 hr is 3640 m3.
Solution
We shall carry out dynamic analysis by the direct method using the lin-
ear theory method based on Q equations.
Steps 1 and 2. The static analysis solution of the network for time t = 0
hr is already obtained in Example 8.6, and is as follows:/^O) = 0.01852;
Rl(0) = 189.0, (2,(0) = 0.1619 m'/s; / 2 (0) = 0.02167, R2(0) = 2238,
Q2(0) = 0.0508 m 3 /s,/ 3 (0) = 0.02325; R3(0) = 10,123, Q3(0) = 0.0311
m 3 /s; / 4 (0) = 0.02322, R4(0) = 3,791, Q4(0) = 0.0324 nrVs; / 5 (0) =
0.02125, /?5(0) = 823.3, Qs(0) = 0.1016 nv7s,/ 6 (0) = 0.02392, R6(0) =
3905, <26(0) = 0.0565 mVs, / 7 (0) = 0.01776, R7(0) = 139.7, Q7(0) =
0.2781 m 3 /s; and hp(0) = 11.44 m.
Step 3. The total outflow from the network during time interval 0-2 hr
as obtained from the average rate of demand = f(0.44 -I- 0.50)/2] X
7,200 = 3,384 m3.
Step 4. The total network demand volume, as obtained from the demand
curve, during the time interval 0-2 hr = 3,640 m3.
Step 5. Correction factor C(0, 2) = 3,640/3,384 = 1.07565.
Step 6. For the sump at node 1, //,(2) = 80.00 m. For the reservoir
at node 2, as q2(t) = -Q7(i), we get H2(2) = 100.00 + 1.07565 X
[-0.2781 - Q7(2)] x 3600 •* 1,000, i.e.,
and
Substituting the values of H2(2) and hp(2) from Equations 1 and 2, re-
spectively, Equation (9). on simplification, becomes
Example 12.4
Using the data given in Example 12.1, carry out the dynamic analysis of
the network for 0-24 hr taking a time step of 4 hr.
Solution
The dynamic analysis for 0-4 hr has been carried out in Example 12.1
and 12.2, by indirect and direct methods, respectively. Following either of
the two methods for the other time steps, the dynamic analysis for 0-24 hr
can be obtained.
For time steps of 4-8 hr and 8-12 hr, reservoir 1 is filled at a constant
rate of 14 nrVmin. Therefore, Equations (12.18) and (12.19) should be used
instead of Equations (12.16) and (12.17) in formulating the desired equa-
tions. Similarly, reservoir 2 is filled at a constant rate of 10 nWmin for time
intervals of 12-16 hr and 16-20 hr, and therefore, the equations stated for
reservoir 1 should be used.
The complete dynamic analysis solution is given in Table: Example 12.4
[162].
12.4 PROBLEMS
Problem 12A: Using the data given in Example 12.3 and Problem 12.3,
carry out the dynamic analysis for the time interval 0-4 hr treating it as
one time step, and compare it with that obtained in Problem 12.3.
Problem 12.5: Carry out the dynamic analysis of the network de-
scribed in Example 12.1 for 0-24 hr if reservoir 1 is filled for 0-8 hr and
reservoir 2 for 16-24 hr at the given rates.
Node Flow Analysis
13.1 INTRODUCTION
(1) The actual nodal demands are more than the anticipated design de-
mands due to an unanticipated accelerated growth over the entire
area. Thus, the system becomes inadequate even during the design
period.
(2) Even though the overall growth is as anticipated, it is an imbalanced
one causing the system to be deficient in some zones and more than
adequate in others.
(3) The withdrawals at some nodes are more than the design ones, due
to large diameter withdrawal connections, installation of pumps, pil-
ferage of water, or excessive leakage. This naturally would create a
deficiency in some other parts.
(4) The distribution system is used beyond its design period and thus
becomes deficient. (The usual design period of a distribution net-
work is 30 years. However, because the cost of replacement of the
entire network is usually enormous, the distribution network is con-
tinued in use beyond the design period by suitably strengthening it
wherever necessary and replacing only those pipes that have exces-
sive corrosion and leakage, thereby creating maintenance problems.)
(5) Some pipes are closed for cleaning, repairs, or replacement. This
may cause a temporary deficiency at some nodes.
353
354 NODE FLOW ANALYSIS
In all these situations either the nodal demands have changed or the dis-
tribution network has been deficient, and therefore it is unable to meet the
demands. In practice, therefore, it would be strengthened or expanded by
a trial-and-error procedure or by using one of the several available op-
timization techniques [85,86,166-1751. It would, however, be worthwhile
to know how the system would actually behave until the necessary modifi-
cations are carried out.
In the distribution network analysis methods that we studied in the ear-
lier chapters, whenever the nodal demands were specified, we presumed
that the network was capable of supplying them, and therefore it actually
supplied them. The primary aim of such analysis is to determine the heads
at the demand nodes and check that they are at least equal to the minimum
required residual heads, and the network is satisfactory. Therefore, Bhave
[163-165] termed this type of analysis as Node Head Analysis (NHA). The
NHA helps us in locating nodes that are deficient in heads and also in de-
termining the head deficiency at such nodes. We can, therefore, decide on
providing additional heads at source nodes, boosting pressures within the
network, strengthening in certain localities, providing check valves and
other similar modifications.
In the NHA it is presumed that the modifications are available, and
therefore, the actual nodal outflows are the same as the nodal demands.
However, such modifications are not available especially when the change
in nodal demands or the deficiency in the network is temporary. In these
situations the network fails to satisfy the demands at all nodes. It may be
able to satisfy demands at some nodes, partially satisfy demands at some
other nodes, while failing to provide even a drop of water at the remaining
nodes. To analyze the network for these situations, it is necessary to con-
sider the network as it is and determine the actual flows that would be
available at different nodes. Because this type of analysis would help in de-
termining the actual available nodal flows, Bhave [163-165] termed it
Node Flow Analysis (NFA).
Because NHA is a steady state analysis, the nodal demands are fixed and
therefore, even though more water may be available at a node, the with-
drawal remains fixed at the specified value. In NFA also, even though more
water may be available at a node, the withdrawal does not exceed the spec-
ified value; however, it may be less or even zero. Thus, the outflow at a de-
mand node in NFA lies between the specified value (the network is capable
of and therefore actually supplies the specified demand) and zero (the net-
work is unable to provide any water).
Before the general NFA principle is developed, we shall study the be-
haviour of a simple two-demand-node serial network under different con-
ditions.
Two-Node Serial Network 355
Figure 13.1 Elevated service reservoir and two-demand node serial network.
356 NODE FLOW ANALYSIS
Datum
Figure 13.2 Hydraulic grade lines for different positions of downstream node for increase in
upstream nodal flow.
Two-Node Serial Network 357
Datum
Figure 13.3 Hydraulic grade lines for different positions of downstream node for increase in
downstream nodal flow.
B3, qB may be affected by increase in qA but does not become zero, i.e.,
qp > qB > 0.
When B coincides with or lies below B3, qB ( = q'Bs) is not affected at
all by an increase in qA. When qA attains its maximum value, qT*, qB is
still equal to qf*, and SAB3 is the limiting HG line.
HG line is SABS. Thus, for B between B4 and B5, qA is between qACh and
zero.
When B is lower than Bs, for a further increase in qB the HG line goes
below SABS. This stage continues until the HG line becomes SAjB6 where
the maximum discharge carrying capacity of pipe 1 is reached. Thus, for
B between Bs and B6, qB can increase further, although qA remains zero.
(Actually in practice, when the HG line approaches SAjB6, air entrainment
takes place and the flow becomes discontinuous, however, we shall neglect
this aspect.)
When B is lower than B6, because subatmospheric pressure at Aj is not
feasible, there is no further increase in qB.
This two-node serial network is shown in Figure 13.4. In Figure 13.4(a),
the two outlets A and B are nearly at the same level and therefore any in-
crease in qA beyond qAcs may affect qB ( = qBa) and qB may even become
zero if outlet B is above or at the same level as that of A. However, any
feasible increase in qB will not affect qA ( = q^). In Figure 13.4(b), outlet
A is much higher than outlet B. Therefore, any increase in qB beyond ^^,es
may affect qA ( = qAcs) and qA may even become zero. However, any feasi-
ble increase in qA will not affect qB ( = qdBcs).
Water distribution networks usually resemble Figure 13.4(a) and there-
fore outflow at downstream nodes may get affected by an increase in out-
flow at upstream nodes. However, for distribution networks in uneven
terrain, for portions resembling Figure 13.4(b), downstream nodes in low-
lying areas may affect the outflow at upstream high-lying nodes. House
plumbing with overhead storage tanks resembles the situation of Figure
13.4(b).
When both qA and qB increase beyond qAK and qp, respectively, de-
pending upon the increase and the relative positions of A and B, different
(Q) (b)
Figure 13.4 Two-node serial system: (a) outlets A and Bat nearly same levels; (b) outlet A much
higher than outlet B.
Two-Node Serial Network 359
Example 13.1
Solution
For the design demands, the pipe discharges are: Qt = 1.8 -I- 1.2 = 3
nrVmin, and Q2 = 1.2 nrVmin.
.-. ft, = RtQi = 0.6 x ¥ = 5.40 m, and h2 = R2Q1 = 2.5 x
1.22 = 3.60 m.
.-. Hdr = 100.00 - 5.40 = 94.60 m > 94.00 m ( = HA), and Hp =
94.60 - 3.60 = 91.00 m > 90.00 m ( = // B ). Thus, an increase in qA and
qB beyond their design values is feasible.
Case 1. Let AqA be the maximum permissible increase in qA without af-
fecting qBe%. Considering path S4, for outflow to occur at A,
For the maximum value of AqA, the condition given by Inequality (1) is
critical. Therefore, 100.00 - 0.6(3 + AqA)2 = 94.00 gives the maximum
value of AqA = 0.1623 nWmin. Thus, qA can increase to 1.8 + 0.1623 =
1.9623 nWmin for which HAl = 94.00 m, and HBl = 94.00 - 3.60 =
90.40 m > 90.00 m (=HB). For the maximum permissible increase in qA,
because Inequality 1 for node A is critical, an increase in qA is governed
by node A itself. Therefore, an increase in qA cannot affect q^.
Case 2. Let AqB be the maximum permissible increase in qB without af-
fecting qAts. From path S4,
From Inequality (3), AqB = 0.1623 rrrVmin, and from Inequality (4),
AqB = 0.1009 m3/min. Therefore, the maximum permissible AqB =
0.1009 nrVmin, and the maximum permissible value of qB = 1.2 +
0.1009 = 1.3009 nrVmin. HA, = 100.00 - 0.6(3 + 0.1009)2 = 94.23 m
(>94.00 m), while HB, = 90.00 m. As Inequality (4) for node B is criti-
cal, the maximum permissible increase in qB is governed by node B itself
and thus an increase in qB cannot affect qAcs.
Both nodes are independent in this example and the increase in outflow
at one cannot affect the design outflow at another.
Example 13.2
Solution
Example 13.3
Solution
From Inequality (1), AqB = 0.1623 nrVmin, and from Inequality (2),
AqB = 1.880 m3/min. Thus, the maximum permissible increase in qB
without affecting qAa is 0.1623 nrVmin. Therefore, qA = 1.8 nrVmin;
maximum qB = 1.2 + 0.1623 = 1.3623 nrVmin; HAl = 94.00 m; and
HBl = 100 - 0.6(3.1623)2 - 2.5(1.3623)2 = 89.36 m.
Herein, node A is critical and the maximum permissible increase in qB
without affecting qAcs is governed by node A.
Case 2. As seen in case 1, when the withdrawal at B increases beyond
qBCi, HA, decreases and just becomes 94.00 m when qB = 1.3623 nrVmin.
When the withdrawal at B increases further, HAt remains at 94.00 m but
qA decreases. Therefore, a stage is reached when HA, is 94.00 m and qA
just becomes zero. Thus, for this situation, HAt — 94.00 m, qA = 0, and
Qi = Q2 = qB- From path SA, 100.00 - O.602 = 94.00 m, which gives
Qi = 3.1623 nWmin = Q2 = qB. HereL the entire withdrawal of
qAcs = 1.8 m3/min at node A is transferred to node B. HB, = 100.00 -
0.6(3.1623)2 - 2.5(3.1623)2 = 69.00 m. Thus, qA = 0, qB = 3.1623
nrVmin, HAx = 94.00 m, and HB, = 69.00 m.
Case 3. The maximum feasible value of qB is obtained when either
(1) HBl = 62.00 m and HA, < HAj (=92.00 m), for flow to continue
362 NODE FLOW ANALYSIS
Example 13.4
Solution
1. The NHA of the network can be obtained by any one of the methods
described earlier. The NHA solution, thus obtained, is shown in columns
TABLE: EXAMPLE 13.4.
(Values Used for Carrying Out Analysis Are Shown in Parentheses.)
4 and 5 of Table: Example 13.4. The values used for carrying out the NHA
are shown within parentheses.
2. The available excess heads (//,• - Hf'n) at the demand nodes for the
NHA solution obtained earlier are shown in column 6. The least excess
head is 0.55 m, available at node 5.
When the outflow at node 3 increases, the available excess heads at all
nodes will decrease. As node 5 has the minimum available excess head, let
us assume that it would govern the maximum possible outflow at node 3
without affecting the outflows at all other nodes. Therefore, we take Hs =
H%m = 88.80 m as a known value and q3 as unknown. Thus, the un-
knowns are qu H2, q3, H3, H4, and H6. (Note that all network-solvability
rules are satisfied and the NHA solution is feasible.) The NHA solution is
obtained and shown in columns 7 and 8. As //, = Hf" at node 5, and
Hj > Hfm at nodes 2, 3, 4, and 6, the obtained solution is admissible.
Thus, the maximum possible outflow at node 3, without affecting the out-
flows at all other nodes, is 1.353 m3/min.
Example 13.5
Obtain the node-flow analysis of the serial network shown in Figure: Ex-
ample 13.5. The pipe head loss relationship is h = RQl 85 in which h is in
Node.O 1 2 3 4
a > • • • » m * *
Pipe, 1 2 3 4
meters and Q in cubic meters per minute. The source HGL is 100.00 m,
and the data for pipes and demand nodes are given in rows 1-3 of Table:
Example 13.5.
Solution
The complete NFA solution is given in Table: Example 13.5. Rows 1-3
show the given data. The corresponding NHA solution is given in rows
4-7 in which assumed values are shown in parentheses. As //, < Hfm at
nodes 2, 3, and 4, the solution is inadmissible and the NFA solution is nec-
essary.
After several trial-and-error iterations, the finally obtained admissible
solution is shown in rows 8-11^ in which the assumed values are shown in
parentheses.
Herein, the NFA solution is obtained by trial and error. The systematic
procedure for carrying out NFA is developed in the following pages.
TABLE: EXAMPLE 13.5.
(Values Used in Carrying Out Analysis Are Shown in Parentheses.)
Row Node Pipe Node Pipe Node Pipe Node Pipe Node
No. Parameter 0 1 1 2 2 3 3 4 4
1. Critical Node (CRIT). When the available HGL at a node, //, is equal
to the minimum required one, Hf" (//, = Hfn), the node is a critical
node.
2. Supercritical Node (SUPC). When the available HGL at a node is
more than the minimum required one (//,• > H]"n), the node is a super-
critical node.
3. Subcritical Node (SUBC). When the available HGL at a node is less
than the minimum required one (//,• < Hjm), the node is a subcritical
node.
The category of a node depends upon its HGL value and the flow. There-
fore, the nature of a node can be completely described by combining both
these categories. However, in practice, some HGL-dependent and flow-
dependent categories are not compatible with each other. For example, a
supercritical node cannot be a partial-flow node or a no-flow node. Simi-
larly, an adequate-flow node cannot be a subcritical node. Thus, the two
categories of a node can be combined only if they are compatible with
each other. Such compatibility of node categories is given in Table 13.1.
Note that a zero-flow node is compatible with any of the HGL-dependent
qj = 0, if Hj < Hf (13.3c)
The loop-head loss relationship must be satisfied for all loops of the net-
work. Therefore,
Subject to:
or
or
the total number of the nodes. In the NFA problem described herein, the
resistance constants of all pipes are known. At all source nodes, the nodal
heads are known but the nodal flows are unknown, the number of these un-
knowns being M. In the case of demand nodes, except for zero-flow de-
mand nodes (qp = 0), the nodal flow will lie between qf* and 0 depend-
ing upon which of the alternate constraints prevails. The number of these
unknowns, therefore, must lie between 1 (the network must have at least
one nonzero-flow demand node) and N (all demand nodes are nonzero-
flow demand nodes). Further, at all the Ndemand nodes, the HGL values
are also unknown and, therefore, the number of these unknowns will be N.
Thus, the total number of the unknowns for the entire network lies be-
tween M + N + 1 and M + IN, while the total number of the nodes and
therefore the total number of the available equations is M + N.
In order to bring down the total number of unknowns to M + N and
make NHA feasible, we shall treat either the demand or the HGL as
known at each demand node. Thus, if we choose qj = q'p as given by
constraint (13.6c) or q} = 0 as given by constraint (13.6e), the nodal flow
will be treated as known and the nodal head as unknown. On the other
hand, if we choose the constraint (13.6d), we shall take //, = Hf" and
therefore the nodal head will be treated as known and the nodal flow as
unknown. Thus, the total number of unknowns at any instant will be M +
N, equal to the number of nodes in the network.
Example 13.6
Solution
Assumed Convert to
Node Category Obtained Node Category Node Category
0) (2) (3)
Row Node Pipe Node Pipe Node Pipe Node Pipe Node
No. Item 0 1 1 2 2 3 3 4 4
Data
0.05 0.10 0.20 0.30
m3/min 2 4 2 3
,m 100.00 92.00 94.00 91.00 88.00
NFA Iteration 1
4 Assumed node category — ADQF — ADQF ADQF ADQF
5 Qy, m3/min (2) - (4) (2) (3)
6 O,, m3/min 11 - 9 - 5 3
7 hx, m 4.222 — 5.826 — 3.928 2.290
8 Hr m (100.00) — 95.778 - 89.952 86.024 83.734
9 Obtained node category — SUPC — SUBC SUBC SUBC
10 Is node category
compatibility achieved? — Yes — No No No
(continued)
TABLE: EXAMPLE 13.6. (continued).
Row Node Pipe Node Pipe Node Pipe Node Pipe Node
No. Item 0 1 1 2 2 3 3 4 4
NFA Iteration 2
11 Assumed node category — — ADQF — CRIT — CRIT — CRIT
12 qr m3/min — — (2) — 2.305 0.850 — 3.472
13 Qx, m.3/min — 8.627 — 6.627 — 4.322 — 3.472
14 h~, m — 2.693 — 3.307 — 3 — 3 —
15 Hr m (100.00) — 97.307 — (94.000) — (91.000) — (88.000)
16 Obtained node category — — SUPC — PRTF — PRTF — SURF
17 Is node category
compatibility achieved? — — Yes — Yes — Yes — No
NFA Iteration 3
18 Assumed node category ADQF — CRIT — CRIT — ADQF
19 qr m3/min — — (2) — 2.305 — 1 322 — (3)
20 Ox, m3/min — 8.627 — 6.627 — 4.322 — 3 —
21 hx, m — 2.693 — 3.307 — 3 — 2.290 —
22 H,, m (100.00) — 97.307 — (94.000) — (91.000) — 88.710
23 Obtained node category — — SUPC — PRTF — PRTF — SUPC
24 Is node category
compatibility achieved? — — Yes — Yes — Yes — Yes
NFA Theory 375
achieved for nodes 2-4. Therefore, as the assumed node category is ade-
quate-flow (row 4) and the obtained node category is subcritical (row 9),
the node category is converted, by referring Table 13.2, to critical (row 11).
Step 4. The next NFA iteration is carried out (rows 11-17). Herein, the
node category compatibility is achieved for nodes 1-3, but it is not
achieved for node 4.
Steps 2-4 are repeated and iteration 3 is carried out. Herein, node cate-
gory compatibility is achieved at all nodes. Thus, the NFA solution is
reached, and it is as shown in rows 19-22.
It will be observed that nodes 1 and 4 can meet the demand completely,
while nodes 2 and 3 can meet the demand only partially.
Example 13.7
Carry out the NFA of the serial network shown in Figure 1: Example
13.7. The pipe head loss relationship is h = RQl8S in which h is in meters
2m)'min 3 4 3
0 M5 1
1 000
? j2 3 3t CK30
0£00
f t
0 3
100-00m 95-00m 96-00m 88-00m 89-00m
and Q in cubic meters per minute. The source HGL is 100.00 m, and data
for pipes and demand nodes are shown in Figure 1: Example 13.7 and also
given in rows 1-3 of Table: Example 13.7.
Solution
The NFA is carried out as shown in Table: Example 13.7. The solution
is obtained in 5 NFA iterations, and is shown in Figure 2: Example 13.7.
2 0 4 1-831
183
41 2
* t t3 1 1 4
10000m 97748m 95-138m 89-918 89- 00m
Row Node Pipe Node Pipe Node Pipe Node Pipe Node
No. Item 0 1 1 2 2 3 3 4 4
Data
1 0.05 0.10 — 0.20 — 0.30 —
2 t7;'?a, rrt3/min — — 2 — 3 — 4 — 3
3 H; ; ", m 100.00 — 95.00 — 96.00 — 88.00 — 89.00
NFA Iteration 1
4 Assumed node category — ADQF — ADQF — ADQF — ADQF
5 qr m3/min — — (2) _ (3) — (4) — (3)
6 Qx, m3/min — 12 — 10 — 7 — 3 —
7 hK, m 4.960 7.079 _ 7.319 — 2.290 _
8 Hr m (100.00) — 95.040 — 87.961 — 80.642 _ 78.352
9 Obtained node category — — SUPC — SUBC — SUBC — SUBC
10 Is node category
compatibility achieved? — — Yes — No — No No
TABLE: EXAMPLE 13.7. (continued).
Row Node Pipe Node Pipe Node Pipe Node Pipe Node
No. Item 0 1 1 2 2 3 3 4 4
NFA Iteration 2
11 Assumed node category — — ADQF CRIT CRIT — CRIT
12 qr m3/min — — (2) — -2.176 — 9.262 _ -1.917
13 Ox, m3/min — 7.169 5.169 — 7.345 — -1.917 _
14 hK, m — 1.912 — 2.088 — 8.000 — -1.000 —
15 H,, m (100.000) — 98.088 — (96.000) — (88.000) — (89.000)
16 Obtained node category — — SUPC — NEGF — SURF — NEGF
17 Is node category
compatibility achieved? — — Yes — No — No — No
NFA Iteration 3
18 Assumed node category — — ADQF NOFL — ADQF — NOFL
19 q;, m3/min — — (2) — (0) — (4) — (0)
20 Ox, m3/min — 6 — 4 4 — 0 —
21 hx, m — 1.376 — 1.300 — 2.599 — 0 —
22 Hr m (100.00) — 98.624 — 97.324 — 94.725 — 94.725
23 Obtained node category — — SUPC — SUPC — SUPC — SUPC
24 Is node-category
compatibility achieved? — — Yes — No — Yes — No
(continued)
3
oo
T A B L E : E X A M P L E 13.7. ( c o n t i n u e d ) .
Row Node Pipe Node Pipe Node Pipe Node Pipe Node
No. Item 0 1 1 2 2 3 3 4 4
NFA Iteration 4
NFA Iteration 5
4-16
Solution
The NHA for the design demands for the entire network (pipe 1-2 not
closed) is carried out, and the available HGL values at all nodes for the
design demands are given in column 4. All nodes are supercritical and the
design is satisfactory.
Pipe 1-2 is closed for repairs and the NFA solution of the network for
this condition is required. The NFA iterations are given in Table: Example
13.8. The values assumed for the iteration are shown in parentheses. All
demand nodes are assumed as adequate-flow for iteration 1 (column 5),
but all of them are found to be subcritical (column 8), and therefore, they
are converted to critical (column 10). At the end of iteration 2, nodes 2 and
5 are negative-flow, nodes 6 and 7 are partial-flow, while nodes 3 and 4 are
surplus-flow as shown by column 13. The negative-flow nodes 2 and 5 are
converted to no-flow, and the surplus-flow nodes 3 and 4 are converted to
adequate-flow. Critical nodes 6 and 7 (obtained as partial-flow) are re-
tained as critical (column 15) because node-category compatibility is
achieved for them. At the end of iteration 3, the obtained node categories
(column 18) are compatible with the assumed node categories (column 15)
for all demand nodes (column 19) and thus the NFA solution is reached
(columns 16 and 17). The complete NFA solution is also shown in Figure
2: Example 13.8. Herein, nodes 3 and 4 completely satisfy the demand,
TABLE: EXAMPLE 13.8.
(Values Used for NFA Iterations Are Shown in Parentheses.)
Pipe 1-2 Closed. NFA Iteration
(continued)
TABLE: EXAMPLE 13.8. (continued).
(92-27)
1-381
The following general observations can be made on the theory and pro-
cedure of the NFA.
1. Each NFA iteration forms one NHA solution under the assumed node
category combination. Therefore, the NFA solution can be obtained by any
one of the NHA methods described earlier. However, when some demand
nodes are assumed as critical nodes, they behave as source nodes (HGL
known). Therefore, it is preferable to use methods based on H or AH
equations for NHA solutions in the NFA iterative procedure. In methods
based on Q or AQ equations, the formation of pseudo loops is necessary.
As pseudo loops may change from iteration to iteration, and as some de-
mand nodes become source nodes, these methods are not well-suited for
carrying out NFA of distribution networks.
2. The NFA procedure starts (iteration 1) with all nodes as adequate-flow
nodes. The total network outflow is maximum but the solution may be in-
feasible because certain constraints from the set of alternate constraints
may not be properly satisfied. (The exception is when the NFA and NHA
are identical.)
3. At successive NFA iterations, the number of demand nodes at which
node category compatibility is achieved increases, and the solution,
though still infeasible, successively approaches feasibility.
4. When the node category compatibility is achieved for all demand
nodes, the procedure ends. The NHA solution obtained in the last iteration
under the assumed node category combinations, satisfies the alternate con-
straints properly at all demand nodes, and thus provides the NFA solution.
384 NODE FLOW ANALYSIS
The total network outflow is less compared to the earlier infeasible solu-
tions, but is maximum for a feasible solution and yields the optimal solu-
tion to the optimization problem of Equations (13.6a) . . . (13.6h).
5. In the usual optimization procedures, the constraints are satisfied and
the solution is feasible at each iteration. The objective function is op-
timized successively, and the procedure ends when no further optimization
is feasible. In the NFA procedure suggested herein, at least some con-
straints are initially violated (an exception occurs when the NHA and NFA
are identical) and the solution is infeasible. In the iterative procedure,
however, the number of violated constraints decreases, and the iterative
procedure ends when the constraints are properly satisfied at all nodes.
6. No proof is provided herein to show that the NFA solution is the opti-
mal solution. However, it can be observed that a slight variation in the
NFA solution results either in a feasible but suboptimal solution or in an
infeasible solution wherein one or more constraints are violated.
13.7 PROBLEMS
Problem 13.1: Solve Example 13.1 if Rol = 0.5 and Ro2 = 3.0.
Problem 13.2: In Example 13.1, nodes A and B are independent of each
other. Find (a) the limiting value of HA beyond which node A will start af-
fecting the design demand at B and (b) the limiting value of HB beyond
which node B will start aflecting the design demand at A.
Problem 13.3: Using the data of Example 13.1, determine the HGL
values of points fi, . . . fi6 described in Figures 13.2 and 13.3.
Problem 13.4: The minimum permissible HGL values at nodes 3 . . .
6 of the network of Example 8.3 are 94.20, 94.10, 90.00, and 92.50 m, re-
spectively. Find (a) the maximum possible outflow at node 5 without af-
fecting the outflows at all the other nodes and (b) the maximum possible
increase in the resistance constant of pipe 2 without affecting the outflow
at any demand node.
Problem 13.5: Obtain the NFA of the serial network of Example 13.5
when the nodal demands at nodes 1 . . . 4 change to 3, 2, 4, and 5 nrVmin,
respectively.
Problem 13.6: Obtain the NFA of the serial network of Example 13.6
when the minimum required nodal heads at nodes 1 . . . 4 are 93.00,
90.00, 92.00, and 91.00 m, respectively.
Problem 13.7: Carry out the NFA of the network of Example 13.8 when
(a) pipe 1-4 is closed for repairs and (b) pipes 1-2, 2-3, and 4-7 are simulta-
neously closed for repairs.
Calibration
14.1 INTRODUCTION
14.1.1 DEFINITION
The heads at source nodes can be determined fairly accurately. The wa-
ter elevations in the reservoirs and the heads supplied by the pumps can be
measured with appropriate gauges.
The rate of water supply at the source nodes can be either measured or
estimated with a reasonable degree of accuracy. The supply rate from a
reservoir can be estimated by observing the depletion in reservoir water
level and thereby in volume over a small time interval over which the sup-
ply rate can be assumed to be constant. Similarly, the supply rate from a
pumped source node can be measured with appropriate meters or esti-
mated from the head-discharge curves of the pumps, which should prefera-
bly be freshly determined prior to network observations.
The length and diameter of the pipes and the location of the various
valves can be obtained from the maps. However, proper selection of the
pipe head loss coefficients such as the Hazen-Wiliams coefficients or the
Darcy-Weisbach friction factors is rather difficult.
390 CAUBRATION
on nodal consumptions may be more reliable than that on pipe head loss
coefficients, or vice versa. For example, in a study on the Boston water
distribution network, Shamir and Howard [90] preferred to adjust pipe
head loss coefficients, because the data on nodal consumptions were more
reliable since they were highly correlated with the basic flow pattern pro-
vided by approximately twenty supply points. On the other hand, in a
study on the Minneapolis water distribution network, many pipe resis-
tances were fairly accurately known and therefore fine tuning of the model
was carried out by adjusting the resistance of a few important pipes.
It would always be possible to carry out localized adjustment in the
nodal consumptions and/or the pipe head loss coefficients and thus force
the model to fit the observations under a single demand pattern. However,
such a procedure would be neither logical nor systematic and would give
widely different calibrated models with different users, and such models
would fail to simulate the field network under other demand patterns.
Therefore, for the reliable calibration of distribution networks:
(1) As wide a range as possible of nodal consumptions should be consid-
ered. This may include minimum-flow, average-flow, peak-flow con-
ditions and fire-flow tests.
(2) Nodal flow consumptions and pipe head loss coefficients should both
be considered for adjustment.
(3) Calibration should not be localized but preferably be spread over the
entire network.
A good calibration technique therefore should preferably be able to con-
sider both nodal consumptions and pipe head loss coefficients for adjust-
ment, systematically decide the extent to which adjustment is necessary,
provide adjustment to all nodal consumptions and/or pipe head loss coeffi-
cients thereby spreading the calibration over the entire network, should
preferably be simple, and, if iterative, should yield calibration in a few
iterations.
The nodal demands and the pipe head loss coefficients (or the pipe resis-
tance constants) are initially assumed during calibration. When these
"values are assumed, the number of the unknown parameters—supplies at
source nodes and heads at demand nodes—is the same as the nodes in the
network and therefore a unique analysis of the network can be obtained.
This yields the values of the unknown parameters. However, when any one
or more of these unknown parameters such as the flow from a reservoir or
a pumping facility in a multisource network, the discharge or head sup-
292 CALIBRATION
plied by a booster pump in the network, or the head at a demand node are
measured, the measured values are compared with the corresponding
values obtained from the model analysis. If the discrepancy between the
measured and obtained values is more than a predetermined permissible
value, calibration of the model becomes necessary.
The calibration of the model is carried out by following one of the fol-
lowing two approaches:
(1) In the first approach of model calibration, one or more of the as-
sumed, i.e., predicted parameters, are adjusted, model analysis is
carried out, the obtained values are compared with the measured
values, and the procedure is repeated until satisfactory calibration is
achieved. In this approach the number of the unknown parameters is
the same as the number of the nodes in the network and the parame-
ters measured in the field are used explicitly, only to compare them
with the corresponding obtained values. Herein the measured values
are not an integral part of the model analysis.
(2) In the second approach of model calibration, the measured parame-
ters are treated as known parameters and are directly used in the
model analysis. This approach provides an additional equation for
each measured value and thus permits the model user to introduce an
additional unknown parameter for each measured value. The addi-
tional unknown parameter can be the head loss coefficient of a par-
ticular pipe in the model or a multiplying factor to correct some or
all predicted pipe head loss coefficients. In this approach the mea-
sured values are used implicitly in the model calibration.
Several methods, based either on the implicit approach or the explicit
approach are available for calibration of network models. These methods
can be classified into two categories: (1) methods in which only the pipe
head loss coefficients are adjusted and (2) methods in which both parame-
ters, pipe head loss coefficients and nodal demands, are adjusted. These
methods are now described.
Herein, only the pipe head loss coefficients are adjusted and it is pre-
sumed that the assumed nodal demands are fairly accurate and do not need
any adjustment.
Rahal et al. [182] proposed a calibration technique for adjusting pipe re-
sistance constants when discrepancy exists between the measured and ob-
Calibration Methods 393
Example 14.1
4OL/S
80
2-7 are demand nodes with estimated water use as shown in the figure.
The length, diameter, and the actual and predicted HW coefficients for the
pipes are given in Table: Example 14.1. The heads are measured at nodes
4 and 7 and the values are 54.10 m and 48.88 m respectively. Calibrate the
model assuming that the predicted demands at all nodes and the HW coef-
ficient of pipe 7 are reasonably accurate.
Solution
Because the heads are measured at two nodes, two additional conditions
are available and therefore the network is divided into two zones. Let these
zones be zone 1 consisting of pipes 1-6 with adjustment factor A; and zone
2 consisting of pipes 8 and 9 with adjustment factor B. (Pipe 7 is excluded
from calibration.) The model has three basic loops, loops I . . . III. Two
pseudo loops IV and V are introduced as shown in the figure. Let the dis-
charges in pipes 1 . . . 9 be Q1 . . . Q9 respectively, in the directions
shown in the figure. Let the estimated resistance constants of pipes 1 . . .
9 for the predicted HW values be Rol . . . Ro9, respectively. The Q equa-
tions can now be written as
Q4 + Q6 - 20 = 0 (3)
Qs - Q* + Qi - 80 = 0 (4)
Q* ~ Qi ~ Q9 - 70 = 0 (5)
Qg - 40 = 0 (6)
Herein both the pipe head loss coefficients and nodal demands are ad-
justed.
14.3.2.1 Trial-and-Error Method
One of the methods used for simultaneously adjusting the predicted pipe
head loss coefficients and nodal demands is to use the conventional trial-
396 CALIBRATION
and-error approach. The model user relies on his experience and judge-
ment for adjustment. Several runs of the model might be necessary to ar-
rive at a reasonably acceptable calibration.
When the model is calibrated under two loading conditions—high and
low flows—qualitative guidelines proposed by Walski [178] can be used.
These include the following:
(1) Adjust the predicted pipe coefficients if there is a similar error at
high and low flows.
(2) Adjust the predicted pipe coefficients and nodal demands by similar
amounts if the model is accurate at low flow but inaccurate at high
flow.
(3) Adjust the predicted nodal flows if the model is accurate at high flow
but inaccurate at low flow.
14.3.2.2 Method Proposed by Walski
A °L fl4n
A U4
~ b(S + F) - aS ' '
B = (14-2)
b(S + F) - aS ~ ^
Calibration Methods 397
in which
ft - H, \"
a = (14.3)
ft - H,F
and
b =
ft/ - ft, \°- (14.4)
Using the adjustment factors A and 5, the adjusted nodal demands qja
and pipe coefficients CIU are obtained from
= BCr, (14.6)
Walski's method is explicit and can be used for multiple loadings. It can
be used iteratively to improve calibration results. The method is illustrated
in Example 14.2.
Example 14.2
Test node
5
node with HGL of 60.00 m and nodes 2 . . . 7 are demand nodes. Nodes
4 and 7 are used as test nodes with fire flows of 160 L/s and 75 L/s, respec-
tively. The observed and predicted HGL values at test nodes are given in
Table 1: Example 14.2. The length, diameter, and actual and predicted
398 CALIBRATION
HW coefficients for the pipes are the same as in Table: Example 14.1. The
actual and predicted nodal demands are given in Table 2: Example 14.2.
Calibrate the model using Walski's method.
Solution
As the fire flow tests are carried out at two test nodes, the network is
divided into two zones: zone 1 comprised of nodes 2, 3,4 and 5, and pipes
1 . . . 6, governed and therefore adjusted by the test results at node 4 and
zone 2 comprised of nodes 6 and 7, and pipes 8 and 9, adjusted by the test
results at node 7. The HW coefficient of pipe 7 is not considered for adjust-
ment as in Example 14.1.
As the test results for low and high flows at node 4 are in agreement, no
adjustment is proposed in zone 1. However, the test results at node 7 give
higher residual heads for predicted conditions and therefore nodal flows
and pipe resistance constants require adjustment in zone 2. Considering
test node 7,
60 - 48.88 0.54
= 1.091
60 - 50.54
60 — 7 15 \ 0 5 4
(
S7 = qo6 + qol = 70 + 40 = 110 L/s
F7 = 75 L/s
1.091 X 75
A7 = = 0.865
1.160(110 + 75) - 1.091 X 110
B .
1.091 ~ °793
Calibration Methods 399
Using the value of B7 the CHW values for pipes 8 and 9 are adjusted as
91.2. The demands at nodes 6 and 7 are also adjusted and these values are
shown in Table 2: Example 14.2.
The adjusted values are used and the pressure heads at nodes 4 and 7 are
obtained as shown in Table 1: Example 14.2. The predicted HGL values at
node 4 have slightly deteriorated, but those at node 7 have considerably
improved and, in general, the calibration is satisfactory.
Bhave [186] proposed a method that also simultaneously adjusts the pipe
coefficients and the nodal demands, and is thus similar to Walski's method.
However, in Walski's method, the nodal demands and therefore the inflow
at the source node increases or decreases depending upon the value of the
adjustment factor A. In Bhave's method it is presumed that the inflow at the
source node can be measured or estimated fairly accurately, and therefore
it remains fixed during calibration. When the nodal demands are initially
assumed, they are adjusted, if necessary, so that their sum equals the esti-
mated inflow. (This adjustment, incidentally, would take care of the water
losses in the network.) In the calibration, however, redistribution of nodal
demands may take place so that some predicted nodal demands increase
and others decrease, the total increase being equal to the total decrease for
the entire network.
Consider a skeletonized network shown in Figure 14.1, in which s is a
source node and ty and t2 are demand nodes where heads are measured si-
multaneously under a particular demand pattern. Divide the network into
two zones, one for each test node as shown in the figure. Let AFi and AF 2
be the demand adjustments for zones 1 and 2, respectively. Naturally,
AF 2 = —AFi so that the demand adjustment for the entire network is
zero.
2 30 40 40.0
3 120 100 100.0
4a 30 20 20.0
5 100 80 80.0
6 50 70 60.6
7b 20 40 34.6
Total 350 350 335.2
a
Rreflowof 160 L/s.
b
Fire flow of 75 L/s.
400 CALIBRATION
Zone 2
Figure 14.1 Skeletonized network.
Consider a path from 5 to f, as shown by thick lines in the figure. Let the
pipes on this path be labelled 1, 3, . . . , x, . . . . For this path,
(14.7)
in which Rip, R3p, . . . Rxp, . . . - predicted resistance constants of pipes
1, 3, . . . , x, . . . ; Qlp, Q3p, . . . . Qxp, . . . = discharges in pipes 1, 3,
. . . , JC, . . . for the predicted nodal demands and pipe resistance con-
stants; Hs = head at source node 5; and H,,p = predicted head at test
node ti.
Assume that the entire correction AF, for zone 1 is applied to pipe 1 on
path s-tu and also that this correction reduces along the path s-tt in pro-
portion to the discharge carried by a pipe so that the adjustment for pipe
x, i.e., AQX = (Qxp/QiP)AFj. If Bis the global adjustment factor for pipe
resistance constants, then for path s-t,,
BRlp(Qlp BR3P(Q3P
n{Hs Htv)
B(HS - Hup) - ~ AF, = ft - ft, (14.11)
Sl
=
C'° D0.54 C*P (14.13)
The zonal demand adjustments AF, and AF2( = - AF,) are distributed
in proportion to the predicted nodal demands in zones 1 and 2, respec-
tively. Therefore, the adjusted nodal demand q]a is given by
AFZ
1+ (14.14)
in which qjp = predicted demand at node;'; AFZ = total nodal flow adjust-
ment for zone z\ and Nz = set of demand nodes in zone z.
This method is explicit and can be used iteratively so that B — 1 and
AF —• 0, and is stopped when reasonable agreement is reached between
the observed and adjusted results.
As the nodal demands and pipe coefficients are simultaneously adjusted,
at least two observations are necessary. These can be obtained by measur-
ing heads at two demand nodes under a single loading condition or at a de-
mand node under normal and fire flow conditions or at a demand node
under two different loading conditions. When more than two observations
are taken, the network is divided into several zones and the number of Bs
and AFs is taken such that their sum exceeds the number of available ob-
servations by one. (This is because the last AF value equals the negative
sum of the other AF values.) The paths are selected such that all the paths
402 CALIBRATION
(1) together cover all zones, (2) lie in the general direction of flow, and (3)
contain pipes that carry the majority of flow towards the test nodes as seen
from the model run under predicted condition.
Example 14.3
Zone,3
Zone.i
Solution
The model is run for the predicted condition and the analysis results are
given in Tables 1, 2 and 3: Example 14.3. As the heads are measured at two
test nodes under normal and fire flow conditions, four additional equations
are available. The network is divided into three zones as shown in Figure:
Example 14.3. These zones are zone 1 comprising pipes 1, 2 and 3; zone
2 comprising pipes 4, 5 and 6; and zone 3 comprising pipes 8 and 9. Let
the pipe adjustment factors for these zones be Bu B2, and B3, respectively.
Pipe 7 is excluded from adjustment as in Example 14.2. Nodes 2, 3, 4, and
5 are affected by flow at test node 4, and therefore they are grouped to-
gether. Similarly, nodes 6 and 7 are grouped together. Let the total adjust-
ment for nodes 2 . . . 5 be AF so that adjustment for nodes 6 and 7 is
— AF. Thus, we have now four unknowns—B t , B2, B3, and AF—for the
four known conditions.
We shall now consider path 1-3-2-4, comprising pipes 2 and 3 (zone 1)
and 4 (zone 2) and path 1-6-7 comprising pipes 8 and 9 (zone 3). Note that
the paths are selected such that all the three zones are covered, they lie in
the general direction of flow and contain pipes that carry the majority of
flow towards the test nodes as seen from model run results for the pre-
dicted condition as given in Table 1: Example 14.3.
TABLE 1: EXAMPLE 14.3.
Predicted
After
Node Actual Predicted Iteration I Final
2 30 40 43.4 40.8
3 120 100 108.4 101.6
4a 30 20 21.7 20.4
5 100 80 86.8 81.2
6 50 70 57.1 67.5
7b 20 40 32.6 38.5
Total 350 350 350 350
a
Fire flow o1 160 L/s.
D
Fire flow o' 75 L/s.
Practical Considerations 405
Considering path 1-3-2-4 and using the pipe flows and head losses from
Table 1: Example 14.3 for the predicted low and high flow conditions, we
get, respectively,
(4.89 + 0.78)3. +
Similarly, for path 1-6-7 comprising pipes 8 and 9, we get for the low and
high flows, respectively,
(.7.57 + 22.57)6, -
Her. The method of Ormsbee and Wood calibrates the model implicitly in
only one iteration. Other methods consider the test results explicitly. If
necessary, several iterations can be carried out to finetune the model so
that the predicted results satisfactorily match the observed ones. However,
it should be remembered that the calibration methods provide adjustments
for groups of pipes and clusters of nodes, and therefore the adjusted and
observed results would not match for water use at individual nodes and for
the resistance constant of individual pipes. Therefore, in practice, a few
iterations would be sufficient to obtain calibration commensurate with the
accuracy of the collected data.
14.4.2 SOURCES OF ERROR IN FIELD DATA COLLECTION
(1) If the adjusted pipe coefficient is much lower than the initially pre-
dicted value (Cxa < C p ), it is implied that the water carrying capac-
ity of the pipes is less than anticipated. While this may be due to an
increase in pipe roughness, it may also be due to a closed or a partly
closed valve or a pipe that is installed and shown on the maps but not
actually commissioned.
(2) If the corrected pipe coefficient is much greater than the initial pre-
diction (Cxa > Cxp), it is implied that the water carrying capacity of
the pipes is more than anticipated. While this may be due to a lesser
increase in pipe roughness than anticipated, it may be also due to the
presence of a rich mesh of smaller diameter pipelines excluded dur-
ing skeletonization, or the presence of new large pipelines placed in
service in the area but not shown on the maps and therefore excluded
in the model.
Practical Considerations 407
(3) If the adjusted water use is much less than the initial prediction
(qja -4 qjp), the model user should attempt to identify water uses not
in operation when the data were collected. For example, the field
data might have been obtained on a school holiday, or on a day when
the workers near the test node may be on strike. It may also be due
to some water users getting water from neighbouring pressure zones.
(4) If the adjusted water use is much greater than the initial prediction
(qJa > qjp) the model user should look for unusual water uses such
as the filling up of a municipal swimming pool, excessive lawn
watering, or new industries excluded during the survey. The error
may also be due to an illegal connection or a large main break. A
pressure reducing valve in operation but excluded in the model may
also show higher head loss indicating more water use than antici-
pated.
(5) If the pressure drops considerably during a fire flow test, nearby
water uses would get reduced during testing, especially in low pres-
sure zones. In such cases it would not be proper to simply add higher
discharge to the normal water use, but the predicted normal water
use should be lowered, in the fire flow test.
(6) When a model consistently predicts much higher heads than those
observed in the field, it is worthwhile to run fire hydrant flow tests
several times with several residual pressure gauges and reconstruct
the hydraulic grade line in the test area. This may reveal that the
water reaches the hydrant through a path not thought to be significant
earlier.
(7) A remote part of a distribution system may also be supplied by an-
other system. If the two systems are modeled as separate, the pre-
dicted heads might be too low, and if they are connected, the heads
may be too high. The results may agree better if the systems are con-
nected at some places in the model, or the water uses are redistrib-
uted between the two systems. In large zones, field testing may con-
tinue for a whole day. Because of the diurnal variation in water use,
it might be necessary to adopt different water use multipliers during
the testing to match the predicted and observed results.
(8) In dynamic model calibration, the reservoir water levels may not
match well. This may be because the pumps are turned on and off at
the beginning of a time step in the model while in reality they may
turn on and off during the time period. It might be necessary to
reduce the time step.
It is thus seen that there are several sources of error in field data collec-
tion and adjustment and therefore proper model calibration is not easy.
408 CALIBRATION
The methods described earlier help the model user in calibration but con-
siderable judgement and intuition is required on the part of the model user
for satisfactory calibration. Calibration involves description of the whole
from observations on parts and therefore as much data as possible should
be collected. Fire flow tests with several residual pressure gauges would be
useful in resolving the valve status and topology of the distribution system.
However, some discrepancies between model results and field observa-
tions are better resolved by posing specific questions to valve crew super-
visors, construction inspectors, and operations personnel. Many problems
arise in calibration due to the lag between installation of new pipes and
showing them on system maps. Accurate and updated maps of distribution
systems are a prerequisite for good calibration.
Walski [188] has given an exhaustive list of the equipment and described
their utility in collection of data in the field. Such equipment includes the
following:
(1) Regularly calibrated pressure gauges, preferably of the same type
and size, ranging from 100 kPa to 2,000 kPa (To obtain good results,
the selected range should barely cover the range of the pressures ex-
pected in the system.)
(2) Pitot gauges for measuring hydrant discharges having a usual range
of 2,500-5,000 L/min, and vernier calipers to measure the inside di-
ameter of the hydrant outlets
(3) Hydrant-cap and pitot fittings
(4) Hydrant wrenches and pliers
(5) Safety equipment such as safety shoes, reflecting vests, hard hats,
traffic cones, and flags
(6) Recording equipment including clipboard, papers, data sheets, ordi-
nary and video cameras
(7) Street maps and distribution system maps at several scales
(8) Communication equipment
(9) Survey equipment for elevation data collection
(10) Equipment for differential pressure measurement for determination
of the head loss coefficients of pipes
1
Miscellaneous Topics
Nodal Heads
Sr. Balancing Balancing Nodal Pipe (Finite Element .
No. Item Heads Rows Heads Loops Discharges Method)
1 Unknown parameters AQ AH H AQ 0 H
and equations used
2 Number of unknowns Minimum Medium Medium Minimum Maximum Medium
4 Mode of analysis Analysis in parts, considers Analysis as a whole, all loops or nodes are considered simulta-
one loop or node at a time neously
5 Suitability for hand Suitable. Networks with Not suitable except for small networks
calculation several nodes can be ana-
lyzed if sufficient time and
manpower are available
TABLE 15.1. (continued).
Nodal Heads
Sr. Balancing Balancing Nodal Pipe (Finite Element
No. Item Heads Rows Heads Loops Discharges Method)
6 Number of iterations Usually large Small and usually not more than 10
7 Effect of network size Increases as the size of the Generally independent of the size of the network
on the number of network increases
iterations
8 Time per iteration Small Large
9 Effect of network size Increases gradually as the Increases rapidly as the size of the network increases
on time per iteration size of the network in-
creases
10 Problem of conver-
May be present Reported to be absent
gence
11 Computer storage Small and increases grad- Large and increases rapidly with the size of the network
capacity requirement ually with the size of the
network
412 MISCELLANEOUS TOPICS
Fj = 0,j = 1 . . .J (15.1)
= 0,j = 1 . . .J (15.3)
(15.4)
dv dv* dyk
changes in J selected free variables vf. The results from sensitivity analy-
sis are approximate because they are derivatives at a point rather than
finite changes. However, sensitivity analysis is a convenient and rapid way
of comparing the effect of all possible changes.
Example 15.1
The network shown in Figure: Example 15.1 has node 1 as source node
with HGL = 100.00 m and nodes 2 . . . 6 as demand nodes with demands
WO-OOrn
19
and minimum required HGL values as given in Table: Example 15.1. The
head loss in pipes is given by h = RQl S5 in which /; = head loss in
meters, R = pipe resistance constant (shown underlined in the figure),
and Q = pipe discharge in cubic meters per minute. The NHA solution of
the network is given in the table. Using sensitivity analysis, determine the
maximum possible outflow at node 3 without affecting the outflows at all
other nodes.
Solution
When the outflow at node 3 increases, the residual heads at all nodes de-
crease. Therefore, for the maximum possible outflow at node 3 without af-
fecting the outflows at the other nodes, the available HGL values at all
nodes, including node 3, should not be less than the minimum permissible
values. Thus, we are interested in knowing the impact of an increase in q3
on the nodal HGL values. Therefore, qu H2 . . . H6 are selected as the six
free variables. The node-flow continuity equations at nodes 1 . . . 6 are,
respectively,
- // 3 \ 0
<?3 = 0
2.5 4.0 6.0
(3)
(4)
1.00 = 0 (5)
2.0 4.0
3X -°-non9 0 = (6)
dF1
dqt 9// 2 u dH4 dq3
dF2 3F 2 9// 2
0 0
dH2 dH3 dq3
9^3
dF, 9F, dF
dH2 dH4 dH6 dq3
_ _ (7)
9^3
dF4
0U o
U ^ ^ ^ 0
9// 3 dH4 dHs dq3
9tf4
1FS dFs dFs dFs
0 0 0
oH4 9<?3 dq3
9// 5
dFb dF6 BHb
0 0
dH3 dHs dH6 dq3 dq3
416 MISCELLANEOUS TOPICS
1 0.1454 0 0.2163 0
0 -0.3672 0.2218 0
dH.,
-1
x (8)
8HA
0
dq3
dHs
0
dq3
dH,
0
dq3
The solution of Equation (8) yields the values given in column 8 of Ta-
ble: Example 15.1. These values are the changes in the values of the free
Sensitivity Analysis 417
Example 15.2
Find the available heads at all demand nodes if the resistance constant
of pipe 2 in Figure: Example 15.1 is 1.5 instead of 1. Use the NHA solution
given in Table: Example 15.1 and adopt the sensitivity analysis.
Solution
Value of v, from
3
<7i. m /min 0 0.5 0 5.00 5.00
H2, m -2.516 0.5 -1.258 91.06 91.15
H3, m -4.165 0.5 -2.082 87.88 88.01
HA. m -5.618 0.5 -2.809 89.88 90.11
H5, m -4.999 0.5 -2.500 86.85 87.04
H6, m -4.627 0.5 -2.313 86.55 86.70
given in the table. Some difference is observed between the two results be-
cause the result from the sensitivity analysis is approximate, as stated
earlier.
As stated in Chapter 5, certain pumps may have more than one flow for
the same head. Therefore, in practical networks having several pumps, the
possibility of numerous solutions cannot be ruled out.
For control elements such as check valves and pressure reducing valves,
the head loss characteristics are different depending upon their mode of
operation. The common approach for considering such elements in the
analysis is to assume a particular operating mode for each element, solve
the network, and check whether the assumed operating modes of different
elements are correct. If they are correct, the network is presumed to be
solved; if not, the operating modes are changed, the network is resolved,
and the process is terminated when the assumptions for the operating
modes are correct for all elements. Even though this procedure is straight-
forward, a problem arises in deciding how accurate the solution must be
before one uses the solution to test the validity of the assumption, regard-
ing the operating mode of the control elements. Naturally, if the solution
being tested is not sufficiently accurate, a wrong conclusion regarding the
validity of the operating modes may be drawn. It is thus seen that the crite-
ria used for terminating the iterative procedure could be critical.
the quantitative criterion for the overall performance of the New York City
water supply tunnel system. They expressed the overall performance as
fV\V
= 0 (15.8)
c dt dt ID
(15.9)
422 MISCELLANEOUS TOPICS
Hp {QP
~ HR ~ QR) QR
I QR - tR) = 0
(15.10)
Hs = Qs) Qs
~ ^ ~ ~
(15.12)
These four equations containing four unknowns xp, tp, Hp, and Qp can be
solved by the method of specified time intervals [154]. The procedure is
continued for successive time intervals. Reservoirs, pumps, and valves can
also be included in unsteady flow analysis.
The emphasis throughout the text has been on understanding the funda-
mentals of the different network analysis methods. Therefore, the networks
considered for the illustrative examples were small so that they could be
tackled by hand calculation. However, as the analytical methods require
t
^ A x * -Ax*
y
t+At
R
Figure 15.1 C* and C~ characteristics on jtf-plane.
Problems 423
15.7 PROBLEMS
Problem 15.1: In Example 15.1, using sensitivity analysis, find the max-
imum possible outflow at node 6 without affecting the outflows at all other
nodes.
Nomenclature
Le =
length of equivalent pipe
Lx —
length of pipe x
In =
natural logarithm
log =
logarithm to base 10
M =
mass flow rate, number of source nodes in net-
work
m = dimensionless form factor for pipe roughness;
node; pipe
N = Manning roughness coefficient; number of de-
mand nodes in network
N. = Manning roughness coefficient of equivalent
pipe
Nx = Manning roughness coefficient of pipe x
Nz = set of demand nodes in zone z
NFA = node flow analysis
NHA = node head analysis
NOFL = no-flow node
n = exponent of discharge in general pipe head loss
relationship
P = wetted perimeter; power delivered by pump;
solution point in jtf-plane of characteristic grid
in unsteady flow analysis
PRTF = partial-flow node
PRV = pressure reducing valve
p = pressure
PJ = average pressure above the minimum standard
at node j
Q = volumetric flow rate; pipe discharge
Qd = discharge at downstream end of pipe
Qe = discharge in equivalent pipe
Qij = discharge in pipe ij
Qp — discharge supplied by pump
Qu = discharge at upstream end of pipe
Qx = discharge in pipe at distance x; discharge in
pipe x
Qxp = predicted discharge in pipe x
QxU) = discharge in pipe x obtained after f* iteration
,QX = assumed discharge in pipe x for /'* iteration
Qo = constant discharge in head-discharge relation-
ship of pumps
Qox = known value of discharge in pipe x
QA, <IB = flow at upstream node A and downstream node
B, respectively
JS <7«es = design demands at nodes A and B, respectively
™\ ^B"" = maximum possible flows at nodes A and B, re-
spectively
qb(t) = flow at boundary element at time /
Nomenclature 431
81 Hazen, R. "Pumps and Pumping Stations," J. New England Viaer Works Assoc.,
67:121 (1953).
82 Fair, G. ML, J. C. Geyer and D. A. Okun. Elements of Water Supply and Waste
Water Disposal. New York, NY:John Wiley and Sons, Inc. (1971).
83 Jeppson, R. W. and A. L. Davis. "Pressure Reducing Valves in Pipe Network
Analysis," J. Hydraul. Div., Am. Soc. of Civil Engrs., 102(HY7):987-1001
(1976).
84 Jacoby, S. L. S. "Design of Optimal Hydraulic Networks,"/ Hydraul. Div., Am.
Soc. of Civil Engrs., 94(HY3):641-661 (1968).
85 Bhave, P. R. "Noncomputer Optimization of Single-Source Networks" J. Envi-
ron. Engrg. Div., Am. Soc. of Civil Engrs., 104(EE4):799-814 (1978).
86 Bhave, P. R. "Optimization of Gravity-Fed Water Distribution Systems: Theory,"
J. Environ. Engrg., Am. Soc. of Civil Engrs., 109(1): 189-205 (1983).
87 Harary, F. Graph Theory. Reading, MA:Addison-Wesley Publishing Co. (1969).
88 Shamir, U. and C. D. D. Howard. "Water Distribution Systems Analysis," J. Hy-
draul. Div., Am. Soc. of Civil Engrs., 94(HYl):219-234 (1968).
89 Shamir, U. and C. D. D. Howard. Closure to "Water Distribution Systems Analy-
sis," J. Hydraul. Div., Am. Soc. of Civil Engrs., 96(HY2):577-578 (1970).
90 Shamir, U. and C. D. D. Howard. "Engineering Analysis of Water Distribution
Systems," J. Am. Water Works Assoc, 69(9):510-514 (1977).
91 Gofman, E. and M. Rodeh. "Loop Equations with Unknown Pipe Characteris-
tics," J. Hydraul. Div., Am. Soc. of Civil Engrs., 107(9): 1047-1060 (1981).
92 Bhave, P. R. "Rules for Solvability of Pipe Networks" / Indian Water Works
Assoc, 22(1):7-10 (1990).
93 Bhave, P. R. "Analysis of Water Distribution Networks: Part I," J. Indian Water
Works Assoc, 13(2): 149-154 (1981).
94 Bhave, P. R. "Analysis of Water Distribution Networks: Part II," J. Indian Water
Works Assoc, 13(3):245-252 (1981).
95 Bhave, P. R. "Analysis of Water Distribution Networks: Part HI,"/ Indian Water
Works Assoc, 13(4):305-314 (1981).
96 Bhave, P. R. "Unknown Pipe Characteristics in Hardy Cross Method of Network
Analysis," J. Indian Water Works Assoc, 18(2):133-135 (1986).
97 Bhave, P. R. "Unknown Pipe Characteristics in Network Analysis: Formulation
of Equations" (to be published).
98 Bhave, P. R. "Unknown Pipe Characteristics in Network Analysis: Solution of
Equations" (to be published).
99 Saaty, T. L. and J. Bram. Nonlinear Mathematics. New York, NY:McGraw Hill
Book Co. (1964).
100 Cross, Hardy. "Analysis of Flow in Networks of Conduits or Conductors," Bulle-
tin No. 286, Engrg. Experiment Station, Univ. of Illinois, Urbana, IL (1936).
101 Cornish, R. J. "The Analysis of Flow in Networks of Pipes," / Inst. of Civil En-
grs., 13:147-154(1939).
102 Liu, K. T. H. "The Numerical Analysis of Water Supply Networks," International
Assoc. of Hydraul. Research, 13th Congress, Kyoto, Japan, 1:35-42 (1969).
103 Bhave, P. R. "Rapid Convergence in Hardy Cross Method of Network Analysis,"
/ Indian Water Works Assoc, 16(l):l-5 (1985).
104 Borgohain, T. P., A. K. Khare, D. J. Desai, K. P. Singh and S. Banerjee. "Com-
446 References
205 Coulbeck, B. "A Computer Program for Calibration of Water Distribution Sys-
tems," J. Inst. ofmter Engrs. and Scientists, 38(l):79-80 (1984).
206 Wood, D. J. Computer Analysis of Flow in Pipe Networks Including Extended Pe-
riod Simulations (KYPIPE). Lexington, KY:Office of Continuing Education, Uni-
versity of Kentucky (1980).
207 Charles Howard and Associates, Ltd. Water Distribution Network Analysis: SPP8
User's Manual. Victoria, B.C., Canada:Charles Howard and Associates, Ltd.
(1984).
Index
453
454 Index
461