Engineering › Fluid Mechanics › Pipes › Water Hammer ›

Incompressible Flow

An anlysis of Water Hammer based on the assumption that water is incompressible

View other versions (6)

Contents

Overview
The Gradual Valve Closure
Valve Closed To Produce Constant Acceleration
Valve Closed To Produce Uniform Rate Of Pressure Increase
The Valve Closes So That The Area Decreases Uniformly With Time
Page Comments

Key Facts

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

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

Overview

Using The Rigid Column Theory the compressibility of the fluid is ignored and it is assumed that pressure changes caused by opening or closing a valve are felt instantaneously through out the pipe. In effect the water column is a solid column, which can accelerate or decelerated as an entity.

The Gradual Valve Closure

A valve is a device that regulates, directs or controls the flow of a fluid by opening, closing, or partially obstructing various passageways.

As the valve closes the pressure at the valve rises decelerating the water column.

Let $\inline p$ be the pressure rise above the initial steady flow pressure, $\inline \alpha$ the deceleration of the water column, $\inline p$ the rise in pressure, $\inline a$ the cross sectional area and $\inline l$ the length of the pipe.

MISSING IMAGE!

23287/Water-Hammer-IF-001.png cannot be found in /users/23287/Water-Hammer-IF-001.png. Please contact the submission author.

The increase in Pressure Head $\inline H$ is given by:

$H=\frac{p}{w}=-\frac{l}{g}\times \alpha=-\frac{l}{g}\left ( \frac{dv}{dt} \right ) /) \begin{quote} By Newton's Second Law: Pressure at Valve=Mass of water multiplied to Deceleration \(\therefore\;\;\;pa=\frac{w\times a\times l}{g}\times\alpha$

$\therefore\;\;\;p=\frac{w\timesl}{g}$

i.e.

$H=\frac{p}{w}=-\frac{l}{g}\times\frac{dv}{dt}$

(2)

Command [\end] not found.

{quote}

Valve Closed To Produce Constant Acceleration

Acceleration is the rate of change of velocity over time.

$T=\frac{w}{g}\times \frac{l}{p}\times v_0$

Let the time taken to close the valve be $\inline T$ seconds.

Then:

$\alpha=\frac{v}{T}$

(3)

$\therefore\;\;\;\;p=\frac{wl}{g}\times\frac{v_0}{T}\times H=+\frac{l\times v_0}{gT}$

So for any given pipeline, if the maximum pressure rise is specified, the time $\inline T$ taken to close the valve is directly proportional to the initial water velocity.

$T=\frac{w}{g }\times \frac{l}{p}}\times v_0$

Valve Closed To Produce Uniform Rate Of Pressure Increase

Pressure is the force per unit area applied in a direction perpendicular to the surface of an object.

The maximum pressure rise

$= \frac{wl}{g}\times \alpha_{max}=\frac{wl}{g}\times KT$

Where $\inline K$ is a constant and $\inline T$ is the total time of valve closure.

At any instant during closure

$\alpha=-K \times t$

(4)

$\frac{dv}{dt}= -K\times t$

$\therefore\:\;\;v= -K\times t^2+C$

where $\inline C$ is Constant.
The values of $\inline K$ and $\inline C$ can be found from the known initial and final conditions $\inline v=v_0$ at $\inline t=0$ and $\inline v=0$ when $\inline t=T$ .

then:

$\alpha_{max}=K\times T_{max}$

And the maximum pressure rise

$=\frac{wl}{g}\times\alpha_{max}=\frac{w\times l}{g}\times K\:T$

The Valve Closes So That The Area Decreases Uniformly With Time

MISSING IMAGE!

23287/Water-Hammer-IF-002.png cannot be found in /users/23287/Water-Hammer-IF-002.png. Please contact the submission author.

$\left(\frac{h_m}{h_0} \right) = \frac{K}{2}\;+\:\sqrt[]{\frac{K^2}{4}\:+\:K}$

(5)

This equation again relates the maximum pressure rise directly to the time of valve closure. The difference between this and the first case is that the rate of increase in $\inline p$ has been controlled. The equation gives for a known pipe, the maximum pressure rise given any closure time. Note: The suffix "o" refers to the initial conditions when the head behind the valve and is equal to the gross supply head if the friction in the pipe is neglected and $\inline h$ is the Inertia head at a time $\inline t$ .

For Valve opening

$\frac{h_m}{h_0} = \frac{K}{2} - \sqrt[]{\frac{K^2}{4}\:+\:K}$

A nozzle is a device designed to control the direction or characteristics of a fluid flow as it exits an enclosed chamber or pipe via an orifice.

Consider the valve as a nozzle at the end of the pipe with a $\inline C_d$ which is constant.

Initial flow

$Q_0=av_o = C_d a_{v_{0}}\sqrt{2gh_0}$

At a time $\inline t$ the discharge

$Q = a v = C_d\:a_v\:\sqrt[]{2\,g\,(h_0\:-\:h)}$

$\therefore\;\;\;\;\frac{Q}{Q_0}\:= \frac{v}{v_0}\:= \frac{a_v}{a_{v_0}}}$

But:

$a_v = a_{v_0}}\times left(1-\frac{t}{T}) \right)$

$\therefore\;\;\;\;\frac{v}{v_0} = \left(1-\frac{t}{T} \right)\sqrt[]{1+\frac{h}{h_0}}$

$\therefore\;\;\;\;{v} = v_0\left(1-\frac{t}{T} \right)\sqrt[]{1+\frac{h}{h_0}}$

$\therefore\;\;\;\frac{dv}{dt}= - \frac{v_0}{T}\sqrt[]{1+\frac{h}{h_0}}+ v_0\left(1-\frac{t}{T} \right)\frac{1}{2}\left(1+\frac{h}{h_0} \right)}^{- \frac{1}{2}}\frac{dh}{dt}$

When $\inline h$ is a maximum:

$\frac{dh}{dt}\ =\ 0$

$\therefore\;\;\;\;\;\left\(\frac{dv}{dt} \right)_{h\,=\,h_{max}}=-\frac{v_0}{T}\sqrt[]{1+\frac{h}{h_0}}$

The Maximum head

$h_m = - \left(\frac{dv}{dt} \right)_{h=h_{max}} = \frac{l}{g}\frac{v_0}{T}\sqrt[]{1+\frac{h_m}{h_0}}$

$\therefore\;\;\;\frac{h_m}{h_0}= \frac{l}{g}\frac{v_0}{Th_0}\sqrt[]{1+\frac{h_m}{h_0}}$

Putting:

$K = \left( \frac{l}{g}\frac{v_0}{Th_0} \right)^2$

$\therefore\;\;\;\left(\frac{h_m}{h_0} \right)^2 = K\left(1 + \frac{h_m}{h_0} \right)$

(6)

$\therefore\;\;\;\left(\frac{h_m}{h_0} \right)^2 - K\left(\frac{h_m}{h_0} \right)-K = 0$

Solving the quadratic:

$\left(\frac{h_m}{h_0} \right) = \frac{K}{2}+\sqrt[]{\frac{K^2}{4}+K}$

(7)

Which must be the positive root for valve closure.

For valve opening.

$\frac{h_m}{h_0} = \frac{K}{2} - \sqrt[]{\frac{K^2}{4}\:+\:K}$

MISSING IMAGE!

23287/Water-Hammer-IF-003.png cannot be found in /users/23287/Water-Hammer-IF-003.png. Please contact the submission author.

For pipe 1

$h_1=-\frac{1}{g}\left ( l_1+\frac{a_1}{a_2}l_2 \right )\;\frac{dv_1}{dt}$

Note: The above expression is the same as equation (2) with an "equivalent" pipe length of:

$l= l_1+\frac{a_1}{a_2}l_2$

And for a number of pipe sections of $\inline n$ different cross sectional areas:

$l= l_1+\frac{a_1}{a_2}l_2+\frac{a_1}{a_3}l_3+.........\frac{a_1}{a_n}l_n$

These equations allow the calculation of an "overall length" which can be used to calculate pressure rises under differing valve opening regimes.

Applying Newton's Second Law for each pipe:

MISSING IMAGE!

23287/Water-Hammer-IF-004.png cannot be found in /users/23287/Water-Hammer-IF-004.png. Please contact the submission author.

$wh_0+wh_2a_2-wh_0a_2=\frac{wa_2l_2}{g}(-\frac{dv_2}{dt})$

$\therefore\;\;\;\;\;h_2 = -\frac{l_2}{g}\,\frac{dv_2}{dt}$

(8)

MISSING IMAGE!

23287/Water-Hammer-IF-005.png cannot be found in /users/23287/Water-Hammer-IF-005.png. Please contact the submission author.

$\therefore\;\;\;\;\;w(h_0+h_1)a_1 w(h_0-h_2)\times a_1$

$=\frac{wa_1l_1}{g}\times \left (- \frac{dv_1}{dt} \right )$

$\therefore\;\;\;\;\;h_1-h_2 = -\frac{l}{g}\times \frac{dv_1}{dt}$

(9)

By continuity:

$a_1v_1 = a_2v_2$

$\therefore\;\;\;\;\;a_1\times \frac{dv_1}{dt} = a_2\frac{dv_2}{dt}$

(10)

Putting equation (10) into equation (8)

$h_2 = -\frac{l_2}{g}\times \frac{a_1}{a_2}\times \frac{dv_1}{dt}$

(11)

Equation (11) can now be inserted into equation (9)

$\therefore\;\;\;\;\;h_1 -\left ( -\frac{l_1}{g}\times \frac{dv_1}{dt} \right ) -\frac{l_2}{g}\times \frac{a_1}{a_2}\times \frac{dv_1}{dt}$

$= -\:\frac{1}{g}\:\left(l_1\:+\:\frac{a_1}{a_2}\:l_2 \right)\:\frac{dv_1}{dt}$

Note: The above expression is the same as equation (2) with an "equivalent" pipe length of:

$l= l_1+\frac{a_1}{a_2}l_2$

And for a number of pipe sections of $\inline n$ different cross sectional areas:

$l= l_1+\frac{a_1}{a_2}l_2+\frac{a_1}{a_3}l_3+.........\frac{a_1}{a_n}l_n$

These equations allow the calculation of an "overall length" which can be used to calculate pressure rises under differing valve opening regimes.

N.B. The following two examples are based on the assumption that water is not compressible.

Example:

[imperial]

Example - Example 1

Problem

Find the maximum increase in pressure head if the valve closes in 12 seconds.

a) So that the velocity decreases uniformly with time.
b) So that the pressure rises uniformly with time.
c) So that the valve discharge area decreases uniformly with time.
d) The last 1000 ft. of pipe is 8 ft. diam. and the velocity decreases uniformly with time.

Workings

Initial velocity $V_0$

$= \frac{1100}{\displaystyle\frac{\pi}{4}\times 100} = 14.01\:ft./sec.$

a) Velocity decreases uniformly with time.

Using equation (2) The Inertia Head

$h = \frac{l}{g}\frac{v_0}{T}$

$= \frac{3000}{32.2}\times \frac{14.01}{12} = 108.8\:ft.$

b) Pressure rises uniformly with time.

Using equation (3)

$v= - \frac{K}{2}\:t^2\:+\:C$

When $T = 0$ $v = v_0$ which equals 14.01 ft./sec.

When $t$ = 12sec. $v = 0$

$\therefore\;\;\;\;\;0 = 1\frac{K}{2}\times 144+14.01$

$\therefore\;\;\;\;\;K = \frac{14.01}{72}$

Thus:

$\alpha= \frac{14.01}{72}\times t$

And:

$\alpha_{max}=\frac{14.01\times 12}{72} = \frac{14.01}{6}\:ft,/sec^2$

But from equation (4)

$h_{max} = \frac{l}{g}\;\alpha_{max}$

$= \frac{3000}{32.2}\;X\;\frac{14.01}{6} = 217.6\:ft.$

c) Valve Discharge Area decreases Uniformly with time

Using equation (5)

$\frac{h_m}{h_0} = \frac{K}{2} + \sqrt[]{\frac{K^2}{4} + K}$

Where:

$K = \left(\frac{lv_0}{gTh_0} \right)^2 = \left(\frac{3000\times 14.01}{32.2\times 12\times 500} \right)^2 = \:0.04731$

Thus from equation (4)

$\frac{h_m}{h_0} = \frac{0.04731}{2} = \sqrt[]{\frac{0.04731^2}{4}+0.04731}=0.2425$

Thus:

$h_m = 0.2425 \times 500 = 121.3\; ft.$

d) The last 1000 ft. of pipe is 8 ft. diam. and the velocity decreases uniformly with time.

Using equation (6) The initial velocity in the 8" pipe

$= \left(\frac{10}{8} \right)^2\times 14.01=$

$= 21.9\:ft/sec.$

And from equation (7) The Equivalent Length

$l = l_1 + \frac{a_1}{a_2}\:l_2 = 1000+\left(\frac{8}{10} \right)^2\times 2000 =$

$=2280\. ft.$

From equation (2) The Pressure on the Valve

$= \frac{l}{g}\;\frac{v_{1.0}}{T} = \frac{2280}{32.2}\times \frac{21.9}{12} = 129.2\: ft.$

And the Pressure rise at the junction

$h_1 = \frac{l_2}{g}\times\frac{v_2}{T}$

$=\frac{2000}{32.2}\times \frac{14.01}{12} = 72.5\:ft.$

Notes:

1. Pipe friction is small compared to the pressure rise and has been neglected
2. The time of closure $T$ is greater than $\displaystyle\frac{L}{1000}$

Solution

a) When Velocity decreases uniformly with time the Inertia Head

is $108.8\;ft.$

b) When Pressure rises uniformly with time, $217.6\;ft.$
c) When Valve Discharge Area decreases Uniformly with time, $121.3\; ft.$
d) When the last $1000 ft.$ of pipe is $8 ft.$ diam. and the velocity decreases uniformly with time, $72.5\;ft.$