Fluids Handling

40 www.aiche.org/cep December 2006 CEP

W
ater hammer is generally defined as a pressure
surge or wave caused by the kinetic energy of a
fluid in motion when it is forced to stop or
change direction suddenly, such as in the slow or abrupt
startup or shutdown of a pump system. It can also occur
because of other operating conditions, such as turbine fail-
ure, pipe breakage or electric power interruption to the
pumps motor.
Today, informatics methods are commonly used to per-
form complex water hammer calculations. However, the
lack of information about certain physical flow properties,
as well as the intricacy in handling complicated equations
make this phenomena analysis difficult. The objective of
this paper is to provide a practical and simplified method-
ology to calculate four main phenomenon parameters of
water hammer:
velocity of the pressure wave or celerity
phenomena critical time
maximum head developed in the maximum pressure time
minimum head developed in the critical time.
Equations and basic considerations
One of the most basic equations for the maximum pres-
sure gradient calculation is Joukowskys equation:
h
w max
= (a)(v)/g (1)
where: h
w max
is the maximum fluid elevation head to water
hammer, ft; a is the wave celerity, ft/s; v is the flow veloci-
ty, ft/s; and g is the gravitational acceleration constant, ft/s
2
.
The volumetric flow is given as:
Q = vA (2)
where A is the flow area, ft
2
. A can be defined as:
A = D
i
2
/4 (3)
where D
i
is the internal diameter of the pipe, in. Substituting
Eqs. 2 and 3 into Eq. 1 and solving for a, yields:
where is the pipe thickness, in.
Eq. 4 can then be rewritten as:
where is the density, lb/ft
3
; K is the liquid compressibility
volume factor, lb/in
2
; C
1
is Poissons ratio; and E is the
maximum yield stress, lb/in
2
. These equations were used to
generate Figures 13.
When the valve that stops the fluid flow is closed in a time
slower than the critical time (this reduces the effect of water
hammer), the Allievi equation is used:
Use this graphical method to quickly
and reliably determine the main data
wave celerity, critical time, maximum
head developed in the maximum
pressure time and the minimum head
developed in the critical time
produced by water hammer.
Alejandro Anaya Durand
Mauricio Marquez Lucero
Maria del Carmen Rojas Ocampo
Carlos David Ramos Vilchis
Gonzalez Vargas Maria de Lourdes
National Autonomous University of Mexico
Simplified Analysis of
Water Hammer
a
g h
Q
D
w i
=

4
4
2
2
max
( )
a
g K
D C
E
=

1
1
1

1
2
5 ( )
CEP December 2006 www.aiche.org/cep 41
where: C is a valve constant; h
0
is the head pump, ft; L is
the length, ft; V is the velocity, ft; and t
c
is the time to close
the valve. Substituting Eqs. 7 and 8 into Eq. 6 results in:
Methodology
The simplified graphic
methodology is based on
the simple nomograms
developed for pipeline
transportation of hydrocar-
bons (piping API-5L-X52
of carbonated steel).
Figures 13 plot the wave
celerity on the y-axis and
the diameter divided by the
pipe thickness in the x-axis
for commercial LPG, crude
oil and water, respectively.
Also plotted in Figures 13
are the different values of
density in degrees (API)
and the quotient value,
h
w max

2
/Q. Using Figures 1,
2 or 3, the maximum head
at the instant the valve clos-
es or the pump stops can be
calculated.
Given the values of D/
and the density in API, the
wave celerity, h
w max

2
/Q
and t
c
/L can be determined.
This is done by drawing a
horizontal line from the
intersection to the left side
of the figure to determine
the wave celerity, and to
the right side to determine the quotient value t
c
/L. The
value of h
w max

2
/Q is determined by observing the value
of the vertical curve at which D/ and the density in
API meets. With the knowledge of the value of h
w
max

2
/Q, the maximum pressure with instant valve closing
time can be obtained. This value should be added to the
systems pressure to determine the system overpressure.
The maximum and the minimum head, with different
valves closing can be calculated using Figure 4, which
plots the pump head against the maximum pump head to
water hammer. It also shows a family of curves with dif-
ferent LQ/D
2
t
c
values. Knowing the values
for the pump head and LQ/D
2
t
c
, the maxi-
mum and minimum pressure that the system
can handle can be determined by drawing a
horizontal line from the point that the pump
head and LQ/D
2
t
c
intersects.
a = wave celerity, ft/s
A = flow area, ft
2
C = valve constant as defined by the Allievi Eq.
C
1
= Poissons ratio
D
i
= internal diameter, in.
E = maximum yield stress, lb/in
2
g = gravitational acceleration constant, ft/s
2
h
o
= pump head, ft
h
w
= fluid elevation head to water hammer, ft
K = liquid compressibility volume factor, lb/in.
2
L = length, ft
Q = volumetric flow, ft
3
/s
t
c
= close time valve, s
v = flow velocity, ft/s
V = velocity, ft/s
X = value of the quotient obtained in Figure 4
Greek Letters
= pipe thickness, in.
= density, lb/ft
3
Nomenclature
h
h
C C C
o
= + ( )
( )
2
4 6
2 2
0 5 .
( )
C
LV
gh t
c
=
0
7 ( )
v
Q
D
i
=
4
8
2

( )
h
h
h g
LQ
D t
o
o i c
=

2
1 4
2

22
2
1 4
4
1

+
h g
LQ
D t
o i c
hh g
LQ
D t
o i c

4
2
2

0 5
9
.
( )
I Figure 1. Wave celerity vs. D/ vs. valve close time for commercial liquefied propane gas (LPG; volume
factor = 67,000 lb/in.
2
) in carbon steel pipe.
2,000.0
2,100.0
2,200.0
2,300.0
2,400.0
2,500.0
2,600.0
2,700.0
2,800.0
10.0 15.0 20.0 25.0 30.0 35.0
0
.
7
0
0
.
6
0
0
.
5
0
0
.
4
0
0
.
3
0
0
.
2
5
0
.
2
0
0
.
1
5
h
w

2 /
Q

=

0
.
1
0
10 API
20 API
30 API
40 API
50 API
60 API
70 API
80 API
90 API
1.0x10
-3
9.1x10
-4
8.7x10
-4
8.3x10
-4
8.0x10
-4
7.7x10
-4
7.4x10
-4
7.1x10
-4
9.5x10
-4
Diameter (D) / Thickness ()
W
a
v
e

C
e
l
e
r
i
t
y

(
a
)
,

f
t
/
s
V
a
l
v
e

C
l
o
s
e

T
i
m
e

(
t
c
)

/

L
e
n
g
t
h

(
L
)
,

s
/
f
t
0
.
8
0
42 www.aiche.org/cep December 2006 CEP
Sample calculation
A 12-in. pipe used to
transport hydrocarbons
(commercial LPG with a
density of 51.4 lb/ft
3
) with
40API has a length of 12
km and a flow of 25,000
bbl/d. The pipe is con-
structed from carbonated
steel API-5L-X-52 and has
a constant thickness of
0.406 in. The pumps dis-
charge pressure is 400
lb/in.
2
and the valve section
closing time is 60 s.
Calculate the: maximum
pressure; maximum pres-
sure critical time; maxi-
mum pressure with a valve
closing time of 30 s; and
overpressure.
Calculate the maxi-
mum pressure. First,
determine D/. Start by
calculating the internal
diameter:
D = 12 in. ((2)(0.406 in.))
= 11.18 in. = 0.93 ft
Therefore:
D/ = (11.18 in.)/(0.406 in.)
= 27.53
The maximum system
pressure is:
System pressure =
(400 lb/in.
2
)(144)/(51.4 lb/ft
3
)
= 1,120.6 ft
Calculate the maxi-
mum pressure critical
time. Using Figure 1, start-
ing at D/ = 27.53, draw a
vertical line until it inter-
cepts the 40API line. At
this intersection, the fol-
lowing information is
obtained:
Fluids Handling
3,000.0
3,100.0
3,200.0
3,300.0
3,400.0
3,500.0
3,600.0
3,700.0
3,800.0
10.0 15.0 20.0 25.0 30.0 40.0
6.7x10
-4
6.5x10
-4
6.1x10
-4
5.9x10
-4
5.7x10
-4
5.6x10
-4
5.4x10
-4
5.3x10
-4
6.3x10
-4
Diameter (D) / Thickness ()
W
a
v
e

C
e
l
e
r
i
t
y

(
a
)
,

f
t
/
s
V
a
l
v
e

C
l
o
s
e

T
i
m
e

(
t
c
)

/

L
e
n
g
t
h

(
L
)
,

s
/
f
t
35.0
3,900.0
4,000.0
4,100.0
5.1x10
-4
5.0x10
-4
4.9x10
-4
10 API
20 API
30 API
40 API
50 API
60 API
70 API
80 API
90 API
h
w

2 /
Q

=

0
.
1
0
0
.
7
0
0
.
6
0
0
.
5
0
0
.
4
0
0
.
3
0
0
.
2
5
0
.
2
0
0
.
1
5
0
.
8
0
I Figure 2. Wave celerity vs. D/ vs. valve close time for crude oil (volume factor = 150,000 lb/in.
2
) in
carbon steel pipe.
I Figure 3. Wave celerity vs. D/ vs. valve close time for water (volume factor = 300,000 lb/in.
2
) in
carbon steel pipe.
12.0 17.0 22.0 27.0 32.0
Diameter (D) / Thickness ()
37.0
3,900.0
4,000.0
4,100.0
4,200.0
4,300.0
4,400.0
4,500.0
4,600.0
4,700.0
W
a
v
e

C
e
l
e
r
i
t
y

(
a
)
,

f
t
/
s
4,800.0
4,900.0
5,000.0
5.13x10
-4
5.00x10
-4
4.76x10
-4
4.65x10
-4
4.55x10
-4
4.44x10
-4
4.35x10
-4
4.26x10
-4
4.88x10
-4
V
a
l
v
e

C
l
o
s
e

T
i
m
e

(
t
c
)

/

L
e
n
g
t
h

(
L
)
,

s
/
f
t
4.17x10
-4
4.08x10
-4
4.00x10
-4
5,100.0
5,200.0
5,300.0
5,400.0
5,500.0
5,600.0
5,700.0
3.92x10
-4
3.85x10
-4
3.77x10
-4
3.70x10
-4
3.64x10
-4
3.57x10
-4
3.51x10
-4
0
.
7
0
0
.
6
0
0
.
5
0
0
.
4
0
0
.
3
0
0
.
2
5
0
.
2
0
0
.
8
0
10 API
20 API
30 API
40 API
50 API
60 API
70 API
80 API
90 API
h
w

2
/
Q

=

0
.
1
0
CEP December 2006 www.aiche.org/cep 43
a = 2.389 ft/s
h
w max

2
/Q = 0.124
t
c
/L = 0.00083 s/ft
Solve for h
w max
, given Q =
25,000 bbl/d = 1.78 ft
3
/s and
= 51.4 lb/ft
3
= 0.029 lb/in.
3
:
h
w max
= XQ/
2
=
0.124(1.78ft
3
/s)/(0.406 in./
12 in./ft)
2
= 192.8 ft
The maximum over-
pressure is:
1,120.6 ft + 192.8 ft =
1,313.4 ft
Calculate the maximum
pressure critical time:
t
c
= (t
c
/L)L =
(0.00083 s/ft)(39,370 ft) =
32.6 s
Calculate the maximum pressure with a valve closing
time of 30 s.
LQ/D
i
2
t
c
= (39,370 ft)(1.78 ft
3
/s) / (0.93 ft)
2
(30 s) =
2,700.8 ft
2
Calculate the overpressure and subpressure. Knowing
the maximum pump head value system, as well as LQ/D
i
2
t
c
,
the overpressure can be determined by using Figure 4. The
overpressure is 108.8 ft. Therefore, the total systems over-
pressure is 1,120.6 ft + 108.8 ft = 1,229.4 ft. Similarly,
using Figure 4, the subpressure is 104.5 ft. Thus, the total
systems subpressure is 1,120.6 ft 104.5 ft = 1,016.0 ft.
ALEJANDRO ANAYA DURAND is professor of chemical engineering at the
National Autonomous University of Mexico (UNAM; Parque Espaa 15B
Colonia Condesa, Mxico D.F. C.P. 06140, Mxico; Phone and fax:
52110385; E-mail: [email protected]). He has 40 years of
experience as an educator, teaching subjects such as heat transfer,
fluid flow and project engineering. After 30 years of holding top-level
positions in project engineering, Durand retired from the Instituto
Mexicano del Petroleo in 1998. He is also a consultant of several
engineering companies. Durand is a Fellow of AIChE and holds an M.S.
in chemical engineering from UNAM. He has published more than 200
technical articles in local and international magazines, related to
chemical engineering and education.
MAURICIO MARQUEZ LUCERO Is project manager for the Delta Project and
Development Co. (Sur 73 No 311 bis 1 Colonia Sinatel, Del Iztapalapa,
Mxico D.F., C.P. 09470; Phone: 56721237; Fax: 55392883; E-mail:
[email protected]). He is a member of the Instituto Mexicano
de Ingenieros Qumicos. Lucero has published several technical articles
in local and international magazines, in chemical engineering.
MARIA DEL CARMEN ROJAS OCAMPO is a chemical engineering student at
the National University Autonomous of Mexico (Sur 73 No 311 bis 1
Colonia Sinatel, Del Iztapalapa, Mxico D.F., C.P. 09470; Phone:
56721237; Fax: 55392883; E-mail: [email protected]).
She has published several technical articles on local and international
magazines, in chemical engineering.
CARLOS DAVID RAMOS VILCHIS is a chemical engineering student at the
National University Autonomous of Mexico (Cerrada Canal Nacional #6
col. Tejomulco Nativitas C.P. 16510 Delegacion Xochimilco Mexico. D.F.;
Phone: 21578360; E-mail: [email protected]). Presently, he
collaborates in a company dedicated to training groups for work in the
chemical industry (pilot plant).
GONZALEZ VARGAS MARIA DE LOURDES is a chemical engineering
student at the National University Autonomous of Mexico (Mxico
Tacuba 1523 edificio Managua departamento 505 colonia Argentina
Poniente C.P.: 11230 Delegacin Miguel Hidalgo Mxico, D.F.; Phone:
55761849; Fax: 55761855).
CEP
I Figure 4. Maximum fluid elevation head to water hammer.
-300
-200
-100
100
300
100 1,100 2,100 3,100 4,100
Fluid Elevation Head, ft
M
a
x
i
m
u
m

E
l
e
v
a
t
i
o
n

H
e
a
d

t
o

W
a
t
e
r

H
a
m
m
e
r
,

f
t
200
0
5,100
S
u
b
p
r
e
s
s
u
r
e
P
r
e
s
s
u
r
e
LQ/D
i
2
T
c
(max) = 10,000
8,000
6,000
4,000
2,000
2,000
4,000
6,000
8,000
10,000
Literature Cited
1. Perry, J. H., Chemical Engineers Handbook, Mc Graw Hill,
4th edition, pp. 644.
2. UNAM, Selection and Operation of Pump System, National
Autonomous University of Mexico (March 1983).
3. American Society of Civil Engineers, Pipe Line Design for
Water and Wastewater, Pipeline Committee on Pipeline Planning
(1975).
4. Douglas, J. F., Solution of Problems in Fluid Mechanics,
Pitman Paperbacks, pp. 109117 (1967).
5. Potter, M. C. and D. C. Wiggert, Mechanics of Fluids,
Prentice Hall, 2nd edition, p. 18 (1997).