What is the Wave Speed Calculation to Analyze Water Hammer in Piping?
Water hammer is a pressure surge or wave caused when a fluid in motion is forced to stop or change direction suddenly (momentum change).
Each series of articles are written by pipe flow analysis engineers from Applied Flow Technology. As industry leaders in water hammer and surge analysis, AFT has collected models and data from projects around the world to use as reference materials for published technical papers, case studies, and blogs. Visit www.aft.com for more information on analysis tools.
Water Hammer Concepts:
What is Wave Speed?
The wave speed is the propagation speed at which waves will be transmitted through a medium. When performing water hammer analysis, the wave speed through the fluid must be known to calculate how fast pressure waves will propagate through the pipes, as well as the severity of pressure changes due to changes in flow.
How is Wave Speed Calculated?
Several parameters are used for wave speed calculations.
The acoustic velocity, or the velocity of the wave traveling uninterrupted through a medium, is the first factor. The acoustical velocity is calculated as follows:
Where
K = fluid bulk modulus of elasticity
ρ = fluid density
The wave speed in a piping system will be slower than the acoustic velocity as the pressure waves will be distorted by the restriction of the pipe walls, losing energy in the process of propagation. The effect of distortions in the pressure waves due to the pipes can be accounted for by modifying the acoustical velocity as follows:
Where
E = pipe modulus of elasticity
D = pipe inner diameter
e = wall thickness
c1 = pipe support constant
It should be noted that the derivation for Equation 2 makes a few assumptions, including:
- The pipe is rigid
- The pipe is essentially liquid full
- The pipe is circular
For cases that do not meet these requirements Thorley provides extensive resources on how to adjust the wave speed, which is detailed further in the next section.
The equations to calculate the pipe support constant (c1) can be found in Wylie & Streeter [1]. The pipe support constant adjusts the wave speed equation to accurately account for the pipe strains and stresses due to elongation of the material when different pipe supports are used. The pipe material poisson’s ratio and/or the pipe geometry are used to calculate c1, depending on what supports are used. The pipe support constant will be 1 for thin-walled pipe anchored throughout with expansion joints, and will be slightly higher or lower for other types of supports.
Often engineers do not know the exact composition of the pipe materials, so parameters such as the pipe modulus of elasticity and poisson’s ratio must be estimated. Thus calculations for the wave speed inherently have a high level of uncertainty.
A Numerical Example
Consider a piping network with water at 72˚F (22˚C) and an average pressure of 290 psia (20 bar). From the NIST REFPROP fluid database, the density is 62.3 lbm/ft3 (998 kg/m3) and the bulk modulus is about 3.2 E+05 psia (2.2 E+04 bar).
The water is flowing through a Steel-ANSI STD 18-inch pipe. The pipe schedule shows that the inside diameter is 17.25 inches (43.8 cm), and the wall thickness is 0.375 inches (0.952 cm). The modulus of elasticity found in Nayar’s Piping Handbook is about 2.9 E+07 psia (2.0 E+06 bar). The Poisson ratio is listed as 0.3 in Perry’s Chemical Engineer’s Handbook. Assuming the pipe is thick-walled anchored throughout, the pipe support constant c1 is 0.88.
By entering the above information into Equation 2, the wave speed is approximately 4060 ft/sec.
As there is a high level of uncertainty in the wave speed, a sensitivity analysis for wave speed may be recommended depending on the calculation being performed. Conveniently, much of this information is standardized so water hammer analysis software can determine wave speed easily.
References
[1] Wylie, E. B. and Streeter, V. L., 1993, Fluid Transients in Systems, Prentice Hall, Englewood Cliffs, NJ.