• https://me.yahoo.com

Compressible Flow

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

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\;&space;C\;v_0$
Where $\inline&space;\rho$ is the density and $\inline&space;C$ the velocity of Sound for the fluid (Water)

Let the initial pipe velocity be $\inline&space;v_0$

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

Consider the instant shown above, $\inline&space;t$ secs after valve closure. At this moment the length of the column brought to rest is $\inline&space;x$.
$\therefore\;\;\;C&space;=&space;\frac{x}{t}$

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

i.e.
$pa&space;=&space;\frac{wax}{g}\times&space;\frac{v_0}{t}$

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

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

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

Applying Newton's second Law:
$p\:a&space;=&space;\frac{w\,a\,x}{g}\;\left(\frac{v_0&space;-&space;v_1}{t}&space;\right)$

Pressure
$p&space;=&space;\frac{w}{g}\:\frac{x}{t}\,(v_0&space;-&space;v_1)&space;=&space;\rho\,C(v_0&space;-&space;v_1)$

Which is commonly written as:

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

And:

$\delta\;h&space;=&space;\frac{\delta\;p}{w}&space;=&space;\frac{\delta\;p}{\rho}&space;=&space;\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&space;C_d$

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

$\delta\,h&space;=&space;\frac{C}{g}\:\delta\,v$

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

After partial closure:

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

Thus:
$\frac{Q_1}{Q_0}&space;=&space;\frac{v_1}{v_0}&space;=&space;\frac{a_v_1}{a_v_0}\sqrt[]{\frac{h_0&space;+&space;\Delta\,h}{h_0}}$

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

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

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

Hence:
$\delta\,h&space;=&space;\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 $\inline&space;K$ = Increase of Pressure/Volumetric strain $\inline&space;=&space;dp\;\div\frac{dv}{v}$

• Assume a thin cylinder constrained longitudinally and subject to an internal
pressure rise $\inline&space;p$ causing lateral expansion from diameter $\inline&space;&space;d\;\text{&space;to}\;&space;d_1$
$C&space;=&space;\sqrt{\frac{1}{\rho\left(&space;\displaystyle\frac{1}{K}\:+\:&space;\displaystyle\frac{d}{t\,E}&space;\right)}}$
And the equivalent Bulk Modulus $\inline&space;(K')$ allowing for Pipe Expansion is given by :
$\frac{1}{K_1}&space;=&space;\frac{1}{K}\:+\:\frac{d}{t\,E}$

Hoop stress/Hoop strain = Young's Modulus

Hoop strain

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

Also:
$p\,d\,l}\,&space;=&space;{2\;X\;f_H\;X}\;t\,l}$

Or:
$f_H&space;=&space;\frac{p\,d}{2\,t}$
where $\inline&space;t$ is the thickness of the pipe wall

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

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

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

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

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

Then:
$v_0\,a\,T&space;=&space;\frac{p}{k}\:a\,x&space;+&space;\frac{\pi}{4}\:2\,d\:\frac{p\:d^2}{2\,t\,E}$
$\therefore\;\;\;\;v_0\,T&space;=&space;p\,x\left(\frac{1}{K}&space;+&space;\frac{d}{t\,E}&space;\right)$

Combining the above equation with $\inline&space;p&space;=&space;\rho\,v_0\,C$

From the above two equations we can :

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

And the increase in pressure head
$=&space;\frac{p}{w}$
or:
$\frac{p}{\rho\,g}$

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

By ignoring Pipe Expansion Equations (35 ) (38 ) (41) and (46 ) can be simplified so that:
$C&space;=&space;\frac{v_0&space;K}{p}$
$p&space;=&space;v_0\;\sqrt[]{K&space;\rho}$
$=\frac{v_0}{g}\sqrt{\frac{K}{\rho}}$

And:
$C&space;=&space;\sqrt[]{\frac{K}{\rho}}$

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

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

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&space;h_0$.

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

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

from equation (46)

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

And from equation (16)

also:
$h_1&space;-&space;h_0&space;=&space;\Delta\,h&space;=&space;\frac{C}{g}\,\Delta&space;=&space;C(v_0&space;-&space;v_1)$

Equations ( 54 ) ( 55 ) and ( 56 ) can be solved for $\inline&space;\delta\;v$ and hence $\inline&space;v_1$ and $\inline&space;\delta\;h$ and hence $\inline&space;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&space;P$ to $\inline&space;R$.

By continuity:
$AV_1=av_2$

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

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

Similarly for the pressure wave in the small pipe:

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

Note: $\inline&space;C_*$ is the velocity of sound in the small pipe.

Adding equations ( 59 ) and(60 )

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

Equate to equation ( 56 )

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

substitute from ( ) and ( ) for V_0 & V_1

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

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

Notes:

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

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

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

• Putting $\inline&space;A$ = infinity i.e.the pipe entering the reservoir

$v_2&space;=&space;2v_1&space;-&space;v_0&space;=&space;2v_1&space;-&space;2v_0&space;+&space;v_0$

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

Where $\inline&space;\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,$\inline&space;K$, is $\inline&space;2\times&space;10^5\;lb./in.^2$
• a) If the initial velocity in a pipeline is $\inline&space;V_0$ show that an instantaneous closure of the valve will result in a pressure rise at the valve of $\inline&space;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.
Workings
a) Using equation (39)

$p&space;=&space;V_0\;\sqrt[]{K\,\rho}$

$\therefore\;\;\;p&space;=&space;V_0\;\sqrt[]{2\times&space;10^5\times&space;12^2\times&space;\frac{62.4}{32.2}}$

$=&space;7470\;V_0\;lb/sq.ft$

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

$Q_0&space;=&space;a_p&space;V_0&space;=&space;Cd\;a_{v_0}\,\;\sqrt[]{2g\,h_0}$

And when the valve is closed by 1/5 th

$Q_1&space;=&space;a_p\;V_1&space;=&space;Cd\;a_{v_1}\;\;\sqrt[]{2g\:h_1}$

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

$\therefore\;\;\;\frac{V_0&space;-&space;\delta\,V}{V_0}&space;=&space;\frac{a_{v_1}}{a_{v_0}}\;\sqrt[]{\frac{h_0&space;+&space;\delta\,h}{h_0}}$

Thus:
$1\:-&space;\frac{\delta\,V}{V}&space;=&space;\frac{a_{v_1}}{a_{v_0}}\sqrt{1&space;+&space;\frac{\delta\,h}{h_0}}$

$=&space;\frac{a_{v_1}}{a_{v_0}}\;\sqrt[]{1&space;+&space;\frac{\delta\,v}{g}\;X\;\frac{1}{h_0}\;\sqrt[]{\frac{K}{\rho}}}$

$1&space;-&space;\frac{\delta\,v}{12}&space;=&space;\frac{4}{5}\;\sqrt[]{1&space;+&space;\frac{\delta\,v}{32.2\;X\;600}\;\sqrt[]{\frac{K}{\rho}}}$

Solving the quadratic in $\inline&space;\delta\,v$

$\,v&space;=&space;1.26\:ft/sec$

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

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