FM Unit3
FM Unit3
FM Unit3
Prof. C. M. Sewatkar
Method of Study
Lagrangian Eulerian
Z
C. M. Sewatkar, Faculty, Mech. Engg. Dept., COE, Pune
Velocity of Fluid Particles
V f ( x, y , z , t )
u f ( x, y , z , t )
v f ( x, y , z , t )
w f ( x, y, z , t ),
Vector notation V iu jv kw
dV du dv dw dp d
Mathematically, =0, =0, =0, =0, =0, =0,
dt dt dt dt dt dt
Unsteady flow - – if at any point in the flow system of a fluid the various
characteristics such as velocity, pressure, density, temperature etc.
which describe the behavior of the fluid flow, change with time then the
flow is called steady flow
dV du dv dw dp d
Mathematically, 0, 0, 0, 0, 0, 0,
dt dt dt dt dt dt
dV du dv dw
Mathematically, =0, =0, =0, =0
ds ds ds ds
dV du dv dw
Mathematically, 0, 0, 0, 0
ds ds ds ds
Q = Constant
Steady uniform flow
Q Constant
Unsteady uniform flow
Q = Constant
Q Constant
Unsteady non-uniform flow
Steady flow V f ( x, y , z ) V f ( x, y ) V f ( x)
unsteady flow V f ( x, y , z , t ) V f ( x, y , t ) V f ( x, t )
Laminar flow – The flow is said to be laminar when the various fluid
particles move in layers (or laminae) with one layer of fluid sliding
smoothly over an adjacent layer
v V
θ
u
P(x,y)
V
0 X
Streamline is an imaginary curve drawn through a flowing fluid
such that a tangent to it at any point gives the direction of the
velocity of flow at that point.
C. M. Sewatkar, Faculty, Mech. Engg. Dept., COE, Pune
Streamlines
Mathematical representation of streamline
v dy
tan dy and dx are the y and x components
of the differential displacement along
u dx the streamline in the vicinity of P
dx dy
; or (udy vdx) 0
u v
dx dy dz
u v w
C. M. Sewatkar, Faculty, Mech. Engg. Dept., COE, Pune
Characteristics Streamline
Path line is the line traced by a single fluid particle as it moves over
a period of time
Streak line is the line traced by the fluid particles passing through a
fixed point in the flow system
Wake
Fluid Mechanics
Heat Transfer
Fixed region
Rate of increase or
decrease of fluid mass Rate of mass flow Rate of mass flow
at the entrance at the exit
within fixed region
C C’
ρ = density of fluid
B B’
P(x,y,z)
D u
δy D’
δz
A
δx A’
( u ) ( v) ( w)
x y z Equation A
x y z
( x y z )
Rate of increase of mass of fluid in the parallelepiped with time
( x y z ) ( x y z ) Equation B
t t
Equate Equation A and B
( u ) ( v) ( w)
( x y z ) x y z
t x y z
The most general form
( u ) ( v) ( w)
0 of continuity equation
t x y z applicable to all types
of flows
( V ) 0
( u ) ( v)
0 Generalized 2D continuity equation
t x y
( u ) ( v)
0 2D continuity equation for steady flow
x y
( u )
0 Generalized 1D continuity equation
t x
( u)
0 1D continuity equation for steady flow
x
Stream tube
N’
M’
Let,
A = area at central plane
V = velocity at mid plane
ρ = density of fluid
δs
s
AV ( AV )
s 2
Mass of fluid entering through plane N’M’ per unit time
s
AV ( AV )
s 2
Net mass of fluid that has remained in the fluid element per unit time
( AV ) s
s
Mass of fluid element A s
Rate of increase of mass of fluid element ( A) s
t
Continuity Equation for 1D flow
Thus
( AV ) s ( A) s
s t
( A) ( AV ) 0
t s
For steady flow ( AV ) 0 AV Constant
s
1 AV
1 1 2 A2V2 3 A3V2
Continuity Equation for 1D flow
Further,
M
dA1, dA2, dA3
v1, v2 and v3
Q vdA
a) u Cx ; v Cy
b) u 3x y; v 2 x 3 y
c) u x y; v x 2 y
d) u A sin xy; v A sin xy
e) u 2 x 2 3 y 2 ; v 3xy
Tutorials
a) u 2 x 2 ; v xyz; w ?
b) u 2 x 2 2 xy; w z 3 4 xz 2 yz; v ?
Tutorials
V (5 x)i (5 y ) j (10 z )k
du dv dw
ax lim ; a y lim ; az lim ;
dt 0 dt dt 0 dt dt 0 dt
We known u f ( x, y , z , t )
Total derivative of u using partial derivative is
du u dx u dy u dz u dt
dt x dt y dt z dt t dt
dx dy dz
lim u , lim v , lim w
dt 0 dt dt 0 dt dt 0 dt
Acceleration of Fluid Particle
Thus du u u u u
lim u v w
dt 0 dt x y z t
u u u u
ax u v w
x y z t
Similarly v v v v
ay u v w
x y z t
w w w w
az u v w
x y z t
In vector notation the acceleration of a fluid particle is written as:
a V .V
Acceleration of Fluid Particle
u v w
t t t
u u u
u ;v ;w ;
x y z
Increase in velocity due to change in
v v v
u ;v ;w ; position of particle
x y z
w w w Convective of spatial acceleration
u ;v ;w ;
x y z
Tangential and Normal Accelerations
Like velocity acceleration is also a vector
s O
Streamline
as dθ
V+δV
a
δV δVn
V
an
B
δs δVs
n
V
A
dVs dVn
as lim and an lim
dt 0 dt dt 0 dt
Vs Vs
as =Vs +
s t
Vn Vn
an Vs +
s t
Vn
Note that though Vn = 0 is not equal to zero
s
Vn is zero at any point on the streamline but at any
other point on the streamline the component of the
velocity in the direction parallel to the of Vn need not be
always zero
Tangential and Normal Accelerations
Further
s Vn
d
r V
Vn Vn V Vs
or (since V Vs )
s s r r
Vs2 Vn
an +
r t
Tangential and Normal Accelerations
Vs
Local tangential acceleration Zero
t for
steady
Vn Local normal acceleration flow
t
Vs
Vs Convective tangential acceleration
s
Vn Vs2 Convective normal acceleration
Vs
s r
Vs
as =Vs
s
For steady flow
Vs2
an
r
Tangential and Normal Accelerations
Y
u
u y
y
B’
B
δθ2 A’
δy
v
v v x
x
δθ1
P(x,yz)
A
u δx
X
Rotational and Irrotational Flow
We can write
1
PA lim θ = s/r
t 0 t
v
v x x v t
PA lim
t 0 x t
v
Thus, PA
x
Rotational and Irrotational Flow
Also
2
PB lim θ = s/r
t 0 t
u
u y u t
y
PB lim
t 0 y t
u
Thus, PB
y
Rotational and Irrotational Flow
The component of rotation – the average angular velocity of two
infinitesimally linear elements in the particle that are perpendicular to
each other and to the axis of rotation (In this case Z axis)
1
z (PA PB )
2
1 v u 1 w v 1 u w
z x y
2 x y 2 y z 2 z x
Rotational and Irrotational Flow
If the component of rotation for all the axes is zero the flow is said
to be irrotational else it is rotational
Thus, for flow to be irrotational
x 0 w v
y z
y 0 u w
z x
z 0 v u
x y
V
ds V cos ds
α C
Vcosα
(udx vdy wdz )
C
Circulation and Vorticity
Circulation around an elementary rectangle
u y
u
y 2
D C
v
v x v x
u v
v x 2
x 2
A B
u y
u
y 2
u y
Circulation along AB = u x
y 2
Circulation and Vorticity
v x
Circulation along BC = v y
x 2
u y
Circulation along CD = u x
y 2
v x
v
Circulation along DA = y
x 2
What ever may be the shape of the curve the circulation must be
equal to the sum of the circulation around the elementary surfaces
of which it consists, provided the boundary of the curve is wholly
in the fluid
Circulation and Vorticity
Thus,
AB BC CD DA
u y v x u y v x
u x v y u x v y
y 2 x 2 y 2 x 2
v u
x y
x y
Circulation
Area
v u
x y
2 z
Some facts about Vorticity
V curl V
Simple Bluff Body Flow
Problem
Velocity Potential
Mathematically, if Φ = f(x, y, z, t)
u ; v and w
x y w
1 w v 1 2 2
x x
2 y z 2 yz zy
1 u w 1 2 2
y y
2 z x 2 zx xz
1 v u 1 2 2
z z
2 x y
2 xy yx
Velocity Potential
If Φ is a continuous function
2 2 2 2 2 2
yz zy zx xz xy yx
B
(x+δx, y+δy)
C
v Flow across the curve ACB in x-direction
D
u y
(x, y) A u
Flow across the curve ACB in y-direction
v x
Stream Function
If dΨ is assumed to represent the total flow across ACB
d u y v x Equation A
d x y Equation B
x y
Stream Function
Compare A and B
v and u
x y
Thus, the assumption that dΨ is the flow across two points is valid
and stream function can be used for determination of flow between
two points if the stream functions at these points is known.
Compare the above equations with the equations for velocity potential
u and v
x y y x
Stream Function
and - Cauchy-Riemann Equations
x y y x
1 v u
Further, we know that z
2 x y
1 1 2 2
z 2 2
Poisson’s
2 x x y y 2 x y equation
0
x y y x
2 2
xy yx
Thus if Ψ is the continuous function and its second derivative exists it
can be a possible case of flow since it satisfies the continuity
equation
Streamlines, Equipotential Lines and Flow Net
Equipotential Φ = C6
lines Φ = C5
Φ = C4 Ψ = C1
Φ = C3 Ψ = C2 Streamlines
Φ = C2 Ψ = C3
Φ = C1 Ψ = C4
Ψ = C5
Vs Ψ = C6
Vn
Streamlines, Equipotential Lines and Flow Net
We know
Vn Vn 0
n Φ is constant along n-direction
Vs The flow is along streamline as Φ is
s not constant along s-direction
Similarly
Vn 0 No flow in the direction
s normal to streamline but the
flow is always along
streamline
Vs
n
Assignment – Prepare a
descriptive note on the
methods of drawing the flow
net and attach it in your lab
practice journal.
Tutorials
Two large circular plates are kept at a distance y0 apart
and contain an incompressible fluid in between. If the
bottom plate is fixed and top plate is moved downward at
a constant velocity of V0, estimate the velocity at which
the fluid moves at a radial distance of r.
Tutorials
u um 1 (r / R )
n