Chapter 5
Chapter 5
Dynamics
CHAPTER 5
FLUID DYNAMICS
5.1 Introduction
In discussing about hydrostatics we were concerned with forces (pressure forces) which
are acting on an object for a liquid at rest and when we deal with kinematics of fluid flow
phenomena related with space time variation (velocity and acceleration) with out
considering the effect of force. But in dealing with dynamics of fluid flow, all forces that
affect the phenomenon are considered.
The dynamics of fluid motion deals with kinetics, which relates the kinematics with the
forces responsible for causing the motion. This relation ship of fluid motion is established
by the use of laws of nature.
i) The principle of conservation of mass (the continuity relation ship)
ii) Newton’s laws of motion
iii) The 1st and second lows of thermodynamics
The forces which may be present in fluid flow problems are gravity forces F g, pressure
force Fp, force due to viscosity Fv, force due to turbulence Ft, Surface tension Fs, and
force due to compressibility of fluid Fc.
Gravity forces (Fg.) is due to the weight of the fluid. Its component in flow direction
results in acceleration.
Pressure force (Fp): It is equal to the product of pressure intensity and cross sectional area
of the flowing fluid. Acts normal to the surface under consideration and produces
acceleration in the given direction.
Viscous forces (Fv): - Exists in real fluids. It is the shearing resistance generated when
there is relative motion between two layers of fluids. It acts opposite to the direction of
motion, and retards flow.
Surface tension (Fs): This force is important when the depth of flow is extremely small.
Force due to compressibility (Fc): for incompressible fluids, this becomes significant in
problems of unsteady flow like water hammer. In most of flow problems, F s and Fc are
neglected.
1
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
Force due to turbulence (Ft): the continuous momentum transfer between layers in highly
turbulent flow results in normal and shear stresses known as Reynolds's stress.
If the changes for the change in forces are small the forces can be taken negligible.
max = (Fg)x+(Fp)x+(Fv)x+(Ft)x
The presence of such a complex system of forces in real fluid flow problems makes the
analysis very complicated. Therefore, mathematical analysis of problems is generally
possible only if certain simplifying assumptions are made.
1. Potential energy /elevation /positional energy/ because of its elevation above datum
level. If a weight w of liquid is at a height of z above datum
Potential energy = Wz
Potential energy per unit weight = z (meters) = Potential head
2. Pressure energy: When a fluid flows in a continuous stream under pressure it can do
work. If the area of cross – section of the stream of fluid is a, then force due to
pressure p on cross- section is Pa.
W = F* S =
3. Kinetic energy
If a mass of fluid (m) moves at some velocity (v),
Kinetic energy = ½ mv2 = ½ W/g v2
Kinetic energy pr unit weight = = kinetic head
2
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
Bernoulli’s assumptions:
dP
[ P ( )ds ]dA
ds ds The flow is Steady an
incompressible type of flow,
Frictionless and non viscous flow,
No shaft work - no pump or
PdA turbines on the streamline,
No transfer of heat - either added
Z or removed.
dsdA
Fig: 5.1 Pressure and gravity forces on a cylindrical element along a streamline
The normal forces on the side faces are in equilibrium and as the fluid is assumed non-
viscous, there is no shear stress. The velocity varies along the streamline and there is
acceleration. It is necessary to take into account force due to acceleration when
considering the longitudinal balance of force. Z
But in the ease of steady flow the velocity doesn’t vary at a point so that local
acceleration will be zero but for velocity variation with position convective
acceleration will be different from zero . The forces tending to accelerate the
fluid mass are pressure force on the two ends of the element,
3
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
For the case of an incompressible fluid may be treated as constant, the integration gives
[Bernoulli’s Equation]
Under special conditions the assumption underlying Bernoulli’s equations can be waived.
The Bernoulli equation is the basis for the solution of a wide range of hydraulics
problems. For two points along a streamline, the Bernoulli equation may be expressed in
the form of:
The increase in temperature of the fluid causes an increase in the internal energy. The
increase in internal energy and the heat transfer from the fluid represent a loss of useful
energy. The total loss per unit mass of fluid is (u2-u1-q).
Energy loss per unit weight in over coming resistance (head loss)
4
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
The total energy of flow decreases in the flow direction, and consequently the energy line
has a down ward slope.
Horizontal
h1
V12/2g The modified Bernoulli’s equation for up stream
V2 /2g
2 section (1) and downstream section (2)
P1/
P2/
1
2 hl = head loss between the two sections.
Z1
Datum Z2
2
v1
EGL HL = head loss
2g 2
p1 HGL v2
2g
p2
Z1
Z2
Datum
The analysis of flow problems is usually based on the one-dimensional approach. The
entire flow is considered to be taking through stream tube with average velocity V at
center of each cross-section. The velocity distribution at any x-section in real fluid flow is
non uniform, on account of the boundary resistance and consequently the kinetic energy
per unit weight given by V2/2g doesn’t represent the kinetic energy across the section.
The kinetic energy per unit time passing through on elemental area dA is ½ (dAu)u2
u - velocity at that point
Total kinetic energy passing the section
And the actual kinetic energy passed on average velocity V passing
Laminar flow is purely a viscous flow; the value of is maximum and equals 2.0. But in
case of fully developed turbulent flow in pipes, is independent of Reynolds number and
may be considered to have almost constant value (1.01 to 1.15) depending on surface
roughness and Reynolds number. Lower value is appreciable for velocity rough surface
and high Reynolds number.
For an identical velocity distribution at two sections, 1 = 2, and if accuracy is not
required in non-uniform velocity distribution = 1.
It is a device used for measuring rate of flow in a pipeline and it consists of three
components:
i) A converging entrance cone of angle of about 200
ii) A cylindrical portion of short length called the “throat”
6
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
The entrance tube and exit tube diameter are the same as that of the pipe line in to which
it is inserted and the length of throat is equal to the throat diameter.
7
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
Theoretical Discharge
For vertical or inclined Venturi meter, the actual discharge can be computed similarly.
Pitot tube is a device used for measuring velocity of flow at any point in a pipe or a
channel. In its elementary form a Pitot tube consists of an L-shaped tube with open ends.
Its may be aligned in open channel or pipe flow measurement as indicated below.
8
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
Stagnation pressure
For a figure below it has been seen that the central streamline terminates
at B the entrance to the Pitot tube. This is on account of the inability of the
streamline to take a sudden turn. The fluid flowing along the central
streamline, therefore, stops moving as it reaches the point B. Hence the
velocity of flow at this point is zero. This point is known as stagnation
point.
If the measurement is made on an open channel flow the surface will be exposed to the
air and there is no static head from the surface, and if measurement is made on pipe flow
there will be static head at A.
Fig.5.5 Simple Pitot tube in pipe flow and open channel flow
Applying Bernoulli’s equation to point A in the undisturbed flow region and the
stagnation point B we have
9
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
A perfect Pitot tube should obey this equation exactly, but all actual instruments must be
calibrated and a correction factor applied to make allowance for the small effects of nose
shape and other characteristics.
Practically it is difficult to read h from a free surface. To overcome this difficulty, the
static tube and the Pitot tube may combine in one instrument (differential U-tube).
For finding the velocity at any point in a pipe by Pitot tube, the following arrangements
are adopted.
1. Pitot tube along with a vertical piezometer tube as shown in Fig.5.6
2. Pitot tube connected with piezometer tube as shown in Fig.5.7
3. Pitot static tube, which consists of two circular concentric tubes one inside the
other with some annular space in between as shown in Fig.5.8. The outlets of
these two tubes are connected to the differential manometer where the difference
of pressure head is measured by knowing the difference of the levels of the
manometer liquid hm.
h
Piezometer
Tube Pitot
tube
Fig.5.6
Fig.5.7
10
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
Fig.5.9
The Pitot tube measures the velocity of only a filament of liquid, and hence it can be used
for exploring the velocity distribution across the pipe cross-section. If, however, it is
desired to measure the total flow of fluid through the pipe, the velocity must be measured
at various distances from the walls and the results integrated. The total flow rate can be
calculated from a single reading only if the velocity distribution across the cross-section
is already known.
The static tube measures the static pressure, since there is no velocity component
perpendicular to its opening and the impact tube measures both the static and impact
pressure (due to kinetic energy). Impact tube head = pressure head + velocity head.
3. Orifices
An orifice is an opening (usually circular) in the wall of a tank or in a plate normal to the
axis of a pipe, the plate being either at the end of the pipe or in some intermediate
location and used for measuring rate of flow out of a reservoir (tank) or through a pipe.
The orifice meter consists of a flat orifice plate with a circular hole drilled in it. There is a
pressure tap upstream from the orifice plate and another just downstream.
11
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
Applying the Bernoulli equation between at 1 (upstream of plate) and 2 (at the orifice)
Q = A1V1 =A1
The actual discharge will be less than the theoretical since the effective flow area near P 2
tapping will be less than A2, the fluid forming a neck or vena contracta. In addition,
there will be some loss of energy between point 1 and 2.
The actual discharge can be determined by determining coefficient of discharge.
12
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
Applying Bernoulli’s equation and neglecting losses between A and B (taking the datum
level at the center of the orifice),
0+0+h =
V= (Theoretical velocity). Where h is the depth of center of
orifice below the free surface. Evidently, the velocity of efflux is equal to the velocity of
free fall from the surface of reservoir. This is known as Torricelli’s theorem. In real
fluid flows, there is always some loss of energy due to viscous effects and accordingly
the actual velocity will be less than the theoretical velocity.
1) Coefficient of contraction, Cc =
2) Coefficient of velocity, Cv =
3) Coefficient of discharge, Cd =
Theoretical discharge, Qt = Ao
Where Ao = (Area of jet at C)
Actual velocity, Va = Cv
As shown in the figure, the paths of the particles of the liquid converge on the orifice so
that the area of the issuing jet is less than the area of the orifice. In the plane of the
orifice, the particles have a component of velocity towards the center so that at C the
pressure is greater than atmospheric pressure. It is only at B a little outside the orifice that
the paths of the particles become parallel. The section through B, which is the section of
least cross-section and hence of maximum contraction, is called the vena contracta.
Area of vena contracta, AB = CcAo (The actual area of the jet is the area of
the vena contracta not the area of the orifice.)
13
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
Determination of Cc and Cv
1) Trajectory method
The equation for the trajectory may be obtained by applying Newton’s equation of
uniformly accelerated motion to a particle of the liquid passing from the nozzle to point
P, whose coordinates are (x, y) in time t. Then x = V xot and z = Vzot - 1/2gt2. Evaluating t
from the first equation and substituting in the second gives
If the jet is initially horizontal, as in the flow from a vertical orifice, Vx0 = V0 and Vz0 = 0,
the above equation is reduced in to;
V0 =
Then Cv =
Pitot tube can be set at the vena contracta so that actual velocity Va is determined.
3) Calipers method
The diameter of the jet at the vena contracta can be approximately measured using an
outside caliper. However, this method is not precise and is less satisfactory than other
methods.
14
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
Types of Orifice
The following figure shows common types of orifice with their coefficient of discharge.
The orifices are classified based on of their size, shape, nature of discharge and shape of
the upstream edge.
1. Depending upon their size: - small orifice and large orifice; If the head of liquid
from the center of the orifice is more than five times the depth of the orifice, the
orifice is called small orifice. If it less than five times it is known as large orifice.
2. Depending upon shape: - as circular, triangular, rectangular and trapezoidal
3. Depending upon shape of edge: - as sharp edged and round or bell mouthed
orifice. (Fig.5.14)
4. Depending up on the nature of discharge : - as free discharging & drowned or
submerged orifice; the submerged orifices are further classified as fully submerged
and partially submerged orifice
y
15
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
4. Weirs
Open channel flow may be measured by a weir or obstruction in the channel that causes
the liquid to back up behind it and flow over it or through it. There may be sharp crested
or broad crested based on their length along the channel section.
a) Rectangular weir
The following figure shows rectangular notch of crest length (L) and working under a
head H.
Total discharge
16
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
Since the velocity of flow through the notch varies from top to bottom, consider a strip of
thickness h at a depth h below the surface. If the velocity of approach is small:
17
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
It is often important to determine the force produced on a solid body by fluid flowing
steadily over it. For example, the force on a pipe bend caused by the fluid flowing
through it; the force exerted by jet of fluid striking against a solid surface; thrust on a
propeller. All these forces are hydrodynamic forces and they are associated with a
change in the momentum of the fluid.
F=m*a
F=m
F * dt = m * dv
dFx =
A = area of jet
v1 = velocity of jet
19
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
Continuity equation Q = Q1 + Q2 Q2 = Q – Q1
Substitute in (*) Q 1 = Q - Q2
Q1 – Q + Q1 – Qcos = 0
2Q1 - Q (1+ cos) = 0 Q – Q2 – Q2 – Qcos = 0
Q1 = Q/2 (1 + cos) Q (1 - cos) – 2 Q2 = 0
Q2 = Q/2 (1- cos)
For a vertical plate, = 90o
Q1 = Q/2 = Q2
21
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
The force exerted by the impact of jet can be fruitfully utilized if the flat plates are
mounted on the periphery of a wheel as shown below. The force exerted by the jet causes
the rotation of the wheel. The flat plates thus occupy the bottom most position according
to their turn.
22
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
The number and location of the plates is so arranged that no portion of jet goes waste
without doing work on the plate.
Efficiency of jet =
Consider a symmetrical curved vane having smooth surface. The jet strikes at the center
and after impact it is deflected equally along the vane surface.
Efficiency,
Let the wheel rotate with a tangential velocity U and a jet moving with velocity v strikes
the wheel.
24
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
Let Fx and Fy be the components of the force exerted on the fluid by the pipe bend.
Then momentum equation in the x–direction can be written as:
P1A1 - P2A2cos - Fx = Q(V2cos - V1)
Fx = Q (V1 -V2cos) - P1A1 + P2A2cos
Momentum equation in the y–direction can be written as:
Fy - P2 A2sin = pQ (V2sin - V1) (V1 = 0)
Fy = P2A2sin + pQ(V2sin)
A reducing bend with deviation in the vertical plane is shown in Fig. below. Due to the
hydrostatic and dynamic pressures a force is exerted by the fluid on the bend which has to
be resisted by a thrust block or other suitable means. This force could be evaluated by
plotting the stream lines and thus determining the pressure distribution. However, by a
simple application of the momentum equation, and quite independently of any energy
losses associated with turbulent eddying (real fluid), we obtain.
25
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
3. Force in nozzle
26
Lecture Note 2007 AMU/WTI
Hydraulics I Chapter 5 Fluid
Dynamics
BY SALMAN ABDU
27
Lecture Note 2007 AMU/WTI