Forgotten Password?

Or login with:

  • Facebook
  • Google
  • Yahoo

Compressible Flow

This analysis of water hammer allows for the compressibility of water and the expansion of the pipe


Although it is usual to consider water as incompressible and of a uniform density, this is clearly not true. It would, for example, be impossible to control the depth of a submersible vessel unless the density of water increased with depth.

Here we consider the impact of compressibility for the purpose of computing water hammer and the subsequent pipe expansion.

Pressure Rise Following Instantaneous Closure


The pressure rise at the valve
=\rho\; C\;v_0
Where \inline \rho is the density and \inline C the velocity of Sound for the fluid (Water)

Let the initial pipe velocity be \inline v_0

At the instant that the valve is closed a pressure wave, moving at the speed of sound \inline C sets off up the pipe bringing the water to rest.

Consider the instant shown above, \inline t secs after valve closure. At this moment the length of the column brought to rest is \inline x.

By Newton's second law, the pressure force equals the change of momentum per second.

pa = \frac{wax}{g}\times \frac{v_0}{t}

Pressure rise at valve
p = \frac{w}{g}\frac{x}{t}\;v_0;= \rho\,C\,v_0
C = \frac{x}{T}

  • The pressure at the valve remains at \inline p above normal from the instant the valve is closed until the pressure wave is reflected from the open end of the pipe as a wave of normal pressure and velocity \inline (= - v_0) reaching the valve at time \inline \displaystyle\frac{2L}{C}.
  • A wave of reduced pressure, \inline p below normal will then set off up the pipe.
  • The time \inline \displaystyle\frac{2L}{C} is the period of the pipe.
  • The above proof still applies if the valve is closed in a time of closure \inline T\leq\;\displaystyle\frac{2L}{C} i.e. the valve is closed before the pressure wave returns to the valve.

Pressure Rise For Instantaneous Partial Reduction Of Flow.


The Pressure rise
\delta\,p = \rho\,C\,\delta \,v
Where \inline \delta\;v is the change of velocity.
The increase in Pressure head \inline \delta\,h = \displaystyle\frac{\delta\,p}{w} =  \displaystyle\frac{\delta\,p}{\rho} =  \displaystyle\frac{C}{g}\delta\,v

The initial velocity \inline v_0 is reduced to \inline v_1 by the sudden partial valve closure. \inline t seconds later the pressure wave will have traveled a distance \inline x up the pipe as before \inline C = \frac{x}{t}

Applying Newton's second Law:

p = \frac{w}{g}\:\frac{x}{t}\,(v_0 - v_1) = \rho\,C(v_0 - v_1)

Which is commonly written as:

Pressure rise
\delta\,p = \rho\,C\delta\,v
Where \inline \delta\;v is the change of velocity.


The increase in Pressure head
\delta\;h = \frac{\delta\;p}{w} = \frac{\delta\;p}{\rho} = \frac{C}{g}\delta\;v

Instantaneous Partial Closure Of Valve

If the reduction in area of the valve is known rather than the reduction of flow or pipe velocity,then the pipe may be treated as a nozzle with a coefficient of discharge of \inline C_d

From the Initial Conditions and from the change in Valve Area:

The change of pressure head
\delta\,h = \frac{C}{g}\:\delta\,v

The initial flow through the valve
= Q_0 a_p v_0=C_d\:a_{v_0}\:\sqrt[]{2\,g\,h_0}


After partial closure:

\frac{Q_1}{Q_0} = \frac{v_1}{v_0} = \frac{a_v_1}{a_v_0}\sqrt[]{\frac{h_0 + \Delta\,h}{h_0}}

v_1 = v_0 - \Delta\,h

\Delta\,h = \frac{C}{g}\:\Delta\,v

For known initial conditions and change in the valve area this equation can be solved as a quadratic in \inline \delta\;v

\delta\,h = \frac{C}{g}\:\delta\,v

The Relationship Between Pressure Rise, Speed Of Sound And Bulk Modulus Allowing For Lateral Expansion Of The The Pipe


Bulk modulus \inline K = Increase of Pressure/Volumetric strain \inline = dp\;\div\frac{dv}{v}

  • Assume a thin cylinder constrained longitudinally and subject to an internal
pressure rise \inline p causing lateral expansion from diameter \inline  d\;\text{ to}\; d_1 And the equivalent Bulk Modulus \inline (K') allowing for Pipe Expansion is given by :


Hoop stress/Hoop strain = Young's Modulus

Hoop strain

p\,d\,l}\, = {2\;X\;f_H\;X}\;t\,l}

f_H = \frac{p\,d}{2\,t}
where \inline t is the thickness of the pipe wall


Consider the situation in the pipe at a time \inline T sec.after an instantaneous valve closure. The pressure wave will have reached a point \inline P at a distance \inline x from the valve where \inline C = \displaystyle\frac{x}{T} The water in the length \inline x is now at rest and at a pressure \inline p above normal.

During time \inline T after closure of the valve a volume \inline  v_0\, a \,T has passed point \inline P and this volume must be accommodated in addition to the initial volume \inline (a x). This is possible due to the compression of the water and the expansion of the pipe.

Assuming that \inline 2d\approx(d_1+d) and substituting from equation 29

v_0\,a\,T = \frac{p}{k}\:a\,x + \frac{\pi}{4}\:2\,d\:\frac{p\:d^2}{2\,t\,E}

Combining the above equation with \inline p = \rho\,v_0\,C

From the above two equations we can :

Eliminate \inline C by putting
\;p = \rho\,v_0X\,\frac{v_0}{p\left(\displaystyle\frac{1}{k}\:+\:\displaystyle\frac{d}{t\,E} \right)}

And the increase in pressure head
= \frac{p}{w}

To eliminate \inline p:

By ignoring Pipe Expansion Equations (35 ) (38 ) (41) and (46 ) can be simplified so that:
p = v_0\;\sqrt[]{K \rho}
And The head rise

C = \sqrt[]{\frac{K}{\rho}}

Thus the equivalent bulk modulus allowing for pipe expansion is given by:

Pipe With Change Of Section


At the initial steady flow conditions it is assumed that friction; velocity head and contraction losses can be neglected. Therefore the pressure throughout both pipes is \inline h_0.


The above diagram shows some instant when the pressure wave is between \inline V & P

from equation (46)

Flow through valve
=Q_1 = a\,v_1 = Cd\,a_{v1}\:\sqrt[]{2g\,h_1}

And from equation (16)

h_1 - h_0 = \Delta\,h = \frac{C}{g}\,\Delta = C(v_0 - v_1)

Equations ( 54 ) ( 55 ) and ( 56 ) can be solved for \inline \delta\;v and hence \inline v_1 and \inline \delta\;h and hence \inline h_1.

Now consider an instant in time after the pressure wave has passed through the junction and a reflected wave has set off up the pipe from \inline P to \inline R.


By continuity:

Using equation ( 16 ) for the pressure wave in the large pipe:

Similarly for the pressure wave in the small pipe:

Note: \inline C_* is the velocity of sound in the small pipe.

Adding equations ( 59 ) and(60 )

Equate to equation ( 56 )

substitute from ( ) and ( ) for V_0 & V_1

Using the value for \inline v_1 found from equations ( 14 ) ( 27 ) and( ) this equation can be solved for \inline  v_2 and hence \inline  V_1 and \inline  h_2.


  • If the velocity of sound is the same for both pipes then \inline  C = C_*

  • Putting \inline A = infinity i.e.the pipe entering the reservoir

Where \inline \delta\,v is the initial reduction of the pipe velocity due to a partial closure of the valve. i.e. When the pressure wave reaches the reservoir entrance the velocity is reduce still further by an amount equal to the initial velocity reduction at the valve.

The solution to the first worked example has been hidden. It can be seen by clicking on the red button. The solution to the other examples can only be seen by Members of Codecogs. Membership is free and registering is quick and easy.

Example - Example 1
A pipe carrying water has a valve at the discharge end. The bulk modulus,\inline K, is \inline 2\times 10^5\;lb./in.^2
  • a) If the initial velocity in a pipeline is \inline V_0 show that an instantaneous closure of the valve will result in a pressure rise at the valve of \inline 7470\;v_0lb/sq,ft.
  • b) Given that the initial velocity is 12 ft/sec. and the pressure head in the pipe is 600 ft, find the increase in pressure head which results when the valve is instantaneously closed by 1/5 th.
a) Using equation (39)

b) When the valve is fully open using equation (42)

And when the valve is closed by 1/5 th

Combining the above two equations and putting \inline h_1=h_0 + \delta\;h

1\:- \frac{\delta\,V}{V} = \frac{a_{v_1}}{a_{v_0}}\sqrt{1 + \frac{\delta\,h}{h_0}}

Solving the quadratic in \inline \delta\,v

The increase in the pressure head
=\frac{\delta v}{g}\times \sqrt{\frac{K}{\rho}}=151\;ft.\text{of water}

The increase in the pressure head is \inline 151\;ft.\text{of water}

Last Modified: 10 Jan 12 @ 17:02     Page Rendered: 2022-03-14 15:56:28