Research Article  Open Access
A Water Hammer Protection Method for Mine Drainage System Based on Velocity Adjustment of Hydraulic Control Valve
Abstract
Water hammer analysis is a fundamental work of pipeline systems design process for water distribution networks. The main characteristics for mine drainage system are the limited space and high cost of equipment and pipeline changing. In order to solve the protection problem of valveclosing water hammer for mine drainage system, a water hammer protection method for mine drainage system based on velocity adjustment of HCV (Hydraulic Control Valve) is proposed in this paper. The mathematic model of water hammer fluctuations is established based on the characteristic line method. Then, boundary conditions of water hammer controlling for mine drainage system are determined and its simplex model is established. The optimization adjustment strategy is solved from the mathematic model of multistage valveclosing. Taking a mine drainage system as an example, compared results between simulations and experiments show that the proposed method and the optimized valveclosing strategy are effective.
1. Introduction
Mine drainage system is an important part in the safety production of coal mine [1]. Due to the limited space and the high cost of equipment changing, water hammer is a common phenomenon in mine drainage system and its harm is inestimable [2, 3]. The minor injuries of water hammer can cause severe shock or even pipes breaking, and the major injuries can cause equipment damaging, pumping station flooding, or even injuries to underground staff. There are several traditional water hammer protection methods such as installing vacuum valves, exhaust valve, and pressure tank. However, when these water hammer protection methods are used to mine drainage system, the original piping arrangement has to be changed because additional equipment is needed. It is unenforceable to do this in this limited mine space and the cost is too high. In this situation, controlling the time and velocity of valveclosing is an effective means to protect water hammer in mine drainage system. A water hammer is easily to be formed in pipeline if the valveclosing is fast. On the contrary, the capability of pumps is inefficient if the valveclosing is too slow because pumps in mine drainage system are centrifugal. This inefficient operation has damage for pumps, because power provided by pumps is converted into heat. Therefore, it is important for us to research the critical valveclosing velocity for Hydraulic Control Valve (HCV). Researches of theoretical system and engineering applications for water hammer protection of pipeline fluid delivery areas are increasingly improved [4–6]. But it is still a difficult problem in mine drainage system because of the limited mine space [7, 8].
At present, the theories and methods of water hammer and its protection means become more and more mature. In summary, researches about this problem are mainly in three aspects: hydraulic transient modeling, the calculation, and protection methods of water hammer.
(1) Hydraulic Transient Modeling. For the research on hydraulic transient, from the mathematical derivation of the 18th century to the graphical analysis of the mid20th century and to the current computer digital simulation, scholars have already made a lot of research results. The major achievements are getting the relationship between multiphase and multicomponent transient flow equations, water hammer equations, and the control equations, such as Joukowsky equation [9]. Based on the transient flow simulation theory, Colombo et al. [10] proposed an aqueducts fault detection technology, Lee et al. [11] proposed the pipe network leak and deterioration over time detection technology by the time domain reflectometry (TDR), Arbon et al. [12] proposed pipeline corrosion and blockage detection technology, Gong et al. [13] proposed a detection technology for pipe friction, wall thickness, velocity, position, and the length of the pipes, and Ferrante et al. [14] presented a leak detection method with coupling wavelet analysis and a Lagrangian model techniques. Meniconi et al. [15] presented a pipe system diagnosis method with the small amplitude sharp pressure waves.
(2) The Calculation Methods of Hydraulic Transient for Water Hammer. The calculation methods of hydraulic transient for water hammer include arithmetic method, graphic method, and numerical method.
(i) Arithmetic Method. Before the 1930s, the hydraulic transient calculation of water hammer used Allievi equations mostly [16]. Allievi equations can be called Arithmetic method which is used to solve the problems of water hammer that with simple boundary conditions and its workload is very large.
(ii) Graphic Method. Graphic method is developed in 1930s to 1960s. Bergeron, Parmakian, and so forth [17] are committed to develop this method. Boundary conditions and the process of water hammer fluctuation are expressed through coordinate graphics of and according to this method. Due to the graphics, it is simple and intuitive for the hydraulic transient calculation of water hammer. However, the accuracy is not high because this method is restricted by calculating means and assumptions.
(iii) Numerical Method. From 1960s, some numerical methods appeared that can be aided by computers, such as Characteristic (MOC) [18], Wave Characteristic Method (WCM) [19], Implicit Method [20, 21], and Finite Element Method (FEM) [22, 23]. The WCM can solve water hammer problems of complex piping systems and boundary conditions. It is the most common method because of the high accuracy and computing. The Implicit Method divides pipeline into several segments and solves equations of the entire pipeline system simultaneously in each segment. The advantage of Implicit Method can be described in a way that a longer segment is selected and the number of calculations is reduced. However, there is more time needed for calculation in large and complex pipeline network system [24]. FEM with flexibility is used in pipe network system which have complex boundary conditions. However, it has a limitation in solving hydraulic transient problems.
(3) Water Hammer Protection. Prevention and controlling means for water hammer are researched with the development of their theory and calculation. Wylie [25–27] has researched several protective devices for water hammer, such as air valves, check valves, pressure tank, and surge tank. Lee [28] and Stephenson [29] discussed the performances of air valves in water hammer protection. However, these researches have not given a quantitative calculation for the problem of water hammer protection.
This paper proposed a water hammer protection method based on velocity adjustment of HCV to deal with the problem of valveclosing water hammer in mine drainage system. The mathematic model of water hammer fluctuations is founded based on MOC according to the hydraulic transient. Then, boundary conditions of water hammer controlling for mine drainage system are determined and the simplex model is established. Finally, the optimization adjustment strategy is solved with simulation and experiment.
The remainder of this paper is organized as follows. Section 2 provides the mathematic model of the propagation and superposition model for water hammer fluctuations. Section 3 provides an optimization method to determine the adjustment strategy for HCV. Section 4 presents a case study. Concluding remarks are offered in the last section of this paper.
2. The Process of Water Hammer Fluctuations
2.1. Mathematic Model
(1) Foundation Equations. The momentum equations of transient flow can be expressed aswhere is longitudinal mean velocity, is pressure head, is friction factor, is diameter of pipeline, is gravitational acceleration, is propagation distance along pipe, and is time.
The left of the equation represents the timevarying inertia force in a unit volume. In the right, the first item represents the pressure in a unit volume of fluid. The second item is friction loss pressure in a unit length of pipe, and the third item is pressure of transient flow.
The continuity equation is used to describe transient flow and it can be expressed as follows:where is propagation speed of water hammer wave and is the angle between the axis of pipe and the horizontal line.
(2) The Process of Water Hammer Fluctuations. Pressures that are caused by friction loss and transient flow velocity are ignored. Furthermore, is considered for continuity equations. Then, the foundation equations can be simplified as
The secondorder partial differential equations about and are obtained through the partial derivative and are taken for variables and in (3). Therefore, the momentum equations for water hammer wave in pipeline drainage system are expressed as
The functions and are introduced along with characteristic lines . Then, the general solution of (4) can be described aswhere is the initial pressure head, is the initial flow velocity, and and are direct and reflection fluctuation functions.
2.2. Propagation and Superposition Model of Water Hammer Wave Based on Characteristic Line Method
Pumps and valves are the generating sources of water hammer wave in hydraulic transition process of pipeline fluid delivery system. These two kinds of water hammer wave start at the same time and propagation directions are superimposed. The superposition leads to strengthening of fluctuations for water hammer. In order to control the water hammer pressure effectively, the superposition effects of water hammer wave need to be weakened as much as possible.
The fundamental equations of water hammer can be rewritten as
These equations are combined linearly using an unknown multiplier . Let be the linear combination. The coefficient can be determined by
In the constraint of characteristic line equation , (6) converts to ordinary differential equation as follows:
and are negligible when and the tilt angle of pipeline is less than 25°. The average flow velocity is replaced by flow . The characteristic equations are changed intowhere is the characteristic line of forward wave in axis, is the characteristic line of reflected wave in axis, and is crosssectional area of pipe.
The integral operation and differential conversion are introduced to (9). The discrete characteristic equations of water hammer are obtained as follows. Equation (10) and the characteristic line grid are shown in Figure 1:
By solving the above equations, transient variables and in the node of the pipe over the time are expressed aswhere .
Consider
The initial condition of water hammer is that of a steady flow, and in (11) are known, and the parameters on each time step can be solved.
Hydraulic characteristics of HCV can be described through loss coefficient and capacity coefficient . Consider where is crosssectional area of the upstream of HCV and is the pressure difference of HCV in upper and lower.
Let be the relative opening displacement of HCV. The dimensionless flow through HCV can be described aswhere is pressure difference on HCV. It can be solved as . The subscript “0” represents HCV being open completely.
3. Optimization Method for Velocity Adjustment of HCV
3.1. Boundary Conditions
Boundary conditions for pump of outlet are shown in Figure 2. If the inlet loss of pump is ignored, the inlet of pool, pump, and HCV can be combined and dealt with as a boundary of characteristic line in foundational equations of water hammer.
The equilibrium condition for water head is established in this boundary. It can be described aswhere is the depth of pool inlet, is the work lift of the pump, is the transient resistance loss of HCV, and is the pressure head of HCV outlet.
We take the balance equations of water head which in the characteristic line and inertia equations of pump unit together. And then the boundary conditions can be expressed as follows based on the characteristic lines of pumps:where the pumplift is . The flow of pump is . The rotational speed is . The torque is . The subscript represents rated working conditions. Then, (16) is generated as follows:where ; ; and is rotational inertia. and are rotational speed and torque in the initial state. Then, (17) can be rewritten as follows when variables and are introduced:The problem is solved as follows:in which and are solved through (18).
The iterative errors should meet (19). Consider
3.2. Optimization Method of HCV Adjustment
The boundary conditions of multistage valveclosing are introduced in the characteristic line calculation of water hammer for which the best procedure for multistage valveclosing is considered. The displacement of valveclosing at any time can be determined through linear procedure in multistage valve closure.
Through a large number of practical tests, the valve closure with a uniform or uniform variable speed is found to have a small influence on the water hammer in the front 2/3 of valveclosing stroke, while there is a large influence on the later 1/3 of valveclosing stroke. Therefore, the valve closure process is divided into multistage. By doing so, the uniform and uniform variable speeds are both introduced to valve closure process. In the following section we will discuss uniform or uniform variable speed calculation method in valve closure process.
Assuming that the times of valveclosing in multistage are , the corresponding displacements of valveclosing are . Then, the displacements of valveclosing at any time are determined as follows:where is displacement of valveclosing at and is the total stroke of valveclosing.
For the valveclosing displacement at any time , we can calculate the dimensionless opening degree using (21) through linear interpolation according to equidistant and its input value:
In the scope of velocity adjustment for HCV, the velocity, displacement, and time of valveclosing have various combinations. This paper adopts a multistage adjustment strategy according to simulation results of water hammer.
The valve closure time for the first stage is equal or close to the time of the flow of pump that is zero. The valve closure time for the following stages is one in four of the previous stage. Consequently, the initial valve closure procedure is determined. Then, search computing is taken to determine the optimization valve closure procedure.
4. Case Study and Experiment
4.1. Working Conditions
The mine with normal water quality was considered. The water normal inflow was 50 m^{3}/h and the maximum inflow was 250 m^{3}/h. Water was discharged into ground by two drain pipes of which the specifications were D273 × 12, respectively (one was used and another was spare). The pipes in this existing and antiquated mine drainage system were seamless steel pipe. The allowed pressure was as lower at 6.4 MPa. The wellhead elevation was 997.991 m and the borehole inclination was 20°. There were three multistage centrifugal pumps and two HCVs were used in Figure 3. And the specification of pumps was MD28065 × 6 ( m^{3}/h, m). One pump was used and others were maintenance when the mine was with a normal inflow, or two pumps were used and one was spare when the mine was with a maximum inflow.
(a)
(b)
The automatic drainage system was chosen in this mine. Check valves were installed in pump outlet. There are HCVs in pipeline shown in Figure 3(a) where the stroke was 220 mm with the structure of HCV shown in Figure 4. During normal pump stopping, HCVs were closed first, and then pumps were stopped. The water hammer phenomenon was significant in drains pipeline during this process. Therefore, we should take some measures to prevent the water hammer phenomenon.
4.2. Simulations
An algorithm was developed using commercial simulation software Fluent6.3 based on the water hammer theory described in Section 2. The acceleration, amplitude, pipeline pressure, and valve closure displacement were obtained in the cases of uniform and uniform variable speed through transient simulation for water hammer. The correctness of simulation and the proposed theory were then verified against experimental results. We simulated the valveclosing in singlestage with a constant velocity and in multistage with a variable velocity based on the proposed theory.
(1) The Simulation of ValveClosing in SingleStage with a Constant Velocity. The velocity of valveclosing in singlestage was determined through repeated simulation. Then, times of equipment in this situation were obtained using (20) to (21) as shown in Table 1.

The simulation results of valveclosing parameters in singlestage with are shown in Figures 5–12.
For Figures 5–12, the static pressure of pipeline was at 3.06 MPa for pump stopping. After the vacuum was completed, the pump was started and HCV was open with a constant velocity. At this moment, pipeline pressure was increased along with the valveopening displacement of HCV. Then, pressure fluctuations appeared. The valveopening displacement of HCV was 42.6 mm and the maximum pressure was 3.84 MPa when s. After this moment, the pressure was steady. The valveopening displacement of HCV was 70.2 mm and the maximum pressure was 3.25 MPa when s. The drainage system pressure reached a steady state along with the increasing of valveopening displacement. The pressure of pipeline began to drop when s and valveclosing displacement was 53.4 mm. Water hammer phenomenon appeared when s. When s, the maximum boost of water hammer was m and the pressure of pipeline at this moment was 4.77 Mpa. It is 1.47 times the normal discharge pressure (3.25 MPa) and the pressure increased by 50% suddenly. The acceleration and amplitude of pipeline in directions of , , and were changed significantly when water hammer occurred. Maximum acceleration in directions of , , and was 4.83 m/s^{2}, 6.46 m/s^{2}, and 5.45 m/s^{2} and the maximum amplitude was 15.8 mm, 25.8 mm, and 15.3 mm, respectively. Additionally, the total stroke of HCV was assumed to be . Valveclosing with rapid and constant velocity were taken in the former and a variable velocity was taken in the latter according to the simulation results of valveclosing water hammer. Due to the velocity of valveclosing that has a large impact on water hammer parameters, it is an effective approach to protect water hammer by choosing a reasonable velocity for valveclosing.
The allowed pressure of pipe is lower and its value was 6.4 MPa. However, the pressure of pipeline reached 4.8 MPa in the case of 22 s valveclosing from the above analysis. This value is near to the allowed pressure of pipe with the pressure being not high, and the water hammer phenomenon was not significant compared with other mine drainage systems. In order to reduce the pressure of pipeline to protection pipeline, it is necessary to develop water hammer protection.
(2) The Simulation of ValveClosing in Multistage with a Variable Velocity. The optimization strategy of valveclosing adjustment for HCV is generated using boundary conditions according to optimization method of valveclosing adjustment that is proposed in this paper. It is a multistage valveclosing with a variable velocity. In the front 2/3 of valveclosing stroke, the velocity is constant and its value is . In the later 1/3 of valveclosing stroke, the velocity is a variable (the deceleration is constant) and the initial value is . The total time of valveclosing is 43.823 s. Times of equipment in multistage valveclosing are shown in Table 2.

The simulation results of valveclosing parameters in multistage are shown in Figures 13–20.
The pipeline pressuring, acceleration, and amplitude in the directions of , , and are significantly reduced in multistage valveclosing with a variable velocity. Valveclosing with a constant velocity is started at 46.007 s, and valveclosing with a variable velocity is started at 60.6 s. The HCVs are closed and pumps are stopped at 89.83 s. In this process, the maximum pressure of water hammer is 3.49 MPa. No effect of water hammer occurred under this situation.
4.3. Experiments
The experimental devices include multifunction data acquisition instrument, vibration sensor, pressure sensors, and cable displacement encoder. Multifunction data acquisition instrument has 16 channels and its acquisition frequency is 102.4 kHz. Data of acceleration, amplitude, flow, pressure, and displacement of valveclosing were acquired by this multifunction data acquisition instrument. Vibration sensor was used to measure piping vibration under the influence of water hammer. It was disposed along with the , , and direction of the pipe. Pressure sensors were used to measure pipeline pressure. Its measurement accuracy and range are 0–10 MPa and ±0.5 FS. Cable displacement encoder was used to posit the displacement of valveclosing. The arrangement diagram of sensors and flow chart of pumps are shown in Figures 21 and 22.
HCVs control system was designed in this experiment which is shown in Figure 23. Sensors were used to collect signals of pressure and acceleration in pipeline. PLC unit was used to analyze pipeline pressure and vibration for water hammer on the basis of pressure and acceleration signals. An expected running speed was given to HCVs. And then, voltage, current, power, and fault protection information of pumps, pressure, flow, temperature, vacuum, and HCVs status information were monitored at real time. These multiparameter data were integrated using of fuzzy control theory with calculus and control strategy. The HCVs control system can control starting and stopping of pumps automatically and can monitor the signals in system running process timely.
For the total stroke of valve , the experimental conditions are as follows. In the front 2/3 of the stroke, the velocity was constant and its value was . In the later 1/3, the velocity is variable (the deceleration is constant) and the initial value was . The full valveclosing times were 44 s. Measured values of water hammer parameters in multistage are shown in Table 3.

5. Results and Discussion
The reliability and effectiveness are illustrated through the comparison between the results of simulation and experiment.
(1) Results of ValveClosing in SingleStage with a Constant Velocity and in Multistage with a Variable Velocity Are Compared. Pipeline pressures and displacements of valveopening under these two conditions are shown in Figures 24 and 25, respectively. The red curves represent the variations in singlestage with a constant velocity while the black curves represent the variations for multistage with a variable velocity.
Figures 24 and 25 show a significant water hammer effect in singlestage valveclosing with a constant velocity, and the valveclosing time is shorter (22 s). On the other hand, water hammer phenomenon was found to be considerably improved for multistage valveclosing condition with a variable velocity and the valveclosing time was extended (44 s). It shows that valveclosing in multistage with a variable velocity can control the water hammer phenomenon effective in mine drainage system.
(2) Results of the Simulation and Experiment Are Also Compared. Results for valveclosing with a variable velocity are shown in Table 4.

It can be seen from Table 4 that the simulation for valveclosing with a variable velocity shows good agreement with the experiment. The errors between simulation and experiment are less than 5%. Therefore, the valveclosing in multistage with a variable velocity can reduce the impact from water hammer in mine drainage system.
6. Conclusions and Future Work
Limited space and high cost of equipment changing are the main characteristics for mine drainage system. For the problem of valveclosing water hammer protection in mine drainage system, this paper proposed a protection method based on velocity adjustment of HCVs. The mathematic model of valveclosing water hammer is established based on fundamental equations. The strategy of valveclosing in multistage with a variable velocity is obtained by velocity adjustment strategy of HCVs, and it is effective in protecting the valveclosing water hammer without change piping arrangement or addition equipment. The results of simulations and experiments in different adjustment velocity are compared. We obtained that the water hammer phenomenon is not obvious with valveclosing velocity in the front 2/3 of its stroke, and it is very seriously within the later 1/3. Therefore, we take a twostage valveclosing strategy. In the first stage, the velocity is constant and its value is . In the second stage, the velocity is variable (the deceleration is constant) and the initial value is .
Extending errors in the theoretical wave speed have an important influence on prediction of the closure rate before a water hammer occurs in the actual pipeline. A very challenging mathematical task is variable. This is beyond the scope of this paper and may be considered as a future work.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work was supported by the National 12.5 Technology Support Program of China under Grant no. SQ2012BAJY3504.
References
 M. Sun, X.R. Ru, and L.F. Zhai, “Insitu fabrication of supported iron oxides from synthetic acid mine drainage: high catalytic activities and good stabilities towards electroFenton reaction,” Applied Catalysis B: Environmental, vol. 165, pp. 103–110, 2015. View at: Publisher Site  Google Scholar
 A. Kaliatka, M. Vaišnoras, and M. Valinčius, “Modelling of valve induced water hammer phenomena in a district heating system,” Computers & Fluids, vol. 94, pp. 30–36, 2014. View at: Publisher Site  Google Scholar
 N. Jeong, D. Chi, and J. Yoon, “Analytic estimation of the impact pressure in a beamtube due to waterhammer effect,” Journal of Nuclear Science and Technology, vol. 50, no. 12, pp. 1139–1149, 2013. View at: Publisher Site  Google Scholar
 W. S. Yi, J. Jiang, D. D. Li, G. Lan, and Z. Zhao, “Pumpstopping water hammer simulation based on RELAP5,” IOP Conference Series: Materials Science and Engineering, vol. 52, no. 7, Article ID 072009, 2013. View at: Publisher Site  Google Scholar
 Q. W. Yong, M. Jiang, and Y. Tang, “Study on liquid pipeline water hammer and protective device with gas concentration,” Advanced Materials Research, vol. 616–618, pp. 774–777, 2013. View at: Publisher Site  Google Scholar
 J. Zhang, J. Gao, M. Diao, W. Wu, T. Wang, and S. Qi, “A case study on risk assessment of long distance water supply system,” Procedia Engineering, vol. 70, pp. 1762–1771, 2014. View at: Publisher Site  Google Scholar
 K. I. M. SangGyun, L. E. E. KyeBock, and K. I. M. KyungYup, “Water hammer in the pumprising pipeline system with an air chamber,” Journal of Hydrodynamics, Series B, vol. 26, no. 6, pp. 960–964, 2015. View at: Publisher Site  Google Scholar
 M. A. Bouaziz, M. A. Guidara, C. Schmitt, E. HadjTaïeb, and Z. Azari, “Water hammer effects on a gray cast iron water network after adding pumps,” Engineering Failure Analysis, vol. 44, pp. 1–16, 2014. View at: Publisher Site  Google Scholar
 M. S. Ghidaoui, M. Zhao, D. A. McInnis, and D. H. Axworthy, “A review of water hammer theory and practice,” Applied Mechanics Reviews, vol. 58, no. 1–6, pp. 49–76, 2005. View at: Publisher Site  Google Scholar
 A. F. Colombo, P. Lee, and B. W. Karney, “A selective literature review of transientbased leak detection methods,” Journal of HydroEnvironment Research, vol. 2, no. 4, pp. 212–227, 2009. View at: Publisher Site  Google Scholar
 P. J. Lee, M. F. Lambert, A. R. Simpson, J. P. Vítkovský, and D. Misiunas, “Leak location in single pipelines using transient reflections,” Australian Journal of Water Resources, vol. 11, no. 1, pp. 53–65, 2007. View at: Google Scholar
 N. S. Arbon, M. F. Lambert, A. R. Simpson, and M. L. Stephens, “Field test investigations into distributed fault modeling in water distribution systems using transient testing,” in Proceedings of the World Environmental and Water Resources Congress, Tampa, Fla, USA, May 2007. View at: Google Scholar
 J. Gong, A. R. Simpson, M. F. Lambert, A. C. Zecchin, Y.I. Kim, and A. S. Tijsseling, “Detection of distributed deterioration in single pipes using transient reflections,” Journal of Pipeline Systems Engineering and Practice, vol. 4, no. 1, pp. 32–40, 2013. View at: Publisher Site  Google Scholar
 M. Ferrante, B. Brunone, and S. Meniconi, “Leak detection in branched pipe systems coupling wavelet analysis and a Lagrangian model,” Journal of Water Supply: Research and Technology, vol. 58, no. 2, pp. 95–106, 2009. View at: Publisher Site  Google Scholar
 S. Meniconi, B. Brunone, and M. Ferrante, “Waterhammer pressure waves interaction at crosssection changes in series in viscoelastic pipes,” Journal of Fluids and Structures, vol. 33, pp. 44–58, 2012. View at: Publisher Site  Google Scholar
 W. H. Hager, “Swiss contributions to water hammer theory,” Journal of Hydraulic Research, vol. 39, no. 1, pp. 3–7, 2001. View at: Publisher Site  Google Scholar
 M. S. Ghidaoui and B. W. Karney, “Modified transformation and integration of 1D wave equations,” Journal of Hydraulic Engineering, vol. 121, no. 10, pp. 758–760, 1995. View at: Publisher Site  Google Scholar
 A. Adamkowski and M. Lewandowski, “Cavitation characteristics of shutoff valves in numerical modeling of transients in pipelines with column separation,” Journal of Hydraulic Engineering, vol. 141, no. 2, 2015. View at: Publisher Site  Google Scholar
 E. Yao, G. Kember, and D. Hansen, “Analysis of water hammer attenuation in applications with varying valve closure times,” Journal of Engineering Mechanics, vol. 141, no. 1, Article ID 04014107, 2014. View at: Publisher Site  Google Scholar
 R. Korbar, Z. Virag, and M. Šavar, “Truncated method of characteristics for quasitwodimensional water hammer model,” Journal of Hydraulic Engineering, vol. 140, no. 6, Article ID 04014013, 2014. View at: Publisher Site  Google Scholar
 M. Rohani and M. H. Afshar, “Simulation of transient flow caused by pump failure: pointimplicit method of characteristics,” Annals of Nuclear Energy, vol. 37, no. 12, pp. 1742–1750, 2010. View at: Publisher Site  Google Scholar
 A. Triki, “Multiplegrid finite element solution of the shallow water equations: water hammer phenomenon,” Computers & Fluids, vol. 90, pp. 65–71, 2014. View at: Publisher Site  Google Scholar
 A. Keramat, A. S. Tijsseling, Q. Hou, and A. Ahmadi, “Fluidstructure interaction with pipewall viscoelasticity during water hammer,” Journal of Fluids and Structures, vol. 28, pp. 434–455, 2012. View at: Publisher Site  Google Scholar
 S.G. Kim, K.B. Lee, and K.Y. Kim, “Water hammer in the pumprising pipeline system with an air chamber,” Journal of Hydrodynamics, vol. 26, no. 6, pp. 960–964, 2014. View at: Publisher Site  Google Scholar
 E. B. Wylie and V. L. Streeter, Fluid Transients, McGrawHill, New York, NY, USA, 1978.
 W. Wan, W. Huang, and C. Li, “Sensitivity analysis for the resistance on the performance of a pressure vessel for water hammer protection,” Journal of Pressure Vessel Technology—Transactions of the ASME, vol. 136, no. 1, Article ID 011303, 2014. View at: Publisher Site  Google Scholar
 G. De Martino, F. De Paola, N. Fontana, and M. Giugni, “Discussion of ‘simple guide for design of air vessels for water hammer protection of pumping lines’ by D. Stephenson,” Journal of Hydraulic Engineering, vol. 130, no. 3, pp. 273–275, 2004. View at: Publisher Site  Google Scholar
 T. S. Lee, “Air influence on hydraulic transients on fluid system with air valves,” Journal of Fluids Engineering, vol. 121, no. 3, pp. 646–650, 1999. View at: Publisher Site  Google Scholar
 D. Stephenson, “Effects of air valves and pipework on water hammer pressures,” Journal of Transportation Engineering, vol. 123, no. 2, pp. 101–106, 1997. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Yanfei Kou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.