Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Lecture 1427495313

Download as docx, pdf, or txt
Download as docx, pdf, or txt
You are on page 1of 102

TABLE OF CONTENT

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

be expressed as :𝑉̅→ = 𝑉̅→(x,y,z,t)


The velocity vector can be written in terms of its three scalar components , i.e
𝑉̅→ =uı^+v𝑗^+w𝑘^
In general ; u = u(x,y,z,t) , v=v(x,y,z,t) and w=w(x,y,z,t)
If properties at any point in the flow field do not change with time , the flow is termed as
steady. Mathematically , the definition of steady flow is 6𝑦 =0 ; where η represents any fluid
6𝑡

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.

One, two and three dimensional flows :


A flow is classified as one two or three dimensional based on the number of space
coordinates required to specify the velocity field. Although most flow fields are inherently
three dimensional, analysis based on fewer dimensions are meaningful.
Consider for example the steady flow through a long pipe of constant cross section
(refer Fig1.6a). Far from the entrance of the pipe the velocity distribution for a laminar flow
𝑢 𝑟
can be described as: = [1 − ( )2]. The velocity field is a function of r only. It is
𝑢𝑚𝑎𝑥
𝑅

independent of r and  .Thus the flow is one dimensional.


Fig1.6a and Fig1.6b
An example of a two-dimensional flow is illustrated in Fig1.6b.The velocity distribution is
depicted for a flow between two diverging straight walls that are infinitely large in z
direction. Since the channel is considered to be infinitely large in z the direction, the velocity
will be identical in all planes perpendicular to z axis. Thus the velocity field will be only
function of x and y and the flow can be classified as two dimensional. Fig 1.7

For the purpose of analysis often it is convenient


to introduce the notion of uniform flow at a given
cross-section. Under this situation the two
dimensional flow of Fig 1.6 b is modelled as one
dimensional flow as shown in Fig1.7, i.e. velocity
field is a function of x only. However,
convenience alone does not justify the assumption such as a uniform flow assumption at a
cross section, unless the results of acceptable accuracy are obtained.
Stress Field:
Surface and body forces are encountered in the study of continuum fluid mechanics. Surface
forces act on the boundaries of a medium through direct contact. Forces developed without
physical contact and distributed over the volume of the fluid, are termed as body forces .
Gravitational and electromagnetic forces are examples of body forces .
Consider an area 𝐴̅ , that passes through ‘c’
.Consider a force 𝐹̅ acting on an area 𝐴̅

through point ‘c’ .The normal stress 𝜎𝑛 and shear

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

Fig 1.9 shows the forces components acting on the area  𝐴 .


The stress components are defined as ;
𝜎 = lim ð𝐹𝑥
𝑥𝑥 𝐴𝑥0 ð𝐴
𝑥

𝜎 = 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 ,𝐹𝑥 .

The shear stress , 𝜎𝑦𝑥 , applied to the fluid element is given by :


𝜎 = lim
𝐹𝑥 𝑑𝐹𝑥
𝑦𝑥 𝐴𝑦0 =
𝑑𝐴𝑦
𝐴𝑦

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 𝑡
= 

𝑑𝑡

To calculate the shear stress, 𝜎𝑦𝑥 𝑑


, it is desirable to express in terms of readily measurable
𝑑𝑡
quantity. l = u t

Also for small angles , l = y 


Equating these two expressions , we have
 𝑢
𝑡 = 𝑦
𝑑
Taking limit of both sides of the expression , we obtain ; = 𝑑𝑢
𝑑𝑡 𝑑𝑦

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𝑒𝑙𝑑



� + μ( )
𝑑𝑦

Examples are : Clay suspension , drilling muds & tooth paste.


Causes of Viscosity:
The causes of viscosity in a fluid are possibly due to two factors (i) intermolecular force of
cohesion (ii) molecular momentum exchange.

#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: 𝑑𝑢
𝑦𝑥 =μ𝑑𝑦

Since the velocity profile is linear ;we have

𝑟𝑦𝑥 −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:

Assumption : linear velocity profile, laminar flow.u = Ω r;


𝑟 =μ𝑑𝑢 = μ 𝗇𝑟; dF= τ dA
𝑦𝑥 𝑑𝑦 ℎ

dF= μ (𝗇𝑟 )2Πr dr



𝑅 2𝑀𝜇𝗇 𝑅 3 𝑀𝜇𝗇𝑅4
T= ∫ 𝑑𝑇= 𝑟 𝑑𝐹 = 𝑟ℎ ∫0 dr =
0 2ℎ
Vapor Pressure:
Vapor pressure of a liquid is the partial pressure of the vapour in contacts with the saturated
liquid at a given temperature. When the pressure in a liquid is reduced to less than vapour
pressure, the liquid may change phase suddenly and flash.

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

(2)Infinitesimal system or differential


analysis
(3) Experimental or dimensional analysis.
In all cases the flow must satisfy three basic laws with a thermodynamic state relation
and associated boundary condition.
1. Conservation of mass (Continuity)
2. Balance of momentum (Newton’s 2nd law)
3. First law of thermodynamics (Conservation of energy)

4. A state relation like ρ=ρ (P, T)

5. Appropriate boundary conditions at solid surface, interfaces, inlets and exits.

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 .

For stream lines : d𝑟̅ ×𝑉̅ = 0


i j 𝑘
|𝑑𝑥 𝑑𝑦 𝑑𝑧| = 0
𝑢 𝑣 𝑤
 ı^( w dy-v dz ) - 𝑗^(w dx –u dz ) + ^𝑘(v dx – u dy ) = 0
 𝑤 𝑑𝑦=v dz ; w dx = u dz & v dx = u dy.
𝑑𝑥 𝑑𝑦 𝑑𝑧
So ; 𝑢 = 𝑣 = w

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

(b)Stream lime plot through ( 𝑥0 , 𝑦0 ,0)

 𝑥0𝑦0 = C
 𝐶 = 16
 xy = 16

(𝑐) 𝑉̅→ = 0.6 ı^– 0.6 𝑗^


𝑑𝑥 𝑥 𝑑𝑥 𝑡
(d) u = Ax , = Ax , ∫ = A∫ 𝑑𝑡
𝑑𝑡 𝑥0 0
𝑥
 ln( ) = At , = 𝑒𝐴𝑡
𝑥

𝑥0 𝑥0
𝑑𝑦 𝑦 𝑑𝑦 𝑡
v = - Ay , = -Ay , ∫ =- A∫ 𝑑𝑡
𝑑𝑡 𝑦0 0
𝑦
 ln( ) = -At ,
𝑦 = 𝑒−𝐴𝑡
𝑦0 𝑦0

At t = 6s ; x=2 𝑒0.3×6 = 12.1 m

; y=8 𝑒−0.3×6 = 1.32 m

(e) 𝑉̅→ = 0.3 ×12.1 ı^ - 0.3 × 1.32 𝑗^ = 3.63 ı^– 0.396 𝑗^


(f)To determine the equation of the path line , we use the parametric equation :

x = 𝑥0 𝑒𝐴𝑡 and y = 𝑦0 𝑒−𝐴𝑡 and eliminate ‘t’


 xy = 𝑥0𝑦0
Remarks :

(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

 x = 1 + 2 (3-2) + 0.5 { (27−8


2 3
) - 4(1 ) }= 3.58 m

(c) The streak line at any given ‘t’ may be obtained by varying ‘𝑡0’ .

# part (a) : path line of particle located at (𝑥0 , 𝑦0) at 𝑡0= 0 s.

𝑡0(s) t X(m) Y(m)


0 0 1 2
0 1 3.08 2.25
0 2 5.67 3.00
0 3 9.25 4.25

#part (b): path lines of a particle located at (𝑥0 , 𝑦0) at 𝑡0 = 2s

𝑡0(s) t(s) X Y
2 2 1 2
2 3 3.58 3.25
2 4 7.67 5.0

#part (c) : 𝑑𝑥𝑢 = 𝑑𝑦𝑣


𝑎𝑦
 dx = ( ) dy
𝑏𝑡
 y dy = 𝑏𝑡 dx
𝑎
 𝑦2 = (2𝑏𝑡) x + c
𝑎

Thus , c = 𝑦 2 – (2𝑏𝑡) 𝑥
0 𝑎 0

For (𝑥0 , 𝑦0) = (1,2) , for different value of ‘t’ .

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.

𝑡0(s) t(s) X(m) Y(m)


0 3 9.25 4.25
1 3 6.67 4.00
2 3 3.58 3.25
3 3 1.0 2.0
CHAPTER – 2
FLUID STATICS
In the previous chapter , we defined as well as demonstrated that fluid at rest cannot sustain
shear stress , how small it may be. The same is true for fluids in “ rigid body” motion.
Therefore, fluids either at rest or in “rigid body” motion are able to sustain only normal
stresses. Analysis of hydrostatic cases is thus appreciably simpler than that for fluids
undergoing angular deformation.

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.

The basic equations of fluid statics :

For a differential fluid element , the body force is 𝑑̅𝐹̅ ̅ = 𝑔̅ dm = 𝑔̅ ρ d


(here , gravity is the only body force considered)where, 𝑔̅ is the local gravity vector ,ρ is
the density & d is the volume of the fluid element. In Cartesian coordinates, d= dx dy dz
.In a static fluid no shear stress can be present. Thus the only surface force is the pressure
force. Pressure is a scalar field, p = p(x,y,z) ; the pressure varies with position within the
fluid.
Pressure at the left face : 𝑃 = ( p - 6𝑝 𝑑𝑦
)
𝐿 6 2

Pressure at the right face : 𝑃 = ( p + 6𝑝 𝑑𝑦


)
𝑅 6 2

Pressure force at the left face : = ( p - 6𝑝 𝑑𝑦


)dx dz 𝑗
𝐿 6 2

Pressure force at the right face : = ( p + 6𝑝 𝑑𝑦


)dx dz (-𝑗^)
𝑅 6 2

Similarly writing for all the surfaces , we have


d̅𝐹→ = ı^ (p - 6 𝑥)dy dz + (p + 6𝑝 𝑑𝑥
)dy dz (-ı^) + ( p - 6𝑝 𝑑𝑦
)dx dz 𝑗
𝑠 6 2 6𝑥 2 6𝑦 2

+ ( p + 6𝑝 𝑑𝑦
)dx dz (-𝑗^) + ( p + 6𝑝 𝑑𝑧
)dx dy (𝑘^ ) +( p + 6𝑝 𝑑𝑧
)dx dy (-𝑘^ )
6 2 6𝑧 2 6𝑧 2

Collecting and concealing terms , we obtain :


6 6𝑝 6𝑝
d̅𝐹→ = - (ı^ + 𝑗^ + 𝑘^ ) dx dy dz
𝑠 6𝑥 6𝑦 6𝑧

 d̅𝐹𝑠→ = - (∇p) dx dy

dz Thus

Net force acting on the body:

 d𝐹→ = d̅𝐹𝑠→ + d̅𝐹̅̅̅ → = ( - ∇p + ρ𝑔̅ ) dx dy dz


 d𝐹→ = ( - ∇p + ρ𝑔̅ )d

or, in a per unit volume basis:


d𝐹
→ = ( - ∇p + ρ𝑔̅ )  (2.1)
𝑑

For a fluid particle , Newton’s second law can be expressed as : d𝐹→ = 𝑎̅ dm = 𝑎̅ ρ dv

Or = 𝑎̅ ρ (2.2)
d𝐹→
d

Comparing 2.1 & 2.2 , we have

- ∇p + ρ𝑔̅ = 𝑎̅ ρ

For a static fluid , 𝑎̅ = 0 ; Thus we obtain ; - ∇p + ρ𝑔̅ =0

The component equations are ; 𝑔̅ = -g 𝑘^


-6𝑝 + ρ𝑔 = 0 𝑔𝑥=0=𝑔𝑦
6𝑥 𝑥
-6𝑝 + ρ𝑔 = 0
6𝑦 𝑦

-6𝑝 + ρ𝑔 = 0
6𝑧 𝑧

Using the value of 𝑔 , 𝑔𝑦 & 𝑔𝑧we have

6𝑝 6𝑝
6𝑝 =0, =0 & = −𝜌 ; since P=P(Z)
6𝑥 6𝑦 6𝑧

𝑑𝑝
We can write = −𝜌𝑔
𝑑𝑧

Restrictions: (i) Static fluid

(ii) gravity is the only body force

(iii) z axis is vertical upward

𝑃g𝑎𝑢g𝑒
𝑃𝑎𝑡𝑚
𝑃𝑎𝑏𝑠

P=0

#Pressure variation in a static fluid :


𝑑𝑃
= -ρg = constant
𝑑𝑧
𝑃 Z

∫ 𝑑𝑃 = - ρg ∫ 𝑑𝑍
𝑃0 Z0

𝑃 − 𝑃0 = - ρg(Z-𝑍0)

𝑃 − 𝑃0 = - ρg(𝑍0 – Z) = ρgh

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:

𝐹𝑧 = σ𝑀Dcos - ρg = 0

Neglecting the volume of the liquid above h , we obtain

 = 𝑀 𝐷2 h
4

Thus ; σ𝑀Dcos - ρg 𝑀 𝐷2 h = 0
4


4𝜎𝑐𝑜𝑠
h = 𝜌g𝐷

Multi Fluid Manometer:


Ex2.2 Find the pressure at ‘A’.

Soln:𝐴 + 𝜌𝑎g ×0.15 - 𝜌𝑚g×0.15 + 𝜌𝑎g ×0.15 - 𝜌wg×0.3 = 𝑃0

#Inclined Tube manometer:


Ex2.3 Given : Inclined–tube reservoir manometer .

Find : Expression for ‘L’ in terms of P.

#General expression for manometer sensitivity

#parameter values that give maximum sensitivity


Soln:

Equating pressures on either side of Level -2 , we have; P = 𝜌𝑙 g (h+H)

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 ; 𝑀 𝐷2𝐻 = 𝑀 𝑑2𝐿


4 4
𝑑

H = L( )2
𝐷

P = 𝜌 g [Lsin + L
𝑑 2
]= 𝜌 gL[ sin + 𝑑 2
] 1

𝑙 (𝐷 ) 𝑙 (𝐷 )
𝑃
Thus, L= 𝑑 2
𝜌𝑙g[sin + ) ]
(
𝐷

To obtain an expression for sensitivity , express P in terms of an equivalent water column


height , ℎ𝑒
2
P= 𝜌wgℎ𝑒

Combining equation 1 &2 , we have

𝜌𝑙gL[ sin + ( 𝑑 )2] = 𝜌w gℎ𝑒


𝐷

𝐿
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
𝐴

Thus , |̅𝐹̅𝑅̅| = 𝑝0 𝐴 + ρg𝑦𝑐 A sin = (𝑝0 + ρg𝑦𝑐 sin)A


Where ℎ𝑐 is the vertical distance between free surface and centroid of the area .
# To evaluate the centre of pressure (c.p) or the point of application of the resultant force
The point of application of the resultant force must be such that the moment of the resultant
force about any axis is equal to the sum of the moments of the distributed force about the
same axis.
If 𝑟̅* is the position vector of centre pressure from the arbitrary origin , then

𝑟̅ * ×̅𝐹̅𝑅̅ = ∫ 𝑟̅ ×d𝐹̅ = - ∫ 𝑟̅ × p d𝐴̅


Referring to fig 2.3 , we can express

𝑟̅* = ı^𝑥* + 𝑗^𝑦*

𝑟̅ = xı^+ y𝑗^ ; d𝐴̅ = - dA 𝑘^ and 𝐹̅𝑅̅ = 𝐹𝑅 𝑘^

Substituting into equation , we obtain

(ı^𝑥* + 𝑗^𝑦* )× 𝑘^ =∫( xı^ + y𝑗^ ) × 𝑑𝐹̅ =∫. ( xı^ + y𝑗^ )× p dA 𝑘^


𝐴
𝐹𝑅
Evaluating the cross product , we get

𝑗^ 𝑥*
𝐹𝑅 + = ∫. ( − 𝑗^ x p + ı^ yp) dA
𝐴
ı^𝑦*𝐹𝑅
Equating the components in each direction ,

𝑦* 𝐹𝑅
= ∫. 𝑦𝑝𝑑𝐴 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𝑛

𝑦* = 𝐼𝑥𝑥

𝐴𝑦𝑐

But from parallel axis theorem , 𝐼𝑥𝑥 = 𝐼𝑥^𝑥^ + A 𝑦𝑐2

Where 𝑥^𝑥^ is the second moment of the area about the centroid al ‘𝑥^’ axis . Thus

𝐼^𝑥^𝑥
𝑦 * = 𝑦𝑐 + 𝐴𝑦𝑐

Or , 𝑦* = ( ℎ𝑐 ) + 𝐼𝑥^𝑥^ 𝑠i𝑛
𝑠i𝑛 𝐴ℎ𝑐

Similarly taking moment about ‘y’ axis ;

𝑥* 𝐹𝑅 = ∫ 𝑥𝑝𝑑𝐴
. .

𝑥* ρg sin 𝑦 A = ∫ 𝑥𝜌𝑔ℎ 𝑑𝐴 = ρgsin∫ 𝑥𝑦𝑑𝐴
𝑐 𝐴 𝐴
.
∫ 𝑥𝑦𝑑𝐴
 *
𝑥 = 𝐴
= 𝐼𝑥𝑦
𝐴𝑦𝑐 𝐴𝑦𝑐

From the parallel axis theorem , 𝐼𝑥𝑥 = 𝐼𝑥^𝑦^ + A𝑥𝑐𝑦𝑐

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.

To find the line of action :

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

Alternative way: By directly using equations:

𝐹 𝑅 = ρ 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

Concept of pressure prism:

𝐹 = 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

If ‘𝐹𝑅 ‘ is the force


acting at a distance 𝑦*
for

the bottom , we have ; 𝐹𝑅 𝑦* = 𝐹1 ×0.3 + 𝐹2 ×0.2 and 𝑦* = 0.296m

Ex-2.6
Soln: Basic equations :
𝑑𝑝
= ρg ; |̅𝐹̅̅| = ∫ 𝑝 𝑑𝐴 ;
𝑑ℎ 𝑅

𝑀̅ =0;Taking moment about the hinge ‘B’ , we have

𝐹𝐴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)

Ex-2.7 :- Repeat the example problem


2.4 if the C.S area of the inclined
surface is circular one , with radius R=2.

Soln: Using integration;


. .
𝐹𝑅 = ∫𝐴 𝑑𝐹 = ∫𝐴 𝜌𝑔ℎ𝑑𝐴 =
𝜌𝑔𝑦𝑠i𝑛 𝑑𝑟 𝑟𝑑

+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

=𝜌g [12×2Π – 0] = 12ρgΠ = 369.458kN


2

Similarly for 𝑦* we can write


Fundamentals of Fluid Mechanics

2𝑀
𝑦* . 𝐹 = ∫ 𝑦 𝑑𝐹 = ∫ ∫( 6 − 𝑟𝑆i𝑛)2 ρgsin dr r d
𝑅 0 0

By using formula : 𝐹𝑅 = ρgℎ𝑐 A = ρg ( 2+2sin30) Π𝑅2 = 369.458kN


44
( )
𝑦 =𝑦 *
=6+ 64
×1
+
𝐼𝑥^𝑥^
𝑐 𝐴𝑦𝑐 44 6
( 4 )

𝑦* = 6.166m

# Find 𝐼𝑥^𝑥^ for a circular C.S

dA = dr rd
2𝑀 𝑅
𝐼 = ∫ 𝑟2dA = ∫ ∫ 𝑟3dr d
Z^ Z^ 0 0

 𝐼 𝑅4
Z^ Z^ 4
×2𝜋
=

But , 𝐼𝑥^𝑥^ + 𝐼𝑦^𝑦^ = 𝐼Z^Z^ (perpendicular axis theorem)


4
2𝐼 = 2𝑀𝑅
𝑥^𝑥^ 4

𝐼 = 𝑀𝑅4
𝑥^𝑥^ 4

# Find 𝐼𝑥^𝑥^ for a semi-circle:


𝑅
∫ 𝑦𝑑𝐴
𝑦𝑐 = ∫ 𝑑𝐴 = ∫0 ∫0 𝑟 𝑠i𝑛 𝑟 𝑑𝑟 𝑑
( 𝑅2)
2

𝑅3
( )[− cos  4𝑅
3
]=0 = 3𝑀
2
(2)
4

𝐼𝑥 = 𝑀𝑅 ( half of the circle)


8
𝑥

𝐼𝑥𝑥 = 𝐼𝑥^𝑥^ +
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𝐴̅
𝐴

We can write; ̅𝐹̅𝑅̅ = ı^𝐹𝑅𝑥 + 𝑗^𝐹𝑅𝑦 + 𝑘^ 𝐹𝑅𝑧

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.

Thus we can write: 𝐹𝑅𝑙 .


= ∫ 𝑝𝑑𝐴𝑙
𝐴
𝑙

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.

𝐹𝑅𝑦 = ∫ 𝑝𝑑𝐴𝑦= ∫ 𝜌𝑔ℎ 𝑑𝐴𝑦 = ∫ 𝜌𝑔𝑑 = ρg

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:-

𝐹𝑅𝐻 = ρgℎ𝑐 (WD) = ρg(0.5) WD = 392kNℎ*


(w𝐷3)
12
h* = ℎ + 𝐼𝑥^𝑥^ = 0.5D +
𝑐 𝐷
𝐴𝑦𝑐 (w𝐷× )
2

 0.5D  =2.67m
D
6

Vertical component:
𝐷2
𝐷2 𝐷2
𝐹𝑣 = ∫0𝑎 𝑝𝑤𝑑𝑥 = ∫0𝑎 𝜌𝑔ℎ𝑤𝑑𝑥 = ρgw∫0𝑎 ℎ𝑑𝑥
𝐷2 1 1

 𝐹𝑣 = 𝜌𝑔𝑤 ∫0 (𝐷 − 𝑎 , (where h+y =D, h = D-y = D-(ax)1/2 )


𝑎2𝑥2)𝑑𝑥
𝐷2
12 3
 𝐹 = 𝜌𝑔𝑤 [ Dx - 𝑎 2 𝑥 ] = (ρgw𝐷3 /3a)
2 𝑎
𝑣 3 0
 𝐹𝑣= 261kN

. 𝐷2
𝑥* = 𝑥𝑝𝑑𝐴𝑦 = ∫0𝑎 𝑥𝜌𝑔ℎ𝑤𝑑𝑥
𝐹𝑣 ∫𝐴
𝑦

𝐷2 1 1
𝜌gw𝐷5
 𝑥* 𝐹𝑣 =0∫ 𝑥(𝐷 − 𝑎 𝑥 )𝑑𝑥 𝑎 2 2
10 𝑎2
=
5
 𝑥*= 1
(𝜌gw𝐷 )= 1.2m
𝐹𝑣 10 𝑎2

Summing moments about ‘O’

𝑀0 = 𝑥* 𝐹𝑣 + 𝐹(𝐷 − ℎ*) − 𝑙𝐹𝑎 =0

32
 𝐹𝑎= 167kN.

33
Fluids in Rigid-Body Motion:-
Basic equation: −𝑝 + 𝜌𝑔̅ = ρ𝑎̅

A fish tank 30cm×60cm×30cm is partially filled with water to be transported in an


automobile. Find allowable depth of water for reasonable assurance that it will not spill
during the trip.

Soln: b=d=30cm= 0.3m


6𝑝 6𝑝
−( ı^+6𝑝 𝑗^ + 𝑘^ ) + ρ +𝑗^ + 𝑘^ )= ρ + + ı^ )
(ı^ 𝑔 (ı^𝑎 ı^𝑎
6𝑥 6𝑦 6𝑧 𝑥 𝑦 𝑧 𝑥 𝑦 𝑧

But; 𝑔𝑥 =0=𝑔𝑧 & 𝑎𝑥 =0=𝑎𝑧


6p
 =0
6z
 p = p(x,y)
6p
− = 𝜌𝑎
6X 𝑥

6p
− = 𝜌𝑔 (gy= -g) 𝑔̅=-g𝑗^
6y

Now we have to find an expression for p(x,y).


6p
dp = 6p 𝑑𝑥 + 𝑑𝑦
6X 6y

But since the force surface is at constant pressure , we have to;


6p
0= 6p 𝑑𝑥 + 𝑑𝑦
6X 6y

𝑎𝑥
 (𝑑𝑦) 𝑠𝑢𝑟ƒ𝑎𝑐𝑒 = − ( the free surface is a plane)
𝑑𝑥 g
(𝑏/2)
 tan = = 𝑏(𝑎𝑥)
𝑒 2 g
 e = ( ) = 0.15(𝑎𝑥)
𝑏 𝑎𝑥
{as b=0.3m}
2 g g

The minimum allowable value of ‘e’ = (0.3 - d )m

34
Thus; 0.3 – d = 0.15 (𝑎𝑥)
g

Hence , = 0.3 – 0.15 (𝑎𝑥)


𝑑 𝑚𝑎𝑥 g

#Liquid in rigid body motion with constant angular speed:


A cylindrical container , partially filled with liquid , is rotated at a constant angular speed ,,
about its axis. After a short time there is no relative motion; the liquid rotates with the
cylinder as if the system were a rigid body .Determine the shape of the free surface.

Soln: In cylindrical co-ordinate;

p =
𝑒 6𝑝 𝑒 6𝑝 6𝑝
𝑟 6𝑟 + 𝑟 6 + 𝑒𝑧 6𝑧

& p + ρg = ρ𝑎̅


−(𝑒 6𝑝
6𝑝 𝑒 6𝑝 ) + ρ(𝑒 𝑔 +
𝑔𝑒 + 𝑒 ) = ρ(𝑒 𝑎 +𝑒 𝑎 +𝑒 𝑎 )
𝑟 6𝑟 + + 𝑒𝑧 𝑟   𝑧 𝑟   𝑧 𝑧
𝑟 6 6𝑧 𝑟 𝑧 𝑟

For the given problem ; 𝑔𝑟 = 𝑔 = 0 & 𝑔𝑧 = −𝑔

and 𝑎 = 𝑎𝑧 = 0 𝑎𝑛𝑑 𝑎𝑟 = −2𝑟

The component equations are:


6𝑝
= 𝜌𝑟r ; 6𝑝=0 and 6𝑝 = ρg
6𝑟 6 6𝑧

Hence , p(r,z) only

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

Volume of the liquid remain constant . Hence  = Π𝑅2ℎ0 ( without rotation)

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𝐴𝑧

d𝐹𝑉2 = (𝑝𝑎𝑡𝑚 + 𝑝2)d𝐴𝑧 = (𝑝𝑎𝑡𝑚 + 𝜌𝑔ℎ2)d𝐴𝑧


The buoyant force (the net force acting vertically upward) acting on the elemental prism is

d𝐹𝐵= (d𝐹𝑉2 − d𝐹𝑉1)= ρg(ℎ2-ℎ1)d𝐴Z = ρgd

Where, d =volume of the prism

Hence, the buoyant force 𝐹𝐵 on the entire submerged body is obtained as :


.
𝐹𝐵 = ∫ pgd , i.e 𝐹𝐵 = ρg

Consider a body of arbitrary shape, having a volume  , is immersed in a fluid. We enclose


the body in a parallelepiped and draw a free body diagram of the parallelepiped with the body
removed as shown in fig. The forces𝐹1, 𝐹2, 𝐹3 &𝐹4 are simply the forces acting on the
parallelepiped, 𝑤ƒ is the weight of the fluid volume (dotted region); 𝐹𝐵 is the force the body
is exerting on the fluid.

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 − 𝑤ƒ

Also; 𝐹1 = ƒ gℎ1 A , 𝐹2 = 𝜌ƒ gℎ2 A and 𝑤ƒ = 𝜌ƒ g[A(ℎ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.

GM=BM – BG , where Metacentric Height


If GM>0 (M is above G) Stable equilibrium
GM=0 (M coincides with G )Neutral Equilibrium
GM<0 (M is below G) Unstable equilibrium

# Theoritical Determination of Metacentric Height:


Before Displacement

𝑥𝐵 = ∫ 𝑥𝑑 = ∫ (𝑧𝑑𝐴) (1)

After Displacement, depth of elemental volume immersed is (z+xtan) and the new centre of
Buoyancy 𝑥𝐵 can be expressed as :

𝑥𝐵  = ∫ (𝑧 + 𝑥 𝑡𝑎𝑛)dA (2)


Subtracting eq.1 from eq.2 , we
have

(𝑥𝐵 − 𝑥𝐵) = ∫ 𝑥2 tan dA = tan ∫

𝑥2dA But ∫ 𝑥2 dA = 𝐼𝑦𝑦

Also, for small angular displacement ;  =tan

41
𝑥𝐵 − 𝑥𝐵= BM tan (as 𝑥𝐵 - 𝑥𝐵 = BM  )

Since ,  BM tan = tan 𝐼𝑦𝑦

 BM =  #Notice that i𝑠 𝑡ℎ𝑒 𝑀. 𝐼 𝑎𝑡 𝑡ℎ𝑒 𝑝𝑙𝑎i𝑛 𝑜ƒ ƒ𝑙𝑜𝑎𝑡𝑎𝑡i𝑜𝑛


𝐼𝑦𝑦 𝐼𝑦𝑦

 GM+BG=

#Notice that  is the immersed volume
𝐼𝑦𝑦

 GM =  - BG
𝐼𝑦𝑦

Fig:Theoritical Determination of Metacentric


Height: #Floating Bodies Containing Liquid:-
If a floating body carrying liquid with free surface undergoes an angular displacement,
the liquid will move to keep the free surface horizontal. Thus not only the centre of
buoyancy moves , but also the centre of gravity ‘G’ moves , in the direction of the
movement of ‘B’.
Thus , the stability of the body is reduced. For this reason, liquid which has to be carried in a
ship is put into a number of separate compartments so as to minimize its movement within
the ship.

#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 𝐼𝑀

Eqn (3) represents an


SHM. 2𝑀 1
𝐼𝑀
2𝑀
The = = 2Π ( )2
= time period, T
1
w W.𝘎𝑀 𝖶.𝐺𝑀
( 𝐼 )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

Also from geometry

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 𝐻

For stable equilibrium of the barge; MG≥ 0


𝐻 𝑏
1
{ ( ) 2 − 1} ≥ 0
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.

Conservation of mass/continuity equation:


The assumption that a fluid could be treated as a continuous distribution of matter – led
directly to a field representation of fluid properties. The property fields are defined by
continuous functions of the space coordinates and time. The density and velocity fields are
related by conservation of mass.

Continuity equation in rectangular co-ordinate system:-

Let us consider a differential control volume of size x, y and z.

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 right face: ρ u y z + 6(𝜌 𝑢 𝑦 𝑧)x


6𝑥

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𝑧

Applying equation (1):

6(𝜌𝑥 𝑦 ∆𝑧) 6(𝜌𝑢) 6(𝜌𝑣) 6(𝜌w)


6𝑡 = 6𝑥 𝑥 𝑦 ∆𝑧  6𝑦 𝑥 𝑦 ∆𝑧  6𝑧 𝑥 𝑦 ∆𝑧

6𝜌 6(𝜌𝑢) 6(𝜌𝑣) 6(𝜌w)


=> 6𝑡 + 6𝑥 + 6𝑦 + 6𝑧 =0
6𝜌
=> +  · (𝜌𝑢̅→ ) = 0 (2)
6𝑡

45
To find the expression for an incompressible flow:
6𝜌
+ 𝜌 · 𝑢̅→ + 𝑢̅→ · ❑𝜌 = 0
6𝑡

=> ( 6𝜌 + 𝑢̅→ · ❑𝜌 )+ ρ· 𝑢̅→ = 0


6𝑡
𝐷
=> + ρ · 𝑢̅→ =0 (3)
𝐷𝑡

Let us define; 𝑢̅→* = 𝑢̅→


; *
= 𝑥i

𝑢r𝑒f i L

𝑢r𝑒f 𝑢r𝑒f 6𝑢i *


. 𝑢̅→ = (❑ * . 𝑢* ) [Since ❑ · 𝑢̅ = 6𝑢i = ]
𝐿 6𝑥i 𝐿 6𝑥i *

𝑢𝑟𝑒ƒ 1 𝐷𝜌
=> (❑*. 𝑢*) = −
𝐿 𝜌 𝐷𝑡
1
=> (❑*. 𝑢* ) = −
ur𝑒f 1 𝐷𝜌 ( 4)
( 𝐿 ) · 𝜌 𝐷𝑡

Eqn (4) may be approximated as (∇*. 𝑢*) = 0


1
If [
ur𝑒f
( 𝐿 )
1 𝐷𝜌
· 𝜌 𝐷𝑡 ]«1 (5)

The velocity field is approximately solenoidal if condition (5) is satisfied.

For incompressible flow, ρ = constant is a wrong statement.(unfortunately such statements


appear in standard books).

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

Hence, for incompressible flow;

· 𝑢̅→ =0, doesn’t matter whether the flow is steady or unsteady.

# 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.

Newton’s 2nd law for a finite system is given by:


𝐹 →
𝑑𝑃̅→ = (1)
)
𝑑𝑡 system

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:

̅6̅𝑉̅→ ̅6̅𝑉̅→ ̅6̅𝑉̅→


u + 6̅𝑉̅→ +w +
v
6𝑥 6𝑦 6𝑧 6𝑡

Hence;

d𝐹→ = 𝑑𝑚 [𝑢 ̅6̅𝑉̅→ ̅6̅𝑉̅


̅6̅𝑉̅→ + + ] (4)
6𝑦 → +̅𝑉̅→
̅6
6𝑥
𝑣 𝑤 6𝑧 6𝑡

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𝑦 𝑦𝑥

(𝜎𝑧𝑥 6𝜎z𝑥 ) dx dy - 𝜎𝑧𝑥 dx dy


+ 6𝑧

On simplifying , we obtain ;

d𝐹𝑠𝑥= 6𝜎𝑥 6𝜎z𝑥


6𝜎𝑦𝑥 + ) dx dy dz
(𝑥6𝑥 + 6𝑦 6𝑧

d𝐹 =d𝐹 6𝜎z𝑥
+ d𝐹 =
𝜌g 6𝜎𝑥𝑥 6𝜎𝑦𝑥 + ) dx dy dz (5)
+ +
𝑥 𝑠𝑥 𝑏𝑥 𝑥 (
6𝑥 6𝑦 6𝑧

Similar expression for the force components in y & z direction are:

6𝜎z
d𝐹 6𝜎𝑥𝑦
+ ) dx dy dz (6)
= 𝜌g𝑦 +( 6𝜎𝑦𝑦
𝑦 6𝑥 + 6𝑦 6𝑧

d𝐹𝑥= 𝜌g𝑧 + 6𝜎𝑥z


( 6𝜎𝑦z ) dx dy dz (7)
6𝑥 + 6𝑦 +6𝜎zz
6𝑧

Now writing the differential form of equation of motion:

(𝜌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𝑧

Newtonian fluid :- Navier-stokes equation:

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.

Stream function for two dimensional incompressible flow:


It is convenient to have a means of describing mathematically any particular pattern of flow.
A mathematical device that serves this purpose is the stream function,  . The stream function
is formulated as a relation between the streamlines and the statement of conservation of mass.
The stream function  (x, is a single mathematical function that replaces two velocity
y,t)
components, u(x, y,t) and v(x, y,t) .
For a two dimensional incompressible flow in the xy plane, conservation of mass can be
u v
written as :  0.
x y

If a continuous function  (x, 


y,t) called stream function is defined such that u  and
y

v , then the continuity equation is satisfied exactly.
x
u v  2  2
Then     0 and the continuity equation is satisfied exactly.
x x y
y x
y


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 kudy  vdx

Thus equation of streamline in a two dimensional flow is: udy  vdx 0

 
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 .

Example problem: Stream Function flow in a corner:


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

Integration with respect to y yields:



 dy  f x
=  Axy  f x ;
y

where f(x) is an arbitrary function of x.

f(x) can be evaluated using the expression for v.


Thus we can write,

 df
v  Ay .
x dx

But from the velocity field description, v  Ay.Hence df


 0 or f(x) =constant.
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

has a velocity given by: ̅𝑉̅𝑝→]t+dt=𝑉̅→(x+dx,y+dy,z+dz,t+dt).

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

The total acceleration of the particle is given by :

̅𝑑̅𝑉̅𝑝→ ̅6̅𝑉̅→ 𝑑𝑥 ̅̅̅6̅̅̅𝑉̅→ 𝑑𝑦𝑝 ̅6̅𝑉̅→ ̅6̅𝑉̅→


𝑎̅̅̅̅→𝑝𝑝=
𝑑𝑧 = + + +
𝑑𝑡 6𝑥 𝑑𝑡 6𝑦 𝑑𝑡 6 𝑑𝑡 6𝑡
𝑧
𝑑𝑥 𝑑𝑦𝑝 𝑑𝑧𝑝
Since = u, = v and = w,
𝑑𝑡 𝑑𝑡 𝑑𝑡

̅𝑑̅𝑉̅𝑝→ ̅6̅𝑉̅→ ̅̅̅6̅̅̅𝑉̅→ ̅6̅𝑉̅→


𝑎̅̅̅̅𝑉̅→𝑝̅→=
̅6 =u +v +w +
𝑑𝑡 6𝑥 6𝑦 6𝑧 6𝑡

̅𝐷̅𝑉̅→ 𝑑̅𝑉̅𝑝→ ̅6̅𝑉̅→ ̅̅̅6̅̅̅𝑉̅→ ̅6̅𝑉̅


̅6̅𝑉̅→ = ̅𝑎𝑝̅̅̅→= =u +v +w +→ (4.1)
𝐷𝑡 𝑑𝑡 6𝑥 6𝑦 6 6𝑡
𝑧
̅𝐷̅𝑉̅→
The derivative is commonly called substantial derivative to remind us that it is computed
𝐷𝑡

for a particle of substance. It is often called material derivative or particle derivative.


From equation 4.1 we recognize that a fluid particle moving in a flow field may undergo
acceleration for either of the two reasons. It may be accelerated because it is convected into
a region of higher (lower) velocity. For example, the steady flow through a nozzle, in which
by definition, the velocity field is not a function of time, a fluid particle will accelerate as it
moves through the nozzle. The particle is convected into a region of higher velocity. If a flow
field is unsteady the fluid particle will undergo an additional “local” acceleration, because
the velocity field is a function of time.
The physical significance of the terms in the equation 4.1 is :
̅6̅𝑉̅→ ̅̅̅6̅̅̅𝑉̅→ ̅6̅𝑉̅→
u +v +w = convective acceleration
6𝑥 6𝑦 6𝑧
̅6̅𝑉̅→
= local acceleration.
6𝑡

Therefore equation 4.1 can be written as:


̅𝐷̅𝑉̅→ ̅6̅𝑉̅→
𝑎̅̅̅̅→ = = ̅(̅𝑉̅→·❑) 𝑉̅→ +
𝑝 𝐷𝑡 6𝑡
54
For a steady and three dimensional flow the equation 4.1 becomes:

55
̅𝐷̅𝑉̅→ ̅6̅𝑉̅→ ̅̅̅6̅̅̅𝑉̅→ ̅6̅𝑉̅→
=u +v +w ; which is not necessarily zero.
𝐷𝑡 6𝑥 6𝑦 6𝑧

Equation 4.1 may be written in scalar component equation as:

𝑎𝑥𝑝 = 𝐷𝑢 = 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𝑦 ]

By considering the rotation about other axes :


6w 6𝑣
ω =1[ − ] and ω
1 6𝑢 6w
x
2 6𝑦 6𝑧
y=
2
[6𝑧 − 6𝑥 ]
6w 6𝑣 6𝑢 6w 6𝑣 6𝑢
Then ̅→ =1 [( − ) ı^ + ( − ) 𝑗^ + ( − ) 𝑘^ ] ; which can be written in
2 6𝑦 6 6𝑧 6𝑥 6𝑥 6𝑦
𝑧
57
vector notation as :

𝜔̅→ =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.)

A term vorticity is defined as twice of the rotation as:

Ç→ =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.

Momentum equation for frictionless flow: Euler’s equations:


The equations of motion for frictionless flow, called Euler’s equations, can be obtained from
the general equations of motion, by putting μ = 0 and σnn = -p.
6𝑃 6
6𝑢 −
ρg =ρ( + 𝑤 )
6𝑢 6𝑢
x
+𝑢 +𝑣
6𝑥 6𝑡 6𝑥 6𝑦 6𝑧
6𝑃 6𝑣 6
ρg − =ρ( + 𝑤 )
6𝑣 6𝑣
y +𝑢 +𝑣
6𝑦 6𝑡 6𝑥 6𝑦 6𝑧
6𝑃 6w 6w
ρg − =ρ( + 𝑤 )
6w 6w
z +𝑢 +𝑣
6𝑧 6𝑡 6𝑥 6𝑦 6𝑧

In vector form it can be written as:


6𝑉̅→
ρg̅→ − ∇P = ρ + )
̅→ ( 6𝑉̅→𝑢
6𝑉̅→
6𝑉 + +𝑣
6𝑡 6𝑥 6𝑦 6𝑧

=> ρg̅→ − ∇P= ρ


̅→ ( + (𝑉̅→ ∇)𝑉̅→ )
6𝑉

6𝑡

=> ρg̅→ − ∇P= ρ 𝐷𝑉̅→


𝐷𝑡

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 𝑠

Simplifying the equation we have:

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𝑠

To obtain Euler’s equation in a direction normal to the streamlines, we apply Newton’s 2 nd


law in the ‘n’ direction to the fluid element. Again, neglecting viscous forces; we obtain:
6 𝑛 6𝑃
(𝑃 − ) ds dx − (𝑃 + 𝑑𝑛
)ds dx – 𝜌 g 𝑐𝑜𝑠𝛽 𝑑𝑛 𝑑𝑥 𝑑𝑠 = 𝜌 𝑎 𝑑𝑛 𝑑𝑥 𝑑𝑠
6𝑛 2 6𝑛 2 𝑛

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.

Bernoulli’s equation: Integration of Euler’s equation along a stream line for


steady flow( Derivation using stream line co-ordinates):
Euler’s equation for steady flow will be:

1 6𝑃 6𝑧 6𝑉
− 6𝑠 −g 6𝑠 = 𝑉 6𝑠

If a fluid particle moves a distance ‘ds’ along a streamline, then

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

* In general the constant has different values along different streamlines.

* For derivation using rectangular co-ordinates, refer page-7.

Unsteady Bernoulli’s equation( Integration of Euler’s equation along a stream line):


1
− ❑P− g̅→ 𝐷𝑉̅→
= 𝐷𝑡 or
𝜌

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𝑡

For an incompressible flow, the above equation reduces to :


𝑃1 2 2
+ �1 + g z1 = 𝑃+2 �2 + g z2 + 2 6𝑉
� �∫ ds
𝜌 2 𝜌 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.

Ans: Applying Bernoulli”s equation between 1 and 2 we have:


𝑃1 2 2
+ �1 + g z1 = 𝑃+2 �2 +g z2 + 2 6𝑉
� �∫ ds
𝜌 2 𝜌 2 1 6𝑡

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𝑡

In view of assumption ‘viii’, the integral becomes


2 6𝑉 𝐿 6𝑉
∫ ds ≈ ∫ ds
1 6𝑡 0 6𝑡

In the tube, V = V2, everywhere, so that


𝐿 6𝑉 𝐿 𝑑𝑉2 𝑑𝑉2
∫ ds = ∫ ds = L
0 6𝑡 0 𝑑𝑡 𝑑𝑡

67
Substituting in the equation (1),

68
2
g ℎ = 𝑉2 𝑑𝑉2
+𝐿 𝑑𝑡
2

Separating the variables we obtain:

𝑑𝑉2 𝑑𝑡
2 =
2gℎ−𝑉2 2𝐿

Integrating between limits V = 0 at t = 0 and V = V2 at t = t,

𝑉2 𝑑𝑉2
∫0 2gℎ−𝑉22
1 𝑉 𝑉2
=[ 𝑡𝑎𝑛ℎ−1 ( )] = 𝑡
√2gℎ √2gℎ 2𝐿

Since 𝑡𝑎𝑛ℎ−1(0) = 0, we obtain

1
𝑡𝑎𝑛ℎ−1 ( )=𝑡
√2gℎ √2gℎ 2𝐿

2 𝑡
=> = 𝑡𝑎𝑛ℎ ( √2gℎ)
√2gℎ 2𝐿

Bernoulli’s equation using rectangular coordinates:


1
− ∇P – g 𝑘^ = (𝑉̅→ ∙ ∇)𝑉̅→
𝜌

Using the vector identity:


1
(𝑉̅→ ∙ ∇)𝑉̅→ = ∇ (𝑉̅→ ∙ 𝑉̅→) − 𝑉̅→ x (∇ x 𝑉̅→)
2

For irrotational flow: ∇ x 𝑉̅→ = 0


1
So (𝑉̅→ ∙ ∇)𝑉̅→ = ∇ (𝑉̅→ ∙ 𝑉̅→ )
2
1
− ∇P – g 𝑘^ = 1
∇ (𝑉̅→ ∙ 𝑉̅→ ) = 1
∇ (V 2 )
𝜌 2 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

will be zero as 𝑉̅→ x (∇ x 𝑉̅→) is perpendicular to V and hence perpendicular to ds.

# A fluid that is initially irrotational may become rotational if:-

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).

HGL and EGL:


Hydraulic Grade Line (HGL) corresponds to the pressure head and elevation head i.e. Energy
Grade Line(EGL) minus the velocity head.
𝑃
EGL = 𝑉2
𝜌g
+ + 𝑧 =H (Total Bernoulli’s constant)
2g

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.

Applying the Bernoulli’s equation between point A and B, we can write

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 ) = 2gZ

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.

Applying the Bernoulli’s equation between point A and B, we can write

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)

Consider a differential annular control volume:

Applying the force balance in Z-direction, we have

73
(P 2πr Δr)z + (2πr Δz τ)r+Δr – (P 2πr Δr)z+Δz – (2πr Δz τ)r = 0

⇒ (Pz – Pz+Δz) 2πr Δr + 2π Δz [(τ r) r+Δr - (τ r) 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

The RHS of the equation will be finite only if C1 = 0.


× 𝑟2
Thus; 𝑣 = + 𝑐
𝜇 4 2

× 𝑟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

The SFEE between (1) & (2) gives;

𝑃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

Applying the momentum relation to the control volume

(𝑃1𝜋𝑅2 − 𝑃2𝜋𝑅2) − 𝑟w 2𝜋𝑅𝐿 + 𝜌𝑔 (𝜋𝑅2𝐿) sin 𝜃 = 𝑚˙ (𝑣̅ ̅ 2̅) = 0


̅ 1̅ = 𝑣̅
𝑃1−𝑃2
⇒ +
(𝑧 − 𝑧 2𝑐w 4𝑐w 𝐿 ------------------------------------ (2)
1 2
)= =
𝜌g 𝜌g 𝑅 𝜌g 𝐷

Comparing eqn. (1) & (2), we have;

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

Rewriting Eqn. (3) in form of Eqn. (4), we have

𝐿𝑣̅2 𝜇 64 𝐿𝑣̅2 64
ℎƒ ( )= ( )
= 2𝑔𝐷 𝜌𝑣̅𝐷 2𝑔𝐷 𝑅𝑒𝐷

(𝑃𝑥 − 𝑃𝑥+Δ𝑥)2𝜋𝑟 Δ𝑟 + 2𝜋 Δ𝑥 [(𝑟𝑟)𝑟+Δ𝑟 − (𝑟𝑟)𝑟] + 𝜌(2𝜋𝑟 Δ𝑟 Δ𝑥)𝑔 sin 𝜃 = 0


77
Dividing by 2𝜋𝑟 Δ𝑟 Δ𝑥
⇒ lim 𝑃𝑥− 𝑃𝑥+Δ𝑥 1 6 −Δ𝑧
)+ ( ) + lim 𝜌𝑔 ( ) = 0
Δ𝑥→0 ( Δ𝑥 6𝑟
𝑟𝑟 Δ𝑥
Δ𝑥→0

78
𝑑𝑃 𝑑𝑧 1 6
⇒− − 𝜌𝑔 + (𝑟𝑟) = 0
𝑑𝑥 𝑑𝑥 6𝑟

𝑑𝑃 𝑑𝑧
⇒(
1
+ 𝜌𝑔 ) =
6
(𝑟𝑟) ------------------------------------------------- (1)
𝑑𝑥 𝑑𝑥 6𝑟

Let 𝑃 = 𝑝 + 𝜌𝑔𝑧 where P→ modified pressure


𝑑𝑃
= 𝑑𝑝 + 𝜌𝑔 𝑑𝑧
𝑑𝑥 𝑑𝑥 𝑑𝑥
Putting in Eqn. (1), we can write;
𝑑𝑃
= 1 ∂ (𝑟𝑟)
𝑑𝑥 𝑟 ∂𝑟
Thus; 𝑑𝑃 2 𝑟2
𝑣 = − ( ) (𝑅 ) [1 − ]
𝑑𝑥 4𝜇 𝑅2

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 = 𝜌𝑄𝑔ℎƒ

The head lost due to friction is called major loss.

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.

The minor head losses may be expressed as ̅2


= where K is determined
ℎƒ 2
experimentally.

(a) Sudden contraction & Enlargements

(b) Entry & Exit losses (c) Pipe bends (d) Valve & fittings

Total loss = 𝑙ƒ + ℎw

Four cases for solving pipe problems:-

a) L, Q & D known, ΔP unknown


b) ΔP, Q & D known; L unknown
c) ΔP, L & D known; Q unknown
d) ΔP, L & Q known, D unknown

Flow through Branched pipes:-

(1) Pipes in series:

80
𝑄𝐴 = 𝑄𝐵 = 𝑄𝐶

ℎ𝑡 = ℎƒ𝐴 + ℎ𝑒𝑛 + ℎƒ𝐵 + ℎ𝑐𝑜𝑛𝑡. + ℎƒ𝑐


ƒ𝑙𝐴𝑣2
Where ℎƒ𝐴 = 2g𝐷 𝐴 & so on, other pipes.
𝐴

(𝑉𝐴−𝑉𝐵)2
ℎ = & ℎ = 𝑉 2 1
𝑐 ( − 1)
𝑒𝑛 𝑐𝑜𝑛𝑡. 2g 𝐶𝑐
2g

(2) Pipes in parallel:-

𝑄 = 𝑄𝐴 + 𝑄𝐵

ℎ𝑡 = 𝐻 1 − 𝐻 2 = 𝐿𝐴 𝑣 2 𝐿𝐵 𝑣 2
𝐴 = 𝐵
ƒ𝐴 𝐷 2𝑔 ƒ 𝐷𝐵 2𝑔
𝐵
𝐴
Sudden Enlargement :

𝑝1𝐴1 + 𝑝′(𝐴2 − 𝐴1) − 𝑝2𝐴2 = 𝜌𝑄(𝑣2 − 𝑣1)

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 = 𝜌𝐴2𝑣2(𝑣1 − 𝑣2)

⇒ 𝑝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 = (𝑣1 − 𝑣2)2


(𝑣1−𝑣2)2
2 2
⇒ ℎ2 = = 𝑣1 [1 − ( 𝐴1 )]
2g 2g 𝐴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

(i) Venturimeter (ii) Orificemeter (iii) Flow nozzle

1. Venturimeter:

𝖺2<𝖺1

Figure shows a venturimeter inserted in a inclined pipe to measure the flow rate through pipe.

Let us consider a steady, ideal and one dimensional flow of fluid.

Applying Bernoulli’s Equation.:

𝑝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

From pressure balance at section 0-0:-

𝑝1 + 𝜌(𝑧1 − ℎ𝑜) = 𝑝2 + 𝜌𝑔(𝑧2 − ℎ − ℎ𝑜) + 𝜌𝑚𝑔ℎ


𝑝1 𝑝2 1
⇒( +𝑧 )= ( + 𝑧 ) + (𝜌 − 𝜌)ℎ ×
𝜌g 1 𝜌g 2 𝑚 𝜌
𝑝1 𝑝2 1
⇒( + 𝑧1 ) − ( ( )ℎ ×-----------------------------------------(2)
𝜌g + 𝑧2 ) = 𝜌𝑚 − 𝜌 𝜌
𝜌g

Putting the above value in Eqn. (1);


2
𝑣 − 𝑣2
2 1 = (𝜌 − 𝜌 )ℎ × 1
2𝑔 𝑚
𝜌

From continuity; 𝐴1𝑣1 = 𝐴2𝑣2


𝐴2𝑣2
𝑣1 =
𝐴1
𝜌𝑚
𝐴2 2 − 1) ℎ
Thus; 2 2
𝑉2 − (𝐴 ) 𝑉 = 2𝑔 ( 𝜌
1
2
𝜌𝑚
[2g( −1)ℎ]
𝜌
⇒ 𝑉22 = 𝐴2 2
[1−( ) ]
𝐴1

⇒ 𝐴
√2𝑔 (𝜌𝜌𝑚 − 1) ℎ
𝑉 1
2 = √𝐴2− 𝐴2
1 2

1 𝐴2 𝜌𝑚
Qth =𝐴 𝑣 = √2𝑔 ( − 1) ℎ--------------------------------------(3)
2 2 𝜌
√𝐴2− 𝐴2
1 2

The above value is the theoretical discharge/flow rate.

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:

Cc = 𝐴 ; where Ao is the area of the orifice.


𝐴𝑜

Applying Bernoulli’s Equation between 1 and c.:

𝑝 𝑝 2
+ 𝑣12 + 𝑧 = 𝑐 + 𝑣𝑐 + 𝑧
1
𝜌 2𝑔 1
𝜌𝑔 2𝑔 𝑐

𝑔
𝑣2−𝑣2 𝑝 𝑝
⇒ 1
=( 1
+ 𝑧1) − ( 𝑐
+ 𝑧𝑐 )-------------------------------------- (1)
2g 𝜌g 𝜌g

From pressure balance at section 0-0:-

𝑝1 + 𝜌(𝑧1 − ℎ𝑜) = 𝑝2 + 𝜌𝑔(𝑧2 − ℎ − 𝑧𝑜) + 𝜌𝑚𝑔ℎ


𝑝1 𝑝2 1
⇒( +𝑧 )= ( + 𝑧 ) + (𝜌 − 𝜌)ℎ ×
𝜌g 1 𝜌g 2 𝑚 𝜌
𝑝1 𝑝2 1
⇒( + 𝑧1 ) − ( ( )ℎ ×-----------------------------------------(2)
𝜌g + 𝑧2 ) = 𝜌𝑚 − 𝜌 𝜌
𝜌g

Putting the above value in Eqn. (1);


𝑣 2
− 𝑣2
𝑐 1 = (𝜌 − 𝜌 )ℎ × 1
2𝑔 𝑚
𝜌

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

Orificemeters are less accurate than venturimeters.

87
CHAPTER-5

FUNDAMENTALS OF IRROTATIONAL MOTION


Irrotational motion is a concept central to the study of fluid dynamics, describing the flow of a
fluid where the fluid particles move without causing any net rotation about any point in the
fluid. In this section, we will delve into the definition and characteristics of irrotational
motion, explore the derivation and explanation of the velocity potential function, and discuss
the relationship between irrotational motion and conservative force fields.

Definition and Characteristics of Irrotational Motion

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:

 ∇× denotes the curl operator,

 v represents the velocity vector field.

Key characteristics of irrotational motion include:

1. Conservation of Angular Momentum: In irrotational flow, the angular momentum of


fluid particles about any point remains constant. This property simplifies the analysis of
fluid motion, as it eliminates the need to consider the effects of rotation.

2. Potential Existence: In regions of irrotational flow, it is possible to define a scalar field


known as the velocity potential function (𝜙ϕ). The gradient of this potential function
yields the velocity vector field, providing a convenient way to represent and analyze fluid
flow.

3. Streamline Behavior: Irrotational flow is characterized by smooth streamline patterns,


with the streamlines representing paths along which fluid particles move. These
streamlines do not intersect, and the flow remains continuous and free from vortices.

Derivation and Explanation of the Velocity Potential Function

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:

 ∇∇ represents the gradient operator,

 𝜙ϕ denotes the velocity potential function.

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.

Relationship between Irrotational Motion and Conservative Force Fields

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.

Furthermore, in potential flow theory, which is based on the assumption of irrotational


motion, the fluid is considered inviscid and incompressible, and the flow is governed by
Laplace's equation. This simplification allows for the application of principles from potential
theory, such as the superposition of elementary solutions, to analyze complex fluid flow
89
problems.

In summary, irrotational motion is intimately connected to conservative force fields through


the velocity potential function, providing a mathematical framework for analyzing and
understanding fluid flow phenomena.

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.

Navier-Stokes Equations and their Simplification for 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:

 𝜌ρ is the fluid density,

 𝑣v is the velocity vector field,

 𝑝p is the pressure,

 𝜇μ is the dynamic viscosity,

 𝑓f represents external body forces.

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.

Introduction to Laplace's Equation and its Role in Describing Irrotational Flow

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:

 𝜙ϕ is the scalar field (e.g., velocity potential),

 ∇2∇2 denotes the Laplacian operator.

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 and Constraints for 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.

3. Prescribed Velocity or Pressure Boundary Conditions: In some cases, the velocity or


pressure may be specified at certain boundaries or interfaces within the fluid domain.
These boundary conditions provide additional constraints on the flow behavior and are
often used in conjunction with other boundary conditions.

4. Compatibility with Laplace's Equation: Boundary conditions must be compatible with


Laplace's equation to ensure a physically meaningful solution. This typically involves

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.

By imposing appropriate boundary conditions, we can define the behavior of irrotational


motion within a given fluid domain and solve for the velocity potential function and flow
field.

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.

As we move forward, it is essential to remain vigilant of the assumptions and limitations of


potential flow theory while leveraging its benefits for practical applications. By integrating
theoretical insights with experimental validation and computational modeling, we can
continue to push the boundaries of fluid dynamics research and contribute to advancements in
various fields of engineering, science, and technology.

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.

3. Introduction to Fluid Mechanics and Fluid Machines by Som and Biswas,2 nd


edition,Tata- McGrawhill.

4. Batchelor, G. K. (2000). An Introduction to Fluid Dynamics. Cambridge


University Press.

5. Kundu, P. K., & Cohen, I. M. (2008). Fluid Mechanics. Academic Press.

6. Lamb, H. (1993). Hydrodynamics (6th ed.). Cambridge University Press.

7. White, F. M. (2011). Fluid Mechanics (7th ed.). McGraw-Hill Education.

8. Anderson, J. D. (2010). Fundamentals of Aerodynamics. McGraw-Hill


Education.

9. Milne-Thomson, L. M. (2012). Theoretical Hydrodynamics. Dover


Publications.

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.

13. Moore, D. W. (2019). Introduction to the Practice of Statistics. W. H.


Freeman.

95

You might also like