Use of Equivalent Celerity to Estimate Maximum Pressure Increase in Serial Pipes During Water Hammer - Numerical Simulations in MATLAB

Use of Equivalent Celerity to Estimate Maximum Pressure Increase in Serial Pipes During Water Hammer - Numerical Simulations in MATLAB

Agnieszka Malesińska Mariusz Rogulski Pierfabrizio Puntorieri Giuseppe Barbaro Beata Kowalska

Warsaw University of Technology Faculty of Building Services, Hydro and Environmental Engineering, Warsaw, Poland

Department of Civil Engineering, Energy, Environmental and Materials, Mediterranea University, RC, Italy

Lublin University of Technology, Faculty of Environmental Engineering, Lublin, Poland

Available online: 
| Citation



Pipe lines are useful for transporting water for drinking, irrigation and for fire-ing over long distances, this pipe lines are called “Transmission line” and are used to carry conveying raw or treated water from a well field or remote storage (large lake, reservoir, etc.,) facility to a treatment plant and/or distribution storage tank. In water-carrying piping systems, dangerous phenomena may occur. One such phenomenon is water hammer.

The water hammer has always been an area of study, which has captivated the minds of research- ers due to its complex and challenging phenomena. Modeling the phenomenon in real conditions is extremely difficult. Due to the dimensions of the piping systems, conducting research at real scales is impossible. However, thanks to the development of numerical methods, the study of water hammer and its effects can be performed using simulation programs. Unfortunately, the simulation results are not always consistent with the actual course of the phenomenon.

One of the parameters that describes the nature of the course of a water hammer is the velocity of propagation of the pressure wave, c, which is called celerity. The transient surge pressure, p, may be calculated from the pressure celerity c, and the sudden change in fluid flow velocity, Δv. In a piping system, the value of the pressure wave celerity is not equal to the individual celerity, c, for a single pipeline. Therefore for piping systems for Δp calculations the equivalent celerity shell be used.

This article presents value of the equivalent celerity calculated from equations derived using linear analysis of natural vibrations of the system. For implement of the equations, an algorithm in MAT- LAB has been developed that allows one to easily calculate the equivalent celerity, ce, for N pipelines connected in series with varying diameter, length and material composition.


equivalent celerity, long distance, numerical methods, water hammer


[1] Shamloo, H. & Mousavifard, M., Numerical simulation of turbulent pipe flow for water hammer. ASME Journal of Fluids Engineering, 137(11), 2015.

[2] Libraga, J., Ribeiro, R., Covas, D. & Ramos, H., Experimental and Numerical Analysis of Water Hammer in a Metal Pipe Rig. 11th International Conference on Computing and Control for the Water Industry, Exeter, UK, 2011.

[3] Bruce, S. & David, S., Water Hammer: Practical Solutions, Butterworth-Heinemann, London, 1995.

[4] Meniconi, S., Brunone, B. & Ferrante, M., Water-hammer pressure waves interaction at cross-section changes in series in viscoelastic pipes. Journal of Fluids and Structures, 33, pp. 44–58, 2012.

[5] Wylie, E.B. & Streeter, V.L., Fluid Transient in Systems, Prential Hall. Inc., Englewood Cliffs, New York, 1993.

[6] Puntorieri, P., Barbaro, G., Martins, N., Covas, D. & Vincenzo, F., Hydraulic Transient Experimental Study in a Copper Pipe. 9th International Conference on Computational & Experimental Methods in Multiphase & Complex Flow, (In Press), 2017.

[7] Malesińska, A., Experimental study of water hammer-induced forces and deformations in dry pipe fire protections sytems. Fire Safety Journal, 72, pp. 16–24, 2015.

[8] Wylie, E.B. & Streeter, V.L., Fluid Transient In Systems. Prential Hall. Inc., Englewood Cliffs, New York, 1993.