APPH 4200 Physics of Fluids - Columbia...

35
APPH 4200 Physics of Fluids Kinematics: Describing Fluid Flow (Ch. 3) 1. Quick Review 2. Lagrangian (material) & Eulerian (field) 3. Streamlines & Pathlines 4. Point Deformations: Stretching, Pinching, and Rotating 5. Principal Axes (Cauchy-Stokes Decomposition) 1

Transcript of APPH 4200 Physics of Fluids - Columbia...

Page 1: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

APPH 4200 Physics of Fluids

Kinematics: Describing Fluid Flow (Ch. 3)

1. Quick Review2. Lagrangian (material) & Eulerian (field)3. Streamlines & Pathlines4. Point Deformations: Stretching, Pinching, and Rotating5. Principal Axes (Cauchy-Stokes Decomposition)

1

Page 2: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Vectors: U = {ux, uy, uz}L = r × p

Not Vector: list = {ρ, η, T}

Pseudovector: L = r × p

2

Page 3: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Ch 2: Problem 4

.:.:~~--~~-

50 Cartesian Tensors

The indicial notation avoids all the problems mentioned in the preceding. The

algebraic manipulations are especially simple. The ordering of terms is unneces-

sary because A¡jBkl means the same thing as BkiAij. In this notation we deal with

components only, which are scalars. Another major advantage is that one does not

have to remember formulas except for the product eijkeklm, which is given by equa-

tion (2. 19). The disadvantage of the indicial notation is that the physical meaning of a

term becomes clear only after an examination of the indices. A second disadvantage

is that the cross product involves the introduction of the cumbersome e¡jk. This, how-ever, can frequently be avoided by writing the i-component of the vector product of 0

and v as (0 x v)¡ using a mixture of boldface and indicial notations. In this book we

shall use boldface, indicial and mixed notations in order to take advantage of each. As

the reader might have guessed, the algebraic manipulations wil be performed mostly

in the indicial notation, sometimes using the comma notation.

Exercises1. Using indicial notation, show that

.Í a x (b x c) = (a. c)b - (a. b)c.

(Hint: Call d == b x c. Then (a x d)m = epqmapdq = epqmape¡jqb¡cj. Using

equation (2.19), show that (a x d)m = (a. c)bm - (a. b)cm.)

2. Show that the condition for the vectors a, b, and c to be coplanar is

e¡jkaibjCk = O.

3. Prove the following relationships:

OijOij = 3epqrepqr = 6epq¡epqj = 2oij.

4. Show thatc . CT = CT . C = 8,

,.;

where C is the direction cosine matrx and 8 is the matrx of the Kronecker delta.

Any matrx obeying spch a relationship is called an orthogonal matri because it

represents transformation of one set of ortogonal axes into another.

5. Show that for a second-order tensor A, the following thee quantities areinvarant under the rotation of axes:

II = A¡¡

/z = I Aii A12 I + I A.2 A23/ + I Aii Al31A21 A22 A32 A33 A31 A33

13 = det(A¡j).

3

Page 4: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Ch 2: Problem 4

.:.:~~--~~-

50 Cartesian Tensors

The indicial notation avoids all the problems mentioned in the preceding. The

algebraic manipulations are especially simple. The ordering of terms is unneces-

sary because A¡jBkl means the same thing as BkiAij. In this notation we deal with

components only, which are scalars. Another major advantage is that one does not

have to remember formulas except for the product eijkeklm, which is given by equa-

tion (2. 19). The disadvantage of the indicial notation is that the physical meaning of a

term becomes clear only after an examination of the indices. A second disadvantage

is that the cross product involves the introduction of the cumbersome e¡jk. This, how-ever, can frequently be avoided by writing the i-component of the vector product of 0

and v as (0 x v)¡ using a mixture of boldface and indicial notations. In this book we

shall use boldface, indicial and mixed notations in order to take advantage of each. As

the reader might have guessed, the algebraic manipulations wil be performed mostly

in the indicial notation, sometimes using the comma notation.

Exercises1. Using indicial notation, show that

.Í a x (b x c) = (a. c)b - (a. b)c.

(Hint: Call d == b x c. Then (a x d)m = epqmapdq = epqmape¡jqb¡cj. Using

equation (2.19), show that (a x d)m = (a. c)bm - (a. b)cm.)

2. Show that the condition for the vectors a, b, and c to be coplanar is

e¡jkaibjCk = O.

3. Prove the following relationships:

OijOij = 3epqrepqr = 6epq¡epqj = 2oij.

4. Show thatc . CT = CT . C = 8,

,.;

where C is the direction cosine matrx and 8 is the matrx of the Kronecker delta.

Any matrx obeying spch a relationship is called an orthogonal matri because it

represents transformation of one set of ortogonal axes into another.

5. Show that for a second-order tensor A, the following thee quantities areinvarant under the rotation of axes:

II = A¡¡

/z = I Aii A12 I + I A.2 A23/ + I Aii Al31A21 A22 A32 A33 A31 A33

13 = det(A¡j).

Multiplication and Contraction of Tensors

4

Page 5: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Ch 2: Problem 4

.:.:~~--~~-

50 Cartesian Tensors

The indicial notation avoids all the problems mentioned in the preceding. The

algebraic manipulations are especially simple. The ordering of terms is unneces-

sary because A¡jBkl means the same thing as BkiAij. In this notation we deal with

components only, which are scalars. Another major advantage is that one does not

have to remember formulas except for the product eijkeklm, which is given by equa-

tion (2. 19). The disadvantage of the indicial notation is that the physical meaning of a

term becomes clear only after an examination of the indices. A second disadvantage

is that the cross product involves the introduction of the cumbersome e¡jk. This, how-ever, can frequently be avoided by writing the i-component of the vector product of 0

and v as (0 x v)¡ using a mixture of boldface and indicial notations. In this book we

shall use boldface, indicial and mixed notations in order to take advantage of each. As

the reader might have guessed, the algebraic manipulations wil be performed mostly

in the indicial notation, sometimes using the comma notation.

Exercises1. Using indicial notation, show that

.Í a x (b x c) = (a. c)b - (a. b)c.

(Hint: Call d == b x c. Then (a x d)m = epqmapdq = epqmape¡jqb¡cj. Using

equation (2.19), show that (a x d)m = (a. c)bm - (a. b)cm.)

2. Show that the condition for the vectors a, b, and c to be coplanar is

e¡jkaibjCk = O.

3. Prove the following relationships:

OijOij = 3epqrepqr = 6epq¡epqj = 2oij.

4. Show thatc . CT = CT . C = 8,

,.;

where C is the direction cosine matrx and 8 is the matrx of the Kronecker delta.

Any matrx obeying spch a relationship is called an orthogonal matri because it

represents transformation of one set of ortogonal axes into another.

5. Show that for a second-order tensor A, the following thee quantities areinvarant under the rotation of axes:

II = A¡¡

/z = I Aii A12 I + I A.2 A23/ + I Aii Al31A21 A22 A32 A33 A31 A33

13 = det(A¡j).

5

Page 6: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Ch 2: Problem 5.:.:~~--~~-

50 Cartesian Tensors

The indicial notation avoids all the problems mentioned in the preceding. The

algebraic manipulations are especially simple. The ordering of terms is unneces-

sary because A¡jBkl means the same thing as BkiAij. In this notation we deal with

components only, which are scalars. Another major advantage is that one does not

have to remember formulas except for the product eijkeklm, which is given by equa-

tion (2. 19). The disadvantage of the indicial notation is that the physical meaning of a

term becomes clear only after an examination of the indices. A second disadvantage

is that the cross product involves the introduction of the cumbersome e¡jk. This, how-ever, can frequently be avoided by writing the i-component of the vector product of 0

and v as (0 x v)¡ using a mixture of boldface and indicial notations. In this book we

shall use boldface, indicial and mixed notations in order to take advantage of each. As

the reader might have guessed, the algebraic manipulations wil be performed mostly

in the indicial notation, sometimes using the comma notation.

Exercises1. Using indicial notation, show that

.Í a x (b x c) = (a. c)b - (a. b)c.

(Hint: Call d == b x c. Then (a x d)m = epqmapdq = epqmape¡jqb¡cj. Using

equation (2.19), show that (a x d)m = (a. c)bm - (a. b)cm.)

2. Show that the condition for the vectors a, b, and c to be coplanar is

e¡jkaibjCk = O.

3. Prove the following relationships:

OijOij = 3epqrepqr = 6epq¡epqj = 2oij.

4. Show thatc . CT = CT . C = 8,

,.;

where C is the direction cosine matrx and 8 is the matrx of the Kronecker delta.

Any matrx obeying spch a relationship is called an orthogonal matri because it

represents transformation of one set of ortogonal axes into another.

5. Show that for a second-order tensor A, the following thee quantities areinvarant under the rotation of axes:

II = A¡¡

/z = I Aii A12 I + I A.2 A23/ + I Aii Al31A21 A22 A32 A33 A31 A33

13 = det(A¡j).

Cartesian Tensors

~ preceding. The:rms is unneces-on we deal withiat one does nots given by equa-cal meaning of and disadvantageeijk. This, how-;tor product of uIn this book wetage of each. As:rformed mostly

?ijqbiCj. Using

mar is

Jnecker delta.ri because it

quantities 'are'w

I

Supplemental Readitig 51

(Hint; Use the result of Exercise 4 and the transformation rule (2.12) to show thati; = A;i = Aii = ii. Then show that AijAji and AijAjkAki are also inva

rants. In

fact, all contracted scalars of the form Aij A jk . . . Ami are invariants. Finally, verify

that

/z = l(ll - AijAji)

h = AijAjkAki - !¡AijAji + /zAii.

Because the right-hand sides are invarant, so are /z and h)

6. If u and v are vectors, show that the products Ui v j obey the transformationrule (2.12), and therefore represent a second-order tensor.

7. Show that oij is an isotropic tensor. That is, show that O;j = Oij under rotation

of the coordinate system. (Hint: Use the transformation rule (2.12) and the results of

Exercise 4.)

8. Obtain the recipe for the gradient of a scalar function in cylindrical polarcoordinates from the integral definition.

9. Obtan the recipe for the divergence of a vector in spherical polar coordinates

from the integral definition.

10. Proye that div(~urlii) = 0 for any vector u regardless of the coordinate

system. (Hit: ùse thë veCtòr integral theorems.)

i 1. Prove that curl (grad Ø) = 0 for any single-valued scalar ø regardless of thecoordinate system. (Hint: use Stokes' theorem.)

Literature Cited

Sommedeld, A. (1964). Mechanics of Deformble Bodies, New York: Academic Press. (Chapter I contansbrief but useful coverage of Caresian tensors.)

Supplemental ReadingArs, R. (1962). Vectors, Tensors, and the Basic Equations of Fluid Mechanics, Englewood Cliffs, NJ:

Prentice-Hal. (Ts book gives a clear and easy treatment of tensors in Caresian and non-Caresian

coordinates, with applications to fluid mechancs.)Prager, W. (1961). Introduction to Mechanics.of Continua, New York: Dover Publications. (Chapters 1

and 2 contan brief but useful coverage of Caresian tensors.)

6

Page 7: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Ch 2: Problem 5

.:.:~~--~~-

50 Cartesian Tensors

The indicial notation avoids all the problems mentioned in the preceding. The

algebraic manipulations are especially simple. The ordering of terms is unneces-

sary because A¡jBkl means the same thing as BkiAij. In this notation we deal with

components only, which are scalars. Another major advantage is that one does not

have to remember formulas except for the product eijkeklm, which is given by equa-

tion (2. 19). The disadvantage of the indicial notation is that the physical meaning of a

term becomes clear only after an examination of the indices. A second disadvantage

is that the cross product involves the introduction of the cumbersome e¡jk. This, how-ever, can frequently be avoided by writing the i-component of the vector product of 0

and v as (0 x v)¡ using a mixture of boldface and indicial notations. In this book we

shall use boldface, indicial and mixed notations in order to take advantage of each. As

the reader might have guessed, the algebraic manipulations wil be performed mostly

in the indicial notation, sometimes using the comma notation.

Exercises1. Using indicial notation, show that

.Í a x (b x c) = (a. c)b - (a. b)c.

(Hint: Call d == b x c. Then (a x d)m = epqmapdq = epqmape¡jqb¡cj. Using

equation (2.19), show that (a x d)m = (a. c)bm - (a. b)cm.)

2. Show that the condition for the vectors a, b, and c to be coplanar is

e¡jkaibjCk = O.

3. Prove the following relationships:

OijOij = 3epqrepqr = 6epq¡epqj = 2oij.

4. Show thatc . CT = CT . C = 8,

,.;

where C is the direction cosine matrx and 8 is the matrx of the Kronecker delta.

Any matrx obeying spch a relationship is called an orthogonal matri because it

represents transformation of one set of ortogonal axes into another.

5. Show that for a second-order tensor A, the following thee quantities areinvarant under the rotation of axes:

II = A¡¡

/z = I Aii A12 I + I A.2 A23/ + I Aii Al31A21 A22 A32 A33 A31 A33

13 = det(A¡j).

7

Page 8: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Ch. 2 Question 10

8

Page 9: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Ch. 2 Question 10

9

Page 10: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Ch. 2 Question 10

10

Page 11: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Ch. 2 Question 11

11

Page 12: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Ch. 2 Question 11

12

Page 13: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Ch. 2 Question 11

13

Page 14: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Kinematics of Fluid Flow (Ch. 3)

1. Streamlines, pathlines, and convective (material) derivative

2.Translations, Deformation, and Rotation of a fluid element (Cauchy-Stokes Decomposition)

14

Page 15: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Kinematics of Fluid Flow (Ch. 3)

15

Page 16: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Eulerian Description (Field Variables)

16

Page 17: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Time Derivatives

“advective” + “unsteady”

17

Page 18: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Streamlines and Pathlines

18

Page 19: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Streamlines and Pathlines (cont.)

2 4 6 8 10 12 14x

1

2

3

4

5

6

7

8

yStreamlines !Blue, t!4" and Pathlines !Red, 0"t"4"

19

Page 20: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Streaklines or “Smoke Lines”~

1

J-\!~

0~L

't

F-t

\(

"~

~'l \ ~

~~

..l~

~~

i.2

,.~

~1\

~~

~l

l~~

~

Q

--

,- .. ~ --.l

'-Li

'l~

\

:j

~i

1t

t()

L/"

.j..

il.J

J'\

)-.0a

' iJ-

~"'

lL\)

)~

(J.i)(0

~

~.j

rJ~

'-"

.."'

~V

\\-

~

~('

~~

II

-.'"

/'~~

luV

\'t

. .'- '"

\J

~

1Jt-

r:.r

4:~

I\

()

~ \-/

4-t

l ~

l-~

tJ

çt

..~

--

Qt

t'"

~

:co

It,$

+-~

V

\f ~:\

l--

.-.-

()

.j

VI

'- \.lI

~~

t~

-L

\I~u

0~

VI

~'v

il,..

y.)~

l-Ù~

lLlLl

\t'-

~-

V1

1l~

~l\

(j3

--l

~V

I

1\"

L \

.l.J

~t;

~..

\J

~lU

Jl!j

\iJ-4~

,-

~l.

..LL

-l~

t--

~

~\)~~

C'

ilyJ

0x:

)-lU

lU1-

M't

l~

~'t

j"

~V

lQ

:V

1C

-~

..2-

J.;

V'

20

Page 21: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Deformation and Flow: Translation, Stretching, Pinching, and Rotating

21

Page 22: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Simple Comments about Velocity Gradient Tensor

22

Page 23: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Simple Example:

Expansion/Contraction

Vector Notation

23

Page 24: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Stretching along one Axis

24

Page 25: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Rotation

25

Page 26: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Bending, Distorting, Shearing

26

Page 27: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Example Shearing

27

Page 28: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Simple Example #1

28

Page 29: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Simple Example #1(from a different point of view)

29

Page 30: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Simple Example #2

30

Page 31: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Simple Example #2(from a different point of view)

-

31

Page 32: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Principal Axes for

Symmetric Tensors

32

Page 33: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Example Principal Axes #2

33

Page 34: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Example Principal Axes #1

34

Page 35: APPH 4200 Physics of Fluids - Columbia Universitysites.apam.columbia.edu/courses/apph4200x/Lecture-4.pdf · APPH 4200 Physics of Fluids ... If u and v are vectors, ... Tensors, and

Summary: Cauchy-Stokes Decomposition

35