Engineering › Fluid Mechanics › Pipes › Water Hammer ›

Compressible Flow

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

View other versions (7)

Contents

Overview
Pressure Rise Following Instantaneous Closure
Pressure Rise For Instantaneous Partial Reduction Of Flow.
Instantaneous Partial Closure Of Valve
The Relationship Between Pressure Rise, Speed Of Sound And Bulk Modulus Allowing For Lateral Expansion Of The The Pipe
Pipe With Change Of Section
Page Comments

Key Facts

Gyroscopic Couple: The rate of change of angular momentum ( $\tau$ ) = $I\omega\Omega$ (In the limit).

= Moment of Inertia.
$\omega$ = Angular velocity
$\Omega$ = Angular velocity of precession.

Overview

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 $\rho$ is the density and

the velocity of Sound for the fluid (Water)

Let the initial pipe velocity be

At the instant that the valve is closed a pressure wave, moving at the speed of sound

sets off up the pipe bringing the water to rest.

Consider the instant shown above,

secs after valve closure. At this moment the length of the column brought to rest is

$\therefore\;\;\;C = \frac{x}{t}$

(2)

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

i.e.

$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$

where

$C = \frac{x}{T}$

Notes:

The pressure at the valve remains at 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 reaching the valve at time $\displaystyle\frac{2L}{C}$ .
A wave of reduced pressure, below normal will then set off up the pipe.
The time $\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 $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 $\delta\;v$ is the change of velocity.
The increase in Pressure head $\delta\,h = \displaystyle\frac{\delta\,p}{w} = \displaystyle\frac{\delta\,p}{\rho} = \displaystyle\frac{C}{g}\delta\,v$

The initial velocity

is reduced to

by the sudden partial valve closure.

seconds later the pressure wave will have traveled a distance

up the pipe as before $C = \frac{x}{t}$

Applying Newton's second Law:

$p\:a = \frac{w\,a\,x}{g}\;\left(\frac{v_0 - v_1}{t} \right)$

(3)

Pressure

$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 $\delta\;v$ is the change of velocity.

And:

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

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:

$Q_1 = a_p\;v_1 = C_d\;a_v_1\;\sqrt[]{2g(h_0 - \Delta\,h)}$

(4)

Thus:

$\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}}$

But:

$v_1 = v_0 - \Delta\,h$

$\therefore\;\;\;\frac{v_0 - \Delta\,v}{v_0} = \frac{a_v_1}{a_v_0}\sqrt[]{1 + \frac{\Delta\,h}{h_0}}$

(5)

But:

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

$\therefore\;\;\;1 - \frac{\Delta\,v}{v_0} = \frac{a_v_1}{a_v_0}\sqrt[]{1\;+\frac{C\;\Delta\,v}{g\:h_0}}$

(6)

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

Hence:

$\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

Notes:

Bulk modulus

= Increase of Pressure/Volumetric strain $= dp\;\div\frac{dv}{v}$

Assume a thin cylinder constrained longitudinally and subject to an internal

pressure rise

causing lateral expansion from diameter $d\;\text{ to}\; d_1$

$C = \sqrt{\frac{1}{\rho\left( \displaystyle\frac{1}{K}\:+\: \displaystyle\frac{d}{t\,E} \right)}}$

(7)

And the equivalent Bulk Modulus

allowing for Pipe Expansion is given by :

$\frac{1}{K_1} = \frac{1}{K}\:+\:\frac{d}{t\,E}$

(8)

Hoop stress/Hoop strain = Young's Modulus

Hoop strain

$=\frac{f_H}{E} = \frac{\pi(d_1\:-\:d)}{\pi\,d} = \frac{d_1 - d}{d}$

(9)

Also:

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

Or:

$f_H = \frac{p\,d}{2\,t}$

where

is the thickness of the pipe wall

$\therefore\;\;\;\frac{p\,d}{2\,t\,E} = \frac{d_1 - d}{d}$

(10)

Consider the situation in the pipe at a time

sec.after an instantaneous valve closure. The pressure wave will have reached a point

at a distance

from the valve where $C = \displaystyle\frac{x}{T}$ The water in the length

is now at rest and at a pressure

above normal.

During time

after closure of the valve a volume $v_0\, a \,T$ has passed point

and this volume must be accommodated in addition to the initial volume

. This is possible due to the compression of the water and the expansion of the pipe.

$v_0\,a\,T = \frac{p}{k}\;X\;a\:x + \frac{\pi}{4}\left(d_1^2 - D^2 \right)x}$

(11)

$= \frac{p}{K}\;X\;ax + \frac{\pi}{4}\:(d_1 + d)(d_1 - d)$

(12)

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

Then:

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

$\therefore\;\;\;\;v_0\,T = p\,x\left(\frac{1}{K} + \frac{d}{t\,E} \right)$

(13)

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

From the above two equations we can :

Eliminate

by putting

$\;p = \rho\,v_0X\,\frac{v_0}{p\left(\displaystyle\frac{1}{k}\:+\:\displaystyle\frac{d}{t\,E} \right)}$

$\therefore\;\;\;p = v_0\;\sqrt{\frac{\rho}{\left(\displaystyle\frac{1}{K}\:+\:\displaystyle\frac{d}{t\,E} \right)}}$

(14)

And the increase in pressure head

$= \frac{p}{w}$

or:

$\frac{p}{\rho\,g}$

To eliminate

$C = \frac{v_0}{\left(\displaystyle\frac{1}{K}\:+\:\displaystyle\frac{d}{t\,E} \right)\rho\:v_0\:C}$

(15)

$\therefore\;\;\;C = \sqrt{\frac{1}{\rho\left(\displaystyle\frac{1}{K}\:+\:\displaystyle\frac{d}{t\,E} \right)}}$

(16)

By ignoring Pipe Expansion Equations (35 ) (38 ) (41) and (46 ) can be simplified so that:

$C = \frac{v_0 K}{p}$

(17)

$p = v_0\;\sqrt[]{K \rho}$

And The head rise

$=\frac{v_0}{g}\sqrt{\frac{K}{\rho}}$

And:

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

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

$\frac{1}{K_1} = \frac{1}{K}\:+\:\frac{d}{t\,E}$

(18)

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

$\therefore\;\;\;Q_0 = A\,V_0\:=\:a\,v_0\:=\:Cd\,a_v_0\:\sqrt[]{2g\,h_0}$

(19)

The above diagram shows some instant when the pressure wave is between

from equation (46)

Flow through valve

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

And from equation (16)

also:

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

Equations ( 54 ) ( 55 ) and ( 56 ) can be solved for $\delta\;v$ and hence

and $\delta\;h$ and hence

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

By continuity:

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

$h_2 - h_0 = \frac{C}{g}\:(V_0 - V_1)$

(20)

Similarly for the pressure wave in the small pipe:

$h_1 - h_0 = \frac{C_*}{g}\;(v_1 - v_2)$

(21)

Note:

is the velocity of sound in the small pipe.

Adding equations ( 59 ) and(60 )

$h_1 - h_0 = \frac{C}{g}(V_0 - V_1) + \frac{C_*}{g}}\;(v_1 - v_2)$

(22)

Equate to equation ( 56 )

$\frac{C}{g}(V_0 - V_1) + \frac{C_*}{g}(v_1 - v_2) = \frac{C_*}{g}(v_0 - v_1)}$

(23)

substitute from ( ) and ( ) for V_0 & V_1

$\frac{C}{g}\left(\frac{a}{A}(v_0 - v_2) \right) = \frac{C_*}{g}(v_0 - v_1) - \frac{C_*}{g}(v_0 - v_2)$

(24)

Using the value for

found from equations ( 14 ) ( 27 ) and( ) this equation can be solved for

and hence

and

Notes:

If the velocity of sound is the same for both pipes then

$\frac{a}{A}(v_0 - v_2) = v_0 - 2v_1 + v_2$

(25)

$\therefore\;\;\;\left(1 + \frac{a}{A} \right) = 2v_1 - v_0\left(1 - \frac{a}{A} \right)$

(26)

Putting = infinity i.e.the pipe entering the reservoir

(27)

$\therefore\;\;\;v_2 = v_0 - 2(v_0 - v_1) = v_0 - 2\,\Delta\,v$

(28)

Where $\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:

[imperial]

Example - Example 1

Problem

A pipe carrying water has a valve at the discharge end. The bulk modulus, $K$ , is $2\times 10^5\;lb./in.^2$

a) If the initial velocity in a pipeline is $V_0$ show that an instantaneous closure of the valve will result in a pressure rise at the valve of $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.