Lecture 1427495313
Lecture 1427495313
Lecture 1427495313
TABLE OF CONTENT 1
Definition of a fluid:- 2
Fluid as a continuum :- 2
ρ= 3
Velocity field: 4
One, two and three dimensional flows : 4
Newtonian Fluids : 8
Non-Newtonian Fluids : 9
Vapor Pressure: 11
Surface Tension: 11
Basic flow Analysis Techniques: 12
Flow patterns: 12
CHAPTER – 2 FLUID STATICS 17
Buoyancy: 36
CHAPTER – 3 42
INCOMPRESSIBLE INVISCID FLOW 55
CHAPTER–4 64
Laminar flow through a pipe. 64
Assumptions: 64
CHAPTER-5 76
FUNDAMENTALS OF IRROTATIONAL MOTION 76
CHAPTER 6 79
MATHEMATICAL FORMULATION OF IRROTATIONAL MOTION 79
CONCLUSION 82
REFERENCES: 83
CHAPTER -1
Definition of a
fluid:-
Fluid mechanics deals with the behaviour of fluids at rest and in motion. It is logical to begin
with a definition of fluid. Fluid is a substance that deforms continuously under the application
of shear (tangential) stress no matter how small the stress may be. Alternatively, we may
define a fluid as a substance that cannot sustain a shear stress when at rest.
A solid deforms when a shear stress is applied , but its deformation doesn’t continue to
increase with time.
Fig 1.1(a) shows and 1.1(b) shows the deformation the deformation of solid and fluid under
the action of constant shear force. The deformation in case of solid doesn’t increase with
time i.e t1 t 2 ....... tn .
From solid mechanics we know that the deformation is directly proportional to applied
shear stress ( τ = F/A ),provided the elastic limit of the material is not exceeded.
To repeat the experiment with a fluid between the plates , lets us use a dye marker to outline
a fluid element. When the shear force ‘F’ , is applied to the upper plate , the deformation of
the fluid element continues to increase as long as the force is applied , i.e t 2 t1 .
Fluid as a continuum :-
Fluids are composed of molecules. However, in most engineering applications we are
interested in average or macroscopic effect of many molecules. It is the macroscopic effect
that we ordinarily perceive and measure. We thus treat a fluid as infinitely divisible substance
, i.e continuum and do not concern ourselves with the behaviour of individual molecules.
The concept of continuum is the basis of classical fluid mechanics .The continuum
assumption is valid under normal conditions .However it breaks down whenever the mean
free path of the molecules becomes the same order of magnitude as the smallest significant
characteristic dimension of the problem .In the problems such as rarefied gas flow (as
encountered in flights into the upper reaches of the atmosphere ) , we must abandon the
concept of a continuum in favour of microscopic and statistical point of view.
As a consequence of the continuum assumption, each fluid property is assumed to have a
definite value at every point in the space .Thus fluid properties such as density , temperature ,
velocity and so on are considered to be continuous functions of position and time .
Consider a region of fluid as shown in fig 1.5. We are interested in determining the density at
the point ‘c’, whose coordinates are , and . Thus the mean density V would be given
by ρ= . In general, this will not be the value of the density at point ‘c’ . To determine the
𝑉
density at point ‘c’, we must select a small volume , 𝑉, surrounding point ‘c’ and determine
the ratio ð𝑚 and allowing the volume to shrink continuously in size.
ð𝑉
Assuming that volume 𝑉 is initially relatively larger (but still small compared with volume ,
V) a typical plot might appear as shown in fig 1.5 (b) . When 𝑉 becomes so small that it
contains only a small number of molecules , it becomes impossible to fix a definite value for
; the value will vary erratically as molecules cross into and out of the volume. Thus there
𝑉
is a lower limiting value of 𝑉, designated 𝑉ꞌ . The density at a point is thus defined as
ρ=
𝑚
𝑉𝑉ꞌ . 𝑉
Since point ‘c’ was arbitrary , the density at any other point in the fluid could be determined
in a like manner. If density determinations were made simultaneously at an infinite number of
points in the fluid , we would obtain an expression for the density distribution as function of
the space co-ordinates , ρ = ρ(x,y,z,) , at the given instant.
Clearly , the density at a point may vary with time as a result of work done on or by the fluid
and /or heat transfer to or from the fluid. Thus , the complete representation(the field
representation) is given by :ρ = ρ(x,y,z,t)
Velocity field:
In a manner similar to the density , the velocity field ; assuming fluid to be a continuum , can
property.
Thus for steady flow is 6𝜌 = 0 or ρ = ρ(x,y,z)
6𝑡
6𝑉̅`
=0 or 𝑉̅→ = ̅→(x,y,z)
6𝑡
Thus in steady flow ,any property may vary from point to point in the field , but all properties
, but all properties remain constant with time at every point.
stress 𝑟𝑛
are
:𝜎 then defined as = lim ð𝐹𝑛
𝑛 𝐴𝑛0 ð𝐴
𝑛
ð𝐹𝑡
𝑟 = lim ;Subscript ‘n’ on the stress is included as a reminder that the stresses are
𝑛 𝐴𝑛0 ð𝐴
𝑛
associated with the surface 𝐴̅ , through ‘c’ , having an outward normal in 𝑛^
direction .For any other surface through ‘c’ the values of stresses will be different .
Consider a rectangular co-ordinate system , where stresses act on planes whose normal are in
x,y and z directions.
Fig 1.9
𝜎 = lim
ð𝐹𝑦
𝑥𝑦 𝐴𝑥0 ð𝐴
𝑥
𝜎 = lim ð𝐹z
𝑥𝑧 𝐴𝑥0 ð𝐴
𝑥
A double subscript notation is used to label the stresses. The first subscript indicates the plane
on which the stress acts and the second subscript represents the direction in which the stress
acts, i.e 𝜎𝑥𝑦 represents a stress that acts on x- plane (i.e the normal to the plane is in x
direction ) and acts in ‘y’ direction .
Consideration of area element 𝐴𝑦 would lead to the definition of the stresses , 𝜎𝑦𝑥 , 𝜎𝑦𝑦 and
𝜎𝑦𝑧 . Use of an area element 𝐴𝑧 would similarly lead to the definition 𝜎𝑧 , 𝜎𝑧𝑦 and 𝜎𝑧𝑧 .
An infinite number of planes can be passed through point ‘c’ , resulting in an infinite number
of stresses associated with planes through that point. Fortunately , the state of stress at a point
can be completely described by specifying the stresses acting on three mutually
perpendicular planes through the point.
Thus , the stress at a point is specified by nine components and given by :
𝜎𝑥𝑥 𝜎𝑥𝑦 𝜎𝑥𝑧
𝜎" = 𝜎𝑦𝑦 𝜎𝑦𝑧]
[𝜎𝑦𝑥
𝜎𝑧𝑥 𝜎𝑧𝑦 𝜎𝑧𝑧
Fig 1.10
Viscosity:
In the absence of a shear stress , there will be no deformation. Fluids may be broadly
classified according to the relation between applied shear stress and rate of deformation.
Consider the behaviour of a fluid element between the two infinite plates shown in fig 1.11 .
The upper plate moves at constant velocity , u , under the influence of a constant applied
force ,𝐹𝑥 .
Where , 𝐴𝑦 is the area of contact of a fluid element with the plate. During the interval t ,
the fluid element is deformed from position MNOP to the position 𝑀'𝑁𝑂𝑃' . The rate of
deformation of the fluid element is given by:
Deformation rate = lim
𝑑
𝑡0 𝑡
=
𝑑𝑡
Thus the fluid element when subjected to shear stress , 𝜎𝑦𝑥 , experiences a deformation rate ,
𝑑𝑢
given by .
𝑑𝑦
#Fluids in which shear stress is directly proportional to the rate of deformation are
“Newtonian fluids “ .
# The term Non –Newtonian is used to classify in which shear stress is not directly
proportional to the rate of deformation .
Newtonian Fluids :
Most common fluids i.e Air , water and gasoline are Newtonian fluids under normal
conditions. Mathematically for Newtonian fluid we can write :
𝑑𝑢
𝜎𝑦𝑥 𝖺
𝑑𝑦
If one considers the deformation of two different Newtonian fluids , say Glycerin and water
,one recognizes that they will deform at different rates under the action of same applied
stress. Glycerin exhibits much more resistance to deformation than water . Thus we say it is
more viscous. The constant of proportionality is called , ‘μ’ .
Thus , 𝑑𝑢
𝜎𝑦𝑥 =μ
𝑑𝑦
Non-Newtonian Fluids :
𝑑𝑢
𝜎 =k ( )𝑛
𝑦𝑥 𝑑𝑦 , ‘n’ is flow behaviour index and ‘k’ is consistency index .
To ensure that 𝑑𝑢
has the same sign as that ( , we can express
𝜎𝑦𝑥 )
of 𝑑𝑦
𝑑𝑢
𝜎 =k 𝑛−1 𝑑𝑢 𝑑𝑢
𝑦𝑥 │(𝑑𝑦)│ ( ) = η( )
𝑑𝑦 𝑑𝑦
𝑑 𝑛−1
Where ‘η’ = 𝑢 )│ is referred as apparent viscosity.
│(
k 𝑑𝑦
# The fluids in which the apparent viscosity decreases with increasing deformation rate
(n<1) are called pseudoplastic (shear thining) fluids . Most Non –Newtonian fluids fall into
this category . Examples include : polymer solutions , colloidal suspensions and paper pulp
in water.
# If the apparent viscosity increases with increasing deformation rate (n>1) the fluid is
termed as dilatant( shear thickening). Suspension of starch and sand are examples of dilatant
fluids . # A fluid that behaves as a solid until a minimum yield stress is exceeded and
subsequently exhibits a linear relation between stress and deformation rate .
𝜎𝑦 𝑥
= 𝜎𝑦i𝑒𝑙𝑑
�
�
�
� + μ( )
𝑑𝑦
#Due to strong cohesive forces between the molecules, any layer in a moving fluid tries to
drag the adjacent layer to move with an equal speed and thus produces the effect of viscosity.
#The individual molecules of a fluid are continuously in motion and this motion makes a
possible process of momentum exchange between layers. Such migration of molecules causes
forces of acceleration or deceleration to drag the layers and produces the effect of viscosity.
Although the process of molecular momentum exchange occurs in liquids, the intermolecular
cohesion is the dominant cause of viscosity in a liquid. Since cohesion decreases with
increase in temperature, the liquid viscosity decreases with increase in temperature.
In gases the intermolecular cohesive forces are very small and the viscosity is dictated by
molecular momentum exchange. As the random molecular motion increases wit a rise in
temperature, the viscosity also increases accordingly.
Example-1An infinite plate is moved over a second plate on a layer of liquid. For small
gap width ,d, a linear velocity distribution is assumed in the liquid . Determine :
(i) The shear stress on the upper and lower plate .
(ii) The directions of each shear stresses calculated in (i).
Soln: 𝑑𝑢
𝑦𝑥 =μ𝑑𝑦
𝑟𝑦𝑥 −0 0
=μ ( 0 =μ
) 𝑑
𝑑−𝑜
Hence;
𝑟 𝑦𝑥 𝑦=𝑑 = 𝑦𝑥 𝑦=0
= μ 0 = constant
𝑑
𝑟
Example-2
An oil film of viscosity μ & thickness h<<R lies between a solid wall and a circular disc as
shown in fig E .1.2. The disc is rotated steadily at an angular velocity Ω. Noting that both the
velocity and shear stress vary with radius ‘r’ , derive an expression for the torque ‘T’ required
to rotate the disk.
Soln:
Surface Tension:
Surface tension is the apparent interfacial tensile stress (force per unit length of interface) that
acts whenever a liquid has a density interface, such as when the liquid contacts a gas, vapour,
second liquid, or a solid. The liquid surface appears to act as stretched elastic membrane as
seen by nearly spherical shapes of small droplets and soap bubbles. With some care it may be
possible to place a needle on the water surface and find it supported by surface tension.
A force balance on a segment of interface shows that there is a pressure jump across
the imagined elastic membrane whenever the interface is curved. For a water droplet in air,
the pressure in the water is higher than ambient; the same is true for a gas bubble in liquid.
Surface tension also leads to the phenomenon of capillary waves on a liquid surface and
capillary rise or depression as shown in the figure below.
Basic flow Analysis Techniques:
There are three basic ways to attack a fluid flow problem. They are equally important for
a student learning the subject.
(1) Control–volume or integral analysis
Flow patterns:
Fluid mechanics is a highly visual subject. The pattern of flow can be visualized in a dozen of
different ways . Four basic type of patterns are :
1. Stream line- A streamline is a line drawn in the flow field so that it is tangent to the
line velocity field at a given instant.
2. Path line- Actual path traversed by a fluid particle.
3. Streak line- Streak line is the locus of the particles that have earlier passed through
a prescribed point.
4. Time line – Time line is a set of fluid particles that form a line at a given instant .
Ex: A velocity field given by 𝑉̅→ = A x ı^ – A y 𝑗^ . x, y are in meters . units of velocity in m/s.
A = 0.3 𝑠−1
(a) obtain an equation for stream line in the x,y plane.
(b) Stream line plot through (2,8,0)
(c) Velocity of a particle at a point (2,8,0)
(d) Position at t = 6s of particle located at (2,8,0)
(e) Velocity of particle at position found in (d)
(f) Equation of path line of particle located at (2,8,0) at
t=0 Soln:
𝑑𝑥 𝑑𝑦
(a) For stream lines ; =
𝑢 𝑣
𝑑𝑥 𝑑𝑦
= −𝐴𝑦
𝐴𝑥𝑑𝑥 𝑑𝑦
∫ =-∫
𝑥 𝑦
ln 𝑥 = - ln 𝑦 + C
ln 𝑥𝑦 = C
xy =C
𝑥0𝑦0 = C
𝐶 = 16
xy = 16
𝑥0 𝑥0
𝑑𝑦 𝑦 𝑑𝑦 𝑡
v = - Ay , = -Ay , ∫ =- A∫ 𝑑𝑡
𝑑𝑡 𝑦0 0
𝑦
ln( ) = -At ,
𝑦 = 𝑒−𝐴𝑡
𝑦0 𝑦0
(a) The equation of stream line through (𝑥0, 𝑦0) and equation of the path line traced out by
particle passing through (𝑥0𝑦0)are same as the flow is steady.
(b) In following a particle ( Lagrangian method of description ) , both the coordinates of
the particle (x,y) and the component ( 𝑢𝑝 & 𝑣𝑝) are functions of time.
Example -2:
A flow is described by velocity field, 𝑉̅ =ay ı^ + bt 𝑗^ , where a = 1 𝑠 −1 , b= 0.5 m/𝑠 2 . At
t=2s ,what are the coordinates of the particle that passed through (1,2) at t=0 ? At t=3s , what
are the coordinates of the particle that passed through the point (1,2) at t= 2s .
Plot the path line and streak line through point (1,2) and compare with the stream lines
through the same point ( 1,2) at instant , t = 0,1,2 & 3 s .
Soln:
Path line and streak line are based on parametric equations for a particle .
𝑑𝑦
v = 𝑑𝑡 = bt , so, dy = bt dt
y - 𝑦 = 𝑏 (𝑡2− 𝑡 2)
0 2 0
& u = 𝑑𝑥 = ay = a [
+ 𝑏 (𝑡2− 𝑡 2) ]
𝑑𝑡 0 2 0
𝑥 𝑡
∫ 𝑑𝑥 = ∫ {𝑎[ 𝑦
+ 𝑏 (𝑡2− 𝑡 2)]}dt
𝑥0 𝑡0 0 2 0
3
(𝑥 − 𝑥 ) = a 𝑦 (t- 𝑡 ) + 𝑏 ( 𝑡 − 𝑡 2t )
0 0 0 2 0 𝑡0
3
𝑥 = 𝑥 + a 𝑦 (t - 𝑡 ) + 𝑎𝑏 { (𝑡 3−𝑡0 )3 - 𝑡 2(t- 𝑡 ) }
0 0 0 2 3 0 0
(a) For 𝑡0 = 0 and (𝑥0 , 𝑦0) = (1,2) , at t = 2s , we have
y-2 = 𝑏(4)
2
y =3 m
x = 1 + 2 (2-0) + 0.5 [8 – 0] = 5.67 m
2 3
(b)For 𝑡0 = 2s and (𝑥0 , 𝑦0) = ( 1,2) . Thus at t = 3s
We have , y -2 = 𝑏(𝑡2− 𝑡 2) = 0.5 (9-4) = 1.25
2 0 2
y = 3.25 m
3
−23
& x = 1+ 2 (3-2) + 0.5 { (3 ) - 22(3-2) }
2 3
(c) The streak line at any given ‘t’ may be obtained by varying ‘𝑡0’ .
𝑡0(s) t(s) X Y
2 2 1 2
2 3 3.58 3.25
2 4 7.67 5.0
Thus , c = 𝑦 2 – (2𝑏𝑡) 𝑥
0 𝑎 0
For t =0 ; c = (2)2 = 4
t = 1 ;c = 4 – ( 1 )1 = 3
1
t = 2 ;c = 4 – ( 2 )1 = 2
1
t =3 ;c = 4 – ( 3 )1 = 1
1
t(s) 0 1 2 3
C= 4 3 2 1
X Y Y Y Y
1 2 2 2 2
2 2 2.24 2.45 2.65
3 2 2.45 2.83 3.16
4 2 2.65 3.16 3.61
5 2 2.53 3.46 4.0
6 2 3.0 3.74 4.36
7 2 3.16 4.00 4.69
# Streak line of particles that passed through point (𝑥0 , 𝑦0) at t = 3s.
Mere simplicity doesn’t justify our study of subject . Normal forces transmitted by fluids are
important in many practical situations. Using the principles of hydrostatics we can compute
forces on submerged objects, developed instruments for measuring pressure, forces
developed by hydraulic systems in applications such as industrial press or automobile brakes.
In a static fluid or in a fluid undergoing rigid-body motion, a fluid particle retains its identity
for all time and fluid elements do not deform. Thus we shall apply Newton’s second law of
motion to evaluate the forces acting on the particle.
+ ( p + 6𝑝 𝑑𝑦
)dx dz (-𝑗^) + ( p + 6𝑝 𝑑𝑧
)dx dy (𝑘^ ) +( p + 6𝑝 𝑑𝑧
)dx dy (-𝑘^ )
6 2 6𝑧 2 6𝑧 2
d̅𝐹𝑠→ = - (∇p) dx dy
dz Thus
Or = 𝑎̅ ρ (2.2)
d𝐹→
d
- ∇p + ρ𝑔̅ = 𝑎̅ ρ
-6𝑝 + ρ𝑔 = 0
6𝑧 𝑧
6𝑝 6𝑝
6𝑝 =0, =0 & = −𝜌 ; since P=P(Z)
6𝑥 6𝑦 6𝑧
𝑑𝑝
We can write = −𝜌𝑔
𝑑𝑧
𝑃g𝑎𝑢g𝑒
𝑃𝑎𝑡𝑚
𝑃𝑎𝑏𝑠
P=0
Ex:2.1 A tube of small diameter is dipped into a liquid in an open container. Obtain an
expression for the change in the liquid level within the tube caused by the surface tension.
Soln:
= 𝑀 𝐷2 h
4
Thus ; σ𝑀Dcos - ρg 𝑀 𝐷2 h = 0
4
4𝜎𝑐𝑜𝑠
h = 𝜌g𝐷
To eliminate ‘H’ , we recognise that the volume of manometer liquid remains constant i.e the
volume displaced from the reservoir must be equal to the volume rise in the tube.
𝑙 (𝐷 ) 𝑙 (𝐷 )
𝑃
Thus, L= 𝑑 2
𝜌𝑙g[sin + ) ]
(
𝐷
𝐿
Thus , S =
ℎ𝑒 1
= 𝑆[𝑠i𝑛 𝑑)2
+( ]
𝐷
𝜌𝑒
Where , SG =
𝜌w
The expression ‘S’ for sensitivity shows that to increase sensitivity SG , sin and 𝑑 should be
𝐷
made as small as possible.
Hydrostatic Force on the plane surface which is inclined at an angle ‘’ to
horizontal free surface:
We wish to determine the resultant hydrostatic force on the plane surface which is inclined at
angle ‘’ to the horizontal free surface.
Since there can be no shear stresses in a static fluid , the hydrostatic force on any element of
the surface must act normal to the surface .The pressure force acting on an element d𝐴̅ of the
upper surface is given by d𝐹̅ = - p d𝐴̅ .
The negative sign indicates that the pressure force acts against the surface i.e in the
.
direction opposite to the area d𝐴̅ .̅𝐹̅̅ = ∫ −𝑝𝑑𝐴̅
𝑅 𝐴
If the free
.
surface is at a pressure ( 0 = 𝑃𝑎𝑡𝑚), then , p = 𝑝0 + ρgh
|̅𝐹̅̅| = ∫ (𝑝 + pgh)dA = 𝑝 𝐴+ ∫. 𝜌𝑔𝑦 𝑠i𝑛 𝑑𝐴
𝑅 𝐴 0 0 𝐴
|̅𝐹̅̅| = 𝑝 𝐴 + ρg sin ∫ 𝑦𝑑𝐴
.
. 𝑅 0 𝐴
But ∫ 𝑦𝑑𝐴 = 𝑦𝑐 dA
𝐴
𝑦* 𝐹𝑅
= ∫. 𝑦𝑝𝑑𝐴 and 𝑥* 𝐹𝑅 = ∫𝐴 𝑥𝑝𝑑𝐴 #when the ambient (atmospheric) pressure , 𝑝0 ,
.
𝐴
acts on both sides of the surface , then 𝑝0 makes no contribution to the net hydrostatic force
on the surface and it may be dropped . If the free surface is at a different pressure from the
ambient, then ‘ 0′ should be stated as
gauge pressure , while calculating the
net force .
. .
∫ 𝑝𝑦𝑑𝐴 ∫ 𝜌g𝑥2𝑠i𝑛 𝑑𝐴
𝑦* = 𝐴
= 𝐴
𝐹𝑅 𝜌g𝑦𝑐𝐴 𝑠i𝑛
𝜌g𝑠i𝑛 ∫ 𝑦2 𝑑𝐴
𝑦* = 𝜌g𝑦𝑐𝐴 𝑠i𝑛
𝑦* = 𝐼𝑥𝑥
𝐴𝑦𝑐
Where 𝑥^𝑥^ is the second moment of the area about the centroid al ‘𝑥^’ axis . Thus
𝐼^𝑥^𝑥
𝑦 * = 𝑦𝑐 + 𝐴𝑦𝑐
Or , 𝑦* = ( ℎ𝑐 ) + 𝐼𝑥^𝑥^ 𝑠i𝑛
𝑠i𝑛 𝐴ℎ𝑐
𝑥* 𝐹𝑅 = ∫ 𝑥𝑝𝑑𝐴
. .
𝑥* ρg sin 𝑦 A = ∫ 𝑥𝜌𝑔ℎ 𝑑𝐴 = ρgsin∫ 𝑥𝑦𝑑𝐴
𝑐 𝐴 𝐴
.
∫ 𝑥𝑦𝑑𝐴
*
𝑥 = 𝐴
= 𝐼𝑥𝑦
𝐴𝑦𝑐 𝐴𝑦𝑐
Where 𝐼𝑥^𝑦^ is the area product of inertia w.r.t centroid al 𝑥^𝑦^ axis.
So, 𝑥* = 𝐼𝑥^𝑦^
𝑥𝑐 + 𝐴𝑦𝑐
For surface that is symmetric about ‘y’ axis , 𝑥* = 𝑥𝑐 and hence usually not asked to evaluate.
Example Problem:
Ex 2.4:Rectangular gate , hinged at ‘A’ , w=5m . Find the resultant force , ̅𝐹̅𝑅̅ , of the
water and the air on the gate .The inclined surface shown , hinged along edge ‘A’ , is 5m
wide .
Determine the resultant force ̅𝐹̅𝑅̅ , of the water and air on the inclined surface.
,
Soln:-
̅𝐹̅ = ∫. 𝑝 𝑑𝐴̅ = − ∫8 𝜌g y sin30 w dy 𝑘^
4
𝑅 𝐴 2
𝐹̅̅̅ = -𝜌gw 𝑘^ [ ]8 = - 999×9.81×5 [64-16]𝑘^
𝑅 2 2 4 4
̅𝐹̅𝑅̅ = -588.01 KN
Force acts in negative ‘z’ direction.
Taking moment about x axis through point ‘ O ’ on the free surface , we obtain :
. 8
𝑦* 𝐹 = ∫ 𝑦 𝑝𝑑𝐴 = ∫ 𝑦 𝜌 𝑔 𝑠i𝑛30 𝑤 𝑑𝑦
𝑅 𝐴 4
𝑦* 𝐹 𝜌gw 𝑦3 8 5×999×9.81
𝑅 =( 2
) [ 3] 4 = 6
[ 83 - 43]
𝑦* ×(588.01 ×103) = 3658.73×103
𝑦* = 6.22 m
#To find 𝑥* ; we can take moment about y axis through point ‘o’.
𝑥* 𝐹𝑅 = . 𝑥 𝑝 𝑑𝐴 = w 8
∫𝐴 ∫0 ∫4 𝑥 𝜌𝑦𝑔 𝑠i𝑛30 𝑑𝑥 𝑑𝑦
𝑥* 𝐹 w 8 w 8
= ∫ 𝑥 𝑑𝑥 ∫ 𝜌𝑔𝑦𝑠i𝑛30 𝑑𝑦 = ∫ 𝜌𝑔𝑦 𝑠i𝑛30 . 𝑤𝑑𝑦
𝑅 0 4 2 4
𝑥* 𝐹 = w 𝐹
𝑅 2 𝑅
𝑥* = w = 2.5 m
2
𝐹 𝑅 = ρ g ℎ𝑐 A = ρ g (
2+2sin30) ×4×5
𝑦* = 𝑦
𝑥^ 𝑥^ w 𝑙3 /12
𝑐 + 𝐴𝑦 = 6 + 20×6
𝑐
= 6.22m
𝑥* = 𝐼𝑥^𝑦^
𝑥𝑐 + 𝐴𝑦𝑐
.
𝐼𝑥^𝑦^ = ∫𝐴 𝑥^𝑦^ 𝑑𝐴 =
w 𝑙
∫ w ∫ 𝑙 𝑥^𝑦^𝑑𝑥^𝑑𝑦^ = 0
2 2
− −
2 2
Thus , 𝑥* = 𝑥𝑐 = 2.5 m
𝐹 = volume = 1 (ρgh)hb
𝑅 2
Ex2.5: A pressurised tank contains oil (SG=0.9) and has a square , 0.6 m by 0.6m plate bolted
to its side as shown in fig . The pressure gage on the top of the tank reads 50kpa and the
outside tank is at atmospheric pressure. Find the magnitude & location of the resultant force
on the attached plate .
Soln :1 = (𝑃𝑠 +
ρgℎ1)×0.36 = 24.4 kN
𝐹 =1ρg(ℎ -ℎ )×0.36 =
2 2 2 1
0.954kN
𝐹𝑅 = 𝐹1 + 𝐹2 = 25.4 kN
Ex-2.6
Soln: Basic equations :
𝑑𝑝
= ρg ; |̅𝐹̅̅| = ∫ 𝑝 𝑑𝐴 ;
𝑑ℎ 𝑅
𝐹𝐴R = ∫ 𝑦 𝑑𝐹 = ∫ 𝜌𝑔ℎ𝑦 𝑑𝐴
dA = r d dr ;
y= rsin ; h = H-y
𝐹𝐴 =1 ∫ 𝑀 𝑅 𝑟 𝑠i𝑛2 𝜌𝑔3 (𝐻 − 𝑟𝑠i𝑛) r dr d
∫0 0
𝑅 ( 𝐻𝑟 -𝑟 sin)dr d
𝐹 𝜌g 𝑀
𝐴= ∫0 ∫0
𝑅
𝜌g 𝐻𝑟3 𝑟4
= ∫[ - sin ]𝑅 sin d
0
𝑅 0 3 4
𝜌g 𝐻𝑅3 𝑅4
= ∫( - sin )sin d
𝑅 0 3 4
𝜌g 𝑀 𝐻𝑅3 𝑀 𝑅4
= [∫ sin d - ∫ 𝑠i𝑛2 d ]
𝑅 0 3 0 4
𝜌g 𝐻𝑅3 𝜌g 𝑅4 1 𝑀
= [ −𝑐𝑜𝑠 ]𝑀 - × [ ∫ ( 1 − 𝑐𝑜𝑠2)𝑑 ]
𝑅 3 0 𝑅 4 2 0
3 3
=- 𝜌g 𝐻𝑅
[-1-1] - 𝜌g𝑅 [ - 𝑠i𝑛2 ]𝑀
𝑅 3 8 2 0
2𝜌g𝐻𝑅2 𝜌g𝑅3
= 3 - 8 [ Π]
𝐹 = ρg [ 2𝐻𝑅 2- 𝑀𝑅 3]
𝐴 3 8
𝐹𝐴 = 366 kN . (Ans)
+y = 6m
y = 6 - = 6 - rsin
2𝑀
𝐹 = ρgsin30 ∫ ∫(6 − 𝑟𝑠i𝑛)𝑟 𝑑𝑟 𝑑
𝑅 0 0
𝜌g 2𝑀
= ∫ ∫(6𝑟 − 𝑟2𝑠i𝑛)𝑑𝑟 𝑑
2 0 0
𝜌g 2𝑀 𝑟2 𝑟3 2 𝑅3
𝐹
𝜌g = ∫ [ (6 − 𝑠i𝑛)]𝑅 d = ∫ (3𝑅2 − 𝑠i𝑛) d
𝑅 2 0 2 3 0 2 0 3
3
= 𝜌g [ 3𝑅2 - (-cos) ]2𝑀 𝑅
2 3 0
2𝑀
𝑦* . 𝐹 = ∫ 𝑦 𝑑𝐹 = ∫ ∫( 6 − 𝑟𝑆i𝑛)2 ρgsin dr r d
𝑅 0 0
𝑦* = 6.166m
dA = dr rd
2𝑀 𝑅
𝐼 = ∫ 𝑟2dA = ∫ ∫ 𝑟3dr d
Z^ Z^ 0 0
𝐼 𝑅4
Z^ Z^ 4
×2𝜋
=
𝐼 = 𝑀𝑅4
𝑥^𝑥^ 4
𝑅3
( )[− cos 4𝑅
3
]=0 = 3𝑀
2
(2)
4
𝐼𝑥𝑥 = 𝐼𝑥^𝑥^ +
A𝑦𝑐 2
𝑀𝑅2 4𝑅
𝑀𝑅
4
+ ( )2
=𝐼
𝑥^𝑥^
8 2 3𝑀
Fundamentals of Fluid Mechanics
𝐼𝑥^𝑥^ = 0.1098 𝑅4
30
#Hydrostatic Force on a curved submerged surface:
Consider the curved surface as shown in fig. The pressure force acting on the element of area
, d𝐴̅ is given by
d𝐹̅ = −𝑝d𝐴̅
𝐹̅ = - ∫. pd𝐴̅
𝐴
Where, 𝐹𝑅𝑥 , 𝐹𝑅𝑦 &𝐹𝑅𝑧 are the components of 𝐹̅𝑅̅ in x, y & z directly respectively.
.
𝐹 . = 𝑘^ ̅𝐹̅ = ∫ 𝑑𝐹̅ . 𝑘^ = − ∫ pd𝐴̅ 𝑘^ = 𝑝𝑑𝐴
-∫
𝑅𝑧 𝑅 𝐴 𝐴z 𝑧
Since the direction of the force component can be found by inspection, the use of vectors is
not necessary.
Where d𝐴𝑙 is the projection of the element dA on a plane perpendicular to the ‘l’ direction.
With the free surface at atmospheric pressure, the vertical component of the resultant
hydrostatic force on a curved submerged surface is equal to the total weight of the liquid
above the surface.
Ex:2.9:The gate shown is hinged at ‘O’ and has a constant width w = 5m . The equation of
𝑦2
the surface is x= /𝑎 , where a= 4m . The depth of water to the right of the is D= 4m.Find
the magnitude of the force , 𝐹𝑎 , applied as shown, required to maintain the gate
in equilibrium if the weight of the gate is neglected.
31
Soln: Horizontal Component of force:-
0.5D =2.67m
D
6
Vertical component:
𝐷2
𝐷2 𝐷2
𝐹𝑣 = ∫0𝑎 𝑝𝑤𝑑𝑥 = ∫0𝑎 𝜌𝑔ℎ𝑤𝑑𝑥 = ρgw∫0𝑎 ℎ𝑑𝑥
𝐷2 1 1
. 𝐷2
𝑥* = 𝑥𝑝𝑑𝐴𝑦 = ∫0𝑎 𝑥𝜌𝑔ℎ𝑤𝑑𝑥
𝐹𝑣 ∫𝐴
𝑦
𝐷2 1 1
𝜌gw𝐷5
𝑥* 𝐹𝑣 =0∫ 𝑥(𝐷 − 𝑎 𝑥 )𝑑𝑥 𝑎 2 2
10 𝑎2
=
5
𝑥*= 1
(𝜌gw𝐷 )= 1.2m
𝐹𝑣 10 𝑎2
32
𝐹𝑎= 167kN.
33
Fluids in Rigid-Body Motion:-
Basic equation: −𝑝 + 𝜌𝑔̅ = ρ𝑎̅
6p
− = 𝜌𝑔 (gy= -g) 𝑔̅=-g𝑗^
6y
𝑎𝑥
(𝑑𝑦) 𝑠𝑢𝑟ƒ𝑎𝑐𝑒 = − ( the free surface is a plane)
𝑑𝑥 g
(𝑏/2)
tan = = 𝑏(𝑎𝑥)
𝑒 2 g
e = ( ) = 0.15(𝑎𝑥)
𝑏 𝑎𝑥
{as b=0.3m}
2 g g
34
Thus; 0.3 – d = 0.15 (𝑎𝑥)
g
p =
𝑒 6𝑝 𝑒 6𝑝 6𝑝
𝑟 6𝑟 + 𝑟 6 + 𝑒𝑧 6𝑧
dp = 6𝑝 |
6𝑟 𝑧 dr + 6𝑝 | 𝑟 dz
6𝑧
Taking (𝑟1, 𝑧1) as reference point , where the pressure is 𝑝1 and the arbitrary point (r,z)
where the pressure is p, we can obtain the pressure difference as ;
𝑝 𝑟 6𝑝 6𝑝
∫ 𝑑𝑝 = ∫ 𝑑𝑟 + ∫ dz
𝑝1 𝑟1 6𝑟 6𝑧
2
p𝑝 = ρ (𝑟2 − 𝑟 2) ρg(z-𝑧 )
35
1 2 1 1
36
If we take the reference point at the free surface on the cylinder axis , then;
𝑝1=𝑝𝑎𝑡 ; 𝑟1 =0 and 𝑧1 = ℎ1
2
p 𝑝 = ρ 𝑟2 ρg(zℎ )
𝑎𝑡𝑚 2 1
Since the free surface is a surface of constant pressure (p= 𝑝𝑎𝑡𝑚) , the equation of the free
surface is given by :
2
0 = ρ 𝑟2 ρg(zℎ )
2 1
2
z= + 𝑟2 = ℎ (𝑟)2
ℎ +
1 2g 1 2g
With rotation :
2
= 𝑅
𝑧 2𝜋r (ℎ +
𝑟2) r.dr
∫
∫00
1 2g
2 4
✯= 𝜋[ ℎ1 𝜔 𝑅
𝑅2+ ]
4g
𝑎𝑛𝑑 ℎ1 𝜔2𝑅2
= ℎ0
− 4g
Finally: z =
ℎ
(𝑟)2 1 𝑟
[ − ( ) 2]
0 2g 2 𝑅
Note that this expression is valid only for ℎ1>0 . Hence the maximum value of is given by
𝑚𝑎𝑥 [2gℎ0]1/2
= 𝑅 .
{ (R)2= (ℎ ℎ ) ×4g and 𝜔2 = 1
(ℎ ℎ ) ×4g
0 1 𝑅2 0 1
For ,𝑚𝑎𝑥 ; ℎ1 ≅ 0 }
37
Buoyancy:
When a stationary body is completely submerged in a fluid or partially immersed in a fluid,
the resultant fluid force acting on the body is called the ‘Buoyancy’ force. Consider a solid
body of arbitrary shape completely submerged in a homogeneous liquid.
d𝐹̅1 =p̅𝑑̅𝐴̅
d𝐹𝑉1 = (𝑝𝑎𝑡𝑚 + 𝑝1)d𝐴𝑧 = (𝑝𝑎𝑡𝑚 + 𝜌𝑔ℎ1)d𝐴𝑧
Alternate approach:-
The forces on vertical surfaces are equal and opposite in direction and cancel,
i.e , 𝐹3 − 𝐹4 = 0.
𝐹1 + 𝐹𝐵 + 𝑤ƒ = 𝐹2 or 𝐹𝐵 = 𝐹2 − 𝐹1 − 𝑤ƒ
38
𝐹𝐵 = 𝜌ƒ gℎ2 A - 𝜌ƒ gℎ1 A - 𝜌ƒ g[A(ℎ2 − ℎ1 )-]
𝐹𝐵 = 𝜌ƒ g , where is volume of the body
The direction of the buoyant force, which is the force of the fluid on the body, will be opposite to
that of ‘𝐹𝐵’ shown in fig (FBD of fluid). Therefore, the buoyant force has a magnitude equal to the
weight of the fluid displaced by the body and is directed vertically upward. The line of action of
the buoyant force can be determined by summing moments of the forces w.r.t some
convenient axis. Summing the moments about an axis perpendicular to paper through
point’A’ we have:
𝐹 𝐵 𝑥 = 𝐹 2 𝑥1 − 𝐹 1 𝑥1 − W ƒ 𝑥2
Substituting the forces; we have
𝑥𝐵 =𝑇𝑥1 − (𝑇 − )2
Where 𝑇=A(ℎ2 − ℎ1). The right hand side is the first moment of the displaced volume
and is equal to the centroid of the volume .Similarly it can be shown that the ‘Z’ co-ordinate
of buoyant force coincides with ‘Z’ co-ordinate of the centroid.
𝑥𝐵 =𝑇𝑥1−(𝑇−)𝑥2
39
Stability:-
Another interesting and important problem associated with submerged as well as floating
body is concerned with the stability of the bodies.
When a body is submerged , the equilibrium requires that the weight of the body acting
through its C.G should be collinear with the buoyancy force .However in general, if the body
is not homogeneous in distribution of mass over the entire volume, the location of centre of
gravity ‘G’ don’t coincide with the centre of volume i.e centre of buoyancy, ‘B’ .Depending
upon the relative location of G & B , a floating or submerged body attains different states of
equilibrium , namely (i) Stable equilibrium (ii) Unstable equilibrium (iii) Neutral equilibrium.
40
Stability of submerged Bodies
#Stability problem is more complicated for floating bodies, since as the body rotates the
location of centre of Buoyancy (centroid of displaced volume) may change.
After Displacement, depth of elemental volume immersed is (z+xtan) and the new centre of
Buoyancy 𝑥𝐵 can be expressed as :
41
𝑥𝐵 − 𝑥𝐵= BM tan (as 𝑥𝐵 - 𝑥𝐵 = BM )
GM+BG=
#Notice that is the immersed volume
𝐼𝑦𝑦
GM = - BG
𝐼𝑦𝑦
#Period of oscillation:
From previous discussion we know that restoring couple to bring back the body to its original
equilibrium position is : WGM sin
42
Since the torque is equal to mass moment of inertia ; we can write
WGM sin = - 𝐼 (𝑑 2), where 𝐼
mass M.I of the body about its of rotation.
𝑀 𝑑𝑡2 𝑀
43
2
If ‘’ is small, sin = , and equation can be written as, 𝑑 + 𝖶𝐺.𝑀 = 0 (3)
𝑑𝑡2 𝐼𝑀
Here time period is the time taken for a complete oscillation from one
side to other and back again. The oscillation of the body results in a flow
of the liquid around it and this flow has been neglected here.
Ex-1
A rectangular barge of width b and a submerged depth of H has its centre
of gravity at its waterline. Find the metacentric height in terms of 𝑏 &
𝐻
hence show that for stable equilibrium of the barge 𝑏 ≥ √6 .
𝐻
Soln:
Given that OG = H
OB = , BG = OG-OB = H- 𝐻 = 𝐻
2 2 2
3
BM= 𝐼 = 𝐿𝐵
( Notice that , is the
12×𝐿𝑏𝐻
immersed volume)
2
𝑏
BM=
12𝐻
2
𝑏 1 𝑏
GM=BM-BG= 𝐻 { ( )2 − 1}
−
𝐻 =
12𝐻 2 2 6 𝐻
( 𝑏 ) ≥ √6 proved.
𝐻
44
CHAPTER – 3
INTRODUCTION TO DIFFERENTIAL ANALYSIS OF
FLUID MOTION
Differential analysis of fluid motion:
Integral equations are useful when we are mattered on the gross behaviour of a flow field and
its effect on various devices .However the integral approach doesn’t enable us to obtain
detailed point by point knowledge of flow field.
To obtain this detailed knowledge, we must apply the equations of fluid motion in differential
form.
Rate of change of mass inside the control volume = mass flux in – mass flux out (1)
Mass fluxes:
At left face: ρ u y z
At bottom face: ρ v x z
(𝜌 𝑣 𝑥 𝑧)
At top face: ρ v x z + 6 y
6𝑦
(𝜌 w 𝑥 𝑦)
At back face: ρ w x y + 6 z
6𝑧
45
To find the expression for an incompressible flow:
6𝜌
+ 𝜌 · 𝑢̅→ + 𝑢̅→ · ❑𝜌 = 0
6𝑡
𝑢r𝑒f i L
𝑢𝑟𝑒ƒ 1 𝐷𝜌
=> (❑*. 𝑢*) = −
𝐿 𝜌 𝐷𝑡
1
=> (❑*. 𝑢* ) = −
ur𝑒f 1 𝐷𝜌 ( 4)
( 𝐿 ) · 𝜌 𝐷𝑡
For example: Sea water or stratified air where density varies from layer to layer but the flow
is essentially incompressible as the density of the particles along its path line don’t change.
𝐷𝜌
𝐷𝑡 = 0 , doesn’t necessarily mean that ρ = constant
# If ρ = constant then the flow is incompressible, but the converse is not true, i.e.
Incompressible flow, the density may or may not be constant.
MOMENTUM EQUATION:
A dynamic equation describing fluid motion may be obtained by applying Newton’s 2nd law
to a particle.
46
where the linear momentum ‘P’ is given by:
𝑃̅→𝑠𝑦𝑠𝑡𝑒𝑚 = ∫𝑚𝑎𝑠𝑠𝑉̅→𝑑𝑚
(2)
Then, for an infinitesimal system of mass ‘dm’, Newton’s 2nd law can be written as:
𝑑𝑉̅→
d𝐹→ = 𝑑𝑚 ( ) (3)
𝑑
𝑡
̅→
𝑉
The total derivative
𝑑𝑡 in equation (3) can be expressed as:
Hence;
Now the force d𝐹→ acting on the fluid element can be expressed as sum of the surface
forces ( both Normal forces and tangential forces) and body forces (includes gravity field,
electric field or magnetic fields) .
To obtain the surface forces in x- direction we must sum the forces in x direction. Thus,
47
6𝜎𝑥𝑥 6𝜎𝑦
𝑑𝐹 = (𝜎 + 𝑑𝑥) 𝑑𝑦 𝑑𝑧 − 𝜎 𝑑𝑦 𝑑𝑧 + + ) 𝑑𝑥 𝑑𝑧 − 𝜎 𝑑𝑥 𝑑𝑧 +
(𝜎
𝑠𝑥 𝑥𝑥 6𝑥 𝑥𝑥 𝑦𝑥 6𝑦 𝑦𝑥
On simplifying , we obtain ;
d𝐹 =d𝐹 6𝜎z𝑥
+ d𝐹 =
𝜌g 6𝜎𝑥𝑥 6𝜎𝑦𝑥 + ) dx dy dz (5)
+ +
𝑥 𝑠𝑥 𝑏𝑥 𝑥 (
6𝑥 6𝑦 6𝑧
6𝜎z
d𝐹 6𝜎𝑥𝑦
+ ) dx dy dz (6)
= 𝜌g𝑦 +( 6𝜎𝑦𝑦
𝑦 6𝑥 + 6𝑦 6𝑧
(𝜌g 6𝜎z𝑥
6𝜎𝑥𝑥 + ) = 𝜌 ( 6𝑢 + u 6𝑢
+v 6𝑢 + w 6𝑢 ) (8)
𝑥
+ 𝑦𝑥
6𝜎 +
6𝑥 6𝑦 6𝑧 6𝑡 6𝑦 6𝑦 6𝑧
6𝜎𝑥 6𝜎𝑦𝑦 6𝜎z𝑦
(𝜌g + + + ) = 𝜌 ( 6𝑣 + u 6𝑣 + v 6𝑣 + w 6𝑣
) (9)
𝑦 6𝑥 6𝑦 6𝑧 6𝑡 6𝑥 6𝑦 6𝑧
(𝜌g 6𝜎𝑥z
6𝜎𝑦z
+ + + 6𝜎zz) = 𝜌 (6w + u 6w + v 6w + w 6w ) (10)
𝑧 6𝑥 6𝑦 6𝑧 6𝑡 6𝑥 6𝑦 6𝑧
The stresses may be expressed in terms of velocity gradients & fluid properties in rectangular
co-ordinates as follows :
𝜎 =
𝜎 = 𝜇 ( 6𝑣 + 6𝑢 )
𝑥𝑦 𝑦𝑥 6𝑥 6𝑦
𝜎 =
𝜎 6w 6𝑣
𝑦𝑧 𝑧𝑦 = 𝜇 ( 6𝑦 + 6𝑧 )
𝜎 =
𝜎 = 𝜇 ( 6𝑢 + 6w )
𝑧𝑥 𝑥𝑧 6𝑧 6𝑥
6𝑢
𝜎 = -P - 2 𝜇 ❑ · 𝑉̅→ + 2
𝑥𝑥
𝜎𝑦𝑦
48
𝜎𝑧𝑧 3 6𝑥
6𝑣
= -P - 2 𝜇 ❑ · 𝑉̅→ + 2
3 6𝑦
6w
= -P - 2 𝜇 ❑ · 𝑉̅→ + 2
3 6𝑧
1
𝜎𝑎𝑣 = (𝜎𝑥𝑥 + 𝜎𝑦𝑦 + 𝜎𝑧𝑧 )
3
49
𝜎𝑎 = -P – 2 𝜇 ❑ · 𝑉̅→ + 2 𝜇 ❑ · 𝑉̅→
𝑃𝑚 = 𝑃 − 𝑌 (❑ · 𝑉̅→)
Where ‘P’ is the local thermodynamic pressure, and ‘𝑌’ is co-efficient of bulk viscosity.
If ds is an element of length along the stream line, the equation of streamline is given by:
V ds 0 = iu jv idx jdy kudy vdx
Then we can write: dx dy 0
----------- (1)
x y
Since x, y,t then at any instant t0 , x, y,t0 . Thus at a given instant a change in may be
evaluated as x, y .
Thus at any instant, d dx dy
----------- (2)
x y
50
Comparing Eqn.1 and 2, we see that along an instantaneous streamline
d 0 or is constant
along a streamline. Since differential of is exact, the integral of
d between any two points in a
flow field depends on the end points only, i.e. 2 1 .
The velocity field for a steady, incompressible flow is given as:V Axi Ayj with A=0.3s-1
Determine the stream function that will yield this velocity field. Plot and interpret the streamlines in
the first quadrant of xy plane:
Solution: u Ax
y
df
v Ay .
x dx
Thus, Axy
. The c is arbitrary constant and can be chosen as zero without any loss in
c
generality. With c=0 and A=0.3s-1, we have, Axy . The streamlines in the 1 st quadrant is shown
in Fig.Regions of high speed flow occur where the streamlines are close together. Lower-speed flow
occurs near the origin, where the streamline spacing is wider. The flow looks like flow in a corner
formed by a pair of walls.
51
Before formulating the effects of force on fluid motion (dynamics), let us consider first the
motion (kinematics) of a fluid element on a flow field. For convenience, we follow a
infinitesimal element of a fixed identity (mass)
As the infinitesimal element of mass ‘dm’ moves in a flow field, several things may happen
to it. Certainly the element translates, it undergoes a linear displacement from x,y,z to x1,y1,z1.
The element may also rotate (no change in the included angle in adjacent sides). In addition
the element may deform i.e. it may undergo linear and angular deformation. Linear
deformation involves a deformation in which planes of element that were originally
perpendicular remain perpendicular. Angular deformation involves a distortion of the element
in which planes that were originally perpendicular do not remain perpendicular. In general a
fluid element may undergo a combination of translation, rotation, linear deformation and
angular deformation during the course of its motion.
For pure translation or rotation, the fluid element retains its shape, there is no deformation.
Thus shear stress doesn’t arise as a result of pure translation or rotation (since for a
52
Newtonian fluid the shear stress is directly proportional to the rate of angular deformation).
We shall consider fluid translation, rotation and deformation in turn.
Fluid translation: Acceleration of a fluid particle in a velocity field. A general
description of a particle acceleration can be obtained by considering a particle moving in a
velocity field. The basic hypothesis of continuum fluid mechanics has led us to a field
description of fluid flow in which the properties of flow field are defined by continuous
functions of space and time. In particular, the velocity field is given by
𝑉̅→=𝑉̅→(x,y,z,t). The field description is very powerful, since information for the entire flow is
given by one equation.
The problem, then is to retain the field description for the fluid properties and obtain an
expression for acceleration of a fluid particle as it moves in a flow field. Stated simply,
the problem is:
Given the velocity field 𝑉̅→ =𝑉̅→(x,y,z,t), find the acceleration of a fluid particle, 𝑎̅̅𝑝→ .
Consider the particle moving in a velocity field. At time ‘t’, the particle is at the position x,y,z
and has velocity corresponding to velocity at that point in space at time ‘t’, i.e.
𝑉̅̅̅𝑝→]t=𝑉̅→(x,y,z,t).
At ‘t+dt’, the particle has moved to a new position with co-ordinates x+dx, y+dy, z+dz and
Fig4.1
53
This is shown in pictorial fig 4.1
𝑑̅̅̅̅𝑉𝑝̅→ ,the change in velocity of the particle , in moving from location 𝑟̅̅→ to ̅𝑟̅→ +̅𝑑̅𝑟→ , is given by:
𝑑̅𝑉̅̅̅̅𝑉̅→̅→=
̅6 ̅̅̅6̅̅̅𝑉̅→ ̅6̅𝑉̅→
dx + ̅6̅𝑉̅→
dy + dt
𝑝 dz +
6𝑥 6𝑦 6 p
6𝑡
p
𝑧
p
55
̅𝐷̅𝑉̅→ ̅6̅𝑉̅→ ̅̅̅6̅̅̅𝑉̅→ ̅6̅𝑉̅→
=u +v +w ; which is not necessarily zero.
𝐷𝑡 6𝑥 6𝑦 6𝑧
𝑎𝑥𝑝 = 𝐷𝑢 = u 6𝑢 + v 6𝑢 + w 6𝑢 + 6𝑢 (4.2 a)
𝐷𝑡 6𝑥 6𝑦 6𝑧6𝑡
𝑎𝑦𝑝 = 𝐷𝑣 = u 6𝑣 + v 6𝑣 + w 6𝑣 + 6𝑣 (4.2 b )
𝐷𝑡 6𝑥6𝑦 6𝑧6𝑡
𝑎𝑧𝑝 = 𝐷w = u 6w + v 6w + w 6w + 6w (4.2 c)
𝐷𝑡 6𝑥 6𝑦 6𝑧6𝑡
We have obtained an expression for the acceleration of a particle anywhere in the flow field;
this is the Eularian method of description. One substitutes the coordinates of the point into the
field expression for acceleration.
In the Lagrangian method of description, the motion (position, velocity and acceleration) of a
fluid particle is described as a function of time.
Fluid rotation: A fluid particle moving in a general three dimensional flow field may
rotate about all three coordinate axes. The particle rotation is a vector quantity and in
general
̅m̅→ = ı^ ωx + ^ ωy + 𝑘^ ωz ; where ωx is the rotation about x axis.
To evaluate the components of particle rotation vector, we define the angular velocity
about an axis as the average angular velocity of two initially perpendicular differential line
segments in a plane perpendicular to the axis of rotation.
To obtain a mathematical expression for ωz , the component of fluid rotation about the z axis,
consider motion of fluid in x-y plane. The components of velocity at every point in the field
56
are given by u(x,y) and v(x,y). Consider first the rotation of line segment oa of length Δx.
Rotation of this line is due to the variation of ‘y’ component of velocity. If the ‘y’ component
of the velocity at point ‘o’ is taken as Vo , then the ‘y’ component velocity at point ‘a’ can be
written using Taylor expansion series as:
V = Vo
+6𝑉Δx
6𝑥
ωoa =𝑙i𝑚 ∆𝛼
∆𝑦
∆𝑡→0 Δ𝑡
= 𝑙i𝑚 ∆𝑡→0 Δ𝑥
Δ𝑡
since Δη= ( Va 6𝑣
- Vo ) Δt = ΔxΔt
6𝑥
6𝑣
ωoa =𝑙i𝑚∆𝑡→0 ( )(Δ𝑥Δ𝑡)
6𝑣
6𝑥
Δ𝑥Δ𝑡
= 6𝑥
The angular velocity of ‘ob’ is obtained similarly. If the x- component of velocity at point ‘b’
6𝑢
is uo + ·Δy
6𝑦
∆£
∆𝛽
ωob =𝑙i𝑚 = 𝑙i𝑚 Δ𝑦
∆𝑡→0 Δ𝑡 ∆𝑡→0 Δ𝑡
6𝑢
u = Δy; which will rotate the fluid element in clock-wise direction, thus –ve sign is
b 6𝑦
multiplied to make it counter clock-wise direction.
6𝑢
But Δ£ = − ΔyΔt (-ve sign is used to give +ve value of ω ob )
6𝑦
6u
Thus ω =𝑙i𝑚 −( 6𝑦)(Δ𝑦Δ𝑡) 6𝑢
=−
ob ∆𝑡→0 Δ𝑦Δ𝑡 6𝑦
The rotation of fluid element about z- axis is the average angular velocity of the two mutually
perpendicular line segments, oa and ob, in the x-y plane.
Thus ω
1 6𝑣 6𝑢
z= [
2 6𝑥
− 6𝑦 ]
𝜔̅→ =1 ❑ x 𝑉̅→
2
58
Under what conditions might we expect to have a flow without rotation ( irrotational flow ) ?
A fluid particle moving, without any rotation, in a flow field cannot develop rotation under
the action of body force or normal surface forces. Development of rotation in fluid particle,
initially without rotation, requires the action of shear stresses on the surface of the particle.
Since shear stress is proportional to the rate of angular deformation, then a particle that is
initially without rotation will not develop a rotation without simultaneous angular
deformation. The shear stress is related to the rate of angular deformation through viscosity.
The presence of viscous force means the flow is rotation.
The condition of irrotationality may be a valid assumption for those regions of a flow in
which viscous forces are negligible. (For example , such a region exists outside the boundary
layer in the flow over a solid surface.)
Ç→ =2 𝜔̅→ = ❑ x 𝑉̅→
The circulation, is defined as the line integral of the tangential velocity component about a
closed curve fixed in the flow ; = 𝑉̅→ 𝑑 ̅ 𝑆
̅̅̅ →
𝑐
where 𝑑̅̅̅̅𝑆→ elemental vector tangent to the curve , a positive sense corresponds to a
counter clock-wise path of integration around the curve. A relation between circulation and
vorticity can be obtained by considering the fluid element as shown:
𝑑𝑣 𝑑𝑢
ΔΓ =uΔx +(𝑣 + Δ𝑥) Δ − (𝑢 + Δ𝑦) Δ𝑥 –vΔy
𝑑𝑥 𝑑𝑦
𝑑𝑣 𝑑𝑢
=( − ) Δ𝑥Δ𝑦 = 2ωz Δ𝑥Δ𝑦
𝑑𝑥 𝑑𝑦
Γ = ∆ = 𝑉̅→ ̅𝑑̅̅̅𝑆→
𝑐
=∫𝐴 2𝜔Z dA
=> Γ =∫ (∇x̅ →) dA
𝐴 Z
59
Angular deformation: Angular deformation of a fluid element involves changes in the
perpendicular line segments on the fluid.
We see that the rate of angular deformation of the fluid element in the xy plane is the rated of
decrease of angle “γ” between the line oa and ob. Since during interval Δt,
Δ γ = γ-90 = - ( Δ α+Δ β )
𝑑𝛾 𝑑𝛼 𝑑𝛽
=>- 𝑑𝑡 = 𝑑 + 𝑑𝑡
𝑡
Now;
𝑑𝛼 𝑑𝑣 𝑑𝛽 𝑑𝑢
𝑑 = 𝑑𝑥 and 𝑑 = 𝑑𝑦
𝑡 𝑡
60
INCOMPRESSIBLE INVISCID FLOW
All real fluids posses viscosity. However, in many flow cases it is reasonable to neglect the
effect of viscosity. It is useful to investigate the dynamics of an ideal fluid that is
incompressible and has zero viscosity. The analysis of ideal fluid motion is simpler because
no shear stresses are present in inviscid flow. Normal stresses are the only stresses that must
be considered in the analysis. For a non viscous fluid in motion, the normal stress at a point is
same in all directions (scalar quantity) and equals to the negative of the thermodynamic
pressure, σnn = −P.
6𝑡
In cylindrical co-ordinates:
r: ρ gr 6𝑃
6𝑉 − 𝜃2
− =ρ 6𝑉 + 6𝑉r + 6𝑉r 𝑉
(
r 𝑉𝑟 𝑉𝜃 + r
𝑉𝑧 )
6𝑟 6𝑡 6𝑟 𝑟 6𝜃 6𝑧 𝑟
:ρg
+ 𝑉𝜃 6𝑉𝜃 +
−𝜃
6𝑉
1 6𝑃
=ρ( 𝑉 6𝑉
+ 𝑉 6𝑉𝜃
+
𝑉𝜃 𝑉r
)
61
θ
𝑟 6𝜃 𝑟 6𝑟 6𝑧 𝑟
6𝑡 𝑟 6𝜃
6𝑃 6𝑉z
𝑧: ρ g − =ρ( + 6𝑉 +
𝑉 6𝑉z
𝑉𝜃 + z
𝑉 6𝑉z
)
z
6𝑧 𝑟 6𝑟 6𝑧
6𝑡 𝑟 6𝜃
62
Euler’s equations in streamline co-ordinates:
Applying Newton’s 2nd law in streamwise (the ‘s’) direction to the fluid element of volume
ds x dn x dx, and neglecting viscous forces we obtain:
6𝑃 𝑑𝑠 6𝑃
(𝑃 − ) dn dx − (𝑃 + 𝑑𝑠
)dn dx – 𝜌 g 𝑠i𝑛𝛽 𝑑𝑠 𝑑𝑛 𝑑𝑥 = 𝜌 𝑎 𝑑𝑠 𝑑𝑛 𝑑𝑥
6𝑠 2 6𝑠 2 𝑠
6𝑃
− 6𝑠 − 𝜌 g 𝑠i𝑛𝛽 = 𝜌 𝑎𝑠
6
Since 𝑠i𝑛𝛽 = , we can write:
6𝑠
6 𝐷𝑉̅→ 6𝑉 6𝑉
6𝑃 −𝜌g = 𝜌 =𝜌 ( +𝑉 )
−
6𝑠 6𝑠 𝐷𝑡 6𝑡 6𝑠
=> − 1 6𝑃 −g 6𝑧 = 6𝑉 +𝑉 6𝑉
𝜌 6𝑠 6𝑠6𝑡 6𝑠
where ‘β’ is the angle between ‘n’ direction and vertical and ‘an’ is the acceleration of the
fluid particle in ‘n’ direction.
63
6𝑃
− 6𝑛 − 𝜌 g 𝑐𝑜𝑠𝛽 = 𝜌 𝑎𝑛
6
Since 𝑐𝑜𝑠𝛽 = , we can write:
6𝑛
1 6𝑃 6𝑧
− 6𝑛
−g 6𝑛
= 𝑎𝑛
The normal acceleration of the fluid element is towards the centre of curvature of the
2
streamline; in the negative ‘n’ direction. Thus 𝑎𝑛 = − 𝑉
𝑅
2
=> 1 6𝑃 +g 6𝑧 = 𝑉
6𝑛 6𝑛𝑅
For steady flow on a horizontal plane, Euler’s equation normal to the streamline can be
written as:
2
=> 1 6𝑃 = 𝑉
6𝑛𝑅
Above equation indicates that pressure increases in the direction outward from the centre
of curvature of streamlines.
1 6𝑃 6𝑧 6𝑉
− 6𝑠 −g 6𝑠 = 𝑉 6𝑠
6𝑃
𝑑𝑠 = 𝑑𝑝 (the change in pressure along ‘s’)
6𝑠
6𝑧
𝑑𝑠 = 𝑑𝑧 (the change in elevation along ‘s’)
6𝑠
6𝑉
𝑑𝑠 = 𝑑𝑉 (the change in velocity along ‘s’)
6𝑠
𝑑𝑃
Thus; − − g 𝑑𝑧 = 𝑉𝑑𝑉
𝜌
𝑑𝑃
=> +𝑉𝑑𝑉 + g 𝑑𝑧 = 0
𝜌
64
=> ∫
𝑑𝑃 𝑉2 g𝑧 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛(𝑎𝑙𝑜𝑛𝑔 ′𝑠′) (5.1)
𝜌 + 2 +
For an incompressible flow, i.e. ′𝑃′ is not a function of ′𝜌′; we can write:
𝑃 𝑉2
+ + g𝑧 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛(𝑎𝑙𝑜𝑛𝑔 ′𝑠′)
𝜌 2
Restrictions:
i. Steady flow
ii. Incompressible flow
iii. Inviscid
iv. Flow along a stream line
1 6𝑃 6𝑧 6𝑉 6𝑉
− 6𝑠 −g 6𝑠 = 6𝑡 + 𝑉 6𝑠
Multiplying ds and integrating along a stream line between two points ‘1’ and ‘2’,
2 𝑑𝑝
∫ + 2
𝑉2
2
1
+ g (z2 − z1) + 2 6𝑉 ds =0
−𝑉
∫
1 𝜌 2 1 6𝑡
Restrictions:
i. Incompressible flow
65
ii. Frictionless flow
iii. Flow along a stream line
66
Ex: A long pipe is connected to a large reservoir that initially is filled with water to a depth
of 3 m. The pipe is 150 mm in diameter and 6 m long. Determine the flow velocity leaving
the pipe as a function of time after a cap is removed from its free end.
Assumptions:
i. Incompressible flow
ii. Frictionless flow
iii. Flow along a stream line for ‘1’ and ‘2’
iv. P1 = P2 = Patm
v. V1 =0
vi. Z2=0
vii. Z1=h
viii. Neglect velocity in reservoir, except for small region near the inlet to the tube.
2
Then; g z1 = g ℎ = � 2 + 2 6𝑉 ds (1)
∫ �
2 1 6𝑡
67
Substituting in the equation (1),
68
2
g ℎ = 𝑉2 𝑑𝑉2
+𝐿 𝑑𝑡
2
𝑑𝑉2 𝑑𝑡
2 =
2gℎ−𝑉2 2𝐿
𝑉2 𝑑𝑉2
∫0 2gℎ−𝑉22
1 𝑉 𝑉2
=[ 𝑡𝑎𝑛ℎ−1 ( )] = 𝑡
√2gℎ √2gℎ 2𝐿
1
𝑡𝑎𝑛ℎ−1 ( )=𝑡
√2gℎ √2gℎ 2𝐿
2 𝑡
=> = 𝑡𝑎𝑛ℎ ( √2gℎ)
√2gℎ 2𝐿
69
Consider a displacement in the flow field from position ‘𝑟→’ to ‘ → +d𝑟→’, the displacement ‘d𝑟→’
being an arbitrary infinitesimal displacement in any direction . Taking the dot product of
𝑑𝑟→ =dx ı^ + dy 𝑗^ + dz 𝑘^ with each of the terms, we have
1
− ∇P· d𝑟→ – g 𝑘^ ∙ 𝑑𝑟→ = 1
∇ (V 2 ) · d𝑟→
𝜌 2
𝑑𝑃
And hence − − g 𝑑𝑧 = 1 (𝑉2)
𝜌 2
𝑑𝑃 1
=> + (𝑉2) + g 𝑑𝑧 =0
𝜌 2
2
=> 𝑃 + +
𝑉 g 𝑧 = constant (5.2)
𝜌 2
Since ‘dr→’ was an arbitrary displacement, equation ‘5.2’ is valid between any two points in
a steady, incompressible and inviscid flow that is irrotational.
If ‘dr→’ = ‘ds→’ i.e. the integration is to be performed along a stream line, then taking the
dot product of 𝑑𝑠, we get:
1
(𝑉̅→ ∙ ❑)𝑉̅→ ∙ 𝑑𝑠 = ❑ (𝑉̅→ ∙ 𝑉̅→ ) ∙ 𝑑𝑠 − 𝑉̅→ x (❑ 𝑥 𝑉̅→) · ds
2
Here even though (∇ x 𝑉̅→) is not zero, the product 𝑉̅→ x (∇ x 𝑉̅→) · ds
1. There are significant viscous forces induced by jets, wakes or solid boundaries.
In these cases Bernoulli’s equation will not be valid in such viscous regions.
2. There are entropy gradients caused by shock waves.
3. There are density gradients caused by stratification (uneven heating) rather than by
pressure gradients.
4. There are significant non inertial effects such as earth’s rotation (The
Coriolis component).
70
Principles of a hydraulic Siphon: Consider a container T containing some liquid. If
one end of the pipe S completely filled with same liquid, is dipped into the container with the
other end being open and vertically below the free surface of the liquid in the container T,
then liquid will continuously flow from the container T through pipe S and get discharged at
the end B. This is known as siphonic action and the justification of flow can be explained by
applying the Bernoulli’s equation.
PA 2
PB VB
0 ZA ZB
g g 2g
The pressure at A and B are same and equal to atmospheric pressure. Velocity at A is
negligible compared to velocity at B, since the area of the tank T is very large compared to
that of the tube S. Hence we get,
VB 2g(ZA ZB ) = 2gZ
71
The above expression shows that a velocity head at B is created at the expenses of the
potential head difference between A and B.
PA 2
PC VC
0 ZA ZC
ρg ρg
2g
Considering the pipe cross section to be uniform, we have, from continuity, VB=VC
2
Thus we can PC Patm VB
write; h
ρg ρg 2g
Therefore pressure at C is below atmospheric and pressure at D is the lowest as the potential
head is maximum here. The pressure at D should not fall below the vapor pressure of the
liquid, as this may create vapor pockets and may stop the flow.
72
CHAPTER–4
Laminar flow through a pipe.
Assumptions:
a) Steady
b) Parallel flow in Z- direction Vr = 0 and Vz = u ≠ 0
c) Constant property fluid (ρ & μ are constant)
6
d) Axisymmetric; = 0 Vθ = 0
6𝜃
Continuity Equation:
1 ∂(𝑟𝑉𝑟 ) 1 ∂𝑉𝜃 ∂𝑉𝑧
𝑟 ∂𝑟 + 𝑟 ∂𝜃 + ∂𝑧 = 0
Since 𝑉𝑟 = 0 = 𝑉𝜃; we
have;
∂𝑉𝑧
= 0 ⇒ 𝑉 = 𝑉 (𝑟, 𝜃)
𝑧 𝑧
∂𝑧
6𝑉z
But = 0 (Axisymmetric)
6𝜃
⇒ Vz = Vz(r) = V(r)
73
(P 2πr Δr)z + (2πr Δz τ)r+Δr – (P 2πr Δr)z+Δz – (2πr Δz τ)r = 0
⇒ −( 1 ∂(𝜏 r)
Pz+Δz− Pz =0
Δz ) + 𝑟 ∂𝑟
1
⇒ 𝑑𝑃
(𝑟𝜇 𝑑𝑣 6𝑃 = × ( a constant)
𝑟 𝑑𝑟 𝑑𝑟) = 6𝑧 = 𝑑𝑧
𝑑
𝑑𝑣
⇒ 𝑑𝑟 [𝑟𝜇 ]=𝑟×
𝑑𝑟
𝑑𝑣 2
⇒ 𝑟𝜇 =×𝑟 +𝑐
𝑑𝑟 2 1
𝑣 𝑟 𝑐1
⇒ =× +
𝑑𝑟 2 𝑟
𝑑𝑣 ×𝑟 𝑐1 1
⇒ 𝑑𝑟 = 𝜇 2+ 𝜇 𝑟
× 𝑟2 𝑐1
⇒𝑣 = + ln r + 𝑐
4 𝜇 2
At r = R; V = 0 (No slip B C)
At r = 0; V = finite
× 𝑟2
At r = R; 0 = +
𝜇 4 2
⇒ × 𝑟2
𝑐 =−
2 4
2 ×
2 ×𝑅
⇒𝑣 = ×𝑟 − = (𝑟2 − 𝑅2)
𝜇 4 𝜇 4 4𝜇
𝑑𝑃 2
𝑅
2
⇒ 𝑣 = − ( 𝑑𝑧 ) ( ) [1 − ( ) ]
4𝜇 𝑅
𝑅
H.W- Evaluate 𝑄 = ∫ 𝑣̅ − 𝑑𝐴̅ = ∫0 𝑣 2𝜋𝑟 𝑑𝑟
𝜋𝑅4
⇒𝑄= − 𝑑𝑃
( )
8𝜇 𝑑𝑧
74
Head loss- the friction factor
𝑃1
( ̅𝑣̅1̅2 𝑃2 ̅𝑣̅2̅2
𝜌 + 𝛼 1 + 𝑧 1 ) = ( + 𝛼 2 + 2 ) + ℎƒ
2𝑔 𝜌 2𝑔
𝑔
𝑔
𝛼1 = 𝛼2 and 𝑣̅
̅ 1̅ = 𝑣̅
̅ 2̅ [velocity profile is not changing from 1 to 2 & c.s area is constant]
Thus ℎ = (𝑧 − 𝑧 𝑃2
)+ (1 − ) ----------------------------------------- (1)
ƒ 1 2 𝜌g 𝜌g
4𝑟w 𝐿
ℎƒ
= 𝜌𝑔 𝐷
𝑑𝑃
2 𝑟2
𝑣= −( ) 𝑅 [1
− ]
𝑑𝑧 4𝜇 𝑅2
1 1
𝑣 = 𝑣𝑑𝐴 =
75
2
𝑑𝑃 − (𝑅 ) ( ) 𝑅 2𝜋
(1 −
𝑑𝑟 𝑟 )
2
𝑎𝑣 𝐴
∫ 𝜋𝑅2 𝑑𝑧 4𝜇 ∫0 𝑅2
76
1 𝑑𝑃
𝑅 𝑅4
2
⇒ 𝑣𝑎𝑣 = (− ) [ ]
− 4𝑅2
2𝜇 𝑑𝑧 2
1 𝑑𝑃
𝑅2
⇒ 𝑣̅ = (− )( )
2 𝑑𝑧 4
𝑑𝑣
𝑑𝑃 𝑅 2𝑟 𝑑𝑃 𝑅
2
𝑟w = 𝜇 ] = 𝜇 (− 𝑑𝑧) ( 4 ) [𝑅2] = 𝜇 (− )( )
𝑑𝑟 =𝑅 𝑟=𝑅 𝑑𝑧 2𝜇
1 𝑑𝑃 4𝜇
= 𝑅2 ) ( )𝜇 8𝜇𝑣̅
2𝜇 (− 𝑑𝑧 ) 𝑣̅
4 =
4 = 𝐷
𝑅
𝑅
4𝑟w 𝐿 4 𝐿 8𝜇𝑣̅
ℎƒ =
𝜌𝑔 𝐷 𝜌𝑔 𝐷 𝐷
=
⇒ ℎƒ 32 𝜇 𝐿
𝑣̅ ------------------------------------------------------ (3)
=
𝜌g𝐷2
Long back, Julius Weisback, a German Prof. In 1850, had shown that ℎƒ . Hagen in his
𝐷
experiment had found that ℎƒ 𝖺 𝑣 (approx.). H. Darey a French engineer prposed a
2
𝗌
dimensionless parameter, ‘f’ which is a function of (Red, , duct shape).
𝑑
𝐿
ℎ = ƒ 𝑣̅ 2
ƒ 𝐷 (4)
2g
𝐿𝑣̅2 𝜇 64 𝐿𝑣̅2 64
ℎƒ ( )= ( )
= 2𝑔𝐷 𝜌𝑣̅𝐷 2𝑔𝐷 𝑅𝑒𝐷
78
𝑑𝑃 𝑑𝑧 1 6
⇒− − 𝜌𝑔 + (𝑟𝑟) = 0
𝑑𝑥 𝑑𝑥 6𝑟
𝑑𝑃 𝑑𝑧
⇒(
1
+ 𝜌𝑔 ) =
6
(𝑟𝑟) ------------------------------------------------- (1)
𝑑𝑥 𝑑𝑥 6𝑟
Because of the gravity the local or/and average velocity increases for the above situation i.e
𝑑𝑝 𝑑𝑧
both & are negative.
𝑑𝑥 𝑑𝑥
𝗌
→ Relative roughness
𝑑
Example:- Determine the head loss in friction when water flows at 15 0C through a 300 mm
long galvanized pipe d = 150 mm & Q = 0.05 m 3/s. ν = 1.14 ×10-6 m2/s, = 0.15 mm. Also
find the pumping power required.
79
𝑐𝑣𝐷 𝑣𝐷
Solution:- 𝑅 = = = 3.72×105 f = 0.02
𝑒 𝜇 𝑣
Power = 𝜌𝑄𝑔ℎƒ
Minor losses:- Due to abrupt changes in geometry, shape of the pipes (i.e sudden expansion,
contraction etc.), loss in mechanical energy occurs. In long ducts these losses are very small
compared to the frictional loss, & hence they are termed as minor losses.
(b) Entry & Exit losses (c) Pipe bends (d) Valve & fittings
Total loss = 𝑙ƒ + ℎw
80
𝑄𝐴 = 𝑄𝐵 = 𝑄𝐶
(𝑉𝐴−𝑉𝐵)2
ℎ = & ℎ = 𝑉 2 1
𝑐 ( − 1)
𝑒𝑛 𝑐𝑜𝑛𝑡. 2g 𝐶𝑐
2g
𝑄 = 𝑄𝐴 + 𝑄𝐵
ℎ𝑡 = 𝐻 1 − 𝐻 2 = 𝐿𝐴 𝑣 2 𝐿𝐵 𝑣 2
𝐴 = 𝐵
ƒ𝐴 𝐷 2𝑔 ƒ 𝐷𝐵 2𝑔
𝐵
𝐴
Sudden Enlargement :
From experimental evidence𝑝′ = 𝑝1; where 𝑝′ is the mean pressure of the eddying fluid over
the annular face g-d.
Thus;
𝑝1𝐴1 + 𝑝1(𝐴2 − 𝐴1) − 𝑝2𝐴2 = 𝜌𝑄(𝑣2 − 𝑣1)
81
But; 𝑄 = 𝐴1𝑣1 = 𝐴2𝑣2 (from continuity)
⇒ 𝑝2 − 𝑝1 = 𝜌𝑣2(𝑣1 − 𝑣2)
From SFEE;
𝑝 2 𝑝2 𝑣 2
+𝑣 1 = + 2 + 𝑔ℎ
1
𝜌 2 𝜌 2 2
𝑝2−𝑝1 𝑣
2
⇒ 𝜌 = 1
−𝑣22
− 𝑔ℎ
2 2
2 2
⇒ 𝑣2(𝑣1 − 𝑣2) = 𝑣1 2−𝑣2 − 2
𝑔ℎ
⇒ 2𝑣1𝑣2 − 2𝑣2 = 𝑣2 − 𝑣2 − 2𝑔ℎ2
2 1 2
𝑣2 𝐴 2 𝐴 2 𝑣2 𝐴 2
1 ×
𝑐 2 2
[1 − 𝐴2)] =
𝑐
ℎ2 = 2𝑔
2𝑔 [1 − 𝐴 ]
𝑐
𝐴2
(
𝑐 2
2 2
𝑣22
⇒ ℎ = 𝑣2 (2 ) − 1] [ 1 − 1]
22g [ 𝐴𝑐 2g 𝐶
𝑐
=
82
Where 𝐶 = 𝐴𝑐
= Coefficient of contraction
𝑐 𝐴2
83
# MEASUREMENT OF FLOW RATE THROUGH PIPE:
Flow rates in a pipe are usually measured by providing a co-axial area contraction within the
pipe & by recording the pressure drop across the contraction. The flow rate can be
determined from the pressure drop by straight forward application o0f Bernoulli’s Eqn. Three
such flow meters operate on this principle i.e
1. Venturimeter:
𝖺2<𝖺1
Figure shows a venturimeter inserted in a inclined pipe to measure the flow rate through pipe.
𝑝1 2 𝑝2 2
+ 𝑣1 + 𝑧 = + 𝑣2 + 𝑧
𝜌𝑔 2𝑔 1
𝜌𝑔 2𝑔 2
84
𝑣2−𝑣2 𝑝 −𝑝
𝑝1 𝑝2
⇒ 2 1
= 1 2
+ (𝑧1 − 𝑧2) = ( + 𝑧1) − ( + 𝑧2 ) ----------------------- (1)
2g 𝜌𝑔 𝜌g 𝜌g
⇒ 𝐴
√2𝑔 (𝜌𝜌𝑚 − 1) ℎ
𝑉 1
2 = √𝐴2− 𝐴2
1 2
1 𝐴2 𝜌𝑚
Qth =𝐴 𝑣 = √2𝑔 ( − 1) ℎ--------------------------------------(3)
2 2 𝜌
√𝐴2− 𝐴2
1 2
Measured value of ‘h’, in actual situation will always be greater than that assumed in case of
ideal case due to friction. Thus overestimates the flow rate. To take this into account, a
multiplying factor Cd , is incorporated in equation (3), i.e.,
Qact = 𝐴 𝐴
1 2 √2𝑔 (𝜌𝜌𝑚 − 1) ℎ
𝑑 √𝐴2− 𝐴2
1 2
Value of Cd for venturimeter usually lies between 0.95 to 0.98. It is interesting to note that
‘Q’ remains same whether the pipe is inclined or horizontal.
85
2. Orificemeter:
𝑝 𝑝 2
+ 𝑣12 + 𝑧 = 𝑐 + 𝑣𝑐 + 𝑧
1
𝜌 2𝑔 1
𝜌𝑔 2𝑔 𝑐
𝑔
𝑣2−𝑣2 𝑝 𝑝
⇒ 1
=( 1
+ 𝑧1) − ( 𝑐
+ 𝑧𝑐 )-------------------------------------- (1)
2g 𝜌g 𝜌g
86
2
1
𝖥 𝜌𝑚
[2𝑔 ( − 1) ℎ]1
𝑉𝑐𝑡ℎ =
I [1 − (𝐴 ) 2] I
L 𝐴1
𝑉𝑐𝑎𝑐𝑡 = 𝑉𝑐𝑡ℎ × Cv
1
𝜌 2
[2g( 𝑚−1)ℎ]
𝑄 = Ac =C A =C C A
𝑉 𝑉 [
𝜌
]
𝑎𝑐𝑡 𝑐𝑎𝑐𝑡 c o 𝑐𝑎𝑐𝑡 c v o 𝐴 𝑜𝐶 𝑐 2
[1−( ) ]
𝐴1
1
𝜌 2
[2g( 𝑚−1)ℎ]
⇒ 𝑄𝑎𝑐𝑡 = CdAo [ 𝜌
𝐴 𝐶 2 ]
[1−( 𝑜 𝑐 ) ]
𝐴1
Where Cd = CcCv
87
CHAPTER-5
Irrotational motion refers to fluid flow in which the curl of the velocity field is zero, indicating
that the fluid particles move in such a way that there is no net rotation about any point in the
fluid. Mathematically, this condition can be expressed as:
∇×v=0
Where:
88
The velocity potential function (𝜙ϕ) is a scalar field used to describe irrotational flow. It
represents the velocity potential energy per unit mass associated with fluid motion. The
gradient of the velocity potential function yields the velocity vector field (𝑣v), allowing us to
express the flow velocity in terms of the potential function.
The relationship between the velocity potential function and the velocity vector field is given
by:
v=∇ϕ
Where:
The velocity potential function satisfies Laplace's equation in regions of irrotational flow:
∇2𝜙=0∇2ϕ=0
This equation describes the distribution of the velocity potential throughout the fluid domain
and serves as a fundamental equation in the analysis of irrotational motion.
Irrotational motion is closely related to conservative force fields, particularly in the context of
fluid dynamics. In conservative fields, the work done by a force on an object moving along a
closed path is zero. Similarly, in irrotational flow, the circulation around any closed path
within the fluid is zero, indicating the absence of rotational effects.
Mathematically, the relationship between irrotational motion and conservative force fields can
be understood through the concept of the velocity potential function. In regions where
irrotational flow prevails, the velocity vector field can be expressed as the gradient of the
velocity potential function (𝑣=∇𝜙v=∇ϕ). This gradient field is analogous to the force field
derived from a scalar potential in conservative systems.
Conclusion
In this section, we have explored the fundamentals of irrotational motion in fluid dynamics,
including its definition, characteristics, the derivation and explanation of the velocity potential
function, and its relationship with conservative force fields. Understanding irrotational motion
is essential for simplifying the analysis of fluid flow problems and applying theoretical
principles to practical engineering and scientific applications.
90
CHAPTER 6
MATHEMATICAL FORMULATION OF IRROTATIONAL MOTION
In this section, we will delve into the mathematical formulation of irrotational motion in fluid
dynamics. This will involve examining the Navier-Stokes equations and their simplification
for irrotational motion, introducing Laplace's equation and its role in describing irrotational
flow, and discussing boundary conditions and constraints specific to irrotational motion.
The Navier-Stokes equations are fundamental equations governing the motion of fluid flow.
They describe the conservation of momentum for a fluid element and are expressed in vector
form as:
ρ(∂t∂v+v⋅∇v)=−∇p+μ∇2v+f
Where:
𝑝p is the pressure,
For irrotational motion, the curl of the velocity field (∇×𝑣∇×v) is zero, simplifying the Navier-
Stokes equations. By neglecting the viscous term and external body forces (assuming inviscid
and irrotational flow), the simplified form of the Navier-Stokes equations reduces to:
(∂𝑣∂𝑡+𝑣⋅∇𝑣)=−∇𝑝ρ(∂t∂v+v⋅∇v)=−∇p
This equation is often referred to as the Euler equation and governs the motion of inviscid,
irrotational fluids.
Laplace's equation is a partial differential equation that arises in many areas of science and
engineering, including fluid dynamics. In the context of fluid flow, Laplace's equation
91
describes the distribution of a scalar field (such as the velocity potential) in regions of
irrotational motion. Mathematically, Laplace's equation is given by:
∇2ϕ=0
Where:
Solving Laplace's equation yields the distribution of the scalar field throughout the fluid
domain, providing insight into the behavior of irrotational flow. In particular, the velocity
potential function (𝜙ϕ) can be derived from Laplace's equation, allowing for the
representation and analysis of irrotational motion.
Boundary conditions play a crucial role in defining the behavior of fluid flow, including
irrotational motion. The choice of boundary conditions depends on the specific problem and
the physical characteristics of the system under consideration. In the context of irrotational
motion, some common boundary conditions include:
1. No-slip Boundary Condition: At solid boundaries, the fluid velocity is zero relative to
the boundary (no-slip condition). This boundary condition applies to viscous flows but is
typically not relevant for inviscid, irrotational flows.
2. Far-field Boundary Conditions: At infinity or far from any boundaries, the flow
conditions approach uniform or free-stream conditions. These boundary conditions are
essential for defining the behavior of the flow at large distances from obstacles or
boundaries.
92
ensuring that the scalar field (e.g., velocity potential) and its derivatives satisfy Laplace's
equation throughout the fluid domain, including at the boundaries.
In this section, we have explored the mathematical formulation of irrotational motion in fluid
dynamics. This involved examining the Navier-Stokes equations and their simplification for
irrotational flow, introducing Laplace's equation and its role in describing irrotational flow,
and discussing boundary conditions and constraints specific to irrotational motion.
Understanding these mathematical concepts is crucial for analyzing and solving problems
related to irrotational flow in various engineering and scientific applications.
93
CONCLUSION
In conclusion, the study of irrotational motion and potential flow theory provides valuable
insights into fluid dynamics and its applications across various fields of engineering and
science. By understanding the fundamental principles, mathematical formulations, and
limitations of irrotational motion and potential flow theory, researchers and practitioners can
effectively analyze fluid behavior, design engineering systems, and develop innovative
solutions to real-world problems.
Throughout this research project, we have explored the definition and characteristics of
irrotational motion, the derivation and explanation of the velocity potential function, the
mathematical formulation of irrotational motion from Laplace's equation, and the application
of potential flow theory in diverse areas such as aeronautics, hydrodynamics, acoustics, and
civil engineering. We have also discussed the limitations and assumptions inherent in
potential flow theory, emphasizing the importance of recognizing these constraints for
accurate and reliable analysis.
Overall, the study of irrotational motion and potential flow theory contributes to our
understanding of fluid dynamics and provides a valuable toolset for engineers, scientists, and
researchers to address complex fluid flow problems. By continuing to explore and refine these
theories, we can further advance our capabilities in fluid mechanics and enhance our ability to
design efficient and innovative engineering solutions.
In summary, the study of irrotational motion and potential flow theory is an ongoing endeavor
with vast potential for innovation and discovery. By building upon the foundations laid out in
this research project, we can continue to make significant strides in understanding fluid
behavior and its implications for society and the environment.
REFERENCES:
1. Introduction to Fluid Mechanics by Fox and Mc Donald ,5th edition, Wiley.
94
2. Fluid Mechanics by F.M White, McGrawhill.
10. Landau, L. D., & Lifshitz, E. M. (1987). Fluid Mechanics (2nd ed.).
Pergamon Press.
11. Van Dyke, M. (1975). An Album of Fluid Motion. The Parabolic Press.
12. Panton, R. L. (2013). Incompressible Flow (4th ed.). John Wiley & Sons.
95