Research Article  Open Access
Water Hammer Modelling and Simulation by GIS
Abstract
This work defined an Eulerianbased computational model compared with regression of the relationship between the dependent and independent variables for water hammer surge wave in transmission pipeline. The work also mentioned control of UnaccountedforWater (UFW) based on the Geography Information System (GIS) for water transmission pipeline. The experimental results of laboratory model and the field test results showed the validity of prediction achieved by computational model.
1. Introduction
Water hammer phenomena occurring during water hammer are explained on the basis of compressibility of liquid. Many researchers have made significant contributions in this area. Zhukovsky introduced the concept of the effective sound speed. He mentioned reducing the motion of a compressible fluid in an elastic cylindrical pipe to the motion of a compressible fluid in a rigid pipe, but with a lower modulus of elasticity of the liquid. Subjects of transients in liquids are still growing fast around the world. Scientists have developed various methods of investigation of transient pipe flow. These ranges of methods are included by approximate equations to numerical solutions of the nonlinear NavierStokes equations. They obtained the differential equations of motion of inviscid fluid forming the basis for further development of the theory of pressure and pressure flow of viscous fluid. By helping of this theory, it became possible to explain of the physical phenomenon, known as water hammer. They introduced the concept of the effective sound speed. Therefore transient flow was solved for the pipeline in the range of approximate equations. These approximate equations are solved by numerical solutions of the nonlinear NavierStokes equations in a method of characteristics (MOC). So, experiences are ensured for the reliable water transmission pipeline. Numerical modeling and simulation which are defined by method of characteristics (MOC) provides a set of results. The (MOC) approaches transform the water hammer partial differential equations into the ordinary differential equations along the characteristic lines defined as the continuity equation, and the momentum equation are needed to determine and in a onedimensional flow system. Solving these two equations produces a theoretical result that usually corresponds quite closely to actual system measurements based on Geography Information System (GIS) if the data and assumptions used to build the (GIS) already are valid.
Curve estimation for the experimental results of laboratory model and the field test results is the most appropriate when the relationship between the dependent variable and the independent variable is not necessarily linear. Linear regression is used to model the value of a dependent scale variable based on its linear relationship to one or more predictors. Nonlinear regression is appropriate when the relationship between the dependent and independent variables is not intrinsically linear. Binary logistic regression is most useful in modeling of the event probability for a categorical response variable with two outcomes. The autoregression procedure is an extension of ordinary leastsquares regression analysis specifically designed for time series. One of the assumptions underlying ordinary leastsquares regression is the absence of autocorrelation in the model residuals. Time series, however, often exhibit firstorder autocorrelation of the residuals. In the presence of autocorrelated residuals, the linear regression procedure gives inaccurate estimates of how much of the series variability is accounted for by the chosen predictors. This can adversely affect the choice of predictors, and hence the validity of the model [1].
2. Materials and Methods
A GIS ready model for liquidvapor flows illustrates the numerical techniques for solving the resulting equations. Hence field test model was chosen for experimental presentation of water hammer phenomenon at the water pipeline: Partial differential equation (1) are solved by method of characteristics MOC.
The method of characteristics is a finite difference technique in which pressures were computed along the pipe for each time step.
Calculation automatically subdivided the pipe into sections (intervals) and selected a time interval for computations; equations are the characteristic equations (2), (4). If the pressure at the inlet of the pipe and along its length is equal to , then slugging pressure undergoes a sharp increase: The Zhukousky formula is as flows: The speed of the shock wave is calculated by the formula: For the velocity of surge or pressure wave in an elastic case with low value of free water bubble, the equation (6) would be valid: The velocity of pressure wave (7) in an elastic case with the high value of free water bubble is presented by the flowing equation [2, 3]:
2.1. Regression Model due to Field Test for Surge Wave Velocity
The autoregression procedure accounts for firstorder autocorrelated residuals. It provides reliable estimates of both goodnessoffit measures and significant levels of chosen predictor variables.
The variables are as follows: โvelocity of surge wave (m/s) as a dependent variable with nomenclature โ,โ independent variable with nomenclature โโ such as โpercent of air volume (m^{3}). The curve estimation procedure allows quick estimating regression statistics, and producing related plots for different models. Hence the autoregression procedure by regression software โSPSS 10.0.5โ was selected for the curve estimation procedure in the present work. The regression model has been built based on field test data.
Regression software โSPSSโ has fitted the function curve and provided regression analysis. So, the regression model has been found in the final procedure. By this model, field test results have been compared by computational model results. The main practical aim of the present work was concentrated on the definition of a condition base maintenance (CM) method for all water transmission systems. In this work the data collection procedures were as follows.
At fast transients, down to 1 second, surge pressure and velocity of surge wave were recorded. These data were used for curve estimation procedure. They were detected on actual systems (field tests). Also, flow and pressure were computed by computational model. Those data have been compared with flow and pressure data which have been collected from actual systems (field tests). The model is calibrated using one set of data, without changing parameter values. It is used to match a different set of results. Curve estimation procedure is illustrated in Figures 1 and 2 for surge wave velocity was formed by estimating regression statistics which are listed in Tables 1 and 2. However, related plots for the field test model were produced [4โ6].
 
^{
a}Predictors: constant, air volume percent. 

3. Results and Discussion
In this work conclusions were drawn on the basis of experiments and calculations for the pipeline with a local leak. Hence, the most important effects that were observed are as follows. The pressure wave speed generated by water hammer phenomenon was influenced by some additional factors. Therefore the ratio of local leakage and discharge from the leak location was mentioned. The effect of total discharge from the pipeline and its effect on the values of wave oscillations period were studied. The outflow to the surge tank from the leak affected the value of wave celerity. The pipeline was equipped with the valve at the end of the main pipe, which was joined with the closure time register. The water hammer pressure characteristics were measured by extensometers.
Power functions are illustrated in Figure 3; however, a variable base is raised to a fixed exponent. The parameter serves as a simple scaling factor, moving the values of up or down as increases or decreases, respectively, and the parameter , called either the exponent or the power, determines the functionโs rates of growth or decay [7, 8].
Chaudhry [9] obtained pressure heads by the steady model which is illustrated in Figure 4. Comparison showed similarity in present work and work of Chaudhry.
(a)
(b)
4. Conclusions
This work focused on the effects of the penetrated air on the surge wave velocity in water pipeline. It showed that Eulerianbased computational model is more accurate than the regression model. Hence in order to present importance of penetrated air on water hammer phenomenon, it was compared the models for laboratory; computational and field tests experiments. At these procedures, it was showed that the Eulerian based model for water transmission line. It was compared with the regression model. On the other hand, this idea were included the proper analysis to provide a dynamic response to the shortcomings of the system. It also performed the design protection equipments to manage the transition energy and determine the operational procedures to avoid transients. Consequently, the results of this work will help to reduce the risk of system damage or failure at the water pipeline.
Nomenclatures
Velocity of surge wave (m/s)  
:  Percent of air volume (m) 
:  Wall thickness of pipe (mm) 
:  Acceleration of gravity (m/s) 
:  Head of the liquid (water) column,โโ(m) 
:  Modulus of elasticity for pipeline material Steelโโ, (kg/m), 
:  Outer diameter of the pipeโโ(mm) 
:  Densityโโ(kg/m) 
:  Surge pressure 
:  Module of elasticity of waterโโ,โโ(kg/m) 
:  Time (S) 
:  Velocity (m/s) 
:  Length (m) 
:  Diameter of each pipe (mm) 
:  Elevationhead (m) 
:  DarcyWeishachcoefficient. 
Acknowledgments
The authors thank all specialists for their valuable observations and advice and the referees for recommendations that improved the quality of this paper.
References
 K. Hariri Asli, F. B. Nagiyev, and A. K. Haghi, โInterpenetration of two fluids at parallel between plates and turbulent moving in pipe,โ in Computational Methods in Applied Science and Engineering, chapter 7, pp. 115โ128, Nova Science, New York, NY, USA, 2009. View at: Google Scholar
 T. S. Lee and S. Pejovic, โAir influence on similarity of hydraulic transients and vibrations,โ Journal of Fluids Engineering, Transactions of the ASME, vol. 118, no. 4, pp. 706โ709, 1996. View at: Google Scholar
 K. Hariri Asli,, F. B. Nagiyev, and A. K. Haghi, โSome aspects of physical and numerical modeling of water hammer in pipelines,โ International Journal of Nonlinear Dynamics and Chaos in Engineering Systems, vol. 60, no. 4, pp. 677โ701, 2010. View at: Google Scholar
 K. Hariri Asli, F. B. Nagiyev, and A. K. Haghi, Physical Modeling of Fluid Movement in Pipelines, Nanomaterials Yearbook, New York, NY, USA, 2009.
 A. Kodura and K. Weinerowska, โThe influence of the local pipeline leak on water hammer properties,โ Materials of the II Polish Congress of Environmental Engineering, Lublin, Poland, 2005. View at: Google Scholar
 M. S. Ghidaoui, S. G. S. Mansour, and M. Zhao, โApplicability of quasisteady and axisymmetric turbulence models in water hammer,โ Journal of Hydraulic Engineering, vol. 128, no. 10, pp. 917โ924, 2002. View at: Publisher Site  Google Scholar
 K. Hariri Asli, โGIS and water hammer disaster at earthquake in rasht water pipeline,โ in Proceedings of the 3rd International Conference on Integrated Natural Disaster Management (INDM '08), Tehran, Iran, 2008, http://www.civilica.com/PaperINDM03INDM03_001.html. View at: Google Scholar
 K. Hariri Asli, โGIS and nonlinear dynamics model: some computational aspects and practical hints,โ International Journal on Technical and Physical Problems of Engineering, vol. 2, no. 5, pp. 1โ5, 2010. View at: Google Scholar
 M. H. Chaudhry, Applied Hydraulic Transients. Van Nostrand Reinhold, New York, NY, USA, 1987.
Copyright
Copyright © 2012 K. Hariri Asli 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.