UNIVERSITY OF NOTRE DAME

DEPARTMENT OF AEROSPACE AND MECHANICAL ENGINEERING

Professor H.M. Atassi AME-60612

113 Hessert Center Mathematical Methods II
Tel: 631-5736
Email: [email protected]

Orthogonal Curvilinear Coordinates

1 Definitions
Let x = (x1 , x2 , x3 ) be the Cartesian coordinates of a point M with respect to a frame
of reference defined by the unit vectors e1 , e2 , e3 , . We introduce three functions defined
by
uj = uj (x1 , x2 , x3 ), j = 1, 3, (1)
in a region R. The equation uj = cj , where cj is a constant, represents a surface. The
system of two equations u2 = c2 and u3 = c3 represent a line γ1 where the two surfaces
intersect. Along γ1 , only u1 varies. The system of three equations u1 = c1 , u2 = c2 and
u3 = c3 represent a point where the three surfaces intersect. At every point M ∈ R,
there are three lines γi (ui ). For Cartesian coordinates, these surfaces are planes. For
cylindrical coordinates, we define

1
u1 = r = (x21 + x22 ) 2 , (2)
x2
u2 = θ = tan−1 ( ), (3)
x1
u 3 = x3 . (4)
(5)

Here r = c1 represents a circular cylinder of radius c1 , θ = c2 represents a vertical plane,

and x3 = c3 represents a horizontal plane. The two equations r = c1 and x3 = c3 repre-
sent a circle in a horizontal plane, only θ varies as we move along the circle.

The position vector of a point M can be expressed in the Cartesian system as

−−→
OM = x = xi ei , (6)

1
where the repeated index implies summation, i.e., xi ei = x1 e1 + x2 e2 + x3 e3 . Note that
∂x
= ei . (7)
∂xi
We now want to use uj as a new coordinate system. We assume that the Cartesian
coordinates xi are given in terms of the new coordinates uj ,

xi = xi (u1 , u2 , u3 ), i = 1, 3.. (8)

Differentiating x with respect to uj , we get

∂x ∂xi
dx = duj = ei duj . (9)
∂uj ∂uj

The vector
∂x ∂xi
Êj = = ei (10)
∂uj ∂uj
is tangent to γj . Note that
∂x ∂sj
Êj = , (11)
∂sj duj
∂x
where ∂sj is the elementary arc length along γj . We also note that the vector Ej = ∂sj
∂sj
is a unit vector. Thus if hj = ∂uj
, along γj , dsj = hj duj and Êj = hj Ej . Hence, using
(10), we get
∂xi
hj Ej = ei . (12)
∂uj
Since both ei and Ej are orthonormal vectors,
i=3
∂xi 2
h2j =
X
( ). (13)
i=1 ∂uj
Equation(13) defines the three scales associated with the new coordinates system.

2 Elementary Quantities
2.1 Elementary Arc Length
The elementary arc length of a line, not coinciding with the three lines defining the
coordinate system at a point M, is obtained by taking the magnitude of (9),

(ds)2 = h2j (duj )2 . (14)

2
2.2 Elementary Surface
The elementary surface of u1 = c1 which contains γ2 and γ3 is

dσ1 = ds2 ds3 = h2 h3 du2 du3 . (15)

2.3 Elementary Volume

The elementary volume

dV = ds1 ds2 ds3 = h1 h2 h3 du1 du2 du3 (16)

3 Differential Operators
3.1 Gradient
The gradient is defined by
df = ∇f · dx. (17)
We can also express df as
∂f
df = duj . (18)
∂uj
Using (9), we get
Ej ∂f
∇f = . (19)
hj ∂uj
Or
Ej ∂
∇= . (20)
hj ∂uj

3.1.1 Useful Results

1.
Ej
∇uj = (21)
hj

2. Equation(21) implies that

Ej
∇× =0 (22)
hj
Since
∇ × (f a) ≡ f ∇ × a + ∇f × a, (23)
then,
Ej 1 1
∇× ≡ ∇ × Ej + ∇ × Ej , (24)
hj hj hj

3
 we deduce
∇hj × Ej
∇ × Ej = (25)
hj

3.2 Divergence
Note that
E1 E2 E3
= × .
h2 h3 h2 h3
Using (21),
E1
= ∇u2 × ∇u3 .
h2 h3
Taking the divergence of both sides and noting that

∇ · A × B =≡ B · ∇A − A · ∇B,

we arrive at
E1
∇· = 0. (26)
h2 h3
Or
Ei
∇· = 0, (27)
hj hk
where i 6= j 6= k.

∇ · F = ∇ · (Fi Ei )
Ei
= ∇·( (hj hk Fi ))
hj hk
Ei
= · ∇(hj hk Fi ).
hj hk
(28)

This gives the expression for the divergence

1 ∂
∇·F= (hj hk Fi ), (29)
h1 h2 h3 ∂ui
where i 6= j 6= k.

4
3.3 curl
Using (25, 23), it is readily shown that
1
∇ × F k Ek = ∇(hk Fk ) × Ek , (30)
hk
or

1 ∂(hk Fk )
∇ × F k Ek = (Ej × Ek ). (31)
hj hk ∂uj
Noting that Ej × Ek = ijk Ei , where the permutation symbol ijk = 1 for i, j, k in order
but i 6= j 6= k, ijk = −1 for i, j, k not in order but i 6= j 6= k, and ijk = 0 when two
indices are equal, we obtain,
hi Ei ∂
∇ × F = ijk (hk Fk ). (32)
h1 h2 h3 ∂uj
The expression (32) for the curl can be cast in the familiar form,
 

 h1 E1 h2 E2 h3 E3 
 
1 
 
∂ ∂ ∂ 
∇×F= ∂u1 ∂u2 ∂u3 
 (33)
h1 h2 h3  
 
 
 h1 F1 h2 F2 h3 F3