Vector and Tensor Analysis
Vector and Tensor Analysis
ANALYSIS
by
G.E.HAY
LIBRARY
Published in Canada by General Publishing Company, Ltd., 30 Lesmill Road, Don Mills, Toronto,
Ontario.
Published in the United Kingdom by Constable
and Company, Ltd., 10 Orange Street, London
WC 2.
TABLE OF CONTENTS
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
.. . . . . . . . .
Addition of vectors. . . . . . .
Multiplication of a vector by a scalar
Subtraction of vectors. . . . . .
Linear functions . . . . . . . .
Rectangular cartesian coordinates
The scalar product . . . .
The vector product
Multiple products of vectors
Moment of a vector about a point
Moment of a vector about a directed line
Differentiation with respect to a scalar variable.
Integration with respect to a scalar variable
Linear vector differential equations
Problems . . . . . . . . . . . . . . .
Page
1
2
4'
6
6
7
10
11
15
18
20
22
25
26
28
15. Introduction. . . . . . . . .
16. Some theorms of plane geometry
37
17.
18.
19.
20.
38
39
40
41
43
46
v
Page
Differential Geometry
Introduction. . . . . . .
The principal triad. . . .
The Serret-Frenet formulas
Curvature and torsion
Problems . . . . . . . .
26.
27.
28.
29.
49
50
51
52
53
55
58
Motion of a Particle
30. Kinematics of a particle. . . . . . . . . . . . . . . .
31. Newton's laws . . . . . . . . . . . . . . . . . . . .
32. Motion of a particle acted upon by a force which is a given
function of the time. . .
33. Simple harmonic motion .
34. Central orbits . . . . . .
35.
36.
37.
38.
39.
40.
41.
42.
62
66
68
69
70
72
73
77
79
80
83
84
94
97
43.
44.
45.
46.
vi
102
102
105
107
Page
47.
48.
49.
50.
CHAPTER
V.
61.
62.
63.
64.
65.
66.
67.
68.
69.
70.
71.
124
127
INTEGRATION
III
117
120
VI.
130
134
138
139
140
143
145
146
149
151
152
154
TENSOR ANALYSIS
Introduction. . . . . . .
Transformation of coordinates . .
Contravariant vectors and tensors
Covariant vectors and tensors . .
Mixed tensors. Invariants . . . .
Addition and multiplication of tensors.
Some properties of tensors.
Tests for tensor character
The metric tensor . . . .
The conjugate tensor . . .
Lowering and raising of suffixes
157
157
159
160
161
162
163
164
166
167
169
Vll
Page
72.
73.
74.
75.
76.
77.
78.
79.
80.
81.
82.
viii
170
170
173
174
177
177
178
180
181
184
186
189
CHAPTER
ELEMENTARY OPERATIONS
1. Definitions. Quantities which have magnitude only are called
scalars. The following are examples: mass, distance, area, volume. A
scalar can be represented by a number with an associated sign, which
indicates its magnitude to some convenient scale.
There are quantities which have not only magnitude but also direction. The following are examples: force, displacement of a point, velocity
of a point, acceleration of a point. Su<;h quantities are called vectors if
they obey a certain law of addition set forth in 2 below. A vector can
be represented by an arrow. The direction of the arrow indicates the
direction of the vector, and the length of the arrow indicates the
magnitude of the vector to some convenient scale.
Let us consider a vector represented by an arrow running from a
point P to a point Q, as shown in Figure 1. The straight line through
P and Q)s called the line of action of the vector, the point P is called the
origin of the vector, and the point Q.is called the terminus of the vector.
To denote a vector we write the letter indicating its origin followed
by the letter indicating its terminus, and place a bar over the two letters. The vector represented in Figure 1 is then represented by the
symbols PQ. In this book the superimposed bar will not be used in
any capacity other than the above, and hence its presence can always
Q
Figure 1
2. Addition of vectors. In I it was stated that vectors are quantities with magnitude and direction, and which obey a certain law of
addition. This law, which is called the law ofvector addition, is as follows.
R
Figure 2
a+b = b+a.
Proof.
(2.1)
a+b = PRo
We now construct a vector equal to b, with its origin at P. Its terminus falls at a point S. A vector equal to a is then constructed with
its origin at S. The terminus of this vector will fall at R, and Figure
3 results. Hence
R
Figure 3
b+a= PRo
(2.2)
From (2.1) and (2.2) it follows that a+b = b+a.
Theorem 2. Vectors satisfy the associative law of addition; that is,
(a+b) +c = a+(b+c).
Proof. Let us construct the polygon in Figure 4 having the vectors
a, b, c as consecutive sides. The corners of this polygon are labelled P,
Q, Rand S. It then appears that
(a+b) +c =-PR+c
PS,
a+(b+c) = a+QS
= PS.
Hence the theorem is true.
According to Theorem 2 the sum of three vectors a, b, and c is
p~,
Figure 4
//
Figure 5
= ma+na,
(3.1)
(m+n)a
(3.2)
m(a+b) = ma+mb.
Figure 6
ma
Figure 7
(3.3)
m(a+b)
mPR,
ma+mb
SUo
The two triangles PQ,R and STU are similar. Corresponding sides
are then proportional, the constant of proportionality being m. Thus
(3.4)
mPR=SU.
rna
Figure 8
a+(-b),
Figure 9
(5.1)
Figure 10
c =PQ+QR.
QR
= nb.
6. Rectangular cartesian coordinates. In much ofthe theory and application of vectors it is convenient to introduce a set of rectangular
cartesian coordinates. We shall not denote these by the usual symbols
X,Y and .e, however, but shall use instead the symbols Xu X 2 and X3
These coordinates are said to have "right-handed orientation" or to be
"right-handed" if when the thumb of the right hand is made to point
in the direction of the positive X3 axis, the fingers point in the direction
of the 90 rotation which carries the positive Xl axis into coincidence
0
T--
J--,~l'-
a,
.L -
-!E--------..-Y
/Q
/ R
",
Figure 11
parallelepiped whose edges have lengths lall, la 21 and laal. Hence the
magnitude a of the vector a is given by the relation
(6.1)
= Va12+a22+aa2.
a+b = alil+a2ia+a~ia+blil+b~+baia.
Now the sum of a number of vectors is independent of the order in
which the vectors are added, by Theorem 1 of 2. Hence we may
write the above equation in the form
a+b = lZtil+blil+a2ia+b2ia+aaia+baia.
By the theorem in 3 we may then write this in the form
a+ b = (al +bl)il + (~+b2)ia+ (a3+lh3)ia.
Hence the components of a+b are ~+bl' a2+b2 and aa+ba. This
proves the theorem when two vectors are added. The proof is similar
when more than two vectors are added.
9
Figure 12
The scalar ab cos ex. arises quite frequently, and hence it is convenient
to give it a name. It is called the scalar product of a and b. It is also
denoted by the symbols a b, and hence we have
a b = ab cos ex..
(7.1 )
al
bl
a b
+~
a
b2 + aa ha
b
a b
a' b = albl+a2b2+aaba.
ab =ba.
Proof.
a b = ~bl +a2b2+aab3,
ba = blal+b2a2+bsaa.
10
Since albl = blal , etc., the truth of the theorem follows immediately.
Theorem 2. The scalar product is distributive', that is ,
a (b+c)
a.b+a.c.
ab =0.
However; if it is given that a b = 0, it does not necessarily follow
that a is perpendicular to b. It can be said only that at least one of the
following must be true: a = 0; b = 0; a is perpendicular to b. Similarly, if it is given that
ab = ac,
it does not necessarily follow that b = c. For this relation can be
written in the form a (b - c) = 0, and hence it can be said only that
at least one of the following is true: a = 0; b = c; a is perpendicular
to the vector b-c.
We note the following expressions, in which a is any vector and
iI, ~ and ia are the unit vectors introduced in 6:
(7.3)
i l i l = 1,
i 2 il =0,
i 3 i l = 0,
a.a
il ~
~.~
ia ~
= a2,
= 0,
= 1,
= 0,
ia = 0,
ia = 0,
iaia=l.
il
~.
(i)
(ii)
Figure 13
laxbl
Figure 11:
12
+ +
Cl
and
C2
in terms of Ca , it is found
C2
a2ba-aab2 = aabl-~ba
Cj = K(a 2ba-aab2),
C2 = K(aabcalba),
Ca = K(a l b2-a2bj ).
Now c2 = C12+C22+C32. Hence
c2 = K2[(a 2ba-aab2)2+ (aa bl-al ba)2+ (a l b2-a2bl )2]
= K2[ a12(b22+ba2) +a22(b32+b12) +aa2 ( b12 + b22)
-2 (a2b2aaba+ a3baal bl + al bl a2b2)] .
(8.3)
The first term inside the square brackets can be written in the form
a12( b2-b 12). If the second and third terms are treated similarly it is
found that
C2 = K2[( a12+a22+aa2) b2- (alb l +a2b2+aaba)2]
= K2[a 2b2-(ab cos 0:)2]
= K 2a2b2(1-cos 2 0:)
= K 2a2 b2 sin 2 0:.
But by condition (iii) above, c2 = a2b2 sin 2 IX. Thus K = 1. If
these two values of K are inserted in (8.3) two vectors c result with
the same magnitude but pointing in opposite directions. Only one of
these vectors satisfies condition (ii) above. Now both values of K are
numerical, and are hence independent of a and h. Thus the same
value of K will satisfy condition (ii) for all vectors a and h. Hence it
13
CI
C3
~b2-~bl'
~b3 ,
axb =
Theorem 2.
is
it2
as
bs
b2
axb = -bxa.
(8.7)
Proof.
il
al
bl
bXa
= bl
itl
is
b2
a2
b3
a3
(8.8)
Proo
aX (b+e) = axb+axe.
Let us write
d
14
= aX (b+e), e = axb,
f = aXe.
Then
dl =
(8.9)
iIxil = 0,
i 2 xil = -ia,
ia xiI = i 2 ,
i l Xi2
= is,
~x~
= 0,
ia xi 2
= -iI,
iIxia = -i2'
i2 xia = iI,
iaxia=O.
axb = axe,
a =
b = c;
a is parallel to b-e.
a (b X e) =
b2cl ),
or
(9.1)
a (b X e) =
a1
d2
aa
bi
C1
b2
C2
b3
Ca
15
= b (c X a) = c . (a X b)
=
-c (b X a)
-a (c X b)
-b (a Xc).
(9.2)
V=
la. (bXc) I,
Figure 15
I ad cos (J I
lad I
= I a (bxc) I.
=
The expression
aX (bxc)
is a vector, and is called a veclor triple product of a, band c. Let us write
d
b X c,
a X d.
Hence e
(9.3)
aX(bxc), we have
= b(a.c)-c(a.b).
(9.4)
(9.5)
(axb) (cxd).
It is a scalar. If we consider it as the scalar triple product of a X b,
c and d, and subject these three vectors to two permutations, then
according to Theorem I of 9, we have
(axb).(cxd) =c.[dx(axb)].
If the vector triple product on the right-hand side of this equation is
expanded by the identity in (9.3), we obtain
(axb).(cxd) = (ca)(d.b)-(c.b)(d.a),
or in a form more easily recalled,
(9.6)
"
a
o
Figure 16
P
Theorem 1.
xXa.
P =pa,
where p is the perpendicular distance from 0 to the line of action ofa.
Proof. Now P = xa sin 0(, where 0( is the angle between x and a.
But p = x sin 0(. Hence P = pa.
Theorem 2. The moments. about a point of any two equal vectors
with the same line of action are equal.
Proof. Let a and a' be two equal vectors with the same line of ac-
o
Figure 17
19
tion, and with origins X and X', as shown in Figure 17. Let 0 be any
point, and let ox = x, ox' = x'. Let P and P' be the moments of
a and a' about o. Then
P = xXa,
P' =x'Xa'.
a we
11. Moment of a vector about a directed line. Each line defines two
directions which are opposite. A line is said to be directed when one of
these directions is labelled the positive direction and the other the negative direction.
Let us consider a directed line L, and let b denote a unit vector
pointing in the positive direction of the line, as shown in Figure 18. We
Figure 18
20
where <I> is the angle between P and the unit vector b. Hence Q, = b. P.
But P = x X a, where x = OR. Thus
(ILl)
Q,=b.(xxa).
--~------~~------~~------------L
o
0'
Figure 19
let P' be the moment of a about 0', and let Q; be the corresponding
moment of a about L. Then
Q, = b (x X a),
Q;
= b (x' X a) ,
Q; = b[(O"O+x) xa]
=b(O"Oxa)+Q,.
Since band 0"0 have the same line of action L, then b (0"0 X x) = 0
by Theorem 2 of 9. Thus Q; = Q" which proves the theorem.
Theorem 2. IfP denotes the moment ofa vector a about the origin
of the coordinates, then the three components of P are equal respectively to the moments of a about the three coordinate axes.
Proof. The truth of this theorem follows immediately from the
above definitions of the moments of a vector about a point and about
a line.
21
o
Figure 20
s~en
. " ~a
1ImltLl.u .... 0 ~ U
is the rate of change of a with respect to u. It is also called the derivative of a with respect to u, and is denoted by the symbol
da
du
22
1"Imlt
. ~a
.
u
Ll.u -+ 0 ~
~:,
so that
~:.
da dm
m du + du a,
(12.2)
d
du (ma)
(12.3)
d
db da
d)a. b) = a du + du' b
(12.4)
Proof of (12.1).
the sum a+b is
( 12.5)
d
du (a X b)
db da
du + du
X b.
According to Theorem 2 of 2 ):he sum of a number of vectors is independent of the order of summ~tion. Thus the right side of (12.5) can
be written in the form a - a+b - b+Ila+/lb, which reduces to
/la+/lb. Thus
/l(a+b) = /la+/lb.
Ifboth sides of this equation are divided by /lu, and if /lu is then made
to approach zero, (12.1) is obtained.
Proof of (12.2).
ma is
(12.6)
and (3.2). Because of the law exemplified by (3.1) we can then write
(12.6) in the form
(12.7)
and because of the law exemplified by (3.2), we can then write (12.7)
in the form
Ll(ma) = ma+mLla+Llm a+Llm Lla - ma
= m Lla+Llm a+Llm Lla.
If both sides of this equation are divided by Llu, and if Llu is then made
to approach zero, (12.2) results.
Proof of (12.3). When u increases by an amount Ll u, the change in
ab is
Ll(a. b) = (a+Lla). (b+Llb) - a b.
Since the scalar product is distributive, this equation may be written
in the form
Ll(a.b)
a b+aLlb+Lla b+LlaLlb - a b
If both sides of this equation are divided by the scalar .6.u, we have
alil+a2i2+a3i3'
d ( .)
du aliI
+dud
( .)
a 212 +
d ( .)
du aal 3
_
q~+ di2+ dia+da l
- al du a2 du a3 du du
24
II
+da2 +daa
du 12 du 13
Now aI' a2 and a3 are scalar functions of u. Also it> i2 and i3 are unit
vectors pointing in the directions of the positive coordinate axis. If they
are the same for all values of u, then
dil = di2 _ dia _ 0
du
du-du-'
and so
da _ dal +da2 +daa
du-du
11
dU
12
dU
13
25
Proof.
We have
a(u) = pf(u) +q g(u) + ... ,
are given
pk+ql+. .
The equation
(14.1 )
in which a and Po, PI, "', Pn are given functions of the scalar u and x
is an unknown vector, is a linear vector differential equation of order'n.
Let F denote the differential operator in (14.1). Then (14.1) can be
expressed in the form
( 14.2)
F[x] = a.
26
... , Cn being arbitrary constant vectors and)'1' Y2, Ya, ... , Yn being
n linearly independent solutions of the homogeneous scalar differential
equation F[y] = 0.
Proof. Let us introduce the unit vectors iI' i2 and i3 with directions
independent of u. Then
C;, C2, C a,
lit
F[x]
Q1i1+a2i2+aai3'
xliI + x2i 2+ x3ia ,
= F[X1]i1+F[X2]i2+F[X3]ia.
the vectors Ct, C2' C 3 , .. , C n being arbitrary constant vectors. Equation (14.4) gives the general solution of Equation (14.1). We note that
Y is the general solution of the homogeneous equation F[ x] = 0, and
that A is a particular solution of Equation (14.1). The particular solution A can be found by procedures very similar to those used to find
particular solutions of linear scalar differential equations. This is demonstrated below.
27
(14.6)
d y _ dy
du 2
du
_~ =
o.
m2 -m-2 =
o.
It has roots -1,2, whence the required solutions of (14.6) are eOu and
e2u . Thus
Y = cleou+c2e2u.
We now use the method of undetermined coefficients to find a particular solution A of Equation (14.5). The function on the right side of
(14.5), and the derivatives of this function, are linear functions of sin u,
cos u, u, 1, none of which are particular solutions of (14.6). Hence, we
look for a particular solution A in the form
A = b sin u+c cos u+du+e,
c =p,
d =-q,
e =0.
Problems
1. The vectors a, b, c and d all lie in a horizontal plane. Their
28
2a- b = n,
a
Ca
ga
b
Cb
g.b
Cc
gc
n.(bxc)
n.(cxa)
n.(axb)
a+---b+
c.
a.(bxc)
b(cxa)
c.(axb)
This formula can be us~d to express any vector n as a linear function of
any three vectors a, band c not lying in the same plane.
26. Express the vector n = i 1 +2i 2+3ia as a linear function of the
vectors a, band c defined in Problem 13. (See Problem 25.)
27. Express the vector n = 2il - 2i2 - 3ia as a linear function of the
vectors a, band c defined in Problem 13. (See Problem 25.)
28. Show that
ax[(fxb) X (gxc)]+b)<[(fXc) X (gxa)]
+cx[(fXa) X (gxb)] = O.
29. If 0 is the origin of the coordinates and A, Band C are three
points such that
OA = 2i1 +2i 2-ia ,
AB = i1 - i 2+2i a ,
BC = -2il +2i2 - 3ia ,
find (i) the moment of BC about A, (ii) the moment of CB about 0,
(iii) the moment of BC about the directed line OA, (iv) the moments
of BC about the coordinate axes.
30. If a and b are two vectors, prove that a times the moment of b
about the line of action of a is equal to b times the moment of a about
the line of action ofb.
3l. If a(u) has a constant magnitude, show that
da
a-=O.
du
32. Ifa = p cos u+q sin u, where p and q are constant vectors and
u is a variable, show that
n=
da). d3~
.
du
31
(i)
(ii)
a = i l cos
... )
(III
a =
.
2il
4-u2
+i2 sin u - i3 e
2i2 .
4+u 2
(i)
(ii)
(iii)
f (p+qu2)du,
f (p cos u +q sec2u)du,
d
f p+qu
4 -u
U.
36. Find the vector x(u) in each of the following cases, given that
p, q and r are constant vectors:
dx
(i)
- = p u2 +q e2u
(ii)
du
d2 x
du 2
2
(iii)
d x
dw
'
(p
SIn
U,
(iii)
(iv)
d2 x
dx
du 2 - 2 du
d'x _
du'
+ 5x = 0,
6 d x + 11 d x _ 6 dx =
2
3
du 3
du
O.
du
32
(i)
dx
du - 3x = p(3u2 + 1),
(ii)
d2 x
du 2 -4x
(iii)
dx
du2
2 dx
du
39. Find the general solutions of the following differential equations, given that p and q are constant vectors:
d2x
dx
(i)
du 2 - 2 du + x = 2pe",
(ii)
(iii)
2
x
dx
u
- 2 - 3 - = 2pe
du
d2 x
dx
u2 - _
U - - 3x
du 2
du
au
l8qu2
6p.
'
33
CHAPTER
II
APPLICATIONS TO GEOMETRY
34
c
~--------~------~B
o
Figure 21
d =ub.
AD =vAC.
Hence we can write (16.2) in the form
(16.3)
= a+vAC.
We now equate the above two expressions given for d in (16.1) and
(16.3), obtaining
(16.4)
a+vAC = u b.
The next step is to express all vectors in this equation as linear functions
of any two vectors in the plane, say a and c. From the figure we see
that
b = a+c,
AC =-a+c,
whence (16.4) becomes
a+v (-a+c)
= u (a+c) ,
(l-u-v) a
= (u-v)
or
c.
35
Since a and c do not have the same line of action it then follows that
1 - u - v = 0,
u- v = 0.
We solve these equations foru and v, obtaining u = v = l. Thus Equation (16.1) becomes d = i b, from which it follows that D is the middle
point of OB.
We have now proved that the point of intersection D of the diagonals is the middle point of one of these diagonals. From symmetry, D
must also be the middle point of the other diagonal.
Theorem 2. The medians of a triangle meet in a single point which
trisects each of them.
Proof. Let us consider the triangle OAB in Figure 22. The points
B
O~-----~------">IA
E
Figure 22
C, D and E are the middle points of the sides, and F is the point of
intersection of the medians AD and BE. We must prove that F is a
point of trisection of each of the three medians AD, BE and OC.
For convenience we denote the vectors drawn from a to the points
A, B, C, D,E and F by a, b, c, d, e and f, respectively. Now F lies on
the median AD. Thus there exists a scalar u such that
(16.5)
DF=uDA
OD = t b,
DA =-t b+a,
OE
EB
= ta,
=-ta+b.
(t-u-tv)a = (t-tu-v)b.
Since a and b do not have the same line of action, it follows that
t - u - tv = 0,
t-tu-v =0.
37
~-------.------------.~
x.
Figure 23
38
Figure 24
(18.l )
AC
CB
-=-.
m
n
If we now denote the position-vectors of A, Band C by a, band c,
respectively, then
AC = c - a,
CB = b - c,
and (18.1) can then be written in the form
n(c-a) = m(b-c).
Solving this equation for c, we obtain
mb+na
m+n
This formula expresses the position-vector c of the desired point C in
terms of the known quantities a, b, m and n.
In books on Analytic Geometry, formulas are usually given which
express the coordinates of C in terms of m, n and the coordinates of
A and B. It should be noted that (18.2) is entirely equivalent to these
formulas, for these formulas can be deduced from (18.2) simply by
equating the components of the left side of (18.2) to the components of
the right side of(18.2).
( 18.2)
c=;=
19. The distance between two points. Let us suppose that A and Bare
two given points, and that it is desired to find the distance d between
A and B in terms of the position-vectors a and b of A and B. Figure
25 illustrates the problem. Now
39
d= IABI.
But AB = b - a.
d2
Thus
= IABI2 =AB'AB = (b-a).(b-a)
whence
= V (b -
a) . (b - a) .
B
Figim 25
20. The area oj a triangle. Let us suppose that A, Band C are three
given points, and that it is desired to find the area 6.ahc of the triangle
ABC in terms of the position-vectors a, band c of A, Band C. Figure
26 illustraties the problem.
B
Figure 26
40
6.ahc = t<P,
where
(20.2)
Now AB
cp
b-a, AC
ABXAC.
c-a. Thus
cp
(b-a)x(c-a).
This simplifies to
(20.3)
cp
bxc+cxa+axb.
The required area of the triangle is thus given by (20.1), qJ being the
magnitude of the vector given by (20.3).
A property of the vector cp will now be recorded, for future use.
Since this vector is equal to AB X AC we conclude that the vector cp is
perpendicular to the plane of the triangle ABC, and its direction is that indicated
by the thumb of the right hand when the fingers are placed to indicate the direction
of the passage around the triangle from A to B to C.
21. The equation of a plane. There are several ways in which a plane
can be specified. For example, three points which are on the plane and
do not lie on a single straight line can be given, or a line in the plane
and a point on the plane but not on the line can be given. In each of
several such cases we sh~ll now deduce the equation which must be
satisfied by the positiori-vector x of every point X on the plane. This
equation will be referred to simply as the equation of the plane. In books
on analytic geometry the equation of a plane usually appears as an
equation which involves scalars only, and is satisfied only by the coordinates of points on the plane. We shall refer to this latter equation
as the cartesian form of the equation of the plane.
(i) To find the equation of the plane through a given point and perpendicular
to a given vector. Let A be the given point and b be the given vector.
Figure 27 illustrates the problem, the plane P being the plane in
question.
Let X be a general point on P, and let a and x denote the positionvectors of A and X, respectively. Now AXis perpendicular to b. Thus
AXb =0.
41
(x-a)b
(21.1 )
=0.
o~~
_____________
x,
x,
Figure 27
42
cp = bxc+cxa+axb
is perpendicular to the plane P. Hence, by Problem (i) above the
equation of P is
(21.2)
(x-a).cp
= O.
(21.3)
xcp = a.(b
c).
~:o-
_______
%.
Figure 28
22. The vector-perpendicular from a point to a plane. The vector-perpendicular from a point D to a plane P is the vector with origin at D and
terminus at the point on P nearest D.
(i) To find the vector-perpendicular from a point D to a plane P through a
given point and perpendicular to a given vector. Let A be the given point
and let b be the given vector. Figure 29 illustrates the problem. We
denote the position-vectors of the points A and [) by a and d, respectively. If the point E is the foot of the perpendicular from the point D
to the plane P, then DE is the vector-perpendicular from the point D to
the plane P. We shall denote it by the symbol p. It is required to find
p in terms of a and b.
43
Figure 29
(22.1 )
From the figure we see that
or
(22.2)
d+Kb+EA - a
O.
(d-a).b+Kb2 =
Thus
K=
44
(a-d)b
b2
o.
(22.4)
(a-d).bb
b2
(22.5)
(ii) To find the vector-perpendicular from a point D to a plane P through three
given points. Let A, Band C be the three given points, with positionvectors a, band c, respectively. Figure 30 illustrates the problem. It
is desired to find the vector-perpendicular p in terms of a, b, c and d.
~-------::Oi.
B
O~~~
_______I
f)
Figure 30
(j'
=bxc+cxa+axb
45
(a-d) cp
p=---cp.
cp2
a (b X c) - d cp
cp2
cpo
23. The equation of a line. There are several ways in which a line in
space can be specified. For example, two points on the line can be
given, or two planes through the line can be given. In each of several
such cases we shall now deduce the equation which must be satisfied
by the position-vector lit of every point X on the line. This equation will
be referred to simply as the equation of the line.
(i) To find the equation of the line through a given point and parallel to a
given vector. Let A be the given point, with position-vector a, and let b
be the given vector. Figure 31 illustrates the problem, L being the line
in question.
o~'~
/
x,
_______________
Figure 31
lit
AX X b = O.
46
But AX
(23.1 )
(x-a)Xb =0.
0/E~
x,
_____________
Figure 32
(x - a) X (b - a)
O.
(iii) To find the equation of the line through a given point and perpendicular
to two given vectors. Let A be the given point with position-vector a, and
let band c be the given vectors. Figure 33 illustrates the problem, L
being the line in question. Now L is parallel to the vector b X c. Hence,
by Problem (i) above, the desired equation of Lis
(23.3)
(x-a)X(bXc) =0.
(iv) To find the equation of the line through a given point and perpendicular
to the plane through three given points. Let A be the given point on the line,
and let B, C and D be the given points on the plane. Figure 34 illustrates
47
o
Figure 33
the problem, Land P being the line and plane in question. We denote
the position-vectors of A, B, C and D in the usual manner. Let us consider the vector (j) given by the relation
(j) = cXd+dXb+bXc.
Figure 34
(23.4)
(x-a)xcp=O.
But
ex =
x-c.
Thus
(24.1 )
(x-c)(x-c) =a2
o
We shall now prove the following well-known property of a sphere:
the angle at the surface of a sphere subtended by a diameter of the
sphere is a right angle. For convenience, the origin of the coordinate
system is chosen at the center of the sphere, as shown in Figure 36. Let
D .and E be points at the ends of a diameter, and let X be a general
point on the sphere. We denote the position-vectors of these points in
the usual manner. From the figure
DX
x-d,
EX
x-e.
49
Thus
DXEX= (x-d)(x-e) =x.x-x(d+e)+d.e.
But d+e=O, and if a denotes the radius of the sphere then xx = ai,
de = _ai. Hence it follows that
DXEX =0,
j)
Figure 36
o
Figure 37
50
or
y-c =CX+XY.
But CX
CX(y-c) =a2
(x - c) .(y - c)
= a2
Differential Geometry
26. Introduction. "'Ie shall consider here only a small portion of the
differential geometry of curves in space. Rectangular cartesian co-ordinates Xu X 2 and Xs are introduced, with origin at a point O. The quantities Xl' x2 and Xs denote the coordinates of a general point X with
position-vector x. If iI' i2 and is are unit vectors in the directions of
the positive coordinate axes, then as before,
(26.1)
= xlil+X2i2+X3i3'
positive for points on one side of A, and negative for points on the
other side of A. The equation of the curve may then take the form
x = x(s).
The derivatives with respect to s of the function x(s) will be denoted by
x', x", x"', etc.
o
Figure 38
way by the symbols j/, j2 I and j3'. They are shown in Figure 38, and
are defined by the conditions:
(i) jl is tangent to the curve C, and points in the direction of s
increasing;
(ii) j2lies in the plane of the vectors jl and jl', and makes an acute
angle with jl' ;
(iii) js is such that the vectors jl' j2 and j3 form a right-handed triad 1.
1 At points on the curve where j/ is equal to zero, these conditions are not sufficient,for a unique determination of j. and j.. We exclude such points from consideration here.
52
The straight line through the point X and parallel to j2 is called the
principal normal to the curve. The straight line through X and parallel
to ja is called the binormal to the curve. The vectors jl, j2 and ja are
called the unit tangent vector, unit normal vector, and unit binormal vector,
respectively. The triad formed by these vectors is called the principal
triad. The plane through X and perpendicular to jl is called the normal
plane. The plane through X and perpendicular to ja is called the
osculating plane.
limit
!J.s ..... 0
ilx
ils
Now
limit lilxl
!J.s ..... 0
ils
1.
Thus x' is a unit vector. Further, the vector ilxjils lies along XY, and
its direction then tends to that of jl as ils tends to zero. Since jl is a unit
vector, we can then write
(28.1 )
The vectors jl" j2' and ja' can each be expressed as ~ linear function of any three non-coplaner vectors. In particular, they can be expressed as linear functions of the vectors jl, j2 and ja, and we then
have relations of the form
(28.2)
53
(28.3)
jljl=l,
j2-j1 = 0,
ja -jl = 0,
jd2= 0,
jd2= 1,
ja-j2=0,
jd3=0,
j2-ja=0,
ja-ja= 1.
(28.2), and then make use of Equations (28.3), we find that a12 +a 21 = 0.
Similarly we can find two similar relations, and we have altogether
(28.5)
. a I2 +a 21 =0,
a23 +aa2=0,
a31 +a I3 =0.
So far, only conditions (i) and (iii) above have been used. By condition (ii) the vector-jl' is to be in the plane of jl and j2. This can be
true only if
(28.6)
ala =
By condition (ii), the vector jl' is to make an acute angle with j2. If
this angle is denoted by oc, then cos oc must be positive. But
Ijl'l cos oc=jdl'
If we substitute here for jl' frrun Equations (28.2) and then use Equa-
54
Thus
(28.7)
> 0,
and because of Equations (28.4), (28.5), (28.6) and (28.8), we can now
express Equations (28.2) in the form
j1' = Xj2,
j2' = "ti3 - Xj1'
js' = - "ti2'
These are the Serret-Frenet formulas. They were given originally by
Serret (1851) and Frenet (1852) in an equivalent form which did not
involve vectors. The quantities x and ", which are functions of the arc
length s of the curve C, will be considered in some detail in the next
section.
(28.10)
. M
= 11m1t-,
/1, -+ 0
.:ls
where .:l6 is the angle between the tangents to the curve C at the points
X and Y in Figure 38. Thus x is the rate at which the tangent at the
pomt X rotates as X moves along the curve. The reciprocal of x is
called the radiuJ of curvature, and will be denoted by the symbol p.
The quantity" appearing in Equations (28.10) is called the torsiop.
of the curve. It can be shown that
" =
. M)
11m1t-,
/1t -+ 0
.:ls
where .:leI> is the angle between the binormals to the curve C at the
55
points X and Y in Figure 38. Thus l' is the rate at which the unit binormal at the point X rotates as X moves along the curve. The reciprocal of l' is called the radius of torsion, and will be denoted by cr.
To find x, we note from Equation (28.1) that
jl' = x".
jl = x',
(29.1)
Substitution from the second of these relations for jl' in the first of the
Serret-Frenet formulas then yields the equation
x" = xj2.
(29.2)
x" . ,,".
x. =
l'
vi,," .,,".
We now substitute for j2' from the second Serret-Frenet formula, obtaining
. . )+'.
(29.4)
x ", =x (1'Ja-XJl
XJ2.
From Equations (29.1), (29.2) and (29.4) it now follows that
x' (x" xx"') = jl[Xj2X (x-rja - X2jl+Xj2)]
= jl . [x 21' jl + x3jaJ
= X 21'.
Thus
(29.5)
't'
,,'" =
56
xl'''il-t-X2'lfi2+Xa'''ia
It
= y;-,,2
+ x----:-:2+X
,,2
1
2
3,
(29.7)
Since It can now be found, we can obtain the unit tangent vector jl
and the unit normal vector j2 by use of Equations (29.1) and (29.2).
The unit binormal vector j3 can then be found easily, since it is equal
to j 1 X j 2. We have the collected results
(29.8)
Jl
= X,
J2
=;;:1 X
II
J3
= -It x' X
x" .
Let us now find the equation of the tangent to the curve at the point
(y- x) Xjl = O.
Because of Equation (29.8), this can be written in the form
(29.9)
(y- x) Xx'
O.
In the same way, the equations of the principal normal and. binormal
can be found in the forms
(29.10)
(29.11)
(y - x) . jl
= o.
(y-x)x'=O.
In the same way we can find the equation of the osculating plane in the
form
(y - x) . (x' X x") = o.
(29.13)
57
Problems
1. Prove that the line joining the middle points of any two sides of a
triangle is parallel to the third side, and is equal in length to one half
the length of the third side .
.2. Prove that the lines joining the middle points of the sides of a
quadrilateral form a parallelogram.
3. If 0 is a point in space, ABC is a triangle, and D, E and Fare
the middle points of the sides, prove that
OA+OB+OC = OD+OE+OF.
= (l+m+n)
OD,
58
and AB are denoted by D, E and F respectively, G is the point of intersection of the perpendiculars from the vertices to the opposite sides,
and H is the center of the circumscribed circle, prove that
iG
= 2 HD,
BG = 2 HE,
GG = 2 HF.
In problems 13-22 the points A, B, C and D have the following coordinates: A(-l, 2, 3), B(2, 5, -3), C(4, 1, -1), D(l, 3, -3).
13. Find the position-vectors of the points of trisection of the line
segment AB.
14. Find the distance between the points A and B.
15. Find the area of the triangle ABC.
16. Find the cartesian form of the equation of the plane through A
and perpendicular to OB.
17. Find the cartesian form of the equation of the plane through
(i) the origin and the points A and B, (ii) the points A, Band C.
18. Find the distance from the point D to the plane through A and
perpendicular to OB.
19. Find the distance from D to the plane ABC.
20. Find the cartesian form of the equation of the line through A
and parallel to BG.
21. Find the cartesian form of the equation of the line through
A andB.
22. Find the cartesian form of the equation of the line through
D and (i) perpendicular to the plane through A, Band C, (ii) perpendicular to BG and OG.
23. A plane passes through a given point A with position-vector a,
and is parallel to each of two given vectors band c. Derive the equation of this plane in the form
(x-a).(bXc) =0.
24. A straight line L passes through a point A with position-vector
a, and is parallel to a vector b. A vector p has its origin at a point C
59
=a-
(a-e)h
e - --b2-- h .
,. +("
-2)"J2 +(2 x't'
'+ X't' ')"J3'
x "" =-3XXJl
X -x3 -x,
28. If the position-vector x of a general point on a curve is given as
a function of a parameter t, and if primes denote differentiations with
respect to t prove that
1
S"2
x'
JI = s' '
Also, derive the equations of the tangent, principal normal, binormal,
normal plane and osculating plane in the following forms:
60
tangent,
(y-x)xx' =0;
principal normal,
(y- x) X (x" -
binormal,
normal plane,
osculating plane,
I;s x')
= 0;
where a and rJ. are constants, and t is a parameter. Find p, (j and the
principal triad. Answer: p = a cosec 2 rJ., (j = 2a cosec 2rJ., jl = sin rJ.
(-il sin t+i2 cos t+ia cot rJ.), j2 = -il cos t- i2 sin t, ja = cos rJ. (i] sin t
- i2 cos t+ia tan rJ.).
30. The position-vector x of a general point on a curve is given by
the relation
x = a(3t - ta)i 1 +3at2i 2+a(3t+t3)ia,
where a is a constant and t is a parameter. Find p, (j and the principal
triad. Answer: p = (1 = 3a yS, v2j] = y-l(rJ. i 1+~i2+yia), j2 = y-l(-~il
+et i 2), v'Z""ja = y-l(_et i] - ~ i 2 +y is), where rJ. = 1 - t2 , ~ = 2t,
Y = 1+t2.
61
CHAPTER
III
~---"'-------'"
o
"
62
Figure 39
Let the curve C in this figure be the path of the particle, and let the
point X denote the position of the particle at time t. The vector OX
is the position-vector of the particle. We denote this vector also by x.
It is a function of the time t.
The velocity v of the particle relative to the frame of reference, and
the acceleration a of the particle relative to the frame of reference,
are defined by the relations
(30.1)
dv d2x
a= dt = dt 2 '
dx
v= dt'
xlil+X2i2+Xsis'
If we now adopt the convention that a single superimposed dot denotes a first time derivative, and a pair of superimposed dots denotes
a second time derivative, then
v =
dx
.
'.
'.
dt
= Xlll+X212+X313,
a =
de =
dv
..
".
..
XIII +X212+X31 3
(30.2)
Xl, 2, X3;
Xl' X2, X3
(ii) The components of the velocity and acceleration in the directions of the
principal triad of the curve traced out by the particle. The curve C in Figure
39 is the path of the particle. Let jl' j2 and j3 denote the principal
triad at the general point X on C. The principal triad was discussed
in 27. If s denotes the arc length of C, then from Equations (28.1)
and (28.10) we have
dx
(30.3)
Jl = ds'
63
dx
v = dt = ds s,
and because of the first equation in (30.~) we then have
(30.4)
v = S jl'
Thus the velocity of the particle is directed along the tangent to its path, and the
speed is v = S.
Because of (30.4) we have
dv
djl
a
= dt =
00.
SJ1+ S
dt'
djl
Tt
djl .
ds s =
o.
XSJ2 ,
and hence
(30.5 )
Thus the acceleration a lies in the osculating plane of C. Also, the components
in the directions of the tangent, normal and binormal are
of a
(30.6)
dv
s=v=vds'
'0
v2
XS2
= xv2 = p ,
0,
64
-"3
~~-----
---- -".
Figure 40
If iI, i2 and i3 are the usual unit vectors associated with the rectangular cartesian coordinate axes in Figure 40, then
dkl = dkl
dt
de
(30o 7)
dk2
dt
(.
6
6) 6 = k 26 ,
-11 sm +12 COS
= dk2
de 6
(.
6
-11 cos
- 12
sin 6) 6 = - k 16 ,
dk3 -0
dt
From Figure 40 it follows that the position-vector x of the particle is
given by the relation x = rk1+x3k30 Thus
dx
= de
:1_
dkl
dk3
Tt .
65
(30.8)
Thus the desired components of v are
(30.9)
From (30.8) we find that
dv..
....
..
a = di = rk1 +(r6+r6)k 2+x}t3
. dk 1 6 dk2 . dk3
+r (ft+r (ft+X3"dt'
Substitution in this equation from (30.7) then yields
(30.10)
a = (r - r( 2)kl+(2rB+r6)k2+X3k 3'
r- ,62,
2;6+r6
=;
;t (r 6),
2
.,
X3 .
= 0,
r = x, and we find from (30.9) and (30.11) that the components ofv
and a in the directions of the parametric lines of the plane polar coordinates x, 6 are
.
. 1 d
.
x, x6', x - x6 2, -x -dt (x 26) .
(30.12)
31. Newton's laws. The concept of force is intuitive. We can define a unit force as that force which produces a standard deflection of
a standard spring. Hence we can assign a numerical value to the magnitude of any force.
We know that forces have magnitude and direction. It has been
verified experimentally to within the limits of experimental error that
forces obey. the law of vector addition. Hence we shall assume that
forces are vectors. The sum of two or more forces is sometimes called
the resultant of the forces.
The term "mass of a body" refers to the quantity of matter present
in the body. We can define a unit mass as that mass which, when sus66
k m a,
67
units. Thus, for example, when a force of one pdl. acts on a particle
with a mass of one lb., the acceleration of the particle is one ft./sec. 2
The systems of units in the second and third columns of Table I are
called foot-pound-second sytems, or simply f.p.s. systems. The sy~
stem of units in the fourth column is called the centimeter-gramsecond system, or simply the c.g.s. system.
e.g.s.
f.p.s.
Unit
Unit
Unit
Unit
of length
of mass
of time
of force
foot (ft.)
pound (lb.)
second (sec.)
poundal (pdl.)
TABLE
foot (ft.)
slug
second (sec.)
pound-weight (Ib.wt.)
centimeter (em.)
gram (gm.)
second (sec.)
dyne
The lb. wt. is the force exerted on a mass of one lb. by the earth's
gravitational field. If G denotes the acceleration due to gravity, expressed in ft./sec. 2, then
.
I lb. wt. = G pdl. ,
1 slug
G lb.
dv
F=m-
dt '
68
where p and q are given constant vectors. Because of the first equation in (31.2) we then have
dv
m dt = 12p+q cos t,
whence
mv =
mx = f(l2pt+qsint+r) dt
where s is an arbitrary constant vector. The arbitrary constant vectors rand s can be found if the initial values of x and v are known.
If these initial values are Xo and vo, it is readily found that
r = mvo ,
or
d2 x
-d2
k
+ m-x
=0.
69
C1
cos
11k
V;; t +
C2
. 11k
sm
V;; t ,
where C 1 and C 2 are arbitrary constant vectors. These arbitrary constant vectors can be found if the initial values of x and v are known.
If these initial values are Xo and vo, it is readily found that
.
C2
V~vo,
whence we have
x = Xo cos
. 1VII; t .
V;; t+vo 1VI;;I sm
11k
F =ma.
Let Xo and Vo be the initial values of x and the velocity v. The entire
path of the particle will be in the plane P containing the vectors
Xo and yo. Let x and 6 be polar coordinates in this plane. The components ofF in the directions of the parametric lines of these coordi70
nates are -F, O. Also, the components of a in these directions are given
in Equation (30.12). Hence we have
(34.1)
-F=m(x-xe2 ),
(34.2)
md
0 = -;
dt (x 26) .
y = 1Jx.
const.
whence
(34.4)
Then
(34.5)
By substitution in (34.1) for x, 6 and
we finally obtain
x from
(34.6)
Now F is a function ofy alone. Once the form of this function has been
assigned, we can find the path of the particle by solving Equation
(34.6).
Let us now consider the special case when F varies inversely as the
square of x. Then we can write
= ymy2,
d2y
d6 2
+y
= h2
'
(34.7)
x1 = hY +
2
C1
Il'
6
cos v+c 2 sm ,
m=:Emj.
(35.1)
j= !
(35.2)
Xc =
~
j=!
mj Xj.
=f
p dV,
=f
p x dV.
(35.4)
mxc
L~
Figure 41
Let p and q denote the lengths of the perpendiculars from the particle
to two perpendicular planes P and Q" as shown in Figure 42. The
product of inertia of the particle with respect to these two planes is
defined to be the scalar K given by the relation
K=mpq.
73
J:
Figure 42
11 =
~ mj (X2j2+X~'3) ,
j = I
(36.1)
12
= ~ mj (X2j3+X2it) ,
13
j= I
~ mj (X2i! +X2i2)
j= I
(36.2)
Kl
=j =~ I mj Xh Xj3
K2
=j =~ I mj Xia Xi!
K3
=. ~
N
J = I
74
mj xiI Xh .
Xa
mj
o"....=-________ X.
x,
Figure 43
1=
~ mjPi2
j = 1
Xj
sin () =
lb X Xjl ,
b 3 x jI -
bi
Xj3,
bi
Xi2 - b 2 XiI'
and hence p/ is equal to the sum of the squares of these three components. Thus
N
1= ~
j
mj [( b 2xj3 - b 3 Xj2)2
+ (b 3 XjI -
bi
Xj3)
+ (b
Xj2 - b 2 XjI)2]
75
N
+
ba2 L
N
mj (Xj12+Xj22) -
2 b2ba L
mj Xj2 Xja
j=l
j= 1
- 2 ba b1 L
j=l
N
mj Xj3 Xjl -
2 b1 b2 L
mj Xjl Xj2
j=l
2 Ka bi b2
It will be noted that bl , b2 and ba, which are the components of the
unit vector b on the line L, are also the direction cosines of L. Equation (36.3) is the desired equation which permits a simple determination of the moment of inertia I of a system about any line L through
the origin, once 11,12 , 1a., KI , K2 and Ka have been found.
If it happens that KI = K2 = Ka = 0, the coordinate axes are said
to be principal axes of inertia at the point O. It can be pr.oved that at
every point there is at least one set of principal axes of inertia. I In
many cases, principal axes of inertia can be deduced readily by considerations of symmetry of the system of particles. For example, at
the center of a rectangular parallelepiped the principal axes of inertia
are parallel to the edges of the body.
We shall now state without proof two theorems the proofs of which
are very simple and may be found in almost any text book on calculus.
The theorem of perpendicular axes. If a system of particles lies entirely
in a plane P, the moment of inertia of the system with respect to a line
L perpendicular to the plane P is equal to the sum of the moments of
inertia of the system with respect to any two perpendicular lines
intersecting L and lying in P.
The theorem of parallel axes. The moment of inertia I of a system of
particles about a line L satisfies the relation
1= 1'+m 12 ,
where l' is the moment of inertia of the system about a line L' parallel
to L and through the center of mass of the system, m is the total mass
1 See J. L. Synge and B. A. Griffith, Principles of Mechanics, McGraw-Hill Book
Co., 1942, pp. 311-321.
76
Figure 44
Let e denote the angle between wand x, and let p denote the length
of the perpendicular from X to the line of action of w. The displacement dx of the point X in time dt has the following properties:
(i)
(ii)
when the fingers are placed to indicate the sense of the rotation
6 from w to x;
(iii) its magnitude is pwdt, which is equal to xwdt sin 6, x being the
magnitude of the vector x.
In view of the definition in 8 of the vector product of two vectors,
it then appears that d x = w X x dt. Thus
dx
Tt=wxx.
If the point 0 is fixed in a frame of reference S, the velocity v of the
point X relative to S is then
v = wxx.
(37.1)
v=u+wXx.
We shall now prove that angular velocity obeys the law of vector
addition. Let us consider a body which is rotating simultaneously
about two lines Land L' which intersect at a point 0 fixed in a frame
of reference S. These angular velocities can be represented by the
arrows wand w' in Figure 45. Let X be a general point in the body,
LEz
/0
'"
Figure 45
= wXx+~'Xx.
To complete the proof we must show that (37.3) can be written in the
form v = w" X x, where w" is an arrow obtained by the application
78
of the law of vector addition to the arrows wand w'. Even though
angular velocity has not been assumed to satisfy the law of vector
addition, Equation (8.5) may be applied to the two products in (37.3)
to yield
VI
W 2 X3 - W3 X 2 +W'2 X3 - W'3 X 2
(W 2 +W'2) X3 -
V2
(w 3 +w'3) X 2 ,
v=w"XX,
where w"is an arrow having components wl+wl', w2+w'2, W 3 +W'3.
But these are the components of the vector obtained by applicai:ion
of the law of vector addition to the arrows wand w'. Hence w" is
equal to the vector sum of wand w', and so angular velocity is a
vector.
It will be noted that two angular velocities can be added only when
their lines of action have a point of intersectiol1, and that the line of
action of the sum passes through this point of intersection.
38. The time derivative of a vector. Let us consider a set of rectangular cartesian coordinate axes with origin 0 fixed in a frame of reference
S, and with axes rotating relative to S with angular velocity w. Then
the line of action of w passes though o.
If iI, i2 and i3 are the usual unit vectors associated with these coordinate axes, then the velocities relative to S of the terminuses of these
vectors are
WXi2,
W X i3
Hence
(38.1 )
di3
Tt =
WX13
79
a = a l i l +a 2 i 2+a3 i3
Because of Equations (38.1) we can write the iast three terms in the
form
which reduces to w X a.
(38'.2)
da
aa
- =-+wXa
dt
at
'
where
_ dali + da 2i + da 3i
at - dt I dt 2 dt 3'
aa
(38.3)
Equation (38.2) expresses dajdt as the sum of two parts. The part
aajat is the time derivative of a relative to the moving coordinate
system. The part w X a is the time derivative of a relative to S when
a is fixed relative to the moving coordinate system.
When the origin of the coordinate system is not at rest relative to S
but has a velocity u, Equations (38.1) still hold, and hence also does
Equation (38.2).
39. Linear and angular momentum. Let us consider a particle of mass
m, with a position-vector x relative to a point 0 fixed in a frame of
m
o
Figure 46
80
M=mv.
The angular momentum of the particle about the point 0 is by definition the moment of M about O. We shall denote it by the symbol h.
Hence, by 10 where the moment of a vector about a point is considered, we have
h
xxM
xXmv
= mxXv.
L mjVj,
j=1
L mjxjXvj.
j=1
m dxc
dt
mj dXj.
dt
j= 1
But
dxc
Tt =
where
Vc
dXj
Vc,
dt =
Vj,
(39.2)
mvc = L mjvj = M.
j=1
81
h = 2:
(39.3)
mj
Xj X (00 X Xj)
j= !
h = 2: mj[oox/-Xj(Xj. 00)].
j=!
Now let us introduce coordinate axes with origin at the point 0 fixed
in S. The directions of these coordinate axes need not be fixed in S.
As before we denote the coordinates of the j-th particle by (Xji' xj2 ,
Xja). The component hi of h then has the value
N
hI = 2:
mj
j=!
[<ill
(XjI2+Xj22+Xja2) - Xjl
<ill 2:
<il2 2:
mj (Xj22+Xja2) -
N
mj Xjl Xj2 -
j=!
j=!
<ila 2:
j=l
mj Xja Xjl
(39.4 )
h2
ha
82
vc+oo X "i.
mj XjX (WXXj).
j=!
By Equation (35.2) the first sum is equal to mxc. Since the origin of
the coordinate system and the center of mass C of the body coincide,
Xc = O. Thus
N
= j=!
~ m- xJ J
(w xx-)
.
J
The right side of this equation is the same as the right side of Equation
(39.3). Hence in the present case the components ofh are also given
by Equations (39.4).
We have then the important result: Equations (39.4) may be usedfor
the determination of the components of the angular momentum h of a rigid
body about either a fixed point 0 in the body or the center of mass C of the body.
In the two cases the origin of the coordinates is at 0 and C, respectively, the
directions of the coordinate axes being quite general. Equations (39.4)
cannot be used in the case of the angular momentum of a rigid body
about a moving point which is not the center of mass of the body.
(40.2)
~
j = I
dv'
mj _ 1
dt
~ Fj+ ~
j= I
~ Fjk.
j= I k = I
dM
dt
'
where M is the linear momentum of the system and F is the sum of
the external forces acting on the system.
Theorem 2. The center of mass of a system of particles moves like
a particle with a mass equal to the total mass of the system acted upon
by a force equal to the sum of the external forces acting on the system.
Proof. In 39 we saw that M = mvc, where m is the total mass of
the system, and Vc is the velocity of the center of mass of the system.
Thus Equation (40.3) can be written in the form
(40.3)
--=F
(40.4)
This completes the proof.
41. The motion of a rigid botty with a fixed point. Let us now consider a system of particles which constitutes a rigid body with a point
o fixed relative to a Newtonian frame of reference.
Theorem 1. The rate of change of the angular momentum of the
body about 0 is equal to the total moment about 0 of the external
forces.
Proof. Let us introduce coordinates with origin at O. Then
(41.1)
84
v.
J
dXj,
dt
where vi> Xj and t have the usual meanings. By Equation (39.1), the
angular momentum h of the body about the fixed point 0 is
N
h=Lmxxv.
j=IJ]
J'
and so
dh
-=A+B
dt
'
(41.2)
where
N
h.
A = L mj -d1 X Vj ,
j= I
t
B =
N
h.
L m X X _ J
j= I J J
dt
A = L mjvjXvj = O.
j=1
Equation (40.1) gives an expression for mivj/dt. Because of this we
have
N
G = j=
L IxxF
]
J'
(41.3)
H = L L xjXFjk .
i=1 k=1
It will be recalled that Fj is the external force acting on the j-th particle and Fjk is the internal force exerted on the j-th particle by k-th
particle. We note that G is the sum of the moments about 0 of the
external forces. Now
N
H = L L XjXF;k = L L XkxFkj.
j=1 k=1
k=1 j=1
Thus
N
(41.4)
85
But F kj = -Fjk
2H
J{
J{
1:
(Xj- Xk) X
Fjk
j= 1 k=l
Since the lines of action of the vectors Xj - Xk and Fjk coincide, their
vector product vanishes. Hence H = 0, and B = G, so Equation
(41.2) reduces to the form
(41.5)
where h is the angular momentum of the system about the fixed
point 0, and G is the total moment about 0 of the external forces.
This completes the proof.
We have placed the origin of the coordinate system at the fixed
point O. Let us now choose as coordinate axes a set of principal axes
of inertia of the body at O. (Principal axes of inertia are defined in
36.) Then the products of inertia K1, K 2, Ka all vanish, and from
Equations (39.4) we obtain for the components of the angular momentum h ofthe body about 0 the expressions
hI
(41.6)
= l 1 w1 ,
h2
= 12w2 ,
ha
= 13 w3'
where 11 ,12 ,13 are the moments of inertia ofthe body about the coordinate axes, and WI' W 2 , W3 are the components of the angular velocity W of the body about O.
In most cases the coordinate axes will be fixed in the body and will
hence have an angular velocity w about O. However, in a few special
cases when the body has a certain symmetry it will be found possible
and desirable to choose coordinate axes not fixed in the body. To
include such special cases we denote the angular velocity of the axes
about 0 by 02, which mayor may not differ from w. According to
Equation (38.3) we then have
dd~ = ~~ +o2xh
or
(41. 7)
dd~
iz1i1
+iz2i2+iz3i:i+(02h3 - 03h2) i l
86
I'rom this equation we can read off the components of the vector
dh/dt. According to Equation (41.5) these components are equal to
the components of G. Hence we have the equations
~1+02ha - Oah2 = G1 ,
h2+0ahl - 0lha = G2 ,
Ga ,
(41.8)
In the case when the coordinate axes are fixed in the rigid body,
then Q. = wand so (12.8) reduce to the form
11 w1 - (12 - 1a) W 2W a = G1 ,
12 2 - (la - 11) wawl = G2,
1awa - (11 - 12) W 1W 2 = Ga.
(41.9)
2: mx = m Xc
j-I
87
mgj
Figure 47
= OCxmgj,
= ai3, j = il sin
('J. -
i3 cos
('J..
Thus
=
('J. -
i3 cos ex)
so
(41.10)
G1 =0,
Ga =0.
(41.11)
The angular velocity w of the sphere about
is
w = <1> i 2 ,
WI
= 0,
W2
= <1>,
W3
= 0.
The moments of inertia of the sphere about the coordinate axes are
(41.13)
II = 12 = tma2, 13 = tma2.
We now substitute in Euler's equations (41.8) from Equations
(41.10), (41.11), (41.12) and (41.13) to obtain the relation
i> _
-
5g sin ('J.
7a
'
89
14a
t,
z = !4 gt2 sin
(41.14)
QI;.
I'
-R
Figure 48
plate as the origin of the coordinate system. We also choose the unit
vectors il and i2 parallel to edges of the plate, as shown. The unit
vector i3 is then perpendicular to the plane of the plate. The directions
of these three vectors are principal directions of inertia of the plate
at o. The moments of inertia of the plate about the coordinate axes
are
90
WI
cos IX,
w2
= w sin
IX,
W3
= 0.
Since the coordinate axes are fixed in the body, Euler's Equations
(41.9) apply. We substitute in these equations from (41.15) and
(41.17), obtaining the relations
G1 =0,
G2 =0,
IX.
2cR
= G3
By Newton's third law ( 31), the forces exerted on the right and
left bearings are -R and R, with magnitudes R given in Equation
(41.18) above.
axis of symmetry is large. For example, the disc and shaft in Figure
49 constitute a gyroscope.
Figure 49
In Figure 50 the line OA is the axis of the gyroscope, the fixed point
being at 0 and the center of mass being at C. We introduce a fixed
unit vector j pointing up from 0, and a set of moving orthogonal
A
o
Figure 50
G = (lia) X (-mgj).
92
But
(41.19)
whence we get
G = -lmg sin a i 2.
Thus
(41.20)
G2 = -lmg sin 6, Ga = O.
The angular velocity w of the gyroscope is
Gl = 0,
= sia+pj.
whence
(41.21)
<ill =psin6,
<il2 =0,
<ila =s+pcos6.
The angular velocity .Q. of the coordinate axes is given by the relation
.Q.
= pj
= P (il sin 6+ia cos a).
Thus
(41.22)
0 1 = p sin a,
O 2 = 0,
Oa =
P cos a.
11P sin a = 0,
(41.25 )
(41.26)
la (S -
Pcos a) = o.
P = constant,
= constant,
(41.27)
O.
lmg
- laP
S---
42. The general motion ofa rigid body. We now consider a rigid body
moving in a general manner. It mayor may not have a fixed point.
The motion of its mass center can be determined from Theorem 2 of
40, which applies to the motion of any system of particles. This
theorem yields
m dvc = F
dt
'
where m is the total mass of the body, Vc is the velocity of its center of
mass, and F is the sum of the external forces acting on the body.
Integration of (42.1) gives the position-vector Xc of the center of mass
C of the body as a function of the time t.
To determine the complete motion of the body it is then only necessary to find the angular velocity of the body about its center of mass.
(42.1 )
94
d/
= j=1
2: mj"j X Vj'
(42.3)
where
(42.4)
N
d,,'
= j=1
2: m _ 1 XV
1 dt
},
B =.2: mj"j
J=I
dv.
-d} .
t
= j=1
2: m.(v1 J
vd XvJ
}=I
= O-vc xM ,
dh
dt -- G
We have thus the result: the rate qf change qf the angular momentum qf
a body about its center qf mass is equal to the total moment qf the external
forces about the center qf mass.
95
Figure 51
96
+ -27t2/11
W =S13
7taV
Thus
(42.6)
W3
S.
v.
all
whence
(42.7)
Also
(42.8)
We now substitute in Euler's equations (41.8) from Equations
(42.6), (42.7) and (42.8) to obtain the relations
SV
G2 = -13 -,
a
G3
= 0.
whence
Problems
1. A particle moves on the curve x 2 = h tan kx 1 , X3 = 0, where h
and k are constants. The X 2 component of the velocity is constantFind the acceleration.
97
98
where p and q are constant vectors and t is the time. Find the velocity v and position-vector x of the particle in terms of t, given that v
= 0 and x = 0 when t = O.
12. A particle of mass m is acted upon by two forces P and Q.
The force P acts in the direction of the Xl axis. The force Q makes
angles of 45 with the axes of X 2 and X 3 Also P = P sin kt and Q. =
q cos kt, where p, q and k are constants and t is the time. At time
t = the particle has coordinates (b,O, 0) and is moving towards the
origin with a speed pjmk. Find the position-vector x of the particle.
Prove that the particle moves on an ellipse, and find the center and
lengths of the axes of the ellipse.
13. A particle of mass m moves under the action of a force pe-qt
and a resistance -lv, where p is a constant vector, q and l are positive
constants, t is the time, and v is the velocity of the particle. Prove that
xoo - x.
rq1 (p+ m q u)
where u is the velocity when t = 0, and lito and Xoo are respectively the
position-vectors of the particle when t = 0 and when t becomes infinite. Is the above result true when l = mq?
14_ A particle of mass m moves under the action of a force p cos qt
- kx, where p is a constant vector, q and k are positive constants,
t is the time, and x is the position-vector of the particle relative to a
fixed point O. When t = 0, the particle is at 0 and has a velocity u.
Find x in terms of t when (i) k =1= mq2, (ii) k = mq2.
IS. A particle of mass m is acted upon by a single force y mjxZ
directed towards a fixed point 0, where y is a constant and X is the
99
mgt2
h
. hI
d
turne d t h rough t h e ang1e a(m+2m')
' were g IS t e acce eratlOn ue
to gravity.
101
CHAPTER
IV
PARTIAL DIFFERENTIATION
43. Scalar and vector fields. Let V denote a region in space, and let
X be a general point in V. Let Xl, X 2 , X3 denote the rectangular cartesian coordinates of X, and let x denote the position-vector of X. Then
x = xlil+X2i2+X3i3'
Let us now consider the case when there is associated with each
point in the region Va value of a scalar j Then we write f = f(x l ,
X 2 , X 3 ), or more compactly
f=f(x).
The values off associated with all the points in V constitute a scalar
field.
Let us now consider the case when there is associated with each
point in the region V a value of a vector a. The:r;t we write a = a (Xl'
X 2 , x a), or more compactly
a = a(x).
f=f(x, t),
= a(x, t).
;,
Figure 52
dj = dJ dX l + dJ dX 2 + dJ dXa.
ds
dX 1 ds
dX 2 ds
dX a ds
Let us now consider all curves through X having the same tangent at
X as C has there. For each of these curves the value at X of the right
side of Equation (44.1) is the same. Thus at each point X there is
(4tl)
a unique value of
(44.2)-
"(7.d
v = 11 -::;vX 1
od
+ 12od-::;+ 1a-::;-
vX
vXa
2
\l
or
1I"
\lJ =
~f
vX 1
11 -::;-
dJ
~f
+ 12 -::;+la -::;-.
dX
vXa
0
103
dJ =
ds
Vf dx,
ds
or
(44.4)
-dJ
ds = vrt
J'
/VJ/
cos 6 =
Vft
or
(44.5)
/VJ/
df
cos 6 = ds'
by Equation (44.4).
(45.1)
V =2: i r
r=1
0
-
oX r
Then
Vf=
(45.2)
r=1
ir
of.
OXr
Proof.
= Vf+Vg.
V (f+g)
(45.3)
()
V(f+g) =r~1 ir
oXr (f+g)
~f + ~ ir {)e
oX r r= I oXr
= ; ir
r= 1
= Vf + Vg.
Theorem 2.
Vf=J~Vu.
(45.4)
Proof.
of
r=l
OX T
Vf=2: i r -
Butf is a function of u, which is a function of X r Thus
of dj ou
~ = -d"""l
(r = 1,2,3),
oX
U OX
r
so
105
V']=~jrdjdU
r=\
dj
du
du
dX r
~ i.
r=\
dU
dX r
dj
= du V'u.
Theorem 3.
U1 ,
U2 , , Un,
then
(45.5)
U1 , U2 , , Un
which
Hence
106
dXr
Theorem 4.
v (fg) = fv g+gV f
(45.6)
Proof.
We have
3
d
V(fg) =r~li,. dX (fg)
r
g
f
ir(fd +gd )
r= 1
dX r
dX r
=fVg+ gVf
=f
(46.1)
aV = (alil+a2i2+aaia).(il
dX 1
d
d
d
= a 1 -::j + a2 -::j + aa ~
uX 1
UX 2
UXa
3
+ i2
uX 2
+ ia
UXa
=:E ar -
r=1 dX r
This operator is a scalar, and can be applied to a scalar field or to a
vector field. Thus, iff and b are two fields, then
3
df
3
db
(46.2)
(aV)f=:E ar-::j, (aV)b=:E a,-::j'
r= 1 UX r
r=1 UX r
We note that
c
or
or
107
it
at
i3
i2
a3
a2
(46.5)
ax\7 =
d
d
d
dX t
dX 2 dX 3
The operator a X \7 is a vector operator. It can be applied to scalar
fields. Thus, ifJis a scalar field, then
(46.6)
(aX\7)J=
We note that
it
at
i2
a2
i3
a3
dJ
dX t
d.f
dXz
d.f
dX 3
(46.7)
\7 . b = (it
~ + i2 dX2
~ + i3 dX~ 3 ). b .
uX t
(46.8)
\7b
=(~ 1 i.~).b
=~ 1 i.. dX
db
dX r
r
r=
r=
But
i.. b = br
Thus
(46.9)
108
(r = 1,2,3).
=~ ~(iTb).
dX
r= 1
v X
J)
d+
o12d
0
= (011 ~
~ + 13 -;;OX
oX
oX
l
Vxb =
d)
3
db
3 d
3
:E Ir-;;- xb =:E lrX -;;- = :E -;;- (Irxb).
0
( r=1
OXr
r=1
OXr
r=10Xr
But
ilxb
ilx (blil+b2i2+b3i3)
= b2ia - b3i2
by Equations (8.9). Similarly we have
i 2xb = bail - bli 3, i3Xb
Hence (46.11) becomes
= bli 2 - b2ip
(46.12)
This can be written conveniently in the form
il
i2
i3
d
d
d
(46.13)
Vxb =
dX l
dX 2
dX 3
b2
bl
b3
The expression V X b is often called the curl of b or the rotation of
b, and is written curl b or rot b.
Theorem. 1. If a and b are two vector fields, then
(46.14)
V (a+b)
V . a+V . b,
Vx(a+b) = Vxa+Vxb.
Proof of Equation (46.14). From Equation (46.9), we have
3
d
V (a+b) =r~1 dX (ar+br)
r
_ ~ (dar + db r)
-r=1 ,dXr
dX r
= Va+Vb.
(46.15)
109
\7 X (a+b) =
il
d
dX l
a1 +b l
i1
d
dX l
a1
i2
i3
d
J
dX 3
dX 2
a2+b 2 a3+b3
i2
i3
d
d
dX 3
dX 2
a3
a2
il
i2
d
d
+
dXl
dX2
b2
bl
i3
d
dX3
b3
\7 xa+ \7xb.
Theorem 2.
(46.16)
Proof.
(46.17)
(a 2;X3-a3;X2) (iIb)
(a ~-a ~) (~b)
3
dX 1
dX 3
it
Figure 53
quantities will be denoted by Xl" X 2', xs' and i 1 ', i 2', is'. Figure 53
shows the coordinate axes and the associated unit vectors.
We now introduce the two operators
'\7
'\7,
oJ
+ 12oJ-::;- + ISoJ-::;-,
J +
J +
J
:1--;
12
11 -::;OX 1
0,
= 11
OX 2
uX 1
OX 2
uXs
IS~'
dX s
LetJbe a scalar field, and let b be a vector field. We shall now consider proofs of the relations
(47.1)
v1= VJ,
(47.2)
V' . b = V . b,
(47.3)
V' xb = V xb.
iI'
au i l +a 12 i 2 +a13 ia
= s=1
~ als is
There are similar expressions for i 2 ' and i a'. The three expressions
can be written in the compact form
(47.4)
i,.'
(r
au is
1, 2, 3) .
s=1
We note that the constant a,s is the cosine of the angle between the
two vectors i,.' and is.
The components of the vector i l relative to the system S' are seen
to be all , a 2l and aal. Thus
3
il =
~
s=1
asl is'.
i,. =
~ asr
s= 1
is'
(r
1,2,3).
Throughout the remainder of this section we will adopt the convention that latin subscripts range over the values 1, 2 and 3, as in the
above equations.
We now introduce a set of nine quantities called the Kronecker
delta. This set is denoted by the symbol 8", and is defined by the
equation
(47.6)
112
8,s = 1 ifr = s
=Oifr=s.
Thus &11 = &22 = &33 = 1 and &23 = &31 = &12 = &32 = &13 = &21 = O.
Since iI, i2 and i3 are unit orthogonal vectors, they satisfy the nine
relations (7.3). These nine relations can now be written compactly
in the form
(47.7)
We also have the nine relations
(47.8)
Let us now substitute in Equation (47.8) from (47.4). This gives
the relation
&rs
= (
art
t= 1
it) .
(~anl i.)
.=1
= t=1
~ ~ art ani iti..
u= 1
But by (47.7) we have iti.
= &tu
so
(47.9)
&rs=
t=1 u=1
art a su &tu
= t=1
~ art ~ ani &t.
.=1
(47.lO)
~ ani tu
u=1
a st ,
(47.11)
113
(47.12)
3rs
= 1=1
1: a"
als'
= ,=1
1: bs' is' .
b = 1: b, i,.
(47.13)
,=1
We now substitute here for i,. from (47.5), obtaining the relation
3
1: b/ is'
r=1
= r=l
1: b,
3
= 1:
1: asr is'
s=l
1: b, asr is' .
,=1,=1
(47.14)
b,
= r=l
1: ars b/.
Xs '
= L asr x"
,=1
x,
= ,=1
1: ars x/.
(47.17)
{)X',
()X s
""::)"7
oX r
-:s-= an ,
oX
r
aTS
\11= \1f
Now
\I1=~i'r{)!'
1
{)X
(47.18)
r=
and
Because of this relation and (47.4) we can write (47.18) in the form
(47.19)
\11 = ~ ~ ~
r=1
,=1 '=1
{)Xt {)X r
art
\11 = T=l
~ ~
s=1
3
t=1
3.
a"
{)j
{)j
art Is '""l
ox,
3
= s=1
~ ~ 1, '""l ~ a" aTt
t=1
oX r=1
t
\11 =
~ Is'""l
,=1 t=1
3
oX,
a,t
{)j
= s=l
~ Is ~ '""l a" .
t=1 oXt
The last summation here reduces to {)jj {)x1 when s
when s = 2, and to {)jj {)xa when s = 3. Thus
1, to {)jj {)X'I.
115
(47.20)
and so
3
df
\I'j=L is s=1
= \If.
dx,
\Ib.
Now
3 db'
\I'b =L - ' .
,= 1 dx's
\I'b =L L -
=~ ~
.1=
1 r= 1
(a b)
.rr
a db,
ST
()x's
= L
db
oX I
L ~ L a" ast
\7'.b =
~ ~ ob, ~'I
,=1 1=1
OXt
~ ob"
OX,
by a repetition of the arguments leading up to Equations (47.10) and
(47.20). This completes the proof.
Proof of Equation (47.3). This is left as an exercise for the reader
(Problem 13 at the end of the present chapter).
.=1
(48.1 )
\7.(fa) =f(\7a)+(\7f).a,
(48.2)
\7 X (fa) =f(\7xa)+(\7f)xa,
(48.3)
\7.(axb) =b(\7xa)-a(\7xb),
(48.4)
(48.5 )
\7(a.b) = (a. \7) b+(b. \7) a+a <(\7 Xb)+bx (\7 xa),
(48.6)
\7 X (\7f) = 0,
(48.7)
\7(\7xa) =0,
(48.8)
\7X(\7xa) =\7(\7.a)-(\7.\7)a,
(48.9)
\7.x=3,
(48.10)
\7xx=O.
(48.11)
(a. \7)x = a.
117
\7.(ja) = L
,=1
=L3
,=1
1,- (fa)
dx,
1,'
f a)
~Jda +ddx,
dx,
3 . da
1,'-
=fL
,=1
dx,
3.d!
+ ,=1
L 1 , - " 'a
dX,
= f(\7a) +(\71) a.
Proof of Equation (48.3). Because of Equation (46.8) we have
3
d
\7 . (~ X b) =,~ 1 i,.. dx, (a X b)
= L3 i, (da
- Xb )
dX,
,=1
+,=1
L3 i,.. ( a X-db) .
dx,
Because of the permutation theorem for scalar triple products (Theorem I of 9), we can write this last equation in the form
da'
3
(
db)
\7(axb) =L3 b ( i,.X-)La'
i,.xdx,
,=1
dx,
,=1
da
3.
db
=bL 1,x--a1: I,X3.
dx,
,=1
,=1
dx,
uX,
We now apply to the right side of this equation the identity (9.3)
for vector triple products, obtaining
(48.12)
d
ax(\7xb) =L3 i,. (db)
a-- - L3 - b (ai,)
,=1
dx,
,=1 dX,
3.(a 'db)
- - L3 a -db
=L
1
,=1'
dX,
,=1' dx,
=L3 i, ( a 'db)
- -(a\7) b.
,=1
118
dX,
(da )
3
b X (\7 xa) =r~l
iT dxT',b - (b \7) a.
()
(48.14)
aX (bxe)
b(a.e) - e(a.b).
\7 (\7.a),
\7 (a.\7),
(\7.a)\7,
we find tha,t the term corresponding to the second term on the right
side of (48.14) is (V, V) a. Thus
V V =
{)2
(()X )2
1
{)2
V(aXb)
Va(axb)+Vdaxb).
V(axb) =b(Vaxa)-a(Vbxb).
But
Vaxa = V Xa,
Thus Equation (483) is proved.
Zl =11 (Xl'
X 2 , X 3 ),' Z2
=12 (Xl'
X 2 , X 3 ),
Z3
=13 (Xl' X 2 ,
X 3 ),
coordinates Xl' x 2 and Xa a new set of coordinates Zl, Z2 and Z3' These
new coordinates are called curvilinear coordinates, and Equations
(49.1) are equations of transformation of coordinates.
We now propose to compute, in terms of quantities pertaining to
curvilinear coordinates only, the expressions \7J, \7. b and \} xb,
where 1 and b denote respectively a scalar and a vector field. With
this goal in mind, we shall devote the rest of this section to some preliminary considerations, and shall complete the final computations
in the next section.
The Jacobian of the transformation (49.1) is the determinant
(49.2)
]' =
dZ l
dX 1
dZ l
dX 2
dZ l
dX a
dZ 2
dX 1
dZ 2
dX 2
dZ 2
dXa
dZ a
dX 1
dZ 3
dX2
dZ a
dX a
We shall consider only the case when ], does not vanish anywhere
in V, so that Equations (49.1) may be solved 1 for Xl, x2 and Xa to
yield the relations
(49.3)
Xa =ga (Zl,Z2,Za)'
121
iarly, we have parametric lines of Z2 and ZS' For example, in the case
of cylindrical coordinates r, e, xs , the parametric lines of r, e and Xs
are respectively horizontal straight lines cutting the Xa axis, horizontal
circles with centers on the Xs axis, and vertical straight lines.
Let X be a general point in a region V. Through X there passes a
parametric line of each of the curvilinear coordinates ZI, Z2' Zs' We
now introduce three unit vectors kl' k2 and ks with origins at X,
defined as follows: kl is tangent at X to the parametric line of ~l'
and points in the direction of Zl increasing; k2 and ks are defined analogously with respect to the par_ametric lines of Z2 and Zs' Figure 54
k,
X~_""""=-+
Figure 54
(49.4)
T=l
122
(r
= 1,2,3).
Because of (ii), (iii) and Equation (49.6) it follows that l'i7z11 = ljh 1 ;
and because of (i) we then have 'i7z 1 = k 1 /h 1 Similar observations
regarding Z2 and Z3 then permit us to write
.
(49.7)
From these equations we get the relation
k 1 (k2 X k s)
But the left side of this equation is equal to one since kl' k2 and ks
form a right-handed orthogonal triad of unit vectors. Thus
123
(49.8)
+ df
CJ VZa
uZa
3 1 df
Vf = ~ -h -;:;- k,.
,-1
uz,
k1 = k 2 xka ,
k2 = kaxk1,
ka = k 1 xk 2 ,
k1 = h2h a V Z2 X V Za ,
k2 = hah1 VZa X VZ 1,
ka = h1h2 VZ1 X Vz 2
We now substitute from these expressions in Equation (50.3), and
operate on the resultant expression with the operator V to obtain
the relation
124
(50.4)
For the first of the three terms on the right side of (50.4) we can write,
because of relation (48.1),
(50.5)
d
\7 (b 1h2 ha) =""'l (b 1h2ha) \7 Zl
~
d
+ ")
(b 1h2ha)
~
\7Z2
d
+ ")
(b 1h2ha)
~
\7 Za'
perpendicular, we obtain
\7 zd\7 Z2 X \7 za)
d
")(b 1h2ha)'
a oZI
\7. b
=.
h h1h [d
-=j (b 1h2 ha)
I 2 a oZI
d (b 2hahl) + Td (b ahlh2) ]
+ ToZ2
oZa.
or
(50.9)
For the first term on the right side we can then write
(50.10)
\l ZI + b1hl (\l
\l ZI) ,
d
";) (b hl ) \l ZI
[ uZI 1
\lZI =
+ lza (b hl)
1
d
+ ";)
(b hl) \l Z2
uZ2
1
\lza]
\lZI'
\lZ2 X \lZI
\lZ3 X \lZI
0,
1
h2hl k2 Xkl
=-
1
hahl ka X kl
= hahl k 2
h2hl ka,
1
This relation and two similar relations involving the second and third
terms on the right side of Equation (50.9) permit us to write Equation
(50.9) in the form
(50.11)
[d
kl
\l X b = h~a dZ 2 (bah a) - dZa (b 2h2 )
d ]
+ hak2hl [d
dZ a (b hl) - dZ I (baha)
1
ka[d
,.
d
]
+ hlh2
dZ I (b2h2) - dZ 2 (b hl)
1
126
Vx = nX'2",
where n is a constant.
7. Find 'Vr and 'V6, where rand 6 are the usual plane polar coordinates. Also, find the magnitudes and directions of 'Vr and 'V6.
8. If f = r3 - cos 2 6, where rand 6 are plane polar coordinates,
find 'V f in terms of r, 6 and the unit vectors i l and i2 associated with
the corresponding rectangular cartesian coordinates.
9. Iff and g are scalar fields, prove that
'V(j/g)
10. Iff = (Xl)2+X3 V(Xl)2+ (X2)2 and g = X l X 2Xa, find at the point
A (3,4,5) the expressions 'V (jg) and 'V (j/g). Note Theorem 4 of
45, and Problem 9 above.
11. Iff = XlX~3' a = xlil - x 2i 2 and b = xsxli2 - xlxJa, compute
the following: (i) (a 'V)f, (ii) (a. 'V) b, (iii) (a X 'V)f, (iv) (a X 'V)' b,
127
t7t
128
t7t
27. If r, e, cp are spherical polar coordinates, describe the parametric surfaces and lines, and show that
(dS)2 = (dr)2+ r2(de)2+r2 sin2e (dcp)2.
28. Prove that for the transformation from rectangular cartesian
coordinates to orthogonal curvilinear coordinates Zl' Z2, Za for which
the metric form is as given in Equation (49.5), the Jacobian I satisfies
the relation hlh2hsl = 1.
29. Express Equation (50.11) in terms ofa determinant.
30. Write out the expressions VI, V b and V xb in the case of
cylindrical coordinates r, e, z.
31. W ri te out the expressions VI, V b and V X b in the case of
spherical polar coordinates.
32. By setting b = Vf in Equation (50.8), deduce an expression
for V2f in terms of general orthogonal curvilinear coordinates.
33. Show that, in cylindrical coordinates r, e, z, we have
2
V2f = d2f !. d2J
df
N.
dr2
r2 de 2
dZ 2 r dr
+1
d;2
d;a
= 0,
0,
~2 = -kl>
~2 = ka cos e,
~~a =
~a =
0,
e.
Note: the corresponding problem in the case of plane polar coordinates is worked out in Chapter III, 30.
129
CHAPTER
lNTEGRATION
51. Line integrals. A curve is called a regular arc it can be represented
in some rectangular cartesian coordinate system by the equation
x
xI(u)il+X2(U)i2+X3(U)i3'
-~B
K.
Q..
Q,
Q.
Figurtl 55
(51.1)
f 1 ds
C
t>.sp .... 0
Figure 56
Let s denote the arc length of C measured from the end A of C. The
vector dx/ds was seen in 28 to be a unit vector tangent to C in the
direction of s increasing. Denoting this vector by t, we have
(51.2)
dx
= ds .
bt = bt.
(51.3)
(51.4)
f b, ds
c
Because of (51.2), we have
= Jb.t ds.
c
(51.5)
Thus
fbtds = fbodx = f(bldxl+b2dx2+b3dx3)'
c
c
c
We maY,also consider line integrals over C with integrands which
are vectors. The integration of vectors was defined in 13. Following
this definition, we have
(51.6)
(51. 7)
1c b ds = i
131
(51.8)
jbxtds = jbXd:x
c
= i1
!c (b dxa-badx
2
2)
jids
(ii)
On C we have
ds
X2
= 2x1,
V1+ (:::Y dX
= (X1)2+8(X1)3,
1
Vs dx
Thus
2
jids
c
= ![(x1)2+8(x1)a] Vs dX 1
= lO;v'5.
132
(i)
(51.9)
We have
bdx = bldxl+b2dx2+badxa
= x2dx I + (xa +XI )2 dx2+xl dxa
=
But on C we have
[X 2+(Xa+X I)2
X2
2x 1, Xa
0, so
jb.dx
C
(ii)
8
23
But on C we have Xl
tX2' Xa = 0, so
4
j bdx
C
= j[tX2+!-(X2)2] dX2 =
8
23
G,a)
"
Figure 57
from the point A to the point B, the curve C being a portion of the
intersection of the cylinder (X I )2+(X2)2 = a2 and the plane xI+xa = a.
133
Now
But on C we have
(X2)2
= a2 -(Xl)2,
Xa =
Xl
dX a
dX l
X2
dX l
-=--,
Thus
a-
dX 2
Xl
= -1.
j bdx = j[a2x l
C
a(a -
Xl)Xl - Xl
(a -
Xl)2]
dX l = - ia4..
Figure 58
(52.1)
lim
2;
N-+oop=l
ASp -+ 0
1--+---',
. . . ----5'---X'
"
OS'
c'-----Figure 59
dS' = nadS.
Let the equation of the surface S be
(52.2)
Denoting the left side of this equati.on by C, we have by Theorems 2
and 3 of 44,
DCI DC dC. dC. dC.
n
Iv
= -
= -::;11 +
oX
I
dXa
Til -
uXI
dXa
-:J
uX 2
~ 12
uX 2
+ oXa
.;:;- la
1 2 +1 a
135
Thus
(52.3)
We then have
(52.4)
(52.5)
fbndS=fb.ndS.
s
s
It is often convenient to introduce an infinitesimal vector dS defined
by the relation
(52.7)
dS
= n. dS,
fbndS = fb.dS.
s
s
It is sometimes necessary to consider surface integrals with integrands which are vectors. For example, we have, following the definition of integration of vectors in 13,
136
(52.9)
jb dS = iljbldS+i2Ib2dS+iajbadS,
s
s
s
jb xndS = ilj (b2na-ban2)dS+i2j (bancblna)dS
s
(52.10)
-,
"
Figure 60
rJxa =
We have
{}xa
dX 1
'dX 2
whence (52.3) yields na = 1/3. Thus dS
= _2
= 3 dx l dx 2, and so
j]dS = j 3] dXI dx 2
S
But on S, we have
Xa
]
= 2 - 2XI - 2x2 ,
= (X I )2+2x2+ (2 - 2x l
2x 2) - 1
= (Xl - 1)2.
137
X IX2
plane, we have
I-x,
S'
1bndS = 21 (x 2+xa)dxI dx 2
But on S we have Xa
= 2 - 2XI -
2x2 , whence
1
I-XII
s'
0 0
53. Triple integrals. Let V be a region in space enclosed by a surface S as shown in Figure 61. Let.f (Xl' X 2 , xa) be a function which is
Figure 61
138
v
v
v
v
To evalute triple integrals, one may divide the region V into elements by means of three systems of planes parallel to the coordinate
planes. The value of the integral is then found by the performance of
three integrations with respect to the rectangular cartesian coordinates
Xu X2 and Xa. Or we may divide V using parametric surfaces of a
curvilinear coordinate system, in which case we evaluate the triple
integral by performing three integrations with respect to the three
curvilinear coordinates. Since most readers will have had considerable
experience in evaluating triple integrals in connection with elementary
calculus, no example need be given here.
Problems
1. If x is the prosition-vector of a general point on a circle C of
radius a, and t is the unit tangent vector to C, evaluate
dx.
ft.
2. If] = (XI)2 - (X2)2, evaluate the line integral of] along the line
= 2 from the point A(O,I) to the point B(2,0).
3. Evaluate the line integral in Problem 2 when the curve C consists
of (i) the two line segments AD and DB, where D has coordinates
(1,1), (ii) the two line segments AO and OB, where 0 is the origin.
4. If ] = 8lxI - 9, evaluate the line integral of] along the curve
(X2)2 = (X I)3 from the origin to the point (1,1).
5. tf] = X2Xa+X3XI +XIX2' evaluate the line.integral of] from the
origin 0 to the point B( 1,2,3) along the path consisting of (i) the line
segment OB; (ii) the three line segments OD, DE and EB, where D
and E have coordinates (1,0,0) and (1,2,0), respectively.
Xl +2X2
139
~.
If J =
xI +X z+X3 ,
(i)fxds,
(ii)fxxtds.
c
c
13. If J = Xl +X a , evaluate the surface integral ofJ over the region
S consisting of the triangle cut from the plane 6x 1 +3x2 +2x3 = 6 by
the three planes Xl = 0, X 2 = 0, Xl +X2 = 1.
14. If b = xliI + (X2) 2ia, evaluate the surface integral ofb over the
region S in Problem 13, the origin being on the negative side of S.
15. If J = (X I)2+(X Z)2, evaluate its surface integral over the region
S consisting of the part of the surface X3 = 2 - (XI)2 - (X2)2 in the: first
octant. Use polar coordinates in the X I X2 plane.
16. Find the area of the region S in Problem 15.
17. If b = (x l)2i 1 +x3i a ,evaluate the surface integral ofb over the
region S of Problem 15, the origin being on the positive side of S.
"
i,
c
. 0 ' - - - - + _ - - x,
i,
Figure 62, and let b be a vector field continuous and with continuous
first derivatives in the region S. Then
!i (\7Xb) dS = c!b.tds,
(54.1)
db2)
j (bldxl~b2dx2)'
j (dbl
dX - dX dS = I
2
c
jsdbdx. dS = - cj bIdxl>
l
in the case when G can be cut by a line parallel to the X 2 axis in two
points at most. Figure 63 illustrates the situation. On G there are
two points D and E where the tangent to G is parallel to the X 2 axis.
Let d and e be the abscissas of D and E, respectively. These points
divide G into two parts G' and Gil. At a general point X(x l , x 2 ) in S
we introduce an element of area lying in a strip parallel to the X 2 axis,
the left edge of the strip cutting G' and Gil at the points X' (Xl' x 2')
and X"(XI, X2/1), as shown. Then
141
Jb1(X1,
X2') dX 1 -
f b1(X1,
= -
X 2") dX 1
f b1(X1,
X2') dX 1 -
Jb1(x'1,
X2") dXl
Let us now consider the case when C can be cut by a line parallel to
the X 2 axis in more than two points, such as the case of the curve C
in Figure 62. Here we have only to join the points F and G where
there are tangents parallel to the X 2 axis by a curve K which is contained
O;M-L---d.~--1------------+I-'"
Figure 63
142
(55.1)
!\l.bdV=!b.ndS,
s
f ;;3
dV =
h3na dS ,
in the case when S can be cut by a line parallel to the Xa axis in two
points at most. Figure 64 illustrates the situation, T being the projection of S on the X 1 X 2 plane. On S there is a curve C consisting of points
where the tangent plane to S is parallel to the X3 axis. The curve C
cuts S into two portions S' and S". At a point X(Xl> X2 , xa) in V we
introduce an element of volume lying in a prism parallel to the X3
axis, the vertical line through X meeting S' and S" at the points X'
(Xl> X2, x3 ') and X" (xl> X2 , xt), as shown. Thus
143
(55.4)
=!
S'
S"
o.
_________ ---x,
--------T
Figure 64
Let n' be the unit outer normal vector at X', and let dS' be the area of
of the element cut from S' by the vertical prism. Let us define n" and
dS" at X" similarly. Then we have
dX 2 dX 1 = n3' dS' = -na" dS",
and so we may write (55.4) in the form
! ~:3 = !
dV
S'
=f b3nadS.
S
Let us now consider the case when S can be cut by a vertical line in
more than two points. In such cases we can always divide V into a
number of regions VI' V2 , by cutting V by a number of surfaces
Xl> K 2 , so chosen that the boundary of each of the regions VI)
144
!~!:dV= !b n
2 2
dS,
!~!3dV= !b n
S 3
dS.
dS,
or
(55.6)
!(\7.\7)]dV = ! t dS,
v
s
where \7. \7 is the familiar Laplacian operator, often denoted by \7 2 ,
and d]jdn is the directional derivative of]in the direction of the outer
normal to the surface S.
56. The symmetric form of green's theorem. Let] and g be scalar fields
with continuous second derivatives in a closed region V bounded by
a surface S. We may then apply Green's theorem as stated in Equation (55.1), but with the vector b replaced by ] \7 g. This yields
(56.1)
But
fj~! dS.
fUV 2g+V!-Vg) dV =
v
s
(56.2)
f (g
vy+ Vg vj) dV =
fg
v
s
Subtraction of (56.3) from (56.2) then yields
dS.
f(i~!
f (jV2g - gVY) dV =
- g Z)dS.
v
s
This equation is called the vmmetric jorm oj Green's theorem.
(56.4)
(57.1)
fb.t ds,
f
S
= f (bIdxI+b2dx2+b3dx3)
f n (V X bIil ) dS
S
fbI dXl>
c
in the case when S is a regular surface element and the positive side of
146
S is that side on which the unit normal vector n points in the direction
of increasing Xa' Figure 65 illustrates the situation, and shows the
S'
C'
Figure 65
(57.4)
X 1X 2
plane into
~!l-ia~!l)dS.
sa
.. 2
Let the equation of the surface Sbe X3 = g(Xl' x 2). Then on S, we have
bl[Xl,X2,xa(XI,X2)]
dC 1
dX 2
(57.5 )
__
db 1
dX 2
C1 (XI,X2) '
db 1 dXa.
dX a dx 2,
(57.6)
j
s
n ('\1 X b1i l )
dS = - j
ni3~:~ dS+ j
s
= - /1 +/2 ,
n'
~!~ dS
where 11 and 12 denote the two integrals on the right side of this
equation.
Let us consider 11' We have n.i a dS = ns dS = dS', where dS' is
the projection of dS on the X 1X 2 plane. Since C1 is a function of Xl and
X 2 only, we can then write
147
11 = - f;~ dS'.
s'
11
=f
C1 (X 1 ,x2 )dx 1
=f
C'
b[x1> x
1
2,
X3 (X1>x 2)]dx 1
C'
!b dx
1
and so
(57.8)
12
dx db
1
n"')"')
dX 2 uX 3
dS.
dx
n-
dX 2
=0.
=f
(58.1 )
f \7.b dV = sf nb dS,
(58.2)
f(nX\7).bdS= ft.bds.
s
f \7JdV = f nJdS,
v
s
(58.4)
f \7 X b dV = f n X b dS.
v
s
(58.5)
f (n X \7)J dS = f tJ ds,
s
c
149
(58.6)
(58.7)
But we have
(58.8)
f \7 . (b
c) dV =
f n (b
c) dS.
v
s
Since e is a constant vector, we have by the permutation theorem for
scalar triple products the relations
The six integration formulas (58.1)-(58.6) may be written compactly in the form
! \7* T dV = s! n* T dS,
(58.9)
(58.10)
!(nX\7)*TdS=!t*Tds,
c
where T can denote a scalar field or a vector field, and the asterisk has
the following meanings: if T is a scalar field, it denotes the multiplication of a vector and a scalar; and if T denotes a vector field, it
denotes either scalar or vector multiplication. Thus, for example,
if T denotes a vector field b and the asterisk denotes scalar multiplication, then Equation (58.10) becomes (58.2).
59. Irrotational vectors. A vector field b (xu x 2 , x 3 ) is said to be
irrotational in a region V in space if everywhere in V we have
(59.1 )
\7xb =0.
Let Cjl be any scalar field with continuous second derivatives; and
let us write b = \7Cjl. Then
!b.tds =! n(\7xb) dS =
c
by Equation (59.1).
151
C'
Xo
Figure 66
above, the line integral of b from Xo to X is the same for paths G' and
Gil and hence has the same value for all paths in V from Xo to X. Thus,
if we write
x
(59.2)
cp = jb.dx,
Xo
drp
ds
Differen-
= b.dx.
ds
152
b - d'Pa _ d'P2 ,
I dX 2 dX a
(60.3)
b - ?'PI _ d'Pa ,
2 dX a dX I
(60.4)
b - d'P2 _ d'Pl ,
a - dX I
dX 2
for 'PI, 'P2 and 'Pa, where bl , b 2 and ba are given functions subject to
the condition
(60.5)
Let us choose 'PI = O. Then we have from Equations (60.3) and
(60.4) by partial integrations with respect to Xl'
(60.6)
a,
(60.7)
a,
_J
x
b =
I
a,
'(db 2 + db 3) dX I
dX 2
dX a
+ d~a _ d~2.
dX 2
dX a
153
h =
jb (al> X
2,
X:i)
0,
dx 2 ,
a,
= 0,
({I2
jbg(Xl> X
Xg) dXl>
2,
a,
({Ig
a1
where all integrations are partial integrations, and al and a2 are constants. The function cp is called a vector potential function.
In the above proof, several arbitrary selections have been made.
This indicates that a given solenoidal vector field b does not possess a
unique vector potential function. In order to see this more clearly,
we let cp be one vector potential function corresponding to the solenoidal vector field b, and let j be any scalar field. Then
Problems
1. Let C be a closed curve in the
X IX2
t f (Xl dX
c
2 -
x 2 dx l )
S, and n is the unit outer normal vector to S, prove that the volume
V of the region enclosed by S is given by the relation
V=tf nxdS.
s
3. If S is a closed surface with a unit outer normal vector n, prove
that
nxxdS =0.
s
13 =
!f
s
fs n(\7xb) dS =0.
6. If b = a l (xI)2iI+aZ(x2)2i2+aS(xs)2is, where al> a2 and as are
constants, evaluate the surface integral of b over the sphere through
the origin with center at the point A (al> a2 , as).
7. If b = (XI)2il+XIXJZ+xsis, evaluate the surface integral of b
over the cube bounded by the planes Xl = 2, X 2 = 2, Xs = 2 and the
coordinate planes.
8. If b = [(XI)2 - xJ i l +[2(x l )z+3xJ i z - 2Xlxaia, evaluate the
surface integral of b over the sphere S with cent~r at the point E
(1,0,2) and passing through the point F(3, -2, 1).
9. If C is any closed curve, prove that dx = o.
= 4,
Xs
0,
the fingers of the right hand when the thumb points in the direction
of the positive Xa axis.
11. Prove Equations (58.5) and (58.6).
12. A vector field b has continuous first derivatives in a closed
region V. On the bounding surface S of V, b is normal to S. Prove
that
!VXbdV=O.
v
13. If j and b are fields with continuous first derivatives in a region
V bounded by a surface S, prove that
!jnxb
! xb
dS = j\7
dV+ \7jxb dV.
s
v
v
14. If band care irrotational vector fields, prove that b X c is
solenoidal.
15. If b = x2xail+X3Xli2+XlX2i3, show that b is solenoidal, and
find its vector potential.
16. Show that the vector field
x1i 1+xJ2
(Xl)2+ (X2)2
Jj2
+-,
JX3
2
156
CHAPTER
VI
TENSOR ANALYSIS
l'
= I ~~' I '
the ranges for rand s being 1, 2, 3. We assume that l' does not vanish,
whence, as mentioned in 49, Equations (62.1) may be solved to
yield
(62.3)
(62.4)
{}
I,
dz" = ~ ~ dz'
s= 1 {}z'
(r=I,2,3).
I,
'+Z
{}Z',
d'"
{}z' (",
(62.6)
=
158
ifr. =1= s.
dZ dZ's
dZ's dZt '
T
8t
(62.7)
which is true because its right side is equal to di/dZI, which, because
of the independence of the coordinates i, satisfies the relations
dZ
dz l = 1
= 0
if r = t,
ifr =1=
t.
159
(63.1)
{} I,
~
AS ,
(}ZS
A"
(63.2)
I, {}
IS
A'rs~ ~ A tu .
- {}Zl {}z"
Contravariant tensors of higher order are defined analogously. Contravariant vectors are often called contravariant tensors of the first order.
A,=
{}ZS
(}Z"
As-
(64.2)
{}<I>
{}Z"
Z2,
z3). We have
{}<I> {}ZS
{}ZS {}Zl"
whence we note that the three expressions {}<I> I {}z' are the components
of a covariant vector.
The definition of a covariant tensor of order two is as follows:
160
65. Mixed tensors. Invariants. We can also define tensors whose law
of transformation involves a combination of both contravariant and
covariant properties. Such tensors are called mixed tensors. Thus,
for example, a set of twenty-seven quantities A:t is said to be the components
of a mixed tensor of the third order with one contravariant suffix and two covariant suffixes (f they transform according to the equation
(65.1 )
It should be noted that (65.1) represents twenty-seven equations,
and that the right side of each of these equations consists of the sum
of twenty-seven terms.
A single quantity A is said to be an invarinat ifit transforms according
to the equation
A' =A.
Invariants may be called contravariant tensors of order zero, or covariant tensors of order zero. An example of an invariant is the temperature in a fluid.
We shall now prove that the Kronecker delta has the tensor character indicated by its suffixes, that is, that it is a mixed tensor of the
second order. We must establish the equation
a:
(65.2)
161
Now
{}z' t
{}Zl t
{}Z" 8. = {}Z" 8,
{}Z 81
+ {}z'S
+ {}Z
81
{}Z'S
3'
When t is given the values 1, 2 and 3, the right side of this relation
reduces respectively to the expressions
{}z'
{}z'S'
{}z'
{}z'S'
{}z'
{}z's'
But this is equal to the left side of Equation (65.2) because of Equation
(62.7).
A:
B:
(66.1)
I t is easily proved that the sum of two tensors is a tensor of the same
order and type as the two tensors added.
There are two products of tensors, called the outer product and the
inner product. The outer product of two tensors is defined to be the
set of quantities obtained by multiplication of each component of the
first tensor by each component of the second tensor. Thus, for example, the outer product of the tensors
and
is the set of 243 quantities C;:. ~iven by the relations
A:
B:,
(66.2)
In writing this equation, we have been careful to keep each particular
162
suffix at the same level on both sides of the equation. When this
convention if followed, it is easily proved that the outer product of
two tensors has the tensor character indicated by the number and
positions of its suffixes. Thus, for example, in Equation (66.2), C;:., is
a mixed tensor of the fifth order, with two contravariant suffixes and
three covariant suffixes.
We now introduce an operation called contraction. It consists in
identifying a superscript and a subscript of a tensor. Thus, for example,
if from the tensor A:t we form the set of quantities B t defined by the
relation
(66.3)
we are performing a contraction. We note that each component
of B t is equal to the sum of three components of A:t It is easily
proved that contraction of a tensor yields a tensor Thus, for example,
B t in Equation (66.3) is a covariant vector. It should be noted that
we do not perform contractions by identifying suffixes at the same level,
since such operations do not yield terisors, in general. For example,
if A:t is a tensor, A:s is not a tensor, in general.
We now define the inner product of two tensors. To obtain it, we
form an outer product of the two tensors and then perform a contraction involving a superscript of one tensor and a subscript of the
other tensor. Thus, for example, an inner product of two tensors A~
and B" is the set of quantities Cst given by the relation
Cst = A:B,t
Of course the inner product of two tensors is not unique. Since outer
multiplication and contraction of tensors both yield tensors, the inner
multiplication of tensors also yields tensors.
A:
163
A:, =0.
(67.1)
Let <." be any other coordinate system, and let A ':, denote the components of this tensor for the coordinate system z". We must show
that
(67.2)
By the laws of tensor transformation we have
"
A sl
A:,
B:
A;, =
B~"
164
dZ'V ,
Xu = dZ u Xv'
or
(68.3)
dz" AtU)X' - 0
(A'TV _ dz'T
dZ t dZ u
v
We now replace the dummy suffixes rand s on the left side by t and u,
whence we obtain the equation
(68.4)
165
where
(68.5)
Equation (68.4) can be written in the form
(68.6)
IS
O.
Since X'is arbitrary, so are the three quantities X". Hence the coefficients in Equation (68.6) must vanish, which leads to the equation
(68.7)
From (68.5) we then have
I
A tu
+ Aut =
{}ZT {};t
{}Z'U
{}z't
{};t {}z'
ATS
In the last term on the right side, we interchange the dummy suffixes
rand s, obtaining the relation
(68.8)
A ' tu
{}z' {};t (A
A )
+ A'ut = {}Z't
{}z'u rs + sr'
A tu
{}z' (}z'
~ ~Ars,
oZ oZ .
(ds)'
(dx,)'+(dx,)"+(dx.)'
dx,dx,.
z: defined by relations
z'
= {}x'd
{}z' ;t,
(69.2)
bll = gll'
b33 = g",
(69.3)
l.b
2 23 =g23 -g
32,
Ib 31
"2
-g
31 -g
13)
t b"
= g12
g., .
(69.4)
(ds)'
= grsdz'dz'.
70. The conjugate tensor. Let g denote the detern;tinant whose elements are the components of the metric tensor. Then
g = Ig,sl
In the expanded form of this determinant, the coefficient of anyone
element grs is called the cofactor of g,S" We denote it by the symbol
f:::.TS. We note in passing that the minor ofgrs is equal to
(-1 y+s /j.rs.
(70.1 )
(70.2)
grs
/j.,t =
g 8!,
(70.3)
gs,
/j.t, =
g 8!.
determinant Igrsl, the subscript r varies over the rows, while the subscript s varies over the columns. We now consider Equation (70.2)
when s = t = 1. The right side is equal to g. The left side is
jj."
g"=-'
g"il
(70.6)
g"i'
= 3~.
3:,
168
gPrYr = gprgsrX' =
~:X'
= XP,
In appearance, this new tensor differs from the original one only in
that the suffix in the second position has dropped from the upper level
to the lower one. We say that we have lowered a suffix. Other tensors
obtained from A" by this lowering process are the following:
B:, = gt Btn
B;'
(72.1 )
A =Vg,)'A'.
A = VA.A' =
Vg" A,A,.
= V(A')'+(A')'+(A3)'.
Let X' and Y' be two unit vectors at a point. The angle 6 between
them is defined by the relation
(72.4)
cos 6 = g"X'Y'.
cos 6
= X'Y'+X'Y'+X'Y'.
(73.1)
z' =
1'(s) ,
the arc length s of the curve being used as the parameter. Because of
Equation (62.5) we see that dz'/ds are the components of a contravariant vector. We denote it by P', so
(73.2)
p' = dz'.
ds
p of
p2 = gmnpmp. = gm.
dz m dz
dS di
(dS)2
(ds)'
l.
Thus p' is a unit vector. It is called the unit tangent vector of the
curve C.
A geodesic may be defined as the curve of shortest length joining
two points. In our three-dimensional space, the geodesics are straight
lines. If we consider surfaces, which are of course two-dimensional
spaces, the geodesics are not necessarily straight lines. For example,
in the case of a spherical surface, the geodesics are the great circles,
that is, those circles on the sphere whose centers coincide with the
center of the sphere.
Let X and Y be two points. The distance L between them, measured
along some curve, is given by the line integral
y
(73.3)
L =
where
(73.4
According to the Calculus of Variations, L is an extremum if the path
joining X and Y is such that
(73.5)
d
ds
(dW\ dW
dP'l- d~ = o.
dW
m n
dp' = gmn Il,'P
+ gmnpm
Il,n
= grnpn
=
171
d ({}w) 2 dp 2 {}gm m.
{}p' = grn ds + dz m P 'P ,
{}w = {}gm. pmp.
ds
{}z'
{}z'
+ {}grn
pm'P. _.1 {}gmn pm'P. =
{}Zm
2 (}Z'
g dp
rn
ds
{}gm
{}zm
pmp. =
{}grm P"P m
(}z".
Substituting for one half of the middle term in Equation (73.6) from
Equation (73.7), we then obtain the relation
(73.8)
gm
dll
ds + [mn, r]pmp- = 0,
where
(73.9)
[mn, r] =
The quantities [mn,r] are the Christoffel symbols of the first kind. We
note that
[mn,r] = [nm,r].
(73.10)
Ifwe multiply Equation (73.8) by g't, we obtain
(73.11)
where
(73.12)
F!.,. = tTmn,r].
The quantities F~. are the Christoffel symbols of the second kind.
We note that
(73.13)
The first term in Equation (73.11) reduces to
t dpdpt
d' z'
8.-=- =-.
ds
ds
ds'
172
(74.1)
and
' IS
d~
ds'
(74.2)
d d'.q
+ F'spqds
~ -~ds
Ip
- 0
-.
But
dz'
ds
dz' dz 'P
dz'PTs'
d'z'
dZ' d'z'P
ds' = dz'P ds'-
d'Z' dz 'P d/ q
dz'P dZ'q dS Ts
dZ' d'z'P
dz'P ds'
d'Z'
dz'P dz ,q
dz m dZ) dz 'P dz,q
+ F'",. dZ'P dz ,q dS dS =
o.
't _
dz's dz m dz,
F pq - dz' d/P dz,qF m
dz's
+ dz'
d'z"
dz'P dz'q
~Ap/~u
Now
(75.3)
dA'p _ ~ ({}z, A)
du - du {}z'P ,
{}i dA,
= {}z'P du
+ {}z'P
Also, because of Equation (74.3), the second term on the right side of
Equation (75.2) satisfies the relation
(75.4)
Substitution from Equations (75.3) and (75.4) in (75.2) then yields
M'p = {};(" (dAm _ F' A dZn )
~u
{}z'P du
mn, du
{}zm Mm
= {}z'P
au'
~A, = 0
~u
3Ar _ dAr
3u - du
+P
m.
Am d::.
du
r.
3u -
Let Am", A~. and Am. be any second order tensors. Their absolute
derivatives along a curve Care
(75.7)
175
( 75.8)
(75.9)
q
+ pq n du - FPnq Am. l'
q
aAmn _ dAmn _ FP A dz _ P A
au - du
mq pn du
nq mp
OA~'n _ dA'!'n
au - ~
Fm AP dz
q
dz
du'
q
dz
du
aA
dA
au = du-
(75.lO)
ag
~
au
= 0'au amn -- 0 ,
agmn
"'u
0
o.
The proofs of these are left as exercises for the reader (Problem 28 at
the end of this chapter).
The rule for the absolute derivative of a product oftwo tensors is the
same as for the ordinary derivative of a product. Thus, for example
(75.12)
!au (ATs B) =
I
aA:s B
au
I
+ AT
's
aB I.
OU
Ap"Iq = {}AP
{}zq
(76 2)
(76.3)
A p'qrls --
+ FP,q A' ,
~
AP'qr + FPts At'qr - Ftqs AP'lr - Ftrs APqt
()ZS
(76.5)
gmnlq
=0,
O~q
=0,
g~~lq
=0.
Because of these relations, the operation of raising and lowering suffixes can be permuted with the operation of taking the covariant
derivative. Thus, for example, we have
Now
A"m
= ~:~-FfmAp.
177
(77.4)
Wrmn =
~m F:" -
;zrF:"
+ P;"Fpm - Ff.,Fpn
From Equation (77.3) it follows by the tests for tensor character that
rmn has the tensor character indicated by its suffixes. It is called the
mixed curvature tensor.
The covariant curvature tensor is
R:
(77.5 )
This tensor plays an important role in mathematical physics. For
any coordinate system in three dimensional space, this tensor vanishes,
since it vanishes for rectangular cartesian coordinates. On the other
hand, if for anyone coordinate system on a surface this tensor vanishes,
then there exists for this surface a curvilinear coordinate system such
that
g" = 1 if r = s
= 0 if r = s.
78. Cartesian tensors. A tensor is said to be cartesian when the transformations involved are from one set of rectangular cartesian coordinates to another. In 47 the special case of such transformations was
considered when the transformation is a rotation about the origin.
In the general case, the transformation is a rotation plus a translation.
If we superimpose a translation on the rotation considered in 47,
the equations of transformation are
(78.1 )
178
where ars are the constants considered in 47 satisfying the orthogonality conditions
(78.2)
Also, as and a's are constants such that a's = - aTS aT' Just as in 47,
we have the relations
(78.3)
dX's
dXT
dx' = aST = dX's'
\~;,J =
laSTI
soI=l.
Theorem 1. For cartesian tensors there is no distinction between
contravariant and covariant character.
Proof. Because of Equation (78.3) it follows that the laws of transformation of contravariant components and covariant components
are the same. Further
so that the raising or lowering of a suffix does not change the values
of the components.
Because of this theorem, when dealing with cartesian tensors we
do not need both superscripts and subscripts, so subscripts will be used
exclusively.
In 6 we introduced the orthogonal projections of a vector on rectangular cartesian coordinate axes, calling these projections the components of a vector. Throughout the rest of this book we shall refer
to these as the physical components of a vector for rectangular cartesian
coordinates.
Theorem 2, The physical components of a vector for rectangular
cartesian coordinates constitute a cartesian tensor of the first order.
179
Proof. The nine constants ar, in Equations (78.1) are the cosines
of the angles between the axes of two rectangular cartesian coordinate
systems. Hence, if b is a vector with physical components br and b'r
for these two systems of coordinates, then just as in 47 we have
b/
= asr br
= uX
i!:~ br.
r
79. Oriented cartesian tensors. Equation (78.4) shows that the Jacobian I of a transformation from one set of rectangular cartesian coordinates to another is equal to-either plus one or minus one. The former
case arises when the transformation is between coordinates whose
axes have the same ori~ntation (both right-handed or both lefthanded), and the latter case arises when the orientations are opposite.
A set of quantities is said to constitute an oriented cartesian tensor if
it is a cartesian tensor when I = 1 and is not a cartesian tensor when
1= -1.
In two-dimensional problems, suffixes have the range 1, 2. For
such problems we introduce a permutation symbol c" defined as follows:
(79.1)
cl l
= c,. = 0,
eu = 1,
Cn
= -1 .
Crs'
Crsl
(79.2)
C312
= 1,
C321
-1,
CU3
= 0.
In Theorem 4 of the next section we shall see that both of these permutation symbols are oriented cartesian tensors.
We can now express in tensor notation many of the formulas and
equations of the earlier chapters of this book. Thus, the scalar and
vector products of two vectors a and b are respectively
a, b" Crsl as b,
We note that this vector product is an oriented cartesian tensor. The
scalar triple product a (b X c) is
The expressioas \l j, \l. band \l X b become
dj
dX, '
db,
dXr '
Crsl
db ,
dXs'
181
"
A" = ~
AS I
{}~
{}z"
{}Z'
{}zP
Am
IW
O'P
g, = !g, ,s! =
{}Z'
{}z
{}z" {}z,;
g,. 1
1;;:,11 ;~I!gvw!
J2g,
Theorem 2. The covariant permutation symbol is a relative covariant tensor of weight -1.
Proof. We must show that
(80.2)
erst =
{}zm
{};C' {}zP
-1
Let us denote the right side of this equation by 1-1 <I> rst. Then
(80.3)
But the expression on the right side here is just the expansion of the
determinant I, so <1>123 = 1. In a similar manner, we see that each
component of <l>rst is equal to a determinant. If the suffixes on any
one component are subjected to a single permutation, two rows or
columns of the corresponding determinant are interchanged, and so
the sign of the component is changed. We thus have the results:
(i) <l>rst = 0 if two suffixes are equal, (ii) <l>rst = 1 if (rst) is an even
permutation of (123), (iii) <l>rst = -lif (rst) is an odd permutation of
(123). Thus we may write
.
and so the proof is complete.
'lJTSt
vg CTSt ,
- tW ( g-tWA"')
A
"'IT -g
... IT'
We note that A::: IT has the same weight as A:::. The absolute derivatives of relative tensors are defined analogously. It can be proved
that both the covariant and absolute derivatives of relative tensors
obey the same product rule as do ordinary derivatives.
(81.1 )
BT
dz:: b
dX s S,
B _ dXs b
T -
dz::
= VgTS BT B S
B = vamnbmb n
184
= vbmb m = b.
aTS
IS
J(,
then by 72 we have
cos 6 = g"
B'
A' .
If A(;) is a unit vector in the direction of the parametric line of Zl, then
(81.3)
A(I)
dz1
= ds = h/
Thus
185
(81.5)
There is a similar procedure for tensors of higher order. For example, let tTS be a set of nine quantities which are the components of a
cartesian tensor. The stress components of elasticity form such a set.
We define the tensorial components of this tensor for curvilinear coordinates z' in terms of trs by the appropriate laws of tensor transformation. Let us denote these components by the symbols Trs, T:"
T;s and T TS. The physical component of this tensor along the directions of two unit vectors ),: and [.tT is defined to be TT'ATfLs. The
physical components of this tensor for the curvilinear coordinates z'
are defined to be the physical components obtained by taking ),: and
[.t' in the directions of the parametric lines of the coordinates. Let us
denote these nine physical components by T(TS). Then
T(TS) = ymnA(T)mA(S)n,
AIs) being the three unit vectors in the directions of the parametric
lines of the coordinates. When the coordinates z' are orthogonal, then
\~) is as given in Equation (81.3), and there are two similar relations
for A(;) andA(;). We then get for T(TS) the expressions
h:T",
hlh, T 12 ,
hlh3 T 13 ,
h3hl T 31 ,
h3h, T 32 ,
h: T 33
(81.6)
There are other expressions for T(TS) in terms of the covariant components TTS' as well as in terms of mixed components.
Physical components of tensors of higher order can be defined
analogously.
186
(82 . 1)
v,
dx,
= dt '
a,
dv,
= dt .
IS
rna, =j,.
Let l be curvilinear coordinates of the particle, V', A' and F' being
contravariant component of velocity, acceleration and force for this
coordinate system. We have by definition
(82.3)
dZ'
dZ' dx,
dz'
= dX, v, = ax, dt = dt
+ F' vm dz n
mn
dt
We could define F' in terms of j" but it is easier to use the relation
(82.5 )
dW
= j,dx, = F,dz' ,
rnA' =F'.
To check these we note that they are tensor equations, and are true
when the coordinates are rectangular cartesian since in this case they
reduce to (82.2). Hence they are true for all coordinate systems.
The mathematical theory of elasticity. Using rectangular cartesian
coordinates x" we have the displacement u" the strain components
ers and the stress components Trs. For the determination of these,
we have the equations
(82.7)
(82.8)
187
(82.9)
JT"
Jx-
(82.lO)
+X
r
pddt
-U '
-
8=e =dUr
rr
dX r
(82.12)
(82.13)
ers
(82.14)
ernlsm +esmlrn
(82.15)
(82.16)
(82.17)
T rs
t(ur!s+uslr)'
ermlsn + esn1rm ,
k,smnemn ,
's
T .. !s
+X =
T
(j2U'
P Jt 2
'
(82.18)
These equations are tensor equations, and reduce to Equations
(82.7)-(82.12) when the coordinates are rectangular cartesian coordinates. Hence Equations (82.13)-(82.18) are true for all coordinate
systems, and are the desired equations. (Note that the term JV/dt"
188
dp
dt
(82.19)
+ dXd r (pur ) =
0,
dV r + dV r Vs = Xr _! dP ,
dt
dx,
p dXr
where p is the density, t is the time, Vr is the velocity, Xr is the external
force per unit volume, and p is the pressure. We convert these to
curvilinear coordinates in a manner analogous to that used for elasticity, obtaining the equations
(82.20)
~~ + (PVr)lr =
(82.21 )
(82.22)
du + r ,
Jt
VI'U =
xr
0,
I rs dP
-pg
dZ:
In a similar manner, we can obtain in terms of curvilinear coordinates the fundamental equations of other branches of mathematical
physics, such as electricity and magnetism, geometrical optics and
heat conduction.
Problems
1. Prove that
(Ars+Asr)lZ: = 2A rslz'.
2. If A = ArsZrz:, where Ars are constants, prove that
8: Art = A,t,
8:8~A~
A~ .
d 2Z"
dZ" dZ'v dZ'w d'Z"
d;:! dl = - dZ" d;:! dZ t dZ'v dZ'w
6. Let b, denote the covariant components of a vector for rectangular cartesian coordinates. Find the covariant components of this
vector for cylindrical coordinates r, e, Z in terms of b, and r , e, z.
7. If z' are curvilinear coordinates and cD is an invariant, do the
.
d2i J > .
expressIOns dZ'd;:! constItute a tensorr
'l
B. Prove that the sum of two tensors A~t and B;t is a tensor.
9. Prove that the outer product of two tensors
and B;t is a tensor.
10. If A;t is a tensor, prove that A;t is a tensor.
11. In three-dimensional space, how many different expressions are
represented by A:B;qC;'? When each such expression is written out
explicitly, how many terms does it contain?
12. Prove that the tensor A;' is transitive.
13. Let A~ be a set of nine quantities such that A~X: is an invariant,
where
is an arbitrary mixed tensor of the second order. Prove that
A; is a mixed tensor of the second order.
14. Let A"t be a set of quantities such that ATSt;rt is a covariant
tensor of the second order, where X t is an arbitrary contravariant
vector. Establish the tensor character of A"t.
15. If g" = 0 for = r s, prove that
1
1
3.
1
gIl = _,
g22 = _,
A:
X:
=-,
g
g23 = (1 =g12 = O.
16. Find the components of g" for cylindrical coordinates r, e, z.
17. Find the components of g" for spherical polar coordinates
gn
r,
e, 'P'
g2.
190
+ [rn, m] =
[Tm, n]
t,:., [mn, r]
g"F~.
T,
[ pq, s]'
dZ m dZ' dZ T
dZ'P dZ'q dZ'S [mn, T]
d 2 Zm
dZ n
Sgm. = 0
S
Su
'Su
sm
n
= ,
SU
.S
SA ~s BAr SBt
= s;;t +
.s &;'.
A~I' =
Jg ~, (Vi A').
n'f
v
1 d (. In
Vg
d,t v g
's dJ)
g d,t .
191
'\7'1 in
~: = 0,
tensor.
40. If Ar is a relative tensor of weight one, prove that
dAr
A: 1r = (Jzr'
41. If Br are the covariant components of a vector for orthogonal
curvilinear coordinates zr whose metric form is given in Equation
(81.2), show that the physical components of this vector for the curvilinear coordinates ;( are as given in (81.5).
42. If T' are the components of a tensor for the coordinates ;( in
Problem 41, find the physical components of T' for these coordinates
in terms of mixed components.
43. Let r, e and <p be the spherical polar coordinates of a particle.
Prove that for these coordinates the contravariant, covariant and
physical components of velocity are respectively
192
dt
dt
dt
JE
= k-,--,
dt
"V
JH
-[L- ,
dt
\JH=O,
\JE = 0,
where H is the magnetic intensity, k is the dielectric constant, E is the
electric intensity, t.J. is the permeability and t is the time. Express these
equations in tensor form for general curvilinear coordinates i.
193
omaNear
46pp.
$6.00
edited by
''Passion,''
Campbell
Pa. $7.50
QA433.H33 c.4
HK POLYTECHNIC LIBRARY
QA
Hay, G. E.
433
.1:.33
c.4
Vector
and
tensor
analaysis.
oJI.eGtiou