Scribd
Upload
112 views

Advanced Fluid Dynamics - Assignment

This document contains an assignment for an advanced fluid dynamics course. It covers several problems involving vector and tensor calculus concepts such as components of velocity vectors, cross products, determinants, and vorticity tensors. Key steps and equations are shown for calculating scalar and cross products of vectors, determinants of matrices, properties under transformations, and the vorticity tensor formula in index notation.

Uploaded by

Aravind
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
112 views

Advanced Fluid Dynamics - Assignment

This document contains an assignment for an advanced fluid dynamics course. It covers several problems involving vector and tensor calculus concepts such as components of velocity vectors, cross products, determinants, and vorticity tensors. Key steps and equations are shown for calculating scalar and cross products of vectors, determinants of matrices, properties under transformations, and the vorticity tensor formula in index notation.

Uploaded by

Aravind
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 10

ME 6303 - Advanced Fluid

Dynamics
Assignment 1

NATIONAL UNIVERSITY OF SINGAPORE


February 7, 2015
Authored by: Aravind Baskar A0136344

Problem 1.1
Components of u,v,w are ui,vi,wi respectively.
U = ui.ei = u1e1 +u2.e2+ u3e3
V = vi.ei = v1e1 +v2e2+ v3e3
U.V = (u1e1 +u2e2+ u3e3).( v1e1 +v2.e2+ v3.e3) = u1v1 +u2v2+ u3v3 = ui.vi ei.ei = 1; ei.ej = 0

(a)

U x V = (u1e1 +u2e2+ u3e3)x( v1e1 +v2e2+ v3e3) = e3u1v2 e2u1v3 - e3u2v1 + e1u2v3 + e2u3v1 e1u3v2
ei x ei = 1; ei x ej = 0;
U x V = eiujvk as =+1 if even permutation of i,j,k;-1 if odd permutation of i,j,k,0 otherwise ( b )
(U x V).W = (e3u1v2 e2u1v3 - e3u2v1 + e1u2v3 + e2u3v1 e1u3v2).( w1e1 +w2e2+ w3e3) = uivjwk ( c )
U.(V x W) = (e3v1w2 e2v1w3 - e3v2w1 + e1v2w3 + e2v3w1 e1v3w2).( u1e1 +u2e2+ u3e3) = uivjwk ( d )
(U x V) x W = (e3u1v2 e2u1v3 - e3u2v1 + e1u2v3 + e2u3v1 e1u3v2) x ( w1e1 +w2e2+ w3e3) = e2u1v2w1 +
e1u1v2w2 e3u1v3w1 + e1u1v3w3 . = (U.W)V (V.W)U

(e)

(U x V)2 = (U x V). (U x V) = (e3u1v2 e2u1v3 - e3u2v1 + e1u2v3 + e2u3v1 e1u3v2 ).(e3u1v2 e2u1v3 - e3u2v1 +
e1u2v3 + e2u3v1 e1u3v2) = (u1v2)2+ (u1v3)2 + (u2v1)2 + (u2v3)2+ (u3ev1)2+ (u3v2)2 - (u1v2 u2v1)
(u1v3 u3v1) - (u2v1 u1v2) - (u2v3 u3v2) - (u3v1 u1v3) - (u3v2 u2v3) = (UV)2 (U.V)2

Problem 1.2

=
[

11
21
31

det [A] = ||

12
22
32

11
21
31

13
23

]
33

12
22
32

13

||
23

33

(f)

=A11(A22A33-A32A23) A12(A21A33-A31A23) + A13(A21A32-A31A22)


= (A1iA2jA3k)

(g)

Det[Almn] =det[A] if l,m,n even permutation +1, -1, odd permutation, 0 otherwise.
det[A] = Ali Amj Ank

(h)

6
.
=

(i)

det[A] = Ali Amj Ank


Multiplying by on both sides we get det[A] = Ali Amj Ank
Det[A] = Ali Amj Ank/6

Problem 1.3
= + +








If i=l then,
= + +








= ( ) + ( ) + (
)






=

=
[

11
21

12
22

13
23

32

]
33

11

12

13

=
21
[
31

22

23

32

]
33

31

(j)


=
[

11

21

31

12

11

12

13

21
] [
33
31

22

23

32

]
33

22
32

13
23


11 11


12 21


13 31

=
12 21
[
12 31


12 22

.. Remaining terms

12 23
]
12 33


12 32

=
[

11

12

13

21

22

31
32

]
33

23

D = [][]

=
[

11

12

13

21

11

12

13

21
]
[

33
31

22

23

32

]
33

22
23

31
32

Problem 1.4

=
[

11
21
31

12
22
32

11

det [A] = ||
21

31

13
23

]
33

12
22
32

13

||
23

33

=A11(A22A33-A32A23) A12(A21A33-A31A23) + A13(A21A32-A31A22)


Aij = ei.ej
ei x ej = ek

(ei x ej ).ek = 1

[A][AT]

=
[

11
21
31

12
22
32

13
23

11
12

] [
33
13

21
22
23

31
32

]
33

=[I]
11
1
[A-1] = || [

]
33

C11 = A22A33-A32A23
11
1
[A-1] = || [

]
33

Problem 1.5

=
1
[
1

]
3

det [U] = u1(v2w3-w3v2) u2(v1w3-w1v3) + u3(v1w2-v1w2) = (uivjwk)

=
[

11
21
31

12
22
32

13
23

]
33

If W is symmetric, W12 = W21, W31 = W13, W32 = W23


ui = Wjk
= W12 - W13 - W21 + W23 + W31 W32
=0
If W is antisymmetric, W12 = - W21, W31 = - W13, W32 = - W23

ui = Wjk
= W12 (-W13) (-W21)+ W23 + W31 (-W32)
= 2(W12 + W23 + W31)
1

Vorticity Tensor = ( )( )
2

)
= (

3
( )

1
1

)
( )
3

( )
1

)
3
( )

= (
1
( )
3

)
1
Vorticity Tensor = (

2
1
( )

3
1

1
2
1
3

)(
)(

1
3
1

(0)

1 | 1

=- ((

2
1

)(

= (

1
3

2
1

1
2

)(

)(

)(

3
1

)(

3
2

2
3

)(

)= xu

Problem 1.6
S = Sij.eiej D = Dlm.elem W = Wpq.epeq
D:S = (Dlm.elem ):( Sij.eiej)

2
1

)(

2
3

1
2

) (

)(

)(

3
2

)(

1
2

(0)

) (

(0)

Determinant = |( ) ( )
2 2
1
(

) (
3
)
( )

(
2

(0)
=

3
2

)(

3
2
3
2

)
)

3
2
3
2

3
3

)(
)(

2
3
2
3

(0)

)(
)(

(
3

2
3
2
3

)
)
)

)
)||

(0)

))+( (

3
2

)(

2
3

)(

1
2

)(

2
1

)(

2
3

)(

3
2

))

= Dlm.el

= DlmSij

Sij.ej

=DlmSml
D:ST = (Dml.emel ):( Sji.ejei)
= Dml.em Sij.ei

= DmlSij

=DlmSml

as Dml=Dlm

D : (1/2) (S+ST) =
(Dml.emel ):(1/2)( Sij.eiej + Sji.ejei)
= (1/2) (DlmSml DlmSml)
= DlmSml
W:S = (Wlm.elem ):( Sij.eiej)
= Wlm.el

= WlmSij

Sij.ej

=WlmSml
W:ST = (Wml.emel ):( Sji.ejei)
= Wml.em Sij.ei

= WmlSij

=-WlmSml
W : (1/2) (W-WT) = 0
D:W = 0

as Wml=-Wlm

Problem 1.10
(
1

1
1

Strain Rate Tensor = (


2
1
( )
3

(
=

1
1
1
2
1
3

)+(
)+(
)+(

)
1
Vorticity Tensor = (

2
3
( )
1

2
1
3
1

)(
)(

1
1
2
1
3
1

3
3

1
3

1
1

2
1
2
2
2
3

)+(
)+(
)+(

1
2
2
2
3
2

1
1

)(
1
( )
( )

1
2
1
3

1
2

)
)

)(

2
1

0
(

3
2

)(

+ ) . ( + + ) =3

(
(

3
1
3
2
3
3

1
3

)
)

3
3

)+(
)+(
)+(

1
3
2
3
3
3

)
)
)

( )
1

)
3
( )
3

(
(

Problem 1.11

. = (

)+(
3
( )

0
=

3
2
3

)(

) (
0

1
3
2

)
)

= (

(a. ) = (

= (

+ ) ( + + ) =0

) . ( + + ) =3a

+ + ) . () = ( + + ) -- to prove

Problem 1.12
S = Srrerer + Srere + Srzerez + See + Sreer + Szeez + Szzezez + Szeze + Szrezer

() =
(
) others are 0
( ) =

. S = (

() + (

+ + ) () = (

) (

) + () + (

) (

) + (

) + ( ) + (

) + (

) + (

) + (

) (

) + (

) + (

)+

)+

Problem 1.13
2u = . (u)
= + +
u = (

+ + ) ( + + )

2u = (
1

+
1

( ) +

+ + ) . ( +

( ) +

+ ( ) +

+ ( ) +

=(

2
2

2
2

2
2

2
2

2
2

+ ( )

( ) +

Problem 1.14
u. (u) = ( + + ).(
1

=(

( )

( ) +

( ) + ( )

+ ( ) +

)
1

( ) ( )

You might also like