Fluid Mechanics Summary Notes PDF
Fluid Mechanics Summary Notes PDF
SUMMARY NOTES
Literature
1. Fluid Mechanics (including Hydraulic Machines) – Dr. A.
K. Jain, Khanna Publishers, Delhi, 2003
2. Fluid Mechanics (6th edition) – frank M. White;
McGraw-Hill 2008
3. Introduction to Engineering Fluid Mechanics.- J. A. Fox
1985
4. Fluid Mechanics:- J. F. Douglas; J. M. Gasiorek; J. A.
Swaffield
5. Hydraulics,Fluid Mechanics and Fluid Machines – S.
Ramamrutham
Literature cont.
6. Essentials of Engineering Hydraulics – J. M. K. Dake,
1992
7. Hydrology and Hydraulic Systems – Ram S. Gupta
Mechanics of Fluids – Bernard Massey, revised by john
Ward- Smith
WHAT IS A FLUID?
Molecules of solids are so closely packed
together that the attractive forces between
the molecules are so large that a solid tends
to retain it’s shape unless compelled by
some external forces to change it. Fluids
are composed of molecules with relatively
larger distances between molecules and
therefore the attractive forces between
molecules are smaller than in solids.
WHAT IS A FLUID?
F
F
Shear τ
Shear τ t1 t3
t0 t2
θ3
θ θ2
θ
Solid Fluid
Definition of fluid
A fluid may be defined in two perspectives:-
fluid dynamics
Fluid Statics
Fluid statics : is the study of the
behaviour of fluids at rest. Since for a fluid
at rest there can be no shearing forces all
forces considered in fluid statics are normal
forces to the planes on which they act.
Fluid Kinematics
Fluid kinematics: deals with the geometry
(streamlines and velocities ) of motion
without consideration of the forces
causin g the motion. Kinematics is
concerned with a description of how
fluid bodies move.
Fluid dynamics
Fluid dynamics: is concerned with the
relations between velocities and
accelerations and the forces causing the
motion.
SYSTEM AND CONTROL
VOLUME
In the study of fluid mechanics, we make
use of the basic laws in physics namely:
The conservation of matter (which is called
A system
Control volume
Control surface
FORCES ACTING ON FLUIDS
(BODY & SURFACE FORCES)
Those forces on a body whose distributions act on
matter without the requirement of direct contact
are called body forces (e.g. gravity, magnetic,
inertia, etc.
Body forces are given on the basis of the force per
unit mass of the material acted on.
Those forces on a body that arise from direct contact
of this body with other surrounding media are called
surface forces eg. pressure force, frictional force,
surface tension
FLUID PROPERTIES
Property :- is a characteristic of a substance
which is invariant when the substance is in a
particular state. In each state the condition
of the substance is unique and is described
by its properties. The properties of a fluid
system uniquely determine the state of the
system.
EXTENSIVE & INTENSIVE
PROPERTIES
Extensive Properties: those properties of the
substance whose measure depends on the
amount of the substance present (weight,
momentum, volume, energy)
Intensive Properties: those properties whose
measure is independent of the amount of
substance present (temperature, pressure,
viscosity, surface tension, mass density etc.
volume per unit mass v and energy per unit mass
e)
PHYSICAL PROPERTIES OF
FLUIDS
Each fluid property is important in a particular
field of application.
Viscosity plays an important role in the
problems of hydraulic friction.
Mass density is important in uniform flow.
Compressibility is a factor in water hammer.
Vapour pressure is a factor in high velocity flow
Mass density & unit
(specific) weight
Mass density and unit weight are the two
important parameters that tend to indicate
heaviness of a substance
Mass density is the mass per unit volume
a d
FLUIDS PROPERTY-
VISCOSITY
At any point at a distance y from the
lower plate, the velocity
U(y) = Uo * (y/Y) Uodt/Y =θ
(du/dy) = (Uo/Y) (θ/dt)=Uo/Y
Experiments show that, other quantities being
held constant F is directly proportional to the
A (area) and the velocity U and inversely
proportional to the distance between the
plates Y
FLUIDS PROPERTY-
VISCOSITY
AU o
F where is the proportionality factor.
Y
The shear stress is defined as Γ (Tau)
F U o du
τ= ; =
A Y dy
4 cos
h
d
Where h – capillary rise
σ – surface tension force per unit length
d – diameter; γ – weight density of fluid and
HYDROSTATICS
Hydrostatic deals with fluid at rest. Hydrostatics
studies the laws governing the behaviour of
fluid at equilibrium when it is subjected to
external and internal forces and bodies at
equilibrium when they are immersed in the fluid.
Shear stress in a fluid at rest is always zero.
Therefore in fluid at rest, the only stress we shall
be dealing with is normal stresses.
WHAT IS HYDROSTATIC
PRESSURE?
HYDROSTATIC PRESSURE
The basic concept of hydrostatics is the
concept of hydrostatic pressure. What is it?
lim it [ p ]
p .
A0 A
Isentropic conditions
dT = Kdz →→ dz = dT/K
Piezometric tube
Pipe
Advantages and Disadvantages
of piezometric tube
Advantages: i) Cheap, easy to install and read
Disadvantages:
i) Requires unusually long tube to measure even
moderate pressures
ii) Cannot measure gas pressures (gases cannot
form free surface)
iii) Cannot measure negative pressures
(atmospheric air will enter the pipe through the
tube).
Manometers:
To overcome the above mentioned limitations of the
piezometer, an improved form of the piezometer
consisting of a bent tube containing one or more
fluids of different specific gravities is used. Such a
tube is called a manometer.
Types of manometers
Simple manometer
Inclined manometer
Micro manometer
Differential manometer
Inverted differential manometer
Manometers:-
Simple manometer
A simple manometer: consists of a tube bent in
U-shape, one end of which is attached to the
gauge point and the other is opened to the
atmosphere. The fluid used in the bent tube
is called the manometric fluid (usually
mercury) and the fluid whose pressure is to
be measure and therefore exerts pressure on
the manometric fluid is referred to as the
working fluid.
Simple manometer
Simple manometer
measuring vacuum pressure
Simple manometer
measuring gauge pressure
Inclined manometer
Working fluid
Manometric fluid
Simple manometer
By using the principle that the pressure on the horizontal and in the same continuous fluid is
the same, we shall state that:
For diagram A
P1=P2
P1=PA + h11
P2=Pa + h22
PA + h11= Pa + h22
PA-Pa= h22- h11
For diagram B
P2=Pa=P1
P1=PB + h11 + h22
PB + h11 + h22=Pa
PB-Pa= -h11 - h22 = vacuum gauge
INCLINED TUBE MANOMETER
This type is more sensitive than the vertical
tube type. Due to the inclination the distance
moved by the manometric fluid in the narrow
tube will be comparatively more and thus give
a higher reading for a given pressure
Micro manometers
Micro manometers
The pressure on level 1 is P1 and pressure on level 2 is
P 2.
PB=P1 + w (h+X-dh)
PD=P2+γw(dh+X) +m.h
But PB=PD------on the same horizontal and in a
continuous fluid.
P1 + w (h+X-dh)= P2+γw(dh+X) +m.h
ΔP=P1-P2 = γw(dh+X) + m.h - w (h+X-dh)= γwdh +
γwX+ m.h - wh - w X + wdh
ΔP=P1-P2 = m.h - wh + 2wdh
Micro manometers
By equation of volumes,
D2dh/4=d2h/(2x4) dh=(d/D)2h/2
ΔP=P1-P2 = mh - wh + w (d/D) 2h
ΔP=P1-P2 = mh - w h[1- (d/D) 2]= w h{SG-
[1-(d/D) 2]
Since d/D is very small, the ratio (d/D) 2 can be
taken as zero
Therefore ΔP=P1-P2 = w h{SG-1}
OTHER TYPES OF MANOMETERS
Differential Manometer :consists of a U-tube
containing the manometric fluid. The two ends of
the tubes are connected to the points, whose
differential pressure is to be measured.
Inverted U-tube Differential Manometer
An inverted U-tube differential manometer is used for
measuring difference of low pressures, where
accuracy is the prime consideration. It consists of
an inverted U-tube containing a light liquid.
MECHANICAL GAUGES
Whenever very high fluid pressures are to
be measured mechanical gauges are
best suited for these purposes. A
mechanical gauge is also used for the
measurement of pressures in boilers or
other pipes, where tube gauges cannot
be conveniently used.
Bourdon’s tube pressure gauge
It can be used to measure both negative (vacuum) and
positive (gauge) pressure. It consists of an elliptical tube
ABC, bent into an arc of a circle. When the gauge tube is
connected to the fluid (whose pressure is to be found) at
C, the fluid under pressure flows into the tube. The
Bourdon tube as a result of the increased pressure tends
to straighten out. With an arrangement of pinion and
sector, the elastic deformation of the Bourdon tube rotates
a pointer, which moves over a calibrated scale to read
directly the pressure of the fluid.
Bourdon’s pressure gauge
Mechanical side with Bourdon
tube
Indicator side with card and
dial
Mechanical Details –
Stationary parts
A: Receiver block. This joins the inlet pipe to the
fixed end of the Bourdon tube (1) and secures the
chassis plate (B). The two holes receive screws
that secure the case.
B: Chassis plate. The face card is attached to this. It
contains bearing holes for the axles.
C: Secondary chassis plate. It supports the outer
ends of the axles.
D: Posts to join and space the two chassis plates
Moving Parts
1. Stationary end of Bourdon tube. This
communicates with the inlet pipe through the
receiver block.
2. Moving end of Bourdon tube. This end is sealed.
3. Pivot and pivot pin.
4. Link joining pivot pin to lever (5) with pins to
allow joint rotation.
5. Lever. This an extension of the sector gear (7).
6. Sector gear axle pin.
Moving Parts
7. Sector gear.
8. Indicator needle axle. This has a spur gear that engages the
sector gear (7) and extends through the face to drive the
indicator needle. Due to the short distance between the
lever arm link boss and the pivot pin and the difference
between the effective radius of the sector gear and that of
the spur gear, any motion of the Bourdon tube is greatly
amplified. A small motion of the tube results in a large
motion of the indicator needle.
9. Hair spring to preload the gear train to eliminate gear lash
and hysteresis.
Diaphragm Pressure Gauge
The principle of work of the diaphragm pressure
gauge is similar to that of the Bourdon tube.
However instead of the tube, this gauges
consists of a corrugated diaphragm. When the
gauge is connected to the fluid whose pressure
is to be measured at C, the pressure in the fluid
causes some deformation of the diaphragm.
With the help of pinion arrangement, the elastic
deformation of the diaphragm rotates the
pointer
Diaphragm Pressure
Gauge
Diaphragm Pressure
Gauge
Dead Weight Pressure Gauge
It is an accurate pressure-measuring instrument and is
generally used for the calibration of other pressure
gauge. A dead weight pressure gauge consists of a
piston and a cylinder of known area and connected to a
fluid by a tube. The pressure on the fluid in the pipe is
calculated by:
p=weight/Area of piston
A pressure gauge to be calibrated is fitted on the other
end of the tube. By changing the weight on the piston
the pressure on the fluid is calculated and marked on
the gauge
Dead Weight Pressure
Gauge
RELATIVE EQUILIBRIUM OF LIQUID
(Liquid under constant acceleration or
constant angular speed)
1
dp a x dx a y dy a z dz ----------------------------( * )
On the accelerating fluid, there are two body forces acting, namely gravity force and inertia
force. From the above equation, we recognise that
ax = -a; ay = 0; az =-g -------------------------------( ** )
-
Po
+
Fig 2-8
Motion in the vertical plane
with constant acceleration
The body forces on such a body are the forces of gravity and inertia. The projections of their
acceleration on the axis are;
1
dp ( g j ) dz 2.41
-j – when ascending
Integrating
p = + (-g + j) Z + C---------------2.42
When Z=0; p=Po
p = (-g + j) Z + Po -------------2.43
p = g(-1 + j/g) Z + Po
p = (-1 + j/g) Z + Po
Motion in the vertical plane
with constant acceleration
Let us represent (-1 + j/g) by k
Then we have
P = -k Z + Po ------------------2.44
Since k is a scalar quantity, we can bring the above expression to the familiar hydrostatic
equation.
Representing -k = 1, we have
p =Po + 1Z ------------------2.45
Though k is a scalar quantity, it can have different values. Let’s look at the different values of k.
1. when j<g, k<1 and becomes small, so the liquid experiences a certain amount of
weightlessness
2. when j = g, k=0 and = 0. Liquid experiences a total weightlessness.
EQUILIBRIUM OF A
ROTATING CONTAINER
EQUILIBRIUM OF A
ROTATING CONTAINER
Consider a cylindrical container filled with a
liquid and rotating with a constant angular
velocity ω about the vertical axis. As a result
of the liquid rotating with the same angular
velocity as the container the liquid is
considered to be at rest relative to the
container. Frictional force (both internal, and
external i.e. friction between particles of
liquid walls) is zero.
EQUILIBRIUM OF A
ROTATING CONTAINER
If the coordinate axis shown on the diagram is
considered fixed to the container, then
relative to the rotating vessel, the liquid will
also be at rest. Therefore the basic
differential equation of hydrostatic of Euler is
applicable in the case of a rotating fluid with
the above conditions.
The body forces acting on the fluid are:
EQUILIBRIUM OF A
ROTATING CONTAINER
1. Gravity dFG = gdM or Z = -g
2. Centripetal force dFCP
v2
dM 2 rdM
r
2
p po x 2
y 2 z Distribution of pressure in the liquid
2
To find lines of constant pressure (isobars) we put the left land side of the equation to zero.
p=constant.
But since Po is atmospheric, we can put p-po=0
Therefore equation of isobars is given by
2
FORCES OF HYDROSTATIC PRESSURE ON
PLANE SURFACES IMMERSED IN FLUIDS
A
FORCES OF HYDROSTATIC PRESSURE ON
PLANE SURFACES IMMERSED IN FLUIDS
The F A p o . A z C A. sin
But .. z C sin hC
F A p o . A hC A
Where paA is the force due to atmospheric
pressure, which is transmitted through out
the fluid onto the planes surface AB and
γhCA is the force due to pressure of the
column of fluid on the surface AB
FORCES OF HYDROSTATIC PRESSURE ON
PLANE SURFACES IMMERSED IN FLUIDS
The moment of the resultant force about the same axis Ox is given by:
Equating equations (2.10) and (2.11) according to the theory of moments, we have
dA
2
z I ox
zD ---------------------
zC .A ( St ) ox
Finding (Centre of
Pressure)
Where Iox = ⌠z2dA – 2nd moment of area or moment of inertia of AB about the
axis ox
and (St)ox = zc.A - 1st moment of area or the static moment of AB about the axis
Ox.
It is also known from the theory of moments that the moment of inertia of a body about a given
axis equals the moment of inertia about an axis parallel to the given axis and passing through the
centre of gravity (centroid) plus the product of the area and the square of the distance between
the axes; i.e.
IC
h D hC Sin 2
hC A
GRAPHICAL METHOD FOR FINDING
HYDROSTATIC FORCE ON PLANE SURFACES
Ox
O x
h
b
m m
H
C C
F D D
B A A
γH
PRESSURE DIAGRAM METHOD
Properties of the pressure diagram
Every ordinate on the pressure diagram gives the
F.H F H
RT . sin .H 1 1 2 . 2
2. 3 2 3
F1 F2 F
RT . sin R B . sin
2 2 2
HYDROSTATIC PRESSURE
FORCES ON CURVED SURFACES
Consider a curved surface ABC with length b.
Let Px and Pz be the horizontal and vertical
components of the force due to hydrostatic
pressure acting on the curved surface. To
find these components lets erect the plane
DE. The plane DE will isolate that volume of
liquid ABCED whose equilibrium we wish to
investigate.
The volume ABCED is acted upon by the ff:
the force Ph acting on the vertical side DE
the force RD-reaction of the base EC
RD=[area (C1CED)] b
the reaction R from the curved surface. Rx, Rz is
the horizontal and vertical components
respectively.
force due to liquid’s own weight G
G=[area (ABCED)] b
Now lets resolve all forces acting on the volume
ABCED onto the x- and z- axis.
Rx=0; Ph – Rx = 0
Ph = Rx = Px
Rz=0; G + Rz - RD = 0
Rz = RD-G
Pz = -Rz = G - RD
Pz= [area ABCED – area C1CED] b
Pz= - [area ABCC1] b
Horizontal component
1. The horizontal component Px of the
force on a curved surface equals the
force of hydrostatic pressure on the
plane vertical figure DE, which is a
projection of the curved surface on the
vertical plane
Vertical Component
2. The vertical component Pz equals the
weight of the imaginary free body of
the fluid ABCC1. This imaginary free
body of the fluid we shall called
"pressure body".
The weight of the pressure body
represent by [area ABCC'] b = Go
Procedure for determining
the horizontal component
1. Place a vertical plane DE behind the
curved surface.
2. Project the curved surface onto the
vertical plane to obtain a plane surface.
3. Determine the horizontal component in a
similar manner as in plane surfaces
immersed in fluid.
Procedure for determining the
vertical component
The cylindrical surface ABC is the surface whose
pressure body is to be found.
1. First fix the extreme ends A and C of the curved
surface;
2. Draw vertical lines from these points to the water
surface;
3. finally note the contour of the pressure body
A'ABCC' ie the body of fluid between the two
vertical lines, the curved surface and the surface of
the fluid.
Procedure for determining the
vertical component cont.-
The cross-section of the pressure body
(positive or negative) is the area between
the two verticals, the cylindrical surface ABC
and the surface of the fluid (or their
continuation). If the pressure body does not
wet the cylindrical surface, then we have
negative body pressure; however if the
pressure body wets the surface, then the
pressure body is positive
Buoyancy
FLOATING BODIES: ARCHIMEDES PRINCIPLE:
the force, which a fluid exerts on a body
immersed in it equals the weight of the fluid
displaced by the body or when a body is
placed (submerged) in a fluid, it experiences
an upward (upthrust) force which is equal to
the weight of fluid the body displaces.
Buoyancy
d
C D
G
Buoyancy
The body AB with volume V completely
submerged in a fluid. The resultant of all
forces due to pressure acting on the surface
element of the body is determined by the
principle of forces on a curved surface.
R R R R
2
x
2
y
2
z
F h2 h1 dA dV V
0 0
-Δh
Vo
V=Vo+aΔh V=Vo-aΔh
The Hydrometer
Let the hydrometer read 1.0 when floating in
distilled water of specific gravity 1. The
corresponding weight of water displaced will be
Voγw; where Vo is the volume of distilled water
displaced. In another liquid of higher (or lower)
density, the hydrometer will pop up (or down)
by an amount Δh. If the stem of the hydrometer
is of cross-sectional area a, the reduction (or
increase) in volume of fluid displaced will be
a.Δh
The Hydrometer
Since the weight of the hydrometer is equal to
the weight of the volume of fluid displaced in
each case
G = γwVo = γf(Vo-a.Δh)
f V o wV o
h
Vo
S .G 1
a f a f
f
KINEMATICS
Kinematics: the study of the geometry of
motion, without considering the forces
causing the motion.
KINEMATICS OF FLUID FLOW
In the 1 8th century, mathematicians sought to
specify fluid motion by mathematical relations.
It must be noted that these relations could be
developed only after certain simplifying
assumptions, notable of which was the concept
of “ideal fluid”, which assumed the fluid as
not having viscosity and not
compressible. The ideal fluid exhibited no
surface tension and could not vaporize if
it was a liquid.
KINEMATICS OF FLUID FLOW
As a result of such assumptions, the relations
obtained for describing the flow of an ideal
fluid may be used to indicate the behaviour
of a real fluid only in certain regions of flow;
e.g. in the regions far removed from
boundaries. The results so obtained may be
only an approximation to the truth, although
in certain cases the theoretical results are
surprisingly close to the actual results.
KINEMATICS OF FLUID FLOW
Irrespective of the way anyone may look at the
relations, they give valuable insight into the
actual behaviour of a real fluid. Therefore in
the forthcoming presentation, we shall only
give an introduction of mathematical
kinematics and its application to a few simple
examples of fluid flow. Attention will be
limited to a steady two-dimensional plane.
TYPES OF FLUID MOTION
(FLOW)
Fluid flow may be classified in a number of ways.
i) Steady and unsteady flows
ii) Uniform and non-uniform flows
filament.
STREAMTUBE AND FILAMENT
Properties of (streamtube)
filament when flow is steady
1). Since streamlines in a steady flow do not change with
time, then filament also does not change its form with
time (i.e. constant form).
2). Since the cross-sectional area of a filament is
elemental, the magnitude of the velocity, U, the
pressure, P, and all fluid properties for all point in a
given cross section of the stream tube are considered
equal. Though U and P are not necessary the same
along the flow.
3). Fluid enclosed in the filament can get out of the tube.
Similarly no particle can enter the stream tube.
The area of a filament normal to streamline
direction is called the filament cross-section dA.
Velocity U, and elementary flow rate dQ are two
fundamental quantities that are used in dealing
with fluid in motion. They give an exact (not
average) differential description of the flow.
Elementary flow rate, dQ
Elementary flow rate
Elementary Flow rate: is the volume of fluid
passing through a given filament cross section
in a unit time (i.e. one second).
The equation of elementary flow rate can be
found by considering fig. 3-5
During the time period dt, all fluid particles from
the section n-n might have moved a distance
dS and might have come to section n'-n’.
ds = Udt
Elementary flow rate
Therefore the volume of fluid passing through
the section n-n during the time dt will be
dV = dA.Udt
Therefore, in a unit time, the volume of fluid
passing through the section n-n will be
dQ = dV/dt = UdA
dQ is the elementary fluid flow rate.
FLOW RATE AND MEAN FLOW
VELOCITY OF A STREAM
A stream consist of numerous filaments. Since flow
velocity of each filament is different from the other
it means the flow velocities at different points in a
Q
given cross section of a stream are different. Since
the velocities are different at different points in a
given cross section, the value of the flow rate of a
stream will be given by the summation of all the
elementary flow over the cross-section
A
Q U .dt
A0
FLOW RATE AND MEAN FLOW
VELOCITY OF A STREAM
Concept of mean flow
velocity
To simplify fluid flow calculations, the concept of mean
flow velocity is introduced. Fluid flowing through a
channel bounded by walls has different velocities at
different points of the cross section. The fluid particles
immediately adjacent to the wall of the tube (channel,
duct etc) adhere to the sides and come to rest. Their
velocities are zero. Filaments in immediate vicinity of the
adhesive particles are dragged because of internal
friction and their velocities are decreased. The farther
the filaments are from the sides of the wall, the greater
their velocity with the maximum at the centre of the
tube
Mean velocity
The mean flow velocity is defined as:
A
Q
U .dA
A0
v
A A
Velocity profile
Let us represent the area of this diagram of A
and lets suppose the stream has a
rectangular cross section with width b
The flow rate is given by:
Q = A.b
CONDITION FOR
CONTINUITY OF FLOW
CONDITION FOR
CONTINUITY OF FLOW
Consider the sections 1-1 and 2-2 of a filament
in a steady flow. We can write that
dQ1 = U1d1
dQ2 = U2d2
It can be seen that
1). dQ1 not greater than dQ2 (because of
incompressibility of fluid)
2). dQ1 not less than dQ2 (because we never
observe a break in the flow)
CONTINUITY EQUATION
Therefore we can write
dQ1 = dQ2 or U1dA1= U2dA2
or dQ = UdA = const
This equation is equally true for a stream i.e.
For any two sections in a stream,
Q1 = Q2 or 1A2 = 2A2
Q=v.A = const. -------- Continuity equation
for a stream
VELOCITY
The velocity of flow for most engineering
problems is of great importance. For flows
past structural or machine parts, knowledge
of the velocity makes it possible to calculate
pressures and forces acting on the structure.
In other cases of engineering as design of
canals and bridge pier, velocity is of interest
from the point of view of its scouring action.
Therefore it is importance to know how to
determine the velocity of flow.
TWO VIEW POINTS ON
FINDING VELOCITY
As particles move in space, their characteristics,
such as velocity, density, etc may change
with space and time. The flow characteristics
are measured with respect to some co-
ordinate system, fixed or moving.
There exist two approaches for finding the
velocity of flow, namely:
i) The Lagrangian Approach
ii) The Eulerian Approach
The Lagrangian Approach
(“Follow that particle”)
When we choose a co-ordinate system attached
to the particles whiles they move. In this
approach, we follow the movement of
individual particles. This means that the
coordinates x, y, z are not fixed but must
vary continuously in such a way as always to
locate the particle. For any particular particle,
x(t), y(t) and z(t) becomes specific time
function which are different for corresponding
time function of other particles
The lagrangian Approach
If the position vector is known, the velocity
could be obtained by differentiating the
position vector with respect to time. For
example if the position vector is expressed in
terms of its components in the x, y, z as:
F(t) = xi + yj + zk
When the equation is differentiated with
respect to time, we obtain the velocity of the
particle as:
The Lagrangian Approach
dx dy dz
u .i .j i
dt dt dt
u u xi u y j u z k
u u 2
x u u
2
y
2
z
The Lagrangian Approach
The difficulty in using this method is that
the motion of one particle is inadequate to
describe an entire flow field. It implies that
the motion of all fluid particles must be
Considered simultaneously which is rather
difficult if not an impossible task.
Eulerian Approach
(“Watch that Space”)
In this approach, choose a co-ordinate system
fixed in space and study the motion of fluid
particles passing through these points.
We fix points in the fluid flow and monitor the
velocity field with time. Hence by this technique,
we express at a fixed positions in space the
velocities of a continuous “string” of fluid
particles moving by this position.In this case the
velocity depends on the point in space and time
Eulerian Approach
ux = f1 (x, y, z, t)
uy = f2 (x, y, z, t)
uz = f3 (x, y, z, t)
u = √(u2x + u2y + u2z)
Since it is almost impossible to keep track of
the position of all the particles in a flow field,
the Eulerian approach is favoured over the
Lagrangian approach.
Velocity as function of
position along a streamline
At times it is useful to express velocity as a
function of position along a streamline and
time as
u = f(s, t)
ACCELERATION
The acceleration of a fluid particle is obtained by
differentiating the velocity with respect to time
uX=f1(x, y, z, t)
uy=f1(x, y, z, t)
uz=f1(x, y, z, t)
dv v dx v dy v dz v dt v v v v
ay . . . . u. v. w.
dt x dt dy dt z dt t dt x y z t
dw w dx w dy w dz w dt w w w w
az . . . . u. v. w.
dt x dt dy dt z dt t dt x y z t
ACCELERATION
a a a a 2
x
2
y
2
z
R
The Continuity Equation
The continuity equation is an expression of the
conservation of mass law and it states that for a
steady flow of fluid in the three-dimensional fluid
element (parallelepiped) of size dx, dy, dz,, the
amount entering the element must be equal to
the amount leaving and for unsteady flow, the
difference between the amount entering and
amount leaving must be stored in the
parallelepiped and this is only possible if density
changes occur in the element.
The Continuity Equation
The Continuity Equation
Let us find the mass of fluid entering the side
ABCD and leaving the side A1B1C1D1 of the
element within a certain interval of time dt
Mass entering the side ABCD
δMe = ρu.dt.dy.dz
And the mass leaving the side A1B1C1D1
δMl = ρ’u’ dt.dy.dz
Note that ρ’ =ρ + (δρ/δx).dx and
u’ = u + (δu/δx).dx
The Continuity Equation
Net mass of fluidbeing retained in the element
= Mass entering and mass leaving within the
time dt is given by:
δMe - δMl = ρu.dt.dy.dz - ρ’u’ dt.dy.dz =
- δ(ρu)/δx).dx.dy.dz.dt
- ( u )
M x dx.dy .dz .dt
x
The Continuity Equation
( u )
Similarly ...... M y dx .dy .dz .dt
y
( u )
.......... .......... M z dx .ddy .dz .dt
z
Therefore the total net gain of mass within the
time dt within the element is given by:
( u ) ( v ) ( w)
M dx.dy .dz .dt
x y z
The Continuity Equation
This gain in mass within the element is only
possible if within the period dt there were
changes in density within the element.
If the density of the fluid within the element
at the time t = 0 was ρ and the density at
the end of the period ie time dy was ρ’ then
the mass of the fluid at the beginning of the
period was
The Continuity Equation
∂Mt=o = ρ.dx.dy.dz
and the mass at the end of the period, dt was
∂Mt=dt = ρ’.dx.dy.dz
But ρ’ = ρ + (δρ/δt).dt
Therefore net change in mass within the element
in time dt due to density changes in the
element is given by
M M t dt M t 0 dx.dy .dz .dt
t
The Continuity Equation
The change in mass due to difference in
volume entering and leaving must be equal
to the change in mass due to density
changes within the element with the same
time period dt. Therefore
( u ) ( v ) ( w)
dx .dy .dz .dt .dx .dy .dz .dt
x y z t
( u ) ( v ) ( w )
0
t x y z
The Continuity Equation
( u ) ( v ) ( w )
0
t x y z
u z
tan .dx.dt
x
The rate of change of α or the rate of rotation
of side oa u z
dt x
Rotation of a fluid
Similarly for the side oc to move to o’c’, the
point o moved through a horizontal distance
of dx1 and point c moved through dx2.
u x
dx1 = uxdt and dx2 = u’xdt = u x dz .dt
z
u x u x
dx2 dx1 dz.dt tan .dz.dt
z z
Rotation of a fluid
Rate of change of β or the rate of rotation of the side
oc. u x
dt z
We adopt a sign convention for rotation. Clockwise
rotation is negative and anticlockwise rotation is
positive.
Rotation of the fluid element about its instantaneous
axis ( in this case about the y-axis) is characterised
by its angular velocity which is defined as the
average rate of deformation of side oa and oc which
Rotation of a fluid
1 1 u z u x
2 dt dt 2 x z
1 u z u x
y
2 x z
Similarly 1 u y u z
x
2 z y
1 u x u y
z
2 y x
Rotational and irrotational
flows
When all the components of rotation, i.e Ωx
Ωy, Ωz are equal to zero, it means rotation is
absent and the fluid flow is referred to as
irrotational flow. On the other hand if
even one component of the angular velocity
is not zero, it means rotation exist and the
flow is called rotational flow
VORTICITY
Vorticity is a concept used in fluid mechanics
to define rotation and it is defined as two
times the angular velocity.
u y u z
x
z y
u z u x
y
x z
u x u y
z
y x
VORTICITY
Just like the angular velocity, the vorticity is
defined in the three axes and characterises
rotation of a fluid element about its
instantaneous axis. If vorticity is zero, then
flow is irrotational
CIRCULATION AND
VORTICITY
Consider a closed curve in a two-dimensional
flow field shown in the diagram below.
Streamlines cut the curve. If P is a point of
intersection of the curve with a streamline,
and θ is the angle which the streamline
makes with the curve, then the component of
the velocity along the closed curve at the
point is equal to v.cos θ. The circulation Γ
(gamma) is defined as the line integral of
velocity around a closed curve in a flow
CIRCULATION
CIRCULATION
Thus the differential circulation dΓ along a
small length ds is given by:
dΓ = (vcosθ)ds.
u.......... .......... .... v.......... .......... ..... w
x y z
VELOCITY POTENTIAL
The total differential of the function φ in a two-
dimensional flow can be written as:
d dx dy udx vdy
x y
Since φ is constant along an equipotential line, we can
write;
d dx dy udx vdy 0
x y
d dx dy
r sin r cos
d x d y d x
v x r sin v y r cos r v x sin v y cos rv
d
rv
d
1 d
v
r d
Radial velocity in terms of the
stream function
d dx dy
r sin r cos
d x d y x y
v y .r sin v x .r cos r v x cos v y sin rv r
d
rv r
d
1 d
vr
r d
Transverse velocity in
terms of stream function
d dx dy
. cos sin
dr x dx y dr x y
v y . cos v x . sin v y cos v x . sin v
d
v
dr
d
v
dr
SOURCE & SINK
A source is a point in space from which fluid
issues uniformly in all directions at a constant
rate. For two-dimensional flow, the flow
pattern is made up of streamlines uniformly
spaced and directed radially outwards from one
point in the reference plane. Continuity
principles shows that the velocity will diminish
as the streamlines spread and the symmetry
will require that all velocities will be the same
at the same radial distance from the origin.
SOURCE
Across all circles, the same discharge will pass.
Therefore the velocity at any point r in the flow field
may be determined by the continuity equation as:
q=Vr x A = Vr x 2πrx1 considering unit depth.
Vr = q/2πr
Where q is the constant rate of flow per unit depth
issuing out of the source.
q – usually called the “strength” of the source per unit
depth of the source.
SOURCE
SOURCE
This implies that any circle drawn to enclose the
source will be discharging the same flow q. The
velocity at any point r in the flow field may be
determined as:
Vr = q/2πr; Vθ = 0
It is customarily for the stream function ψ =0 to be
made coincident with the x-axis. Since it is
convenient to express the source in the polar
coordinate, the stream function is written as a
function or r and θ i. e ψ(r, θ).
STREAM FUNCTION FOR
SOURCE
.dr d v .dr r .v r d
r
Integratin g
q
( r , ) v .dr rv r d 0 r d C
2r
q
(r , ) C
2
Choo sin g ...boundary ...conditions ..such ..that ..when .. 0 o ,... 0,..then ..C 0
q
(r , ) .
2
VELOCITY POTENTIAL FOR
SOURCE
d ( r , ) .dr .d
r
Φ(r,θ) = -(q/2π).lnr +c
Φ(r,θ) = -(q/2π).lnr
SINK
Sink: is the exact opposite of a source i.e. a
point in space to which fluid converges
uniformly and from which fluid is
continuously removed. As a result, a sink is
treated as a negative source flow and the
expression for velocities and the functions ψ
and φ are the same as those for a source but
with the signs reversed.
SINK
Consequently the stream function for a
sink is given by:
ψ(r,θ) = - (q/2π).θ
and the velocity potential of a sink given
as:
Φ(r,θ) = (q/2π).lnr
Combination of a source and
a uniform rectilinear flow
A uniform flow in the x- direction with a stream
function ψu = U.y and a source of uniform
strength of stream function ψs = (q/2π).θ
located at the origin O. When these flows are
brought together, the resulting stream function
is obtained by adding the respective stream
functions. i.e
ψcomb = U.y + (q/2π).θ
and since y= r. sin θ
Ψcomb =U.r.sinθ + (q/2π).θ
Graphical superposition of the
two streamlines
Graphical representation is obtained using
Rankine’s method by superposition of the
two streamlines. For example for the ψcomb =
8 is obtained by adding at the points of
intersections ψ1=8 + ψ0 =0; ψ1=7 + ψ2 =1;
ψ1=6 + ψ2 =2; ψ1=5 + ψ2 =3; ψ1=4 + ψ2
=4; ψ1=3 + ψ2 =5; ψ1=2 + ψ2 =7; ψ1=1 +
ψ2 =7; ψ1=0 + ψ2 =3; etc
Graphical superposition of the
two streamlines
All streamlines of the combined flow are obtained
in such manner. It must be observed that the
resulting streamlines are grouped into two
distinct sets. In one set, all streamlines emerge
from the origin (under ψ=8.0 streamline) and in
the other they approach the rectilinear
asymptotically at some distance upstream of
ψ=9. The two sets are separated by the
streamline ψ=8, which passes through the point
S.
Graphical superposition of
the two streamlines
This point is a stagnation point, where the
velocity from the source equals to the
uniform rectilinear velocity, so that the
resultant is zero at S. The distance OS =a
can be obtained by equating the uniform
velocity to that of the source at a radius a,
from the origin. Thus
U= q/2πa
Or a = q/2πU
Graphical superposition of the
two streamlines
Assignment: Plot graphically the combination
of the ff. streamlines. Uniform flow + source
Group 1 & 8: u=1m/s; q=8m3/s/m
Group 2 & 9: u=2 q=16
Group 3 & 10: u=1.5 q=6
Group 4 & 11: u= 2.5 q=12
Group 5 & 12: u=1.25 q=10
Group 6 & 13: u=1.75 q= 7.5
Group 7 & 14: u=2.75 q= 14