TPI 17 CHP 3 Revised-BTE
TPI 17 CHP 3 Revised-BTE
TPI 17 CHP 3 Revised-BTE
1
Scope of the Chapter
(u z ) z z
y y
(u x ) x (u x ) x x
Flowing fluid
y
z z
z (u z ) z
y
z
z x
control x x x
volume
x
Mass balance: rate of rate of rate of
mass - mass increase (kg/s)
in out of mass 3
Mass balance over a volume element xyz:
A x (u x ) x - (u x ) x x A y (u y ) - (u y )
y y y
A z (u z ) z - (u z ) z z V
t
A x yz, A y xz, A z xy, V xyz
yz (u x ) x - (u x ) x x xz (u y ) - (u y )
y
y y
xy (u z ) z - (u z ) z z xyz
t
By dividing equation throughly with xyz we get:
x y z 4 t
By taking limit as x, y &z go zero we get:
(.v)
t Net rate of mass addition
Rate of increase per unit volume by
of mass per unit convection
volume
(.v) 0
Notice that:
6
The equation of Motion
Fluid flow decribed by cartesian coordinate
y
y
x
Momentum transport
z in the fluid through
flowing fluid x the 3 directions
y
control volume z
z x
y
( yx )
y y
y
( zx ) z z
x
z
x
y
( xx ) x ( xx ) x x
( zx ) z
z y
z x
z ( yx )
y
z x
x
Notice that: Momentum enters across & leaves the surfaces
normal to x-dir. by all mechanisms (both convective & molecular
8
transport).
Momentum balance in the x-dir.:
A x ( xx ) x - ( xx x x
) A y ( yx ) - ( yx A z ( zx ) z - ( zx z z
y y y
(u x V)
(V)g x
t
A x yz, A y xz, A z xy, V xyz
yz( xx ) x - ( xx x x
) xz( yx ) - ( yx xy( zx ) z - ( zx z z
y y y
(u x xyz)
(xyz)g x
t
By dividing equation throughly with xyz we get:
x y z
(u x )
t
Taking limit as x, y &z go zero:
x-component: (u x ) xx yx zx
- g x
t x y z
Similarly:
(u y ) xy yy zy
y-component: - g y
t x y z
z-component: (u z ) xz yz zz
- g z 10
t x y z
By using vector-tensor notations, 3 equations can be written
as follows:
u i -.i g i (i x, y, z)
t
vv p
v -.vv - p - . g (3.2 9)
t
Rate of Rate of mom. Rate of mom. External force
increase of addition by addition by mol. on fluid per
mom. per conv. Per unit transport per unit vol.
unit vol. vol. unit vol.
Differential volume element in cylindrical coordinates
(r, , z)
Directions:
Radial : r
Axial : z
Angular :
x r cos ,
Coordinate transformation between
rectangular & cylindrical coordinates: y r sin,
zz
Differential volume element in spherical coordinates
(r, , )
Directions:
Radial : r
Angular : &
u y u y u y u y p τ xy τ yy τ zy
ux uy uz g y
t x y z y x y z
u z u z u z u z p τ xz τ yz τ zz
ux uy uz g z
t x y z z x y z
Newton’s Law of Viscosity
Cartesian Coordinates:
u x 2 u z u y
τ xx μ 2 μ κ u τ yz τ zy μ
x 3 y z
u y 2 u x u z
τ yy μ 2
μ κ u τ zx τ xz μ
y 3 z x
u z 2 u y u x
τ zz μ 2 μ κ u τ xy τ yx μ
z 3 x y
u x u y u z
in which: u
x y z
15
Newton’s Law of Viscosity
Cartesian Coordinates:
u x 2 u z u y
τ xx μ 2 μ κ u τ yz τ zy μ
x 3 y z
u y 2 u x u z
τ yy μ 2
u τ zx τ xz μ
x
μ κ
y 3 z
u z 2 u y u x
τ zz μ 2 μ κ u τ xy τ yx μ
z 3 x y
u x u y u z
in which: u
x y z
16
Navier-Stokes Eqns in Rectangular Coordinates
Newtonian Fluid with Constant and
u x u x u x u x p 2u x 2u x 2u x
ux uy uz 2 2 2 g x
t x y z x x y z
u y u y u y u y p 2u y 2u y 2u y
ux uy uz 2 2 2 g y
x
t x y z y y z
u z u z u z u z p 2u z 2u z 2u z
ux uy uz 2 2 2 g z
t x y z z x y z
u
ρ u u p μ 2u ρg
t 17
Equations of Change in Cylindrical Coordinates
1 (rur ) 1 ( u ) ( u z )
Continuity: 0
t r r r z
Motion:
ur ur u ur u2 u p 1 r rr 1 r rz
ur u z r g r
t r r r z r r r r 2
r z
u
u u u ur u u
1 p 1 r 2 r 1 z
ur uz 2 g
t r r r z r r r r z
u z u u u z u p 1 r rz 1 z zz
ur z uz z g z
t r r z z r r r z
18
Newton’s Law of Viscosity in Cylindrical Coordinates
ur 2 u z u r
τ rr μ 2 μ κ u τ rz τ zr μ
r 3 r z
u z 2 u 1 u z
τ zz μ 2 μ κ u τ z τz μ
z 3 z r
1 u ur 2
τ μ 2 2 μ κ u
r r 3
u 1 ur
τ r τr μ r
r
r r
1 (rur ) 1 u u z
in which: (.u )
r r r z 19
Navier-Stokes Eqns in Cylindrical Coordinates
Newtonian Fluid with Constant and
ur ur u ur u2 ur
ur uz
t r r r z
p 1 rur 1 2ur 2 u 2ur
2 2 2 g r
r r r r r r z
2
u u u u u u u
ur r u z
t r r r z
1 p 1 ru 1 2u 2 ur 2u
2 2 2 g
r r r r r
2
r z
u z u z u u z u z
ur uz
t r r z
p 1 u z 1 uz uz
2 2
r 2 2 g z
z r r r r z
2 20
Equations of Change in Spherical Coordinates
1 ( r 2ur ) 1 ( u sin ) 1 ( u )
Continuity: 2 0
t r r r sin r sin
ur u u u u u u 2
u 2
p
Motion: t ur r r r sin r r
r r r
2
1 r 2 τ rr
1 ( τ r sin )
1 τ r
g r
r r r sin r sin r
u u u u u u u u u 2
cot 1 p
ur
r
r r sin r
t r r r
2
1 r 2 τ r
1 ( τ sin )
1 τ r cot
g
r r r sin r sin r r
u u u u u u u ur u u cot 1 p
ur
t r r r sin r r r sin
2
1 r 2 τ r 1 ( τ )
1 τ τ r 2 cot
τ g
r r r sin
21
r r r
Newton’s Law of Viscosity in Spherical Coordinates
u 2 u 1 ur
τ rr μ 2 r μ κ u τ r τ r μ r
r 3 r r r
1 u ur 2
τ μ 2 μ κ u
r r 3
1 ur u
τr τ r μ r
r sin r r
1 u ur u cot 2
τ μ 2 μ κ u
r sin r 3
sin u 1 u
τ τ μ
r sin r sin 22
Navier-Stokes Eqns in Spherical Coordinates
Newtonian Fluid with Constant and
ur ur u u r ur
p
2 2
u u u
ur
t r r r sin r r
1 2 ur 1 u r
r
r 2 r r r 2 sin sin
g r
1 ur 2ur 2 u 2u cot
2
2 u
2 2 2 2 2
r sin r r sin
2 2
r r
u u u u u u u u u cot
2
1 p
ur
r
r r sin r
t r r r
1 2 u 1 u
r 2 r r r r 2 sin sin
g
1 u 2 ur
2
u 2 cos u
2 2 2 2 2 2 2
r sin r r sin r sin
2
23
Navier-Stokes Eqns in Spherical Coordinates
Newtonian Fluid with Constant and
2 1 2
r sin 2 r 2 sin 2 r 2 sin r 2 sin 2
24
The Equation of Mechanical Energy
Equation of mechanical energy (ME) may obtained from
equation of motion by some (lengthy) modification &
involves only mechanical terms
1 2
( ) -. v - (.pv) - p(-.v) - .(.v) - ( : v)
2
t 2
( v.g)
Reversible Work done (J/s)
Increase of Increase of Work done
conversion of by viscous
kinetic kinetic pressure of Work done
kinetic energy forces on
energy per energy per surroundings by external
into internal the fluid
unit vol. unit vol. on fluid forces on the
energy
fluid
Reversible
conversion
from kinetic
into internal 25
energy
The Equation of Mechanical Energy
(v.g) ( v. ) (. v) (. v) (. v) ( / t )
Velocity distribution
Pressure distribution
Density distribution
Temperature distribution
Mass distribution
Example: Operation of a Couette Viscometer
The measurement of pressure difference vs. mass flow
through a cylindrical tube is the basis for determination of
viscosity in commercial capillary viscometers. The viscosity
may be determined by measuring the torque required to turn
a solid object in contact with a fluid (Couette instrument, in
Fig3.6-1)
u 2 p
(p = p(r)) ( 4)
r r
-dir.:
() r 0 z 0
SS 0 0
0 0 0
u u u u u u u
ur r u z r 0
t r r r z 0 0 0 ignored
r 0
1 p 0 1 ru 1 2u 2 ur 2u
2 2 2 g
r r r r r
2
r z
p p() (5)
1 ru ( 6)
0 (Velocity profile
r r r equation)
z-dir.:
z 0 z 0 z 0 z 0
0 0 0
u z u z u u z u z
ur uz
t r r z 0 0
0
SS p 1 u z 1 2
u 2
uz
r 2 z
2 g z
z r r r r z
2
z 0 z 0 z 0 (7 )
p (p = p(z))
g z (8)
z
The simplified equations:
u 2 p
(p = p(r)) ( 4)
r r
1 ru ( 6)
0
(Velocity profile
r r r equation)
p (p = p(z))
g z (8)
z
Integration of Eq.6:
1 ru
C1'
r r
ru
C1' r
r
C1' 2
ru r C2
2
C1' C2
u r (C1 C1' )
2 r
BC1 : @r R, 0
BC2 : @r R, R0
C2 0
From BC1 : 0 C1 ( R) C1
R (1 2 )
C2 ( R ) 2 0
From BC2 : R0 C1 (R) C2
R (1 2 )
R r 1 (Velocity profile
u R 0 /( )
r R equation)
0
u 1 u r
r r r
r
r r
u
2
2
R
r r r 0
r r r R r 1
Tz k 0 (k 4R 2 L( 2 / 1 2 ))
Tz
k / 4R 2 L( 2 / 1 2 )
Slop = k
0
Pressure distribution in the viscometer:
p p( r , z )
p p
dp dr dz
r z
p R 0 R r ( R ) 2 r 2 1
2 2
1
A 2 2 2
r
r r r R r ( R )
p ( R ) 2
A 3
2
r 2
A R0 /(
1 2
2
r r r ( R )
p
g z
z
( R ) 2 2 r2
p(r, z) A 3 2
dr (g z )dz
r r ( R )
( R ) 2 r2
p(r, z) A 2 ln r 2
g z z C1
2r
2
2( R )
BC : @ r R, z z 0 ; p patm
C1 patm 2A ln(R ) g z z 0
Example: The shape of the surface of a rotating liquid
Assumptions:
Steady-state
Fluid properties are constant z0
Analysis:
Similar to the previous example, the equation of continuity &
equation of motion may be simplified to obtain:
SS ur 0 u u () uz 0
1 (rur ) 1 ( u ) ( u z )
0
t r r r z
SS ur 0 ur 0 ur 0
ur ur u ur u2 ur
ur uz
u u () u r 0
t r r r z
p 1 rur 1 2ur 2 u 2ur
2 2 2 g r
r r r r r r z
2
45
SS u r 0 u u () u r 0 u z 0
u u u u ur u u
ur uz u r 0 u u (z)
t r r r z
r 0
1 p 1 ru 1 2u 2 ur 2u
2 2 2 g
r r r r r
2
r z
p p()
SS ur 0 uz 0 uz 0
u z u z u u z u z
ur uz
t r r z
p 1 u z 1 uz uz
2 2
r 2 2 g z
z r r r r z
2
46
The simplified equations:
u 2 p
(1)
r r
1 ru ( 2)
0 (Velocity profile equation)
r r r
p (p = p(r, z))
g z (3)
z
C2
u C1r
( 4)
Integration of Eq. (2):
r
C2
u C1r
r
BC1: @r 0, u finite (5)
BC2: @r R, u R ( 6)
C2
From BC1: 0 0 C2 0
0
From BC2: R C1R C1
There is no shear stress, therefore, no deformation.
u r Every fluid particle moves with the same angular
(7 )
velocity.
This equation states that each element of the rotating fluid moves
as an element of rigid body.
Surface shape may be modelled through pressure distribution
in the fluid body.
Pressure distribution:
Method-1:
p u
2 r z0
r r (8)
p (9)
g z g
z
p p
dp dr dz ( 2 r)dr (g)dz (10)
r z
p p 1
r z
r gz C1
2 2
dp dr dz p(r, z) (11)
2
(12)
BC : @ r 0, z z 0 ; p patm
patm 0 gz0 C1 C1 patm gz0
1
p(r, z) 2 r 2 gz p atm gz 0
2
1 (13)
p(r, z) p atm 2 r 2 g(z z 0 )
2
This equation gives the pressure at all points within the liquid.
At air-liquid interface, p = patm
Pressure distribution at the interface (parabola)
2g
Method-2:
u r
p u 2
1 2 2
p(r,z) r f1 (z) C
(15)
r r 2
p
g p(r,z) f 2 (r) gz C (16)
z
These two equations must be identical, thus,
1
f 2 (r) 2 r 2 , f1 (z) gz
2
Therefore,
1 2 2
p(r,z) r gz C (17)
2
At the liquid-gas interface:- Viscous stresses and surface
tension are negligible. Thus,
pliquid pgas BC: @z h(r), p p0
1 2 2 1 2 2
p(r,z) r gz C p0 r gh(r) C (18)
2 2
2 2 C p0
h(r) r (19)
2g g
Recall that:
@r 0, h(r) z 0
To determine C: The total volume of the liquid is,
(20)
dV 2h (r )rdr
R R 2 2 C p0
V V0 2 h(r)rdr 2 r rdr (21)
0 0
2g g
2 R 2 gV0
C p0 (22)
4 R 2
Inserting C into Eq.19 gives ,
R r 1
2 2 2
h(r) h 0 (h 0 V0 / R 2 ) (23)
2g R 2
Initial liquid level before rotation