The focus of this work is on deriving the governing equations for the ‘shock damper’, an instrument used in water flow systems. The shock damper is a mechanical tool with a single degree of freedom system that includes a tank, connecting pipe, mass, spring and damper. This tool serves as a boundary condition for characteristic line equations. Additional equations are included by considering the conservation of mass, momentum and energy. This system of equations is then explicitly solved at each time step. In order to illustrate how the shock damper performs, two gravity feed systems with and without a damper are considered. In the first system, the control valve suddenly shuts, and in the second system, rapid demand change occurs in one of the nodes to impose a water hammer condition. Subsequently, unsteady flow parameters such as minimum/maximum flow velocity and pressure are evaluated and a sensitivity analysis is carried out. The results demonstrate that despite being simple and economically efficient, the shock damper is highly capable of moderating unsteady flow characteristics.

One of the key factors to consider in designing a water distribution system is the effect of possible unsteady flow caused by a water hammer. A water hammer is a pressure surge or wave propagating at the speed of sound, which is created when a fluid (usually a liquid) in motion is forced to stop or change direction suddenly (momentum change). This phenomenon can occur for various reasons, including rapid valve opening and/or shutting, pump, joint or pipe failure, and using check valves. Water hammers cause intensive damage in some cases, making it a significant factor to keep in view while designing water distribution systems.

Two approaches to tackle this problem have been reported: (a) taking the necessary precautions to avoid a sudden change in flow velocity (or direction) that would cause a water hammer; and (b) equipping the water system with efficient tools to dampen the harmful impacts of a water hammer. As water hammers seem to be inevitable, the second approach is widely employed. Such tools (e.g., shock dampers) help depreciate the consequences of unsteady flow due to water hammers by compensating for any changes in flow pressure. Surge tanks and safety valves are the two kinds of water hammer damping equipment (Tijsseling 1996; Tijsseling et al. 1996). In order to simulate the dynamic behaviour of surge tanks, they are introduced as a boundary condition for characteristic line equations. Taking the conservation of mass, momentum and energy into account augments the number of equations in a system to n, where n is the number of unknown surge tank variables. Detailed discussions on this topic have been presented by Larock et al. (1999) and Skulovich et al. (2015, 2016). Safety valves are also applied as a boundary condition, and additional equations are derived to form a system of equations. Niełacny (2004) modelled spring safety valves in a water distribution system with pumps.

High operational cost is a major concern with utilizing surge tanks. Despite being expensive, one-way surge tanks are unable to dampen the positive pressure wave caused by a water hammer. They are only capable of preventing the negative waves and column separation by injecting flow to the pipeline (Larock et al. 1999; Chamani et al. 2013; Fathi-Moghadam et al. 2013). Air surge tanks can handle both positive and negative waves, but they are costly (Larock et al. 1999; Kim 2008, 2010). On the other hand, safety valves have much lower installation and maintenance costs, yet they perform poorly in restraining positive waves and are incapable of hindering negative waves; thus, column separation and consequent damage to the system are inevitable (Bergant & Simpson 1995; Ghidaoui et al. 2005). For this reason, we aim to find a simple and economic alternative for handling both negative and positive waves.

This study focuses on the shock damper in an attempt to derive its governing equations and illustrate its effectiveness. According to the results, shock dampers are highly capable of moderating unsteady flow characteristics. The rest of this paper is structured as follows. The next section defines the problem under study, which includes deriving and solving the equations for the given boundary conditions. Then the performance of the proposed approach is demonstrated by solving a case study. The results are depicted in the following section and the conclusions comprise the final section.

The current study focuses on the proposed instrument, i.e., the shock damper. It is capable of encountering positive and negative waves caused by water hammers. In this section, the shock damper is demonstrated, the governing equation is derived, boundary conditions are considered and the resulting system of equations is solved.

Shock damper

A shock damper consists of a damper, spring, mass, tank and joint (Figure 1). The spring, damper and mass perform together as a vibrant system with a single degree of freedom (vertically) when positive or negative waves occur. The spring here is pre-pressed to capture sufficient force, which is used to balance the force caused by the pressure head. In the case of positive waves, the excessive pressure imposed on the spring compresses it further, storing potential energy. This potential energy protects the rest of the water distribution system from water hammer effects (waves). When a negative wave hits the shock damper, the pre-pressed spring elongates and water stored in the shock damper's tank is injected into the flow to compensate for the negative pressure. This helps prevent column separation and subsequent damage to the system. Theoretically, after a water hammer occurs, an infinite number of periodic waves must be created in the system. However, this is never the case as a result of friction and energy depreciation. The shock damper acts to accelerate the process of water hammer wave depreciation. Figure 1 illustrates a shock damper and its components, including the (1) damper, (2) spring, (3) mass, (4) joint and (5) tank.

Figure 1

Shock damper and its components.

Figure 1

Shock damper and its components.

Close modal

Governing equations

To simulate the condition of unsteady flow in the system, a pair of partial differential equations (Equations (1) and (2)) are used. These equations are obtained by considering the conservation of mass and momentum for a moving control volume as an element of the pipeline (Chaudhry 1979; Wylie 1984)
formula
(1)
formula
(2)
where H is the piezometric head and V is the flow velocity, both of which are functions of time (t) and space (x); a is the speed of sound wave propagation in the flow; g is the gravitational acceleration; D is the pipe diameter; and f is the Darcy–Weisbach coefficient. All dimensions are in standard SI units. One of the most widely-used approaches for solving equations similar to Equations (1) and (2) is the characteristic line method (Chaudhry 1979). In this method, the solution space is discretized to a finite number of nodes including intermediate and boundary nodes.

Intermediate nodes

The discretized forms of the equations for intermediate nodes in unsteady flow condition are represented below (Chaudhry 1979; Larock et al. 1999):
formula
(3)
formula
(4)
In the above equations, and are the piezometric head and flow rate in node i at time step n+1, respectively, and A is the pipe's cross section.

Boundary conditions

To derive the governing equations of a shock damper, a gravity feed system is considered. Such a system includes a pipe, tank, valve and shock damper, each with a specific equation (boundary condition).

Reservoir

The following equations are considered for a reservoir with constant head located at the beginning of a pipeline with i = 1 node (Larock et al. 1999):
formula
(5)
formula
(6)
If the reservoir is located at the downstream end (node i = N) with constant head , the following equations are considered:
formula
(7)
formula
(8)

Valve at node number i = Bv

Assuming the valve is located in the pipeline at node number i=Bv, the boundary condition equations for the valve that suddenly closes at time Tc are as follows (Larock et al. 1999):
formula
(9)
formula
(10)
formula
(11)
formula
(12)
formula
(13)
formula
(14)
formula
(15)
formula
(16)
formula
(17)
formula
(18)
where are the flow rates at the valve's node and the node after the valve at time step n + 1, which can be calculated with Equation (9). Moreover, the pressure head values at the nodes in each time step are calculated with Equations (10) and (11) using the C+ and C characteristic lines concept. Equation (12) calculates the valve head loss coefficient as a function of the valve's opening percentage , assuming the valve closes as a linear function of time (Equation (13)). Equations (14) and (16) calculate the values of that will be used in other equations. In Equation (15), the values of represent the wave speeds in the pipe before and after the valve. The values of are calculated with Equations (17) and (18) and are then used in other equations. The other parameters were defined previously.

Shock damper at node number i = Bsh

As shown in Figure 2, it is assumed the shock damper is located at node number i=Bsh. Expanding the conservation of mass for this node yields:
formula
(19)
where is the flow rate from the main pipe to the shock damper joint at time step n + 1. Whenever flow enters the shock damper it is assumed to be a positive flow rate, otherwise it is negative.
Figure 2

Shock damper as a boundary condition.

Figure 2

Shock damper as a boundary condition.

Close modal
Next, the conservation of energy is expanded for node i. If the space interval used is small enough, the energy depreciation due to friction in the distance between nodes i and i + 1 along the pipeline is deemed negligible. By using the mean of the values in each time step and assuming the junction pipe is a control volume for solving force equilibrium equation, the control volume is influenced by the force from the top and force from the bottom. Friction force also acts along the junction pipe, so by considering the mentioned forces, the force equilibrium is developed as follows:
formula
(20)
(here, only undefined parameters define other previously defined parameters), where Lc is the length of the joint; is the mass level proportional to the pipe axis in each time step; and is the pressure imposed on the moving mass at each time step. By using the continuity principle, can be calculated as follows:
formula
(21)
Figure 3 shows a free body diagram of any element of mass in the tank. Applying Newton's second law of motion along the vertical axis leads to:
formula
(22)
where M is the moving mass; C is the damping coefficient; Ar is the cross section of the shock damper's tank; Ks is the spring stiffness factor; Z0 is the initial mass height of the pipe's axis at t = 0; and is the distance over which the spring has been compressed at t = 0. The forces are defined as follows: gravity force ; spring force ; damping force and inertia force . The above equations represent the vibration of a system with a single degree of freedom (vertically). in the above equation are the acceleration and velocity of the flow mass, respectively, and are calculated as:
formula
(23)
formula
(24)
Figure 3

Free diagram of motion mass.

Figure 3

Free diagram of motion mass.

Close modal
By using Equation (21) and the relations between position, velocity and acceleration, the velocity and acceleration of the mass in the tank can be calculated based on the flow rate in the joint using the continuity principle as follows:
formula
(25)
formula
(26)
The flow rate and head at nodes i and i + 1 can be estimated by considering the characteristic line equations C and C+ and the equations below:
formula
(27)
formula
(28)
formula
(29)
formula
(30)
In the above equations, are derived with Equations (17) and (18) (with node number i=Bsh). In each time step the nine variables are computed with Equations (20)–(22) and (25)–(30).

Solving the system of equations for a shock damper

The system of equations defined in the previous section is used to simulate the behaviour of a shock damper when a water hammer occurs. The next step is to solve this difficult nonlinear and implicit system of equations. To do so, the term in Equation (20) is replaced with . By replacing in Equation (20) with their corresponding values from Equations (21), (22) and (25)–(30), Equation (20) changes into an explicit linear form with as the only variable:
formula
(31)
In the above equation, new parameters are introduced and defined as follows:
formula
(32)
formula
(33)
formula
(34)

The efficiency of the proposed model is examined against two benchmark examples from the literature. The examples are solved in two scenarios: with and without a shock damper. The aim is to demonstrate the effect of installing a shock damper on the flow characteristics of such systems.

First case study (rapid valve closure)

The example involves transients caused by rapid valve closure downstream of a long conduit with an upstream reservoir (solved by Saikia & Sarma (2006) and Afshar & Rohani (2008)). A water hammer condition is created by shutting the downstream valve at TC = 4 seconds. The initial system condition includes a constant flow rate in the pipe.

Scenario 1

In this case, the system is not equipped with a shock damper. The piezometric head in the tank is 600 ft; the main pipe diameter is 2 ft; the pipe length is 12,000 ft and the Darcy–Weisbach coefficient is f = 0.02 (for both steady and unsteady flow). The steady discharge is 12 ft3/sec, the pressure wave velocity after valve closure is a = 3,000 ft/sec and and are assumed to be 3,000 ft and 1 sec, respectively.

Scenario 2

A shock damper is used in this scenario, which is located 50 ft upstream of the closure valve. The parameters related to the shock damper are:
formula

All other parameters are assumed similar to those in scenario 1.

Second case study (rapid demand change)

This example considers transients caused by rapid demand changes in one of the nodes (solved by Larock et al. 1999). This network is supplied by gravity flow from two elevated reservoirs. The pipe and node numbers are shown in Figure 4 with the pipe diameters and lengths. The pressure wave velocity in all pipes is a = 2,800 ft/sec, and are assumed to be 300 ft and 0.1 sec, respectively.

Figure 4

Schematic of gravity network and parameters in case study 2.

Figure 4

Schematic of gravity network and parameters in case study 2.

Close modal

Scenario 1

In this case, the system is not equipped with a shock damper. Nodes 1 through 4 have ground elevations of 860 ft, while nodes 5 and 6 have ground elevations of 980 ft and 960 ft, respectively. The equivalent sand roughness of all pipes is e = 0.002 in. The demand increases from 450 to 900 gal/min at node 2 to meet the sudden need for more water for fire suppression. This effect indicates an increase in demand on the head of node 4.

Scenario 2

A shock damper is used in this scenario located 50 ft upstream of node 4. The parameters related to the shock damper are similar to those in the first case study, except the tank diameter here is Dr= 5 ft.

As seen in Figures 5 and 6, a water hammer caused by sudden valve closure can generate piezometric heads (and corresponding pressure) of up to 1,200 ft and 100 ft when no shock damper is used. This is relatively severe pressure resulting in extensive damage to the system. Maximum pressure (positive values) can cause pipes to explode. Moreover, the next head peaks reach 1,100 ft and 150 ft due to friction in the pipeline and energy depreciation, although the resulting pressure is high. In this scenario, the impact of the water hammer is still visible after around 300 sec. When a shock damper is utilized as described in scenario 2, the maximum and minimum pressure heads calculated are 720 ft and 540 ft, respectively. This is a greater reduction (40% for the maximum and 80% for the minimum head) compared to scenario 1. Moreover, the peak heads significantly depreciate after almost 50 sec, reaching 650 ft and 550 ft. In the second case study and as indicated in Figure 7, the maximum and minimum pressure heads adjust from 260 ft and 0 ft, respectively, to 165 ft and 55 ft. After about 15 sec, the pressure oscillation also disappears due to the shock damper.

Figure 5

Pressure head vs time at valve node calculated by previous researchers and in the current work for result fitting (first case study – scenario 1).

Figure 5

Pressure head vs time at valve node calculated by previous researchers and in the current work for result fitting (first case study – scenario 1).

Close modal
Figure 6

Pressure head vs time with shock damper at valve node (first case study – scenario 2).

Figure 6

Pressure head vs time with shock damper at valve node (first case study – scenario 2).

Close modal
Figure 7

Pressure head vs time with and without shock damper at node 4 (second case study – scenarios 1 and 2).

Figure 7

Pressure head vs time with and without shock damper at node 4 (second case study – scenarios 1 and 2).

Close modal

A sensitivity analysis of the system is described herein. The results indicate that shock damper performance is more sensitive to three parameters: main damping coefficient, tank diameter, and shock damper spring stiffness. Tables 13 summarize the sensitivity analysis outcome.

Table 1

Maximum pressure change with damping coefficient (Dr= 10 ft and Ks = 1.4 × 107 lb/ft)

C (lb.s/ft) 0.5 × 105 1.0 × 105 1.5 × 105 2.0 × 105 2.5 × 105 
PMAX (ft) 1,100 988 720 814 1,152 
C (lb.s/ft) 0.5 × 105 1.0 × 105 1.5 × 105 2.0 × 105 2.5 × 105 
PMAX (ft) 1,100 988 720 814 1,152 
Table 2

Maximum pressure change with tank diameter (C = 1.5 × 108 lb.sec/ft and Ks = 1.4 × 107 lb/sec)

Dr (ft) 10 12 14 
PMAX (ft) 1,011 860 720 690 573 
Dr (ft) 10 12 14 
PMAX (ft) 1,011 860 720 690 573 
Table 3

Maximum pressure change with spring stiffness (Dr= 10 ft and C = 1.5 × 108 lb·sec/ft)

Ks (lb/ft) 0.6 × 107 1.0 × 107 1.4 × 107 1.8 × 107 2.8 × 107 
PMAX (ft) 1,150 1,023 720 812 990 
Ks (lb/ft) 0.6 × 107 1.0 × 107 1.4 × 107 1.8 × 107 2.8 × 107 
PMAX (ft) 1,150 1,023 720 812 990 

Table 1 shows the maximum pressure head obtained at the valve location for damping coefficients of +66% and −66% of its initial value. Evidently, an increase in the damping coefficient decreased the shock damper efficiency, hence a greater piezometric head was generated. Decreasing this parameter improved the shock damper's performance. However, a limitation is that further parameter reduction would deteriorate efficiency, suggesting the necessity to seek an optimal level. The system responses to main tank diameter changes are demonstrated in Table 2. As expected, the shock damper performed better with increasing tank size. However, there was a feasible upper limit to this expansion due to physical and economic constraints. Besides, beyond a certain level, any further increase in size appeared to have no considerable effect. Again, this implies there was an optimal level for this capacity expansion. In a situation where the diameter approaches zero, a system with a shock damper and main tank diameter of zero is expected to perform similarly to a system without a shock damper. Table 3 presents the maximum pressure head changes with spring stiffness. The table indicates that an increase and decrease from the optimum values increased the maximum head, so there were optimum values at which the maximum heads could be minimum.

Water hammers in water distribution systems are inevitable. Due to the extensive damage caused by this phenomenon, its impact on a system must be moderated efficiently. A shock damper is a simple and economically efficient piece of equipment that can be used to handle devastating consequences such as pipe explosions. The present study concentrated on deriving the governing equation of a shock damper and using it to model the system. Results demonstrated that installing a shock damper in a gravity-feed system can effectively dampen the pressure waves caused by a water hammer. A sensitivity analysis revealed that shock damper performance is mostly dependent on the main tank diameter and the shock damper's damping coefficient. This study also suggests obtaining optimal parameter values in order to maximize performance.

Afshar
,
M. H.
&
Rohani
,
M.
2008
Water hammer simulation by implicit method of characteristic
.
International Journal of Pressure Vessels and Piping
85
(
12
),
851
859
.
Bergant
,
A.
&
Simpson
,
A. R.
1995
Water Hammer and Column Separation Measurements in an Experimental Apparatus
.
Report No. R128
,
Dept. of Civil and Environmental Engineering University of Adelaide
,
Adelaide
,
Australia
.
Chamani
,
M.
,
Pourshahabi
,
S.
&
Sheikholeslam
,
F.
2013
Fuzzy genetic algorithm approach for optimization of surge tanks
.
Scientia Iranica
20
(
2
),
278
285
.
Chaudhry
,
H. M.
1979
Applied Hydraulic Transients
.
Springer
,
New York
,
USA
.
Fathi-Moghadam
,
M.
,
Haghighipour
,
S.
&
Vali Samani
,
H. M.
2013
Design-variable optimization of hydropower tunnels and surge tanks using a genetic algorithm
.
Journal of Water Resources Planning and Management
139
(
2
),
200
208
.
Ghidaoui
,
M. S.
,
Zhao
,
M.
,
McInnis
,
D. A.
&
Axworthy
,
D. H.
2005
A review of water hammer theory and practice
.
Applied Mechanics Reviews
58
(
1
),
49
76
.
Larock
,
B. E.
,
Jeppson
,
R. W.
&
Watters
,
G. Z.
1999
Hydraulics of Pipeline Systems
.
CRC Press
,
New York
,
USA
.
Niełacny
,
M.
2004
Model of the water-hammer effect considering a spring safety valve
.
Archives of Hydro-Engineering and Environmental Mechanics
51
(
1
),
25
40
.
Saikia
,
M. D.
&
Sarma
,
A. K.
2006
Numerical modeling of water hammer with variable friction factor
.
APRN Journal of Engineering and Applied Sciences
1
(
4
),
35
40
.
Skulovich
,
O.
,
Bent
,
R.
,
Judi
,
D.
,
Perelman
,
L. S.
&
Ostfeld
,
A.
2015
Piece-wise mixed integer programming for optimal sizing of surge control devices in water distribution systems
.
Water Resources Research
51
(
6
),
4391
4408
.
Skulovich
,
O.
,
Perelman
,
L. S.
&
Ostfeld
,
A.
2016
Optimal closure of system actuators for transient control: an analytical approach
.
Journal of Hydroinformatics
18
(
3
),
393
408
.
Tijsseling
,
A.
1996
Fluid-structure interaction in liquid-filled pipe systems: a review
.
Journal of Fluids and Structures
10
(
2
),
109
146
.
Tijsseling
,
A.
,
Vardy
,
A.
&
Fan
,
D.
1996
Fluid-structure interaction and cavitation in a single-elbow pipe system
.
Journal of Fluids and Structures
10
(
4
),
395
420
.
Wylie
,
E. B.
1984
Fundamental equations of waterhammer
.
Journal of Hydraulic Engineering
110
(
4
),
539
542
.