Vector And Tensor Analysis Notes 2c63c

1

Ivan Avramidi, MATH 332: Vector and Tensor Analysis, Formulas

MATH 332: Vector and Tensor Analysis Vector Algebra

Symmetry Unity Matrix

A · B = |B||A| cos θ

co m

δjk = δkj

Scalar Product



 1 0 0 I = (δik ) =  0 1 0  0 0 1

Commutativity

e.

A·B=B·A Magnitude

re

Transformation of Rectangular Coordinates 2

llo

A × B = −B × A

nl

Anti-commutativity

in

Vector Product A × B = n | A | | B | sin θ

x0i = αi0 k xk + x(0) i

ef

|A| = A · A

.a

A×A=0

w

Scalar Triple Product

w

(A × B) · C

w

Cyclic symmetry

(A × B) · C = (B × C) · A = (C × A) · B

Summation from 1 to 3 is assumed over all repeated indices

Vector Form  x01 x0 = (x0i ) =  x02  , x03 

 x1 x = (xi ) =  x2  x3 

Transformation matrix 

 α10 1 α10 2 α10 3 α = (αi0 k ) =  α20 1 α20 2 α20 3  α30 1 α30 2 α30 3

Matrix Form of the Transformation x0 = αx + x0

Vector Triple Product A × (B × C) = B(A · C) − C(A · B)

Orthonormal Basis i k ,

(k = 1, 2, 3)

ik · ij = δkj Kronecker Symbol δik =

δki

=

δik

=

1 if i = k 0 if i 6= k

Orientation of Orthonormal Basis: right-handed if ( i 1 × i 2 ) · i 3 = +1 left-handed if ( i 1 × i 2 ) · i 3 = −1

2

Ivan Avramidi, MATH 332: Vector and Tensor Analysis, Formulas

i 0j = αj 0 k i k ,

Levi-Civita (Alternating) Symbol  +1 if (i, j, k) = (1, 2, 3),     (2, 3, 1), (3, 1, 2)   −1 if (i, j, k) = (2, 1, 3), εijk =   (3, 2, 1), (1, 3, 2)     0 otherwise

i k = αj 0 k i 0j

αj 0 k = cos( i 0j , i k ) Orthogonality condition αi0 k αj 0 k = δij ,

αk0 i αk0 j = δij

co m

Transformation of Orthonormal Basis

εijk = εijk

Matrix Form of Orthogonality Condition

Anti-symmetry

εijk = −εjik = −εikj = −εkji

e.

T means transposition of a matrix (replacement of rows by columns)

Cyclic symmetry

re

ααT = I

Proper transformation (no change of orientation)

ef

det(αj 0 k ) = 1

εijk = εjki = εkij

Orthonormal basis

llo

Cartesian Vectors

nl

det(αj 0 k ) = −1

in

Improper transformation (changes orientation)

Ak = A · i k

.a

A = Ak i k ,

i j × i k = εjkl i l ( i j × i k ) · i l = εjkl

Vector Product in Cartesian Components (A × B)i = εijk Aj Bk

Scalar Product in Cartesian Components

w

Scalar Triple Product

A · B = Ai Bi

w

(A × B) · C = εijk Ai Bj Ck

w

|A|2 = Ai Ai

Vector Product

Tensor Notation



 i1 i2 i3 A × B = det  A1 A2 A3  B1 B2 B3

δii = δii = 3

Scalar Triple Product 



A1 A2 A3  (A × B) · C = det B1 B2 B3  C1 C2 C3

δij Aj = Ai δ ij Ai Bj = Ai Bi = A · B

3

Ivan Avramidi, MATH 332: Vector and Tensor Analysis, Formulas

A0 = αA εijj = εjij = εjji = 0

Matrix Form of a 2-tensor   A11 A12 A13 A = (Aik ) =  A21 A22 A23  A31 A32 A33

εijk δij = εijk δik = εijk δjk = 0 εijk Aj Ak = εijk Ai Ak = εijk Ai Aj = 0

co m

A0 = αAαT εijk εmnl = δim δjn δkl + δjm δkn δil + δkm δin δjl −

δim δkn δjl

−

δjm δin δkl

−

Stress Tensor pik Stress

δkm δjn δil

pi = pik nk

εijk εmnk = δim δjn − δjm δin

e.

ni unit normal Moment of Inertia Tensor

re

εijk εmjk = 2δim

Iik =

ef

εijk εijk = 6

N X

h i (j) (j) (j) (j) mj δik xl xl − xi xk

j=1

in

Angular Momentum

(tensors in rectangular coordinates)

w

.a

llo

Scalar: 0-tensor ϕ Vector: 1-tensor Ai 2-tensor Aik n-tensor Ai1 ...in

nl

Cartesian Tensors

Transformation Laws

w

w

ϕ0 = ϕ

A0i = αi0 k Ak

A0ij = αi0 k αj 0 l Akl A0i1 ...in

= αi01 j1 · · · αi0n jn Aj1 ...jn

Matrix Form of the Transformation Laws Vector   A1 A = (Ai ) =  A2  A3

Li = Iik ωk

ωk angular velocity Deformation Tensor ∂uk ∂ul ∂ul 1 ∂ui + + uik = 2 ∂xk ∂xi ∂xi ∂xk ui displacement vector Rate of Deformation Tensor 1 ∂vi ∂vk vik = + 2 ∂xk ∂xi vi velocity vector field Isotropic Tensors (built from δik only; no preferred directions) Isotropic 2-tensor Aik = pδik

4

Ivan Avramidi, MATH 332: Vector and Tensor Analysis, Formulas

Isotropic 4-tensor

Antisymmetrization 1 T[ij] = (Tik − Tki ) 2

Aiklm = pδik δlm + ρδil δkm + λδim δkl

Decomposition of 2-Tensor

Tensor Algebra

Tik = T(ik) + T[ij]

Cik = Ai Bk ,

co m

Tensor Product

Duality (Equivalence of antisymmetric 2-tensor to an axial vector)

Cijkl = Aij Bkl

Contraction Djkl = Ciijkl

e.

1 A˜i = εijk Ajk , 2

Trace

re

A = Aii

A˜1 = A23 ,

Ci = Aik Bk ,

Dij = Cijkl Bkl

A˜2 = A31 ,

A˜3 = A12 ,

ef

Inner Product

Aij = εijk A˜k

in

Principal Axes

nl

Symmetric Tensors

.a

llo

Sij = Sji   S11 S12 S13 (Sik ) =  S12 S22 S23  S13 S23 S33

| n (r) | = 1 Characteristic Equation

w

Anti-symmetric Tensors

Eigenvalues λ(r) (characteristic values) and Eigenvectors n (r) (characteristic or principal directions) (r) (r) Tik nk = λ(r) ni

det(Tik − λδik ) = 0

w

Aij = −Aji

 0 A12 A13 0 A23  (Aik ) =  −A12 −A13 −A23 0

w



Contraction of antisymmetric tensor Aii = 0 Symmetrization 1 T(ij) = (Tik + Tki ) 2



 T11 − λ T12 T13 T22 − λ T23  = 0 det  T21 T31 T32 T33 − λ λ3 − I1 λ2 + I2 λ − I3 = 0 Invariants of a 2-Tensor: Ik I1 = Tii = T11 + T22 + T33 T22 T23 T11 T12 I2 = det + det T32 T33 T21 T22

5

Ivan Avramidi, MATH 332: Vector and Tensor Analysis, Formulas

Curvilinear Coordinates



T11 T13 + det T31 T33   T11 T12 T13 I3 = det(Tik ) = det  T21 T22 T23  T31 T32 T33

q i = q i ( r ),

(i = 1, 2, 3)

Cartesian Coordinates xi = xi (q)

co m

The eigenvalues λr , (r = 1, 2, 3), of a symmetric Radius vector 2-tensor are real

r(q) = xk (q) i k = x1 (q) i 1 + x2 (q) i 2 + x3 (q) i 3 A symmetric 2-tensor has three orthogonal principal axes n(r) , (r = 1, 2, 3) Basis (tangent vectors to coordinate curves)

e.

n (r) · n (p) = δrp

∂r ∂q i

re

ei =

0

e k = αj k e 0j 0

αi k0 αk j = δji

nl

in

ef

In the principal axes a symmetric 2-tensor Tik has Transformation of basis diagonal matrix e 0j = αk j 0 e k ,   λ1 0 0 0 (Tik ) =  0 λ2 0  αi k αk j 0 = δji , 0 0 λ3

.a

Tik =

3 X

llo

Orientation: Decomposition in of orthonormal eigenvecright-handed if ( e 1 × e 2 ) · e 3 > 0 tors n (r) and eigenvalues λ(r) and left-handed if ( e 1 × e 2 ) · e 3 < 0 (r) (r)

λ(r) ni nk

w

w

r=1

w

Traceless Tensors (Deviators) 1 T ik = Tik − δik T 3

where T = Tjj (trace) T = T ii = 0 Decomposition

Reciprocal Basis e j · e k = δkj ej =

ek × el ( e 1 × e 2) · e 3

where (j, k, l) = (1, 2, 3), (2, 3, 1), (3, 1, 2) Contravariant Components A = Ai e i ,

Ai = A · e i

Covariant components 1 Tik = T ik + δik Tjj 3

A = Ai e i ,

Ai = A · e i

6

Ivan Avramidi, MATH 332: Vector and Tensor Analysis, Formulas

ds2 = gik dq i dq k

Transformation of components 0

A0i = αi k Ak ,

Volume Element √ dV = G dq 1 dq 2 dq 3 ,

A0i = αk i0 Ak

G = det(gik )

Metric Tensor ∂r ∂r · , ∂q i ∂q k

g ik = e i · e k

Contravariant components Aik Covariant components Aik Mixed components Ai k , Ai k Relations

Symmetry

e.

g ik = g ki

gik = gki ,

Tensors in Curvilinear Coordinate System

co m

gik = e i · e k =

Aik = gin An k = gkn Ai n = gin gkm Anm

re

gik g kn = δin

Aik = g in An k = g kn Ai n = g in g km Anm

det(g ik ) =

1 G

Ai k = g in Ank = gkn Ain

etc

in

G = det(gik ) ,

ef

Determinant of the Metric Tensor

ik

det M ik G

llo

g ik = (−1)ik

nl

The matrix g is the inverse matrix of the matrix Orthogonal Coordinate System gik given by

w

e i · e k = e i · e k = 0 if i 6= k | e i | = hi ,

w

Relations between components Ai = g ik Ak

w

Ai = gik Ak ,

if i 6= k

Basis

.a

where M ik is a 2 × 2 matrix obtained from the 3 × 3 matrix gij by removing the i-th row and the k-th column

∂r ∂r · = 0, ∂q i ∂q k

| e i| =

1 hi

Metric gik = 0,

if i 6= k

gii = h2i (no summation!)

Displacement dr =

∂r i dq ∂q i

Line Element (Arc Length) ds2 = d r · d r =

∂r ∂r i k · dq dq ∂q i ∂q k

g ik = 0, if i 6= k 1 g ii = 2 (no summation!) hi Metric Coefficients ∂r hi = i ∂q

7

Ivan Avramidi, MATH 332: Vector and Tensor Analysis, Formulas

s hi =

∂x1 ∂q i

2

+

∂x2 ∂q i

2

+

∂x3 ∂q i

x1 = ρ cos ϕ,

2

r = ρ cos ϕ i 1 + ρ sin ϕ i 2 + z i 3 Orthonormal basis

(no summation!)

e ρ = cos ϕ i 1 + sin ϕ i 2

co m

Ai =

x3 = z

Radius Vector

Relation between covariant and contravariant components h2i Ai

x2 = ρ sin ϕ,

e ϕ = − sin ϕ i 1 + cos ϕ i 2 Orthonormal Basis

ez = i3

∂r ∂q i e i = e i = ∂r ∂q i

e.

Line Element

re

ds2 = (dρ)2 + ρ2 (dϕ)2 + (dz)2

Line Element ds =

h21

1 2

(dq ) +

h22

2 2

(dq ) +

h23

3 2

(dq )

nl

Volume Element

dV = ρ dρ dϕ dz

Metric Coefficients

in

2

ef

Volume Element

hρ = 1,

hϕ = ρ,

hz = 1

.a

Cartesian Coordinates Line Element

llo

dV = h1 h2 h3 dq 1 dq 2 dq 3

w

Volume Element

w

ds2 = (dx1 )2 + (dx2 )2 + (dx3 )2

w

dV = dx1 dx2 dx3

Spherical Coordinates r ≥ 0,

0 ≤ ϕ < 2π p q x21 + x22 2 2 2 r = x1 + x2 + x3 , tan θ = , x3 x2 tan ϕ = x1 x1 = r sin θ cos ϕ, x2 = r sin θ sin ϕ, x3 = r cos θ

Metric Coefficients h1 = h2 = h3 = 1

0 ≤ θ ≤ π,

Radius Vector r = r sin θ cos ϕ i 1 + r sin θ sin ϕ i 2 + r cos θ i 3

Cylindrical Coordinates ρ ≥ 0, 0 ≤ ϕ < 2π, −∞ < z < ∞ q x2 ρ = x21 + x22 , tan ϕ = , z = x3 x1

Orthonormal Basis e r = sin θ cos ϕ i 1 + sin θ sin ϕ i 2 + cos θ i 3 e θ = cos θ cos ϕ i 1 + cos θ sin ϕ i 2 − sin θ i 3

8

Ivan Avramidi, MATH 332: Vector and Tensor Analysis, Formulas

e ϕ = − sin ϕ i 1 + cos ϕ i 2

Speed dr ds |v| = = dt dt

Line Element ds2 = (dr)2 + r2 (dθ)2 + r2 sin2 θ(dϕ)2

Acceleration

Volume Element

a=

dV = r2 sin θ dr dθ dϕ

co m

Arc Length (Line Element)

Metric Coefficients

ds = |v|dt

hθ = r,

hϕ = r sin θ

e.

hr = 1,

dv d2 r = 2 dt dt

Fields in Cartesian Coordinates

re

Partial derivatives

ef

Vector And Tensor Analysis Functions of Single variable

∂i =

∂ ∂xi

in

Nabla (Del) Operator

nl

Product Rules

llo

dϕ dA d (ϕ A ) = A +ϕ dt dt dt

Gradient grad f = ∇f = ik ∂k f

w

Trajectory

w

w

.a

d dB dA (B · A) = ·A+B· dt dt dt dB dA d (B × A) = ×A+B× dt dt dt

∇ = ik ∂k = i1 ∂1 + i2 ∂2 + i3 ∂3

r(t) = xk (t) ik ,

a≤t≤b

Velocity dr v= dt Unit Tangent Vector dr dr u = = dt ds dr dt

( grad f )i = ∂i f Directional Derivative df dr = · gradf = u · gradf ds ds Flow Lines

dr =βF dt dx1 dx2 dx3 = = F1 F2 F3

Divergence div F = ∇ · F = ∂i Fi

9

Ivan Avramidi, MATH 332: Vector and Tensor Analysis, Formulas

∇×r=0 q r = | r | = x21 + x22 + x23

Curl i1 i2 i3 curl F = ∂1 ∂2 ∂3 F1 F2 F3

r , r

∇r =

(curl F)i = εijk ∂j Fk

|∇r| = 1

df r dr r r dF ∇ · F (r) = · r dr r dF ∇ × F (r) = × r dr

curl F = ∇ × F = il εljk ∂j Fk

co m

∇f (r) =

Laplacian ∆ = div grad = ∇2 = ∇ · ∇

e.

∆ = ∂i ∂i = ∂12 + ∂22 + ∂32

re

(F · ∇)r = F

in

∇ × (∇ × F ) = −∆ F + ∇(∇ · F )

Fields in Orthogonal Coordinate System (in orthonormal basis e i )

ef

Vector identities:

∇ × ∇ϕ = 0,

llo

curl grad = 0

Vector Components

nl

curl curl F = −∆ F + grad div F

Fi = F · e i

∇ · (∇ × F) = 0,

.a

div curl = 0

w

∇(f g) = (∇f )g + f (∇g)

grad f = e 1

1 ∂ 1 ∂ 1 ∂ f + e2 f + e3 f 1 2 h1 ∂q h2 ∂q h3 ∂q 3

w

∇(f F) = (∇f ) · F + f ∇F ∇ × (f F) = (∇f ) × F + f (∇ × F)

w

∇( F × G ) = G · (∇ × F ) − F · (∇ × G ) df ∇ϕ dϕ dF ∇ · F (ϕ) = ∇ϕ · dϕ dF ∇ × F (ϕ) = ∇ϕ × dϕ ∇f (ϕ) =

∇r = 3

1 div F = h1 h2 h3

(

∂ (h2 h3 F1 ) ∂q 1 ) ∂ ∂ + 2 (h3 h1 F2 ) + 3 (h1 h2 F3 ) ∂q ∂q

1 ∂ ∂ curl F = e 1 (h3 F3 ) − [ 3 (h2 F2 ) h2 h3 ∂q 2 ∂q 1 ∂ ∂ +e2 (h1 F1 ) − [ 1 (h3 F3 ) h3 h1 ∂q 3 ∂q

10

Ivan Avramidi, MATH 332: Vector and Tensor Analysis, Formulas

1 ∂ ∂ +e3 (h2 F2 ) − [ 2 (h1 F1 ) h1 h2 ∂q 1 ∂q ∆f = (

∂ h2 h3 ∂ 1 f ∆f = h1 h2 h3 ∂q 1 h1 ∂q 1 ) h2 h3 ∂ h2 h3 ∂ ∂ ∂ + 1 f + 1 f ∂q h1 ∂q 1 ∂q h1 ∂q 1

1 1 ∂θ (sin θ∂θ f ) ∂r (r2 ∂r f ) + 2 2 r r sin θ +

1 ∂2 f sin2 θ ϕ

co m

Integrals

Parametrization of a curve C

Cylindrical Coordinates:

r = r (t),

re

Line Integrals

1 1 div F = ∂ρ (ρFρ ) + ∂ϕ Fϕ + ∂z Fz ρ ρ

Zb F · dr =

F(x(t)) ·

dr dt dt

a

C

in

ef

Z

llo

nl

eρ ρeϕ ez 1 curl F = ∂ρ ∂ϕ ∂z ρ Fρ ρFϕ Fz

a≤t≤b

e.

1 grad f = eρ ∂ρ f + eϕ ∂ϕ f + ez ∂z f ρ

1 1 ∆f = ∂ρ (ρ ∂ρ f ) + 2 ∂ϕ2 f + ∂z2 f ρ ρ

Circulation of vector field along a closed contour I F · dr C

.a

Parametrization of a surface S r = r(u, v),

w

Spherical Coordinates

w

a ≤ u ≤ b,

w

1 1 grad f = er ∂r f + eθ ∂θ f + eϕ ∂ϕ f r r sin θ

div F =

r2

e reθ r sin θ eϕ 1 r ∂r ∂θ ∂ϕ curl F = 2 r sin θ Fr rFθ r sin θFϕ

Unit Normal n=

1 1 ∂r (r2 Fr ) + ∂θ (sin θ Fθ ) 2 r r sin θ 1 + ∂ϕ Fϕ r sin θ

c≤v≤d

∂u r × ∂v r |∂u r × ∂v r|

For a surface given by f( r ) = C the Unit Normal is n=

∇f |∇f |

11

Ivan Avramidi, MATH 332: Vector and Tensor Analysis, Formulas

Surface Element

Circulation of a Gradient along a closed contour I grad ϕ · d r = 0

dS = n dS = ∂u r × ∂v r du dv

C

d S = |∂u r × ∂v r| du dv For a surface given by z = f (x, y) a ≤ x ≤ b, y1 (x) ≤ y ≤ y2 (x)

D

the Surface Element is q dS = 1 + (∂x f )2 + (∂y f )2 dy dx

e.

re

Green’s Theorem ZZ I (∂1 F2 − ∂2 F1 )dx1 dx2 = (F1 dx1 + F2 dx2 )

ef S

Zd Zb

ϕ( r (u, v)) |∂u r × ∂v r| du dv c

a

nl

S

llo

Zb yZ2 (x) = ϕ(x, y, z(x, y))

.a

a y1 (x)

q

w

1 + (∂x f )2 + (∂y f )2 dy dx

w

w

Flux of a Vector Field F through the surface S ZZ ZZ F · dS = F · n dS S

∂S

∂S is a closed plane curve, which is the boundary of the region S in the x1 x2 -plane

in

ϕ dS =

∂D

∂D is a closed surface, which is the boundary of the solid region D

Surface Integral of a scalar field ZZ

co m

Gauss (Divergence) Theorem ZZZ ZZ div F dV = F · dS

S

Stokes’ Theorem ZZ I curl F · dS = F · dr S

∂S

∂S is a closed space curve, which is the boundary of the surface S Flux of a curl through a closed surface S ZZ curl F · dS = 0 S

Line Integral of a Gradient

Tensor Fields

ZQ grad ϕ · d r = ϕ(Q) − ϕ(P ) P

Flux of a Tensor Field ZZ Tik nk dS, S

ZZ Tik ni dS S

12

Ivan Avramidi, MATH 332: Vector and Tensor Analysis, Formulas

Divergence of a Tensor (in Cartesian Coordinates) ∂i Tik , ∂k Tik

co m

Directional Derivative (in Cartesian Coordinates) dTik dxj = ∂j Tik ds ds Analog of curl of an antisymmetric 2-tensor 3εijk ∂i Ajk = ∂1 A23 + ∂2 A31 + ∂3 A12

w

w

w

.a

llo

nl

in

ef

re

e.

˜ = ∂1 A˜1 + ∂2 A˜2 + ∂3 A˜3 = div A