Strojniški vestnik - Journal of Mechanical Engineering 62(2016)10, 614-622 © 2016 Journal of Mechanical Engineering. All rights reserved.
D0l:10.5545/sv-jme.2016.3535	Original Scientific Paper
Received for review: 2016-02-29 Received revised form: 2016-06-28 Accepted for publication: 2016-07-04
Numerical Simulation of Transient Pressure Control in a Pumped Water Supply System Using an Improved Bypass Pipe
Xiaozhou Li- Manlin Zhu* - Jiancang Xie
Xi'an University of Technology, State Key Laboratory Base of Eco-hydraulic Engineering in Arid Area, China
The use of a bypass pipe in a pump system is an economical and simple method of preventing water hammer surges to control the extreme transient pressures induced by pump failure. Currently conventional bypass pipes (CBPs) have certain limitations in engineering applications. To overcome these limitations, in this study, a CBP is modified to create an improved bypass pipe (IBP). Unlike that of a CBP, the control valve in an IBP system is a hydraulic control valve that uses a controller on an uninterruptible power source (UPS) to turn the valve on or off. This modification enables the hydraulic control valve to be precisely controlled and facilitates the prevention of water hammer surges. To study the transient control effect of the IBP using the method of characteristics (MOC), a mathematical model of a complex system that includes an IBP, a check valve, and a water pump is established. Using MATLAB software, a water supply system in an industrial zone is used as an example for numerical simulations of the transient process experienced by this pump after it fails. The extreme pressure variations and the envelope of the pipe pressure in a typical section of the pumped water supply system are analysed and compared. The results show that an IBP can effectively reduce the maximum extreme pressure in a pumped water supply system.
Keywords: improved bypass pipe, hydraulic transient, pumped water supply system, UPS controller, numerical simulation Highlights
•	Based on a conventional bypass pipe, a new, economical, and effective water hammer prevention measure, i.e., an improved bypass pipe, is proposed.
•	A hydraulic control valve with an on/off status determined by an uninterruptible power source controller is installed in the bypass pipe to control the transient pressure.
•	The on/off time for the newly developed improved bypass pipe is accurately controlled to facilitate transient pressure control.
•	Based on the method of characteristics, a mathematical model for a complex system containing an improved bypass pipe is established.
0 BACKGROUND
Pump failures induce rapid changes in the flow rate, which causes extremely high-pressure water hammer surges in the water supply system, including extremely high- and low-pressure water hammer surges. Extreme water hammer pressure damages pipes and hydraulic equipment and results in water supply service interruptions. To eliminate or reduce the impact of extreme water hammer pressure on water delivery systems, water hammer prevention measures are normally implemented; these include bypass pipes, air valves, check valves, air tanks, pressure regulating chambers, unidirectional pressure regulating chambers, pressure discharge valves, water hammer prevention valves, etc. [1] to [3]. An air valve is a common, economical, and effective piece of water hammer prevention equipment that is normally installed at the local peak of a pipe. When negative pressure occurs inside the pipe, the air valve opens to allow air in; when the pressure inside the pipe increases, the air valve releases it by discharging air [4] and [5]. A check valve is an economical and
reasonable water hammer prevention measure that is widely used in pumping stations [6]. Check valves can prevent water pumps from rotating in reverse due to liquid backflow and thus prevent damage to the motor. However, instant check valve closure results in a catastrophically high pressure [7] and [8]. Therefore, when a check valve is installed in a pumping station to control transients, an air tank should also be installed near the pumping station to prevent the increase in pressure that results from closing the check valve [9]. However, air tanks have drawbacks that include their large volume, high maintenance cost and complex operation (the gas in the tank should be replenished). A pressure discharge valve with the proper parameters can prevent a water delivery system from generating excessively high or low pressures; otherwise, the negative impact of water hammer surges on the water delivery system is exacerbated [10]. To ensure proper protection of a water delivery system by a pressure discharge valve, the pressure discharge valve should have so little inertia that it can be opened promptly in response to a rapid change in the pressure to prevent delays in its opening [11]. A water hammer
614
*Corr. Author's Address: Xi'an University of Technology, NO.5 South Jinhua Road, Xi'an 710048, Shaanxi ,China, zhuml@xaut.edu.cn
Strajniski vestnik - Journal of Mechanical Engineering 62(2016)10, 614-622
prevention valve may have inaccurate hydraulic control that causes its opening to be delayed, or that causes it to close excessively fast, which results in water hammer damage. Wylie [12] proposed two types of conventional bypass pipes (CBPs) for preventing water hammer surges in pumping stations. One was to install a check valve in the water pump's outlet and a control valve in the bypass pipe to prevent the pump from rotating in reverse due to water backflow. After pump failure, the bypass pipe control valve opens instantly and then closes gradually to release the high pressure inside the system and eliminate water hammer surges. The other type was designed to prevent low pressure and water column separation in the water pump's outlet from being caused by pump failure. A control valve was installed at the water pump's outlet, and a check valve was installed in the bypass pipe; therefore, the water flowing into the reservoir could be directed to the water pump's outlet through the bypass pipe. Both types of CBPs have certain drawbacks. The first type of bypass pipe cannot accurately control the valve's on/off timing; therefore, it may not be able to prevent water hammer surges as well as is desired. The second type of bypass pipe cannot release the high pressure in the system and is only applicable to pumping station water supply systems that ignore the effect of backflow and have water pump inlets with positive water head pressures.
Compared with experimental studies, numerical simulations take less time and are more economical. With advances in computer technology, numerical simulations have become widely employed in water hammer analyses [13] to [15]. To predict pressure variations in a water supply system and to choose appropriate water hammer preventive measures, numerical simulations of water hammer surges are indispensable, simple, and effective [16]. Methods of numerically analysing water hammer surges include arithmetic, graphical, algebraic and linear analysis methods as well as the method of characteristics (MOC). At present, the MOC is the simplest and most popular method for numerically solving transient flow problems. Izquierdo proposed mathematical models for simulating transient boundary conditions in simple and complex hydraulic systems [17] and [18]. Because the results of hydraulic transient numerical simulations are in good agreement with those of test measurements [19] and [20], the results of numerical simulations of water hammer surges are widely accepted.
A composite bypass pipe and check valve is installed in a pumping station to reduce the extremely high water pressure associated with water hammer
surges in a pumped water supply system. However, this type of CBP has limitations. Therefore, this study has three purposes:
1)	To overcome the limitations of CBPs, an improved bypass pipe (IBP)-based water hammer prevention measure is proposed, and the operational principle of an IBP is explained.
2)	Based on previous studies and the method of characteristics (MOC), a mathematical model of the boundary conditions of an IBP is established, and a solution for this mathematical model is described.
3)	A water supply system in an industrial zone is used as an example in a numerical simulation of the water hammer prevention effect of a pump with an IBP-based protective device on a water supply system and to provide evidence for the use of an IBP-based water hammer prevention measure in a system.
1 THE STRUCTURE AND OPERATIONAL PRINCIPLE OF THE IMPROVED BYPASS PIPE
Fig. 1 shows a diagram of the IBP. It is based on a CBP, but the control valve in the CBP is replaced by a hydraulic control valve. The hydraulic control valve consists of a valve body, a hydraulic pressure system, an oil cylinder, an energy storage tank, and an electromagnetic valve. The IBP's hydraulic control valve leverages an uninterruptible power source (UPS), which provides power to the electromagnetic valve, which controls when and for how long the main valve is open. This improvement can prevent the valve from turning on or off in response to man-made or hydraulic interference and provide more accurate control over the valve's status than a conventional manually or hydraulically controlled valve can.
Fig. 1. Schematic of the IBP
After the pump fails, the UPS-powered controller can predefine a delayed opening time, t1 (the interval between when the pump fails and when the electromagnetic valve powers on), for the electromagnetic valve according to the requirements. After the electromagnetic valve powers on and opens, oil from the energy storage tank enters the
Numerical Simulation of Transient Pressure Control in a Pumped Water Supply System Using an Improved Bypass Pipe
615
Strojniski vestnik - Journal of Mechanical Engineering 62(2016)10, 614-622
rod-less chamber of the oil cylinder and rotates the valve board to open the main valve. The main valve's opening time, t2 (the time it takes the main valve to make the transition from being completely closed to being completely open), is regulated via the pressure at the valve that controls the oil's flow rate according to actual requirements. After the time allotted for the electromagnetic valve to be open expires, the UPS-powered controller closes the electromagnetic valve, and the oil in the energy storage tank enters the rod chamber of the oil cylinder and pushes the piston back to close the main valve. The main valve's closing time, t4 (the time it takes for the main valve to make the transition from completely open to completely closed), is regulated by a valve closure regulation valve according to actual requirements. The difference between the time for which the electromagnetic valve is open and the time required for the main valve to open is defined as the main valve's opening duration, t3 (the time the main valve is completely open). Fig. 2 shows the IBP in operation.
Fig. 2. Flowchart of the IBP's operation
2 CHARACTERISTIC EQUATIONS FOR CALCULATING WATER HAMMER SURGES AND THEIR SOLUTION
equations comprise a pair of non-linear hyperbolic partial differential equations, which is very difficult to solve using conventional methods. Therefore, the MOC is used to convert the two partial differential equations that take pipeline hydraulic friction into account into ordinary differential equations of a specific form. These are called characteristic equations and are given in Eqs. (3) and (4) for C+ and in Eqs. (5) and (6) for C-:
g dh a dt
dv dt
+ — + = 0,
M*
2D
dx dt
_ gdh + dv + fvJH = 0 a dt dt 2D
dx
— = -a, dt
(3)
(4)
(5)
(6)
where, because v sinff is typically much smaller compared with any other item in the equation for the transient process and because wave speed a is far larger than the flow velocity v, v sinff is ignored in Eq. (3) and Eq. (5), and the flow velocity, v, is ignored in Eq. (4) and Eq. (6).
Next, the characteristic equations are integrated to generate a set of simplified finite-difference equations, Eqs. (7) to (10).
C + : hiJ&l = Cp - BQt
i, jAt '
C- : hiJ&l = CM + BQi
i, jAt'
(7)
(8)
CP hi-l,(j-1)At + BQi-1,(j-1)At RQi-1,(j-1)At Qi-1,(j-1)At ' (9)
2.1 Basic Differential Equations for Water Hammer Surges
CM KM j-i)At - BQi+i
,( j-1)A - RQi+i ,(j-1) At Q+i	(10)
The basic differential equations for water hammer surges consist of an equation of motion Eq. (1), and a continuity equation, Eq. (2).
dh 1 dv v dv f vlvl
— +--+--+ ^—i-1^ 0, (1)
dx g dt g dx D 2g
dh dh „ a2 dv — + v(--sm0) +--= 0.	(2)
dt dx	g dx
2.2 Simplified Finite-Difference Equations
The numerical calculation of water hammer surges is based on comprehensive and complete basic differential equations, Eqs. (1) and (2). These two
In Eqs. (7) to (10) B = a/(gA), R =fAx/(2gDA2),
hf-iti-i)^ bin,^)^ Qi-i,(j-\)At and Gin/j^A^ and are known from the previous step; only hij-At and QijAt are unknown.
2.3 Procedure for Solving the Water Hammer Equations
Fig. 3 shows the grid diagram used in the MOC. The solution process normally starts with a stable flow at t=0. Therefore, the initial values of h and Q in each calculation section are known. Based on their values at t=0 in all the sections of a layer, the values of h and Q at t=At in all the sections of the layer are calculated. The boundary section is determined by the boundary condition. Next, the values at t = 2At in all the sections
616
Li, X. - Zhu,	M.- Xie, J.
Strajniski vestnik - Journal of Mechanical Engineering 62(2016)10, 614-622
of the layer are calculated. The same procedure is followed until values at the required time, t =jAt, are calculated.
Fig. 3. A characteristic line in space (x) - time (t)
3 BOUNDARY CONDITIONS AND A MATHEMATICAL MODEL FOR THE IBP
Fig. 4 shows that the check valve and the controllable bypass pipe composite are installed at the pump's outlet to prevent water hammer surges. The check valve is installed at the water pump's outlet to prevent backflow in the water delivery pipe due to accidental power failure and to facilitate maintenance and repairs. However, a check valve by itself may close suddenly and cause damage if there is a water hammer surge resulting from a pump failure or if there is a severe engineering accident. Therefore, installing a controllable bypass pipe can prevent an increase in pressure resulting from a sudden closure of the check valve. To study the water hammer prevention effect of a controllable bypass pipe, the MOC is employed to perform a numerical analysis and calculations. The MOC is a widely used numerical method that is verified in reference [12].
Before establishing a mathematical model, the following assumptions are made:
1.	From when the pump fails to when the flow rate at the water pump's outlet drops to zero, the check valve local coefficient of resistant is a constant, kc1, and when the flow rate at the water pump's outlet is Q6 < 0, the check valve closes instantly.
2.	The local water head losses at T-type nodes upstream and downstream of the water pump are ignored.
3.	The T-type nodes connecting the pipe and the bypass pipe upstream and downstream of the water pump (the dashed box in Fig. 4) are very short, and loss during the flow is ignored; they are analysed as an integral unit.
When a power failure occurs, the water pump in the water supply system suddenly loses power, the pump's outlet pressure drops, and water flows back. When the flow rate at the water pump's outlet is less than or equal to zero (Q6,(j■_1)&t = 0), the check valve closes instantly (Q6, ;At = 0). When the pump fails at time T < t1, the IBP closes (Q4, jM = Q 5, jM = 0). When the pump fails at time t1 < T < /j+t2+t3+t4 , the IBP opens quickly and stays open. At this point, the backflow water column is discharged from the IBP to the reservoir. When the pump fails at time T> 11+t2+t3+t4 , the IBP closes slowly (04, ]A1 = Q 5, ]At = 0).
Fig. 4. A boundary model of the IBP Continuity equation:
Q2, jAt = Q3, jAt + Q4, jAt ' Q7, jAt = Q5, jAt + Q6, jAt, Q3, jAt = Q6, jAt = UQR ' Q4, jAt = Q5, jAt,
jAt	jAt
water balance equation:
hjA, + hr ( )
a0 (n + tan 1 —) + b0 a
Q3,jAt Q3,jAt 1 ,
-kc-i—2—1 = h,
1 2gA2
'6, jAt'
Q4,jAt Q4,jAt ,
h — kc ---- = h
4, jAt /vl2 0 .2	5, jAt
2 gA2
h2,jAt = h3,jAt =h4,jAt, h5,jAt = h6,jAt =h7,jAP (12)
characteristic equation:
h2,jAt CP B1-2Q2,jAt, ^7, jAt CM + B7-8Q7,jAt,
(13)
water pump's unit inertia equation:
(a2 +u2)
«j (n + tan 1 —) + \ a
-P0 -
WRA Nr n , , A --R-(a0-a) = 0.
g Tr 15 A? Eqs. (11) to (14) are consolidated into
(14)
Numerical Simulation of Transient Pressure Control in a Pumped Water Supply System Using an Improved Bypass Pipe
617
Strojniski vestnik - Journal of Mechanical Engineering 62(2016)10, 614-622
F1 = Cp - CM - (( + B7_s){QRv + QaJM ) + +Hr (a2 +u2 )
a0 (n + tan —) + b0 0a0
-kc, Qr , = 0, 1 2 gA '
(15)
F2 = Cp - CM - (( + B7_s )( + Qj ) -
Q4,jAt Q4,jAt „
-kc2-!-7—1 = 0,
2 2g^2
(16)
F 3 = (a2 +u2 )
WR2 NR n
al (n + tan 1 —) + bl a
+P0---R-(a0-a) = 0.
0 g Tr 15A? 0
(17)
Eq. (15), Eq. (16) and (17) can be solved using the Newton-Raphson method [12].
dF1
F1„= — = -( B-2 + B-zQ + du
+Hr \2u
a01 n + tan 1 — | + b0
+ a0a kc1 Qu ,(18)
F1„=dF1 = H\2a da
a0 \n+ tan1 — | + b0 - a0v\, (19)
F1

dF1
"da
= -( B1-2 + B7-S)'
4,jAt
dF 2
F 2U = — =-( BI-2 + B7_S)QR , du
F 2..ÊFÎ = 0,
da
(20)
(21) (22)
F 2r
dF 2
- = -( Bx_2 + B7_g) - kc2 Qjr, (23)
4, jAt
T73 9F3 2
F 3U =-= 2u
du
al | n + tan 1 — | + bl
-ala, (24)
w dF3 i F 3a =-= 2a
da
al \k+ tan 1 — | + bl
F 3r
WR Nr _ g Tr 15 At '
dF 3- = 0.
(25)
(26)
dQj
The Newton-Raphson equations are represented
as matrices:
(
Flu	Fla	Flo Q4, jAt
F 2V	F 2a	F 2o Q4, jAt
F 3u	F 3a	F 3o °4, jAt
(Au ^		(-Fl 1
Aa	=	- F 2
Uâ,jAt		F 3,
(27)
The solutions of the above equations are represented as matrices,
^Au Ï Aa 1
AQj
(F1 F2 F3Q + F1 F2Q F3
v u a Qa.ja,	a QaJA, u
+F1Qaja,F2uF3a- F1Qija,F2a F3u -F1uF2Q4ja,F3a- F1a F2uF3QtJa, )
(F2aF3Q - F2q F3a F1Q F3a- F1a F3Q
a	Q4,jAt	Q4,jAt	a	Q4,jAt	a	a	Q4 ,jAt
F2q F3u - F2uF3Q
Q4, jAt	u	u	Q4,jAt
v F 2u F 3a- F 2a F 3u
F1 F2Q -F1Q F2 *)(-F1\
a Qaja,	QaJA, a F 1
F1 F3n - F1n F3 F1n F2 - FI F2
U Qa F1a F3u- Flu F3a
u Qa, j
Flu F2a- Fla F2u ,
-F2 -F3
(28)
2
Initial values are assigned to u and a as follows:
U = 2U(j_i)At — U(j—2)At, a = j-1)At — a(j-2)At,
Q4,jAt = 2Q4,(j-1)At — Q4,( j-2)At.	(29)
Eqs. (15) to (26) are substituted into Eq. (28) and solved for Aw, Aa and AQ4,jM . Then, Eq. (29) is rewritten as:
618	Li, X. - Zhu,
v=v + ^v, a=a+^a, Q4JAt = Q4JAt + Aq4ja,. (30)
The above procedure is repeated until Aw, Aa and AQ4 ;Ai are within the allowable deviation range:
|Ao| + |Aa| + |AQ4, ;Ai| < e.
The allowable deviation e is set to approximately 0.0002 [12].
M.- Xie, J.
Strajniski vestnik - Journal of Mechanical Engineering 62(2016)10, 614-622
4 USE OF AN IBP IN A PUMPED WATER SUPPLY SYSTEM
4.1 System Overview
coefficient^ , and its local resistance coefficient, kc, is shown in Eq. (31).
kc = 1/C2 .	(31)
As shown in Fig. 5, a pumping station water supply system consists of three identical water pumps operating in parallel, two parallel pipes, and upstream and downstream reservoirs. The pumping station's static head is 142 m, the water pump's rated head is Hr = 149.5 m, the rated flow rate is QR = 0.37 m3/s and the water pump's rated rotation speed is Nr = 1450 rpm. The pump's characteristic curve is based on data for the specific rotation speed ns=25 (SI), according to reference [12]. Each pump unit's extreme moment of inertia is WR2 = 450Nm2. The pipe length is L = 1390 m, and the pipe's inner diameter is D = 781 mm. The wave speed is a = 1213 m/s, and the pipe's coefficient of friction is f= 0.02. When the pump stops and backflow starts, the check valve behind the pump closes immediately, and the bypass pipe control valve opens or closes according to the predefined procedure. According to Appendix B of reference [12], when the check valve is completely open (i.e., the valve's opening level is t = 1), its discharge coefficient is Cd = 1.4, and when the check valve is completely closed (i.e., the valve's opening level is t = 0), its discharge coefficient is Cd = 0. Fig. 6 shows the relationship between the hydraulic control valve's discharge coefficient and opening level. The relationship between the valve's discharge
Pump Check valve
Upstream reservoir	-U-w-i Ott- -o-
	I \ BP
Pipeline
Downstream reservoir
Pipeline
Fig. 5. A diagram of a pumped water supply system
1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 Opening ratio
Fig. 6. The hydraulic valve's discharge coefficient
4.2 Numerical Simulations
Numerical simulations and calculations for two postpump failure scenarios, without an IBP and with an IBP, are performed for the pumped water supply system shown in Fig. 5.
4.2.1 Case 1: Without an IBP
When all three parallel water pumps fail due to an accident and backflow occurs, the check valves behind the pumps immediately close, and the IBP is not involved. The transient process that occurs in the pumped water supply system in this situation is numerically simulated. Fig. 7 shows the flow rate at each water pump's outlet (Q6) and the typical section pressure (h6).
Fig. 7. Variations in the flow rate and pressure without an IBP
Fig. 7 shows that when the three parallel water pumps fail simultaneously after an accident, the flow rate at each water pump's outlet (the pipe between the water pump and the main pipe) drops rapidly until T = 1.99 s, when it becomes negative (i.e., reaches its minimum value), and the typical section pressure at the water head also drops to its minimum value. At this point, the check valve closes instantly to prevent a high volume of backflow, and the flow rate at each water pump's outlet stabilizes at zero. However, the closure of the check valve also stops the transient high-pressure wave from returning and being discharged, which may lead to an extremely high-pressure water hammer surge in the system.
4.2.2 Case 2: With an IBP
Fig. 8 shows that the hydraulic control valve's on/ off program opens the IBP 2.5 s after pump failure (t1 = 2.5 s). It takes 3 s (t2 = 3.0 s) to make the transition
Numerical Simulation of Transient Pressure Control in a Pumped Water Supply System Using an Improved Bypass Pipe
619
Strojniski vestnik - Journal of Mechanical Engineering 62(2016)10, 614-622
from being completely closed to being completely open; it remains completely open for 5 s (t3 = 5 s), and then it takes 15 s (t4 = 15 s) to make the transition from being completely open to being completely closed.
25 Time [s]
Fig. 8. Operation of the hydraulic valve
Fig. 9 shows the flow rate at the outlet of a water pump with an IBP and typical variations in the section pressure. Fig. 9 shows that when the three parallel water pumps fail simultaneously, the flow rate at each water pump's outlet (Q6) also drops rapidly until T = 1.99 s, when it becomes negative (i.e., reaches its minimum value), and the typical section pressure at the water head drops to its minimum value as well. At this point, the check valve closes instantly to prevent a high volume of backflow, and the flow rate at each water pump's outlet stabilizes at zero. Because the IBP opens at a predefined time, the returned transient high-pressure wave is not blocked by the closed check valve. Instead, it is discharged by the IBP, eliminating the extremely high-pressure water hammer surge from the system. During this period, the maximum flow rate in the IBP reaches Q5 = -1.062 m3/s. In addition, the IBP gradually closes at a predefined time to prevent valve closure-induced water hammer surges.
10 15 20 25 30 35 40 45 50 Time [s]
Fig. 9. The flow rate and pressure variation at the IBP
4.3 Comparison and Analysis of Simulation Results
The mathematical model of the IBP presented in this study is used in numerical simulations and
calculations of a transient process in a pumped water supply system. The typical variations in the section pressure determined by the simulations are shown in Fig. 10. In the graph, the red dashed line represents typical variations in the section pressure when there is no IBP, and the solid blue line represents typical variations in the section pressure when there is an IBP. The simulation results show that without an IBP, the water head of the maximum transient pressure in a typical section reaches 223.2 m, which is 1.51 times the initial pressure water head (147.5 m). With an IBP, the water head of the maximum transient pressure in a typical section drops to 171.8 m, which is 1.16 times the initial pressure water head (147.5 m). These results show that the presence of an IBP can significantly reduce the water head of a typical section's maximum transient pressure.
0 5 10 15 20 25 30 35 40 45 50 Time [s]
Fig. 10. Pressure variations in a typical section
600 800 Distance [m]
Fig. 11. Envelopes of extreme pressure
Fig. 11 shows the envelope of the extreme pipe pressure during the transient process in a pumped water supply system. In the graph, the blue dashed-dotted line represents the water supply pipe's central line, the blue solid line represents the water supply system's initial hydraulic grade line, the red dotted line represents the envelope of the extreme pressure during the transient process in a water supply system without an IBP, and the green dashed line represents the envelope of the extreme pressure during the transient process in a water supply system with an IBP.
620
Li, X. - Zhu,	M.- Xie, J.
Strajniski vestnik - Journal of Mechanical Engineering 62(2016)10, 614-622
This graph shows that an IBP significantly reduces the maximum extreme pressure in a pumped water supply system, which ensures that the water supply system operates safely and reliably.
5 CONCLUSIONS AND FUTURE WORK
In a pumped water supply system, pump failure reduces the pressure in the water supply pipe and results in subsequent liquid backflow in the pipe. An enormous backflow not only results in water column separation but also impacts the water pump's impeller and makes it rotate in reverse. In severe cases, it threatens the pumped water supply system's safe operation. Check valves are widely used in pumped water supply systems to prevent pump failure-induced water backflow and reverse rotation of the water pump. However, rapid check valve closure also results in extremely high-pressure water hammer surges. Previous researchers combined a CBP with a check valve to release the high pressure in the system and eliminate water hammer surges. However, because the timing of the CBP valve's on/off setting cannot be accurately controlled, it may not be possible to achieve the desired water hammer prevention effect. To overcome this limitation of CBPs, this study improved the CBP and proposed a new water hammer prevention measure, the IBP. An IBP is a bypass pipe with a hydraulic control valve controller whose on/ off procedure is predefined by a UPS. Therefore, the valve's setting can be accurately controlled. Using the MOC, a mathematical model of the IBP is established, and calculations are performed. This mathematical model simulates and calculates typical variations in the transient section pressure and the distribution of the pressure in the water supply pipe for a pumped water supply system in an industrial zone. A comparison of the simulation results shows that the IBP (whose hydraulic control valve's setting is changed at t1 = 2.5 s, t2 = 3.0 s, t3 = 5 s and t4 = 15 s) water hammer prevention measure in a pumped water supply system can significantly reduce the maximum extreme pressure of water hammer surges in the pumped water supply system. However, some of the problems encountered in this study require further research: 1) The mathematical model of the IBP does not take into account the bypass pipe's length and diameter; it only considers the bypass pipe valve's local coefficient of resistance, so the next step is to incorporate the bypass pipe's length, l, and diameter, d, into the mathematical model to study the effect of the bypass pipe's length and diameter on water hammer prevention. 2) The numerical simulation and calculation of the pumped
water supply system's transient process shows that an IBP can protect the water supply system from damage due to water hammer surges; however, the sensitivity of key parameters affecting the IBP's water hammer prevention effect (the bypass pipe's length, l, and diameter, d, and the bypass pipe hydraulic control valve's delayed opening time, t1, opening time, t2, opening duration, t3 and closing time, t4) will analyzed to provide reference parameters for engineering applications.
6 ACKNOWLEDGMENTS
This research program is sponsored by the National Natural Science Foundation of China (Grant Nos. 41471451 and 51479160).
7 NOMENCLATURE
t T h
v g
a x 0
f
D
A
C+
C-
p
C cm
B R Ax At
i j
N
A:
A2
kc1
kC2
Qr Hr
V
a
ao, b{) a1 b1
fi '
WR2
time [s] time [s]
pressure head [m]
flow velocity [m/s]
acceleration of gravity [m/s2]
wave speed of water hammer surge [m/s]
distance along pipe [m]
pipe slope
Darcy-Weisbach friction factor pipe inner diameter [m] area of pipe section [m2] name for forward characteristic line name for reverse characteristic line known constant in MOC known constant in MOC known constant in MOC known constant in MOC spatial step size [m] time step size [s]
subscript denoting the position of the section
number of time segments
number of pipe segments
area of check valve section [m2]
area of hydraulic valve section [m2]
minor loss coefficient of check valve
minor loss coefficient of hydraulic valve
rated discharge of pump [m3/s]
rated head of pump [m]
dimensionless pump discharge
dimensionless pump speed ratio
parameters for pump head characteristics
parameters for pump torque characteristics
dimensionless torque ratio
combined polar moment of inertia of pump
Numerical Simulation of Transient Pressure Control in a Pumped Water Supply System Using an Improved Bypass Pipe
621
Strojniski vestnik - Journal of Mechanical Engineering 62(2016)10, 614-622
Nr	rated rotation speed of pump [rpm]
Tr	rated torque of pump [Nm]
n	constant
F	function symbol
e	assigned precision control
t	valve opening ratio
L	pipe length [m]
ns	specific speed
Cd	discharge coefficient of valve
l	bypass pipe length [m]
d	bypass pipe inner diameter [m]
8 REFERENCES
[1]	Boulos, P.F., Karney, B.W., Wood, D.J., Lingireddy, S. (2005). Hydraulic transient guidelines for protecting water distribution systems. Journal American Water Works Association, vol. 97, no. 5, p. 111-124.
[2]	Bergant, A., Kruisbrink, A., Arregui, F. (2012). Dynamic behaviour of air valves in a large-scale pipeline apparatus. Strojniški vestnik - Journal of Mechanical Engineering, vol. 58, no. 4, p. 225-237, DOI:10.5545/sv-jme.2011.032.
[3]	Ramezani, L. (2015). An exploration of transient protection of pressurized pipelines using air valves, PhD thesis, University of Toronto, Toronto.
[4]	Karadžič, U., Bulatovic, V., Bergant, A. (2014). Valve-induced water hammer and column separation in a pipeline apparatus. Strojniški vestnik - Journal of Mechanical Engineering, vol. 60, no. 11, p. 742-754, DOI:10.5545/sv-jme.2014.1882.
[5]	Carlos, M., Arregui, F.J., Cabrera, E., Palau, C.V. (2011). Understanding air release through air valves. Journal of Hydraulic Engineering, vol. 137, no. 4, p. 461-469, DOI:10.1061/(ASCE)HY.1943-7900.0000324.
[6]	Thorley, A.R.D. (2004). Fluid Transients in Pipeline Systems: A Guide to the Control and Suppression of Fluid Transients in Liquids in Closed Conduits. ASME Press, New York.
[7]	Sugiyama, H., Mizobata, K., Ohtani, K., Ohishi, T., Sasaki, Y., Musha, H., Miwa, T. (2003). Dynamic characteristics of a check valve for drain pumps. ASME/JSME 4 Joint Fluids Summer Engineering Conference, p. 2851-2856, DOI:10.1115/ FEDSM2003-45256.
[8]	Thorley, A.R.D. (1989). Check valve behavior under transient flow conditions: A state-of-the-art review.
Journal of Fluids Engineering, vol. 111, no. 2, p. 178-183, DOI:10.1115/1.3243620.
[9]	Purcell, P. (1997). Case study of check-valve slam In rising main protected by air vessel. Journal of Hydraulic Engineering, vol. 123, no. 12, p. 1166-1168, DOI:10.1061/(ASCE)0733-9429(1997)123:12(1166).
[10]	Zhang, K.Q., Karney, B.W., McPherson, D.L. (2008). Pressure-relief valve selection and transient pressure control. Journal American Water Works Association, vol. 100, no. 8, p. 62-69.
[11]	Chaudhry, M.H. (1987). Applied Hydraulic Transients, 2nd ed. Van Nostrand Reinhold, New York.
[12]	Wylie, E.B., Streeter, V.L., Suo, L. (1993). Fluid Transients in Systems. Prentice Hall, New Jersey.
[13]	Kaliatka, A., Vaisnoras, M., Valincius, M. (2014). Modelling of valve Induced water hammer phenomeno In a district heating system. Computers&Fluids, vol. 94, p. 30-36, DOI:10.1016/j. compfluid.2014.01.035.
[14]	Skulovich, O., Perelman, L., Ostfeld, A. (2014). Bi-level optimization of closed surge tanks placement and sizing in water distribution system subjected to transient events. Procedia Engineering, vol. 89, p. 1329-1335, DOI:10.1016/j. proeng.2014.11.449.
[15]	Wan, W., Huang, W., Li, C. (2014). Sensitivity analysis for the resistance on the performance of a pressure vessel for water hammer protection. Journal of Pressure Vessel Technology, vol. 136, no. 1, p. 011303, DOI:10.1115/1.4025829.
[16]	Wood, D.J. (2005). Waterhammer analysis-essential and easy (and efficient). Journal of Environmental Engineering, vol. 131, no. 8, p. 1123-1131, DOI:10.1061/(ASCE)0733-9372(2005)131:8(1123).
[17]	Izquierdo, J., Iglesias, P.L. (2004). Mathematical modelling of hydraulic transients in complex systems. Mathematical and Computer Modelling, vol. 39, no. 4-5, p. 529-540, DOI:10.1016/S0895-7177(04)90524-9.
[18]	Izquierdo, J., Iglesias, P.L. (2002). Mathematical modelling of hydraulic transients in simple systems. Mathematical and Computer Modelling, vol. 35, no. 7-8, p. 801-812, DOI:10.1016/S0895-7177(02)00051-1.
[19]	Kendir, T.E., Ozdamar, A. (2013). Numerical and experimental investigation of optimum surge tank forms in hydroelectric power plants. Renewable Energy, vol. 60, p. 323-331, DOI:10.1016/j.renene.2013.05.016.
[20]	Kwon, H., Lee, J. (2008). Computer and experimental models of transient flow in a pipe involving backflow preventers. Journal of Hydraulic Engineering, vol. 134, no. 4, p. 426-434, DOI:10.1061/(ASCE)0733-9429(2008)134:4(426).
622
Li, X. - Zhu,	M.- Xie, J.