Pressure Transients in Water Engineering
Pressure Transients in Water Engineering
Pressure Transients in Water Engineering
water engineering
A guide to analysis and interpretation
of behaviour
John Ellis
Independent Hydraulics Consultant and Honorary Research Fellow,
University of Glasgow, UK
Published by Thomas Telford Publishing, Thomas Telford Ltd, 1 Heron Quay, London E14 4JD.
www.thomastelford.com
A catalogue record for this book is available from the British Library
ISBN: 978-0-7277-3592-8
Whilst every reasonable effort has been undertaken by the author and the publisher to
acknowledge copyright on material reproduced, if there has been an oversight the publishers
will endeavour to correct this upon a reprint.
The right of John Ellis to be identified as the author of this work has been asserted by him in
accordance with the Copyright, Designs and Patents Act 1988.
All rights, including translation, reserved. Except as permitted by the Copyright, Designs and
Patents Act 1988, no part of this publication may be reproduced, stored in a retrieval system or
transmitted in any form or by any means, electronic, mechanical, photocopying or otherwise,
without the prior written permission of the Publishing Director, Thomas Telford Publishing,
Thomas Telford Ltd, 1 Heron Quay, London E14 4JD.
This book is published on the understanding that the author is solely responsible for the statements
made and opinions expressed in it and that its publication does not necessarily imply that such
statements and/or opinions are or reflect the views or opinions of the publishers. While every
effort has been made to ensure that the statements made and the opinions expressed in this
publication provide a safe and accurate guide, no liability or responsibility can be accepted in
this respect by the author or publishers.
Acknowledgements xv
Notation xvi
Introduction 1
iii
1.10.1 Thermosetting plastics, 14
1.10.2 Thermoplastics, 15
1.11 Failure modes of pipes, 17
1.12 Maximum pressure and allowable amplitude of surge
in plastic pipes, 18
1.13 Minimum pressures, 18
1.14 When is analysis necessary?, 19
Chapter 3 Interpretation of a 27
3.1 Fluid properties, 27
3.2 Influence of the conduit wall, 28
3.3 Simple expression for a, 29
3.4 Variation of a with conduit shape, 32
3.5 Influence of gas on a, 32
3.6 The effect of sewage, 38
iv
Chapter 6 Boundaries 60
6.1 Types of boundary, 61
6.2 Reservoirs and tanks, 62
6.3 Branches and changes in pipe properties, 64
6.3.1 Specific cases — number of pipes ¼ 1, 67
6.3.2 Specific cases — change of cross-sectional area, 67
6.4 Response of a large pipe or trunk main, 69
6.5 Actuated valves and pipeline fittings, 71
6.5.1 Terminal valves, 73
6.5.2 In-line valve, 74
6.5.3 Automatic control valves, 75
6.6 Use of more than one time step, 77
6.7 Non-reflecting boundary, 78
6.8 Other bifurcation conditions, 82
6.8.1 Bifurcation with operating valves, 82
6.8.2 Isolating valves, 83
6.9 Continuous drawoff, 84
v
9.5.2 Actuated valve, 119
9.6 Control of transient pressures and estimation of valve
operating time, 122
vi
Chapter 11 Flywheels 159
11.1 Moment of inertia, 159
11.2 Flywheels, 160
11.3 Limitations on flywheel size, 161
11.4 Pipeline limitations, 162
11.5 Case study with different pump speed options, 163
11.6 Flywheels on a larger system, 167
11.7 Booster pump installations, 170
11.8 Multi-pump installations, 170
11.9 Advantages of flywheels, 171
Appendix Moment of inertia, 171
vii
Chapter 13 Further aspects of pressure vessels 215
13.1 Pressure vessel types and their fittings, 215
13.2 Vessels having an air—water interface, 215
13.2.1 Air compressors, 215
13.2.2 Control of gas charge/liquid level, 217
13.2.3 Other vessel fittings, 218
13.3 Bladder vessels, 219
13.4 Positioning a pressure vessel, 221
13.5 Installation with air valves, 225
viii
16.4.1 Above-ground storage tanks, 284
16.4.2 Vacuum disconnecting valves, 286
16.5 Air valve operation, 292
16.6 Summary of influence of discharge arrangements, 293
ix
19.2 Initial valve sizing, 362
19.3 Valve positioning, 363
19.4 Analysis of behaviour, 363
19.5 Automatic surge control valve, 365
19.6 Surge anticipation valve, 366
19.7 Pumping station pressure relief valve, 366
19.8 Grove regulator, 369
19.9 High head relief valves, 371
19.10 Bursting disk, 375
x
20.11.6 Backflow check valve, 407
xi
Chapter 23 Amplification of transient pressures 474
23.1 Transmission of pressure waves through a branch
connection, 474
23.2 Pressure wave transmission through a change of cross-
section, 476
23.3 Meeting of opposing pressure waves, 478
23.4 Pressure waves in a suction main, 479
23.4.1 Protection of the rising main, 479
23.4.2 Conditions in the gravity main, 479
23.5 Amplification within a network, 481
23.5.1 Kirkleatham Lane Pumping Station, 481
23.5.2 System modelling, 484
23.5.3 Recorded pipe bursts and pressure extremes, 485
23.6 Wellfield transients, 488
23.6.1 Collector pipeline system, 488
23.6.2 Borehole and wellhouse configuration, 490
23.6.3 Wellfield operating conditions, 490
23.6.4 Pumpset inertia, 492
23.6.5 Sequenced pump operation, 493
23.6.6 Pumping failure, 494
23.6.7 Air valve operation, 496
23.7 Potential for amplification, 497
xii
24.4.2 Comments on observed behaviour, 515
24.4.3 Modelling behaviour, 516
24.4.4 Remedial measures, 518
24.4.5 Final comments, 519
References 520
Index 527
xiii
Acknowledgements
The author would like to express his appreciation for assistance received
from a number of sources. To Chris Bone of Morton & Bone Services,
UK agents for Charlatte bladder pressure vessels, who provided much
useful material on this form of surge protection. To Sam Gilbert of
Glenfield Valves Limited, Kilmarnock, Scotland, for many interesting
conversations regarding valve performance; also for providing infor-
mation on a wide variety of valve types. This book would not have
been possible without the work of engineers over many years who
laboured to improve understanding of this subject. From among now
retired engineers I would like to thank George A. Milne, formerly of
Crouch & Hogg, Glasgow, for taking time to discuss the subject at
length.
xv
Notation
xvi
Fro ¼ FVo2 =z
f friction factor in equation hf ¼ fLV 2 =ð2gDÞ
fc ovality correction factor P
fe equivalent overall friction factor ¼ KL D=L þ f
fos factor of safety
fs pipe support factor
G shear modulus of pipe material (N/m2 )
g acceleration due to gravity (m/s2 )
Hb inertial head rise in a branch pipe
Hbs base outlet head at a PRV
Hi inertial head rise ¼ aV=g (m)
Hj inertial head rise at a junction of pipes
Ho initial head or piezometric level
Hp head developed by pump
Hr liquid level in reservoir or large tank
Hv constant outlet head at a PRV or head rise at a valve
Hx inertial head rise at a point ‘x’ along a pipeline
h pressure head ¼ p=ðgÞ at a point in the pipeline (m)
also absolute pressure head in vessel
habs absolute pressure head
hatm atmospheric pressure head
hf head loss due to pipeline flow resistance (m)
hvap vapour pressure head (gauge)
I moment of inertia of a pumpset (kg/m2 )
Is ¼ s3 =12
i pipe number
i.l. invert level of pipeline
Jþ invariant or quasi-invarient quantity propagating along
Cþ path (m/s)
J invariant or quasi-invarient quantity propagating along
C path (m/s)
Jo initial value of invariant or quasi-invarient
K bulk modulus of fluid (N/m2 )
Kg bulk modulus of gas (N/m2 )
Ks bulk modulus of matter in sewage (N/m2 )
ki þ1 if isolating valve open and 0 if valve closed
ko equivalent uniform sand grain roughness height of new
pipe (mm)
ks equivalent uniform sand grain roughness height (mm)
L pipe length measured along axis (m)
lnðxÞ natural logarithm (to base e) of x
xvii
log10 ðxÞ common logarithm (to base 10) of x
m number of waves in buckling calculation
mg mass of gas dissolved/unit volume of liquid
N rotation speed of a pump (rev/min)
n polytropic exponent of gas law
Pr power output of a pump
p gauge pressure ¼ gh at a point in the pipeline (N/m2 )
pabs absolute pressure (N/m2 )
patm atmospheric pressure (N/m2 )
pb critical buckling pressure
pex external pressure
pH percentage hydrogen ion concentration in water (%)
pk critical pressure
ps set pressure at which regulator starts to open
Q flow rate at a pipeline section at some time (m3 /s)
Qair air inflow rate at air valve
Qe flow rate from an ejector system
Qm mixture discharge rate from a jet pump
Qo initial steady flow
Qout outflow to a gravity sewer
Qp flow rate of pumped medium through a jet pump
Qr discharge through relief valve or regulator
q ratio Qp =Qm for a jet pump
also vertical load on pipe/metre axial length
qi discharge at a demand point
R reactive force between sliding surfaces
Rc collapse rate of flexible pipe ¼ minutes/% decrease in
diameter
Re Reynold’s number of flow at a section ¼ VD=
rm distance from pivot to centre of gravity of check valve
door
rp radius of pivot pin
S gradient of head loss due to pipeline resistance
So initial steady flow gradient of head loss due to resistance
Sr ring stiffness/unit length of pipe
s pipe wall thickness (m)
sf sliming factor
sj þ1 for Cþ characteristic and 1 for C characteristic
sv Vo =jVo j provided jVo j > 0
T torque
Tabs absolute temperature
xviii
Teff effective valve closure time
Tr resultant thrust on valve hinge
t time (s)
V mean velocity of flow at a cross-section at some time
(m/s)
Vo initial velocity
Vol volume of liquid þ volume of gas ¼ Voll þ Volg (m3 )
Volf volume of free or undissolved gas (m3 )
Volg volume of gas dissolved in liquid (m3 )
Volin inlet volume of sleeve valve
VolL volume of liquid (m3 )
Volp volume of a pipeline
Vols volume of matter in sewage (m3 )
Volt total volume (m3 )
Ws buoyant weight of check valve door
x axial distance along pipeline from some datum (m3 )
Z head relationship for jet pump
z pipe centreline level measured from a common
horizontal datum (m3 )
also d/spreservoir level pressure vessel level þ habs
z ¼ Vo = ½LA=ðgAs Þ
zw weir crest level
¼ y=Dm
k rate of increase in roughness height with time
x distance increment in numerical scheme
xp small distance increment
t time increment in numerical scheme
tp small time increment
y vertical change in pipe diameter
! rotation speed change over a time increment
relief valve diameter
void fraction ¼ ð1 þ f Þðpatm =pabs Þð1=Þ
¼ Dmin =Dm ¼ 1
f void fraction of free gas ¼ Volf =ðVolf þ Voll Þ
s sewage matter void fraction ¼ Vols =Volt
1 void fraction of free gas
2 void fraction of dissolved gas ¼ Volg =ðVolg þ Voll Þ
ratio of specific heats for gas ¼ Cp =Cv
h horizontal radial deformation of pipe
i initial pipe deflection
xix
y vertical radial deformation of pipe
" maximum strain
"cr critical strain for a perfectly round pipe
"f effective maximum strain
pump efficiency
angle between flats of a polygonal duct (rad)
also angle of pipeline with the horizontal
also angle of check valve door from the vertical
Poisson’s ration for pipe wall material
also coefficient of sliding friction
kinematic viscosity of fluid (m2 /s)
ratio hm =h
liquid density (kg/m3 )
d density of disk material
g density of gas (kg/m3 )
m mean density of mixture ¼ g þ ð1 Þ (kg/m3 )
also mean density of mixture exiting from a jet pump
p density of pumped medium through a jet pump
s density of sewage matter (kg/m3 )
xx
Introduction
1
Pressure transients in water engineering
2
Introduction
desirable that he should have on his team someone familiar with aspects
of transient behaviour. The better informed the client’s representative
the less scope for misunderstanding.
However, until now this type of analysis remains an activity which is
not encountered on a very regular basis within consultancy organisa-
tions generally and a considerable time may elapse between one
analysis and the next. Coupled with turnover of staff, this can make
the process of gaining experience somewhat protracted. The new
graduate or novice engineer may be provided with a quite sophisticated
modelling tool but still be unsure of what design cases should be consid-
ered and for how long a computer simulation is required to be continued
to ensure that all transient activity of significance has been predicted.
Being new to surge analysis, the novice may be uncertain whether or
not the predictions which he has obtained from the modelling
exercise are realistic, or possibly contain the consequences of some
mistake in the information provided to the computer or even the
result of using the program inappropriately. These concerns are not
unique to pressure transient analysis and have also arisen in connection
with the use of computer software in other areas of engineering such as
structural analysis where similar concerns have been raised.
However, the frequency with which transient investigations are
undertaken may increase generally if a novel approach to mains con-
dition assessment using controlled surges becomes widespread. As
described by Misiunas et al. (2007), this technique is in the early
stages of development but offers hope of a low-cost means of mains
assessment.
In a recent article, Gordon (2004) described how a computer
program written by a mathematician without adequate understanding
of engineering had resulted in expensive repairs and several months
delay in commissioning of a scheme. This article highlighted the
dangers of subcontracting the hydraulic transient study based on
lowest cost, without determining whether the program was adequate
for the task and had a proven track record on other projects. The
client had not worked on a project of this kind for a few years and
experienced personnel were no longer available. In another instance,
Gordon (2006a) described how a high head turbine runner was fabri-
cated with an even number of blades with the corresponding number
of wicket gates also being even. The basic principle laid down back in
1941, that an even number of wicket gates should be accompanied
by an odd number of runner blades, had been lost. The result was
that a large number of blades were cracked after a couple of months
3
Pressure transients in water engineering
4
Introduction
5
Pressure transients in water engineering
6
1
Motivation for hydraulic
transient analysis
In times past when methods of analysis and materials were perhaps less
clearly understood than now or quality control was less exacting, a
relatively large factor of safety would be applied to any design. With
increasingly sophisticated computer modelling capabilities available
and with new materials available and a better understanding of these,
there has been a movement towards more detailed analysis and better
use of materials. Strengths of pipes and pipeline fittings are of para-
mount importance in determining whether a design is safe and so this
chapter examines some aspects of pipe materials and allowable internal
pressures before proceeding to consider analysis of pressure transients in
pipeline systems. Also included is some mention of pipe linings as these
can have a bearing upon allowable pressures. Other aspects such as
trench conditions and flexibility which can influence allowable pres-
sures will be considered at a later stage.
7
Pressure transients in water engineering
8
Motivation for hydraulic transient analysis
systems comprising a mix of old and new pipes and the properties of
both types are important as well as any deterioration in these over
time. Since pipeline systems will usually be in service for many years
and be subject to extension and upgrading, to suit changes in
demand for instance, an understanding of properties of older or even
historic materials must still be retained as some reanalysis is often
necessary to accommodate changes to a network.
Table 1.1. Ratings and test pressures for various classes of cast iron
9
Pressure transients in water engineering
Description Class
15 20 25
10
Motivation for hydraulic transient analysis
Test pressure 50 40 32 25
(bar g)
Maximum working 40 25 25 25
pressure (bar g)
11
Pressure transients in water engineering
DN 100 1.6
DN 150 2.1
DN 200 2.4
DN 250 2.7
DN 300 3.0
DN 350 3.1
DN 400 3.2
DN 500 3.4
DN 600 3.6
DN 700 3.8
DN 800 and above 4.0
12
Motivation for hydraulic transient analysis
1.9.1 Bitumen
Bitumen linings comprise a mixture of 80% bitumen and 20% dry lime.
The bitumen deteriorates with age. Burstall (1997) reported ‘puffing-
up’ of bitumen to give a rough tuberculated appearance. Deterioration
can also result in loss of lining in raw water pipelines.
13
Pressure transients in water engineering
treated water mains a long life can be anticipated and with correct
chlorine residuals gives a clean surface with little or no change in
carrying capacity of the main. It is a brittle material which can fail
catastrophically. Maximum pipe deflection (change in diameter/original
diameter) should be 2% to avoid risk of cracking. Failure of lining has
also been attributed to pressure transients.
14
Motivation for hydraulic transient analysis
1.10.2 Thermoplastics
Thermoplastics are synthetic materials softened by application of heat
and are capable of repeated softening by subsequent heating. The
principal types of thermoplastics used in the water industry are polyvinyl
chloride (PVC) and polyethylene (PE). Other materials such as acrylo-
nitride butadience styrene (ABS) and polypropylene (PP) may be used
in special circumstances.
PVC was developed by German scientists shortly before World War
II. Meanwhile in the UK, ICI scientists had discovered a means of
producing polyethylene. Being lightweight, flexible and virtually free
from chemical attack made these materials contenders for pipeline
materials.
Initial poor performance of uPVC was attributed to a combination of
poor installation, interference damage and unsuitable operating condi-
tions including high transient pressures. From about 1973, when the
industry fully appreciated the limitations of the material, the failure
rate of uPVC has been below that of spun grey iron (cast iron).
Pressure pipes for water services may be obtainable in four classes as
indicated in Table 1.6.
Pressure class
B C D E
Pressure, bar g 6 9 12 15
15
Pressure transients in water engineering
Table 1.7. The five series for HDPE defined in DIN 8074
Series
1 2 3 4 5
16
Motivation for hydraulic transient analysis
Class
A B C D E
Pressure, bar 3 6 9 12 15
17
Pressure transients in water engineering
18
Motivation for hydraulic transient analysis
19
Pressure transients in water engineering
20
2
Derivation of basic equations
Two basic methods of analysis have been used extensively in the inves-
tigation of unsteady flows in pressure conduit networks. The simplest of
these has the assumption that the flowing liquid is incompressible and
that the conduit walls are rigid or inelastic. This method is sometimes
known as the rigid-column analogy and has been widely employed in
the analysis of longer-period oscillations such as occur in low-pressure
tunnels of hydropower networks.
The second approach considers the liquid to be compressible and the
conduit walls to be deformable resulting in a variation of flow along the
pipeline or tunnel at any instant. This ‘elastic’ approach coupled with
the large computing power and memory of modern desktop computers
has the potential to reproduce both high-frequency and longer-period
aspects of hydraulic transients in a way that was hitherto not possible.
21
Pressure transients in water engineering
Hi
Diameter = D M
Velocity = V
Cross-section A
22
Derivation of basic equations
23
Pressure transients in water engineering
Piezometric line
p/(rg) + ∂[p/(rg)]/∂x dx
h = p/(rg)
rAV
H dx
z z + ∂z/∂x dx
Horizontal datum
24
Derivation of basic equations
becomes:
@V=@t þ V @V=@x þ 1=ðgÞð@H=@x @z=@xÞ þ gS þ g @z=@x ¼ 0
or finally:
@V=@t þ V @V=@x þ g @H=@x þ gS ¼ 0 ð2:5Þ
In the mass balance equation:
@V=@x þ 1=ð@=@p @p=@t þ V @=@p @p=@xÞ
þ 1=Að@A=@p @p=@t þ V @A=@p @p=@xÞ ¼ 0 ð2:2Þ
Writing,
25
Pressure transients in water engineering
or,
p
a ¼ 1= fð1= d=dp þ 1=A dA=dpÞg ð2:6Þ
then the mass balance equation finally becomes:
a @V=@x þ g=af@H=@t þ V @H=@xg g=aV dz=dx ¼ 0 ð2:7Þ
26
3
Interpretation of a
The value of a is dependent upon the properties of the fluid and d=dp
and upon the cross-sectional area of the conduit and its deformation
characteristics A and dA=dp.
For water, the value of K varies only slowly with temperature with a
decrease of the order of 2% between 08C and 508C. A value of
2.137 GN/m2 at 218C is given by Knapp et al. (1970) with Thorley
and Enever (1979) providing a slightly higher value of 2.193 GN/m2
at 208C. Seawater has a value of 2.4 GN/m2 at 258C. For a computer
simulation conducted assuming an almost uniform temperature, K
can be considered constant and independent of flow conditions. If it
27
Pressure transients in water engineering
Diameter = D
ds
Change of internal
pressure = dp
28
Interpretation of a
also
Youngs modulus E ¼ d=ðdc=cÞ ¼ dpD=ð2s dD=DÞ
¼ D2 =ð2sÞ dp=dD
assuming that the pipe is entirely free to contract in length. Re-
arranging,
dD=dp ¼ D2 =ð2sEÞ
Cross-sectional area A ¼ =4D2 and A þ dA ¼ =4ðD þ dDÞ2 , or
neglecting second-order effects:
dA ¼ =2D dD or dA=dD ¼ =2D
and
1=A dA=dD ¼ 2=D
then,
1=A dA=dp ¼ 1=A dA=dD dD=dp ¼ 2=D D2 =ð2sEÞ
¼ D=ðsEÞ ð3:3Þ
29
Pressure transients in water engineering
Steel
a = 1350 1300 1200 1100 1000 900 800 700 600
Hard cast iron
500
300
Concrete
Asbestos cement
200
Polyester
(glass fabric)
Polyester 10
(glass fibre)
100
Perspex
PVC (hard)
PVC 50
(blow resisting)
PE (hard) 9
PE (soft)
log10 (E )
8 1 2 3
log10 (Dc 1/s)
a = ÷(K/r)/÷[1 + DKc 1/(sE)]
30
Interpretation of a
a circular pipe section. A wave speed value for wood stave pipe was
given by Fedosoff and Szpak (1979) as 762.0 m/s. The influence of
constraints has been considered by various authors such as Streeter
and Wylie (1967), Thorley and Enever (1979) and Jaeger (1977).
Expressions for the constraint factor c1 for different circumstances,
including thick-walled conduits defined by D=s < 25, and tunnels
have been included in Table 3.1.
Expansion joints
c1 ¼ 1 =2
Expansion joints
throughout
c1 ¼ 1 c1 ¼ D=ðD þ sÞ þ 2s=Dð1 þ Þ
c1 ¼ 2ð1 þ Þ
31
Pressure transients in water engineering
32
Interpretation of a
33
Pressure transients in water engineering
34
Interpretation of a
35
Pressure transients in water engineering
36
Interpretation of a
1200
1000
800
Wave speed (m/s)
600
Vf = 1%
400 Vf = 0.1%
Vf = 0.01%
Vf = 0.001%
Vf = 0.0001%
200
0
0.1
0.5
0.9
1.3
1.7
2.1
2.5
2.9
3.3
3.7
4.1
4.5
4.9
5.3
5.7
6.1
6.5
6.9
7.3
7.7
8.1
8.5
8.9
9.3
9.7
10.1
Pressure (bar(abs))
Fig. 3.3. Acoustic wave speed plotted against pressure for various void fractions
solution at this gas release head and that equilibrium concentration was
achieved with no delay. Other modellers have adopted an exponential
decay of rate of gas release.
Figure 3.3 shows how wave speed changes with pressure for a range of
free gas concentrations. The gas fractions are those at atmospheric
pressure and range from 0.0001% to 1.0% by volume.
Adoption of a wave speed which varies with changing pressure
conditions within each pipeline of a network has important implications
with regard to the way in which computations are conducted. In the
presence of free gas the liquid in the pipeline system over most of the
range in free gas concentration has effectively become a fluid which
retains a density dominated by the weight of the liquid component
but with compressibility more like that of the gas. Only at the extremi-
ties of concentration range does this alter, with liquid properties
dominant when only very small amounts of gas are present, and gas
properties dominant if only quite small amounts of liquid are available.
Evolution of gas from solution and subsequent expansion of the free
gas has the effect of introducing an additional source of fluid to the
system. Suppose that upstream of a region in which free gas is present
and expanding, the flow is decelerating at some rate dV=dt. Down-
stream of the zone of gas expansion the flow will be decelerating at a
lesser rate as the flow from upstream will be augmented by liquid
released from the length of pipe containing the expanding gas in
37
Pressure transients in water engineering
p=ðm gÞ ¼ KL V 2 =ð2gÞ
38
Interpretation of a
39
Pressure transients in water engineering
40
4
Characteristic equations
41
Pressure transients in water engineering
along dx=dt ¼ V þ a
and
ð ð ð ð
dV g 1=a dH þ g S dt þ g V=½aðV aÞ dz ¼ constant
along dx=dt ¼ V a
or if a is assumed constant:
ð ð
V þ g=aH ¼ constant g S dt þ g V=½aðV þ aÞ dz ð4:5aÞ
42
Characteristic equations
and
ð ð
V g=aH ¼ constant g S dt g V=½aðV aÞ dz ð4:6bÞ
43
Pressure transients in water engineering
Ignoring any small changes in kinetic energy head then in the absence
of pipeline resistance and losses generally H remains essentially
constant along a path dx=dt ¼ V a so that along some length of
conduit over which velocity changes by an amount dV and over
which pipeline elevation changes by an amount dz:
ð ð
dV ¼ g V=½aðV aÞ dz or since V a; dV g V=a2 dz
44
Characteristic equations
and substituting:
1=K þ Dc1 =ðsEÞ ¼ 1=V dV=dp
Since pressure change dp is due entirely to elevation change dz then
dp ¼ g dz, or,
½1=K þ Dc1 =ðsEÞ ¼ 1=V dV=ðg dzÞ
Substituting a2 and rearranging,
gV=a2 dz ¼ dV
Between limits 0 and 1 then:
ð
V1 V0 ¼ g V=a2 dz ð4:7Þ
45
Pressure transients in water engineering
46
Characteristic equations
4.4.1 Corrosion
Metallic pipes such as steel and cast iron are liable to corrode. Hydraulic
Research Station paper No. 4 gives some data from studies by Cole-
brooke and White (Colebrooke, 1939) whereby the increase in k may
be estimated for older asphalted cast-iron pipes. The rate of increase
in roughness height with time represented by can be obtained from:
2 log10 ð12Þ ¼ 3:8 pH ð4:14Þ
where ¼ ft=year for cast iron asphalted pipes. The current value of
k is then found from the equation,
k ¼ ko þ t ð4:15Þ
t being pipe age in years and ko is the new pipe roughness height.
4.4.2 Sliming
Wastewater mains are liable to sliming. Recent research has shown that
k is primarily a function of velocity:
k ¼ V ð2:34Þ ð4:16Þ
where
¼ 0:054 (lower bound)
¼ 0:446 (mean) and
¼ 3:660 (upper bound)
Ninety-five per cent of observations lie between these upper and lower
limits.
From equation (4.16) it will be noted that as V ! 0 then k ! 0. To
avoid unreasonably small values of k it is recommended that a minimum
value of k be used equal to the value ko for the clean pipe.
For all types of wastewater rising mains, the following has also been
proposed as in the Biwater Manual (1988).
47
Pressure transients in water engineering
f=ð2=DÞðV 2 Þm ðt1 t0 Þ
where ðV 2 Þm may be represented as 12 ðV12 p
þ V02 Þ; Vm may be represented
1
as 2 ðV1 þ V0 Þ the arithmetic mean, or as ðV1 V0 Þ the geometric mean,
Smith (1969). p 2
The representation
p of Vm as ðV1 V0 Þ has some advantage as Vm
becomes ½ ðV1 V0 Þ2 or simply V1 V0 . Preserving the sign of resistance
gradient then, this can be written V1 jV0 j. In the equations for the
quasi-invariant, ignoring the second small integral for the effect of
changing pipe elevation:
V g=aH ¼ J f=ð2DÞVjV0 jðt t0 Þ
or rearranging,
f1 ðÞf=ð2DÞjV0 jðt t0 ÞgV g=aH ¼ J ð4:17Þ
This formulation also preserves the linear nature of the equations.
48
5
Application of characteristic
equations
49
Pressure transients in water engineering
dx/dt = V ± a
Wavelets
Velocity of flow = V
Horizontal datum
50
Application of characteristic equations
C– characteristics
dx/dt = V – a
C+ characteristics
dx/dt = V + a
+t
Limits of region
influenced by P
Typical intersections
+x
51
Pressure transients in water engineering
Representative characteristics
C– C+
+t
+x
52
Application of characteristic equations
C– paths
C+ paths
+t
+x
significant difficulties and has been used effectively by Larsen (1976) for
example, in his investigation of undersea outfalls. Considerable
complexity is introduced where a looped or branching system is
encountered or where internal changes of pipe material or diameter
are present for instance. Although not insurmountable, the computa-
tional difficulties of following progress of characteristics through more
complicated systems are sufficient to make this approach unattractive
for general use.
The irregular distribution of intersection points at which solutions
are obtained also means that some post-processing is required to
make the predictions easy to understand. Typically the end-user
wishes to see time variations of pressure, head, flow, etc., at specified
locations and this requires interpolation within the irregular mesh of
solutions.
Almost all computational schemes designed for general use employ a
regular or ‘fixed’ mesh arrangement which provides predictions at
predefined locations throughout the pipeline system and at regular
increments of time. Characteristic paths are selected, from the 1
53
Pressure transients in water engineering
54
Application of characteristic equations
Dt
+t C– characteristics C+ characteristics
dx/dt = V – a dx/dt = V + a
Dt
Dx Dx
+x
55
Pressure transients in water engineering
Dt
+t
C+ characteristics
C– characteristics
Dt
Dx Dx
+x
propagation of the pressure waves but contains the same risks of disper-
sion of the wave front.
By setting:
x=t ¼ dx=dt ¼ constant
56
Application of characteristic equations
Vu/s Vd/s
Dx/2 Dx/2
Dx Dx
+x
Vu/s Vd/s
C+ characteristic C– characteristic
Dt
Vu/s(o) Vd/s(o)
the need for interpolation is avoided and with it the effects of any
numerical dispersion.
In this scheme, free air or gas released from solution, which can influ-
ence the wavespeed a, is handled by ‘collecting’ all of the gas in the
neighbourhood of a computing point, into a single ‘bubble’ occupying
the cross-section of the pipe at a local computing point (Fig. 5.7). A
discontinuity of velocity (and flow) occurs across the bubble which
expands and contracts as pressure changes. Behaviour of the gas mass
is represented using a gas law.
If the gas void fraction in the neighbourhood of a solution point and
at the start of a time increment t is o at absolute pressure pabsðoÞ and
57
Pressure transients in water engineering
58
Application of characteristic equations
each pipe, based upon the composite pipe and liquid system. Gas release
is modelled as a separate phase which may be either distributed along
the top of the pipe in the form of a continuous gaseous phase, or
could be represented as a discrete bubble where it occurs at a well-
defined summit on the pipeline. A flexible model will include the
option of representing free gas or vapour in one or both of these forms.
59
6
Boundaries
In earlier chapters it was described how the task of determining flow and
pressure changes along a pipeline could be achieved by using the
network of characteristics and the values propagating along these
paths. The challenge remains of finding corresponding values of flow
parameters at the many and various features and fittings which are to
be found in most pipeline networks. Each feature will act to transmit
and/or reflect transient pressure waves in a manner which is particular
to that type of fitting or element of the pipeline. Each of these elements
is referred to as a boundary. Analysis of response of these features
requires that an additional ‘auxiliary’ equation or equations has to be
introduced to complete the process of solution and this chapter intro-
duces just some of the wide variety of possible pipeline elements and
the corresponding equations.
At each computing point within a pipeline, values of velocity V and
total head H were obtained through solving a pair of simultaneous equa-
tions at the intersection point of characteristics travelling in opposing
directions within the pipeline. At the upstream and downstream extremi-
ties of a pipeline only a single characteristic intersects the boundary at any
instant of time, as shown in Fig. 6.1. Since values of both V and H are
required it is necessary to introduce additional equation(s), sometimes
referred to as auxiliary equations, to provide a complete description of
conditions at the extremity of the pipeline. The nature of boundaries is
varied and it is worthwhile spending a little time in considering some
of the types and common assumptions which are made regarding these.
As a rule the steady flow characteristics of a boundary such as a valve
or pump are assumed applicable to unsteady conditions.
Boundaries usually represent some actual physical feature of a
network but may also be numerical devices introduced for convenience
60
Boundaries
Characteristics Dt
Single characteristic Dt
at boundary
Dx Dx
61
Pressure transients in water engineering
V g=aH ¼ J or V ¼ J ðÞg=aHr
62
Boundaries
Reservoir level = Hr
Pipeline
Dx
Boundary node
Internal node
Dt
Characteristics
Single characteristic
arrives at boundary
When flow is into the pipeline from the reservoir, kinetic energy is
developed and losses are incurred so that Hr ¼ H þ ð1 þ KL ÞV2 =ð2gÞ,
where KL is the entry loss coefficient. The relationships at the pipe
inlet at any time t become:
63
Pressure transients in water engineering
64
Boundaries
Pipe No. 1
Pipe No. 2
Pipe No. 3
Positive direction of flow may be
defined as towards the junction or
away from the junction in any pipe
Area = A3
(a)
C2
Dt
C1
1 C3
Dx1
Dx2
Dx3
(b)
65
Pressure transients in water engineering
Common head = H
Area = A2
Area = A1
V1 V2
Wave speed = a1
Wave speed = a2
+x
C+ characteristic C– characteristic
Dt
Dx1 Dx2
66
Boundaries
67
Pressure transients in water engineering
Internal node
Boundary node
Dt
Single characteristic
at boundary
Dx
68
Boundaries
since
A1 =a1 A2 =a2 and jA1 C1 j jA2 g=a2 Hr j
H ¼ 1=gðA2 g=a2 Hr Þ=ðA2 =a2 Þ
or
H ¼ Hr
69
Pressure transients in water engineering
Aa
+x
Ab
+x
+x
J+ = VoAb/Aa
J– = 0
J– = –Vo
present:
Jþ ¼ Vo Ab =Aa þ g=a0 ¼ Vo Ab =Aa ð6:12Þ
And from the downstream section of aqueduct where zero flow is taking
place:
J ¼ 0 þ g=a0 ¼ 0 ð6:13Þ
Due to valve closure in the branch main, the head change Hj at the
junction of the branch main with the aqueduct is given by equation
(6.4). Substituting in this equation then:
Hj ¼ 1=g½ðVo ÞAb þ Vo Ab =Aa Aa 0Aa =ðAb =ab þ Aa =aa þ Aa =aa Þ
or
Hj ¼ 2=gAb Vo =ðAb =ab þ 2Aa =aa Þ ð6:14Þ
Head rise at the junction as a proportion of the initial head rise at the
valve in the branch is then:
Hj =Hb ¼ 2=½1 þ 2ðAa =aa Þ=ðAb =ab Þ
70
Boundaries
or for
a ¼ constant Hj =Hb ¼ 2=ð1 þ 2Aa =Ab Þ
Inserting values for area and wave speed into the above equation it
can be shown that as aqueduct area/branch area increases the response
at the connection tends more and more towards that of a simple reser-
voir. This relationship allows the investigator to assess the implication
of replacing the aqueduct connection by a simpler boundary such as a
reservoir.
71
Pressure transients in water engineering
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
1000 n
Valve type* D (in.) ao/A Port shape
A Gate 12 1.0 Circular
B Butterfly 24 0.65 –
C Globe 6 1.0 Shaped (iii)
D1 'In-line' regulator 12 1.0 Rectangular
D2 'In-line' regulator 9 Shaped (i)
100 i ii n iii
n
n
Port shapes
* See Fig. 4
Valve headloss coefficient, Kn
D2 † Model
E
Valve full open
C
B F
10 10
D1
G
A
Ko
1 1
Fig. 6.7. (a) Head-loss coefficients for different valve types; (b) control valve type
72
Boundaries
Discharge level = Hr
Piezometric level = H
Pipeline Valve
Internal node
Boundary
Dt Single characteristic
at boundary
Dx
73
Pressure transients in water engineering
Piezometric line
Hu/s
Valve Hd/s
C+ characteristic C– characteristic
Dt
Dx Dx
74
Boundaries
75
Pressure transients in water engineering
76
Boundaries
For a viable result Hu=s Hd=s Ko V2 =ð2gÞ where Ko is the fully open
valve head loss coefficient.
Pipeline
Dtp
Dtp
Dt
Dtp
Dtp
Dtp
77
Pressure transients in water engineering
78
Boundaries
Non-reflecting boundary
Dx
Non-reflecting boundary
Jp = Vp ± g/aHp
Characteristics without
the area being modelled
Dt
Jo = Vo ± g/aHo
Jo = Vo ± g/aHo
Jo = Vo ± g/aHo
Characteristics within the
area being modelled
79
Pressure transients in water engineering
given by:
Jo ¼ Vo g=aHo ð6:24Þ
If no reflection or transient effects arrive at the boundary from the main
network then the value Jo will remain the same regardless of the time
which has elapsed. From within the modelling area the pressure wave
effect from the transient event will be given by the prevailing values
of velocity Vp and Hp as shown in Fig. 6.12. So that:
Jp ¼ Vp ðÞg=aHp ð6:25Þ
is the quasi-invariant value at the end of the segment of characteristic
leading to the point t at which a solution is to be found. The effects of
friction and any other influences along the segment of path have been
added to yield this final value of Jp . The solution after any time step is
then given by:
V ¼ ð Jp þ Jo Þ=2 and H ¼ ð Jp Jo Þ=ð2g=aÞ ð6:26Þ
This technique is clearly approximate since although no effects will
return to the boundary from features along the pipeline such as
changes of cross-section, bifurcations, valves, etc. there will be a
more gradual and continuous change due to the action of pipeline
resistance. If necessary, some allowance can be made for this effect by
slowly modifying the value of Jo over time, producing a quasi non-
reflecting boundary. A modified technique which includes the effect
of resistance in the system outwith the modelled area can be found in
Chapter 20.
Non-reflecting boundaries can be used in a variety of contexts. Two
such boundaries can be employed on either side of a booster pumping
station for instance, as depicted in Fig. 6.13. The non-reflecting
boundary can be used in many circumstances. For instance, consider
the case of branch pipe network connected to a larger aqueduct
system. If only events within the branch network are of interest then
non-reflecting boundaries could be installed a short distance along
the aqueduct on each side of the connection point. Alternatively, the
junction itself could be made non-reflecting using the same principle,
by introducing the initial steady flow values of velocity and head in
the aqueduct upstream and downstream of the junction, so that head
at the junction Hj becomes:
Hj ¼ ðJb Ab þ ðVo Ab =Aa þ g=aa Ho ÞAa
ð0Ab =Aa g=aa Ho ÞAa Þ=ð2Aa =aa þ Ab =ab Þ
80
Boundaries
Hd/s
Hu/s
Booster
pumping
station
Dx Dx
or
Hj ¼ ðVo Ab þ 2g=aa Ho Aa Jb Ab Þ=ð2Aa =aa þ Ab =ab Þ
¼ ðconstant Jb Ab Þ=ð2Aa =aa þ Ab =ab Þ ð6:27Þ
where constant ¼ Vo Ab þ 2g=aa Ho Aa .
There are limitations to the application of the above approach.
Clearly the period of analysis has to be limited to < 2L=a where this
is the reflection time from the non-reflecting boundary to the first
source of significant reflection in the pipeline system outwith the area
of interest. The use of this type of boundary is also approximate in
the sense that reflections are continually being transmitted back from
81
Pressure transients in water engineering
Piezometric levels
H1
Valve 1
Area = A1
V2
V1 Valve 2 Area = A2
V3
Area = A3
Common horizontal datum H2 H3
82
Boundaries
are:
V1 g=a1 H1 ¼ J1 for pipe No. 1
V2 g=a2 H2 ¼ J2 for pipe No. 2
V3 g=a3 H3 ¼ J3 for pipe No. 3
H1 H2 ¼ KL1 V22 =ð2gÞ for valve No. 1
H1 H3 ¼ KL2 V32 =ð2gÞ for valve No. 2
A1 V1 A2 V2 A3 V3 ¼ 0 at the bifurcation
where KL1 and KL2 are head loss coefficients for valves No. 1 and No. 2
respectively. This set of relationships contains two quadratic equations
and cannot be solved explicitly. An iterative solution is necessary for
each time step. Ellis and Tint (1976) described a method of solving
these equations.
An alternative to solving the complete set of equations together
using a specially constructed routine is to compromise and separate
the valves from the branch so that small increments xi are introduced
between the bifurcation and each valve with a corresponding small time
step t. Beyond each valve and on pipe No. 1, boundaries of the type
described in section 6.6 can be introduced allowing larger distance and
time increments to be used in the rest of the network. In this way
solutions can be found for the six unknowns using existing routines
and without introducing significant inaccuracies.
83
Pressure transients in water engineering
H1
Piezometric line @ open valve
H4 H2 H3
A1 Shut valve
A2
A4
A3
Open valve
with the summation over all pipes. The head at the valve on each pipe is
then:
Hi ¼ H if ki ¼ 1 or Hi ¼ ai =g Ji if ki ¼ 0
Velocity at the valve on each pipeline i is then:
Vi ¼ Ji ðÞg=ai Hi
84
Boundaries
Piezometric line
H Demand points
zi
z Pipeline
85
Pressure transients in water engineering
Area 50R
Pressure vessel
Area 49R
Distribution area
Blank flange
86
7
Valve closure in a simplified
system
Vo þ g=aHo ¼ Jþ
87
Pressure transients in water engineering
Valve
Vo
+X
Length = L
t = 4L/a (seconds)
Time +
J+ = Vo + g/aHo
t = 3L/a (seconds)
J– = Vo – g/aHo
t = 2L/a (seconds)
H = Ho
J+ = –Vo – g/aHo
J– = –Vo – g/aHo
C– characteristics
t = 0 (seconds)
C+ characteristics
J+ = Vo + g/aHo
88
Valve closure in a simplified system
Static ‘layers’
Steady flow velocity
Vo
Additional travel of
first and second layers
89
Pressure transients in water engineering
90
Valve closure in a simplified system
Compression wave
Ho a Hi
M
V = Vo
V=0
(a) 0 £ t £ L/a
Rarefaction wave
Ho a Hi
M
V = –Vo V=0
Rarefaction wave
Ho
M
a V=0 –Hi
V = –Vo
Ho Compression wave
M
a V=0 –Hi
V = +Vo
Fig. 7.3. Pressure wave front at different stages following valve closure
91
Pressure transients in water engineering
value arriving at the valve face at this time remains constant at:
Jþ ¼ Vo þ g=aHo
At the valve face V ¼ 0 so that:
Jþ ¼ Vo þ g=aHo ¼ 0 þ g=aH
or
H ¼ Ho aVo =g ¼ Ho Hi
From initial valve closure at time ¼ 0 until time ¼ 2L=a, when a
reflected pressure wave returns to the valve from the reservoir, piezo-
metric level at the valve remains at Ho þ Hi . After time ¼ 2L=a,
piezometric level at the valve has fallen to Ho Hi . This time 2L=a
from the initiation of the pressure transient is the minimum time
before any ‘relief of pressure’ occurs at the source of the hydraulic
transient. It is sometimes referred to as the ‘critical period’. Remaining
with this simple example, it follows that any closure time which is
2L=a will produce the same maximum head rise as an instantaneous
valve closure. Closure times >2L=a will start to bring benefits of pres-
sure relief with consequent maximum inertial head rise Hi being
<aVo =g. Figure 7.3 shows the position of the wavefront at different
stages as it travels to and fro in the pipeline.
92
Valve closure in a simplified system
U/s of valve
–aVo/g
Mid-point
Ho of pipeline
Time
Fig. 7.4. Time histories of inertial head in pipeline after valve closure
93
Pressure transients in water engineering
the inlet to the pipeline from the reservoir. These effects will pro-
gressively erode the transient event until eventually a new steady
flow regime prevails. In theory there will always be some small residual
vestige of the transient with truly steady flow never actually being
reached. In practice the flow may be declared steady when amplitude
of the residual oscillation has become smaller than some selected
finite value.
It will have been noted that compression waves have travelled in
both the positive and negative directions of motion within the pipeline.
The same applies to rarefaction waves. Direction of wave motion, for
example in the positive x direction, does not imply that the wave
produces an increase in pressure.
94
8
Actual pipelines
8.1 Attenuation
Before proceeding to consider specific examples, some discussion is
appropriate regarding the phenomenon of attenuation which occurs
in pipelines generally but whose existence may not be obviously
apparent except where the pipeline is relatively long. Consider the
circumstance of a valve closed abruptly at the downstream end of a
long pipeline having uniform properties as shown in Fig. 8.1. Steady
flow velocity is Vo , pipe diameter is D, uniform wave speed is a and
dx=dt ¼ a for simplicity. Friction factor f is taken to be constant for
convenience and the effects of pressure change on velocity under
steady flow are ignored. Uniform gradient of piezometric line is So
under the original steady flow. Datum for head measurement is taken
to be at the level of the steady flow piezometric line at the valve as
shown in Fig. 8.1.
Head rise Hi just upstream of the valve when the valve is just closed at
time to ¼ 0 is aVo =g. At some time t after closure, the compression wave
front caused by valve closure will have propagated upstream to a point
95
Pressure transients in water engineering
V = Vo
Transient piezometric line or
hydraulic gradient
Hv
Hx
Compression wave
H
Steady flow
Hi = aVo/g
Piezometric line
Gradient = So
Horizontal datum
Initial steady
velocity = Vo
X
L
C– characteristic
J– = –Vo – g Ú S dt
J+ = Vo + gSot
t
C+ characteristic
J+ = Vo + gSot/2
t/2
x/2
Time to
x
96
Actual pipelines
and
ð
1
H ¼ ð Jþ JÞ=ð2g=aÞ ¼ a=g Vo þ gSo t Vo g S dt
2
ð
¼ aVo =g þ a=2 So t þ S dt ð8:4Þ
Total head H is made up of the initial steady flow head plus the inertia
head rise Hx , so that Hx ¼ H So at or Hx ¼ H So x. Then:
ð
Hx ¼ aVo =g þ a=2 So t þ S dt So at
ð
¼ aVo =g þ a=2 S dt So t ð8:5Þ
97
Pressure transients in water engineering
Using the Darcy equation to evaluate head loss due to friction then,
So ¼ f Vo2 =ð2gDÞ and S ¼ f V 2 =ð2gDÞ. Substituting in the expression
for velocity V at the wavefront:
ð
2 2
V ¼ f=ð4DÞ Vo t V dt ð8:6Þ
giving:
V ¼ ft=ð4DÞðVo2 V2 =3Þ
or, solving for V:
p
V ¼ 6D=ð ftÞð f1 þ 3½Vo ft=ð6DÞ2 g 1Þ ð8:7Þ
A simpler expression for V can be found by assuming that
3½Vo ft=ð6DÞ2 1. Using a series expansion and truncating at the
first term:
98
Actual pipelines
y2 þ 2y þ 1 ¼ 1 þ 3y2
rearranging then,
2y2 2y ¼ 0 or y ¼ 1
giving,
Vo ft=ð6DÞ ¼ 1 or t ¼ 6D=ðVo fÞ
Since x ¼ at therefore,
x ¼ 6aD=ðVo fÞ ð8:10Þ
The inertial head rise Hx at the wavefront is given by the expression,
Hx ¼ a=g½Vo þ ft=ð4DÞðV 2 =3 Vo2 Þ
When V ¼ Vo and t ¼ 6D=ðVo fÞ then,
Hx ¼ a=g½Vo þ ft=ð4DÞðVo2 =3 Vo2 Þ
¼ a=gðVo ft=ð4DÞ 23 Vo2 Þ
or as expected,
Hx ¼ a=g½Vo f6D=ð4DfVo Þ 23 Vo2 ¼ 0 ð8:11Þ
Confirming that the wave front is of zero height when V ¼ Vo .
To gauge the magnitude of x consider the following values: Vo ¼ 1 m/s;
f ¼ 0:02; D ¼ 1 m and (i) a ¼ 1200 m/s and (ii) a ¼ 300 m/s. Sub-
stituting in equation (8.10), then (i) x ¼ 120 km and (ii) x ¼ 30 km.
Only in very long pipelines is it likely that the wave will be completely
attenuated. In shorter lines significant attenuation can still be
observed.
99
Pressure transients in water engineering
100
Actual pipelines
Compression wave
160
140 t = 18.41 s
120
t = 13.81 s
100 t = 9.2 s
t = 4.6 s
Head (mAD)
80
60
Steady flow profile
t=0s
40
Downstream valve
Pipeline profile
20 Um Ghafa Water Project
3 s closure of downstream valve
instantaneous hydraulic profiles 4.6 £ t £ 23.01 s
0
0
790.3
1580.6
2370.8
3161.1
3951.4
4741.7
5532
6322.2
7112.5
7902.8
8693.1
9483.4
10 273.6
11 083.9
11 854.2
12 644.5
13 434.8
14 225
15 015.3
15 805.6
16 595.9
17 386.2
18 176.4
18 986.7
19 757
Chainage (metres from upstream reservoir)
180
140
80
t = 46.02 s t = 41.42 s
60
40 Closed valve
Um Ghafa Water Project
20 3 s closure of downstream valve
instantaneous hydraulic profiles
27.61 £ t £ 46.02 s
0
0
790.3
1580.6
2370.8
3161.1
3951.4
4741.7
5532
6322.2
7112.5
7902.8
8693.1
9483.4
10 273.6
11 083.9
11 854.2
12 644.5
13 434.8
14 225
15 015.3
15 805.6
16 595.9
17 386.2
18 176.4
18 986.7
19 757
101
Pressure transients in water engineering
80
70
t = 78.23 s
Pipeline profile
60
50
Head (mAD)
t = 73.63 s
t = 55.22 s
40
t = 59.83 s
30
t = 69.03 s
20
t = 59.83 s
Um Ghafa Water Project Closed valve
10 3 s closure of downstream valve t = 55.22 s
instantaneous hydraulic profiles
55.22 £ t £ 78.23 s
0
0
790.3
1580.6
2370.8
3161.1
3951.4
4741.7
5532
6322.2
7112.5
7902.8
8693.1
9483.4
10 273.6
11 083.9
11 854.2
12 644.5
13 434.8
14 225
15 015.3
15 805.6
16 595.9
17 386.2
18 176.4
18 986.7
19 757
Chainage (metres from upstream reservoir)
120
t = 96.64 s
100 t = 115.05 s
t = 87.44 s
80
Head (mAD)
t = 119.65 s
60
t = 82.84 s
40 Pipeline profile
102
Actual pipelines
103
Pressure transients in water engineering
now ceases to fall. After 2L=a seconds the reflected wave from the valve
travels upstream, reducing pressure along the majority of the pipeline as
shown in Fig. 8.2b at times 41.42 and 46.02 s. However, head does not
remain at vapour pressure near to the valve. Instead an initial hydraulic
gradient equal to the slope of the pipeline is established along these
parts of the line subject to vapour pressure. A quasi free surface flow
is established towards the valve, causing pressure to increase gradually
at the valve and in adjoining parts of the pipeline as shown in Fig. 8.2c
at times 55.22, 59.83 and 69.03 s.
When the rarefaction wave reaches the upstream reservoir after time
3L=a seconds it is reflected into the pipeline in the form of a compres-
sion wave producing a recovery of head. A situation is developed in
which head is increasing along upstream parts of the pipeline as
shown for times 55.22, 59.83 and 69.03 s. At the same time head is
rising over downstream parts near to the closed valve (Fig. 8.2c),
with vapour bubbles being removed as pressure rises. Eventually these
compression waves meet at a point within the system as the final
vapour cavity is removed and a head rise occurs at this point as
shown in Fig. 8.2c at time 73.63 s.
This head rise then spreads both upstream and downstream, pro-
ducing a more general head rise along the pipeline as illustrated in
180
Um Ghafa Water Project
3 s closure of downstream valve
160 time histories of head
140
100
80
60
8 km from upstream
40 reservoir
20
0
0.354
4.248
8.142
12.036
15.93
19/824
23.718
27.612
31.508
35.4
39.294
43.188
47.082
50.976
54.87
58.764
62.658
66.552
70.446
74.34
78.234
82.128
86.022
89.916
93.81
97.704
101.596
105.492
109.386
113.279
117.173
Time (s)
104
Actual pipelines
140
Um Ghafa Water Project
2 min closure of downstream valve
120 time histories of head
100
Head (mAD)
80
60
8 km from upstream
40 reservoir Upstream of valve
20
0
0.354
7.06
13.806
20.532
27.258
33.964
40.71
47.436
54.162
60.888
67.614
74.34
81.066
87.792
94.518
101.244
107.97
114.695
121.421
128.147
134.873
141.599
148.326
155.052
161.778
166.504
175.23
181.956
188.682
195.400
202.134
208.80
Time (s)
Fig. 8.4. Head plotted against time for 2 min valve closure
Fig. 8.2d at times 82.84, 87.44 and 96.64 s. Thereafter head starts to
decline as shown for times 115.05 and 119.65 s.
In deciding on an appropriate time of valve closure it is important to
consider not only the initial head rise at the valve on closure but also
the potential for unacceptable pressures being developed following
reflection of the initial pressure wave.
Figure 8.3 shows predicted variations of piezometric level at chainage
8 km from the upstream reservoir and just upstream of the valve for the
3 s valve closure. Figure 8.4 shows corresponding predictions for a more
extended valve closure lasting 2 min.
105
9
Valve operations
The term valve encompasses many different types, each with a wide
range of duties. Valves discussed in this chapter are those located
on a pipeline and through which flow passes during normal opera-
tion. These are often used to control discharge and/or pressure.
Many valves are under operator control using actuators to open and
close the valve at a pre-defined rate. Others may be manually oper-
ated using a handwheel. Still other valves are fitted with control
systems which respond to changes in flow and pressure throughout
the day. These valves may be set to maintain a set upstream or
downstream pressure or to produce a set range of downstream
pressure and flow. Each of these valve patterns has the capability to
initiate hydraulic transients during the process of flow adjustment.
Where a valve is not moving it will merely act as a partial trans-
mitter/reflector of surge effects arising elsewhere in the pipeline
system.
This chapter includes several examples of gravity pipeline systems
in which flow is controlled using valves. Valve operation is one of the
primary sources of pressure transient behaviour. In the majority of
cases discharge is controlled by a valve at the downstream end of the
main. A valve in this position may be called a ‘terminal discharge
valve’. Pressure surges will be developed in the line upstream of the
valve as it is opened or closed. In some instances a valve may be
positioned at some point within a pipeline system with appreciable
lengths of line both upstream and downstream of the valve. A valve
in this position may be called an ‘in-line valve’ with transient effects
produced both upstream and downstream of the valve. Other valves
are the discharge valves just downstream of pumps and the isolating
valves distributed along long pipelines.
106
Valve operations
130
TWL Flow rate = 95.47 Mld when
upstream reservoir is at
120 top water level (TWL)
LDO
110
Flow rate = 90.92 Mld when
100 upstream reservoir is at
lowest draw-off (LDO)
90
Head (mAD)
80
70
60
50
40
30
Pipeline profile
20
10
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
Chainage (km)
107
Pressure transients in water engineering
10 000.0
Valve head loss characteristics
1000.0
100.0
Valve head loss coefficient, KL
Submerged discharge
10.0 Butterfly
Gate
Needle
1.0
0.1
0.01
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
Valve position as a proportion of port stroke
The head loss characteristic of the butterfly valve and of some other
valves was as shown in Fig. 9.2, where head loss H ¼ KL V 2 =ð2gÞ. It
will be noted that changes in KL are modest for much of the stroke
from open to shut. Only when the valve has become about 15%
closed does the value of KL start to change more rapidly. Unfortunately
the availability of head loss data is more limited towards the more
interesting closed valve end of the stroke. Typically, data curves
obtained from manufacturers terminate around the 10% open position
and it is often necessary to complete the curve down to the closed
108
Valve operations
125
Start of valve operation
120
115
Head (mAOD)
Mid-point of main
110
Upstream of
butterfly valve
105
95
0.365
3.285
6.205
9.125
12.045
14.965
17.885
20.805
23.725
26.645
29.565
32.485
35.405
38.325
41.245
44.165
47.085
50.005
52.925
55.845
58.785
61.685
64.605
67.525
70.445
73.365
76.285
79.205
82.125
85.045
Time (s)
109
Pressure transients in water engineering
1.4
At butterfly valve
1.2
1.0
0.8
Velocity (m/s)
0.6
–0.2
0.365
9.125
17.885
26.645
35.405
44.165
52.925
61.685
70.445
79.205
87.965
96.725
105.485
114.245
123.005
131.765
140.525
149.285
158.045
166.805
175.585
184.326
193.086
201.846
210.606
219.366
228.126
236.886
245.647
254.407
263.167
271.926
Time (s)
operation of air valves at high points along the pipeline. Figure 9.5
shows the curves of minimum and maximum head during this valve-
opening sequence.
With computer simulation carried out for a sufficient length of time,
essentially steady flow conditions will have been established at the end
140
Maximum head
120
100
Minimum head
Head (mAOD)
80
60
40 Pipeline profile
20
Maximum and minimum head in
DN 1000 main 1 min butterfly valve opening
0
0
835
1671
2508
3342
4177
5012
5848
6401
6968
7818
8669
9235
10 095
10 668
11 529
12 390
13 251
14 112
14 701
15 585
16 469
Chainage (m)
110
Valve operations
200
180
160
120
100
40
DN 100 steel gravity main
20 3 min butterfly valve closure
0
0.365
12.045
23.725
35.405
47.085
58.765
70.445
82.125
93.805
105.485
117.165
128.845
140.525
152.205
163.885
175.585
187.246
196.926
210.606
222.286
233.966
245.647
257.327
269.006
280.686
292.368
304.045
315.725
327.405
339.064
350.764
362.444
Time (s)
111
Pressure transients in water engineering
200
180
Maximum head
160 Maximum
Maximum and minimum head during pressure head
140 3 min butterfly valve closure
Head (mAOD)
120
Minimum head
100
80
60
40
18
+
17
+
16
Pressure (bar g)
15
+
Maximum pressure on
14 float valve closure
+
13
+
12
0 2 4 6 8 10 12 14 16 18
Time (min)
112
Valve operations
similar levels obtained during float valve closure requires a valve stroke
time >13 min. Closure intervals can become quite prolonged. This is a
consequence of the need to slow valve movement during the final 10—
20% or so of closure in order to reduce the deceleration rate. This same
slow valve movement also takes place over the larger part of the stroke
when little flow change occurs.
113
Pressure transients in water engineering
Spillway
Penstocks
Valve house
Draw
off
Downstream
Dam dam face
crest
120 m
Draw-off well
Flow
regulating valve
Emergency Valve
Section through dam
butterfly house
valve
Burst site
114
Valve operations
Weight
Lever
Valve door
Flow
115
Pressure transients in water engineering
30
30 s emergency valve closure
25
20 At base of At emergency
draw-off well butterfly valve
Velocity (m/s)
15
5 Valve shut
–0.5
0.02
1.04
2.06
3.08
4.10
5.12
6.14
7.16
8.18
9.20
10.22
11.24
12.26
13.28
14.30
15.32
16.34
17.36
18.38
19.40
20.42
21.44
22.46
23.48
24.50
25.52
26.54
27.56
28.58
29.60
30.621
31.641
Time (s)
140
Base of
120
draw-off well
100
80
60
30 s emergency valve closure
40 Upstream of butterfly valve
20
0
0.02
1.00
1.98
2.96
3.94
4.92
5.90
6.88
7.86
8.84
9.82
10.80
11.78
12.76
13.74
14.72
15.70
16.68
17.66
18.64
19.62
20.60
21.58
22.56
23.54
24.52
25.50
26.48
27.46
28.44
29.42
Time (s)
116
Valve operations
Actuator
Water surface
Spillway
Flowmeter
Inlet bend
Discharge sump
Vertical
downpipe
Corner fillets
Spindle
Baseplate
117
Pressure transients in water engineering
where the energy of flow is dissipated. The port shape can be modified
to suit the circumstances both of the system in which it is installed and
also the range of flows over which the valve is to exercise control. The
port shapes can be altered over time to suit any changes in flow rate as a
scheme is developed. For example, in the early years of a project, flow
may be modest to suit prevailing demands while in later years flow
may have to be increased to meet anticipated increase in demand.
118
Valve operations
Gradient steepens
post-burst
Isolating valve
119
Pressure transients in water engineering
Inlet main
BPC
Receiving
reservoir Flow direction
M
Flow direction
M
Suction tank
PS (idle)
180
160
140
Head (mAOD)
120
100
80
0 1 2 3
Time (min)
Fig. 9.16. Head upstream of needle valve for 1.8 min closure
120
Valve operations
45
42
39
Head (mAOD)
36
33
30
Head downstream of needle valve –
1.8 min closure at a flow rate of 16 Mld
0 1 2 3
Time (min)
Fig. 9.17. Head downstream of needle valve for 1.8 min closure
type whose head-loss characteristic is shown in Fig. 9.2 and which is
installed in a chamber as shown in Fig. 9.15.
When the valve is closed over a time of 110 s, head upstream of the
valve increases smoothly to a maximum as depicted in Fig. 9.16. After
wave reflection from the BPC, head falls to produce an irregular
oscillation. On the downstream side of the valve, head falls during
valve closure (Fig. 9.17). After valve closure and wave reflection, head
largely recovers and a fairly regular oscillation is developed. The timescale
of oscillations shown in Figs 9.16 and 9.17 are indicative of the respective
lengths of pipeline upstream and downstream of the closed in-valve.
Overall variations in head through the system can be seen in Fig. 9.18.
175
BPC
150
Envelope curves of max. and
125 min. head for 1.8 min closure
Head (mAOD)
75
Downstream
reservoir
50
Needle valve chamber
25
Pipeline profile
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14
Chainage (km)
121
Pressure transients in water engineering
122
Valve operations
Butterfly valve 21
Gate valve 39
Cone valve 49
Globe valve 52
123
Pressure transients in water engineering
124
10
Pumps
Pumps are used to move fluids from point to point where insufficient
piezometric head is available to permit gravity flow. Pumping equip-
ment comes in a wide variety of forms to suit a range of different
functions. The potential for pressure transients to be initiated by each
type of pump is dependent upon the characteristics of the pump. The
majority of pumpsets in service initiate hydraulic transient events
when starting or stopping while other types of pump are themselves
incapable of initiating such events. Some designs of pumping plant
create transient behaviour as a part of their operation and in one
case the pump only functions by creating a hydraulic transient and
using its properties to lift water. In any discussion of hydraulic transients
it is therefore worthwhile to spend some time considering the types of
pump available and their characteristics as these will have a bearing
upon the nature of transient events developed and also on the types
of surge protection equipment appropriate in each instance.
125
Pressure transients in water engineering
126
Pumps
Ns = 20–35 Ns = 35–60
Radial Radial
Ns = 60–90 Ns = 90–120
Francis/semi-axial/diagonal Francis/semi-axial or diagonal
Ns = 120–160 Ns = 160–300
Mixed flow Axial/propeller
127
Pressure transients in water engineering
128
Pumps
100
90
Q > 631 litre/s
Q>
= 631 litre/s
80
Q = 315 litre/s Q = 63.1 litre/s
Q = 31.5 litre/s
Efficiency (%)
Q = 12.6 litre/s
70
Q = 6.3 litre/s
60
50
40
10 20 30 40 50 60 70 80 90100 200 300
Specific speed (rpm ÷(m3/s)/m3/4)
129
Pressure transients in water engineering
passed through fine screens, radial pumps can be used with efficiencies
reaching 90%.
Other considerations other than pumping efficiency can be impor-
tant, for example in sewage pumping. When required to pass solid
material, the impeller shape can be modified to avoid choking and/or
to disintegrate matter present in the sewage. Some efficiency is
sacrificed and such pumps may have peak efficiencies in the 50—60%
range. Sewage pumps will usually operate at relatively low speed,
typically not greater than 960 rpm. This compares with speed of a
clean water pump which will usually be 1460 rpm and upwards. An
unchokable pump may be classified according to the size of solid
object which it can pass, with the normal maximum being specified
as 100 mm. For very small flows, such as from isolated properties, it is
not practical to make a centrifugal pump capable of handling the
H, P H, P
H–Q H–Q
P–Q
P–Q
Q Q
(a) (b)
H, P H, P
H–Q
H–Q
P–Q
P–Q
Q Q
(c) (d)
Fig. 10.3. Performance curve forms: (a) unstable, overloading; (b) stable, non-
overloading; (c) high shut-off head, flat power curve; (d) high shut-off head, high
shut-off power
130
Pumps
Hp
Duty point
Speed N reducing
131
Pressure transients in water engineering
Hp/N 2
Q/N
2
Fig. 10.5. Hp /N against Q/N curve
132
Pumps
where dE=dt ¼ power ðPrÞ is the rate at which energy changes with
time, I is the pumpset ‘moment of inertia’ and ! is the angular velocity
of rotation. Differentiating with respect to time then:
Pr ¼ 12 I d!2 =dt ¼ I! d!=dt ð10:5aÞ
Equation (10.5) can be written in terms of torque (T) where torque is
Pr=! so that:
T ¼ I d!=dt ð10:5bÞ
When the pump is accelerating, say when being started, d!=dt ¼ þve,
and when the pump is decelerating, as after being tripped,
d!=dt ¼ ve. After being tripped power input becomes zero. (Note:
for anyone unfamiliar with the term ‘moment of inertia’, a brief descrip-
tion of this parameter and its determination has been included in the
appendix to Chapter 11.)
The Hp =N2 —Q=N and power P=N3 or torque T=N2 —Q=N curves
can be digitised to yield either linear or parabolic segments for each
zone of pump operation. Using the head Hp , speed N and flow Q
relationship for the normal zone of operation as an example, Fig. 10.6
illustrates how a series of overlapping parabolic segments can be used
to represent the performance curves.
Angular acceleration d!=dt and torque T are related through the
moment of inertia I using the equation:
T ¼ I d!=dt ð10:5bÞ
For a parabolic approximation to a segment of curve:
Hp =N2 ¼ aH ðV=NÞ2 þ bH V=N þ cH ð10:6aÞ
T=N2 ¼ aT ðV=NÞ2 þ bT V=N þ cT ð10:7Þ
where aH , bH , cH and aT , bT , cT are the coefficients of the parabolic
segments relating to Hp , N and V and T, N and V. These coefficients
can be pre-computed before any simulation of pressure transients
commences, with the values stored and used during all studies involving
the particular performance curves.
At the start of any time increment, the velocity V or flow rate Q at the
pump and the speed of the pump N will be known. The value of Q=N can
be used to select the appropriate parabolic segment of curve and hence
the coefficients aT , bT and cT of the parabola. Torque T can then be
found at the start of the time increment. Speed change ! over the
time step t can then be found using:
! ¼ tT=I or ! ¼ !o þ tT=I ð10:5cÞ
133
Pressure transients in water engineering
Parabola centred
Hp/N 2 on I – 1
I–1
(Hp/N 2)I–1
I
(Hp/N 2)I Parabola centred
on I + 1
I+1
(Hp/N 2)I+1
Q/N
(Q/N)I–1 (Q/N )I (Q/N)I+1
using the revised speed ! and initial velocity, the coefficients aH , bH and
cH of the parabola defining Hp can be obtained and hence the value of
pumping head at the close of the time increment.
Total moment of inertia I is made up of the individual inertias for the
motor, shafting, pump impeller and entrained liquid. Values of these
components can usually be obtained from pump and motor manufacturers
when a pumpset has been selected. For purposes of preliminary studies
before a particular pump supplier has been selected it may be necessary
to obtain a preliminary estimate of moment of inertia. Linton (1972)
has compiled data for both pump and electric motor inertias and presented
these in a graphical form covering a wide range of flow, head and speed.
These curves allow initial estimates of inertia to be obtained.
134
Pumps
Hp Piezometric line or
hydraulic gradient
Hr
M
H
Rising main
Tank or reservoir V +
Pump
Horizontal datum
J± = V ± g/aH
Dt
C± characteristic
Dx
Also, from the characteristic path arriving at the pump at the end of
the time increment, the quasi-invariant value J from Chapter 4,
gives the relationship just downstream of the pump:
V g=aH ¼ J or H ¼ ð J VÞ=ðg=aÞ
Substituting for H then:
ðÞa=gðJ VÞ Hr ¼ aH V 2 þ bH VN þ cH N2
or, rearranging:
aH V 2 þ fbH N þ ðÞa=ggV þ fcH N2 þ Hr ðÞa=g Jg ¼ 0
ð10:8Þ
This quadratic equation (10.8) can be solved for V. The corresponding
value of H at the end of the time increment can then be obtained
from:
H ¼ ð J VÞ=ðg=aÞ
135
Pressure transients in water engineering
Hp
Piezometric line or
hydraulic gradient
Hd
Booster pump
Delivery main
Suction main Hu
V +
+
Horizontal datum
J+ = V + g/aHu J– = V – g/aHd
Dt
C+ characteristic C– characteristic
Dxu Dxd
136
Pumps
137
Pressure transients in water engineering
Screens
Upstream river
Intake weir
Bellmouth inlet
Flap gates
Penstocks
4 off concrete
volute pumps
Tunnel length
5.3 km
The same equation for the pumps is used but with the order of
upstream and downstream head reversed, that is:
Hd=s Hu=s ¼ aH V 2 þ bH VN þ cH N2
138
Pumps
139
Pressure transients in water engineering
The highest system curve will embody maximum static head and
maximum system resistance representing a mature main. In carrying
out surge analyses, usually the amplitude of pressure transient is a
function of flow rate and so computer simulation should also employ
the minimum assessment of system resistance to simulate clean main
conditions and also minimum static lift. The pump will deliver a
greater discharge when mains are new and static is low than when
mains are mature and static is high. Variation of system resistance is
particularly important in raw water and sewage lines. Resistance
changes will usually be less significant for treated water mains where
a pipe lining is included.
140
Pumps
motors as in start mode the motor may not develop sufficient torque to
allow the pump to accelerate. An approximate limit of 30 kw has been
quoted as a maximum for Start/Delta starting.
Change-over from Star to Delta can cause problems. Timing is very
important and also the method of switching. Open switching completely
breaks the circuit during transition. This break, and the subsequent
switch to Delta mode can produce high transient currents and large
negative torques. These currents may damage windings and the
torques can damage the shaft between motor and pump. A flexible
coupling will alleviate this effect.
Where Star/Delta is unsuitable, transformer starting may well be an
option although more costly. The supply voltage in the initial stage is
reduced to a percentage of mains supply voltage through use of a trans-
former. Typically three tappings are provided on the transformer to
allow for difference between actual and anticipated values. These
tappings may represent voltages in the range 60—80% of mains
voltage. As the motor approaches its maximum speed the supply is
switched onto mains voltage. Complicated switching arrangements
reduce transient currents during switching.
141
Pressure transients in water engineering
Carlops rising main – start of duty pump. Head at pumping station and mid-point of main
400
390
380
370
Head (mAOD)
360
350
340
330 Ch. –0 km
Ch. –3.4 km
320
0.055
2.255
4.455
6.655
8.855
11.055
13.255
15.455
17.655
19.855
22.055
24.255
26.455
28.655
30.855
33.055
35.255
37.455
39.655
41.855
44.055
46.255
48.455
50.655
52.855
55.055
57.255
59.455
61.655
63.855
Time (s)
pressure which will be close to the pump’s shut-valve delivery head. Initi-
ally, head in the rising main will exceed the design head. Figure 10.10
shows this for a solo duty pump delivering into a single rising main
carrying treated water. After a delay while the startup compression
pressure wave travels along the main, points along the pipeline are
affected by the pump start and experience a head rise similar to that
at the pump. Figure 10.10 also shows head conditions at the mid-
point of the pipeline, chainage 3.4 km from the pumping station.
After time L=a the pressure wave reaches the downstream end of the
pipeline and a rarefaction wave reflection travels back along the main
towards the pumping station, bringing with it a relief of pressure.
Points further from the pumping station experience this relief of
pressure first, with a fall in head occurring at the mid-point of the pipe-
line. When the reflected pressure wave reaches the pumping station a
fall in head occurs to a lesser extent as the pump operating point
moves to a lower delivery head and higher flow. The pressure waves
continue to travel to and fro in the main, with each wave reflection
producing a fall in pump delivery head and corresponding increase in
flow. Effects of developing pipeline resistance as flow increases gradually
diminish the transient effect, with an almost steady flow occurring after
142
Pumps
Carlops rising main. Start of duty pump. Max. and min. head along main
400
350
300
Elevation (mAOD)
250
200
150
100
Invert level (i.l.) (mAD)
h (max.)
50
h (min.)
0
0
280
560
840
1120
1400
1680
1960
2240
2520
2800
3080
3360
3640
3920
4200
4480
4760
5040
5320
5600
5880
6160
6440
6720
7000
Chainage (m)
143
Pressure transients in water engineering
60
50
Head (mASL)
40
30
20
10 PS
Ch. 1525 m
0
0.173
3.114
6.055
8.996
11.937
14.878
17.819
20.760
23.701
26.642
29.583
32.524
35.465
38.406
41.347
44.288
47.229
50.170
53.111
56.052
58.993
61.934
64.875
67.816
70.757
73.698
76.639
79.580
82.521
85.462
Time (s)
pressure wave reflection from the bifurcation into two mains down-
stream of which there is an overall increase in cross-sectional area.
Thereafter, head starts to rise more slowly. About halfway along the
mains, at chainage 1.525 km, head rises when the pressure wave
reaches this location (Fig. 10.12).
Following reflection of this start-up pressure wave from the outfall
end of the system, head falls within the mains. Those points closest
to the outfall experience the effects of pressure relief first, as at the
mid-point of the mains (Fig. 10.12), with the pumping station being
affected after about 21 s. At 30 s the second pump is operated while
pressure is relatively low at the pumping station, and a second upsurge
pressure wave is created. It will be noted that this second upsurge is
smaller than the first. This is common in multi-pump installations, with
successive pump starts having a progressively smaller effect.
Subsequent wave reflections progressively reduce system pressures
with near-steady conditions achieved after about 112 min of starting
the first pump. The difference in head between the pumping station
and the halfway point is indicative of pipeline resistance.
Development of flow in the two mains follows a stepped pattern as
shown in Fig. 10.13. Velocity in each main is initially similar but as
144
Pumps
2.5
Velocity (m/s)
1.5
0.5
Ch. 0450
Ch. 0600
0
0.173
2.941
5.709
8.477
11.245
14.013
16.781
19.549
22.317
25.085
27.853
30.621
33.389
36.157
38.925
41.693
44.461
47.229
49.997
52.765
55.533
58.301
61.069
63.837
66.605
69.373
72.141
74.909
77.677
80.445
83.213
85.981
Time (s)
60
50
40
Elevation (mASL)
30
i.l. (mASL)
20 h (max.)
h (min.)
10
0
0
122
244
366
488
610
732
854
976
1098
1220
1342
1464
1586
1708
1830
1952
2074
2196
2318
2440
2501
2623
2745
2867
2989
3111
3233
–10
Chainage (m)
Fig. 10.14. Envelope curves for direct start of two pumps in sequence
145
Pressure transients in water engineering
(a) Prediction of start-up transient pressures and time taken for these
transient pressures to decay in order to establish essentially steady
flow.
(b) Determining an appropriate time delay to be allowed between
successive pump starts in a multi-pump installation to avoid any
adverse interaction between surging produced by each pump.
(c) Providing a convenient means of establishing steady flow condi-
tions throughout the network.
146
Pumps
147
Pressure transients in water engineering
Az Zour to Wafra Main. Trip of two pumps. Max. and min. head along main
400
i.l. (mPWD)
350
h (max.)
h (min.)
300
Head (mPWD)
250
200
150
100
50
0
0
1153.36
2367.56
3532.11
4751.58
5983.65
7155.67
8448.01
9605.88
10 823.28
11 979.81
13 198.89
14 418.39
15 634.59
16 858.59
18 029.24
19 328.61
20 559.99
21 843.30
23 005.95
24 300.18
25 532.26
26 754.46
28 037.77
29 262.87
30 552.27
31 841.67
33 075.07
34 372.03
35 607.23
Chainage (m)
Fig. 10.15. Envelope curves for long rising main after pump failure
pumps are tripped, delivery head falls steeply until discharge head is
below the suction well level. Under the action of the suction well
level, sewage continues to flow through the pump, keeping the check
valves partially opened. The head downstream of the pumps is largely
Az Zour to Wafra Rising Main. Trip of two duty pumps. Head at PS, ch. 5 km and ch. 10 km
250
D/s pumps
Ch. 5 km
200
Ch. 10 km
Head (mPWD)
150
100
50
0
0.056
1.344
2.632
3.920
5.208
6.496
7.784
9.072
10.360
11.648
12.936
14.224
15.512
16.800
18.088
19.376
20.664
21.952
23.240
24.528
25.816
27.104
28.392
29.680
30.968
32.256
33.544
34.832
36.120
37.408
38.696
Time (s)
Fig. 10.16. Head variations in a long rising main after pump failure
148
Pumps
Egaila, Kuwait. Simultaneous trip of three pumps. Head at PS and ch. 2.6 km
160
140
120
Head (mPWD)
100
80
60
40
Ch. –0 km
20
h –2.6 km
0
0.116
1.160
2.204
3.248
4.292
5.336
6.380
7.424
8.468
9.512
10.556
11.600
12.644
13.688
14.732
15.776
16.820
17.864
18.908
19.952
20.996
22.040
23.084
24.128
25.172
26.216
27.260
28.304
Time (s)
stabilised, showing only a slow and small recovery till around 14 s when
head increases more rapidly. Part way along the mains at chainage
2.6 km the fall in head is similarly halted. Figure 10.18 shows the corre-
sponding variation of velocity at the start of rising mains. Following an
Egaila, Kuwait. Trip of three pumps. Velocity at start of mains
1.6
1.4
1.2
1.0
Velocity (m/s)
0.8
0.6
0.4
0.2
Ch. 0.0 m
0 Time (s)
0.116
0.580
1.044
1.508
1.972
2.436
2.900
3.364
3.828
4.292
4.756
5.220
5.684
6.148
6.612
7.076
7.540
8.004
8.468
8.932
9.396
9.860
10.324
10.788
11.252
11.716
12.180
12.644
13.108
13.572
14.036
Time (s)
149
Pressure transients in water engineering
Egaila, Kuwait. Trip of three pumps. Max. and min. head along mains
160
i.l. (mPWD)
140
h (max.)
h (min.)
120
Elevation (mPWD)
100
80
60
40
20
0
0
200
400
600
800
1000
1100
1300
1500
1700
1900
2000
2200
2400
2600
2800
3000
3100
3300
3500
3700
3900
4000
4200
4400
4600
4800
4900
5100
Chainage (m)
Fig. 10.19. Envelope curves in a sewage rising main after pump failure
150
Pumps
Piston
Discharge
Suction
Time
(a)
Time
(b)
Time
(c)
Fig. 10.21. Flow against time relationships for piston pumps: (a) single cylinder
single-acting; (b) twin cylinder single-acting or single cylinder double-acting; (c)
three cylinder single-acting
151
Pressure transients in water engineering
P and Q
Flow Q against head Hp
Hp
152
Pumps
Air compressor
Compressed Compressed air
air receiver supply pipe
Discharge
pressure pipe
Isolating valve
Check valve
Ejector cylinder
Check valve
Elevation
Check valve
Check valve
Isolating valve
Isolating valve
Ejector cylinder
Isolating valve
Isolating valve
Check valve
Check valve
153
Pressure transients in water engineering
dx/dt = V – a
DH
H
dx/dt = V + a
Steady flow piezometric line
Datum
Vu/s Vd/s +
Qe Collector/rising main
Pumping station
154
Pumps
If the discharge from the ejector is Qe then with reference to Fig. 10.24
the equations required for solution at the ejector connection are:
from characteristics
Jþ ¼ Vu=s þ g=aH and J ¼ Vd=s g=aH
and from conservation of volume
Vu=s A þ Qe ¼ Vd=s A
eliminating Vu=s and Vd=s then,
H ¼ ð Jþ J þ Qe =AÞ=ð2g=aÞ ð10:11Þ
To obtain an idea of the magnitude of the inertial pressure rise caused
by an ejector, assume that there is no time for wave reflections to occur
prior to the maximum discharge from the ejector being achieved. Then
if Jþ and J remain at their initial steady flow values, the initial head
Ho at the ejector connection before discharge Qe occurs will be:
Ho ¼ ð Jþ JÞ=ð2g=aÞ
and when the maximum discharge from the ejector is Qemax , maximum
head will be:
Hmax ¼ Ho þ Qemax =ð2Ag=aÞ
or peak inertial head
H ¼ Qemax =ð2Ag=aÞ ð10:12Þ
Suppose Qemax ¼ 5 litres/s, D ¼ 300 mm (A ¼ 0:0707 m2 ) and a ¼
1000 m/s. The maximum inertial head increase will be:
H ¼ 0:005=ð2 0:0707 9:81=1000:0Þ ¼ 3:6 m
In a smaller main, say D ¼ 200 mm,
H ¼ 0:005=ð2 0:0314 9:81=1000:0Þ ¼ 8:1 m
The smallest sewage rising main may have diameter D ¼ 100 mm
giving: H ¼ 0:005=ð2 0:00785 9:81=1000:0Þ ¼ 32:4 m.
In a plastic rising main these inertia head values would be reduced
proportionately for the smaller wave speed in these pipes.
155
Pressure transients in water engineering
Header tank
Supply
reservoir
M
Air vessel
Riser
Non-return valve
Supply pipeline
Fast-closing valve
156
Pumps
H1
Piezometric levels
H3
H2
Qp Pumped fluid
2
Diffuser 3
Mixing tube
1
Qm Qp + Qm
Motive fluid
Nozzle
157
Pressure transients in water engineering
h, Z
Z–q
h–q
158
11
Flywheels
159
Pressure transients in water engineering
11.2 Flywheels
The flywheel consists of one or more metal disks rotating on the same
shaft as the pump and motor and with the same speed. Figure 11.1
shows a typical installation for a vertical pump.
On the delivery side of a pumping station after pump failure, the rate
at which head decreases and flow decelerates is reduced by adding
additional inertia. In the case of a booster station, on the suction side
of the station, rate of head rise and rate of flow deceleration are also
Flywheel assembly
Lower bearing
Coupling
Shaft extending
downwards to pump
160
Flywheels
reduced by adding a flywheel. After flow has ceased and the check valve
on each pump delivery branch has closed, the pump/flywheel arrange-
ment plays no further part in the transient event in a rising main. The
flywheel may have been added for the purpose of controlling either
minimum or maximum transient pressures in the main. Where peak
pressures are to be controlled then these will arise as a consequence of
deceleration of a reversed flow which has occurred subsequent to the
initial downsurge following pumping failure. The only way in which the
flywheel can limit such maximum pressures is by previously controlling
these initial minimum piezometric levels. The adverse hydraulic gradient
producing flow deceleration and accelerating flow in the reversed
direction in the main is thus limited. So it is only by limiting the
initial downsurge that the subsequent upsurge is controlled.
The flywheel forms an integral part of the pumpset assemblage and it
is speed change of the machine that initiates the transient event. By
adding the flywheel the rates of speed change are reduced and the
pressure transient is controlled at its source. The entire system benefits
from the flywheel solution. This contrasts with some other forms of
surge alleviation in which the hydraulic transient is allowed to
develop and its subsequent spread through the majority of a pumping
system is restricted by additional surge equipment. Reducing rates of
velocity change at the pump itself gives check valves time to respond
thus reducing potential for valve slam. The same applies during pump
start if the flywheel means that a slower acceleration of pump occurs.
The check valve will not then be thrown open.
This form of hydraulic transient suppression also has the advantage of
ease of maintenance.
Where a gearbox is included, relative speeds of pump and motor
require to be considered (Miller, 1978). Resultant inertia is expressed
in terms of pump speed so that
161
Pressure transients in water engineering
162
Flywheels
163
Pressure transients in water engineering
Lawhead/Silverburn rising main. Max. and min. head after trip of Uniglide pump
350
300
250
Elevation (mAOD)
200
150
i.l. (mAOD)
100 hmax/0.0
hmin/0.0
50 hmax/2.7
hmin/2.7
0
0
150
300
450
600
750
900
1050
1200
1350
1500
1650
1800
1950
2100
2250
2400
2550
2700
2850
3000
3150
Chainage (m)
moment of inertia ¼ 2.7 kg.m2 was found adequate for the task. Figure
11.2 shows that both maximum and minimum piezometric levels were
improved with minimum head remaining above the pipeline elevation.
In the case of the DUOGLIDE option it was necessary to have a
flywheel with moment of inertia ¼ 13.5 kg.m2 before a similar degree
Lawhead/Silverburn rising main. Max. and min. head after trip of Duoglide pump
350
300
250
Elevation (mAOD)
200
150
i.l. (mAOD)
hmax/0.0
100
hmin/0.0
hmax/13.5
50 hmin/13.5
0
0
150
300
450
600
750
900
1050
1200
1350
1500
1650
1800
1950
2100
2250
2400
2550
2700
2850
3000
3150
Chainage (m)
164
Flywheels
SDC/0.0
1.2 SDC/2.7
DRA/0.0
1.0 DRA/13.5
0.8
Velocity (m/s)
0.6
0.4
0.2
0
0.041
0.615
1.189
1.763
2.337
2.911
3.485
4.059
4.633
5.207
5.781
6.355
6.929
7.503
8.077
8.651
9.225
9.799
10.373
10.947
11.521
12.095
12.669
13.243
13.817
14.391
14.965
15.539
16.113
–0.2
Time (s)
165
Pressure transients in water engineering
300
250
Head (mAOD)
200
150
100 SDC/0.0
SDC/2.7
50 DRA/0.0
DRA/13.5
0
0.041
1.025
2.009
2.993
3.977
4.961
5.945
6.929
7.913
8.897
9.881
10.865
11.849
12.833
13.817
14.801
15.785
16.769
17.753
18.737
19.721
20.705
21.689
22.673
23.657
24.641
25.625
26.609
27.593
28.577
Time (s)
166
Flywheels
the DUOGLIDE. The pumps do not turbine once the flywheels have
been added, as minimum head is maintained at higher levels than
those within the suction well.
25
20
Head (mASL)
15
10
0
0.103
2.781
5.459
8.137
10.815
13.493
16.171
18.849
21.527
24.205
26.883
29.561
32.239
34.917
37.595
40.273
42.951
45.629
48.307
50.985
53.663
56.341
59.019
61.697
64.375
67.053
69.731
72.409
75.087
77.765
80.443
–5
Time (s)
167
Pressure transients in water engineering
1.5
Velocity (m/s)
1.0
0.5
0
0.103
2.060
4.017
5.974
7.931
9.888
11.845
13.802
15.759
17.716
19.673
21.630
23.587
25.544
27.501
29.458
31.415
33.372
35.329
37.286
39.243
41.200
43.157
45.114
47.071
49.028
50.985
52.942
54.899
–0.5
Time (s)
until it is below the sewage level in the suction well. At this time the
pumps begin to turbine and piezometric level downstream of the
pumps is largely stabilised by continued flow from the suction well.
This continues until around 40 s have elapsed when flow finally
ceases and check valves close. Flow reversal along the mains produces
a modest head rise and a fairly regular oscillation.
Figure 11.7 shows the corresponding velocity variation in the DN
700 line leaving the pumping station. Flow declines in a stepped
manner, with each wave reflection from the treatment works producing
a further reduction in flow until check valve closure.
Notwithstanding the stabilising effect of the suction well and pumps’
turbining, Fig. 11.8 shows that the entire length of pipelines up until the
summit at chainage 4321 m is subject to severe sub-atmospheric pres-
sure. Operation of the air valves at this high point alleviates conditions
between this summit and the treatment works.
The minimum pressures predicted were unacceptable and a flywheel
having inertia ¼ 255 kg.m2 was added to each pump. Subsequent
predictions after failure of the pumps produced a more gradual fall in
piezometric level downstream of the pumps (Fig. 11.6), with
minimum head remaining above suction well level. Flow reversal
occurred around 56 s with pumps tripped at time ¼ 1 s. After check
valves shut, subsequent peak pressures were below the steady
168
Flywheels
Sharjah PS No. 1 max. and min. head following pumping failure DN 450 main
30
i.l. (mASL)
25 hmax/0.0
hmin/0.0
hmax/255
20
hmin/255
Elevation (mASL)
15
10
0
0
301.5
603.0
904.5
1206.0
1507.5
1809.0
2110.5
2411.9
2713.4
3014.9
3316.4
3617.9
3919.4
4220.9
4422.3
4725.2
5028.0
5330.8
5633.6
5936.4
6144.9
-5
–10
Chainage (m)
Fig. 11.8. Envelope curves after pumping failure for DN 450 main
169
Pressure transients in water engineering
Sharjah PS No. 1 air volumes around ch. 4321 m after pumps are tripped
12
10
6
Volume (m3)
4 450/0.0
700/0.0
2 450/255
700/255
0
0.103
13.081
26.059
39.037
52.015
64.993
77.971
90.948
103.926
116.904
129.881
142.859
155.836
168.814
181.792
194.769
207.747
220.724
233.702
246.680
259.657
272.635
285.612
298.590
311.568
324.545
337.523
350.500
363.478
376.456
389.433
–2
Time (s)
will allow, this example has also illustrated how a flywheel can be
considered for larger schemes, acting in conjunction with other forms
of alleviation, to provide an overall package of protection measures.
170
Flywheels
the maximum number of operating pumps but may result from trip of a
lesser number of units.
171
Pressure transients in water engineering
Angle, dq
r dr
D/2
electric motor, the shaft connecting the motor and pump and the
impeller of the pump, possibly including entrained liquid.
Power P is rate of working or,
P ¼ dE=dt ¼ 12 I d!2 =dt ¼ I! d!=dt ðA10:2aÞ
If speed ! is increasing — that is, the machine is running up to full
speed — then d!=dt is þve and rotational kinetic energy is increasing,
and when d!=dt is ve the machine is decelerating and stored energy is
decreasing. The reduced kinetic energy of rotation is partially trans-
ferred to the liquid being pumped and so pumping is continued for a
time after a pump is switched off or tripped.
Torque T is defined as:
T ¼ P=! ¼ I d!=dt ðA10:2bÞ
and is an alternative to using power as a variable.
172
12
Pressure vessels
Pressure vessels form one of the most versatile means of providing surge
protection at pumping stations. A wide range in vessel capacity can be
provided, allowing the vessel or vessels to be optimised to suit any
system. This chapter describes the principles of operation of a pressure
vessel and provides some examples of installations.
All forms of hydraulic transient protection endeavour to limit the
rates of change of flow, acceleration or deceleration within at least
part of a pipeline network. Such flow changes may be introduced to a
system by starting or stopping of pumps. One method of limiting the
rate of acceleration/deceleration of flow is to provide some means of
augmenting flow into the pipeline. This is the principle of operation
of a pressure vessel when it is installed on the pipeline (Fig. 12.1).
A mass of gas, usually air or nitrogen, is contained in the upper part
of the vessel which is usually cylindrical in shape and installed
either nearly horizontally or vertically. Occasionally a sloping vessel
may be included such as at Ross Priory on Scottish Water’s Loch
Lomond Scheme where four vessels inclined at 458 are used
(Fig. 12.2).
A pressure vessel may be used to improve minimum head conditions
and/or maximum transient pressures.
173
Pressure transients in water engineering
Pressure vessel
Gas charge
M
Pipeline
174
Pressure vessels
Auto-sequence start of 24 MGD pumps No. 4 and 5. Surge vessel level from static
1.2
Sluice valve level
1.0
Water level (metres from static level)
0.8
0.6
0.4
0.2
0
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
–0.2
Time (min)
175
Pressure transients in water engineering
176
Pressure vessels
Pressure vessel
M
Small |dV/dt| downstream of
vessel connection
177
Pressure transients in water engineering
178
Pressure vessels
and the maximum expanded gas volume. To this maximum gas volume
is added a factor of safety in order to determine the gross vessel capacity
necessary.
Fixing minimum transient pressures is of particular importance as far
as treated water pipelines are concerned. Ideally the minimum transient
pressure should remain positive throughout the pipeline system for a
treated water application to avoid risk of contamination, but pressure
may be allowed to fall to a lower level in a raw water or sewage
system. Unfortunately it is a sad fact that in many countries economic
constraints make it difficult to find money for surge suppression, let
alone meet the strict criteria of maintaining positive pressures in treated
water supply pipelines. In such countries, circumstances dictate that
minimum transient pressures become sub-atmospheric with depressing
regularity given the often unreliable state of power supplies.
Having estimated a gas charge for a pressure vessel, more detailed
computations can be carried out. These calculations may commence
from static conditions with pumps being operated in sequence, for
multi-pump installations, to allow time for the hydraulic transient
developed from a pump start to dissipate before the next duty pump
is operated.
179
Pressure transients in water engineering
0.5
Velocity (m/s)
0
0.045
3.060
6.075
9.090
12.105
15.120
18.135
21.150
24.165
27.180
30.195
33.210
36.225
39.240
42.255
45.270
48.284
51.299
54.314
57.329
60.344
63.359
66.374
69.389
72.403
75.418
78.433
81.448
84.463
87.478
–0.5
pump 1
pump 2
–1 u/s of connection
riser
ch 1.5 km
–1.5
Time (s)
pipelines the pump output increases once more and some flow leaves
the vessel with pressure in the gas mass falling and the gas volume
expanding. After 30 s the second pump is operated and a further
period of pressurisation of the gas charge in the vessel occurs but to a
lesser extent than when the first pump was operated. After about
Abu Hamour PS 44. Sequenced pump start. Max. and min. head along mains
45
40
35
Elevation (mQNHD)
30
25
10
0
0
150
350
550
750
900
1100
1300
1500
1650
1850
2000
2200
2400
2550
2750
2900
3100
3300
3500
3700
3900
4100
4300
4500
4700
4904
Chainage (m)
Fig. 12.6. Envelope curves along a rising main after pump start
180
Pressure vessels
112 min essentially steady flow has been established in the pumping
system. It will be noted that the initial discharge from the first pump
exceeds its design steady output. This can have implications with
regard to the short-term power requirements for the pump. The overall
variation of piezometric level along the pumping mains is shown in
Fig. 12.6. Typically a vessel installation will produce a fairly uniform
variation of head along a main. As shown in the figure, inclusion of
the pressure vessel has produced a relatively smooth variation in max-
imum head.
1.5
Velocity (m/s)
0.5
0
0.016
1.568
3.120
4.672
6.224
7.776
9.328
10.880
12.432
13.984
15.536
17.088
18.640
20.192
21.744
23.296
24.848
26.400
27.952
29.504
31.057
32.609
34.160
35.712
37.264
38.816
40.368
41.920
43.472
45.024
46.576
–0.5
Time (s)
181
Pressure transients in water engineering
2.5
Air volume (m3)
1.5
1
n = 1.2
0.5 n = 1.0
n = 1.4
0
0.016
1.552
3.088
4.624
6.160
7.696
9.232
10.768
12.304
13.840
15.376
16.912
18.448
19.984
21.520
23.056
24.592
26.128
27.664
29.200
30.737
32.273
33.808
35.344
36.880
38.416
39.952
41.488
43.024
44.560
46.096
47.631
Time (s)
182
Pressure vessels
250
200
Head (mAOD)
150
100
n = 1.2
50 n = 1.0
n = 1.4
0
0.016
1.552
3.088
4.624
6.160
7.696
9.232
10.768
12.304
13.840
15.376
16.912
18.448
19.984
21.520
23.056
24.592
26.128
27.664
29.200
30.737
32.273
33.808
35.344
36.880
38.416
39.952
41.488
43.024
44.560
46.096
47.631
Time (s)
183
Pressure transients in water engineering
250
200
Head (mAOD)
150
100
50 No throttle
DN 150 bp
0
0.016
1.600
3.184
4.768
6.352
7.936
9.520
11.104
12.688
14.272
15.856
17.440
19.024
20.608
22.192
23.776
25.360
26.944
28.528
30.113
31.697
33.281
34.864
36.448
38.032
39.616
41.200
42.784
44.368
45.952
47.535
Time (s)
184
Pressure vessels
Total head
available
Head drop
in throttle
Piezometric line
without throttling
Pumping main
Pressure vessel
Pumping station
185
Pressure transients in water engineering
A throttling effect can be achieved in several ways but in all cases the
head drop H achieved by the throttle can be represented as:
X
H ¼ ðKL =A2 ÞQ2 =ð2gÞ ð12:6Þ
186
Pressure vessels
Do
(3.34 – 5.0) × Do
Do
2 × Do
Orifice may be shaped onto the
vessel base or held between pipe flanges (1.67 – 2.5) × Do
(a) (b)
Outflow from vessel Check valve closed Outflow from vessel Check valve
during inflow
(c) (d)
Helicoidal inflow
Spiral casing
Pressure vessel
Rising main
(e)
187
Pressure transients in water engineering
250
200
Elevation (mAOD)
150
Fig. 12.13. Envelope curves showing throttle effect after a pumping failure
188
Pressure vessels
189
Pressure transients in water engineering
60
50
Head (mBNSD)
40
30
20
Head (TWL)
10
Head (BWL)
0
0.090
8.910
17.730
26.550
35.370
44.190
53.010
61.830
70.650
79.469
88.289
97.109
105.928
114.748
123.568
132.387
141.207
150.027
158.846
167.666
176.486
185.305
194.125
202.944
211.764
220.584
229.403
238.223
247.043
255.862
264.682
Time (s)
0.8
0.6
0.4
0.2
0
0.090
8.910
17.730
26.550
35.370
44.190
53.010
61.830
70.650
79.469
88.289
97.109
105.928
114.748
123.568
132.387
141.207
150.027
158.846
167.666
176.486
185.305
194.125
202.944
211.764
220.584
229.403
238.223
247.043
255.862
264.682
–0.2
Time (s)
190
Pressure vessels
5
Volume (m3)
1 Vol. (BWL)
Vol. (TWL)
0
0.090
8.730
17.370
26.010
34.650
43.290
51.930
60.570
69.210
77.850
86.489
95.129
103.768
112.408
121.048
129.687
138.327
146.967
155.606
164.246
172.886
181.525
190.165
198.805
207.444
216.084
224.724
233.363
242.003
250.643
259.282
267.922
Time (s)
191
Pressure transients in water engineering
60
50
Elevation (mBNSD)
40
30
i.l. (mBNSD)
20 max. (BWL)
min. (BWL)
10 max. (TWL)
min. (TWL)
0
0
204
407
611
814
1018
1221
1425
1629
1832
2036
2239
2443
2646
2850
3053
3257
3461
3664
3868
4071
4275
4478
4682
4886
5089
5293
Chainage (m)
192
Pressure vessels
350
300
Elevation (mAOD)
250
200
150
100
i.l. (mAOD)
50
Hydraulic gradient
0
0
977
1954
2932
3909
4886
5863
6840
7818
8795
9772
10 749
11 726
12 703
13 681
14 658
15 635
16 612
17 589
18 567
19 544
20 195
20 883
21 968
23 053
24 139
25 224
26 309
27 026
28 090
Chainage (m)
193
Pressure transients in water engineering
considerably lower than that of the pumps serving the main to Fossdale
SR. This difference in pump discharge head for the two mains has a sig-
nificant role to play in determining hydraulic transient behaviour.
Pressure vessels are located downstream of both the Fossdale and
Aysgarth pumps at Bainbridge and also upstream and downstream of
Stonehouse pumping station. Figure 12.19 shows these vessels which
are fitted with insulation jackets to avoid freezing in winter.
194
Pressure vessels
350
300
250
Head (mAOD)
200
150
100 Suction
Aysgarth
50 Fossdale
0
0.270
9.180
18.090
27.000
35.910
44.820
53.730
62.640
71.550
80.460
89.370
98.280
107.190
116.099
125.009
133.919
142.830
151.740
160.650
169.560
178.470
187.380
196.290
205.201
214.111
223.021
231.931
240.841
249.751
258.661
267.571
Time (s)
but at a greater rate due to the smaller vessel on this line. When head
on the suction side of the station exceeds head downstream of the
Aysgarth pump, that pump starts to turbine with flow passing through
the pump from upstream. Suction head and downstream head follow
one another closely. Head downstream of the Fossdale pump remains
unaffected by the suction level.
A compression wave travels upstream in the suction main to Yarker
Bank and is reflected back to Bainbridge arriving after about 45 s when a
steep fall in suction head occurs. The fall in head causes flow into the
Aysgarth line to cease and its pump check valve to close. Head down-
stream of the Aysgarth pump then falls smoothly as its vessel air volume
expands. Subsequent further wave reflections in the long suction main,
coupled with falling head downstream of the Aysgarth check valve,
allows this valve to reopen with a further period of flow through the
valve from upstream occurring. This sequence is repeated on one
further occasion.
Figure 12.21 shows the predicted variation of air volume within the
Aysgarth vessel. After trip at time ¼ 10 s, air volume starts to expand
but after a short interval volume is stabilised by the redevelopment of
flow through the pump. Wave reflection in the suction main causes
upstream head to fall and the check valve to close after about 45 s.
The Aysgarth main now becomes entirely reliant upon the vessel for
a continued supply of water. Air volume in the pressure vessel then
195
Pressure transients in water engineering
0.40
0.35
0.30
Volume (m3)
0.25
0.20
0.15
0.10
0.05 Volume
0
0.270
9.180
18.090
27.000
35.910
44.820
53.730
62.640
71.550
80.460
89.370
98.280
107.190
116.099
125.009
133.919
142.830
151.740
160.650
169.560
178.470
187.380
196.290
205.201
214.111
223.021
231.931
240.841
249.751
258.661
267.571
Time (s)
Fig. 12.21. Air volume variations after pump trip at booster station
expands to its maximum. Subsequent reopening of the check valve at
around 80 s produces a reduction in vessel volume as head rises and
flow passes through the Aysgarth pump and check valve. Thereafter
only modest variations in vessel volume were predicted to occur.
The development of pump turbining action has been beneficial in
two respects.
(a) When head upstream of the pump exceeds the downstream head
and the pump commences to turbine, this acts as a form of ‘relief
valve’ allowing water to escape from the pressurised suction main
and preventing further significant head rise on the upstream side
of the pumping station.
(b) The continued flow through the pumps while turbining means that
the downstream pipeline is not entirely dependent upon the pres-
sure vessel for a supply of water to achieve hydraulic transient alle-
viation but some of the water is being provided by flow continuing
to pass through the pump during the intervals of turbining.The
capacity of pressure vessel can be made smaller than if no turbining
action occurred.
196
Pressure vessels
350
300
Head (mAOD)
250
200
150
100
u/s head
50
d/s head
0
59.890
64.395
68.900
73.405
77.910
82.415
86.920
91.425
95.930
100.440
104.940
109.450
113.950
118.460
122.960
127.470
131.970
136.480
140.980
145.490
149.990
154.500
159.000
163.510
168.010
172.520
177.020
181.530
186.030
190.540
Time (s)
Fig. 12.22. Head variations at booster station after all pumps are tripped
197
Pressure transients in water engineering
400
350
300
Head (mAOD)
250
200
150
100
u/s head
50 d/s head
0
0.265
2.385
4.505
6.625
8.745
10.865
12.985
15.105
17.225
19.345
21.465
23.585
25.705
27.825
29.945
32.065
34.185
36.305
38.425
40.545
42.665
44.785
46.905
49.025
51.145
53.265
55.385
57.505
59.625
61.745
63.865
65.985
Time (s)
Fig. 12.23. Head variations at downstream station while upstream pumps are
running
198
Pressure vessels
199
Pressure transients in water engineering
water to flow into the suction vessel rather than continue to the pumps
has the effect of further aggravating rates of flow deceleration through a
failing pump thus imposing more severe closure conditions on the pump
discharge check valves.
As the transient progresses it is not uncommon for the head upstream
of the pumps to exceed the downstream head. Then the check valves
will reopen, allowing flow to re-establish and also largely stabilising air
volumes within the vessels. Following flow reversal upstream and down-
stream of the pumping station, the upstream head will start to fall while
downstream head rises. Check valves will then reclose.
In the case of high-lift pumps, the check valves may not reopen if the
head upstream of the pumps does not rise to the minimum level of fall-
ing head in the downstream vessel. In these circumstances the more
rapid closure of the check valves will mean that more water will
enter the upstream vessel and more water will be withdrawn from the
downstream vessel. There will be a tendency for the necessary vessel
volumes to increase above the volumes required where a vessel was
installed only on one side of a pumping station.
As with a pressure vessel located on the downstream side of a pump-
ing station, throttling of flow can also be applied to the upstream vessel.
This throttling action may be applied either during inflow to the vessel
or to the outflow depending upon a need to suppress the upsurge on the
suction side (throtting of inflow) or the downsurge (throtting of out-
flow).
A pressure vessel may also be used as part of an overall protection
package. It can be used to prolong the time taken for piezometric
level to fall, at the start of a main, to a sufficient extent that no
unacceptably low pressures are developed between downstream points
on the main where other forms of protection are present.
200
Pressure vessels
201
Pressure transients in water engineering
habs
M
z h
Minimum hydraulic gradient:
no vacuum allowed
vacuum allowed
Gas volume = Vol
Pumping station
202
Pressure vessels
203
Pressure transients in water engineering
or
Vp ðe em Þ þ zðVol Volm Þ constant lnðVol=Volm Þ ¼ 0
for n ¼ 1 ðA12:3aÞ
and
Vp ðe em Þ þ zðVol Volm Þ
204
Pressure vessels
and
Vp eo þ Volm zðVolo =Volm 1Þ
Volm hm =ð1 nÞðVolo ho =Volm hm 1Þ ¼ 0 for n 6¼ 1
The equations can also be arranged to allow solution for Volm , noting
that for the frictionless case ho ¼ z, thus,
Volm ¼ Vp eo =fzð1 hm =zÞ þ hm lnðhm =zÞg for n ¼ 1 ðA12:5aÞ
and
Volm ¼ Vp eo =fzð1 ðhm =zÞð1=nÞ Þ þ hm =ð1 nÞððhm =zÞð1=n 1Þ 1Þg
for n 6¼ 1 ðA12:5bÞ
205
Pressure transients in water engineering
ðA12:8aÞ
For n 6¼ 1 the resistance term becomes:
ð ð
F e dVol ¼ F=Vp z ðVolm VolÞ dVol
ð
ð1 nÞ
þ constant=ð1 nÞ ½Vol Volð1
m
nÞ
dVol
Volð1
m
nÞ
Vol þ Volð2
m
nÞ
ð1 nÞ=ð2 nÞg
substituting VolðnÞ
m hm ¼ constant then,
ð
F e dVol ¼ F=Vp fVol2m ðz=2Þð1 Vol=Volm Þ2
206
Pressure vessels
207
Pressure transients in water engineering
or
b ¼ zðð1=nÞ 1Þ hm =ð1 nÞðð1=n 1Þ 1Þ n 6¼ 1
where c ¼ Vp eo .
The value of ¼ hm =ho is obtained with ho being the initial steady
flow absolute pressure head in the vessel and z is the level difference
between the pressure vessel and the downstream reservoir þ atmos-
pheric pressure head. Note that ho ¼ z þ Feo .
Example
Z ¼ 100 m, L ¼ 1000 m, D ¼ 1:0 m, eo ¼ 0:2, F ¼ 20, ho ¼ 104 m and
hm ¼ 20 m. Vp ¼ 785 m3 and ¼ hm =ho ¼ 20=104 ¼ 0:19231
n¼1
a ¼ ð32:618343 þ 9:8128514Þ 20=785 ¼ 0:5810316
b ¼ 80:769231 þ 32:973173 ¼ 47:796058
c ¼ 785 0:2 ¼ 157
p
Volm ¼ f47:796058 ½2284:4632 4 ð0:581036Þ 157g
=½2 ð0:5810316Þ
¼ 3:163 m3
The corresponding result with zero friction is:
Volm ¼ 785 0:2=ð80:0 þ 32:188758Þ ¼ 3:284 m3
n ¼ 1:001
a ¼ ð32:592745 þ 9:807059Þ 20=785 ¼ 0:580527
b ¼ 80:737531 þ 32:968 ¼ 47:769531
c ¼ 157
p
Volm ¼ f47:769531 ½2281:9281 4 ð0:580527Þ 157g
=ð1:161054Þ
¼ 3:165 m3
n ¼ 1:2
a ¼ ð27:891309 þ 8:661267Þ 20=785 ¼ 0:4898135
208
Pressure vessels
209
Pressure transients in water engineering
Project Max. air volume, Max. air volume, Max. air volume, estimated,
detailed analysis estimated, n ¼ 0 n ¼ 0, zero friction
(m3 ) (m3 ) (m3 )
become:
for n ¼ 1:
zð 1Þ hm lnðÞ þ ðFÞVolm =Vp ðz=2ð1 Þ2
þ hm ½1 þ ðlnðÞ 1ÞÞ ¼ 0 ¼ fnðÞ ðA12:10aÞ
and for n 6¼ 1:
zfð1=nÞ 1Þ hm =ð1 nÞðð1=n 1Þ 1Þ
210
Pressure vessels
65 fn(x)
60
55
50
45
40
35
30 x = hm/hp
25
20
15
hmax = hm exp{z/hm}
10 x2 = hm/hmax
–5
211
Pressure transients in water engineering
Example
The following values were used to illustrate the iterative solution which
can be easily programmed.
n ¼ 1, z ¼ 100 m, hm ¼ 20 m, L ¼ 1000 m, D ¼ 1 m, hmax ¼
hm expðz=hm Þ ¼ 20 expð100=20Þ ¼ 2968:3 m, hmin ¼ z ¼ 100 m,
FVolm =Vp ¼ 20:0 3=785 ¼ 0:0764
Since it is only an estimate of peak head that is required, the process
can be truncated when a sufficient degree of convergence has been
achieved.
For the Kielder project, maximum upsurge pressure in the vessels was
estimated at 373.9 m absolute. The detailed analysis gave a maximum of
335 m absolute.
212
Pressure vessels
Example
The desired value of maximum pressure hp in the vessel is selected.
Volm ¼ 3 m3 , Vp ¼ 785 m3 , hm ¼ 20 m, hp ¼ 500 m, z ¼ 100 m.
n¼1
F ¼ f96:0 þ 64:377516g=f46:08 þ 16:624899g 785=3
¼ 280:92
n ¼ 1:001
F ¼ f95:987117 þ 64:416g=f46:067633 þ 16:629134g 785=3
¼ 280:62263
n ¼ 1:2
F ¼ f93:160096 þ 70:9976g=f43:394018 þ 17:219815g 785=3
¼ 221:56
n ¼ 1:41
F ¼ f89:801103 þ 75:596298g=f40:321191 þ 17:373275g 785=3
¼ 161:9722
The throttle at the vessels’ inlet in order to reduce the peak pressure
to 252 m absolute was estimated at 230 for n ¼ 1 while the detailed
analysis gave a throttle value of 125. Use of n > 1 gave higher values
of throttle coefficient, for example n ¼ 1:2 produced an estimated
throttle coefficient of 260.
213
Pressure transients in water engineering
214
13
Further aspects of pressure
vessels
215
Pressure transients in water engineering
Bourdon-type
pressure gauge
Silicon fluid filling point
High-level alarm
Manway
Compressor start
Compressor stop
Low-level alarm
Pressure relief valve Water level indicator
Inflow/outflow
connection
Drain valve
are used three-stage oil removal filters can be incorporated. These filters
will reduce oil carryover to <0.003 ppm where the temperature of air
flowing into the filter has been cooled to <308C after leaving the com-
pressor. Compressor capacity should be a balance between not allowing
the set to run too long and overheating, and too short not allowing the
set to reach a correct working temperature. Unmanned installations
should incorporate low oil level protection. Critical installations may
use duty and standby compressor sets.
While not absolutely necessary, air receivers are often installed in
surge vessel systems. Instead of the compressor delivering direct to
the vessel, a solenoid valve is included downstream of the receiver
and this opens to allow a flow of air into the vessel.
216
Further aspects of pressure vessels
217
Pressure transients in water engineering
218
Further aspects of pressure vessels
should consider site conditions carefully and all failure conditions which
might develop.
A circulation system can be included in those larger installations
containing a substantial volume of water in the vessel. Eventually the
water can become stale unless an arrangement to change the water is
incorporated. A small pump and pipework system providing about
5 litres/min (7 m3 /day) will ensure that water in the vessel remains
sweet.
219
Pressure transients in water engineering
Air or nitrogen
Liquid-filled bladder
Outlet to pipeline
(a)
Anti-extrusion ring
Shell in stainless steel,
polypropylene and PVC Unrestricted
Liquid inlet/outlet inlet–outlet port
(b) (c)
between the bladder and the vessel walls. As pressure drops and the gas
expands, then to prevent possible expansion of the bladder out of the
vessel a poppet valve is included which is pushed shut by the expanding
bladder thus preventing the bladder being damaged. Its chief function is
to store fluid under pressure. Stored fluid is then released to the pipeline
circuit to supplement pump output and to absorb pressure peaks and
shocks generated elsewhere in the system.
220
Further aspects of pressure vessels
45.0 m
21.
4m
Sight Max WL
glass
Pumping main Min WL
221
Pressure transients in water engineering
Anti-vacuum valve
Triple air valve cluster
Varies
Typical arrangement of
anti-vacuum valve Letch House tank
250 Elevation = 223 mAOD
50
16.5 m
P
0
0 930 1400 1760 2160 2455 3050 3815 4230 4865 6207
Chainage (m)
pipeline. During the time delay of around 0.06 s before the vessels are
able to respond, the first component of wave has been travelling
along the rising main shown in Fig. 13.4. The resultant effect is a
short-duration pressure pulse which travels back and forth in the
rising main on top of the main surge controlled by the pressure ves-
sels. Should this short-term pressure fluctuation coincide with the
peak or trough of the primary wave, the effect will be to increase or
decrease the maximum or minimum transient pressure based upon
the main wave only. The predicted and observed head variation at
Riding Mill pumping station following simultaneous trip of six duty
pumps, can be seen in Fig. 13.5. Irregularities in head variation are
present on both traces.
Within the Riding Mill pumping main, this effect was responsible for
operation of an air valve at Penn’s Hill (Fig. 13.6), although a
preliminary analysis which neglected the effect of vessel connection
length had suggested that no such valve operation would occur. Both
the observed and simulated head variations at Penn’s Hill show clearly
the presence of higher-frequency components of wave motion. These
waves travel together with the primary wave which is dominated by
the presence of pressure vessels. Where possible, the pressure vessels
should be placed as close as possible to the pumping main to avoid
unnecessary delay in vessel response and to minimise the development
of these secondary wave components.
222
Further aspects of pressure vessels
320
6 pump trip record at Riding Mill PS
300
d
ve
280
er
ed
bs
260 ict
O
240 ed
Pr
220
Elevation (mAOD)
200
180
160
140
120
100
80
60
40 n = 1.4
20 a = 900 m/s
0 K = 85
?60
Recording failure 6 pump trip record at Penn’s Hill
t ed
e dic
?40 Pr
?20
Elevation (mAOD)
Observe
?00
d
223
Pressure transients in water engineering
the pipeline to the centreline of the vessel may be quite significant. The
smaller the difference in level from the vessel centreline to the down-
stream discharge level the larger will be the required vessel volume.
In systems of low static it is important to employ a realistic assessment
of vessel elevations if a possible underestimation of vessel volume is to
be avoided.
It has been said that the function of the pressure vessel is to limit
changing head or pressure along the rising main. This leaves the rela-
tively short stretch of pipeline between the suction well upstream of
the pump and the connection with the pressure vessel to be considered.
The function of the pressure vessel is to sustain pressure at the vessel
connection and along the rising main. However, on the downstream
side of a failing pump, pressure can fall relatively quickly and a steep
adverse hydraulic gradient can be established within this relatively
short section of pipeline between the pump and the vessel connection
where pressure is being maintained at a relatively high level (Fig. 12.4).
Compared with the situation without a pressure vessel, jdV=dtj
downstream of the vessel is reduced while jdV=dtj upstream of the
vessel is increased.
Flow will tend to reverse more quickly upstream of the vessel connec-
tion than if the vessel were not included. This has important implica-
tions with respect to closure of the check valve downstream of the
pump and will be considered in greater detail when discussing check
valves in Chapters 20 and 21. It may be noted that siting of the
vessel close to the pumps will tend to aggravate check valve closure
conditions and if practical the vessel should be moved further down-
stream rather then being say teed off a manifold directly opposite the
pump branch connections (Fig. 13.7).
Suction well
Pumps
To rising main
Check valves
Discharge manifold
Pressure vessel
224
Further aspects of pressure vessels
The pressure vessel does not suppress the transient at its source as
does the flywheel, but rather the transient is developed at the pump
and the vessel acts to suppress its effects to an acceptable extent as
the pressure wave enters a downstream pipeline for instance. As
described, the positioning of the vessel and configuration of its con-
nection can influence the effectiveness of response.
(a) At static, the air valve is closed and the tube is filled with liquid. A
charge of air is enclosed in the upper part of the vessel above the
lower end of the tube and under the static pressure head.
(b) While pumping, the air valve remains closed with the tube liquid-
filled as before. The air charge is compressed under the pumping
head in the top of the vessel.
(c) When head starts to fall after a pumping failure, the air mass starts
to expand, releasing liquid into the rising main thus suppressing the
transient head drop in the main. Thus this vessel acts in the same
manner as a conventional pressure vessel at this stage. The dip-tube
still remains liquid-filled at this time.
(d) After a time when head has fallen to the air valve operating level,
the valve will open to allow air to flow into the tube and finally into
the vessel itself. This air flow augments the expanded air volume
already within the vessel. Operation of the air valve largely stabi-
lises head at the vessel and the continued flow of liquid out of
the chamber is essentially balanced by air entering the tank via
the air valve.
(e) After flow reversal has occurred, head starts to rise at the vessel
with air flowing out of the air valve and liquid entering the
vessel. When liquid level reaches the lower end of the dip-tube
the air mass contained in the vessel outside the tube is isolated
while the air in the tube exits via the air valve. The air valve
225
Pressure transients in water engineering
Protective cap
Pressure gauge
Flanged outlet
Interchangeable pipeline
connection DN 100 to DN 400
To
pipeline
shuts when liquid level reaches the top of the tube and the air mass
remaining in the vessel is compressed by the return flow in the
vessel.
226
Further aspects of pressure vessels
Static head
(a) (b)
Head shortly after trip
AV closed
(c)
Head rising during reversed flow
227
Pressure transients in water engineering
General layout
To rising mains
228
Further aspects of pressure vessels
In the event of pump shutdown, pressure falls and the air mass in the
lower chamber expands allowing liquid to flow from the vessel thus
providing a measure of surge protection. Only when hydraulic gradient
has fallen to the operating level of the air valves will these valves open
to allow air inflow which augments surge protection.
229
14
Surge tanks and related
structures
230
Surge tanks and related structures
Cap
Ground level
M
2134 j tunnel
1524 j steel pumping main
Pumping level
M
M
M
Pumping flow
Tunnel water level maintained
by throttle valve
Reversed flow in main while downstream = 67.67 mAOD
pressure vessels and
Static level in shaft = 68.58 mAOD
feeder tanks refilling
Fig. 14.1. Simple surge shaft between pumping main and gravity flow tunnel
231
Pressure transients in water engineering
305 M
Surge pond
295
M
Steel-lined penstock to
underground generating station
232
Surge tanks and related structures
Piezometric level in
pipeline at connection
Piezometric level
at water surface M
Ht H
Vt At
+
Ac
Au/s Ad/s
Vu/s Vd/s
+ +
that:
Vu=s Au=s At dHt =dt ¼ Vd=s Ad=s ð14:1Þ
Head difference between the connection and the water surface
Ht H is made up of two components: inertial head and head loss in
the connection, so that:
Ht H ¼ fKL þ 1gfAt =Ac dHt =dtg2 =ð2gÞ þ Lt =g dVt =dt ð14:2Þ
where KL is the overall loss in the connection, At is cross-section of the
chamber, Ac is cross-sectional area of flow in the connection and
Vt ¼ dHt =dt.
In addition the Cþ and C characteristics provide the following rela-
tionships at the time when the wave paths arrive at the tank connection:
Vu=s þ g=au=s H ¼ Jþ and Vd=s g=ad=s H ¼ J
Four equations are available to solve for the unknowns, Vu=s , Vd=s , H
and Ht . The above equations can be solved in a straightforward manner
with dHt =dt ðHt Hto Þ=t.
233
Pressure transients in water engineering
Since cross-sectional areas are the same and wave speed is also
assumed common, this reduces to:
X
Hj ¼ a=g ðJÞ=3 ð14:3Þ
At the free surface (Fig. 14.4), assuming the positive direction of flow to
be from the shaft into the pipeline, then for conservation of volume:
VAp ¼ As dHs =dt ð14:4Þ
or averaging over a time increment t:
ðV þ Vo ÞAp =2 ¼ As ðHs Hso Þ=t
giving:
V þ Vo ¼ 2As =ðAp tÞðHs Hso Þ
Also from the characteristic arriving at the free surface:
V g=aHs ¼ J
234
Surge tanks and related structures
Piezometric level at
the connection point
Piezometric level
at the free surface
Ap
Hj Ht
Dx
Ap
Dx +
Ap
Dx
This neglects the inertia of that part of the water column above the first
computational section. Substituting for V in the conservation equation
at the free surface:
Hs ¼ ð2As =ðAp tÞHso J Vo Þ=ð2As =ðAp tÞ þ g=aÞ ð14:5Þ
Having found Hs , the velocity at the free surface can be obtained from:
V ¼ J þ g=aHs
In this way the response of the shaft is included and the effect of a
pressure transient partially bypassing the shaft can be modelled. Varying
water level in the shaft requires some method of identifying those
sections of shaft which are above water level and only the parts of
235
Pressure transients in water engineering
the system below the water level included in the analysis. Where raw
or treated water is being handled there is also the possibility of allow-
ing some spillage from the top of the surge pipeline provided the loss
of water is acceptable and the spillage can be conveyed safely to a
convenient watercourse.
Dawson (c. 1980) described use of pipes as a form of surge tank in
sewage pumping applications. Clearly any spillage from the top of the
surge pipe has to be avoided and there is some risk of offensive
odours being vented when sewage level in the shaft is rising. Sizing of
the surge shaft requires consideration of possible septic conditions
developing if the sewage is allowed to dwell within the shaft for too
great a time. During downsurge the shaft must essentially empty.
Thus the shaft cannot have too great a capacity. On the other hand,
if too small a shaft volume is provided it may not fulfil its surge suppres-
sion function and violent oscillations may develop which could
adversely affect other parts of the system such as check valves.
236
Surge tanks and related structures
dx/dt = V – a
As Valve
Ls
A A
V = Vo V = 0.0
have wave propagation times much longer than in the surge tank itself
— that is, Lu=s =a Ls =as and Ld=s =a Ls =as . If a valve at the down-
stream end of the pipeline is abruptly closed with an initial steady
flow velocity Vo the inertial head rise Hi will be aVo =g as shown. The
pressure wave thus created will travel upstream to the surge chamber
connection at a rate dx=dt ¼ V a. At the connection the head rise
which occurs as the wave front arrives is given by:
Hj ¼ fVo A ðVo ÞA þ 0As g=f2A þ As g
For simplicity, if a common wavespeed a is assumed for both pipeline
and surge shaft then:
Hj ¼ a=gVo 2A=ð2A þ As Þ ð14:6Þ
This head rise will induce a velocity Vs into the surge chamber given by:
Jþs ¼ 0 ¼ Vs þ g=aHj
or
Vs ¼ 0 g=aHj ¼ g=aa=gVo 2A=ð2A þ As Þ
¼ Vo f2A=ð2A þ As Þg ð14:7Þ
and an invariant value
Js ¼ Vs g=aHj ¼ 2Vo f2A=ð2A þ As Þg ð14:8Þ
237
Pressure transients in water engineering
Table 14.1. Ratio of head rise after reflection to initial head rise at connection
Time from arrival Ratio of inertial head rise Hjn =Hj1 at junction
of pressure wave
0 1
2Ls =a 1 A
4Ls =a 1 2A þ A2
6Ls =a 1 3A þ 3A2 A3
8Ls =a 1 4A þ 6A2 4A3 þ A4
10Ls =a 1 5A þ 10A2 10A3 þ 5A4 A5
12Ls =a 1 6A þ 15A2 20A3 þ 15A4 6A5 þ A6
238
Surge tanks and related structures
Surge shaft
239
Pressure transients in water engineering
240
Surge tanks and related structures
6 5
4 3
Pumping station
2 1 Surge
vessels
NRV
DN 1200 bypass
241
Pressure transients in water engineering
below the prevailing water level in reservoir No. 1 then a check valve on
the bypass line will open to permit flow from reservoir No. 1 into the
rising mains thus augmenting flow from the vessels and inhibiting
further substantial head drop in the mains. The bypass line is some
60 m in length and overall losses in this line should be taken into
account in any simulation.
975 mm
800 mm
750 mm
300 mm
Warren Park PS
375 mm
1000 mm
1500 mm
Viking SR complex
450 mm
242
Surge tanks and related structures
Viking SR, causing outflow from this storage facility. This outflow pro-
vides a source of water for surge suppression in addition to that from the
pressure vessels. The vessels’ volume could be reduced by 50 m3 were
Viking SR relied upon to yield water after all pump trip events.
Use of service reservoirs in this way depends upon the existence of a
backflow connection. Otherwise, instead of outflow from the reservoir,
air valves may open allowing the line to deprime. If there are no air
valves, high vacuum pressures may occur in the branch serving the
reservoir. Even if it is not the intention to utilise an SR for surge
duties in this way, hydraulic transients are such that the reservoir
branch will be affected and should therefore be included in any
modelling exercise. If it is not the intention to allow reversed flow
from the SR and it is necessary to avoid vacuum pressures then a
check valve can be included on the branch main close to its connection
to the trunk main.
243
244
Primary units
Secondary
units
EL +81.0
M
Pressure transients in water engineering
GL +78.80
M GL +78.70
CIP M
f600
f600
f700
EL +75.50
M
Contact tank
WL = 78.5 mAOD
DN 450 DN 600
DN 500
DN 600
UF membrane plant
Permeate line to
contact tanks
245
Pressure transients in water engineering
Steady pumping
Piezometric line –
min. resistance Steady flow – max. resistance
min. resistance
Primary strainers
Surge tank Min. piezometric level
Membrane with surge tank
units
Feed tank Contact tank
WL 78.5 WL 78.5
Sub-atmospheric head
246
Crane 30 T
+13.50
2000
247
Pressure transients in water engineering
0.048
6.166
12.285
18.403
24.521
30.640
36.758
42.876
48.994
55.113
61.231
67.349
73.467
79.585
85.703
91.822
97.940
104.058
110.176
116.294
122.412
128.531
134.650
140.769
146.888
153.007
159.126
165.245
171.365
177.484
183.603
189.722
Level (mASL) –0.5
–1.0
–1.5
–2.0
–2.5
Time (s)
248
Surge tanks and related structures
Qw
Seal weir
H
Inlet pipeline
Zw
Cross-section = A
Tailwater
249
Pressure transients in water engineering
5.4
5.3
5.2
Level (mASL)
5.1
5.0
4.9
4.8
0.048
6.166
12.285
18.403
24.521
30.640
36.758
42.876
48.994
55.113
61.231
67.349
73.467
79.585
85.703
91.822
97.940
104.058
110.176
116.294
122.412
128.531
134.650
140.769
146.888
153.007
159.126
165.245
171.365
177.484
183.603
189.722
Time (s)
Fig. 14.15. Head variation upstream of seal weir following pumping failure
M 153.700 mASL GL
250
Surge tanks and related structures
Al Raiqa ELT
Mizaira’a
ND 400 DI
ND 150
Mizaira’a PS
Chlorination
Surge protection
To Nafeer
ND 400
ND 150
ND 150 DI
Sabkha ELT
ND 400
ND 400 DI
ND 150 DI
Gurmeda ELT
ND 150 DI
ND 400 DI
Shah ELT
Tharwaniya GT ND 150 DI
ND 400 DI
ND 400 DI
ND 150 DI
To Hameem
Hawatheen ELT
Fig. 14.17. Pipeline system serving water towers and ground tanks
251
Pressure transients in water engineering
210
190
170
Elevation (mOD)
150
130
110
4 8 12 16 20 24 28
Chainage (km)
Fig. 14.18. Envelope curves for pump trip with water towers connected
210
190
170
Elevation (mOD)
150
130
110
4 8 12 16 20 24 28
Chainage (km)
Fig. 14.19. Envelope curves for pump trip without water towers
252
Surge tanks and related structures
Pumping station
Pipe bridge
Static hydraulic
gradient
Air valve
Maximum flow
Piezometric level
hydraulic gradient
after pump trip
Maximum flow
Hydraulic gradient at
HWL,
Twin seaward pipelines
LWL
LWL HWL
Diffuser section
253
Pressure transients in water engineering
19.xxx m
2971
Detail P
i.l. 13.450 m Drg No. P2325/5/R16
Detail D 2 No. 152 × 152 mm × 23 kg
Drg. No. P2325/5/R16 1 m UCs
2971
4 No. 365 × 368 mm × 177 kg
Galvanised 76.2 × 3.25 mm dia 1 m UCs
Weldmesh (Ref. no. 310) supported (4 No. 13.362 mm long)
on angles to within 1829 mm of Note Column ends fabricated
base: thereafter mesh is flush for bearing
25.4 × 25.5 × 2.64 mm dia.
(ref. no. 1118) 20 No. 150 × 90 mm angles
Detail F
8.476 m Drg. No. P2325/5/R15
2971
Detail G
See detail L
Drg No/ P2325/5/R15 Drg No. P2325/5/R16
Inner tie 25 dia.
305×127 mm×37 kg m
UB (2 No. 5632 mm long)
16 No. 60 × 60 mm angles
(angles back to back)
2432
Detail No. 1
4.355 m
534 mm dia. flange plain ended
GRP pipes 1000 mm long
I.L. 3.720 m
2000
2445
1829
walls lined with G.R.P. 2.300 m 3.200 m
200
200
Ion
I.L. 1.450
1.100 m
Blass infill
200
200
254
Surge tanks and related structures
the seawall are primed with the air valves shut. Hydraulic gradient over
the air valves is sufficient to keep the valves tight shut. At lesser flow
rates the air valves will be open and the top of the loop represents
the downstream limit of the pumping mains. At low flows and/or low
tide level, the seaward section of the system flows under gravity. This
arrangement has the merits of providing an almost constant static
head against which pumps can operate. As far as pressure transient
behaviour is concerned, the loop of pipework creates a useful adverse
hydraulic gradient along the landward stretch of pipeline, which
produces a more rapid flow deceleration along this part of the system.
Sub-atmospheric pressures are thus minimised and quantities of air
admitted through operation of air valves along the landward stretches
of pipeline are reduced. Demands on pressure vessels located at the
pumping station are also more modest, allowing smaller vessels to be
installed. Without the inverted U of pipework at the seawall, large sec-
tions of the landward pipelines would deprime after a pumping failure,
more substantial vacuum pressures would occur and the size of pressure
vessels would have to be greatly increased. The resulting process of
restarting pumps would involve removal of large volumes of air from
the system before a steady pumping regime could be established.
255
Pressure transients in water engineering
Piezometric level
Air inflow
Air valve
(closed) Zu/s Zd/s
–
M –
M
Hu/s = Hd/s
Pressure
vessel Hu/s
Hd/s
Vu/s Vu/s
Vd/s Vd/s
(a) (b)
Zu/s
M
Zd/s
M
Hu/s Zu/s
M
Hd/s
Zd/s
Hu/s M
M
–
M
Hu/s Hd/s
Hu/s
M
Vu/s Hd/s Vu/s
Vd/s Vd/s
(e) (f)
256
Surge tanks and related structures
257
Pressure transients in water engineering
(c) Air inflow continues and effluent flows over the top of the loop
cascading into the downstream vertical pipe with level falling in
this part of the loop. Air charge within the pressure vessel continues
to expand. Hu=s < Hd=s , Zd=s < Zu=s , Vu=s < Vd=s , Qair is þve and h
is ve.
(d) Effluent level in the downstream vertical section falls below the
vessel connection and the air masses are united. Upstream flow
will finally reverse as effluent moves upstream attempting to refill
the vessel at the pumping station and removing air admitted to
the pipeline through upstream air valves. The effluent level falls
in the vertical upstream leg of the loop. Eventually upstream and
downstream effluent levels are stabilised with the seaward level
equal to the prevailing tide level. Pressure of the air mass becomes
atmospheric.
(e) When pumps are restarted effluent level in the upstream leg rises
and pressure in the air mass his þ ve. The increase in pressure
starts to push air out of the air valves on top of the loop and
volume of the air mass starts to decrease. Increased pressure
causes a small positive velocity to develop in the downstream
pipeline.
(f ) With continued pumping the upstream leg is filled and flow
cascades into the downstream pipe producing a level increase and
isolating the vessel. Under the action of rising downsteam level,
flow develops in the seaward section of pipelines. If flow is below
the maximum or if tide level is low, the loop may remain only
partially primed. When an equilibrium condition is achieved and
Vu=s ¼ Vd=s , air venting ceases with the remaining air mass being
at atmospheric pressure h ¼ 0:0.
258
15
Feeder tanks or volumetric
tanks
259
Pressure transients in water engineering
Service
reservoir
Minimum hydraulic gradient following M
a pumping failure and with pressure
vessel at PS for protection
Pipeline profile
Possible location
of feeder tank
Pressure vessel
Pumping station
Float-operated valve
TWL
Feeder tank
Filling
connection BWL
Upstream pipeline
Isolating valve Downstream pipeline
260
Feeder tanks or volumetric tanks
the water level in this tank. Feeder tanks will usually be located at a
high point on the main where preliminary analysis without local protec-
tion has shown transient sub-atmospheric pressures to occur.
To limit maximum water level in the feeder tank a non-return valve
(NRV) is positioned in the outlet connection. This NRV shuts as flow
attempts to re-enter the tank, forcing flow through the relatively
modest filling connection. As water level in the tank rises, an inlet
valve, often a float valve, progressively closes, shutting off inflow as
level approaches TWL. The filling connection is sized so as to restrict
inflow rates to ensure that no secondary surging of any significance is
generated during filling. Capacity of the feeder tank should be arranged
so that no air is able to enter the pipeline.
261
Pressure transients in water engineering
Hu/s
Hd/s
(a) (b)
M
M M
M
(c) (d)
Hu/s Hd/s
Hd/s
Hu/s
(e) (f)
262
Feeder tanks or volumetric tanks
where
P Ac is the cross-sectional area of the outflow connection and
KL is the sum of loss coefficients for the components making up
this connection — that is, entry loss, NRV loss, tee, isolating valve
and pipe resistance. The feeder tank supplies treated water from the
tank to the main, thus minimising risk of contamination.
Characteristics are selected which arrive at the feeder tank
connection at the time when a solution is required. The quasi-
invariant values propagating along these paths yield a pair of
equations at the feeder connection so that:
Vu=s þ g=aHmain ¼ Jþ and Vd=s g=aHmain ¼ J
These equations can be solved for Vu=s , Vd=s , Hmain and Htank for
normal operation.
(b) After the feeder tank and its connection have drained, initially a
common water surface and piezometric level will exist upstream
and downstream of the feeder connection. Conservation of
volume is given by:
Vu=s A þ Qair ¼ Vd=s A þ As dz=dt
where As is the air—water interface area in the pipe and
z ¼ zu=s ¼ zd=s is the water surface level.
(c) If the downstream water level should fall then a situation may
develop where upstream and downstream water levels become
different so that:
Vu=s A ¼ Asu=s dzu=s =dt þ q
and
Vd=s A ¼ q Asd=s dzd=s =dt
where q is the spillage from upstream to downstream.
(d) If a pressure vessel at an upstream pumping station is refilling, velo-
city upstream of the feeder connection may reverse, causing the
water level to decline on this side of the feeder connection, then:
q ¼ 0:0
(e) If downstream flow should reverse, spillage may occur from down-
stream to upstream of the feeder connection, so that:
Vu=s A ¼ Asu=s dzu=s =dt þ q
and
Vd=s A ¼ q Asd=s dzd=s =dt
263
Pressure transients in water engineering
(f ) During refilling water enters the feeder tank via the modest-
diameter filling connection. Inflow is regulated by a valve which
closes in response to rising water level in the tank. As level
increases, the valve progressively closes so that a quiet shut-off of
inflow is achieved as the tank level approaches top water level
(TWL).
190
Downstream reservoir
Feeder tank location
M M
180
170
Elevation (mAOD)
160
150
Pipeline profile
140
130
1070 mm diameter prestressed concrete pipeline
Stanton steel cylinder type – 15 796 m long
M Pressure vessels
Pumping station
0 2 4 6 8 10 12 14 16 18
Chainage (km)
264
Feeder tanks or volumetric tanks
265
Pressure transients in water engineering
Feeder tank
refilling
Feeder tank
full
1 min.
10 m
Start of ‘WEIR 16’
and Hmain must remain above the minimum acceptable þve pressure
otherwise a vacuum pressure may appear in the main. Note that
Htank may be at a fairly low level at this stage. In an installation
involving large pressure vessels, the demand for refilling water can be
266
Feeder tanks or volumetric tanks
267
Pressure transients in water engineering
160
Trip of 1st pump
Predicted head variation
150
Trip of 2nd pump
140
Trip of 3rd pump
Elevation (mAOD)
130
120
110
Observed head variation
100
90
2 4 6 8 10 12 14
Time (min.)
rising main had a route which was quite favourable, with pressure
vessels at the pumping station alone able to avoid development of
sub-atmospheric pressures. The new main had a less attractive profile
and using only pressure vessels would require a substantially increased
capacity than the straightforward duplication of vessel volume envisaged
for the new main.
Before embarking on prediction of transient behaviour within the
duplicated system, the opportunity was taken to compare transient
predictions for the original main with field observations. This is very
desirable for all systems but often only possible on larger installations.
Figure 15.6 shows comparisons of transient piezometric level at the
pumping station for a normal auto-sequence shutdown of three duty
pumps.
Simply adding two new vessels to the existing pair was not sufficient
to avoid sub-atmospheric pressures at higher points on the main follow-
ing a pumping failure. Figure 15.7 depicts the maximum and minimum
transient head along the new main together with the initial profile of
the pipeline.
From the initial studies, possible sites for feeder tanks were identified at
high points towards the downstream end of the new main. A number of
feeder arrangements were considered with the option of using feeders at
two different sites being investigated. Since the route had not been fina-
lised, the opportunity was taken to explore alternative possible profiles.
Figure 15.8 shows one such route with two feeder tanks at chainages
14.5 km and 18.2 km. The feeder capacities for this configuration were
268
Feeder tanks or volumetric tanks
160
140
Maximum transient
piezometric level
120
Elevation (mAOD)
100
80 Minimum transient
piezometric level
60
Pipeline profile
40
2 4 6 8 10 12 14 16 18
Chainage (km)
Fig. 15.7. Envelope curve following pump failure without feeder tank
160
140
Maximum transient
piezometric level Feeder tank location
120
Elevation (mAOD)
100
80 Pipeline
profile
60 Minimum transient
piezometric level Feeder tank
location
40
2 4 6 8 10 12 14 16 18
Chainage (km)
Fig. 15.8. Envelope curve following pump failure with possible feeder tanks
269
Pressure transients in water engineering
160
piezometric level
100
80 Modified
pipeline
profile
60 Feeder tank
location
40
2 4 6 8 10 12 14 16 18
Chainage (km)
Fig. 15.9. Envelope curves following pump failure with amended profile and
feeder tank
270
Feeder tanks or volumetric tanks
80
Maximum piezometric level
4 8 12 16
Chainage (km)
Fig. 15.10. Envelope curves following pump failure with two feeder tanks
271
272
54.600 Fall (See detail) 54.775 Fall 54.600 54.725
54.400 54.300
Item 6 Soffit 54.400
D15.
53.527
280
53.435 TWL 53.400 2
D11.
3
1
Item C
3
1118
D14. For details of access
Open
3
70 300 Wall thickening around
2621
platform see Record C.L. 51.900
300
T.O.C. H. wheel 500
25 mm grout
under baseplate 1000 50.200 m
D13.
273
Pressure transients in water engineering
53.6
53.4
Predicted level
53.2
Level (mAOD)
53.0
Observed level
52.8
52.6
52.4
52.2
52.0
0.335
31.825
63.315
94.805
126.295
157.785
189.276
220.767
252.257
283.747
315.236
346.725
378.214
409.703
441.193
472.682
504.171
535.662
567.154
598.646
630.138
661.631
693.123
724.615
756.107
787.599
819.091
850.583
882.075
913.567
945.059
Time (s)
1 min.
274
Feeder tanks or volumetric tanks
275
Pressure transients in water engineering
90
Elevation (mOD) 80
70 First upsurge
controlled by feeder Second upsurge while
the feeder is re-filling
60
50
1 2 3 4 5 6 7 8
Time (min.)
controlled by the water level in the feeder and by the inertia of the
relatively short length of main between the pumping station and the
feeder tanks. During the second and much higher upsurge, the feeder
tanks are refilling and the head at their connections is much greater.
Also the inertia of the entire length of main is involved, leading to
the much greater head rise.
276
Feeder tanks or volumetric tanks
277
Pressure transients in water engineering
278
16
Discharge conditions
279
Pressure transients in water engineering
TWL M
As
BWL or LDO
As
A
q
H
Drawoff
V
Piezometric levels
during reversed flow
280
Discharge conditions
While filling
or emptying
M
As
Ap
V Rising main
281
Pressure transients in water engineering
282
Discharge conditions
Piezometric level
while filling
Filling connection
Inlet loss
Loss in reversed
flow connection
Piezometric level
during reversed flow
Bottom water level (BWL) or
Lowest drawoff level (LDL)
M
Outflow connection
or drawoff
283
Pressure transients in water engineering
284
Discharge conditions
TWL
M
TWL
M
WL < breaker
M
285
Pressure transients in water engineering
inflow rate, piezometric level could rise to the level of water at the top of
the filling connection. During reversed flow, piezometric level will clo-
sely follow that of the tank water level while the siphon is primed. If the
siphon is broken by a low tank water level, piezometric level will rise to
the top of the filling connection and then descend within the vertical
pipe outwith the tank according to the relationship:
As dH=dt ¼ AV ð16:1Þ
If vacuum breakers are included then during inflow the siphon will
be broken at a much higher level when tank water level falls below
the level of the vacuum-breaking arrangement. While the siphon is
operating upstream, head H closely follows tank water level. Once
the siphon has been broken, H will tend to stabilise to some extent.
At low rates of inflow, head may reach the water level at the top of
the filling connection. For higher flow rates, head will also depend on
air inflow rates at the vacuum breaker. On flow reversal with tank
water level above the vacuum breaker then the siphon will operate
and upstream head will follow tank water level. If water level is below
the vacuum breaker, water level will start at the top of the inlet pipe
and descend within the vertical line outside the tank.
286
Fig. 16.6. Siphon outfall with vacuum disconnecting valve
Discharge conditions
287
Pressure transients in water engineering
288
Discharge conditions
Tailwater (canal) H
Piezometric level during
effluent level steady pumping and
with siphon primed
H Qu/s
V +ve
M
H
Q
during forward flow. As air is removed, the siphon is finally primed with
piezometric level over the summit, falling to reach tailwater
level þ downstream losses (Fig. 16.8). Sub-atmospheric pressure will
prevail over the summit, thus minimising pumping head.
For a typical pumping system of the Al Ghouta Project, Fig. 16.9
illustrates changing effluent level predicted in the upstream leg of the
discharge siphon. Initially the siphon is fully primed with pumps run-
ning. After pump failure, flow rate decreases and eventually starts to
reverse. On flow reversal over the summit, the vacuum disconnecting
valve opens and the siphon is broken. Under reversed flow, the effluent
level starts to fall, slowly at first and then more steeply as the upstream
leg of the siphon starts to empty in order to supply the liquid necessary
for pressure vessel refilling. When the effluent level reaches the rising
main, piezometric level is largely stabilised for a time by the relatively
flat pipeline profile. When flow becomes positive once more, the
upstream leg of the siphon partially refills. When vessels are filling,
289
Pressure transients in water engineering
622
System primed
621
Valve opens
620
Elevation (mOD)
Water level in
619 discharge pipework
the siphon pipes are emptying to provide the vessels with the required
liquid and when the vessels release effluent, piezometric level in the
upstream part of the siphon is recovering.
Effluent level is again determined by the equation:
As dH=dt ¼ AV ð16:1Þ
with As ¼ Au=s =sinðu=s Þ while H soffit of rising main and As ¼
A=sinðÞ for H < soffit of rising main. The much shallower gradient
of the rising main usually has a stabilising influence upon water level
and piezometric level.
The effect of this oscillation in head at the siphon end of the system
can be seen in Fig. 16.10 which depicts predicted maximum and mini-
mum hydraulic levels along the rising mains. Initial downsurge after
pumps are tripped takes place while flow in the system is positive.
Along the upstream parts of the pumping mains the minimum down-
surge is developed while discharge head is controlled by the tailwater
level during siphoning. After flow reversal, head at the siphon end of
the system first rises as the vacuum disconnecting valve opens, and
then falls as the upstream limb of the siphon deprimes. Maximum
head towards the downstream end of the system is influenced by the
maximum level as the siphon is broken. As the vessels refill, head
at the downstream end of the system decreases and so peak head at
the pumping station is partly determined by the low downstream
head resulting in a reduced maximum pressure. Minimum head along
290
Discharge conditions
630
Maximum piezometric level
Vacuum
breaker valve
625
Elevation (mOD)
620
Minimum piezometric level
615 Pumping
station
610
Pipeline profile
Maximum transient
630 hydraulic gradient
Vacuum breaker
625 valve
Elevation (mOD)
610
Pipeline profile
291
Pressure transients in water engineering
35
30
il (mASL)
25
hmax
hmin
Elevation (mASL)
20
15
10
0
301.5
603.0
904.5
1206.0
1507.5
1809.0
2110.5
2411.9
2713.4
3014.9
3316.4
3617.9
3919.4
4220.9
4422.3
4725.2
5028.0
5330.8
5633.6
5936.4
6144.9
–5
Chainage (m)
Fig. 16.12. Envelope curves for sewage rising main with air valves
292
Discharge conditions
35
d/s vessel Start 450 Mid 450
30
25
20
Head (mASL)
15
10
0
0.276
19.596
38.916
58.236
77.556
96.876
116.196
135.516
154.836
174.156
193.477
212.797
232.117
251.437
270.757
290.077
309.397
328.717
348.037
367.357
386.677
405.997
425.317
444.637
463.957
483.278
502.598
521.918
541.238
560.558
579.878
599.198
–5
Time (s)
293
Pressure transients in water engineering
294
17
Air valves
295
Pressure transients in water engineering
Cowl
Cover
Orifice cover Seat ring
Seal ring
Orifice bracket Float guide
Sealing face
Float
Fulcrum pin
Body
Float and lever
Body
Adjusting
screw
296
Air valves
HGL
DOV DOV
DOV DOV
PS
(d) SOV
SOV
SOV
(e)
SOV
SOV
SOV
(f)
800 m 800 m
DOV
M M
DOV
(g) (h)
297
Pressure transients in water engineering
298
Air valves
299
Pressure transients in water engineering
1.4
Small orifice – outflow only
1.2
0.6
Q = 1.02 ¥ ÷(Dpp) Q = 48 ¥ ÷(Dpp)
0.4
Outflow
0.2
Inflow
Large orifice – inflow and outflow
40 30 20 10
10 20 30 40 50 60 70 80 90
Discharge of free dry air in STD m3/min
0.2
Inflow
0.4
Vacuum gauge
pressure (bar)
0.6
Supercritical flow
0.8
300
Air valves
Air valve
Air inflow Qair
Vu/s Vd/s
+x
C+ characteristic C– characteristic
Dt
Dx Dx
Space << Dx
a more rapid flow reversal in that section of main upstream of the valve.
In a pumping main where air valves have been provided, it is possible
that some of these valves will operate following a pumping failure. This
is particularly likely if no special surge protection has been included at
the pumping station to prevent this from happening. After the pumping
failure, as the negative pressure wave or downsurge propagates into the
rising main, the hydraulic gradient at an air valve may fall to the point
where the valve opens to allow inflow of air.
Where a number of air valves along a pipeline operates during the
downsurge following pumping failure then air pockets will be developed
in the pipeline below the connection to each operating valve.
Successive air valves will each act to further suppress the downsurge
so that the minimum piezometric line may appear as shown in
Fig. 17.6. For each valve which functions, deceleration of flow in the
section of line downstream of that valve is < deceleration rate upstream
of the valve connection. Suppression of the pressure transient
downstream of air valves has been achieved by providing alternative
301
Pressure transients in water engineering
= air valve
Successive positions
of pressure wavefront
Pipeline
profile
Zones of transient
sub-atmospheric pressure
PS
Minimum transient piezometric levels without air valves
(ignoring gas release and possible column separation)
302
Air valves
Air pocket
Head drop through
upper air valve = KinV 2/(2g)
M
Hydraulic grade line (HGL)
using small orifice
at lower air valve
HGL using
large orifice at
lower air valve
V+
Head drop using small orifice DH
of lower air valve = KoutV 2/(2g)
303
Pressure transients in water engineering
air outflow rates are substantial and the head loss through a valve will
remain relatively modest.
While the air valves may be instrumental in alleviating the initial
downsurge, in many instances the ‘secondary’ pressure transients
produced by closure of air valves can be more severe than the initial
downsurge after a pumping failure. This problem arises because of the
substantial hydraulic gradient which can develop between neighbouring
operating air valves where inflow/outflow of air is via a large orifice.
Only a modest differential head occurs across an operating air valve
while using this large orifice and so the air valve has little influence
on the development of reversed flow in the main.
Many air valves are also fitted with a second ‘small’ orifice valve
intended for venting of small accumulations of air under pressure.
Such modest quantities of air may enter the pipeline, through a pump
suction for instance, tending to accumulate at high points on a line
leading to an eventual decrease in overall efficiency of transmission.
Supposing that air outflow is only through the small orifice. Under
the initial adverse hydraulic gradient between two air pockets, flow
starts to develop towards the lower pocket and air begins to evacuate
through the small orifice of the lower air valve. As air outflow develops,
air pressure within the lower pocket increases and the volume of air
decreases due to compression as well as venting.
The rising head at the lower air pocket will cause hydraulic gradient
to flatten between the adjacent air pockets, since:
Hupper Hlower ¼ head loss at valves
þ pipeline resistance between valves
þ inertial head
Head in the pipeline at the upper air valve will not drop significantly
below atmospheric pressure if the air pocket is increasing in volume as a
result of continued air inflow. If a still higher air valve is present further
downstream, it is possible that air may also be flowing out of the middle
air valve, in which case the head at this valve will also rise significantly.
In contrast with the earlier illustration when the large orifice was used
to discharge air at the lower valve, use of the small orifice restricts rates
of outflow and flattens the hydraulic gradient thus reducing head
available to accelerate flow and to overcome pipeline resistance.
Accordingly dV=dt is smaller and the eventual flow velocity when the
valve shuts is also limited. Upsurge after valve closure is correspondingly
lower because the differential velocity is smaller.
304
Air valves
305
Pressure transients in water engineering
or
h X i
2
1 F=HVðmaxÞ ¼ exp VolðmaxÞ =Volp fL=D þ KL
Finally,
pn X
VðmaxÞ ¼ 2gH= fL=D þ KL
h X io
1 exp VolðmaxÞ =Volp fL=D þ KL ð17:5aÞ
If VolðmaxÞ ! 1,
pn X o
VðmaxÞ ! 2gH= fL=D þ KL ð17:5bÞ
306
Air valves
0
3
2.5
1
2
1.5 3
7
1
15
0.5 31
0
0.1 0.316 1.0 3.16 10.0 31.6 100.0 316.0
Maximum air volume, Vol(max) (m3)
V(max) = ÷{2gDH/(fL/D + SKL)(1 – exp[–Vol(max)/Volp(fL/D + SKL)])}
DH = 10 m; L = 1000 m; D = 1000 mm; f = 0.02; Volp = pD 2L/4
Fig. 17.8. Maximum reversed velocity plotted against initial air volume at air
valve
307
Pressure transients in water engineering
308
Air valves
Ho H Ho
V
Air valve
+x
Common horizontal datum
(a)
V+a
Ho
Shut valve, e.g. NRV V=0
+x
Ho Ho
V
V=0
Shut valve,
e.g. NRV +x
Common horizontal datum
(c)
Fig. 17.9. Head conditions (a) and (b) just after air valve closure and (c) as
pressure waves are reflected
309
Pressure transients in water engineering
Suppose the air valve is just downstream of a shut valve, such as the
check valve at a pumping station (Fig. 17.9b). A single characteristic
will arrive at the air valve connection at the moment of closure, then:
J ¼ Vd=s g=aHo
Just after valve closure, velocity ¼ 0 and head will be H, so that:
0 g=aH ¼ J ¼ Vd=s g=aHo
giving
H ¼ Ho aVd=s =g ð17:7Þ
Considering a wavespeed of 1000 m/s and a reversed velocity
Vd=s ¼ 1:0 m/s then the head rise after valve closure will be:
H Ho ¼ 1000:0 ð1:0Þ=9:81 ¼ 102 m
This is double the amplitude of the head rise at an air valve located
along a pipeline. For a valve along the pipeline, the abrupt pressure
rise after air valve closure will propagate both upstream and downstream
from the now closed valve as compression waves (Fig. 17.9c), thus
affecting adjacent sections of main. Subsequently the pressure waves
will be reflected from closed check valves, other still-existing air pockets
and discharge points to tanks and reservoirs. The result can be a
complicated pattern of hydraulic transient behaviour in the pipeline
as a whole.
310
Air valves
Envelope curves of maximum and minimum head along the rising main
70
50
Elevation (mHSD)
40
30
20
Pipeline profile
10
–10
Minimum piezometric level
1 2 3 4
Chainage (km)
valves when the first air mass has been expelled. Air volume variation at
chainage 2.19 km is more complicated with the air valve opening and
closing repeatedly (Fig. 17.12). Transient head variation at the
pumping station (Fig. 17.13) shows the initial fall in head after trip
followed by a long interval during which air first flows into the line
and then is removed from the pipeline, culminating in an abrupt
3500
3000
2500
Air volume (l)
2000
1500
Air valve opens Air valve opens
1000
Air valve shuts
500 Air valve shuts
Air valve opens
Fig. 17.11. Air volume variation at PS air valve after pump failure
311
Pressure transients in water engineering
150
100
Air valve shuts Air valve shuts
50
head rise after around 1.5 min. Subsequent reflection of the pressure
wave produced by this sudden head rise causes the valve to reopen
and a second period of air inflow and outflow to occur. This second
air flow interval is concluded with a further steep head rise, at around
2.5 min, to be followed by a third interval of air inflow. At chainage
2.19 km, head variations are more complicated (Fig. 17.14). Each
70
60
50
Elevation (mHSD)
30 Pump trip
10
0
PS air valve shut PS air valve shut
312
Air valves
70 Effects of PS
air valve closure
60
50
Elevation (mHSD)
Pump trip
40
Effects of air valve closure
at chainage 2.19 km
30
20
10
closure of this air valve produces an abrupt head rise and compression
waves whose reflection causes a further valve reopening. Behaviour at
this location is however dominated by the pumping station air valve
closures which produce the steep rise in head after around 1.5 min
and 2.5 min.
It is the uncontrolled acceleration dV=dt in the pipeline, between the
neighbouring open air valves which is responsible for the magnitude of
the differential velocity present upstream and downstream of the air
valve just as that valve shuts. Subsequent equalisation of velocities in
the main after air valve closure produces the shock pressure rise
predicted.
In view of the beneficial effect of throttling air outflow, it may be
thought that all air valves which may operate during a transient
event should have an air outflow regulating device fitted. A short-
coming of this approach can be the time taken for a large volume of
air to be expelled. During pipeline filling operations, if air outflow is
throttled at every valve in a system, the time taken to fill the main
could become excessive.
Where a number of air valves function during the initial downsurge
after a pumping failure, the subsequent behaviour of hydraulic transi-
ents can become quite complicated. The use of air valves can be
beneficial in relieving negative pressures as a result of pump trip.
Inflow of air via the air valves reduces rates of deceleration downstream
313
Pressure transients in water engineering
314
Air valves
Sharjah PS No. 1
35
AV 450 AV 700
30
25
Air volume (m3)
20
15
10
0
0.828
7.452
14.076
20.700
27.324
33.948
40.572
47.196
53.820
60.444
67.068
73.692
80.316
86.940
93.564
100.188
106.812
113.436
120.060
126.684
133.308
139.932
146.556
153.180
159.804
166.428
173.052
179.676
186.300
192.925
–5
Time (s)
Fig. 17.15. Air removal on pump restart with large orifice in use
of the DN 450 pipeline air valve. At around 2 min the air valve on the
DN 700 main shuts, producing an abrupt head rise and further irregular
fluctuations in head.
Replacing the standard large-orifice air valves with restricted outflow
valves has the effects of prolonging air venting times and reducing the
secondary transient effects following air valve closure. Figure 17.17
Sharjah PS No. 1
40
35
30
25
Head (mASL)
20
15
10
5
strt 700 mid 700
0
0.828
7.452
14.076
20.700
27.324
33.948
40.572
47.196
53.820
60.444
67.068
73.692
80.316
86.940
93.564
100.188
106.812
113.436
120.060
126.684
133.308
139.932
146.556
153.180
159.804
166.428
173.052
179.676
186.300
192.925
Time (s)
315
Pressure transients in water engineering
Sharjah PS No. 1
35
AV 450 AV 700
30
25
Air volume (m3)
20
15
10
0
0.828
79.488
158.148
236.809
315.469
394.129
472.789
551.450
630.110
708.770
787.431
866.091
944.751
1023.411
1102.072
1180.732
1259.392
1338.052
1416.713
1495.373
1574.033
1652.694
1731.354
1810.014
1888.674
1967.335
2045.995
2124.621
2203.247
2281.872
2360.498
2439.123
–5
Time (s)
shows the manner in which air volume changes in the DN 450 and DN
700 pipelines when pumps are started. A relatively rapid initial reduc-
tion in each air volume occurs as the pumped flow compresses the air
masses. Some pressure wave effects are present, resulting in an oscilla-
tion in air volume before an almost steady rate of reduction in volume is
attained at each site. The time to completely remove the air mass in the
DN 450 main is around 7 min, while expulsion of the air pocket in
the DN 700 pipeline takes almost 37 min. The corresponding head
variations at the start and at the mid-point of the DN 450 main are
shown in Fig. 17.18. The start-up transient has dissipated before the
air valves close. At around 7 min when the valve on the DN 450
main shuts, only a modest head rise was predicted. Similar head changes
occur in the DN 700 main. After 7 min only very small effects of the DN
450 pipeline air valve closure were predicted to occur in the larger main.
After 37 min or so, closure of the air valve on the larger main also
produces a quite modest head rise.
Throttled air valve closure on the DN 700 main produces a more
modest transient effect than an identical valve on the DN 450 main.
The same valve will produce a smaller head rise on closure as pipeline
diameter increases. To illustrate this effect, for simplicity assume that
resistance is essentially that of the air valve throttle and this can be
represented as Q ¼ Kout Habs , where Q is flow rate, Habs is absolute
pressure head at the valve connection and Kout is the valve throttle
316
Air valves
Sharjah PS No. 1
35
30
25
Head (mASL)
20
15
10
5
strt of 450 mid 450
0
0.828
19.872
38.916
57.960
77.004
96.048
115.092
134.136
153.180
172.224
191.269
210.313
229.357
248.401
267.445
286.489
305.533
324.577
343.621
362.665
381.709
400.753
419.797
438.841
457.885
476.930
495.974
515.018
534.062
553.106
572.150
Time (s)
317
Pressure transients in water engineering
H
(c)
(b)
Air valves
PS
Pipeline profile
and flow rate from the pumping station to the air valve is < in case (a).
The increased pressure in the air mass produces a downstream hydraulic
gradient which is useful in accelerating the liquid column towards the
discharge point. When the air valve closes, the differential pipeline
velocity across the valve connection will be reduced considerably
from that with the large orifice. After steady flow has been achieved,
the hydraulic gradient will be as curve (c) with the corresponding
pump operating point (c).
318
Air valves
Components
Large orifice
Sealing ring
Small orifice
Small orifice lever
Air valve float
Tappet
Air valve float guide
Air valve unit
Elevator
Guide sleeve
Main cover
Float chamber
Operating float
319
Pressure transients in water engineering
the large float chamber of the Glenfield ‘EPEX’ sewage air valve. Other
suppliers such as Erhard similarly provide a spacious float chamber for
their sewage valves.
17.7.3 Inspection
An adequate programme of inspection and maintenance of these valves
forms an important part of the overall package of protection for any
pipeline system. Figure 17.21 shows a sewage air valve chamber open
for inspection during transient tests.
320
Air valves
vent pipe terminating above the ground surface so that a free inflow/
outflow from the chamber can take place and the air valve can function
as intended.
If air valve covers are not maintained and kept free from debris and
sand for instance, air flow can be inhibited and the performance of the
air valve may not be as modelled. Where a substantial number of air
valves are being used for the purpose of hydraulic transient suppression,
the possibility of malfunction will increase. Analyses conducted should
examine a range of possibilities to ensure that satisfactory transient
conditions are achieved with one or more valves not behaving as
intended.
321
Pressure transients in water engineering
Air valve
Buffer tank
Air vol.
Rising main
level. Negligible head loss is assumed through the air valve during
inflow. It is also assumed that negligible deceleration of the water
column has occurred prior to air valve operation and that the friction
factor f and the elevation difference z remain essentially constant.
Equation (2.1) can be written:
L=g dV=dt þ fLV 2 =ð2gDÞ þ z ¼ 0 ð2:1Þ
or
L=g dV=dVol dVol=dt þ fL=DV 2 =ð2gÞ þ z ¼ 0
but
dVol=dt ¼ Q ¼ VA and LA ¼ Volp
therefore
Volp =gV dV=dVol þ fL=DV 2 =ð2gÞ þ z ¼ 0
V
M
Rising main
Buffer tank Diameter = D
Length = L
Cross-section = A
Pumping station
322
Air valves
or
fL=DVol=Volp þ lnð fL=De þ zÞ ¼ constant ð17:9aÞ
Initial velocity is Vo and e ¼ eo ¼ Vo2 =ð2gÞ while Vol ¼ 0. Substituting
then,
0 þ lnð fL=Deo þ zÞ ¼ constant
When air volume has reached its maximum as flow downstream of
the air volume comes to rest then, e ¼ 0 and Vol ¼ Volmax . Then,
fL=DVolmax =Volp þ lnðzÞ ¼ lnð fL=Deo þ zÞ
rearranging,
fL=DVolmax =Volp ¼ lnfð fL=Deo þ zÞ=zg
and finally,
Volmax =Volp ¼ D=ð fLÞ lnf1 þ hf=zg ð17:9bÞ
Equation (17.9b) was presented by Livingston in 1969. Static head is
a function of volume, and as air volume increases z will increase. Where
overall static head is appreciable, any variation in z will be relatively
unimportant. On the other hand, if say a vertical buffer tank is selected
for a low head installation, the variation in z may have a significant
influence on overall static. In these circumstances it may be worthwhile
to insert extremes of z into the equation to assess the importance of
changes in z. These extremes could include values of static head
measured from the top and bottom of the buffer tank.
Equation (17.9b) can also be used to estimate air volume when no
tank is included but the air valve sits directly on top of the pipeline.
The variation in z is likely to be more modest in this situation, with
the air volume stored within the pipeline itself.
323
Pressure transients in water engineering
324
18
Air and gas
325
Pressure transients in water engineering
Start levels:
pump 3 = –16.025
pump 2 = –16.125
pump 1 = –16.225
WPL SWALLOWGLIDE
SRVS 200DS (4 off)
–17.74 mAOD
326
Air and gas
Fig. 18.2. Air valves and non-return valves at submersible pumping station
327
Pressure transients in water engineering
Single-orifice air
release valve
DN 500 cross-over pipe
Rise
Odour pipe to
wet well
+6.0 mAOD
To rising main
DN 350 line from
foul pump
H
N = 100%
Duty point
N = 75%
Loss through
air valve
Direct start
riser filling
Loss in riser
N = 50%
Z
N = 25%
Speed necessary Q
to just fill riser Duty flow
Flow when air valve closes
Maximum flow – riser empty
328
Air and gas
progressively become filled, with liquid level and static head at the
pump increasing and discharge reducing. The air column above the
water surface in the riser will become smaller, partly through compres-
sion but mainly as a result of venting through the air valve located at
the top of the riser. Preferred location for the air valve is downstream
of the 908 bend. During the venting process, pressure rise at the air
valve should not be too great. An excessive pressure rise upstream of
the check valve may overcome the downstream static head, allowing
the check valve to crack open and air to escape into the downstream
line. Often a typical large-orifice air valve would be installed to
ensure adequate air venting capacity.
When all air or gas has finally been purged, the air valve will shut and
the differential flow rate upstream and downstream of the check valve
will abruptly be eliminated. Downstream of the initially closed check
valve, velocity is zero and the upstream velocity will be given by the
operating point on the pump H—Q curve at the moment of valve
closure (Fig. 18.4). This upstream velocity will still be greater than
design velocity. An indication of the sudden head rise associated with
the process of velocity equalisation after air valve closure can be
obtained as follows.
From the Cþ characteristic arriving upstream of the check valve
(Fig. 18.5):
Vu=s þ g=aHu=s ¼ Jþ
where Vu=s and Hu=s are given essentially by the pump operating point
(Fig. 18.4), ignoring riser losses and head loss through the air valve.
From the C characteristic arriving on the downstream side of the
check valve (Fig. 18.5):
0 g=aHd=s ¼ J
where Hd=s is the prevailing head downstream of the check valve. If no
other pump is already operating, Hd=s will be the static head. Where
pumps are already running, Hd=s will be given by the appropriate
point on the system curve.
When flow is abruptly equalised, common velocity V is given by:
V ¼ ð Jþ þ JÞ=2 ¼ ðVu=s þ g=aðHu=s Hd=s ÞÞ=2
The transient pressure H when upstream and downstream water
columns meet is given by:
H ¼ ð Jþ JÞ=ð2g=aÞ ¼ ½Vu=s þ g=aðHu=s þ Hd=s Þ=ð2g=aÞ
329
Pressure transients in water engineering
Hd/s
Air valve
Hu/s
Discharge
NRV header
Riser
Vu/s Vd/s = 0.0
Dx + Dx
C+ C–
Dt
From pump
or
H ¼ aVu=s =ð2gÞ þ ðHu=s þ Hd=s Þ=2 ð18:1Þ
This head rise H can be very considerable and unacceptable.
Figure 18.6 shows an example of predicted velocity in the filling riser
of Fig. 18.1. The pump reaches design speed in 1.0 s and a maximum
velocity of 3.2 m/s is developed. As the riser fills, head at the pump
delivery increases and pump operation gradually moves back up its
performance curve (Fig. 18.4), to produce a reduced flow rate. Velocity
continues to decrease gradually until the riser is filled and the air valve
shuts. The velocities upstream and downstream of the air valve con-
nection are abruptly equalised and a high-frequency transient occurs.
Eventually a steady velocity of <1.0 m/s is produced. Head changes
at the air valve connection for this 1 s pump start are contained in
Fig. 18.7. Only small changes occur until the air valve shuts when
upstream and downstream velocities are equalised. An abrupt head
rise develops reaching a peak of 180 mAOD. The high-frequency
transient dissipates quite rapidly to leave a steady operating head of
330
Air and gas
2.5
2.0
Velocity (m/s)
1.5
1.0
0.5
0.0
0.003
0.382
0.760
1.139
1.517
1.896
2.274
2.653
3.031
3.410
3.788
4.167
4.545
4.924
5.302
5.681
6.059
6.437
6.816
7.194
7.573
7.951
8.330
8.708
9.087
9.465
9.844
10.223
10.601
10.980
11.358
–0.5
Time (s)
180 AV connection
160
140
120
Head (mAOD)
100
80
60
40
20
0
0.003
0.422
0.840
1.258
1.677
2.095
2.513
2.932
3.350
3.768
4.187
4.605
5.023
5.441
5.860
6.278
6.696
7.115
7.533
7.951
8.370
8.788
9.206
9.625
10.043
10.462
10.880
11.298
11.717
12.135
12.554
–20
Time (s)
331
Pressure transients in water engineering
120
Head (mAOD)
100
80
60
40
20
0
8.297
8.326
8.356
8.386
8.416
8.446
8.476
8.506
8.536
8.566
8.595
8.625
8.655
8.685
8.715
8.745
8.775
8.805
8.835
8.864
8.894
8.924
8.954
8.984
9.014
9.044
9.074
9.103
–20
Time (s)
332
Air and gas
valve body downstream of the door. When the valve is later opened
fully after the air valve shuts, the door will travel into a region partially
filled with air. Damping of the door motion can be reduced by this air
pocket, allowing the door to open more rapidly and causing it to
violently strike its stop, possibly resulting in damage.
333
Pressure transients in water engineering
120
100
80
60
40
20
0
0
9
18
27
36
45
54
63
72
81
90
99
108
117
126
135
144
153
162
171
180
189
198
207
216
225
234
Time to reach design speed (s)
Fig. 18.9. Maximum head during riser filling plotted against ramp slope
334
Modulo 1
Cruce
ro Cla
rines 3
ue
anq
Est rines
Cla
360
Clarines
4+
2
875
4+
2
.
V.C cion
" P. enta lara
∆10
re
40 m ta C
Alimrines = 28 510
na
Cla ro L Sannta dento
1+
U
"Ace Pla tame
io
∆10 0 Tra 476
1
R
6+
10m 1
L = 15 se
ero bal
" Ac Em Clara
∆12 ta
1
6 San
2
3
4
5
0
Tom
∆20" Aduccion
El Yai
335
Pressure transients in water engineering
3.5
3
Velocity (m/s)
2.5
1.5
0.5
0
0.003
0.325
0.647
0.968
1.290
1.611
1.933
2.255
2.576
2.898
3.220
3.541
3.863
4.184
4.506
4.828
5.149
5.471
5.793
6.114
6.436
6.758
7.079
7.401
7.723
8.044
8.366
8.688
9.009
9.331
9.652
–0.5
Time (s)
from Rio Unare and to transfer this water at a rate of 350 litres/s either
to Santa Clara storage reservoir or directly to the nearby Santa Clara
treatment works. Pressure rating of the existing 20 in (508 mm) steel
rising main is 16 bar(g) and the pump duty head is 110 m.
336
Air and gas
10.670 m
3 tonne manual hoist
8.510 m
400 on 400 on TFA
6V AV
400 on
TPI RCV
SW
6.995 m
6.760 m
6.470 m
770 SD
floor opening
min W.L.
0.000 m
Inlet pipe
(by others)
–2.030 m
–3.164 m Strainer
–3.430 m
200
Head (mAD)
150
100
50
0
2.5
2.540
2.580
2.619
2.659
2.699
2.739
2.779
2.818
2.858
2.898
2.938
2.978
3.017
3.057
3.097
3.137
3.176
3.216
3.256
3.296
3.336
3.375
3.415
3.455
3.495
3.535
3.574
3.614
–50
Time (s)
337
Pressure transients in water engineering
n
ai
m
y
er
liv
de
N
D
0
50
g
tin
is
Ex
TP1 TP2
400 DN
300 DN BF Existing
GV 300 DN 300 DN platform
TFA GV
4400 Approx.
Back flushing pipework
300 DN
New access
walkway
750
Pumpset Pumpset
1900
No. 1 No. 2
338
Air and gas
El Yai/duty start
250
Ch. 0.0 km Riser AV
200
150
Head (mAD)
100
50
0
2.5
2.669
2.838
3.007
3.176
3.346
3.515
3.684
3.853
4.022
4.191
4.360
4.529
4.698
4.868
5.037
5.206
5.375
5.544
5.713
5.882
6.051
6.220
6.390
6.559
6.728
6.897
7.066
7.235
7.404
7.573
–50
Time (s)
pressure rise forces the non-return valve open and flow is established in
the downstream rising main.
Following a pump trip, as piezometric level falls, the butterfly valve
automatically reopens allowing air to re-enter the riser.
339
Pressure transients in water engineering
While the liquid column is static, gas may evolve from the column of
sewage upstream of the check valve both under septic conditions and
also because the column is under vacuum pressure. When the pump
is started and flow commences in the riser, the gas volume is
compressed. As pressure within this volume rises to match the piezo-
metric head downstream of the check valve, gas will be pushed into
the downstream line through this valve. When all gas has been removed
from upstream of the non-return valve then the column of liquid will
impact against the valve door, throwing it open violently into the
downstream fluid which is a mixture of sewage and gas. Damage to
valve doors and stops has occurred as a consequence.
340
Air and gas
10 m
DN 300
Gas pocket
Hydraulic gradient at
time of flow reversal
8m
Riser
1000 m
M
of the liquid component. The amount of reversed flow and its duration are
modest in these circumstances. A short duration of reversed flow allows
only a modest time for flow acceleration to occur, resulting in a relatively
small reversed velocity. Inertial head rise is a function of this maximum
reversed velocity and so the peak pressure will also be relatively modest.
When a volume of gas is present, an additional storage component
exists, namely the volume available as the gas mass is compressed by
the increasing pressure produced by the reversed flow. This component
can be written by differentiating equation (12.2) as:
dVol=dpabs ¼ 1=nVol=pabs ð18:3Þ
Allowing reversed flow to develop, the gas volume is compressed with
rising pressure. The hydraulic gradient flattens and the acceleration
towards the gas pocket reduces and finally ceases as a positive gradient
is established. The peak pressure within the gas pocket will depend on
the maximum reversed velocity and the initial volume of gas present.
The duration of reversed flow and the maximum reversed velocity
will increase as the initial volume of gas is increased. Deceleration of
this reversed velocity will occur when the combined effects of rising
pressure in the gas pocket plus developing pipeline resistance produce
a positive hydraulic gradient within the pipeline. Even a relatively
341
Pressure transients in water engineering
160
n = 1.4
140 n = 1.2
n = 1.0
Maximum pressure head (mWG)
120
100
80
60
40
20
0
0 2 4 6 8 10 12
Initial air/gas volume (m3)
342
Air and gas
relatively small initial gas volume and that this peak can be considerably
in excess of the maximum with no gas present.
The maximum amplification of 560% shown in Fig. 18.17 over the case
when no air pocket is present is greater by a factor of about 2 than the
260% amplification reported by Escarameia (2005) for an air pocket
located partway along a pipeline. Explanation for this difference may lie
in the position of the air pocket. The compression wave produced by
amplification of pressure at a point within the pipeline, will propagate
along the main. On encountering a closed valve for example, then
amplitude of the compression wave will be doubled against the valve,
resulting in an overall amplification of 520%. The 560% amplification
was for the case of an air pocket being compressed at a closed valve
while the 260% amplification was at a pocket along the pipeline.
343
Pressure transients in water engineering
20
15
10
0
0.004
0.487
0.970
1.453
1.936
2.420
2.903
3.386
3.869
4.352
4.835
5.318
5.801
6.284
6.767
7.250
7.733
8.216
8.700
9.183
9.666
10.149
10.632
11.115
11.598
12.081
12.564
13.047
13.530
14.013
14.496
14.980
–5
Time (s)
Fig. 18.18. Head plotted against time for small orifice valve
during the initial reversed flow interval. Figure 18.19 illustrates that the
pressure rise on valve closure is a more significant component of the
total pressure rise.
For the system under consideration these curves demonstrate that for
a given flow capacity, peak pressure is at a maximum for a particular
n = 1.0
50
Vol = 10 litres Vol = 80 litres
40 Vol = 100 litres Vol = 150 litres
30
Pressure head (mWG)
20
10
0
0.003
0.297
0.591
0.885
1.179
1.473
1.767
2.061
2.355
2.649
2.943
3.237
3.531
3.825
4.119
4.413
4.707
5.001
5.295
5.589
5.883
6.177
6.471
6.765
7.059
7.353
7.647
7.941
8.235
8.529
8.823
–10
–20
Time (s)
Fig. 18.19. Head plotted against time for larger outflow orifice
344
Air and gas
initial air volume and that the peak pressure increases as air flow
capacity rises.
345
Pressure transients in water engineering
Control panel
Washwater
booster set
8.00
7.10
Air relief
(existing)
0.34
TP. 2
–3.20
TP. 1
TP. 1
2500 N T.F.A.
300/2500 N T.F.A.
Part Section A–A
use and also with both mains in service. In each case two duty pumps
were tripped together. Figures 18.26 and 18.27 show these comparisons.
The precise volume of gas evolved is not possible to predict with any
accuracy nor is the rate of gas reabsorption with increasing pressure.
346
Air and gas
These parameters may vary to some extent from one test to another. A
modest increase in gas release rate is sufficient to produce a peak head
at the pumping station >100 m (Fig. 18.28).
The development of high pressures in the presence of a gas pocket is
essentially a matter of the time available to accelerate flow towards the
gas pocket. If the pocket is small then following flow reversal it does not
347
Pressure transients in water engineering
60
Observation d/s of pump
50
30
20
Prediction halfway along rising mains
10
0
0.008
0.815
1.621
2.428
3.234
4.040
4.847
5.653
6.460
7.266
8.072
8.879
9.685
10.492
11.298
12.104
12.911
13.717
14.524
15.330
16.136
16.943
17.749
18.556
19.362
20.168
20.975
21.781
22.588
23.394
24.200
25.007
–10
Time (s)
take long for the reversed flow to compress the gas and increase pressure
within the gas volume. As pressure increases, the hydraulic gradient
causing flow to accelerate towards the pocket is first flattened and
then starts to decelerate the reversed flow. If the time available to
accelerate flow is small then the maximum reversed velocity will also
be more limited and the peak upsurge pressure will tend to remain
more modest. As gas volume is increased, the rate of pressure rise in
60
Elevation (mAOD)
50
40
30
20
10
0
0 100 200 300 400 500 600 700 800 900
–10
Chainage (m)
348
Head (mAOD) Pressure (bar g)
0
1
2
3
4
5
6
7
8
9
10
0
10
20
30
40
50
60
70
80
90
0.027 07/22/99 14:29:52
07/22/99 14:31:48
0.729
07/22/99 14:33:45
1.431
07/22/99 14:35:41
2.133
07/22/99 14:37:37
2.835 07/22/99 14:47:41
3.537 07/22/99 14:49:38
4.239 07/22/99 14:51:34
4.941 07/22/99 14:53:30
5.643 07/22/99 14:55:44
6.345 07/22/99 14:59:15
7.047 07/22/99 15:01:12
7.749 07/22/99 15:03:08
8.451 07/22/99 15:05:04
9.153 07/22/99 15:07:02
9.855 07/22/99 15:08:58
10.557 07/22/99 15:10:54
Time
11.259 07/22/99 15:12:51
Time (s)
07/22/99 15:14:47
11.961
07/22/99 15:21:38
12.663
07/22/99 15:23:34
13.365
07/22/99 15:25:31
14.067 07/22/99 15:27:27
14.769 07/22/99 15:29:23
15.471 07/22/99 15:31:19
PS head
21.087 07/22/99 15:47:01
07/22/99 15:48:57
Air and gas
flow and with compression of the gas mass, pressure rises increasingly
the gas pocket becomes smaller and the hydraulic gradient accelerating
The gas mass is still not sufficient to effectively cushion the reversed
flow persists for longer, allowing a larger reversed velocity to develop.
349
Pressure transients in water engineering
50
Head (mAOD)
40
30
20
10
0
0.027
0.675
1.323
1.971
2.619
3.267
3.915
4.563
5.211
5.859
6.507
7.155
7.803
8.451
9.099
9.747
10.395
11.043
11.691
12.339
12.987
13.635
14.283
14.931
15.579
16.227
16.875
17.523
18.171
18.819
19.467
Time (s)
steeply to produce the peak pressure. With a large gas pocket there is
ample time to create the maximum reversed velocity. Eventually the
hydraulic gradient flattens and flow starts to decelerate. Gas volume
remains relatively large and rate of pressure rise is modest, producing
a gradual deceleration of flow and a lower peak. With such relatively
80
Head (mAOD)
60
40
20
0
0.027
0.675
1.323
1.971
2.619
3.267
3.915
4.563
5.211
5.859
6.507
7.155
7.803
8.451
9.099
9.747
10.395
11.043
11.691
12.339
12.987
13.635
14.283
14.931
15.579
16.227
16.875
17.523
18.171
18.819
19.467
–20
Time (s)
350
Air and gas
351
Pressure transients in water engineering
140
120
100
Head (mAOD)
80
60
40
20
0
0.002
0.134
0.267
0.400
0.532
0.666
0.798
0.931
1.064
1.197
1.330
1.462
1.595
1.728
1.861
1.994
2.126
2.259
2.392
2.525
2.658
2.791
2.923
3.056
3.189
3.322
3.455
3.587
3.720
3.853
3.986
4.119
–20
–40
Time (s)
closed and velocities upstream and downstream of the air valve connec-
tion were equalised when an abrupt upsurge developed. Peak head was
around 170 mAOD. Actual measurements indicated that maximum
pressure could reach almost 20 bar(g). The resulting elastic oscillation
decays rapidly over about 1 s. Predicted velocity variations during this
event are shown in Fig. 18.30 for the pump delivery (DN 350), at
the elevated air valve (DN 500) and at the start of the rising main
Ayr solo pump op/blockage and restart
3.0
2.5
2.0
Velocity (m/s)
1.5
1.0
0.5
0
0.002
0.139
0.277
0.415
0.553
0.691
0.828
0.966
1.104
1.242
1.379
1.517
1.655
1.793
1.931
2.068
2.206
2.344
2.482
2.620
2.757
2.895
3.033
3.171
3.308
3.446
3.584
3.722
3.860
3.997
4.135
–0.5
Time (s)
352
Air and gas
30
25
Air volume (l)
20
15
10
0
0.002
0.136
0.271
0.405
0.539
0.674
0.808
0.943
1.077
1.212
1.346
1.481
1.615
1.750
1.884
2.019
2.153
2.287
2.422
2.556
2.691
2.825
2.960
3.094
3.229
3.363
3.498
3.632
3.767
3.901
4.036
–5
Time (s)
353
Pressure transients in water engineering
The peak pressures recorded were sufficient to damage air valves and
pressure transducers. Initially 10 bar(g) rated transducers were installed
and destroyed by the pressure transients. Replacement instruments
were rated to 20 bar(g) and survived. Air valve floats were crushed
and Fig. 18.32 shows the buckled elevator between floats of an APEX
valve typical of the damage caused.
75 mm recirculation
line with valve
Elevated section of 300 mm line
Tied bend
300 mm valve
Pump
Strainer
Crossings below drainage channels, DN 450 seaward
water main and main road steel main
Treated effluent lagoon
Diffuser
354
Air and gas
Steel straps
Ground level
450∆ DI pipeline
Concrete footing
300 ¥ 450 concentric taper
355
Pressure transients in water engineering
356
Air and gas
70
Tied joint
Ch –1600 m
60 Ch –3200 m
50
40
Head (mAD)
30
20
10
–10
Time (s)
70
Invert level
Max. head
60 Min. head
50
40
Level (mAD)
30
20
10
–10
–20
Chainage (m)
357
Pressure transients in water engineering
180
Tied joint
160 Ch –1600 m
Ch –3200 m
140
120
Head (mAD)
100
80
60
40
20
0
Time (s)
primed main. In general, design of a pipe bend would be based upon the
equation:
Resultant thrust ðNÞ ¼ 2gA½h þ V 2 =ð2gÞ sinð=2Þ ð18:4Þ
where is the total angle of turn and h is the pressure head (metres). It
may be necessary to add other forces such as the weight of a water
column in some instances. The pressure head and velocity should be
the maximum occurring during the surge event.
If gas and vapour are present initially then following pump start rising
pressure will cause vapour cavities to collapse, gas to be compressed and
progressively transported and reabsorbed into the flow. Pump start
studies were carried out using a wide range of cavity volumes and
other parameters such as the rate at which gas could be reabsorbed
by the flow. It was found only too easy to create pressures substantially
greater than the 7 bar(g) design limit. Figure 18.37 shows typical time
histories of head at the tied joint and along the main. The most
likely cause of failure was considered to be movement of the 908
bend allowing the upper section of pipe at the Viking-Johnson coupling
to be pulled out. The coupling was designed to accommodate normal
thermal movements of up to possible 10 mm. It would seem likely
that the 908 had been subject to shock pressures above the 7 bar(g)
limit on various occasions when the pump was started, leading to
eventual failure of the coupling.
358
19
Relief valves
359
Pressure transients in water engineering
Deflector hood
Upper plate
Spring
Isolating valve
From pipeline
360
Relief valves
Isolating
valve
Connection
Outflow
Pipeline
361
Pressure transients in water engineering
Oil reservoir
Spring
Cylinder
Piston
Discharge to
wet well
Valve giving control on outflow
from lower side of cylinder
From pipeline
relief to occur. As the valve starts to close, the restricted outflow of oil
from the lower part of the cylinder only permits a slow closure. This
damping effect also eliminates any risk of resonance developing.
Where damping is not included in the valve then it is necessary for
the period of oscillation of the valve’s spring/mass system to be much
less than that of the pipeline in order to avoid occurrence of pressure
pulsations. The moving parts should be as light as possible to reduce
inertia and ensure freedom of movement. Lightness is also favourable
for stable operation.
362
Relief valves
The connection from the main to the relief valve should be as short as
possible to reduce any risk of resonance and to minimise head loss.
Hmin
Discharge (l/s)
Fig. 19.4. Head plotted against flow relationship for relief valve
363
Pressure transients in water engineering
exceeds that force exerted by the compressed spring. When the valve
opens fully (Fig. 19.4), outflow Qr is given by:
p
Qr ¼ constant 2 ðH zÞ ð19:1Þ
where Qr is discharge through the relief valve, is valve diameter, H is
head at the valve connection and z is the valve elevation. Also the
upstream and downstream characteristics (Fig. 19.5) yield the following
equations for the quasi-invariants:
Vu=s þ g=aH ¼ Jþ and Vd=s g=aH ¼ J
Wavefront
Hu/s
Hset
Relief valve
Hmax
Qr
Vu/s Vd/s
C+ C–
Dt
Dx Dx
364
Relief valves
Wet well
To wet well
Check valve
To atmosphere
Isolating valve To rising main
365
Pressure transients in water engineering
366
Elevation (mPWD) Head (mPWD)
0
50
100
150
200
250
300
350
0
50
100
150
200
250
300
350
0 0.112
1153.36 11.088
2367.56 22.064
3532.11 33.040
hmin
hmax
4751.58 44.016
5983.65 54.992
il (mPWD)
7155.67 65.968
8448.01 76.944
9605.88 87.920
10823.28 98.896
11979.81 109.872
13198.89 120.848
14418.39 131.823
15634.59 142.799
16858.59 153.775
18029.24 164.751
175.727
Time (s)
186.703
Chainage (m)
20559.99
21843.30 197.679
23005.95 208.655
219.631
30552.27
Ch. 20 km
Ch. 10 km
35607.23
329.391
Relief valves
367
Pressure transients in water engineering
250
200
Head (mPWD)
150
100
50 d/s pmps
Ch. 10 km
Ch. 20 km
0
0.112
11.872
23.632
35.392
47.152
58.912
70.672
82.432
94.192
105.952
117.712
129.471
141.231
152.991
164.751
176.511
188.271
200.031
211.791
223.551
235.311
247.071
258.831
270.591
282.351
294.111
305.871
317.631
329.391
341.151
352.911
Time (s)
Fig. 19.9. Head variations following pump trip with pumping station relief valve
250
200
Elevation (mPWD)
150
100
50 il (mPWD)
hmax
hmin
0
0
1153.36
2367.56
3532.11
4751.58
5983.65
7155.67
8448.01
9605.88
10823.28
11979.81
13198.89
14418.39
15634.59
16858.59
18029.24
19328.61
20559.99
21843.30
23005.95
24300.18
25532.26
26754.46
28037.77
29262.87
30552.27
31841.67
33075.07
34372.03
35607.23
Chainage (m)
Fig. 19.10. Envelope curves following pump trip with relief valve
368
Relief valves
200
150
Discharge (l/s)
100
50
0
123.088
124.544
126.000
127.455
128.911
130.367
131.823
133.279
134.735
136.191
137.647
139.103
140.559
142.015
143.471
144.927
146.383
147.839
149.295
150.751
152.207
153.663
155.119
156.575
158.031
159.487
160.943
162.399
163.855
165.311
166.767
–50
Time (s)
369
Pressure transients in water engineering
Pipeline
set pressure imposed by the pressurised annulus, is the liquid within the
system able to push the rubber sleeve aside, allowing flow to pass
through the slotted openings from the upstream to the downstream
side of the cylinder. In this way the internal pressure in the pipeline
can be limited in the same manner as achieved by the more common
pressure relief valve. Figure 19.13 depicts the sequence of operation
370
Relief valves
1000
800
600
400
200
0
0 5 10 15 20 25 30 35 40 45 50 55
% (pressure rise/set pressure)
of the regulator. Advantages are its quiet rapid operation and positive
shutoff, even if foreign matter should become trapped between the
core and the tube. There is no risk of oscillation because of the
design of the regulator and the absence of any spring-mass system.
Discharge Qr through a regulator can be established using the
formula:
p
Qr ¼ Fg Ps ð19:3Þ
where Fg is a flow factor and Ps is the set pressure at which the regulator
is set to start opening. Valves of flow factor as a function of percentage
rise above the set pressure are shown for one model of Grove regulator
in Fig. 19.14.
371
372
Control shaft
Piston
(closing
Spiral member)
extension pipe Outer
cylinder
Dashpot A
Cam C
Connecting pipe
Inner
body
Control valve B
Dismantling ring Bridge piece
Outer cylinder
373
Pressure transients in water engineering
to pass over the parapet of a road bridge (Fig. 19.17). This situation was
not considered acceptable, as one day a tourist bus stopped on the
bridge with its roof vents open. Two turbines tripped and the resulting
relief valves’ operation soaked the passengers, with water flooding out of
the bus when the door was opened. A longer tailwater channel would
have solved the problem but this was not a practical proposition. The
1200 mm diameter outlet hood was modified to include a truncated
cone containing orifices. This caused outflow to break into a large
number of jets and increased aeration and energy dissipation. While
not completely eliminating the problem of water spilling over the
parapet, some improvement has been achieved.
Other options were tried in an effort to reduce severity of the prob-
lem. One attempt adjusted the relief valve movement so that it could
only open to possibly 70%. A test using a partially opened valve setting,
and without any preliminary analysis of likely effects, produced results
described by the operators as ‘horrifying’. The entire power station
vibrated and shuddered, with penstock pressures reaching 500 psig,
which is above design maximum pipeline pressure. The relief valve
was designed to pass flow without large energy loss within the stream-
lined passages. Part-open operation creates conditions for substantial
head loss within the valve itself. Given the overall amount of energy
involved, around 40 MW/machine, any attempt to dissipate flow
energy within the valve through turbulence will inevitably result in
substantial vibration and must be avoided.
374
Relief valves
This example has illustrated that the solution does not just involve
choosing a relief valve but also involves conditions downstream of
the valve to ensure that adequate provision is made for adequate dissi-
pation of energy associated with the outflow and safe disposal of
discharging water. In the present instance, a greater length of tailwater
channel would have avoided the problems experienced.
Snap-through position
Initial position
Vu/s Vd/s
Pipeline
Dx + Dx
C+ C–
Dt
375
20
Check valve dynamics
376
Check valve dynamics
377
Pressure transients in water engineering
Discharge header
Isolating valve
Isolating valve
Suction
branch
Delivery branch
Check valve
Pump
Suction header
Check valve
Isolating valve
378
Check valve dynamics
250
Pressure head (mWG)
200
50
Static head Pump trip
= 75 mWG
0
0 1 2 3 4 5
Time (s)
–50
379
Pressure transients in water engineering
record shows steady pumping conditions before pump trip and also the
transient pressure head following pumping failure. Total period of the
record was around 5 s. Static pressure head at the check valve was
70 mWG. Steady running pressure head while pumping rose to
78.83 mWG. Subsequent to pump trip, pressure fell to a minimum of
46.61 mWG as flow reversed. Head then started to increase slightly
prior to check valve closure. After closure, transient pressure rose
abruptly to a head of 350 mWG. Elastic oscillations produced further
severe pressure peaks and after 5 s the maximum head still remained
above the height of the building. The high transient pressures caused
burst pipe connections, vibrations and noise, with the result that
occupants had to be evacuated. This is an instance of an unsuitable
check valve selection resulting in delayed valve closure and allowing
development of substantial reversed flow in the system. The remedy
was to install a valve possessing fast-closing attributes or alternatively
to use a valve with a damped closure.
380
Check valve dynamics
381
Pressure transients in water engineering
Head in discharge
manifold during
pump deceleration
Pumping head as M
pump decelerates
Suction well .
Discharge branch
M
Discharge manifold
NRV or header
Failing pump
Steady pumping
Suction branch hydraulic gradient
Hydraulic gradient
L during deceleration
Fig. 20.4. Hydraulic gradient in a pumping station shortly after pump trip
382
Check valve dynamics
check valve may fall below the upstream head in say a wet well or
suction main and forward flow may recommence with the check
valve partially reopening. A check valve installed within a bypass
around the pump(s) would also open under these circumstances. This
can occur several times before the valve finally remains shut. The
following cases illustrate both types of valve reopening.
DN 100 NRV
Intake
WPL Duoglide
DQC 100/125
Non-reflecting boundary
at start of rising main
383
Pressure transients in water engineering
1.2
1.0
0.8
0.6
Velocity (m/s)
0.4
0
0.001
0.015
0.029
0.043
0.057
0.071
0.085
0.100
0.114
0.128
0.142
0.156
0.170
0.184
0.198
0.212
0.227
0.241
0.255
0.269
0.283
0.297
0.311
0.325
0.339
0.354
0.368
0.382
0.396
0.410
–0.2
–0.4
Time (s)
station. Weir Pumps Limited Duoglide units type DQC 100/125 were
employed. A pressure vessel was connected 8 m downstream of the
second pump branch. Figure 20.6 shows the changing velocity just
downstream of the check valve of a failing pump and also downstream
of the valve of the remaining operating pump. After trip at 0.1 s the
velocity falls steeply and reverses at 0.285 s. A reversed velocity of
0.3 m/s is achieved at the time the check valve closes after about
0.37 s. Velocity becomes zero and remains so until pressure wave reflec-
tions cause head across the valve to allow a brief interval of reopening
before the valve finally shuts. Velocity in the operating pump branch
increases until the first check valve shuts and then falls with a
modest oscillation developing. Changing head downstream of the
check valves can be seen in Fig. 20.7, with an abrupt upsurge when
the check valve of pump No. 1 closes. Peak head is almost 120 m
above datum. Steady pumping head is around 84 mAOD. Head down-
stream of the second pump falls to a lesser extent. Wave reflection
causes head downstream of the failing pump to reach a minimum and
it is at this time that the check valve briefly reopens. A smaller head
variation occurs at the operating pump. The reason that the check
valve is able to reopen under such high head conditions is because
the pump is still rotating at around 85% of its design speed at the
time of valve closure. The reduced head developed by the pump is
still sufficient to reopen the check valve when downstream head falls.
384
Check valve dynamics
140
Delivery No. 1
120 Delivery No. 2
100
Head (mAOD)
80
60
40
20
0
0.001
0.015
0.029
0.043
0.057
0.071
0.085
0.100
0.114
0.128
0.142
0.156
0.170
0.184
0.198
0.212
0.227
0.241
0.255
0.269
0.283
0.297
0.311
0.325
0.339
0.354
0.368
0.382
0.396
0.410
Time (s)
385
Pressure transients in water engineering
180
160
140
120
Head (mAD)
100
80
60
40
PS No. 5A
20 Ch. 7 km
0
0.327
11.118
21.909
32.700
43.491
54.282
65.073
75.864
86.655
97.446
108.237
119.029
129.820
140.610
151.401
162.192
172.983
183.774
194.565
205.356
216.147
226.937
237.728
248.519
259.310
270.101
280.892
291.683
302.473
313.264
324.055
Time (s)
386
Check valve dynamics
1.5
Velocity (m/s)
0.5
0
0.327
5.232
10.137
15.042
19.947
24.852
29.757
34.662
39.567
44.472
49.377
54.282
59.187
64.092
68.997
73.902
78.807
83.712
88.617
93.522
98.427
103.332
108.237
113.142
118.048
122.953
127.858
132.763
137.668
142.572
–0.5
–1
–1.5
Time (s)
causes head at the pumping station to rise, with the result that the flow
through the pumps decelerates and the check valves reclose quietly.
The necessary pressure vessel capacity is reduced by reopening of the
check valves. Figure 20.10 shows air volume expanding quite rapidly
40
Air volume (m3)
30
20
10
0
0.327
10.791
21.255
31.719
42.183
52.647
63.111
73.575
84.039
94.503
104.967
115.431
125.896
136.360
146.823
157.287
167.751
178.215
188.679
199.143
209.607
220.070
230.534
240.998
251.462
261.926
272.390
282.854
293.318
303.781
314.245
324.709
Time (s)
387
Pressure transients in water engineering
after pumps were tripped. When the check valves reopen, only a modest
additional expansion occurs, with air volume reaching a maximum of
60 m3 . Effectively, the recommencement of flow through the check
valves has replaced the need for further outflow from the pressure
vessel.
1.6
1.4
Reversed velocity (m/s)
1.2
1.0
0.8
0.6
0.4
0.2
1 2 3 4
Number of operating pumps
388
Check valve dynamics
dx/dt = +a
Differential head across
valve door after closure a/gVo
Vo V=0 V=0 Vo
Door closed
389
Pressure transients in water engineering
390
Check valve dynamics
+t
V, H
C Co
Dt
T/2
Dt Vs, HT/2
Dt
Steady-state region
V = Vs (constant)
T/2
Dt
Cs
391
Pressure transients in water engineering
With reference to Fig. 20.13, over a time interval 0 to T=2, the quasi-
invariant becomes:
Js ¼ Vs þ sg=aHs f T=ð4DÞVs2
and
HT=2 ¼ Hs f T=ð4DÞVs2 ð20:3Þ
Proceeding along the characteristic Co from time T=2 to a solution
point on the boundary at time T, the quasi-invariant value is:
ð
Jo ¼ Vs sg=aHT=2 g S dt ð20:4Þ
with the integral evaluated between the limits T=2 and T. Applying the
D’arcy equation once more p and approximating average velocity, over
the interval T=2 T, by ðVVs Þ then:
S ¼ f VVs =ð2gDÞ
and at the solution point on the model boundary:
V sg=aH ¼ Jo ¼ Vs sg=a½Hs f T=ð4DÞVs2 f T=ð4DÞVVs
or rearranging,
V½1 þ f T=ð4DÞVs sg=aH ¼ Vs sg=a½Hs f T=ð4DÞVs2 ð20:5Þ
From within the region being modelled the characteristic C arrives at
the solution point with a quasi-invariant value:
V þ sg=aH ¼ J
392
Check valve dynamics
Limit of modelling
H(i+1)
i
H(i)
Dt
System curve
u/s system
Q(i) Q(i+1) Q
Dx
393
Pressure transients in water engineering
394
Check valve dynamics
395
Pressure transients in water engineering
396
Check valve dynamics
The value quoted above is sometimes stated as the maximum head rise
possible at a check valve on closure. This is only the case as far as the
initial head rise is concerned, as observations have revealed consider-
ably greater head rise after initial closure as described in section
20.10.2. At the upstream valve face, the corresponding head drop
will be:
Hv ¼ aVo =g
397
Pressure transients in water engineering
2.5
1.5
Velocity (m/s)
0.5
0
0.001
0.109
0.216
0.323
0.431
0.538
0.645
0.753
0.860
0.967
1.075
1.182
1.289
1.397
1.504
1.611
1.718
1.826
1.933
2.040
2.148
2.255
2.362
2.470
2.577
2.684
2.792
2.899
3.006
3.114
3.221
–0.5
Time (s)
–1
–1.5
the cavity collapsed was of a similar order as the initial reversed velocity
at the time of initial closure.
Miller (c. 1970) observed behaviour of a Glenfield and Kennedy
Swing Check Valve subsequent to pump trip from an initial steady
pumping head of 14.33 m. Following check valve closure, a peak pres-
sure head of 405 ft or 123.44 mWG was measured at the inlet to the
check valve while the corresponding peak pressure at the valve outlet
was 270 ft or 82.3 mWG. Figure 20.15 shows predicted coastdown of
velocity at the pump discharge. After closure, the velocity oscillates
as a vapour cavity upstream of the check valve repeatedly opens and
shuts. The corresponding pressure head is shown in Fig. 20.16 together
with the observed values. A more detailed picture of head changes can
be seen in Fig. 20.17. It will be noted that head variations downstream
of the valve occur with higher frequency and dissipate more rapidly.
Modest vapour cavities were predicted to occur on the downstream
side of the valve but these were smaller and relatively short-lived
under the generally higher pressure prevailing in the downstream pipe-
line. The head changes upstream of the valve take place over a longer
time span with larger cavities forming. From this and other studies it has
been noted that transient pressures upstream of a check valve can
exceed the corresponding downstream pressures often by a considerable
margin. Changes in diameter along the pump suction and discharge
branches can aggravate conditions.
398
Check valve dynamics
Pressure head at inlet and outlet of Glenfield and Kennedy Swing Check Valve
following pump trip at time = 0.1 s
160
140
Valve inlet
123.44 metres – peak head
120 measured at valve inlet
Pressure head (mWG)
Valve outlet
100
82.3 metres – peak head
80 measured at valve outlet
60
40
20
0
0.001
0.121
0.240
0.359
0.478
0.597
0.716
0.835
0.954
1.073
1.192
1.311
1.430
1.549
1.668
1.788
1.907
2.026
2.145
2.264
2.383
2.502
2.620
2.700
2.858
2.977
3.096
3.216
3.394
3.4
3.5
–20
Time (s)
120
Valve outlet
Pressure head (mWG)
100
80
60
40
20
0
2.570
2.596
2.623
2.649
2.675
2.792
2.728
2.755
2.781
2.808
2.834
2.861
2.887
2.914
2.940
2.967
2.993
3.019
3.046
3.072
3.099
3.125
3.152
3.178
3.205
3.231
3.258
3.264
3.311
3.367
3.383
–20
Time (s)
Fig. 20.17. Head plotted against time as check valve shuts (detail)
399
Pressure transients in water engineering
400
Check valve dynamics
Pump
NRV
401
Pressure transients in water engineering
402
Check valve dynamics
100
80
Elevation (mASL)
60
40
i.l. mASL
20 hmax
hmin
0
0
575.6
1247.1
1726.7
2398.3
3069.8
3645.3
4316.9
4892.4
5468.0
6043.6
6619.2
7098.8
7578.5
8154.1
8633.7
9209.3
9784.9
10 360.5
10 840.1
11 415.7
12 087.2
12 662.8
13 334.3
13 909.9
14 389.5
14 965.1
15 540.7
16 212.2
Chainage (m)
403
Pressure transients in water engineering
0.6
0.4
0.2
0
0.093
2.139
4.185
6.231
8.277
10.323
12.369
14.415
16.461
18.507
20.553
22.599
24.645
26.691
28.737
30.783
32.829
34.875
36.921
38.967
41.013
43.059
45.105
47.151
49.197
51.243
53.289
55.335
57.381
59.427
–0.2
–0.4
Time (s)
60
50
Air volume (m3)
40
30
20
10
Vessel
0
0.093
3.906
7.719
11.532
15.345
19.158
22.971
26.784
30.597
34.410
38.223
42.036
45.849
49.662
53.475
57.288
61.101
64.914
68.727
72.540
76.353
80.166
83.979
87.792
91.605
95.418
99.231
103.045
106.858
110.671
114.484
118.297
Time (s)
404
Check valve dynamics
80
Head (mASL)
60
40
20
0
0.093
4.092
8.091
12.090
16.089
20.088
24.087
28.086
32.085
36.084
40.083
44.082
48.081
52.080
56.079
60.078
64.077
68.076
72.075
76.074
80.073
84.072
88.071
92.070
96.069
100.068
104.068
108.067
112.066
116.065
120.064
Time (s)
Fig. 20.22. Head variations after pumping failure, upstream part of system
When pumps are restarted they will commence operation against a low
static head and with the start-up compression wave being reflected from
the first closed check valve.
Figure 20.22 shows the predicted head variations at the pumping
station also upstream and downstream of the first NRV, after pumps
are tripped. A smooth fall in head at the pumping station is followed
by a very modest oscillation around an eventual static head of about
45 mASL. Upstream of the first NRV the head drop initially follows
the same form as at the pumping station until a wave reflection from
downstream produces a head rise beginning after about 14 s. A
second reflection causes head to fall once more and at around 37 s
the NRV shuts. The rate of head drop increases upstream of the
valve after closure and head falls to a minimum at about 58 s.
Thereafter an oscillation occurs between the pumping station and
the NRV, with maximum amplitude at the closed valve. Downstream
of the NRV, head conditions are essentially the same as upstream.
After closure a modest high-frequency oscillation occurs in that section
of main between the now closed first and second NRVs.
Figure 20.23 depicts changing head downstream of the second and
third NRVs. After closure of the check valves a modest high-frequency
oscillation is developed between the shut second and third valves.
405
Pressure transients in water engineering
100
80
Head (mASL)
60
40
20 d/s 2 NRV
d/s 3 NRV
0
0.093
4.092
8.091
12.090
16.089
20.088
24.087
28.086
32.085
36.084
40.083
44.082
48.081
52.080
56.079
60.078
64.077
68.076
72.075
76.074
80.073
84.072
88.071
92.070
96.069
100.068
104.068
108.067
112.066
116.065
120.064
Time (s)
Fig. 20.23. Head variations after pumping failure, downstream part of system
406
Check valve dynamics
Downstream
reservoir
Piezometric gradient
without in-line NRV
Reversed flow
Double-orifice air valve
Piezometric gradient
with in-line NRV
Closed NRV
PS
Fig. 20.24. Hydraulic gradients with in-line check valve and bypass
407
Pressure transients in water engineering
408
21
Check valve characteristics
409
Pressure transients in water engineering
15
P Case 3: Amri wafer check valve
(bar(g)) 400 mm NB Series 2000
10 8.24 bar(g) max Fig. No. 1H3T2KF – heavy springs fitted
410
Check valve characteristics
411
Pressure transients in water engineering
Flow
Door
(d) External lever and axial spring (e) External torsional spring
(f) Air chamber + lever and weight (g) Cylinder and oil reservoir
412
Check valve characteristics
m = hinge friction
rm = distance to centre of gravity
rp
Centre of gravity
+V
Seating angle = b
q dq/dt, d2q/dt 2
a = opening angle
Ws = buoyant weight
413
Pressure transients in water engineering
414
Check valve characteristics
415
Pressure transients in water engineering
door moves freely towards its seat until the final stages of closure. The
piston can also be mounted internally to make contact with the door
itself.
A relatively slow-closing valve may allow high reversed velocity to
occur and the damping device may be unable to control the hydrody-
namic forces and can become unstable if not carefully chosen. A
quick-closing free-acting valve may be considered preferable in some
circumstances.
416
Check valve characteristics
Diaphragm Stopper
Air plug
Hinge pin Cover
Outlet body
Hinge pin
Inlet body Outlet body
Door
Flow
Flow
Door
Seat
Seat ring
Face ring
Inlet body Face
417
Pressure transients in water engineering
Flow
Seat ring
Stop
Open
Opening angle
(b) Wafer-type body tilting disk
418
Check valve characteristics
Closing disk
Ductile iron/nitrile coated
+
Flexible hinge
Damped closure
419
Pressure transients in water engineering
door travel and thus closure time. Since no shaft projects through the
valve body, the options of having external lever and weight or springs
to assist closure are not available with these valves.
In some valves, such as shown in Fig. 21.6b, the option of installing a
‘bottom buffer’ is available. This is in the form of a piston which projects
into the valve near to the bottom of the seat. When the closing door
makes contact with the piston over the final stages of closure,
the piston gradually moves back into an oil-filled cylinder, effectively
throttling any reversed flow which develops and controlling decelera-
tion rate dV=dt. This modification produces a damped closure over
the final stages of movement and it is necessary to know characteristics
of the oil cylinder and piston for any modelling exercise involving valve
closure.
Being rubber-coated, door closure will not be accompanied by the
sharp crack which is present with metal-to-metal seating.
420
Check valve characteristics
Door closed
90∞
421
Pressure transients in water engineering
+ +
Open
Spring Open Shut
(a) Demag type DRV-Z nozzle check valve (b) Nozzle check valve
e.g. SOCLA type 402
Spring
Open Shut
(c) Alstom Clasar check valve
422
Check valve characteristics
and pushing the door on to its seat. The distance which the valve door
has to travel is generally small, allowing a relatively prompt closure to be
achieved.
Closure performance is governed by the equation:
X
axial forces ¼ 0
These pertinent forces are:
ð
Hydrodynamic force ¼ p dA
Friction force ¼ R
Inertial force ¼ Ws =g dv=dt
Buoyant weight ¼ Ws (þ sign in upwards flow and sign in down-
wards flow)
Spring force ¼ x (where is in N/mm and x is movement in mm)
Good closure performance can be achieved with a short stroke and a
small disk mass. An adjustable spring force gives flexibility in applica-
tion of the valve.
As with other internal spring-assisted valves it is deemed inadvisable
to use this type of valve when the liquid carried contains material which
may become entangled with the spring and impede closure. Some of
these valves, such as the Demag nozzle type (Fig. 21.9a), have a stream-
lined body designed to maximise recovery of kinetic energy developed
as flow passes through the throat of the valve. This streamlining
tends to produce a longer valve than other types which do not recover
the same amount of energy. There has to be a trade-off between energy
savings in the case of the more expensive longer-body type for which
increased space must be reserved within a pumping station and a
cheaper valve in which greater head loss occurs but which is shorter
and does not require the same amount of space. The globe-shaped
version is an example of a shorter body type (Fig. 21.9b), with a door
travel of about DN/4.
The Alstom Clasar wafer type illustrated in Fig. 21.9c, is a variation
on the axial movement valve. The valve door comprises a set of
concentric streamlined rings having a central shaft and axial spring.
When closed, these moving rings seal against a set of fixed rings. On
opening, the spring is pulled into tension by the hydrodynamic force
acting on the set of moving rings. Travel of the moving rings is typically
<DN/10 and valve length is <DN/2. When a pump is tripped and flow
423
Pressure transients in water engineering
Spring
Cover
Guide bush
+ Flow direction
Body
Disc Seat
424
Check valve characteristics
Valve open
where flow deceleration is relatively low at the valve, for instance solo
duty pump installations with modest static lift, so that the system-
imposed dV=dt is small. Advantages are that the ball is self-cleaning
with virtually no risk of blockage. This valve pattern finds application
in sewage systems and other areas such as mine drainage.
425
Pressure transients in water engineering
Direction of
reversed flow
closure the valve acts much as a free-acting valve in that it does not
greatly affect the rate of flow change. Subsequent to initial closure,
deformation of the sleeve over a significant time produces a more
gradual flow deceleration than say a swing-check valve, and this part
of the process is more akin to that of a damped valve.
This valve is suitable for lines conveying slurry, sludge, lime, chemical
slurries, paper plant waste and raw sewage, with the flexibility of the
rubber sleeve allowing the valve to seal even in the presence of
gravel. When installed in appropriate conditions, the valve will produce
a quiet closure. If subject to excessive transient differential pressure
across the sleeve, the possibility of damage exists through the rubber
sleeve being pushed inside out or ‘inverted’. To prevent damage it is
essential to choose a valve having the ability to resist the maximum
differential pressure likely to be experienced across the valve. A range
426
Check valve characteristics
Hd/s
Valve body
Hu/s
Reversed flow
V = velocity
+
Vol in
DH = Hd/s – Hu/s
427
Pressure transients in water engineering
alter and with it the volume of liquid on the inlet and outlet sides of the
sleeve. The total volume within the valve body remains constant — that
is, inlet volume þ outlet volume ¼ constant. It is necessary to know the
relationship between inlet or outlet volume and differential head. This
information can be obtained from, for example, experimental measure-
ment or finite-element analysis of the sleeve under load. Given the
relationship between inlet volume and differential head and ignoring
inertial effects within the sleeve, then over a time increment of the
analysis:
VA ¼ dVolin =dt ¼ dVolin =d H d H=dt ð21:1Þ
where Volin is the inlet volume, H is the differential head ¼
Hd=s Hu=s , V is the velocity and A is area of flow. Digitising the
relationship between inlet volume and differential head (Fig. 21.13),
and interpolating using initial values and averaging over the time
increment t then with dVolin =d H Volin =ðHÞ ¼ gradient of
Volin =H curve and suffix ‘o’ denoting conditions at the start of the
time increment:
ðV þ Vo ÞA=2 ¼ gradient ½Hd=s Hu=s ðHd=so Hu=so Þ=t
From the characteristic paths:
Jþ ¼ V þ g=aHu=s and J ¼ V g=aHd=s
or
Hd=s Hu=s ¼ a=g½2V ð Jþ þ JÞ
Writing, constant ¼ 2a=g gradient=ðA tÞ and rearranging then:
V ¼ constant½ Jþ þ J þ g=aðHd=so Hu=so Þ=ð2 constant 1Þ Vo
ð21:2Þ
Corresponding values of Hu=s and Hd=s are obtained by substituting the
value of V into the quasi-invariant equations.
Prior to initial sleeve closure and subsequent to any sleeve inversion,
the relationship between differential head and velocity should be
expressed in the form:
H ¼ KL V 2 =ð2gÞ
where KL is a loss coefficient. The sleeve may become inverted if the
differential pressure exceeds a maximum permitted limit. Transient
analysis can be used to establish this limit for any system so that an
appropriate sleeve can be used.
428
Check valve characteristics
1.5
Velocity (m/s)
0.5
0
0.001
0.121
0.240
0.359
0.478
0.597
0.716
0.835
0.954
1.073
1.192
1.311
1.430
1.549
1.668
1.788
1.907
2.026
2.145
2.264
2.383
2.502
2.621
2.740
2.859
2.978
3.097
3.216
3.336
3.455
3.574
–0.5
–1
–1.5
Time (s)
429
Pressure transients in water engineering
14
12
10
Head (mAD)
6 Series 1
Series 2
4
Series 3
0
0.001
0.119
0.237
0.354
0.472
0.589
0.707
0.825
0.942
1.060
1.177
1.295
1.413
1.530
1.648
1.765
1.883
2.001
2.118
2.236
2.353
2.471
2.589
2.706
2.824
2.942
3.059
3.177
3.294
3.412
3.530
3.647
–2
Time (s)
Fig. 21.15. Head upstream and downstream of a sleeve valve after pump trip
of the higher rate of deceleration for the much larger reversed velocity
achieved in the 20 m sleeve.
If the reversed velocity is too great, a sleeve may be deformed to an
extent that it becomes inverted. Under these conditions the valve may
no longer remain fully closed and could allow some reversed flow to
continue. The rate of this continued reversed flow will depend upon
the inverted shape and stiffness of the sleeve. To illustrate conditions,
a simple system was considered. A uniform pipeline of length 980 m was
chosen which rises continuously from invert level þ0.0 mAD down-
stream of a pumping station to an elevation þ30.0 mAD at chainage
980 m. Discharge level was set at þ40 mAD. Pipe diameter was
400 mm. Suction well level was 4.0 mAD.
Pump start produced an eventual steady velocity of 1.9 m/s. Alterna-
tive sleeve stiffnesses were examined, ranging from 60 m to 100 m.
After pump trip, velocity at the pump and in the pipeline decelerates
and reverses after about 4 s as shown in Fig. 21.16 for a 60 m sleeve. An
initial reversed velocity is developed, reaching around 0.3 m/s. As the
sleeve deforms in an endeavour to mobilise sufficient strength to
decelerate the reversed flow, it becomes inverted after around 6 s as
shown in Fig. 21.17.
Figure 21.18 shows the head variations predicted for this event. After
trip, head downstream of the pump falls to 10.0 mAD and thereafter
increases gradually. After around 4 s the sleeve becomes closed and flow
430
Check valve characteristics
1.5
Velocity (m/s)
0.5
0
0.005
0.353
0.701
1.049
1.396
1.744
2.092
2.440
2.788
3.136
3.484
3.832
4.180
4.526
4.875
5.223
5.571
5.919
6.267
6.615
6.963
7.311
7.659
8.007
8.354
8.702
9.050
9.398
9.746
10.094
10.442
–0.5
–1
Time (s)
40
30
Volume (l)
20
10
0
4.091
4.170
4.248
4.327
4.405
4.483
4.562
4.640
4.719
4.797
4.875
4.954
5.032
5.111
5.189
5.267
5.346
5.424
5.503
5.581
5.659
5.738
5.816
5.895
5.973
6.051
6.130
6.208
6.287
6.365
6.443
6.522
–10
–20
–30
Time (s)
431
Pressure transients in water engineering
50
40
30
Head (mAD)
20 Series 1
Series 2
10 Series 3
Series 4
0
0.005
0.475
0.946
1.416
1.886
2.357
2.027
3.218
3.768
4.218
4.749
5.109
5.650
6.120
6.590
7.061
7.581
8.002
8.472
8.942
9.413
9.883
10.354
10.824
11.294
11.765
12.285
12.706
13.176
13.646
14.117
14.587
–10
–20
Time (s)
0.8
Pump speed/design speed
0.6
0.4
0.2
0.0
0.005
0.480
0.955
1.431
1.906
2.381
2.857
3.332
3.807
4.283
4.758
5.233
5.708
6.184
6.659
7.134
7.610
8.085
8.560
9.036
9.511
9.986
10.461
10.937
11.412
11.887
12.363
12.838
13.313
13.789
14.264
–0.2
–0.4
–0.6
–0.8
–1.0
Time (s)
432
Check valve characteristics
0.5
0
0.005
0.603
1.200
1.798
2.396
2.994
3.592
4.189
4.787
5.385
5.983
6.581
7.178
7.776
8.374
8.972
9.570
10.167
10.765
11.363
11.961
12.559
13.156
13.754
14.352
14.950
15.548
16.145
16.743
17.341
17.939
18.537
–0.5
Time (s)
this oscillation taking place around the zero flow condition. A similar
velocity variation occurs at the valve inlet and with no continuous
reversed flow there is no motivation for reversed pump rotation to
occur. The predicted oscillation was due to substantial deformation of
the sleeve. Under the action of transient effects within the pipeline,
only a 10 litre margin remains before inversion could occur. Head
changes are shown in Fig. 21.21, with peak head of þ75 m at the
40
30
20
10
0
0.005
0.612
1.220
1.818
2.430
3.043
3.650
4.258
4.866
5.473
6.081
6.688
7.296
7.904
8.511
9.119
9.726
10.334
10.942
11.549
12.157
12.764
13.372
13.080
14.587
15.195
15.802
16.410
17.018
17.625
18.233
18.840
–10
Time (s)
433
Pressure transients in water engineering
434
Check valve characteristics
435
Pressure transients in water engineering
436
Check valve characteristics
60 60 qm qc 6
40 40 4
30 30 3
Vm Vc
20 20 Hm 2
Hc
10 10 1
0 0 0
0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
70 70 7
qm qc
60 60 6
40 40 4
Vm Vc
30 30 3
20 20 Hm 2
Hc
10 10 1
0 0 0
0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
30 30 3
Vm Vc
20 20 Hm Hc 2
10 10 1
0 0 0
0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
437
Pressure transients in water engineering
2.5
400 300
DN 100 Delft hydraulics data
250
Keystone Fig. 85 DN 200 200
Delft hydraulics data
2.0 150
1
Effect of shifting centre of 2
3
gravity of valve door
4
Magnitude of reversed velocity at closure (m/s)
1.5
Glenfield M1
0.0
0 5 10 15
Deceleration gradient |dV/dt| (m/s2)
438
Check valve characteristics
1.5
Magnitude of reversed velocity at closure (m/s)
Glenfield recoil
series 5100
DN 300
1.0
DN 200
0.5
DN 150
DN 100
0
0 5 10 15
Deceleration gradient |dV/dt| (m/s2)
Fig. 21.25. Performance curves for Glenfield recoil and tilting disk valves
439
Pressure transients in water engineering
DN 100
from Delft hydraulics data
1.5
Magnitude of reversed velocity at closure (m/s)
DN 200
from Delft hydraulics data
1.0
0.5
0
0 5 10
Deceleration gradient |dV/dt| (m/s2)
440
Check valve characteristics
2.4
2.3 DN 450
DN 750
DN 600
2.2
1.8
1.7
Magnitude of reversed velocity at closure (m/s)
1.6
441
Pressure transients in water engineering
2.5
DN 800
DN 700
2.0 Velocity when fully
open = 1.5 m/s
Demag
DRV-B
Magnitude of reversed velocity at closure (m/s)
1.5
DN 600
Measured DN 800
medium spring
DN 500
1.0
DN 400
DN 300
0.5
0.0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28
Deceleration gradient |dV/dt| (m/s2)
442
Check valve characteristics
1.3
1.2
1.1
No spring, Vo = 0.0 m/s
Magnitude of reversed velocity at closure (m/s)
1.0
0.8
0.6
0.5
0.4
0.3
Stiff spring, Vo = 3.0 m/s
0.2
0.1
0.0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28
Deceleration gradient |dV/dt| (m/s2)
443
Pressure transients in water engineering
changes with varying spring stiffness. In this figure the initial velocity is
representative of the linear spring stiffness, with a broken or missing
spring requiring zero velocity to maintain the valve in the open position.
A stiff spring by contrast needs a flow velocity of 3.0 m/s to fully open
the valve. Not surprisingly, the stiffer spring is able to close the valve
quickly but the absence of a spring requires differential pressure
1.2 Gestra RK 56
DN 200 weak spring
1.1
1.0
0.9
Demag DRV-B
DN 300 weak spring
Magnitude of reversed velocity at closure (m/s)
0.8
0.7
0.6
0.5
Demag DRVg
DN 200 weak spring
0.4
0.2
0.1
Demag DRV-B
DN 300 strong spring
0.0
0 2 4 6 8 10 12 14 16 18 20 22 24 26
Deceleration gradient |dV/dt| (m/s2)
444
Check valve characteristics
445
22
Flexible pipe
446
Flexible pipe
FRP 0.3
RPM 0.3
PVC 1.5
HDPE 4.0
447
Pressure transients in water engineering
(a) (b)
448
Flexible pipe
449
Pressure transients in water engineering
or
"f Dm =s ¼ 1 Dm =Dmax
giving,
Dm =Dmax ¼ 1 "f Dm =s ð22:2Þ
where Dm is the original diameter of the round pipe and Dmax is the
maximum deformed diameter.
One of the most difficult tasks is prediction of deformation and strain
over time. A relationship which according to Moser (1981) may be used
for deflections of up to 10% in PVC pipes is:
" ¼ ðs=Dm Þ½3 y=Dm =ð1 2 y=Dm Þ ð22:3Þ
where y is the vertical change in diameter, Dm is mean diameter and "
is the maximum strain. This equation does not include the influence of
internal pressure. If a cylinder is deformed and therefore not perfectly
circular it is inherently weaker and will collapse at a differential pressure
somewhat below that of a circular pipe in the same circumstances.
When a pipe is relatively unconstrained, for example above ground or
when placed in loose soil, the maximum strain may be:
Design value ¼ 0.5%
in emergency situations ¼ 1.5%,
where repeated surging occurs ¼ 1.3%
The ultimate strain is defined as the initial strain at first installation
which would cause failure after 50 years. For PVC this has been
suggested as 1%. Allowable design strain is half the ultimate strain.
Moser could find no basis for this 1% ultimate strain value.
The relationship between pipe ring deformation and tangential strain
in the pipe wall has been expressed in a different form by Molin (1981)
as:
"r ¼ pDm =ð2sEÞ Df ðy=Dm Þðs=Dm Þ ð22:4Þ
where "r is the tangential strain in the pipe wall, y=Dm is the deforma-
tion of the pipe ¼ , s is the wall thickness and Df is the deformation
factor.
This equation represents a linear combination of strains due to
internal pressure pDm =ð2sEÞ and strain due to elliptical deformation
Df ðy=Dm ÞðsDm Þ. If the deformation is elliptical then Df ¼ 3. In
most practical conditions 4 Df 6, although values of Df in some
circumstances may range from 2.85 to 10 or greater. In circumstances
450
Flexible pipe
Table 22.2. Deflection and shape v. embedment (Bishop and Lang (1984))
451
Pressure transients in water engineering
P P
s s Dy
Dy
Dm Dm
(a) (b)
452
Flexible pipe
length of pipe:
Sr ¼ EI=r3 ð22:10Þ
3
where I ¼ s =12.
As regards the allowable buckling pressure within a deformed pipe
Timoshenko (1936) and others have studied the elastic stability of
cylindrical shells and predicted buckling of a round pipe using the
following equation:
pb ¼ 24EI=½ð1 2 ÞoD3 ð22:11Þ
where pb is the critical buckling pressure; E is the modulus of
elasticity of the pipe material; is Poisson’s ratio for the pipe
material, for HDPE ¼ 0:4; I is moment of inertia of the original
circular cross-section ¼ s3 /12; and oD is mean outside pipe diameter.
Where out-of-roundness has occurred, a reduction factor c is applied,
giving critical pressure pc ¼ pb c. Alternative expressions for the reduc-
tion factor c have been developed. These reduction factors may include
corrections both for ellipticity as well as variations in pipe wall thick-
ness. Considering flexible pipes ‘creep’, effects were accounted for by
using the long-term secant modulus and long-term Poisson’s ratio for
pipe material. For a plastic pipe displaying ‘creep’, any calculation of
long-term deflection or buckling strength requires use of the ‘creep’
modulus. For HDPE the creep modulus at 50 years is typically 0.2—
0.3 GN/m2 . When calculating initial deflection and strength for
HDPE, the short-term modulus is 0.8—1.0 GN/m2 . The modulus for
polyethylene may be calculated from equation (22.12):
E ¼ 10:0½5:88 0:067 log10 ðhoursÞ 0:6 log10 ð8FÞ ð22:12Þ
The reduction factor c is given variously by:
c ¼ ðoDmin =oDmax Þ3 ðsmin =smax Þ3 ; that is, c always <1 ð22:13aÞ
Subscripts ‘max’ and ‘min’ denote maximum and minimum values of
diameter and wall thickness.
Ignoring the influence of variations in pipe wall thickness, Allman
(1975) reported an experimentally derived expression for Alydn R ‘D’
PE pipe giving reduction factor ca as:
ca ¼ ðoDmin =oDmax Þ4:62 ð22:13bÞ
Again, neglecting the effect of wall thickness changes, Gaube et al.
(1974) gave a reduction factor cg :
cg ¼ ½ð1 Þ=ð1 þ Þ2 3 ð22:13cÞ
453
Pressure transients in water engineering
454
Flexible pipe
(%) 1 2 3 4 5 6
resistance and r the pipe radius. The modulus of soil reaction E0 is not a
constant but among other factors depends upon the stiffness of the pipe
used. Spangler recommended that E0 values should be based upon
experience and judgement.
Equation (22.14) is acknowledged to be a relatively crude tool for
calculating collapse pressure.
Compacted soil embedment factors in the pipe zone increase pipe
resistance to hydrostatic collapse. For properly compacted embedment
areas a supporting factor fos ¼ 3 may be applied. Taprogge (1981)
calculated critical pressure pc using the relationship:
pc ¼ po fc fos ð22:15aÞ
where po is the buckling pressure for a round pipe, fos is a factor of
safety ¼ 3 for properly compacted embedments, and fc is an ovality
correction factor. Table 22.3 lists values of fc as a function of and
shows how fc varies with deflection.
For polyethylene sewer pipes Taprogge indicates a factor of safety
ð fosÞ ¼ 1:5 so that the allowable differential pressure becomes:
pex p ¼ pc =1:5 ð22:15bÞ
For GRP pipes the vertical deflection is limited to 5%.
Greatorex (1981) has supplied data for GRP/RPM pipes and some of
his comments were as follows.
(a) In the USA GRP pipes are supplied with stiffness EI=D3 ¼ 1280 N/m2
and this is considered close to the optimum for most applications
with a 5% deflection being normal. In Europe GRP pipe is supplied
with a range of stiffness.
(b) Pipes of stiffnesses < 1000 N/m2 cannot be safely installed at any
depth where there is a risk of full vacuum occurring.
(c) Minimum stiffness with 2D cover should be the greater of
1100 N/m2 or that obtained from the equation:
p
pex p ¼ f32E0 rtEI=D3 g1=fos ð22:16Þ
where pex is external pressure; p is the internal pressure (ve if
vacuum); E0 rt is the tangent value of the resultant modulus of
455
Pressure transients in water engineering
456
Flexible pipe
Kurt and Mark (1981) quoted the equation of Slocum for collapse
pressure of thermoplastic elliptical sections as:
pc ¼ pb 9 ð22:18Þ
where ¼ Dmin =Dm ¼ 1 .
At 5% eccentricity ¼ 0:05 and ¼ 0:95, giving p=po ¼ 0:63%,
and at 10% eccentricity ¼ 0:10 and ¼ 0:90, giving p=po ¼ 39%.
These results do not consider the stiffening effect of the reducing
minimum radius of curvature of the deforming section.
Investigations have shown that it is possible to place pipes in clay
surround and the necessary interaction between pipe and soil will be
obtained. However, the support from this clay will be less than if the
trench were refilled with sand. The pipe has to be designed for its
surroundings.
Leaving aside the matter of factor of safety, the equation for long-
term buckling pressure pex p is:
p
pex p ¼ f32E0 trEI=D3max g ð22:16Þ
This may be written in terms of the original round pipe diameter Dm
thus:
p
pex p ¼ ð1 "f Dm =sÞ1:5 f32E0 trEI=D3m g ð22:19Þ
For long-term buckling calculations either the long-term ‘creep’
modulus E50 may be used giving stiffness ¼ E50 s3 =ð12D3m Þ. Alterna-
tively the long-term elastic stiffness ¼ 0:7Eo s3 =ð12D3m Þ can be used.
The factor 0.7 is an adjustment for reduction in modulus over time
and may also include some compensation for out-of-roundness. The
long-term modulus is thus taken as 70% of the short-term modulus.
Equation (22.19) would indicate that if:
1 "fð50Þ Dm =s ¼ 0
then the pipe will buckle under a small differential pressure pex p, giving:
"fð50Þ ¼ s=Dm ¼ Dfð50Þ ð50Þ s=Dm
or
Dfð50Þ ð50Þ ¼ 1; giving ð50Þ ¼ 1=Dfð50Þ ð22:20Þ
In the long term Dfð50Þ may reach values of 10—15, corresponding with
values of ð50Þ in the range 6—10%. Practical experience has shown
that thin-walled pipes have collapsed after some years, even at low
backfill heights — that is, low values of pex p.
457
Pressure transients in water engineering
458
Flexible pipe
Pex
d
s
P
Dm
Pex Pex
P P
d s
Pex
where h is the horizontal radial deviation from the circular shape and
v is the corresponding vertical deviation, then if circumference is of
constant length:
2r ¼ ðr þ h þ r þ vÞ ¼ 2r þ ðh þ vÞ
giving,
h ¼ v ¼ ð22:24Þ
Substituting in the formula for cross-sectional area A then:
A ¼ ab ¼ ðr þ Þðr Þ ¼ ðr2 2 Þ ð22:25Þ
Irrespective of the value and sign of the cross-sectional area is
smaller than for the circular form with ¼ 0.
Differentiating with respect to :
dA=d ¼ 2
and
1=A dA=d ¼ 2=ðr2 Þ ¼ 2=r2 ð22:26Þ
From theory of elastic stability described by Timoshenko (1936), the
relationship between deflection and buckling pressure is given by:
=i ¼ pb =ð pb pex þ pÞ ð22:27Þ
459
Pressure transients in water engineering
where differential pressure normal to the pipe wall is pex p, pex being
the external pressure and p is the internal pressure; is deformation and
i is the initial deformation. pb is the differential pressure which
produces buckling in a circular section and is given by:
pb ¼ 2E=ð1 2 Þðs=Dm Þ3
or
hb ¼ 2E=½ gð1 2 Þðs=Dm Þ3 ð22:11Þ
where hb ¼ pb =ðgÞ. When pex p ¼ pb the pipe collapses completely.
Writing:
¼ i pb =ð pb pex þ pÞ
and differentiating with respect to p:
d=dp ¼ i pb =ð pb pex þ pÞ2
Assuming this is applicable to the approximately elliptical shape
under consideration, then:
1=A dA=dp ¼ 1=A dA=d d=dp ¼ 2=r2 ½i pb =ð pb pex þ pÞ2
or
1=A dA=dp ¼ 2fi pb =rg2 =f pb pex þ pg3 ð22:28Þ
In considering the speed of propagation of a pressure wave in an
elliptically deforming HDPE outfall pipeline, Larsen (1976) derived
the expression:
ðA Ao Þ=Ao ¼ 3fi pb =½rð pb pex þ pÞg2
Differentiating with respect to p gives:
1=A=dA=dp ¼ 6fi pb =rg2 =f pb pex þ pg3 ð22:29Þ
Equations (22.28) and (22.29) have different coefficients, reflecting the
alternative assumptions with respect to the shape of deformed section
but the form of expression is the same.
Whereas the expression for changing area, at constant shape,
due to hoop stress yields a constant gradient 1=A dA=dp ¼ Dc1 =ðsEÞ,
the gradient of area with pressure under changing shape is itself
pressure dependent. The pressure wave propagation rate is not a
system constant but varies with changing flow conditions in a similar
manner to inclusion of gas/vapour content variation with pressure.
460
Flexible pipe
with pressure terms p, pb and pex having been replaced by the corre-
sponding pressure head values h, hb and hex , using h ¼ p=ðgÞ.
The changing shape factor effectively creates more storage within the
pipe than is given by simply using the area change due to hoop stress
variations. The effect is to produce a reduction in wave speed with
falling pressure and therefore a reduced pressure wave amplitude or
‘damping’ of the wave. Figure 22.4 shows the dependence of wave
speed upon ovalisation for HDPE pipes of SDR ratio 17.6 and 26.0
and for initial deformations i of 2% and 4%. Variation of the
coefficient c2 from 2 to 6 was also included. SDR 11 pipes are
sufficiently stiff to display little variation of wave speed with internal
pressure.
Larsen (1976) predicted the propagation of pressure transients
using the irregular grid of characteristics produced by varying wave
speed caused by ovalisation. While quite suitable for the relatively
straightforward case of a single pipeline without internal boundary
conditions or branches, this approach becomes more awkward for
general use. Just as the option exists of separating air or gas from the
water component when predicting wave propagation following air/gas
release under low pressures, so the same approach can be used to
separate the ovalisation factor from the area change produced by
hoop stress changes. This may be achieved by ‘collecting’ the ovalisa-
tion effect at computing sections in the same way as gas or vapour is
numerically collected. Alternatively the ovalisation effect can be
analysed as a separate pressure-dependent component distributed
along the pipeline.
Larsen (1976) was primarily interested in the use of HDPE pipes for
undersea outfalls and his concerns were largely with pipes surrounded
by water. In this context he presented curves showing how theoretical
buckling pressure varied when stiffeners in the form of anchoring
weights were placed at intervals along an HDPE pipeline. Inclusion of
pipe ovalisation would also be applicable to above-ground installations
or buried pipe where the response of the pipe was not significantly
influenced by its surroundings. This implies soft soil surroundings
such as marsh, swamp and loose or slightly compacted fill. Influence
461
Pressure transients in water engineering
190
SDR = 26
180 c2 = 2
Di = 0.02
SDR = 17.6
170
c2 = 6
Di = 0.04
160
150
140
Acoustic wavespeed a (m/s)
130
SDR = 26
120 c2 = 6
Di = 0.02
110
100
90 SDR = 26
c2 = 6
80 Di = 0.04
70
60 HDPE pipes
50
40
30
20
10
0
–10 0 10 20 30 40 50
Internal pressure head (mWG)
462
Flexible pipe
463
Pressure transients in water engineering
Cg
1 1.5 2 3 4
40
30
a
b a
c
d b
e c 2
20
d
e
15 1.5
10 1
9 0.9
8 0.8
7 0.7
H/B
6 0.6
H/B
5 0.5
4 0.4
3 0.3
2.5 0.25
2 0.2
1.5 0.15
0.10
0.1 0.15 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.5
Cg
464
Flexible pipe
Horizontal flexible pipeline. con = 1.0, del/D = 2%, h b = 12.3, gas included
70
Series 1
60 Series 2
Series 3
50
Elevation (mAD)
40
30
20
10
0
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10 000 11 000 12 000
–10
Chainage (m)
Fig. 22.6. Envelope curves for horizontal pipe and 2% deflection ratio
465
Pressure transients in water engineering
30
Head (mAD)
20
10
0
2.46
19.68
36.90
54.12
71.34
88.56
105.78
123.00
140.22
157.44
174.66
191.88
209.10
226.32
243.54
260.76
277.98
295.20
312.42
329.64
346.86
364.08
381.30
398.52
415.74
432.96
450.18
467.40
484.62
–10
Time (s)
Fig. 22.7. Head variation for horizontal pipe after pump failure
40
30
20
10
0
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10 000 11 000 12 000
–10
Chainage (m)
466
Flexible pipe
20
Elevation (mAD)
10
0
0
100
200
300
400
500
600
700
800
900
1000
1100
1200
1300
1400
1500
1600
1700
1800
1900
2000
2100
2200
2300
2400
2500
2600
2700
2800
2900
3000
–10
–20
–30
Chainage (m)
Fig. 22.9. Envelope curves for descending outfall after pump failure and 2%
deflection ratio
Pressure transients in water engineering
Descending flexible pipeline. con = 1.0, del/D = 2%, h b = 12.3, gas c = 0.000 01
40
Series 1
30 Series 2
Series 3
20
Elevation (mAD)
10
0
0
100
200
300
400
500
600
700
800
900
1000
1100
1200
1300
1400
1500
1600
1700
1800
1900
2000
2100
2200
2300
2400
2500
2600
2700
2800
2900
3000
–10
–20
–30
Chainage (m)
Fig. 22.10. Envelope curves for descending outfall after pump failure and with
stiffeners
20
15
10
0
0
100
200
300
400
500
600
700
800
900
1000
1100
1200
1300
1400
1500
1600
1700
1800
1900
2000
2100
2200
2300
2400
2500
2600
2700
2800
2900
3000
–5
Chainage (m)
468
Flexible pipe
20
15
10
0
0
100
200
300
400
500
600
700
800
900
1000
1100
1200
1300
1400
1500
1600
1700
1800
1900
2000
2100
2200
2300
2400
2500
2600
2700
2800
2900
3000
–5
Chainage (m)
20
Elevation (mAD)
15
10
0
0
100
200
300
400
500
600
700
800
900
1000
1100
1200
1300
1400
1500
1600
1700
1800
1900
2000
2100
2200
2300
2400
2500
2600
2700
2800
2900
3000
–5
–10
Chainage (m)
Fig. 22.13. Rising main, envelope curves for rigid pipe and gas release
469
Pressure transients in water engineering
20
15
Elevation (mAD)
10
0
0
100
200
300
400
500
600
700
800
900
1000
1100
1200
1300
1400
1500
1600
1700
1800
1900
2000
2100
2200
2300
2400
2500
2600
2700
2800
2900
3000
–5
–10
–15
–20
Chainage (m)
Fig. 22.14. Rising main, envelope curves for rigid pipe and no gas release
are both ignored, the minimum head falls to a much greater extent, as
depicted in Fig. 22.14.
The previous two examples considered produced predictions closely
in line with Larsen’s results when similar pipe properties were used. This
finding would indicate that it is quite feasible to utilise a fixed grid
computational scheme for such studies as well as a variable wave
speed scheme.
470
Flexible pipe
30
d/s PS
Seawall
25 Outfall
20
15
Head (mAD)
10
–5
–10
Time (s)
471
Pressure transients in water engineering
30
Invert lev.
Max. head
25
Min. head
20
15
10
Head (mAD)
–5
–10
–15
–20 Chainage
472
Flexible pipe
473
23
Amplification of transient
pressures
474
Amplification of transient pressures
Distribution mains
At Ad
J+ = Vo2/(1 + 1/2Ad/At)
+T
J– = Vo(1 – 1/2Ad/At)/(1 + 1/2Ad/At) PS
J– = 0
J+ = +Vo
Rising main
+X
J– = –Vo
+X
or
X X
H ¼ a=g ðA JÞ= ðAÞ ð6:4Þ
where wave speed a is considered constant. Suppose a pump trip occurs
at the PS which reduces the initial velocity Vo to zero. The resultant
invariant travelling upstream from the pumping station will have a
value Vo . The invariant from the upstream conduit has value þVo
and the value from the trunk main is zero. Head at the connection
point is then:
H ¼ a=gVo 2At =ð2At þ Ad Þ
where At is the cross-sectional area of the trunk main and Ad is the
cross-section of a distribution main.
Velocity in the trunk main upstream of the connection becomes:
Vt ¼ Vo Vo 2At =ð2At þ Ad Þ
and the resultant invariant propagating upstream of the connection
is:
Jt ¼ Vo ð1 12 Ad =At Þ=ð1 þ 12 Ad =At Þ ð23:1Þ
475
Pressure transients in water engineering
H1
+a
Ho = 0
A1
A2
V1 V=0
+X
J2– = 0
J1+ = V1 + g/aH1
476
Amplification of transient pressures
477
Pressure transients in water engineering
When A2 < A1 then J > J1 and so the pressure wave travelling into a
network containing pipes of reducing diameter may be amplified. For
example, if A2 ¼ 12 A1 then H ¼ 133:0 m/s and V2 ¼ 1:33 m/s, giving
a transmitted wave invariant value ¼ 2.67 m/s as compared with the
incident wave value of 2.0 m/s.
+a –a
+V V=0 –V
+X
J+ = V + g/aH J– = –V – g/aH
–a aV/g +a
H H
+V V=0 –V
478
Amplification of transient pressures
be obtained from:
V ¼ ð Jþ þ JÞ=2 and H ¼ ð Jþ JÞ=ð2g=aÞ
with velocity reduced to zero and head H ¼ a=gV. The head rise is
increased by the intersection of opposing compression waves. When a
pair of opposing rarefaction waves meet, the head drop at intersection
of the waves is likewise greater, as will be illustrated in section 23.5
concerning wellfield transients.
479
Pressure transients in water engineering
350
300
250
Elevation (mAD)
200
150
100
Invert
50 Max. head
Min. head
0
Chainage
300
250
200
Elevation (mAD)
150
100
Invert
50
Max. head
Min. head
0
Chainage
480
Amplification of transient pressures
350
300
250
Elevation (mAD)
200
150
100
Rawyards
50 d/s pump
u/s pump
0
Time (s)
Fig. 23.6. Head variations in the gravity main after pump start
481
Pressure transients in water engineering
Uplands PS
Tank or reservoir
Pumping station
Kirkleatham Lane PS
Redcar
distribution
zone
AC 600f
Yearby TWL = 85.675 mOD
limit of C62. From this pumping station a rising main extends to Yearby
Reservoir at an elevation of 85.675 mAOD. From Yearby a series of
pumping stations and pumping mains serves substantial areas of the
North York Moors. Between the steelworks branch and Kirkleatham
PS lies the Redcar residential area served by a set of branch connections
482
Amplification of transient pressures
150f
Roseberry Road Race
600f AC to 150f course
Yearby Res.
150f 380f 150f
457f 150f
Coast Rd
Redcar Lane 150f
Low Farm Drive 380f
150f 150f
150f
200f
Laburnum Rd
Kirkleatham Lane West Dyke Road 150f
250f 228f
250f
0 1 2 km
from C62. Figure 23.8 shows the networks of pipelines supplying Redcar
as modelled in this study. Branches from C62 comprise two DN 150
mains along Corporation Road, a DN 457 main laid along Trent
Road and a DN 250 pipeline into Kirkleatham Lane itself.
Following earlier studies, a pressure vessel was installed on the down-
stream side of Kirkleatham Pumping Station to protect the rising main
to Yearby from sub-atmospheric pressures. On pump start-up there is an
upsurge on the downstream side of the pumps and a drawdown on the
suction side. The presence of the pressure vessel on the delivery side of
the pumps allows a relatively smooth and gradual rise in pressure to
occur. In contrast, on the suction side the drawdown is very abrupt.
Similarly, when Kirkleatham pumps are switched off there is a
smooth and controlled decline in the head downstream of the pumps
483
Pressure transients in water engineering
110
100
Head (mAOD)
90
80
70
Fig. 23.9. Head downstream of Kirkleatham Lane pumping station after pump trip
regulated by the pressure vessel (Fig. 23.9), while on the suction side
pressure rises sharply (Fig. 23.10).
90
80
70
Head (mAOD)
60
50
40
Fig. 23.10. Head upstream of Kirkleatham Lane pumping station after pump trip
484
Amplification of transient pressures
Metres
50 * * * * *
40
30
20
10
30 Kirkleatham PS
outflow – MLD
3 bursts 1 burst
20
10
Time
0
Fig. 23.11. Head recordings at Kirkleatham indicating pipe bursts had occurred
485
Pressure transients in water engineering
90
80
70
60
Flow (l/s)
50
1
40
30
20
3 bursts burst
10
0
hr:mn 00:00
MT/94 26/10 27/10 28/10 29/10 30/10 31/10 01/11 02/11 03/11 04/11
Time
pumping station. When head rises and flow falls to zero this indicates
that pumps have been tripped. Three bursts were recorded on one
occasion within the distribution zone and a further single burst after
a subsequent pump trip. Flow recordings with the distribution system
(Fig. 23.12) noted an abrupt flow increase which correlated with the
pipe bursts. Figure 23.11 shows the predicted steep head increase
upstream of the pumping station following pump trip, with a
maximum head of around 86 mAOD. Within the distribution system
along Redcar Lane, the corresponding predicted head variation is
shown in Fig. 23.13, with maximum head exceeding 150 mAOD.
Minimum head was also predicted to fall below atmospheric pressure,
which may interfere with sensitive equipment such as that used for
home dialysis which draws supplies direct from mains.
Predicted peak pressure within the distribution system, following
trip of Kirkleatham pumps, is shown in Fig. 23.14 in the form of
486
Amplification of transient pressures
160
120
Head (mAOD)
100
80
60
40
Time (min)
487
Pressure transients in water engineering
12.0
11.0 13.5
13.0 13.9
10.1
10.6
10.7
14.0
14.5
10.8
10.5 15.0
12.5 10.6
15.5
12.9 10.8
14.0
10.9
13.0 15.7
12.0
11.0
11.0
12.0
11.1
11.2
11.5
11.3
11.4
488
Amplification of transient pressures
W-111
Reservoir compound W-105 W-103
DN 600 DI
W-106
W-113
Motorised valves
Line 2
Line 3
Line 4
Line 1
Line 5
W-49
Line 1
DN 100 DI
W-1 W-124
W-125
W-123
Wellhouse
0 1 2 3 km
489
Pressure transients in water engineering
close to the reservoir was 600 mm while at the extremities of the system,
for instance near to well number 1, diameter was 100 mm.
490
Amplification of transient pressures
Control valve
electrically operated
Compensator
Air valve
Flow meter
Flap valve Gate valve
Pressure
gauge
Ground level
Branch pipeline
Non-return valve
Multi-stage pump
Motor
Aquifer
Aquifer
Sump pipe
Clay (impervious)
Bail bottom
491
Pressure transients in water engineering
Borehole losses
Control valve loss
Collector
main loss
Minimum
Pumping head
pumping head
after some time
Reservoir
Wellhouse
Cone of depression
Cone of depression
Aquifer
Bore hole
to permit flows from one group of pumps to pass into the alternative
system of mains in order to reach the reservoirs. Normally, around
50% of the pumps would be in use at any time, yielding about 20 Mld.
492
Amplification of transient pressures
pump trip and that the wellhead and other air valves may come into
operation.
380
370
360
Elevation (mAOD)
350
340
330
320
310
Fig. 23.18. Head at well No. 1 for sequenced start commencing at well No. 1
493
Pressure transients in water engineering
420
400
Elevation (mAOD)
380
360
340
320
Fig. 23.19. Head at well No. 1 for sequenced start ending with well No. 1
somewhat higher peak pressures were predicted. Figure 23.19 shows the
corresponding head variation at well No. 1. Maximum transient pres-
sure exceeded 420 mAD. Other starting arrangements also produced
higher starting pressures than those of Fig. 23.18 and it was concluded
that the optimum arrangement was to first start those pumps towards
the upstream end of the system and progress downstream towards the
reservoirs.
494
Amplification of transient pressures
(i)
Vo Collector main
Vo
Wellhouse Vo 2Vo
Riser Branch
Pump (a)
Ho
(ii)
(ii) Ho
(ii)
(ii)
Ho
(iii)
(iii)
Vo
2Vo
Ho
(ii)
(b)
Ho
(ii)
(iv)
(iv)
(iii)
(iii)
2Vo
Vacuum pressure
(c)
495
Pressure transients in water engineering
360
340
320
mAOD
300
280
260
0 1 2 3 4 5 6
km
Fig. 23.21. Area blackout envelope curves along line No. 1 (initial stages)
Figure 23.20a shows the initial rarefaction waves (i) travelling from
each wellhead along the branch pipeline towards its connection with
the collector main.
In Fig. 23.20b the initial rarefaction waves have reached the
connection point with the collector. A partial reflection in the form
of a compression wave (iii) travels from the connection back towards
each wellhouse while rarefaction waves (ii) travel both upstream and
downstream along the collector from the connection points. These
rarefaction waves produce a fall in head within the collector main.
In Fig. 23.20c opposing rarefaction waves have met in the collector
main between the branch connections, leading to a doubling in ampli-
tude of head drop (iv). There is a strong possibility that severe vacuum
pressures will occur in the collector between adjacent wells.
Simulation of an area power failure affecting all GROUP I pumps
shows that such vacuum pressures develop along line No. 1 as shown
in Fig. 23.21. The profile of line No. 1 shows relatively low elevations
over the middle parts of the main with higher levels both towards
well No. 1, chainage 0 km, and also at the reservoir end of the pipeline,
chainage 6 km.
496
Amplification of transient pressures
440
420
400
380
360
mAOD
340
320
300
280
260
0 1 2 3 4 5 6
km
Fig. 23.22. Area blackout envelope curves for line No. 1 (later stage)
reservoirs ceases and then reverses while positive flows continue from
the higher parts of the system around well No. 1. Air pockets at
these more elevated wellhead air valves will continue to grow while
air is expelled from the lower parts of the system. Air removal through
these lower air valves is facilitated by continued positive flows from the
higher parts of the system around well No. 1 and also by reversed flow
from the reservoirs.
As lower air valves close, compression waves are developed and a
general rise in head is developed over these lower central parts of the
network as shown in Fig. 23.21.
With reversed flow from the reservoirs increasing and air valves
closing, the compression wave propagating upstream travels into
progressively smaller pipes, being amplified until it finally enters the
DN 100 pipes around well No. 1. Closure of final air valves coupled
with the strength of the compression wave produces a maximum
transient head of over 440 mAD at well No. 1 (Fig. 23.22).
497
Pressure transients in water engineering
498
24
Flow instabilities
499
Pressure transients in water engineering
500
Flow instabilities
Constant head
Static head = Ho
Self-acting valve
any stretching which might occur in a cable connection by using the ball
bearings. Unfortunately, a small amount of slack in the chain of ball
bearings was sufficient to confuse the governor about the precise
valve opening. This caused ‘hunting’ of the valve at a frequency
equalling the period of the penstock and produced a self-exiting
oscillation. Reverting to a cable connection eliminated the problem.
The impedance method of analysis can be used to predict the most
critical periods of oscillation in a pipeline network. It is not intended
to discuss this method but the interested reader can find a description
of the method in books, for example by Streeter and Wylie (1967) or by
Jaeger (1977).
To illustrate circumstances in which it is possible for flow within a
pipeline network to become unstable and fail to settle at a steady
operating condition, three examples have been included in this chapter.
The first illustration concerns a pumping installation and the second an
entirely gravity-driven system.
501
Pressure transients in water engineering
East Main
762f
Cranhill PS
457f
Castlemilk Low Castlemilk PS
TWL 128.6 610f
Burnside PS
Castlemilk High 457f
TWL 171.5 Link Main DN 700
Distribution main
to Cathkin DN 200
Link Main DN 700
546f 546f
244f 244f
Auchentibber Udston
TWL 208.8 TWL 181.4
Water tower
Ground tank
Pumping station 546f 546f
Pressure reducing valve
Isolating or fixed throttle valve
0 1 2 3 4 5
are involved, one being the downstream part of the East Main supplied
originally with water from the Loch Katrine system and the other a
gravity-driven network supplied from Daer reservoir to the south.
These usually independent systems operate at significantly different
piezometric levels, with the Daer gravity system around 100 m higher
than the East Main network. To improve security of supply a link
main has been constructed between the two systems. This link runs
502
Flow instabilities
503
Pressure transients in water engineering
100
90
80
70
Generated head (m)
60 Duty point –
link main pumps off
50
Hmax
20
10
Cut impeller
0
0 100 200 300
Flow rate (l/s)
will be to reduce the suction head at existing booster pumps along the
East Main and in particular at Castlemilk pumping station. The steady
operating head for these pumps will increase. Added to this increased
pumping duty at Castlemilk Low station are transient pressure varia-
tions produced when the new booster pump is operated.
504
Flow instabilities
Piezometric gradient
without link main pump
Piezometric gradient
with link main pump
Pumping head
without link
Transient drawdown main pump
Link Main PS
Additional head
+ with link main
pump running
East Main
Castlemilk PS
75
Start of link main pumps
70
Head u/s of Castlemilk PS
65
Head (mAOD)
60
55
50
Fig. 24.5. Head conditions at Castlemilk after link main pump operated
505
Pressure transients in water engineering
110
90
Head (mAOD)
80
Burnside pumps operating
70
60
0 1 2 3 4
Time (min)
Fig. 24.6. Head conditions at Burnside after link main pump is started
506
Flow instabilities
300
M 283.68 mAOD
957.6 litres/s
260
BPC
240 1078.6 litres/s 957.6 litres/s
220
200
M
180
160
140
100
60
40
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28
Chainage (km)
507
Pressure transients in water engineering
Float chamber
Overspill weir
M
Linkage
Throttled connection
between float
chamber and BPC
838 mm f Float
MS main in
2 No. 457 mm f
float operated valves
838 mm f
MS main out
508
Flow instabilities
Actual phase 2
port shape
Initial phase 2
port shape
Actual phase 2
valve travel = 103.2
closure duty. In the event that these downstream valves are closed,
flow into the pipeline below the BPC will fall, causing level in the
chamber to rise. The floats in the BPC will then move upwards,
progressively closing the twin H1 valves and throttling flow entering
the BPC. Each H1 valve contains a set of four openings or ports
through which flow passes. Shape of the ports can be altered to provide
a range of flow v. head loss characteristics to suit particular require-
ments while providing the necessary maximum design discharge
when fully opened. If discharge required should change over time
the port shape can be altered to suit. Phase I of the project had port
shapes as illustrated in Fig. 24.9, with a total flow area 373.16 cm2
when fully opened and with a valve stroke of 85.7 mm. Originally it
was intended that the port shapes would be altered to permit an
increased maximum flow in Phase II through a total ported area of
625.81 cm2 over the same stroke length. However, to maximise
discharge a considerably increased flow area totalling 1574.19 cm2
was actually created, with stroke length being increased to
103.2 mm. These port forms for Phase II are also shown in Fig. 24.9.
Travel of floats within the BPC for the final Phase II arrangement
was 762 mm — that is, from a minimum elevation of 231.058 mAOD
509
Pressure transients in water engineering
7
BPC inlet valve loss coefficient (kV)
as a function of valve stroke
Phase I porting
5
1
Actual Phase II porting
0
0 10 20 30 40 50 60 70 80 90 100 110
Valve stroke (mm)
510
Flow instabilities
511
Pressure transients in water engineering
380
360
340
320
300
280
260
240
0 1 2 3 4 5 6 7 8 9 10
Time (min)
200
150
d/s reservoir
100
50
Pipeline profile
Helvetica Flush left
0 2 4 6 8 10 12 14 16 18 20 22 24 26
Chainage (km)
512
Flow instabilities
Intake tower
Submerged discharge
valve 800 ¥ 600
Trifurcation
513
Pressure transients in water engineering
measurements were made within the ND 2000 penstock both with and
without the small hydro turbine in operation. When flow through the
SDV was 2.7 m3 /s or higher, the following effects were observed:
(a) a pulsating sound within the small hydrostation
(b) the noise of cavitation from within the SDV in the area of the
DN 800 bend (Fig. 24.15)
(c) increased vibration in the small hydrostation
(d) more pronounced turbulence in the SDV sump evidenced by
surface wave action (Fig. 24.16)
(e) large pressure fluctuations recorded on pressure gauges at a hydrant
close to the penstock at the small hydrostation.
Recordings of pressure fluctuations in the penstock showed a
maximum head range of about 9 mWG. While a range of frequencies
of oscillation was recorded, it was noted that some oscillations had a
period of around 1.6 s which is close to the fundamental period 4L=a
of the conduit system as a whole. Other much higher frequencies
were also present.
514
Flow instabilities
515
Pressure transients in water engineering
516
Flow instabilities
90
80
70
60
50
mAD
40
30
20
Trifurcation
SDV tee
10
SDV bend
0
Time (s)
517
Pressure transients in water engineering
Diaphragm
DN 800
Flow
518
Flow instabilities
60.0
59.5
59.0
58.5
58.0
57.5
57.0
56.5
56.0
55.5
SDV tee
55.0 SDV bend
54.5
Time (s)
519
References
520
References
Ellis, J. and Mualla, W. (1986) Selection of Check Valves. Proceedings of the 5th
International Conference on Pressure Surges, BHRA Fluid Engineering
Centre, Hannover, Germany, 22—24 September, pp. 213—222.
Ellis, J. and Tint, K. M. (1976) Simulation of pressure transients at valves and bifur-
cations. International Water Power and Dam Construction, 8:8, 24—27 August.
Esleek, S. H. and Rosser. (1959) Check Valve Water Hammer Characteristics.
Proceedings of an American Water Works Association Society meeting,
November.
Fedosoff, F. A. and Szpak, E. (1979) Cyclic Overpressure Effects on High Density
PE Pipe. Proceedings of Western Canada Water and Sewage Conference,
Regina, Saskatchewan. Engineering Digest, September, 35—38.
Fok, A. T. K. (1980) Design Charts for Surge Tanks on Pump Discharge Lines.
Proceedings of the 3rd International Conference on Pressure Surges, Eng-
land, March, BHRA, Paper J3, pp. 445—472.
Gaube, E., Müller, W. and Falcke, E. (1974) The Statics of Rigid Polyethylene
Drain Pipes. Proceedings of Kunststoffe, 64, April, 193—196.
Glass, W. L. (1980) Cavitation and Corrosion in a Pumping Main. Proceedings of
the 3rd International Conference on Pressure Surges, Canterbury, 25—28
March, Paper H3, pp. 401—414. Cranfield: BHRA.
Gordon, J. L. (2004) The untested computer program. Hydro Review World-
wide, 12:13, July, 32 and 34.
Gordon, J. L. (2006a) Nintendo engineers. Hydro Review Worldwide, 14:4,
September, 46 and 51.
Gordon, J. L. (2006b) The hunting impulse turbine. Hydro Review Worldwide,
14:1, March, 46 and 54.
Gotham, K. V. and Hitch, M. K. (1975) Design considerations for fatigue in a
PVC pressure pipeline. Pipes and Pipelines International, 20:2, February, 10.
Graze, H. R. (1968) A Rational Thermodynamic Equation for Air Chamber Design.
Proceedings of the 3rd Australian Conference on Hydraulics and Fluid
Mechanics, Institution of Engineers Australia, 25—29 November, p. 57.
Graze, H. R. and Forrest, J. A. (1974) New Design Charts for Air Chambers.
Proceedings of the 5th Australian Conference on Hydraulics and Fluid
Mechanics, Christchurch, New Zealand, December 9—13, Vol. II, p. 34.
Graze, H. R., Schubert, J. and Forrest, J. A. (1976) Analysis of Field Measure-
ments of Air Chamber Installations. Proceedings of the 2nd International
Conference on Pressure Surges, London, Paper K2. Cranfield: BHRA.
Greatorex, C. B. (1981) The Relationship between the Stiffness of a GRP Pipe and
its Performance when Installed. Proceedings of an International Conference
on Underground Plastic Pipe, New Orleans, 30 March—1 April (Schrock,
B. J. (ed.)), pp. 117—129. New York: ASCE.
Heinsbroek, A. G. T. J. and Tijsseling, A. S. (1994) The Influence of Support
Rigidity on Waterhammer Pressures and Pipe Stresses. Proceedings of the 2nd
International Conference on Water Pipeline Systems, Edinburgh, UK,
24—26 May, pp. 17—30. Cranfield: BHR Group.
521
Pressure transients in water engineering
522
References
523
Pressure transients in water engineering
524
Further reading
525
Pressure transients in water engineering
526
Index
527
Pressure transients in water engineering
528
Index
529
Pressure transients in water engineering
compression wave propagation, 100, 101, Demag DRV-Z valves, 422, 423
103 demand points
computer programs, 3 mains, 85
concrete pipes, 1011 demand-sensing pressure-reducing valves,
condensers, 6 7677
conduits, 2829, 32, 34 densities
configurations modified Proctor density, 456
break pressure chamber, 508 depth of cover
pipelines, 137139 pipes, 463, 464
pumping stations, 137139, 138 descending outfall, 467468, 467468
connections, 236, 400, 474476, 475 differential orifices, 186, 187
conservation of force, 2324, 24 dip-tubes see central tubes
constrained conditions direct pump start, 140, 141146, 142145
flexible pipes, 454458, 463465 envelope curves, 143, 145
constraint factors multi-pump operation, 143146,
parameter, 31 144145
construction discharge
feeder tanks, 271, 273 relief valves, 369
continuous drawoff, 8486 submerged valves, 117118, 117
control valves see valves see also outlets
correction factors disconnecting valves, 286291, 287291
buckling pressures, 455 disks
corrosion equations, 47 bursting disks, 375, 375
costs lift-disk check valves, 424
protection, 200201 split disk valves, 420421, 420, 441, 442
coupling tilting disk valves, 417418, 437, 439,
pumped outfall pipeline, 355, 356 440
covers, 320321, 464 displacement pumps see reciprocating pumps
crevices, 34 distribution
cross-section changes, 476478, 476 free gas, 5859
cross-section of pipeline, 214 distribution systems, 85, 86, 482483
cyclical oscillations, 472473 DN 450 mains, 179, 292, 293, 314316,
cylindrically balanced relief valves, 371, 372, 315317
373 DN 600 mains, 179
DN 700 mains, 292, 314316, 315317
Daer network link, 501507, 502 doors
damage check valves, 382388
air valves, 353 recoil valves, 416, 417
Damascus irrigation project, 286291, reopening, 382388
287290 swing check valves, 413
damped check valves, 415416 double orifice valves (DOVs), 295, 296
damped relief valves, 362 air flow, 299, 300
damper installations, 518519 locations, 297298, 297
dams, 114 downstream head variations
Darcy equation, 46 check valves, 378, 379
definition sketches pump trips, 484
connection analysis, 475 sleeve valves, 429, 430
cross-section change, 476 downstream pipework
deformation analysis, 459 check valves, 328, 329
sleeve valves, 427 downstream pumping failure
definitions, 12 heads, 406
deflections downstream pumping stations, 193,
flexible pipes, 447, 449451, 465, 467 196199
pipes, 12 downstream sleeve volume, 431
deformation, 448, 452, 458462, 459 drainage pumping stations, 292293
Delft Hydraulics, 436 drawoff, 8486
Demag DRV-B valves, 442443 duckbill valves see sleeve valves
530
Index
531
Pressure transients in water engineering
532
Index
533
Pressure transients in water engineering
534
Index
535
Pressure transients in water engineering
536
Index
537
Pressure transients in water engineering
538
Index
539
Pressure transients in water engineering
540