1 Description

Velocity Field U = Ux i + Uy j + Uz z = Ur r + Uz z + U = Ur r + U + U
First Viscosity
Second Viscosity

2 Continuity Equation
2.1 General Form

+ (U) = 0
t

2.2 Cartesian Coordinate

(Ux ) (Uy ) (Uz )
+ + + =0
t x y z

2.3 Cylindrical Coordinate

1 (rUr ) 1 (U ) (Uz )
+ + + =0
t r r r z

2.4 Spherical Coordinate

1 (r2 Ur ) 1 (U sin ) 1 (U )
+ 2 + + =0
t r r r sin r sin

2.5 Orthogonal Curvilinear Coordinate

" #
U1 U2 U3 (h2 h3 U1 ) (h1 h3 U2 ) (h1 h2 U3 )
+ + + + + + =0
t h1 x1 h2 x2 h3 x3 h1 h2 h3 x1 x2 x3

1
3 Equation of Motion
3.1 General Form
3.1.1 General From
   
DU U U UU
= + (U )U = + ( ) + U ( U)
Dt t t 2
= f + T = f p + T0

3.1.2 Constant Viscosity

   
DU U U UU
= + (U )U = + ( ) + U ( U)
Dt t t 2
= f p ( U) + ( + 2)( U)
= f p 2 U + ( + )( U)

3.2 Cartesian Coordinate

3.2.1 Variable Viscosity
     
DUx p Ux Ux Uy Ux Uz
= fx + 2 + U + ( + ) + ( + )
Dt x x x y y x z z x
     
DUy p Uy Uy Uz Ux Uy
= fy + 2 + U + ( + ) + ( + )
Dt y y y z z y x y x
     
DUz p Uz Uz Ux Uy Uz
= fz + 2 + U + ( + ) + ( + )
Dt z z z x x z y z y

3.2.2 Constant Viscosity

DUx p
= fx + 2 Ux + ( + ) U
Dt x x
DUy p
= fy + 2 Uy + ( + ) U
Dt y y
DUz p
= fz + 2 Uz + ( + ) U
Dt z z

2
3.3 Cylindrical Coordinate
3.3.1 Variable Viscosity

r-direction
! !
DUr U2 Ur Ur Ur U Ur U2
= + Uz + Ur +
Dt r t z r r r
" # " #
p Ur Uz Ur
=fr + ( + ) + 2 + U
r z z r r r
" # " #
1 1 Ur U U 2 Ur 1 U Ur
+ ( + ) +
r r r r r r r r

z-direction
! !
DUz Uz Uz Uz U Uz
= + Uz + Ur +
Dt t z r r
" # " #
p Uz 1 Ur Uz
=fz + 2 + V + r( )
z z z r r z r
" #
1 1 Uz U
+ ( + )
r r z

-direction
! !
DU Ur U U U U U U Ur U
+ = + Uz + Ur + +
Dt r t z r r r
" #
1 p 1 Uz U
=f + ( + )
r z r z
" #
1 Ur U U
+ ( + )
r r r r
" # " #
1 2 U 2 1 Ur U U
+ + U + +
r r r r r r

3.3.2 Constant Viscosity

r-direction
! !
DUr U2 Ur Ur Ur U Ur U2
= + Uz + Ur +
Dt r t z r r r
" #
p Ur 2 U
=fr + 2 Ur 2 2 + ( + ) ( U)
r r r r

z-direction
! !
DUz Uz Uz Uz U Uz
= + Uz + Ur +
Dt t z r r
p
=fz + 2 Uz + ( + ) ( U)
z z
-direction
! !
DU Ur U U U U U U Ur U
+ = + Uz + Ur + +
Dt r t z r r r
" #
1 p U 2 Ur 1
=f + 2 U 2 + 2 + ( + ) ( U)
r r r r

3
3.4 Spherical Coordinate
3.4.1 Variable Viscosity

r-direction
! " #
DUr U2 + U2 p Ur
=fr + 2 + U
Dt r r r r
" #
1 U Ur
+ r ( ) +
r r r r
" #
1 Ur U
+ + r ( )
r sin r sin r r
" #
Ur 2 U 4Ur 2 U 2U cot U cot Ur
+ 4 + r cot ( ) +
r r r r r sin r r r r

-direction
! " ! #
DU Ur U U2 cot 1 p 1 1 U Ur
+ =f + 2 + + U
Dt r r r r r r
" #
1 sin U U
+ ( )+
r sin r sin r sin
" #
U Ur
+ r ( ) +
r r r r
" #
2 cot U 2 cot U 2U cot2 U 3 Ur
+ + 3r ( ) +
r r r sin r r r r

-direction
! " #
DU U Ur U U cot 1 p 1 1 U Ur U cot
+ + =f + 2( + + ) + U
Dt r r r sin r sin r sin r r
" # " #
Ur U 1 sin U U
+ + r ( ) + ( )+
r r sin r r r r sin r sin
" #
3 Ur U 2 cot sin U 2 cot U
+ + 3r ( ) + ( )+
r r sin r r r sin r sin

3.4.2 Constant Viscosity

r-direction
! !
DUr U2 + U2 p 2Ur 2 U 2U cot 2 U
=fr + 2 Ur 2 2 2 + ( + ) ( U)
Dt r r r r r2 r sin r

-direction
! !
DU Ur U U2 cot 1 p 2 Ur U 2 cos U 1
+ =f + 2 U + 2 2 2 2 2 + ( + ) ( U)
Dt r r r r r sin r sin r

-direction
! !
DU U Ur U U cot 1 p U 2 Ur 2 cos U 1
+ + =f + 2 U 2 2 + 2 + 2 2 + ( + ) ( U)
Dt r r r sin r sin r sin r sin r sin

4
3.5 Orthogonal Curvilinear Coordinate
x1 -direction
" ! !#
U1 U1 U1 U2 U1 U3 U1 U2 h2 U1 h1 U1 h1 U3 h3
+ + + U2 + U3
t h1 x1 h2 x2 h3 x3 h2 h1 x1 h1 h2 x2 h1 h3 x3 h3 h1 x1
"   
1 p 1 1 1 U1 U2 h1 U3 h1
= f1 + ( U) + 2h2 h3 + +
h1 x1 h1 x1 h1 h2 h3 x1 h1 x1 h1 h2 x2 h3 h1 x3
     #
h2 U2 h1 U1 h1 U1 h3 U3
+ h3 h1 ( )+ ( ) + h1 h2 ( )+ ( )
x2 h1 x1 h2 h2 x2 h1 x3 h3 x3 h1 h1 x1 h3
   
h2 U2 h1 U1 h1 h3 U3 h1 U1 h1
+ ( )+ ( ) + ( )+ ( )
h1 h2 h1 x1 h2 h2 x2 h1 x2 h1 h3 h1 x1 h3 h3 x3 h1 x3
   
2 1 U2 U3 h2 U1 h2 h2 2 1 U3 U1 h1 U3 h3 h3
+ + + +
h1 h2 h2 x2 h2 h3 x3 h1 h2 x1 x1 h1 h3 h3 x3 h1 h3 x1 h2 h3 x2 x1

x2 -direction
" ! !#
U2 U1 U2 U2 U2 U3 U2 U3 h3 U2 h2 U2 h2 U1 h1
+ + + U3 + U1
t h1 x1 h2 x2 h3 x3 h3 h2 x2 h2 h3 x3 h2 h1 x1 h1 h2 x2
"   
1 p 1 1 h2 U2 h1 U1
= f2 + ( U) + h2 h3 ( )+ ( )
h2 x2 h2 x2 h1 h2 h3 x1 h1 x1 h2 h2 x2 h1
     #
1 U2 U3 h2 U1 h2 h3 U3 h2 U2
+ 2h3 h1 + ) + + h1 h2 ( )+ ( )
x2 h2 x2 h2 h3 x3 h1 h2 x1 x3 h2 x2 h3 h3 x3 h2
   
h3 U3 h2 U2 h2 h2 U2 h1 U1 h2
+ ( )+ ( ) + ( )+ ( )
h2 h3 h2 x2 h3 h3 x3 h2 x3 h2 h1 h1 x1 h2 h2 x2 h1 x1
   
2 1 U3 U1 h3 U2 h3 h3 2 1 U1 U2 h1 U3 h1 h1
+ + + +
h2 h3 h3 x3 h3 h1 x1 h2 h3 x2 x2 h2 h1 h1 x1 h2 h1 x2 h3 h1 x3 x2

x3 -direction
" ! !#
U3 U1 U3 U2 U3 U3 U3 U1 h1 U3 h3 U3 h2 U2 h2
+ + + U1 + V2
t h1 x1 h2 x2 h3 x3 h1 h3 x3 h3 h1 x1 h3 h2 x1 h2 h3 x3
"   
1 p 1 1 h1 U1 h3 U3
= f3 + ( U) + h2 h3 ( )+ ( )
h3 x3 h3 x3 h1 h2 h3 x1 h3 x3 h1 h1 x1 h3
     #
h3 U3 h2 U2 1 U3 U1 h3 U2 h3
+ h3 h1 ( )+ ( ) + 2h1 h2 + +
x2 h2 x2 h3 h3 x3 h2 x3 h3 x3 h3 h1 x1 h2 h3 x2
   
h1 U1 h3 U3 h3 h3 U3 h2 U2 h3
+ ( )+ ( ) + ( )+ ( )
h1 h3 h3 x3 h1 h1 x1 h3 x1 h2 h3 h2 x2 h3 h3 x3 h2 x2
   
2 1 U1 U2 h1 U3 h1 h1 2 1 U2 U3 h2 U1 h2 h2
+ + + +
h1 h3 h1 x1 h1 h2 x2 h3 h1 x3 x3 h2 h3 h2 x2 h2 h3 x3 h1 h2 x1 x3

5
4 Stress Tensor in various coordinate
4.1 General Form
 
Uk Ui Uj
ij = pij + ij + + (1)
xk xj xi

4.2 Cartensian Coordinate

Ux Uy Uz
xx = p + 2 + U ; yy = p + 2 + U ; zz = p + 2 + U
x y z

     
Ux Uy Uy Uz Uz Ux
xy = yx = + ; yz = zy = + ; xz = zx = +
y x z y x z

4.3 Cylindrical Coordinate

 
Uz Ur 1 U Ur
zz = p + 2 + U ; rr = p + 2 + U ; = p + 2 + + U
z r r r

     
Uz Ur 1 Ur U U U 1 Uz
zr = rz = + ; r = r = + ; z = z = +
r z r r r z r

4.4 Spherical Coordinate

   
Ur 1 U Ur 1 U Ur U cot
rr = p + 2 + U ; = p + 2 + + U ; = p + 2 + + + U
r r r r sin r r

     
 U  1 Ur sin  U  1 U 1 Ur  U 
r = r = r + ; = = + ; r = r = +r
r r r r sin r sin r sin r r

4.5 Curvilinear Coordinate

 
1 U1 U2 h1 U3 h1
11 = p + 2 + + + U
h1 x1 h1 h2 x2 h3 h1 x3
 
1 U2 U3 h2 U1 h2
22 = p + 2 + + + U
h2 x2 h2 h3 x3 h1 h2 x1
 
1 U3 U1 h3 U2 h3
33 = p + 2 + + + U
h3 x3 h3 h1 x1 h2 h3 x2

 
h2 U2
h1 U1

12 = 21 = +
h1 x1 h2 h2 x2 h1
 
h3 U3
h2 U2

23 = 32 = +
h2 x2 h3 h3 x3 h2
 
h1 U1
h3 U3

31 = 13 = +
h3 x3 h1 h1 x1 h3