FM Unit3
FM Unit3
FM Unit3
Prof. C. M. Sewatkar
Department of Mechanical Engineering,
Govt. College of Engineering and Research ,
Avasari (Kh) Tq: Ambegaon, Dist: Pune
Fluid Kinematics
Method of Study
Lagrangian Eulerian
Z
C. M. Sewatkar, Faculty, Mech. Engg. Dept., GCOEARA
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., GCOEARA
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., GCOEARA
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
( u ) ( v) ( w)
For steady flow + + =0
x y z
u v w
For incompressible fluid + + =0
x y z
( V ) = 0
( u ) ( v)
+ + =0 Generalized 2D continuity equation
t x y
( u ) ( v)
+ =0 2D continuity equation for steady flow
x y
(u ) (v)
+ =0 2D continuity equation for steady flow of
x y an incompressible fluid
( u )
+ =0 Generalized 1D continuity equation
t x
( u)
=0 1D continuity equation for steady flow
x
du
=0 1D continuity equation for steady flow of
dx an incompressible fluid
Stream tube
N’
M’
Let,
A = area at central plane
V = velocity at mid plane
ρ = density of fluid
N
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
C. M. Sewatkar, Faculty, Mech. Engg. Dept., GCOEARA
Continuity Equation for 1D flow
Further,
➢However, the above equation can also be applied to flow passages of large
areas, even if velocity is not uniform over a particular section i. e. it varies
from point to point over a section
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
C. M. Sewatkar, Faculty, Mech. Engg. Dept., GCOEARA
Tutorials
a) u = 2 x 2 ; v = xyz; w = ?
b) u = 2 x 2 + 2 xy; w = z 3 − 4 xz − 2 yz; v = ?
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
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
C. M. Sewatkar, Faculty, Mech. Engg. Dept., GCOEARA
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
C. M. Sewatkar, Faculty, Mech. Engg. Dept., GCOEARA
Tangential and Normal Accelerations
➢ Like velocity acceleration is also a vector
s O
Streamline
as dθ
V+δV
a
δV δVn
V
an
B
δVs
n
V
A
Vs = f1 ( s, n, t ) and Vn = f 2 ( s, n, t )
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
Vn
steady
Local normal acceleration flow
t
Vs Convective tangential acceleration
Vs
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
A
P(x,yz) u δx
X
C. M. Sewatkar, Faculty, Mech. Engg. Dept., GCOEARA
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
C. M. Sewatkar, Faculty, Mech. Engg. Dept., GCOEARA
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
C. M. Sewatkar, Faculty, Mech. Engg. Dept., GCOEARA
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
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
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
C. M. Sewatkar, Faculty, Mech. Engg. Dept., GCOEARA
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
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
= 2z
Some facts about Vorticity
= V = curl V
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
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
= 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
= 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
The product of slopes of two lines is -1, which means streamline and
equipotential line intersect orthogonally
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 C. M. Sewatkar, Faculty, Mech. Engg. Dept., GCOEARA
Assignment – Prepare a
descriptive note on the
methods of drawing the flow
net and attach it in your lab
practice journal.
u = um 1 − (r / R)
n