Cu 31924001246341
Cu 31924001246341
CORNELL
UNIVERSITY
LIBRARY
GIFT OF
^ale University
MATHBMAtlCS ;
V Date Due
- NOV 2 1 1H43 HI y
m rr
iif29''«!
\ 5 1983
issr
jggr29_i!^
i^uu
^ww<- ^m:r
-
FEB 2 2 1952
^(R ^
NOV
Hk f
p.
^
8 1994
NOV 4 1957
'OCT 2
i
TOF 063
JAN 6 1565
0£CUJii^er
h
JUL29
_iite^«a
Rir!H977
;LlZi.
FEB 19 979
MAR 19 19-5
'
••j
^[v '
#&-
341
3 1924 001 246
The original of tiiis book is in
http://www.archive.org/cletails/cu31924001246341
pale '^icmtmnial ^hlicatim$
VECTOR ANALYSIS
pale 'Bicentennial publications!
BY
By Yale University.
Published, December, igoi.
K.lSUljoc^
the public.
As, however, the years passed without my finding the
leisure to meet
this want, which seemed a real one, I was
very glad to have one of the hearers of my course on Vector
Analysis in the year 1899-1900 imdertake the preparation of
a text-book on the subject.
I have not desired that Dr. Wilson should aim simply
at the reproduction of my lectures, but rather that he should
use his own judgment in all respects for the production of a
text-book in which the subject should be so illustrated by an
adequate number of examples as to meet the wants of stu-
dents of geometry and physics.
J. WILLARD GIBBS.
GENERAL PREFACE xi
Page
Pkeface by Professor Gibbs vii
General Preface ix
CHAPTER I
CHAPTEK II
CHAPTEE ni
THE DIFFEEENTIAL CALCULUS OE VECTORS
55-56 Derivatives and differentials op vector functions
with respect to a scalar variable 115
57 Curvature and torsion of gauche curves .... 120
58-59 Kinematics op a particle. The hodograph . .. 125
60 The instantaneous axis op rotation 131
61 Integration with applications to kinematics . 133
. .
74-76
73 Laws OF operation of V; V* > VX
The partial application of V- Expansion op a vec-
157
CHAPTER IV
THE INTEGRAL CALCULUS OE VECTORS
Arts. Paqb
79-80 Line integrals of vector functioks with applica-
tions 179
81 Gauss's Theorem 184
82 Stokes's Theorem 187
83 Converse of Stokes's Theorem with applications . 193
84 Transformations of line, surface, and volume in-
tegrals. Green's Theorem 197
85 Remarks on Multiple-valued functions 200
86-87 Potential. The integrating operator " Pot " . . 205
88 Commutative property op Pot and ^ 211
89 Remarks upon the foregoing 215
90 The integrating operators "New," "Lap," "Max" 222
91 Relations between the integrating and differ-
entiating operators 228
92 The potential " Pot " is a solution of Poisson's
Equation 230
93-94 Solenoidal and irrotational parts of a vector
function. Certain operators and their inverse . 234
95 Mutual Potentials, Newtonians, Laplacians, and
Maxwellians 240
96 Certain boundary value theorems 243
Summary of chapter iv 249
Exercises on chapter iv 265
CHAPTER V
LINEAR VECTOR FUNCTIONS
97-98 Linear vector functions defined 260
Dyadics defined
99 264
100 Any linear vector function may be represented
BY A dyadic. Properties of dyadics .... 266
101 The nonion form of a dyadic 269
102 The dyad or indeterminate product of two vec-
tors IS the most general. Functional property
of the scalar and vector products 271
103-104 Products of dyadics 276
105-i07 Degrees of nullity of dyadics 282
108 The idemfactor 288
Akts.
YECTOE ANALYSIS
VECTOR ANALYSIS
CHAPTER I
ues which are wholly determined when the unit of the scale
is fixed. On the other hand the remaining quantities are not
capable of such representation. Force to be sure is said to be
of so many pounds or grams weight; velocity, of so many
feet or centimeters per second. But in addition to this each
of them must be considered as having direction as well as
magnitude. A force points North, South, East, West, up,
down, or in some intermediate direction. The same is true
of displacement, velocity, and acceleration. No scale of num-
bers can represent them adequately. It can represent only
their magnitude, not their direction.
2 VECTOR ANALYSIS
The positive and negative numbers of ordinary/ algebra are the
typical scalars. For this reason the ordinary algebra is called
scalar algebra when necessary to distinguish it from the
vector
points move in straight lines in the same direction and for the
same distance D. This is equivalent to shifting space as a
rigid body in that direction through the distance D without
rotation. Such a displacement is called a translation. It
must be scalar ; m, the mass, and the timCj are also t.,
the point P falls upon the point and P' upon some point T.
Then .
A = PP' =
,
OT =
.
'S.
T~{x',y',z').
When two vectors are equal the three scalars which repre-
sent them must be equal respectively each to each. Hence
one vector equality implies three scalar equalities.
-
6 VECTOR ANALYSIS
Definition : A vector A is said to be equal to zero when its
magnitude A is zero.
A= if ^ = 0.
All null vectors are regarded as equal to each other without
any considerations of direction.
Scalar Multiplication
X {y A) = (^xy') A= y (x A)
a ^— A=— •
A A
8 VECTOR ANALYSIS
7.] Definition : The negative sign, — , prefixed to a vector
reverses its direction but leaves its magnitude unchanged.
For example if A be a displacement for two feet to the right,
—A is a displacement for two feet to the left. Again if the
strokeA £ he A, the stroke S A, which is of the same length
as A B but which is in the direction from B to A instead of
from A to B, will be — A. Another illustration of the use
of the negative sign may be taken from Newton's third law
of motion. If A denote an "action," — A wiU denote the
" reaction." The positive sign, + may, be prefixed to a vec-
tor to call particular attention to the fact that the direction
has not been reversed. The two signs + and — when used
in connection with scalar multiplication of vectors follow the
same laws of operation as in ordinary algebra. These are
symbolically
+ + = +; +- = -; - + --; -- = + ;
— (m A) =m (— A).
8.] The
addition of two vectors or strokes may be treated
most simply by regarding them as defining translations in
space (Art. 2). Let S be one vector and T the other. Let P
be a point of space (Fig. 2). The trans-
lation^ carries P into P' such that the
line PP' is equal to S in magnitude and
direction. The transformation T will then
carry P' into P"— the line IFW being
parallel to T and equal to it in magnitude.
Consequently the result of S followed by
T is to carry the point P into the point
P". If now Q be any other point in space, S will carry
Q
into Q' such that CS' = S and T will then carry Q'
into Q"
ADDITION AND SCALAR MULTIPLICATION 9
That is
R= S -I- T.
by writing
K=S+T=T+ S.
note that in case the vectors lie along the same line vector
addition becomes equivalent to algebraic scalar addition. The
lengths of the two vectors to be added are added if the vectors
have the same direction ; but subtracted if they have oppo-
site directions. In either case the sum has the same direction
as that of the greater vector.
10.] After the definition of the sum of two vectors has
been laid down, the sum of several may be found by adding
together the first two, to this sum the third, to this the fourth,
and so on until all the vectors have been combined into a sin-
ADDITION AND SCALAR MULTIPLICATION 11
gle one. The final result is the same as that obtained by placing
^ ' the origin of each succeeding vector upon the terminus of the
preceding one and then drawing at once the vector from
the origin of the first to the terminus of the last. In case
'
these two points coincide the vectors form a closed polygon
I* and their sum is zero. Interpreted geometrically this states
that if a number of displacements R, S, T are such that the • •
A - B = A + (- B).
By this means subtraction is reduced to addition and needs
12 VECTOR ANALYSIS
no special consideration. however an interesting and
There is
ence of two vectors which are drawn from the same origin is
the vector drawn from the terminus of the vector to be sub-
-tracted to the terminus of the vector from which it is sub-
tracted. Thus the two diagonals of the parallelogram, which
is constructed upon A and B as sides, give the sum and dif-
ference of A and B.
II A B = B + A,
-1-
III,, (m + n) A = m A + nA,
III, m (A -I- B) —mA + mB,
IIL - (A + B) = - A - B.
:
3 A + 4B = E
and 2 A+ 3 B =F
yield by the usual processes the solutions
A=3E-4F
and B = 3 F - 2 E.
Components of Vectors
Tj = ajj +a 2/i + Si c,
b
Ta = ajj +a 2/2 + c,
^ «2!2
Tg = ajg a + 2/3 b + Sg c,
«! + OJg + iTg + • = V
• •
yi + + 2^8 +
2^2
• = ^ r = 0.
• •
(5)
zi + z^ + z^+ ... = 0^
Should the vectors all be coplanar with a and b, all the com-
ponents parallel to c vanish. In this case therefore the above
equations reduce to those given before.
16.] If two equal vectors are expressed in terms of the
same three non-coplanar vectors, the corresponding scalar co-
efficients are equal.
:
18 VECTOR ANALYSIS
Let r = r',
Z Z
Right-handed Left-handed
Fig. 8.
of the tbnmb, first finger, and second finger of the right hand. If the thumb and
first finger be stretched ont from the palm perpendicular to each other, and if the
second finger be bent over toward the palm at right angles to first finger, a right-
handed system is formed by the fingers taken in the order thumb, first finger,
second finger.
20 VECTOR ANALYSIS
the Z^-plane on which rotation from the F-axis to the Z-
axis is counterclockwise. Thus it appears that the relation
between the three axes is perfectly symmetrical so long as the
same cyclic order XYZXY is observed. If a right-handed
a; i,
y j, z k.
Applications
R= i (S + T).
Ii = AB = AE+ EB = LT + x(S-tT:)
Hence ^ T + a; (S -i T) = y (S -|- T)
or a;S-l-i(l-a;)T = yS-f-2/T.
Hence, equating corresponding coefficients (Art. 16),
x = y,
-Q.-x)=y.
ADDITION AND SCALAR MULTIPLICATION 23
From which y
"
= \
3
.
E = (m + w) S +m (T - S}.
Hence (m + w) S is the fraction oi AB which is cut off by the
line through B parallel to B 0. Consequently by similar tri-
(1 - m) -1- (1 - w) -I- (m + ») = 2,
and the theorem is proved.
Example 3: If from any point within a parallelogram lines
E= AB = m S + ra T,
Equating coefficients,
By solution, =»
•^
m+n—1
m
y =
m+n—1
Substituting either of these solutions in the expression for P,
the result is
expressions for the same vector are usually found. The co-
efficients of the corresponding terms are equated. In this way
the equations between three unknown scalars are obtained
from which those scalars may be determined by solution and
then substituted in either of the expressions for the required
vector. The vector method has the same degree of adapta-
bility as the Cartesian method in which oblique axes with
different scales are employed. The following examples like
those in the foregoing section are worked out not so much for
their intrinsic value as for gaining a familiarity with vectors.
Example 1 : Let A B CD be a tetrahedron and P any
point within it. Join the vertices to P and produce the lines
until they intersect the opposite faces in A\ B', C, D'. To
show
PA' PB' PC PD'
DD' ~
"* *"
AA''^ 'BB' GC '
= AP=l'&-\- mC + n'O,
"B
A' = ^ ^' = P = Ai G B m C 4-
;&! -I- 71 D).
if
Hence Ai =
I +m+n
PA'
and -T—p = ^1-1
-^7— =1-
, ,,
(l + m ^+
,
,
7i).
26 VECTOR ANALYSIS
Vi = *2 »•
Hence k^ =
PB' k,-l
BB' = -j;-=^-'
^^•^
PC = mand^^
,PD'
= ».
^^
Adding the four ratios the result is
~''
Then ,-*.,>^t*nr^tt=:2r ?- ^fc ^^ ^
= z[C + y (lA + mB+
nO — C)].
Hence 1 — x = zyl,
X = zym,
= z[l + y(n-l-)l
Hence
1 +
28 VECTOR ANALYSIS
The components of P parallel to A and B are in inverse ratio
and
l-x=\y.
1
ADDITION AND SCALAR MULTIPLICATION 29
nA-\-mB
P=
m+n
and M=\(^A + B+ C).
aA-t-5B-|-cC-|-c?D-| = 0.
30 VECTOR ANALYSIS
Change the origin from to C by adding a constant vector
E = OW' to each of the vectors A, B, C, D The equation
then becomes
a (A + E) + 5 (B + E) + c (C + E) + d (D + R) + = • • •
a +h+ + d-\ = 0.
If P = nA + m'B
m+n
= — m
n
1 +
m+n m+n
If M = J (A + B + C),
'• ~ 3
•"
3 •" 3 •
First let aA + 6B =
and a + b = 0.
First let aA + 6B + cC =
and a + 6 + c = 0.
(B - A) = a; (C - A)
subsists. The sum of the coefficients on the two sides is
the same.
The necessary and sufficient condition that an equation,
in which the sum of the scalar coefficients is zero, subsist
' Yectois which have a common origin and teiminate in one line aie called hj
Hamilton " termino-collinear."
32 VECTOR ANALYSIS
between four vectors, is that if drawn from a common origin
they terminate in one plane.^
^ (B - A) + m (C - A) + 71 (D - A) = 0,
where I, m, and n are certain scalars. The sum of the coeffi-
cients in this equation is zero.
Between any five vectors there exists one equation the sum
of whose coefficients is zero.
Let A, B, C, D, £ be the five given vectors. Form the
differences
E-A, E-B, E-C, E-D.
One of these may be expressed in terms of the other three
— or what amounts to the same thing there must exist an
equation between them.
/fc (E - A) + Z (E - B) + m (E - C) + w (E - D) = 0.
The sum of the coefficients of this equation is zero.
I Vectors which
have a common origin and terminate in one plane are called
by Hamilton " termirwxomplanaT"
ADDITION AND SCALAR MULTIPLICATION 33
«i
—
»
34 VECTOR ANALYSIS
E + wB = D + wC = (l + w)G.
n(B + C^ = (l + n)6-(l- m) A.
This vector cuts BC and AG. It must therefore be a
multiple of F and such a multiple that the sum of the coeffi-
Hence w (B + C) = (1 + m) G - (1 - w) A = 2 wF.
Hence F = B—+-—C
and the theorem has been proved. The proof has covered
considerable space because each detail of the reasoning has
been given. In reality, however, the actual analysis has con-
sisted of just four equations obtained simply from the first.
aA + 6B+^P =
and a + h + p = 0.
Therefore P = ;
a + b
ADDITION AND SCALAR MULTIPLICATION 35
P = xA + yB
x + y
P = xA + yB + zC
x + y + s
BD intersect in G. To show
that FG intersects AB ina, point F' and CD in a point F"
such that the lines AB and CD are divided internally at
F' and F" in the same ratio as they are divided externally
by F, That is to show that the cross ratios
36 VECTOR ANALYSIS
Choose the origin at random. The four vectors A, B, C, D
drawn from it to the points A, B, C, D terminate in one
plane. Hence
aA + 6B + cC + dD =
and a-t-6 + c + = 0.
<^
aA + cC = - (&B + dD),
a + c = — ip + d).
aA + cC 6B + c?D
Divide G=
a + c h + d
aA + dTi 6B + cC
In like manner P=
a + d h + c
(a + c)G — (a + d)F _ cC — dH
Form
(a + c) — (a + d) {a + c) — {a + d)
{a + c)Q,-{a + d)-E
= cG ^ dJi
or
— c
-z
d — ^— = Ji c a
. (.aj
aA + 6B cC + dH =
E. (6)
a + b c + d
Hence E divides AB
in the ratio a : 6 and CD in the ratio
c .• d. But equation (a) shows that H" divides CD in the
ratio — c:d. Hence B and B" divide CD internally and
externally in the same ratio. Which of the two divisions is
internal and which external depends upon the relative signs
of c and d. If they have the same sign the internal point
and a + h + c + d=
possible to obtain
£ = aA
—
+ hB
= =
cC-iLdJ)
*-—r— J
a + +d h , c
F=
aA+ cC SB + dD
a + c b + d
aA + dD 6B + cC
G= — =— e ; J
a + d +
by splitting the equations into two parts and dividing, Next
four vectors such as A, D, E, F may be chosen and the equa-
tion the sum of whose coefficients is zero may be determined.
This would be
— aA + d'D + (a + 6)E+ {a + c)F = 0.
By treating this equation as (c) was treated new points may
be obtained.
38 VECTOR ANALYSIS
_-ak + dJi _ (g + 6)E + (a + c)E
~
J a +c+ ei
Centers of Gravity
a A + 6B
a + b
Then Q = (a +
^ b) "^^^f^ + e C
'^
a + b
a +b + c
Hence G = aA + bB + ,
cC
a +
,
b +
,
c
40 VECTOR ANALYSIS
Evidently the center of gravity of any number of masses
a, h, c, d, ... situated at the points A, B, C, D, ... may
be found in a similar manner. The result is
— a + i + c + d+ ' ^ ^
GBC:GOA:GAB:ABO=a:h:c:a + h + c.
GBO
In a similar manner
BOA _ a J + c -I-
WCA I
and
CAB _ g + 5 + c
6AB~ c
BCDQ:CDAG:DABG:ABCG:ABCD
= a :i : c : d:a + h +c+ d.
In like manner
and
and
;
42 VECTOR ANALYSIS
a, J, c may therefore be looked upon as coordinates of the
points P inside of the triangle ABG. To each set there
corresponds a definite point P, and to each point P there
corresponds an infinite number of sets of quantities, which
however do not differ from one another except for a factor
of proportionality.
To obtain the points P of the plane ABC which lie outside,
of the triangle ABG one may resort to the conception of
negative weights or masses. The center of gravity of the
masses 2 and — 1 situated at the points A and B respectively
would be a point G dividing the line A B externally in the
ratio 1 : 2. That is
GA:GB = 1:2.
Any point of the line AB produced may be represented by
a suitable set of masses a, h which differ in sign. Similarly
any point P of the plane ABG may be represented by a
suitable set of masses a, b, c of which one will differ in sign
from the other two if the point P lies outside of the triangle
ABG. Inasmuch as only the ratios of a, h, and c are im-
portant two of the quantities may always be taken positive.
The idea of employing the masses situated at the vertices
as coordinates of the center of gravity is due to Msbius and
was published by him in his book entitled '' Barycentrische
GalciU" in 1826. This may
be fairly regarded as the starting"^
point of modern analytic geometry.
The conception of negative masses which have no existence
in nature may be avoided by replacing the masses at the
vertices by the areas of the triangles GBG, GGA, and
GAB which they are proportional. The coordinates of
to
a point P would then be three numbers proportional to the
areas of the three triangles of which P is the common vertex
and the sides of a given triangle ABG, the bases. The sign
of these areas is determined by the following definition.
ADDITION AND SCALAR MULTIPLICATION 43
ABG=BGA = GAB.
Interchange of two letters which amounts to a reversal of
the cyclic order changes the sign.
FAB+PBC+FGA=ABG
must hold. The same will also hold if P be outside of the
triangle provided the signs of the areas be taken into con-
sideration. The areas or three quantities proportional to
them may be regarded as coordinates of the point P.
The extension of the idea of " iarycentric " coordinates to
space is immediate. The four points A, B, G, D situated at
the vertices of a tetrahedron are weighted with mass a, 6, c, d
respectively. The center of gravity G is represented by
these quantities or four others proportional to them. To
obtain points outside of the tetrahedron negative masses
may be employed. Or in the light of theorem 2, page 40,
the masses may be replaced by the four tetrahedra which
are proportional to them. Then the idea of negative vol-
umes takes the place of that of negative weights. As this
idea is of considerable importance later, a brief treatment of
it here may not be out of place.
Definition: The volume ABGD of a tetrahedron is said
to be positive when the triangle ABG appears positive to
44 VECTOR ANALYSIS
the eye situated at the point D. The volume is negative
ABGP-BG1)P+ GDAP-DABP=ABGD
holds good. It still is true if P be without the tetrahedron
provided the signs of the volumes be taken into considera-
tion. The equation may be put into a form more symmetri-
cal and more easily remembered by transposing all the terms
to one number. Then
ABGD + BGDP+ GDPA + DPAB + PABG=0.
The proportion in theorem 2, page 40, does not hold true
if the signs of the tetrahedra be regarded. It should read
BGDG:GDGA:DQAB: GABG:ABGD
=a ; 6 ; c ; d ; a + 6 +c+ d.
ADDITION AND SCALAR MULTIPLICATION 45
xA + yB + zC
x+ y + e
aA + &B + cC
a + b + c
x + y + z -{ w
which expresses any vector P drawn from an indeterminate
origin in terms of four given vectors A, B, C, D drawn from
the same origin, it may be seen by comparison with
G = aA + &B + c C + tZD.
a + b +c+d
that the four quantities x, y, z, w are precisely the bary-
centric coordinates of P, the terminus of P, with respect to
the tetrahedron A B CD. Thus the vector methods in which
the origin is undetermined and the methods of the " Bary-
centric Calculus " are practically co-extensive.
the negative ^-axis lies. For this reason the point of view
and the direction of description of the boundary must be kept
clearly in mind.
Another method of stating the definition is as follows : If
A' =A cos X,
4
50 VECTOR ANALYSIS
Theorem 2 : The vector which represents a closed polyhedral
surface is zero.
covered twice, once with positive area and once with negative.
These cancel each other. Hence the total projection of a
zero.
The theorem has been proved for the case in which the
closed surface consists of planes. In case that surface be
' Snch a state of affairs is realized to all practical purposes in the case of a
Summary of Chapter I
r — xi + yi + zTs.. (6)
m+n
The necessary and sufficient condition that a vector equation
represent a relation independent of the origin is that the sum
of the scalar coefficients in the equation be zero. Between
any four vectors there exists an equation with scalar coeffi-
^^aA + 5B + cC+..._
a + h + c-\- •• ^ '
Exercises on Chapter I
1 A complete quadrangle consists of the six straight lines which may be passed
through four points no three of which are collinear. The diagonals are the lines
which join the points of intersection of pairs of sides.
54 VECTOR ANALYSIS
point E and that this point is such that PH passes through
and is bisected by the center of gravity of the tetrahedron.
10. Show that without exception there exists one vector
equation with scalar coefficients between any four given
vectors A, B, C, D.
/ 11. Discuss the conditions imposed upon three, four, or
five vectors if they satisfy two equations the sum of the co-
efficients in each of which is zero.
CHAPTER II
the scalar product owing to the fact that its value is sca-
In case the two vectors A and B are collinear the angle be-
tween them becomes zero or one hundred and eighty degrees
and its cosine is there^^ce equal to unity with the positive or
negative sign. Hence the scalar product of two parallel
vectors is numerically equal to the product of their lengths.
The sign of the product is positive when the directions of the
vectors are the same, negative when they are opposite. The
product of a vector by itself is therefore equal to the square
of its length
K'L=A\ (3)
A B cos (A, B) = 0.
i • i =i •
J =k • k = 1, (4)
i . j = j . k =k . i = 0.
If more generally a and b are any two unit vectors the
product
a • b = cos (a, b).
DIRECT AND SKEW PRODUCTS OF VECTORS 57
A A
—— A = -j^
"R 7i
A a cos (A, H) — B cos (A, B) a, (5)
A* A -^ ^
where a is a unit vector in the direction of A If A is itself a
unit vector the formula reduces to
58 VECTOR ANALYSIS
(A + B) . C =A • C +B • C. (6)
tion. To show
(A + B) . (Ce) = A (Cc) + B. • (Cc)
or (A + B) .0 = A'C + B'C.
A • c is the projection of A upon c ; B • c, that of B upon c
+
The scalar product may be used just as the product in ordi-
nary algebra. It has no peculiar difficulties.
A = ^1 i + ^2 j + ^3 k,
and "&= B^i + B^i-ir B^]s.,
then A.B= {A^i + A^i + A^Ts.) (^Byi . +
B^i + B^Y)
= ^1 ^1 i. i -I- ^1 ^ji. j 4- ^^^3 i.k
+ ^2^1] .i^+ ^2-BJ.i -H ^2^8 j. k
+ ^3 5lk J^+ ^3 ,52 k j ^3 ^3 k k.
. . -1- .
cos (A, B) = cos (A, X) cos (B, X) + cos (A, Y) cos (B, Y)
+ cos (A,^ cos (B,^).
C.C-D.D = 4 A^B
or C^2 _ 2>2 =, 4 ^ ^ cos (A, B).
A = ^1 i + ^2 j + ^3 k,
then L'li. = A'^ = A^^ + A^ + A^. (8)
C = AxB
and read A cross B. For this reason it is often called the cross
the scalar product did in the case of the cosine. The two prod-
ucts are in a certain sense complementary. They have been
denoted by the two common signs of multiplication, the dot
and the cross. In vector analysis they occupy the place held
by the trigonometric functions of scalar analysis. They are
at the same time amenable to algebraic treatment, as will be
seen later. At present a few uses of the vector product may
be cited.
If A and B (Fig. 18) are the two adjacent sides of a parallel-
ogram the vector product
C =A X B = ^ 5 sin (A, B) c
although the sign is not required for the case of the cosine.
25).
ABG = -ACB.
In the same manner this reversal of sign, which occurs
when the order of the factors in a vector product is reversed,
will appear after a little practice and acquaintance just as
when expanding.
(A + B)xC = AxC + BxC. (11)
AX C, B X C, - (A + B) X C.
i (A X B) and - ^ (A X B).
But the sum of all the faces of the prism is zero; for the
prism is a closed surface. Hence
AxC + BxC-(A + B)xC + |(AxB)-i(AxB) = 0,
"'
A X C + B X C - (A + B) X C = 0,
or A X C + B X C = (A + B) X C. (11)
AX (C + D) + B X (C + D) = (A + B) X (C + D).
AxC+AxD+BxC+BxD
= (A + B) X C + (A + B) X D.
But by (11) A>(b + B X D = (A + B) X D.
Hence AxC + BxC=(A + B)xC.
This completes the demonstration. The distributive law holds
for a vector product. The generalization is immediate.
+ {A^B^-A^B,)ls..
1 This follows also from the fact that the sign is changed when the order of
factors is reversed. Hence j X j = — j X j = 0.
6
66 VECTOR ANALYSIS
1 J k
Ax B =
^1 B.
-"2 B,
-"3
The formulae for the sine and cosine of the sum or dif-
= cos + sin
a a; i a; j,
/I b = cos y + sin y 1 j,
= H + mj + mk,
a
a' = Z'i + m'j + 5i'k,
+ Qm' — I'm) k,
(A-B) C.
1 Later (Chap. V.) the product BC, where no sign either dot or cross occurs,
will be defined. But it will be seen there that (A«B) C and A« (B C) are identical
and consequently no ambiguity can arise from the omission of the parenthesis.
68 VECTOR ANALYSIS
37.] The second triple product is the scalar product of
two vectors, of which one is itself a vector product, as
A'(BxC) or (AxB).C.
is negative.
In case A, B, and C are coplanar this volume will be
neither positive nor negative but zero. And conversely if
BxC = — CxB.
Hence : A scalar triple product is not altered by interchanging
the dot or the cross or by permuting cyclicly the order of the
vectors, but it is reversed in sign if the cyclic order be changed.
38.] A word is necessary upon the subject of parentheses
in this triple product. Can they be omitted without am-
biguity ? They can. The expression
A'BxC
can have only the one interpretation
I
A.(BXC).
B = ^1 i + ^2 j + ^3 k,
then [A B C] = A, B, C, + B, C, A, + C, A, B,
-AB,C,-B,C,A,-C,A,B,. ^^^^
This may
be obtained by actually performing the multiplica-
tions which are indicated in the triple product. The result
may be written in the form of a determinant.^
A A ^3
[A B C] = B^ B, B, (18)' .
Ax(BxC) = a; B + ^ C,
(AxB)xC = mA + toB
in general
(AxB)xC is not equal to Ax (BxC).
The parentheses therefore cannot be removed or inter-
changed. It is essential to know which cross product is
72 VECTOR ANALYSIS
formed first and which second. This product is termed the
vector triple product in contrast to the scalar triple product.
The vector triple product may be used to express that com-
ponent of a vector B which is perpendicular to a given vector
A. This geometric use of the product is valuable not only in
itself but for the light it sheds
upon the properties of the product.
AXB
Let A (Fig. 21) be a given vector
axb'
and B another vector whose com-
ponents parallel and perpendicular
to A are to be found. Let the
components of B parallel and per-
pendicular to A
B and B " re-
be '
B'=_A (21)
B = B' + B'-^A-^^^>
A'A A«A
(22) '
This proves the formula in case two vectors are the same.
To prove it in general express A in terms of the three
non-coplanar vectors B, C, and BxC.
A = 6B + cC + a (BxC), (I)
+ a (BxC)x(BxC).
(BxC)x(BxC) =
Ax(BxC) = 6Bx(BxC) + c Cx(BxC). (II)'
Multiply by B.
-A.[Bx(BxC)]=-A.[B.C B-B.BC]
= -B.C A.B + B.B A.C.
Hence n = l and Ax(BxC) = A.C B - A.B C. (24)
(AxB).(CxD).
A-C A.D
(AxB).(CxD) = (25)'
B.C B.D
If A and D be called the extremes ; B and C the means ; A
and C the antecedents ; B and D the consequents in this
product according to the familiar usage in proportions, then
the expansion may be stated in words. The scalar product
of two vector products is equal to the (scalar) product of the
antecedents times the (scalar) product of the consequents
diminished by the (scalar) product of the means times the
(scalar) product of the extremes.
To reduce a vector product of two vectors each of which
is itself a vector product of two vectors, as ,'^
(AxB)x(CxD).
Let CxD = E. The product becomes ' ^ "'^ ^^ ^
' ,^-, 'l
[B C D] A- [C D A] B + [D A B] C - [A B C] D = 0. (27)
[A B ^D| = [B C D] a+ [C a D] B + [A B D] C. (27)'
A£, AG,. etc., give the angles between the vectors A and B,
A and C, etc. The points A,B,C,D determine a quadrilateral
upon the sphere. AG and £D are one
pair of opposite sides ; AD and B G, the
Hence
sin A B sin CD cos x = cos A C cos BD cos AD cos B G.
In words : The product of the cosines of two opposite sides
(AxB)x(AxC) = [ABC]A.
Let the three vectors a, b, c be unit vectors in the direction of
BxC, CxA, AxB respectively. Then
AxB = c sin c, AxC = — b sin 6
(AxB)x(AxC) = — cxb sin c sin 6 = A sin c sin b sin A
[A B C] = (AxB).C = cC sin c = cos (90° -p„) sin c
[A B C] A = sin c sin p^ •
80 VECTOR ANALYSIS
By equating the results and cancelling the common factor,
a^ = h^ + c^-2bccosA
&2=:c2 + a^-2cacoaB
c2 = a2 + 62 _ 2 a J cos C.
2axb=:2bxc = 2cxa =
2 a & sin C=26c sin ^ =^casm B.
If each of the last three equalities be divided by the product
\ ah c, the fundamental relation
Multiply by • b X c.
r.fexc = a a-bxc + b b'bxc + cc-bxc
or [r b e] = a [a b c].
[r c a] = & [b c a]
[r a b] = c [c a b]
Hence r = ^^^ a +
[a be] f^ b+ &^
[be a] [c a b]
c. (28)
bxc »
cxa >
axb
[abe] [abc] [abe]
DIRECT AND SKEW PRODUCTS OF VECTORS 83
a', b', c' which are the normals to the planes of b and c,
[a b c] r-a = [b c a],
a; x = rta
[abcjr.b = y [cab], y = r.b
[a b c] r»c = a [a b c], z = r»c
= b'-b = c'.c = l
a'.a (32)
a'«b = a'»c = b'»a = b'«c = c'»a = c'«b = 0.
a
,
• a = bxc • a = bxc-a = [be a]
=1
[a b c] [a b c ] [a b o]
,
a'.b
^
= bxc,b = bxcb =
. =0
be]
[a [a b cj [a b c]
and so forth.
Conversely if two sets of three vectors each, say A, B, C,
and a, b, c, satisfy the relations
, = B-b = C'e = 1
A.a
A.^ = A-e = B.a = B-c = C^a = C^b =
DIRECT AND SKEW PRODUCTS OF VECTORS 85
r '1,' n r ^
X ° cxa axb l
-*
"^"[[alJc] [abc] [abc]J
1 _„ 1
[abc] 2_
.
Hence [a'b'c'] = .
(33)'
[a b c] ^ ""
-^ [a b c]
P-A
(P.axR)(A.BxC) =
DIRECT AND SKEW PRODUCTS OF VECTORS 87
For tliis reason the primes i', j', k' are not needed to denote
a system of vectors reciprocal to i, j, k. The primes will
therefore be used in the future to denote another set of rect-
angular axes i, j, k ,
just as X', Y', Z' are used to denote a
set of axes different from X, Y, Z.
The only systems of three vectors which aretheir own reciprocals
A-B = A.C = 0.
B.A = B.C = 0.
Hence B is perpendicular to A and C.
C'A = C'B = 0.
Hence C is perpendicular to A and B.
Hence A, B, C must be a system like i, j, k or like i, j, — k.
* 46.] A scalar equation of the first degree in a vector r is
88 VECTOR ANALYSIS
where a, b, c, d, e, f are known vectors ; and a, h, c, d, known
scalars. Obviously any scalar equation of the first degree in
an unknown vector r may be reduced to the form
r«A =a
where A is a known vector ; and a, a known scalar. To ac-
r cos (r, A)
r«A =a
r.B = b (37)
r«C = c,
r-A =a
r.B =b (39)
r-C = c
r.D =d
A a-r + B b^r + C cr = D,
where A, B, C, D, a, b, c are known vectors. No scalar coeffi-
cients are written in the terms, for they may be incorporated in
the vectors. Multiply the equation successively by A', B', C
It is understood of course that A, B, C are non-coplanar.
DIRECT AND SKEW PRODUCTS OF VECTORS 91
a-r = D«A'
bT = D.B'
c.r = D«C'.
in this form.
ni = na.' a,'T + nh'^'T + n c'c»r
Applications to Mechanics
R = f, + f2 + f3 + ...
(41)
taken.
Mo {f}=2 0P<2.
f r f 2' • • •
•
Mo m = dxf (42)
For ix{ = df sin (d, f) e,
dxf = (isin(d,f)/e.
Now <?sin'(d, f) is the perpendicular distance from to f.
Mo {f i} = dxf 1
Mo {i^l = dxfa
Mo {f i} + Mo {f J + • • • = dxf + dxf 2 +
1 (43) • •
•
di = d/ + c, da = da' + c, • • •
R.Mo'{fi,f2,--}=R-Mo {fi.fa,---}
where 0' and are any two points in space. This important
relation follows immediately from the equation
R.Mo' ^Ro} =
and the relation is proved. The variation in the total
S6 VECTOR ANALYSIS
braces { }.
Then
RxMc = RxMo — Rx(cxR) = 0.
The problem is to solve this equation for c.
RxMo = R'R c
:=
RxMo
C •
R«R
If c be chosen equal to this vector the total moment about
the point 0', which is at a vector distance from equal to c,
the total moment and the resultant is constant and since the
resultant itself is constant it is clear that in the case where
they are parallel the numerical value of the total moment
will be a minimum.
The total moment is unchanged by displacing the point
about which it is taken in the direction of the resultant.
zero act upon a rigid body, then there exists in space one
and only one line such that the total moment about any
point of it is parallel to the resultant. This line is itself
V = axr. (45)
's — *3-'^'^8
DIRECT AND SKEW PRODUCTS OF VECTORS 99
where r^, r^, Tj, • • • are the radii vectores drawn from points
on the axis a^, ag, fig, • • • to the same point of the body. Let
the vectors r^, t^, fg, • • • be drawn from the common point of
intersection of the axes. Then
ri = r2 = r3 = ...=r
and
V = Vi + Vj + Vg + • • • = a^xr + sl^xt + agXr + • • •
= (ai + aa + ag + .)xr. . •
V = axr (45)
axr = — Vq.
This may be done as follows. Multiply by xa.
(axr)xa = — VoXa
or a-a r — a«r a = — VoXa.
^ ^- Yq X a
a>a
The point r, thus determined, has the property that its veloc-
Vo = Vo' + Vo"
D.fi, D.f2, . .
On the other hand the work done by f " is equal to the work
done by f during the displacement. For f' being parallel to
a is perpendicular to its line of action. If h be the common
vector perpendicular from the line a to the force f ", the work
done by f" during a rotation of angular velocity a for time
t is approximately
•
W = h f" a t = &'la.xt" t.
Mo'{fi,f2, •••}=0.
Hence a«Mo' {fp ^2> ' "I ^ =
where a any vector whatsoever. But this expression is
is
equal to the work done by the forces when the body is rotated
for a time t with an angular velocity a about the line a
passing through the point 0'. This work is zero.
Applications to Geometry
From this
The scalarsiXj, a^, a^; b^, b^ti^; Cj, o^, Cg are respectively the
direction cosines of i'; j'; k' with respect to i, j, k.
That is
6i
= cos (j', 1) &2 = cos (j',j) 63 = 003 (j', k) (48)
Cj — cos (k', 1) Cg = cos (k', j) Cg = cos (k', k).
]
j'-j' = 1 = Si^ + b^^ + 632 (49)
C k'.k' = 1 = c^a + c^2 + C32
^ i.i = 1 = aj2 + 2 + c^2 J,
i.j = = «! a^ + 61 ^2 + Cj C2
«1
and [ijk] = [i'j'k'] = l =
;
Let be the origin and B the vector 'OB. Let E be the ra-
dius vector from to any point of the required line. Then
E — B is parallel to A. Hence the vector product vanishes.
Ax(E-B) = 0.
This is the desired equation. It is a vector equation in the
imknown vector E. The equation of a plane was seen (page
88) to be a scalar equation such as
E.C =c
in the unknown vector E.
The point of intersection of a line and a plane may be
found at once. The equations are
( - B) =
Ax(E
IH'C = c
AxE = AxB
(AxE)xC = (AxB)xC
A'C E — C-E A = (AxB)xC
A'C E - c A = (AxB)xC
Hence
(AxB)xC +cA ^
A.C
The solution evidently fails when A C • = 0. In this case how-
ever the line is parallel to the plane and there is no solution
or, if it lies in the plane, there are an infinite number of solu-
tions.
r = a;i + yj + «k
p = ui + vj + wk
Hence i-t? = xu + yv + sw = 1.
The quantities u, v, w are the reciprocals of the intercepts of
the plane p upon the axes.
The relation between r and p is symmetrical. It is a relar
r«p = 1
SUMMAKY OF ChAPTEK II
A.B = ^1 ^1 + A^ B^ + ^3 Sg (7)
AxB=-BxA (10)
jxk = — kxj = i
kxi = — ixk = j
(13)
i J k
AxB= A, . (13)'
Bl Ba B.
The dot and the cross in a scalar triple product may be inter-
changed and the order of the letters may be permuted cyclicly
without altering the value of the product ; but a change of
cyclic order changes the sign.
A-^ A^ A^
[ABC] = Bi B^ B^ (18)'
^1 ^2 ^s
= A.C B - A-B C
Ax(BxC) (24)
(AxB)xC = A.CB-C.BA (24)'
(AxB).(CxD) = A'C B.D - A-D B-C (25)
(AxB)x(CxD) = [A C D] B - [B C D] A
= [ABD] C-[ABC]D. (26)
when
- ,
bxc
a'= .-.-, V=
cxa
= f-tT
axb
» ..
^ ,' <* (29)
[a be] [abc] [abc] ^ ^
= b'.b = c'.c = 1
a'.a
r.A = a. (36)
v = Vo + a X r (46)
ExEBCiSES OK Chapter II
(Pv Vv «i)
IS
1 Vx z, 1
^2 1
»8 1
Vi 3. 1
T^
:
CHAPTER III
At
a1
-
-
A
— «
wiU approach
^^ —
dt
which is a vector tangent to the curve at P directed in that
sense in which the variable t increases along the curve.
If r be expressed in terms of i, j, k as
—Ai!
=z ATI i + A!J + A^a k
AiS '
A«
i
Ai!
r = aa + 6b + cc
d"i _d"a d"l d"c
'dT"~'dF^^dT" ""J^"'
THE DIFFERENTIAL CALCULUS OF VECTORS 117
dx
-J- = — a sm i 4- b cos t.
at
Hence —=a cos (« + 90°) -f- b sin (i5 + 90°).
d7^ ~"'
is evidently a differential equation satisfied by the ellipse.
—
dr
dt
= a sinh ^ + b cosh t,
Hence ^=
d^T
r
ordinary calculus.
For
M ~ M Lt A«
Hence in the limit when A < = 0,
_(a.b) = a._ + -.b (3)
dt
^(ax[bxc]) = ax[bx(^^)]+)^;[(^)xc]
s»<'
(7)
120 VECTOR ANALYSIS
Often it is more convenient to use not the derivatives but
r.r = 1.
The locus of the terminus of r is ,a spherical surface of unit
radius described about the origin, r depends upon two vari-
r = f(0,
where some vector function of the scalar t. In most appli-
f is
LiM At dt
C = (9)
As=0 Ai ds
The vector C is coUinear with At and hence perpendicular to
t ; for inasmuch as t is a unit vector At is perpendicular
to t.
P^^^«^ N = txC.
If c be a unit vector collinear with C
n=t X c
n.' di = 0.
Hence c?n • t= 0.
The increment of n is perpendicular to t. But the increment
of n is also perpendicular to n. It is therefore parallel to c.
The tortuosity T is
„
T = —d ^^
d (dx
(txc) = — — X-j-^
IN^ d'^x
,
,.,.
(11)
)
ds ds\ds ds^ VC'C/
„
T = d^T
-T—o X
d s^
d^T
r-s
c? s^
—= + di
1
y"C7C
-T-
ds
X
d^i
ds^ \/c»C
1
dT d^T d
d s ds^ ds a/ Q.Q
The first term of this expression vanishes. T moreover has
been seen to be parallel to C = d^ t/d s^. Consequently the
magnitude of T is T by the unit vec-"!'
the scalar product of
tor c in the direction of C. It is desirable however to have
dx d^T di
— C.T = -C.— 1
X-r-o ^=. - c
ds ds^y^Q.Q
. —
ds
X
d^T d 1
ds^dS'^/Q.Q
But c is parallel to the vector d^ i/d s^. Hence
di d'^r_
ds ds^
Or T= —• (13)
d^i d^
ds^'Ts^
The tortuosity may be obtained by another method which
is somewhat shorter if not quite so straightforward.
r = txc.|-% (14)'
ds
124
de
t X
ds
But
THE DIFFERENTIAL CALCULUS OF VECTORS 125
»! * a2 ~ [^^i ^a]
ai X ag ~ [«! I
ttg]
LiM rAr"] dT
if sbe the length of the arc measured from some fixed point
of the curve. It is found convenient in mechanics to denote
v = s=- (16)
V =v t. (17)
dv d(vi) dv dt
dt dt dt dt
dv d^s
dt'^dt^^^'
dt
= dt
d s
-— _- -— — Cv,
dt ds d t
A= s t + «2 c.
r= a;i + yjH-sk,
v = r = a!i + yj + sk, (15)'
(16)'
(18)'
V x^ + y^ + z^
From these formulae the difference between s, the rate of
change of speed, and A = r, the rate of change of velocity,
is apparent. Just when this difference first became clearly
recognized would be hard to say. But certain it is that
Newton must have had it in mind when he stated his second
law of motion. The rate of change of velocity is proportional
to the impressed force ; but rate of change of speed is not.
r =V=a X r.
'
A = a V — a^ r.
The acceleration A may be given by the formula
f = A = g.
The hodograph of the hodograph reduces to a constant
vector. The curve is merely a point. It is easy to find
the hodograph. Let Vo be the velocity of the projectile
in path at any given instant. At a later instant the velocity
will be
V = Vo + « g.
— r„
r =r
a;
(r — f° )
Multiply by rx. r x r =rx •
consequently
Lm
A« = 02
1 rxCrx^Ar) _ j^^ 1, ^
At = 02"^ ~M ^ -
At~^
^
; v r
r — ro _ Af '
At At
*
A *
*
Af
Hence
Ar
rxr— roXro=-r-XroA« + a;A<rxro + yA«rxn.
But each upon the right-hand side is an
of the three terms
infinitesimal of the secovd order.Hence the rates of descrip-
tion of area at P and P^ differ by an infinitesimal of the
second order with respect to the time. This is true for any
point of the curve. Hence the rates must be exactly equal
at all points. This proves the theorem.
60.] The motion body one point of which is
of a rigid
fixed is at any instant a rotation about an instantaneous axis
passing through the fixed point.
Let i, j, k be three axes fixed in the body but moving in
space. Let the radius vector r be drawn from the fixed point
to any point of the body. Then
= a;i + yjH-2k,
r
dt = xdi+ydi + zdTs.. ,
Moreover i.i=j.j=k»k = l.
Hence i.(ii = j.<Zj = k«cZk = 0.
132 VECTOR ANALYSIS
Substituting these values in the expression for d r.
Let ,. ^j • , . <^^ . , .
'^i ,
Then dr
r
.
= -J-
=aX r.
d t
If r be displaced by a, it becomes
T + d-^T = t + A^xx dt.
If it then be displaced by aj, it becomes
r + tZ r =; r + c^i r + aj X \r + (b.^xx) df] dt.
Hence (^ r = a^ X r d i + ag X r d « + aj X (a^ X r) (d i)^
r =a X r.
a X (a X r).
axr.
This is equal in magnitude to the rate of change of angular
velocity multiplied by the radius of the circle. It does not
depend in any way upon the angular velocity itself but only
upon its rate of change.
61.] The subject of integration of vector equations in which
If then dT = d6,
r =s+ C,
134 VECTOR ANALYSIS
where C is some constant vector. To accomplish the integra-
tion in any particular case may be a matter of some difficulty
just as it is in the case of ordinary integration of scalars.
Example 1 : Integrate the equation of motion of a
projectile.
f = g + b, <
r = f (O-
Since the acceleration is parallel to the radius,
r X if = 0.
But r X if = ^- (r X r).
at
For
j^
(r X f) = r X f + r X if.
Hence (r x f) =
at
and r X r = C,
X C —1 —1
if
^ r"
,
{r.r r
.
- r.r r}.
.-
r • r = r^.
-
J
TT
Hence —X—C =
r
g2
-
r
_
ly
r
-r.
^2
Integrate. Then —— = 6^
- -
r
+ e I,
r .r X C r • r
+ e r -
c
C2 r
But
r-rxC rXr.C C«C
-.2 /.2
136 VECTOR ANALYSIS
I^et
p_ —_ c
and cos v, = cos (r, I).
Or p
r =
1 + e cos tt
P
1 — e^
It is possible to carry the integration further and obtain
the time. So far merely the path has been found.
Aa; = L Aa; J-
This is the partial derivative of Fwith respect to x. Hence
the rates at which V increases in the directions of the three
axes X, Y, Z are respectively
9r 9V 3V
3 as
'
9 y' S z
of V is denoted by V V.
vr = i|r+jf
X
+ k|^.
3 z 9 y 9
(21)
V." How this is to be read is not divnlged. Indeed, for printing no particular
name is necessary, but for lecturing and purposes of instruction something is re-
quired — something too that does not confuse the speaker or hearer even when
often repeated.
THE DIFFERENTIAL CALCULUS OF VECTORS 139
By making use of the formulae (47)' and (47)", Art. 63, page
104, for transformation of axes from and by i, j, k to i', j', k'
actually carrying out the differentiations and finally by
taking into account the identities (49) and (50), may V
actually be transformed into V.
V' = V.
The details of the proof are omitted here, because another
shorter method of demonsti-ation is to be given.
64.] Consider two surfaces (Fig. 30)
lowed when passing from one surface to the other. The rate
of increase will be the quotient of the actual increase d c and
the distance V^r • di between the surfaces at the point
X, y, z in d r. Let n be a unit normal to the
the direction
surfaces and d n the segment of that normal intercepted
between the surfaces, n dn will then be the least value for
d r. The quotient ,
Cb C
^/d I • dT
~
dn
n (23)
VV-dr=^n.dT. (25)
d n
dV
S7V.dx = -~dn = dV. (25)'
_ 5V SV 8V
But dV=yr-dx+^dl,+^dz,
where (d xy + (d yy +(dzy = dr.dT.
lid I takes on successively the values idx, jdy, k rf« the
equation (25)' takes on the values
Vk« idx = 5V
-i^— dx
dx
SV
^^'i^y = 3^dy (26)
5V
VF-'kdz = -p^dz.
dz
SV 3r 9V
dX 3y^ Sz'
VF= (VF. i) + (VF. j) j + (VF.
i k) k.
dV 3V k 5V
Hence by (26) VF=i^- +j + ^. ^ (21)
dT'VV=dV
for all values of d r.
This definition is certainly the most natural and important
from theoretical considerations. But for practical purposes
either of the definitions before given seems to be better.
They are more tangible. The real significance of this last
definition cannot be appreciated until the subject of linear
vector functions has been treated. See Chapter VII.
The computation of the derivative V of a function is most
frequently carried on by means of the ordinary partial
differentiation.
^ . X . V
Vx' + y-^ + a*
1
Hence Vr=
Vx^ + y^ + z^
and Vr = ,
=-
VrTr r
Example 2 : Let „ 1 1
+ ga
r 's/x'^ + y2
V-1 = — 1—
„ . a; .
1
y
*
(a;2 + ya + «2)i
^1 -r -r
and V-=-
r (r«r)' i r*
n— p
Example 3 Vi^ = nr '^
r = »ir" •
(ix+jy).
x^a +I ,,a
and ix + }y = T — k k»r.
VrVr V r«r
d r=dT'a = di'W.
Hence VF=a.
Example 3: Let F= (rxa) • (rxb), where a and b are
constant vectors.
F = r-r a«T) — r«a r«b.
dV = 2dT'i a«b — dr«a r«b — dr»b r«a = di'W.
t = -hVT,
where ft is a constant depending upon the material of the
body. Suppose again that V be the gravitational potential
due to a fixed body. The force acting upon a imit mass at
the point (.», y, z) is in the direction of most rapid increase of
potential and is in magnitude equal to the rate of increase
per unit length in that direction. Let F be the force per unit
mass. Then
F = VF.
As different writers use different conventions as regards the
sign of the gravitational potential, it might be well to state
that the potential Preferred to here has the opposite sign to
the potential energy. If W denoted the potential energy of
a mass m situated at x, y^ z, the force acting upon that mass
would be
T = -VW.
In case V represent the electric or magnetic potential due
to a definite electric charge or to a definite magnetic pole re-
spectively the force F acting upon a unit charge or unit pole
as the case might be is
F = -VF.
The force is in the direction of most rapid decrease of
potential. In dealing with electricity and magnetism poten-
tial and potential energy have the same sign; whereas in
attraction problems they are generally considered to have
opposite signs. The direction of the force in either case is in
10
X
3x 9y 3z
As far as combinations of these are concerned, the formal laws
are precisely what they would be if instead of differentiators
three true scalars
a, i, c
9 9 9 9 ^ ^
9x 9y 9y d
the associative law
9
^ {9 9\ f 9 9 \9 ,..,^''^^'' J.
9 -x\ry9^r\J-x9-y)9^~"'^^'^ =
and the distributive law
9 f 9 9\ 9 9^99 ,_,
^
,^
^
9x\9y 9zJ 9x9y 9x9z ^ >
Sx 3y dz
as a vector. It is not a true vector, for the coefEcients
5 d d
dx 3y 9z
This operator A V is a
• scalar differentiator. When applied
a-vr
without parentheses and no ambiguity can result from the
omission. The two different forms (a . V) Fand a • (V F)
may however be interpreted in an important theorem,
(a V) F is the directional derivative of F in the direction
.
V {x, y, z)
of position in space.
The scalar operator a • V may be applied to a vector func-
tion V to yield another vector function.
/ 5F 5Fi 5FA.
and (a.V)V=(«:^+«,y^ + «3-^jx
,
/ 5F2 SF2 ,
9V^\.
(31)
(a - V) V =a VV .
S
The differentiators —9- ,
d
-r— . ;r-, being scalar operators, pass
dx dy d z
by the dot and the cross. That is
3V 5V SY
VxV=ix^ + jx^ + kX^. (33/
5 3 3 a!
y z
Now
dx 9x 9x 9x
dy''
Then
Hence
Moreover
» J
3V 3V 5V
and VxV = ix^+jx^-fkx^-
doc 9y 9z
(33)'
^ '
V= 1
9x
1-
3y
1- k — 9z
V. V=div V.
The reason for this term is that V-V gives at each point the
rate per unit volume per unit time at which fluid is leaving |
this consider a small cube of matter (Fig. 31). Let the edges
of the cube be dx, dy, and dz respectively. Let
perpendicular to the X
axis. The normal to the
face whose x coordinate is
^JfA
s- j
x4dx , yg.
(fjf t/e
face is
through it is therefore
X
3V
= i.Y (x,y,z) dy dz + i'-—-dx dy dz.
cl OS
The total flux outward from the cube through these two
faces is therefo^ the algebraic sum of these quantities. This
is simply
1 . -;;— dx dy dz = -^ dx dy dz.
dX dX
i .
-:
3y
— dx dy dz and k • -=
OS
— dx dy dz,
The total flux out from the cube is therefore
This is the net quantity of fluid which leaves the cube per
unit time. The quotient of this by the volume dx dy dz of
SV 5V 5V_5F, 5F 5F-3
8 8y 8z 8X dy 8z
Because V V thus
• represents the diminution of density
or the rate at which matter is leaving a point per unit volume
per unit time, it is called the divergence. Maxwell employed
the term convergence to denote the rate at which fluid ap-
proaches a point per unit volume per unit time. This is the
negative of the divergence. In case the fluid is incompressible,
as much matter must leave the cube as enters it. The total
div D = V D = 0.
.
V X V = curl V.
The curl of a vector function V is itself a vector function
of position in space. As the name indicates, it is closely
connected with the angular velocity or spin of the flux at
each point. But the interpretation of the curl is neither so
easily obtained nor so simple as that of the divergence.
Consider as before that V represents the flux of a fluid.
Take at a definite instant an infinitesimal sphere about any
point (x, y, a). At the next instant what has become of the
sphere ? In the first place it may have moved off as a whole
in a certain direction by an amount d r. In other words it
Let a = «! i + 0^2 j + as k
•o is a constant vector.
Vo Hence the term V x Vo vanishes.
V X V = curl V = 2 a, t'--
a = 2Vxv = 2 curl v.
T = v„ + ^ (V X V) X r = Vo + ^ (curl v) x r. (34)
V + v') = Vu + Vv
(u (35)
V.(u + v)=V.u + V.v (36)
Vuv, Vwv, Vm X V.
V« X v = (Vm) X V.
If V is to be applied to more than the one term which follows
it, the terms to which it is applied are enclosed in a paren-
thesis as upon the left-hand side of the above equations.
^ dx
^ dx
The summation extends over x, y, z.
vx(.,) = 2ix(|^v)+2ix(»?,)
9u\ .^A . 9v
r—
.
=201-:) X V + y.ui X
dx
Hence Vx (ttv) = VMXv-|-MVxv.
To demonstrate
V (u • t) =V Vu
. -1- u. VV + V X (V X u) 4- u X (V X v).
U
THE DIFFERENTIAL CALCULUS OF VECTORS 159
V(u.v) = :^iA(u.y)=2ir|^-v
'^ 9x^ -^
+ u.|T'\ -' \9x QxJ
>i
V(u.v)= '^ — .v+T
dx -^ iu.—
3x
Now
IX 5u ,^ Su _- 5u
TX(Vxu) = vx2ix^=2v.^i-2v-i^"
;:^l
ax
= VX(VXU)+>V.i--
^ -^ 9x
5u
0' 2T«T-i = v X (V X u) + vVu.
Hence V (u^v) = v Vn + u^ Vv
+ V X (V X u) + u X (V X v).
The other formulae are demonstrated in a similar manner.
74.] The notation!
V(u.v)„ (44)
Let u= Mj i + ^2 J + '^a ^»
v= + v^} + v^}si.
'yj i
fZ (u . v) = d (u . v)„ + c? (u • v),
or £? (u • v) = u • (Z V + (^ u• V.
= V(uvX + V(uv\
^(y,v) (38)'
V X (u X v)n = u V v--» u V V. • •
u X (V X v)
11
162 VECTOR ANALYSIS
formally as if V, u, v were all real vectors. Then
ux(Vxv) = u«vV — n.Vv.
The second term is capable of interpretation as it stands.
The first term, however, is not. The operator V has nothing
upon which to operate. It thereforemust be transposed so
that it shall have u • v as an operand. But u being outside
of the parenthesis in ux (V x v) is constant for the differen-
tiations. Hence
u•V V = V (u v)n •
V (x, y, z)
/(a3)=/(a;„)+/'(iBo)da3.
That is V (r v)y = v. .
For V (r . v)t = V V r — (V
• X r) x v,
= xi-\-y^ +
x zls.,
Hence V (r • v)^ = v.
V (di- v)y = v,
By(47) v= Vo + V((ir.v)d,+ (Vx v) x<ir
V (rfr. v) = V (cZr.v),,' + V (c?r vV .
Substituting
dT'di = c^-
At each point of this sphere the velocity is
V = Vo + £? r • V V.
In the increment of time B t the points of this sphere will have
moved the distance
(Vo + <Z r • V v) S ^.
The distance between the center and the points that were
upon the sphere of radius di a,t the commencement of the
interval S t has become at the end of that interval B t
dr' = di + dT-VvSt.
To find the locus of the extremity of cZr' it is necessary to
eliminate d r from the equations
dt' = dT + dT'Vv8t,
e^ = di' dr.
dr-dr = c^
has been transformed into an ellipsoid by the motion of the
fluid during the time B t.
T = Vo + |[V(cZr.v)-vJ +|(Vxv)xtZr;
or of the equation (49)
v = Vo+V(rfr.v)^,+ (Vxv)xdr,
v = v„+ [V(dr.v)^,+ i(Vxv)x dr]+|(Vx v) x «?r.
|(Vxv)xdr = | curl v x df r
|v((Zr.v)-v„,
3a! "
9y 9z
dx » dx = d x^ + d y^ + d z^ = c^
goes over into the totality of points upon the ellipsoid of
which the equation is
V-VF, V X VF
THE DIFFERENTIAL CALCULUS OF VECTORS 167
VxVF=curlVF. (52)
the terms cancel out. Later (Art. 83) it will be shown con-
versely that if W possesses no curl,
a vector function i. e. if
SymboUcaUy, -
^ —
V V = 9x^ + 9y^ + ,
9z^
The operator V V is
• therefore the well-known operator of
Laplace. Laplace's Equation
^. VF=^ + ^ + ^ = (53)
V-VF=0. (53)'
dT= pk
-V'{ dt.
dt pk
This is Fourier's equation for the rate of change of tempera-
ture.
Let V be a vector function, and F^, F^, V^ its three com-
ponents. The operator V V of Laplace may be applied to V.
•
Vx(VxV) = V.VV-V-VV.
The terms V •
it operates upon V.
to
Vx Vx V = VV-V-V-VV, (58)
or curl curl V = VdivV — V.VV. (68)'.
V VV = and
ciu:l curl V = V div V.
V div V = V VV. -
To sum up. There are six of the dels of the second order.
V.VF, VxVF,
V.VV, V^-V, V-VxV, VxVxV.
Of these, two vanish identically.
VxVV = 0, V.VxV = 0.
A third may be expressed in terms of two others.
Vx Vx V = VV-V-V-VV. (58)
Then
l^'^^LiM 2
2da;2a=0 a^
l5^u
2 Sa!2
THE DIFFERENTIAL CALCULUS OF VECTORS l7l
2 5s2 a=
Add:
1 /52m 32^ 52«\_ 1 „ „
"1
_LiM r 6
a=0\_ a2 J-
As V and V« are independent of the particular axes chosen,
this expression may be evaluated for a diiferent set of axes,
then for still a different one, etc. By adding together all
these results
1*1 + M2 + • • 6 » terms
'
'^0
lV'Vu=^^ ^
a
Let n become infinite and at the same time let the different
by Mo-
^V.Vm=L^,^^S^°.
6 a^ a==0
vector additions.
d' ,^ da dh
_(a.b)=-.b + a.-,
,
(8)
^(axb)=-xb+ax^, (4)
the unit tangent to the curves directed toward that part of the
curve along which s is supposed to increase.
d s d s^
T=
d^T d^T
ds^' dl^
174 VECTOR ANALYSIS
If r denote the position of a moving particle, t the time,
V the velocity, A the acceleration,
dx . ,^..
d s
(16)
'"^dTt'
dv d^ T .. ,^„
'^ = --Tt=di^ = '- (^«)
9 X " d y Sz
9V 9V 9V
If V is a vector function a • W
of that vector function in the direction a.
is the directional derivative
-9Y SV 9V
V><V=ix|^ + jx|^ + tx|^ (88/
V.V=|5
9
+ |5+«^.,
9 9z
(82)"
cc y
\9x 9yJ
176 VECTOR ANALYSIS
Proof that V V is • the divergence of V and V x V, the cwd
ofV.
V.V = divV,
V X V = curl V.
= Vu + Vv,
V(u + v') (35)
V (u + v) = V u + V V,
. - . (36)
V (« v) = V V + M V T,
. ti • (39)
+ ux(Vxv), (41)
= v,Vxu-u.Vxv,
V.(tixv) (42)
Application to hydrodynamics.
The dels of the second order are six in number.
THE DIFFERENTIAL CALCULUS OF VECTORS 177
VxVr=curlVr=0, (52)
V V X V = div curi V = 0,
•
(56)
sion of a function.
(r-re^')-, (rO+ihe^kx-.
a = (r • r) r, b = (r - a) e, c =r X ^
where d, e, f are constant vectors,
r r^
8. Compute V W, VV . V, and V x V x V when V is
and
(VxV) X W=Vx (Vx W)w + WV.V-V(V.W)w.
CHAPTER rV
THE INTEGEAL CALCTJLTTS OF VECTORS
/W. dT.
If W^W^i + W^i+Wslt,
and di = i dx + J dy + Ts. dz,
That is,
IT
Jv V.dT = V(i)-V(to) = V
o
<ix, y,z)-Y (,x„ y^ ».)•
By definition dr.W=dV
fdV=r(x) - V(To) = V(x, y, z-) - V(x„ y,, z^). (2)
o
F(r) = F(r„).
Given fW'dT = 0.
«/
To show W = VV.
Let To be any fixed point in space and r a variable point.
The line integral
dr
/7
is independent of the path of integration C. For let any two
)
fW'dT+ Jfw-dt = 0,
Jo —C
fW'dT = — JfW',dr.
J —o' o'
rW'di = V {x,y,z).
together.
W • dr= V(x^ y, z).
c?
f = -^k.
The work done by the force when its point of application
moves from the position x^ to the position r is
r=w = g(z„-z').
Or more simply V=— gz.
The force is f = VF= — ^k.
THE INTEGRAL CALCULUS OF VECTORS 183
w=
=/ t I • dv.
And conversely if { = W
the integral vanishes. The force-function and the work done
differ only by a constant.
V= w + const.
In case there is friction no force-function can exist. For the
work done by friction when a particle is moved around in a
closed circuit is never zero.
The by a fixed mass
force of attraction exerted upon M
a unit mass is directed toward the fixed mass and is propor-
tional to the inverse square of the distance between the
masses.
t
M
= -C—,T.
w=— ^
c -^ I'dr.
Jx
^' dr
w
a
184 VECTOR ANALYSIS
By a proper choice of units the constant c may be made
V=- \
-^ + -? + — +••
where r^^^, r-g, are the distances of the attracted unit
• • •
ILW.da. (4)
f da
/xa •
186 VECTOR ANALYSIS
gives the amount of that substance which is passing through
the surface per unit time. It was seen before (Art. 71) thaF
the rate at which matter was leaving a point per unit
volume per unit time was V - f The total amount of mat-
.
V.f = —
dv Jff{.d&.
Js (9)
where the symbol S r has been introduced for the sake of dis-
tinguishingit from d r which is to be used as the element of
Fig. 32.
fW'dT=lfw„'dT+lfv(W'ST^.dt
+ ~ C (V xW) X Si.di.
lfw(W'STy'dT
St X dr
^St X dT = da,
where d a denotes the positive area of the triangular element
of surface. For the infinitesimal triangle therefore the
relation p
JI A
W-dr = 'V X W • da,
holds.
Let the surface S be divided into elementary triangles.
For convenience let the curve which bounds the surface
be made up of the sides of these triangles. Perform the
integration
fW'di
W'dr
?/.
which arise from those sides of the triangles lying within the
surface S cancel out, leaving in the sum only the terms
which arise from those sides which make up the bounding
curve of the surface. Hence the sum reduces to the line in-
tegral of W along the curve which bounds the surface S.
y^ fW'dT = fW'dT.
s'J A Jo
S fW'di= fW'di- Jo
fw-dx,
J Jc
Cd(W'ST-) = 0.
SrW'dT = f\sw-dr-dW'8i\.
3W 5W
But dW = 9W
dX ay
,
-^^ dx+ -^;— dy + -^r-de,
dz
, ,
SW SW^
or dW = 5W ,
-fr—i'dr+ -f;—i'dT + -^r-ls.'dT,
8X ay az
. , ,
and gW = ?^i.Sr+^j.Sr+^k.8r.
dX ay az
X —
192 VECTOR ANALYSIS
Substituting these values
J J ( 3 ax
+ similar terms in y and z. i •
/.
( 1 X
o'W\
—
-T )
•
,^
(5 r X d, r)
,
= SW
-tt— .
,
dr
.
i •
^
3r
9W
— -^— •
^
Sr
.
i • d r.
\ &X J dX dX
or S rw . c? r = fv x W • Sr X <i r.
Srx<ir = £?a,
S ('•W'di = fv xW'da..
/ W'dr
J W'di
THE INTEGRAL CALCULUS OF VECTORS 193
2s/w.dr = 2j'vxW.da.
Or fw-di=f('vxW-d&. (11)
or rr (XT- V x"W)-<Za = 0.
Let the surface S over which the integration is performed be
infinitesimal. The integral reduces to merely a single term
(TJ-Vx W).£?a = 0.
As this equation holds for any element of surface d a, the
first factor vanishes. Hence
IT - V X W = 0.
Hence TJ = V x W.
area d a.
value of d&,
X Exfr.
o
THE INTEGRAL CALCULUS OF VECTORS 196
n.a B>(2a.
Experiment therefore shows that
or f'E-di= f j —i-da.
Jo JJa
V X E = - B, curl E = - B.
It is also a fact of experiment that the work done in carry-
ing a unit positive magnetic pole around a closed circuit is
equal to 47r times the total electric flux through the circuit.
The work done in carrying a unit pole around a circuit is
X "E'dr.
o
VXH=4 TT C.
V V X W = 0,
. div curl W= 0.
VxVF=0, curlVF=0.
According to Gauss's theorem
JJfV'VxWdv= ffvxW'dii.
V'VxWdv= fW'dx = 0,
Jo
V V• X W= 0.
:
JJ a Jo
Apply this to any infinitesimal portion of surface. The curve
bounding this surface is closed. Hence the line integral of
the derivative VF^ vanishes.
V X VV'da = 0.
Vx VF=0.
In a similar manner the converse theorems may be
demonstrated. If the divergence V • TJ of a vector function
TJ is everywhere zero, then TJ is the curl of some vector
function W.
XT =V X W.
TJ = VF.
84.] By making use of the three fundamental relations
between the line, surface, and volume integrals, and the
dels, viz.
I
V (u v) • di = j uV V' diJtrj vVu'dv.
r
The expression [u v]
/uV V'dt — — Jo vV u-
o
I
di. (14)'
Hence
rrr\7-(uv)dv=rrrus7-vdv+ c c ^^ wvdv.
Hence
or
f f i '^U'V dv= j f
uvda— I j I W^'W dv, (17)'
functions of position,
:
=1 I
jV'(vVu)dv— I I IvV'Vudv.
Hence
/ /
/Vm'V«'c?v = / I uVv'da.— I j l uS7''^vdv,
/ (y'V'V>' — 'oS7''^u)dv= I
j {u'^v — VVu)'d2L.
/ /
following
/ / I
wV>i'V'>'dv= I j uwVv'da— I j j
u'^'[wVv']dv,
F=tan"^?^'
^'^— 1.
-2 I
,,2^+ ^2 ,
„2J-
y = r sin 0,
fv-dv=f^-^=lfd0.
Form the line integral from the point ( + 1,0) to the point
lying below that axis. The value of the integral along the
first path is
1 r'
- I de=7r;
- Ir~"
1
along the second path, dd^ir.
From this it appears that the integral does not depend merely
upon the limits of integration, but upon the path chosen.
202 VECTOR ANALYSIS
the negative of the value
the value along one path being
along the other. The integral around the circle
which is a
F= tan 1 = ^*
(—1,0) had not included the origin the values of the integral
would not have differed. In other words the value of the
(^-^
''"^^'^
/ now a diaphragm be inserted reaching from
the surface of the cylinder or column to the
surface of the cube the region thus formed
bounded by the surface of the cylinder, the
cny., surface of the cube, and the two sides of the
diaphragm is acyclic. Owing to the inser-
Fig. 34.
tion of the diaphragm it is no longer possible
to draw a circuit which shall pass completely around the cyl-
the given cyclic region and the two faces of each diaphragm.
In acyclic regions or regions rendered acyclic by the fore-
hold true. For cyclic regions they may or may not hold
true. To enter further into these questions at this point is
+
,
^2
.
J
,
+ «2 ^»
Pj = ajj i -f- yj j + sj k
is the vector drawn from the point (ajj, y^, Zj) to the point
(asj, yj, z^"). As this vector occurs a large number of times
in the sections immediately following, it will be denoted by
This is what is termed the potential at the point (aj, y^, «i)
J>otV=fff^^^^l;^^^dx,dy,dz,. (22)
Potr=/ / /^d«2.
JJJ »-j2
= r ainB cos
x <f),
dv — r^sva.6 dr dO d<^.
Let the point (ajj, y^, z{) which is jQxed for the integration
be chosen at the origin. Then
or simply PotF= I
/ /
f^-sin^ dr dO d^.
r =R
Hence the triple integral taken over all space outside of a
sphere of radius B (where B is supposed to be a large quan-
tity) is less than 4 ir K /B, and consequently converges as far
as regions distant from the origin are concerned.
210 VECTOR ANALYSIS
If the function V remain finite or if it become infinite so
Vr<K
r =R r =R
C f fvrsmedr d0 d^ < f f Cxdr d0 d<f>.
= rzzO
r =R
dO d^ = 'iirKB.
f f flTdr
r =
Hence the triple integral taken over all space inside a sphere
of radius iJ (where B is now supposed to be a small quantity)
is less than ^iir K B and consequently converges as far as
regions near to the origin which is the point (aj^, y-^, Zj) are
concerned.
If at any point (Xg, j^, Zj) not coincident with the origin,
(Xg, y2> Zg) % ^^^ square of the distance of that point from
(Xg, y2» Z2) I'^fi^O'i''''^ finite as that distance approaches zero, then
the integral converges as far as regions near to the point (Xg, j^, Zg)
[PotF]^,,„...
is therefore
Hence ^m \
[Pot n». + A«..y..».- [Pot ^]«.,y.,^ ?^
Aa;i = 0^ AiBj ^
LiM ^
[PotF(a!2 + Aa;2,y2,g2)]^,y„^-[Potr(a;2,ya>g2)]».,v.,^. ?
Aa3i = 0^ Aa;, j'
12
= rff yjo^i+^^yyvH)
^^ J
err V{.x^+^x^,yi,^)
dv.
t As all the following potentials are for the point x^, y^, xy the bracket and
indices have been dropped.
THE INTEGRAL CALCULUS OF VECTORS 213
LiM
Aa;2=0( Aa!2 =
Lm fff n^yy^^H) ^^
when A a?! approaches zero as its limit the regions mand m',
which are at no point thicker than A 05, approach zero ; M'
and M
both approach J!f as a limit.
t There are cases in which this reversal of the order in which the two limits
are taken gives incorrect results. This is a question of double limits and leads to
the mazes of modem mathematical rigor.
X If the derivative of Fis to exist at the snrface hounding T
the values of the
function V must diminish continnonslj to zero npon the surface. If Fchanged
suddenly from a finite value within the surface to a zero value outside the de-
rivative 3VI Q xi
would not exist and the triple integral would be meaningless.
For the same reason V
is supposed to be finite and continuous at every point
within the region T.
214 VECTOR ANALYSIS
Then if it be assumed that the region Tis finite and that V
vanishes upon the surface bounding T
LiM
VPotF=PotVF: (27)
= i Pot —
9V
H- j Pot
SV k SV
;^ + Pot 11-.
THE INTEGRAL CALCULUS OF VECTORS 215
Vi X Pot W = Pot Va X W,
curl Pot W = Pot curl W
(28)
or
of higher osder
V V Pot V= Pot V
• . VF, (30)
VV Pot W = Pot VV W,
^ . (32)
sPotr err i 3V
9 —-fff
u J J
a! J
-V-^H
'
ji ^12 ""a
+ A«>,^O^J dv,
J J„,
Lm ^ r r r
-^-.^ox^JJJ^
—-—
F(a;a. y^ h) dv.
J
+ ^ ^2>
Hence ^fff ^^^^ ^a^^a)
^
J J ri2
Consequently
3 PotF r r r sv "^
, r r v
iBi J J J uTy^dx^ * J Ja ri2
THE INTEGRAL CALCULUS OF VECTORS 217
i«c?a<r2 dd d<f>.
1 9V ,
J J JM ri»
12 5*2
5Potr , 9V
=„
Pot •
9xi 9x,2
gPotF
f''ffLk/^'"'^-fJ''-^J-'"-
3X
1 sv
dv^
ISL
may or may not approach a limit when i? becomes smaller
and smaller. Hence the equation
5PotF
— =„ 5V
^ Pot -z
-1 £r
///,
THE INTEGRAL CALCULUS OF VECTORS 219
5PotF _ Pot5F
IS.
need not become zero and consequently the equation
5PotF „ SV
=Pot
need not hold for any point (a;^, ^j, z^ of the region. But
if V becomes infinite at x^, y^, z^ in such a manner that
ri<K,
then the surface integral will approach zero as its limit and
the equation will hold.
Finally suppose the function V remains finite upon the
surface S bounding the region T, but does not vanish there.
In this case there exists a surface of discontinuities of V.
Within this surface V is finite ; without, it is zero. The
surface integral
V
i'da
fl S 7-12
220 VECTOR ANALYSTS
does not vanish in general. Hence the equation
PPotr
— ji
^ ^3V
= Pot TT
cannot hold.
Similar reasoning may be applied to each of the three
partial derivatives with respect to ajj, y^, z^. By combining
the results it is seen that in general
ViPotF=PotV2F.
da.
a •
12
iir^j- '•12
(37)
5PotF r r rCx2-xi)V
Or
9x^ 7=ffP^^^^' <''^
sPotw r r riXi-xi)VT
and
s,
—
Sail oy^
— + —
5PotF
T
^— Szj =
Hence
In like manner
V Pot F= /// ^ dv^. (39)
V X Pot W =
//J '-1^ dv^, (40)
JJJhlzJL^?ll^l^dx^dy^dz^=^MaxW. (44)
VPotF=NewF
V X Pot W = Lap W (45)
V Pot W = Max W.
-
The first is written Ifeiv V and read " The Newtonian of F."
224 VECTOR ANALYSIS
The reason for calling this integral the Newtonian is that if
in magnitude and has the direction of the vector r^j from the
point (x^, y^, z^)to the point (x^, y^, Sg). Hence the force is
fi2 ^"^a
^ 12
"-^"'///^^^^"-'-^///^ dv.
(^2-^1)^
v„
O U %J '
V!,
° ilx "'y\U U U /^i2
dC^
^ 12
15
226 VECTOR ANALYSIS
Integrating over all space, the total magnetic force acting at
(43)'
<^^-^^^-^^-^"--"^^^
j.LapW = d.,
///
k . Lap W ^JJJ^^l^^^J^^^^lAIdv,
In terms of the potential (if one exists) this may be writt
.
1 .
-r
Lap
^
W = SVotZ
— ?j
SPotF
-z
J
•
. Lap «r
T
W = 5PotX
— ,r
SPot^
^
,
,,„,,,
(43)"
T „ = 3Poty
k.LapW
,
-^:j
SPotX
=
dv„ = Max
/// ^l' I.
Let
I (^2, ^2' ^2) = i^ ('^a. ^2' ^2) + J -^ («'2' Vv ^^ + k C(a;2, y^i ^2)
ri2 • I = (x^ -x{)A + (2/2 - 2/1) 5 + («2 - z^) G.
=
Max I
/// \A{^-x) + B{'n-y)+ C(^-z) \~^dv.
According to the notation employed for the Laplacian
Max W .^fff(-^—^)^+(y2-y^)^+(^.-^l)^,,^,
(44)'
(29)
V-VPotr= V.Pot VF=Pot V-VF.
But by (45) V Pot F = New F,
and V.PotVF=Max VF
Hence V- V PotF= V NewF= Max VF -
or V Pot V X W = 0..
or VxPotVr=0.
Hence V x New V = Lap VV=0. (50)
And by (58), Chap. IH. VxVxW = VV.W-V-VW,
V.VW = VV-W-VxVxW.
Hence V V W = New V W - Lap V x W,
• •
(51)
or V V W = V Max W — V x Lap W.
•
Poisson's Equation
V-VPotr=-47rr, (52)
^^""'IJIvrJ'-
Id X2
Vi Vi Pot F= Vi
. . New F= Max V^V= f f C ^^^'J^^ d v^.
But Vi — = - Va —
and Va - (f V^ — ") = V^ — • V2F+ F Vg • Vg —
\ ^12/ ^12 *'l2
:
Integrate
But V2 . V„ —= 0.
= ff VV, — 'da.
The surface integral is taken over the surface which bounds
the region of integration of the volume integral. This is
w 1 „ 1 r
r 12
|//rv.2....
|<JJ^,,,,=2^.
Hence when B becomes infinite the surface integral over the
large sphere approaches zero as its limit.
1 r
Vi da. = - — r^sme dd d6.
7-3 ^
7-12
-//'Tsin^ de d^.
THE INTEGRAL CALCULUS OF VECTORS 233
(^i» Vv '^i) which has been selected as origin. Then for the
surface integral V is practically constant and equal to its
value
W = Xi + Fj + ^k,
V- VPotX=-4 7rX,
V. VPotF=-4 7r F,
V V Pot ^ = - 4 TT ^.
.
Consequently
V V Pot (Xi +
• Fj + ^k) = - 4 TT (Xi + Fj + -^k).
234 VECTOR ANALYSIS
Theorem : If V and W are such functions of position in space
that their potentials exist in general, then for all points at which
V and W are finite and continuous those potentials satisfy
Poisson's Equation,
V V Pot F= - 4 ttF,
. (52)
V V Pot W = - 4 W.
• TT (52)'
or MaxVF=-4 7rF.
In a similar manner it was seen (51) Art. 91 that
= New V W — Lap V x W.
Hence V Max W — V x Lap W = — 4 W, tt (54)
W = :j—
47r
LapVx
^
W — 47r
:i— NewV'W. (55)
^ '
Let W = Wi + W 2'
V . Wj = J— V
4 TT
. Lap VxW = -i—
4
V'Vx ^^^P ^i-
TT
Hence V Wj = div Wj =
•
(58)
As before
and ^V
47r
X New V.W = 7-VxVPotV.W
47r
= 0.
The divergence of the first part and the curl of the second
part of W are therefore zero.
236 VECTOR ANALYSIS
;i^VxVxPotVxW
-i-VxLapVxW = 47r *^
Air
J- VV PotVxW = -^ VPotV-Vx W =
= 0,
for V V XW= . 0.
Hence ^- V VPot V X W = V X W.
•
Air
-;— Lap V X W
into which W is divided. Hence as the second part has no
curl, the third part can have none. Moreover
4 TT
"^ New V • W.
— LapW2 = 0,VxW2 = 0.
(61)
— New
~ and — 'w • or — div
are inverse operators. That is
- ^ New V
4 TT
•
Wa = - V -^
4
. New Wj = Wj.
'JT
(62)
•
and also — ~.
47r
— Max and V
are inverse operators. That is
_V.ri-NewF=F (63)
-J— Pot
47r
and V xV X or curl curl
are inverse operators. That is
7— Pot and — VV
47r
-. — Pot and — V V •
_PotV.VF=-V.VT^PotF=r
47r Air
—V V . and V x V X
are equivalent
- V VWi = V. X V X Wj. (67)
V V . and VV •
7rF=- V.NewF"
4
and 4 TT W = V X Lap W — V Max W
—
Air
-.— New Max
is an operator which is equivalent to Pot. For solenoidal vector
Pot ( U, F)
=fff^(^v Vv h) Pot ^^^1
(73)
One integration may be performed.
(74)
New ( ^'^=///J//w(a3i,2/i,%) •
^ Vix^,y^,z^)dv^ dv^.
MaxiW,V)=ffJJJJvix^,y„z{)^^.W(x^,y^,z^)dv,dv^.
''
(77)
16
:
=-
f f fuFotVdv
[NewtrPotFJdw
4k7r
f f fV'
+ f f f New U New F . cZ v.
If now W' and W" exist oftly in finite space these surface
Given VF=0.
To show F= const.
Choose a fixed point {x-^, y^, Zj) in the region. By (2) page
180
But JvV-dr=fo.dT = Q.
Hence V{x, y, z) =V {x-^, y^, z^) = const.
Theorem : If V (a;, y, «) be a scalar function of position
in space which possesses in general a definite derivative V V;
if the divergence of that derivative exists and is zero through-
out any region of space,^ finite or infinite but necessarily
continuous; and if furthermore the derivative V F vanishes
at every point of any finite volume or of any finite portion of
surface in that region or bounding it, then the derivative
vanishes throughout all that region and the function F re-
duces to a constant by the preceding theorem.
1 The term throughout any region of space must be regarded as indnding the
boundaries of the legion as well as the region itself.
THE INTEGRAL CALCULUS OF VECTORS 245
Given V V F"=
• for a region T,
To show V— const.
Since VF vanishes for the portion of surface S, Fis certainly
constant in S. Suppose that, upon one S and in the
side of
region T, V were not constant. The derivative V V upon
this side of Shas in the main the direction of the normal to
the surface S. Consider a sphere which lies for the most
part upon the outer side of S but which projects a little
VF.<?a
in the surface integral all have the same sign and cannot
cancel each other out. The surface integral of V F over
that portion of S which is intercepted by the spherical sur-
face vanishes because V F is zero. Consequently the surface
integral of VF taken over the entire surface of the spherical
segment which projects through S is not zero.
Hence C C^V-d& = 0.
the divergence V -
[ f7 V F] of the product of U and VV
exists and is zero throughout and upon the boundaries of T
and at infinity ; and if furthermore V be constant and equal
to c upon all the boundaries of T and at infinity ; then the
function V is constant throughout the entire region T and
is equal to c.
SuMMAUT OP Chapter IV
The line integral of a vector function W along a curve C is
defined as
CvV'dv=0 (3)
or
and if
Jjv • da
=J W ' dr, then TI =V X W. (12)
/ / I uV 'Vdv= j I uv da,—
'
j j
jVu'vdv, (17)
Green's Theorem:
= CI'vVu'd&-jjfvV-Vudv, (19)
Kelvin's generalization:
Pot
V=fff^^^^^:^^ dx, dy, dz,. (22)
VPotF=PotVF; (27)
V Pot W = Pot V W,
. . (29)
and V V Pot r= - 4
. tt W. (52)'
F=7^V.Newr,
4 IT
(53)
— Lap^ V X W — 4Tr
W = 47r 1 — New V ;; • W. (55)
^ ^
;j—
4 TT
Lap V X Wi =V X ;j—
4 TT
Lap Wi = W^. (60)
^ Lap W2 = 0, V X Wj = (61)
--r— New
47r
V . Wa = - V . — New Wj = Wg.
4 TT
(62)
254 VECTOR ANALYSIS
f-V. ^NewF=
47r
r
(63)
--^
47r
Max VF= F.
— Pot VV
47r
- . Wj
Wa^ =-V V—
J-•
.
477
Pot Wa = Wa- (65)
_ Jl
4 TT
Pot V VF^-V.V:;—
.
4 TT
Pot F=F
(^^)
1 1
-
__Pot
47r
V . VW = - V. V-p-PotW=W.
47r
V V Wjj = VV Wg
- . (68)
may be. Formulae (71) to (80) inclusive deal with these inte-
grals. The chapter closes with the enunciation of a number
of theorems of a function-theoretic nature. By means of
tion. The work done in the text has for the most part assumed
that the potential exists. But many of the formuloR connecting
Newtonians, Laplacians, and Maxwellians hold when the poten-
tial does not exist. These are taken up in Exercise 6 referred to.
THE INTEGRAL CALCULUS OF VECTORS 255
ExEBCisEs ON Chapter IV
^=L^ X dx
is a vector quantity. It may be called the shew line integral
of the function W. If c is any constant vector, show that if
1 The first four exercises are taken from Foppl's Einfuhrnng in die Max-
well'sche Theorie der Electricitat where they are worked out.
256 VECTOR ANALYSIS
and H.c = c. I
/ / V. Wc?a- J J^V (W-da)
I
H =J fg[ V. W cZ a - V (W • c? a)j.
S=fJ^Vda.
This is a vector quantity. Show that the surface integral
of F taken over any closed surface is equal to the volume
integral of VF taken throughout the volume bounded by
that surface. That is
//.^'^^ = -///^^^'
Hence conclude that the surface integral over a closed sur-
face vanishes if F be constant and conversely if the surface
1 ^ r,,
^n ^ 12
Then New V
=fff^i Pn ^(^2' Vv H) <^^p (81)
Vli'l3=-V2Pi2 (84)
JJJ'^iPn Vdv,=fffv,(p,,V)dv,
-JJJpn'^.Vdv,.
17
258 VECTOR ANALYSIS
By exercise (3)
/ / /
"^2 O12 ^ '^^a =J J -^12
'^^^•
V X Lap W = Lap V X W,
V New F=Max VF,
=
V Max W = New V W. •
Then prove
or if W be a vector function
M{W)^JJfr,^Wdv^. (90)
V ^ (W) = ^ (V W) = Max
. . (7-2 W), (92)
fi-(F)=-^PotPotF (96)
and vector potential, that is, to V, V«, Vx, and Pot with the
auxiliaries, the Newtonian, the Laplacian, and the Maxwellian.
The treatment of the second topic will likewise introduce
novelty both in concept and in notation — the linear vector
function, the dyad, and the dyadic with their appropriate
symbolization.
The simplest example of a linear vector function is the
product of a scalar constant and a vector. The vector r'
r' = cr (1)
r are
i • r, j • r, k • r.
axes, a point
r = ia;-l-jy + k2
becomes r' = i c^ + j a; Cg y + k Cg «,
or r' = i Cj i . r +j (jj j • r + k Cg k • r.
If T — xa + + zc, where [a b c]
1/Ta ^ 0,
\ m) \m
and /('±')^l/(r).
\m, J m
Hence f (\ n—^=f(
m \'m J
— A= —f
J m
(j).
That is, equation (5) has been proved in case the constant a
is a rational positive number.
To show the relation for negative numbers note that
/(0)=/(0 + 0) = 2/(0).
Hence /(0) = 0.
But /(O) =/(r-r) =/(r +(-r)) =/(r) +/(-r).
Hence / (r) = — / (— •
LiM
X f./C")=/(jf„(.r))
/(a!r) = x/(r).
Hence Lim = a / (r)
05 z^ w
y ^^ ,)
264 VECTOR ANALYSIS
^™ (irr) = ar.
x =a
Hence / (« r) = ^ / (')
known.
Let l=/(a),
m=/(b),
r=/(o).
= a;a + yb + zc.
r
r' = (i cj + 1 j C2 j + k cj k) . r.
(6).
Definition : An expression a b formed by the juxtaposition
of two vectors without the intervention of a dot or a cross is
called a dyad. The symbolic sum of two dyads is called a
dyadic binomial ; of three, a dyadic trinomial ; of any num-
ber, a dyadic polynomial. For the sake of brevity dyadic
binomials, trinomials, and polynomials will be called simply
dyadics. The first vector in a dyad is called the antecedent
and the second vector, the consequent. The antecedents of a
dyadic are the vectors which are the antecedents of the
individual dyads of which the dyadic is composed. In like
manner the consequents of a dyadic are the consequents of
the individual dyads. Thus in the dyadic (7) aj, aj, aj are • • •
T>=0'r. (8)
By definition <? • r = aj b^ • r + ag bg • r 4- ag bg • r + • • •
T'0=0c'r. (9)
r = xa, + i/'b-\-ee,
utive. That is
(a + b)c = ac + bc /-i -1 \
and a (b + c) = a b Xa c.
[(a + b) c] . r = (a + b) c T = a c . r + b c • r = (a c + b c) • r
and
+ cl + cm + CI1+ • •
+
The dyad therefore appears as a product of the two vectors of
which it is composed, inasmuch as it obeys the characteris-
tic law of products — the distributive law. This is a justifi-
+ agikiH-ttgakj + a33kk.
.X=W 'X
for any value of and the dyadics by (10) must be equal.
r
for all values of s and r. Let s and r each take on the values
i,j,k. Then (14)
= a^ al + aj2 am + ^13 an
The expressions (15), (16), (16)' for are unique. Tw'o equal
dyadics which have th,e same three non-coplanar ante-
cedents, a, b, c, have the same consequents A, B, C — these
however need not be. non-coplanar. And two equal dyadics
which have the same three non-coplanar consequents 1, m, n,
associative.
The most general product conceivable ought to have the
property thatwhen the product is known the two factors are
alsoknown. Certainly no product could be more general.
Inasmuch as scalar multiplication is to be associative, that is
Let a = aj +
i ^2 J + «8 ^>
b = Ji + Ja +
i i Jg ^
LINEAR VECTOR FUNCTIONS 273
a'= a/ + a^ + a^ k,
i j
+ <*i' + aa'S/
ji + a2*3 Jlf
jj
But aj 61 = a/ 6j'.
This shows that the product of the lengths (including sign)
are equal and the theorem is proved.
The proof may be carried out geometrically as follows.
Since ab is equal to a'b'
ab • r = a'b' • r
which shows that the products of the lengths are the same.
18
274 VECTOR ANALYSIS
The indeterminate product a b imposes Jive conditions upon
I
the vectors a andb. The directions of a and b are fixed and
likewise the product of their lengths. The scalar product
a • b, being a scalar quantity, imposes only one condition upon
a and b. The vector product a X b, being a vector quantity,
imposes three conditions. The normal to the plane of a and
b is fixed and also the area of the parallelogram of which they
are the side. The nine indeterminate products (12) of i, j, k
into themselves are independent. The nine scalar products
are not independent. Only two of them are different.
i • i=j = l, •
j =k • k
and i.j=:j •i = •k = k'j=k-i = i.k = 0.
j
<Px = aj X bi + aj X bg + ag x bg -1
(19)
Or 0s=i-0-i+i'0'3+'k'0'k, (20)'
<2>^ = (j k - k (? j) i + (k <^ 1 - i
.
(?•
. . . • • • <? • k) j
1 A Bubscript dot might be used for the scalar of * if it were snfiScientlj distinct
and free from liability to misinterpretation.
276 VECTOR ANALYSIS
Products of Dyadics
That is, the antecedent of the first and the consequent of the
second dyad are taken for the antecedent and consequent
respectively of the product and the whole is multiplied by
the scalar product of the consequent of the first and the
antecedent of the second.
Thus the two vectors which stand together in the product
(ab).(cd)
are multiplied as they stand. The other two are left to form
a new dyad. The direct product of two dyadics may be
defined as the formal expansion (according to the distributive
law) of the product into a sum of products of dyads. Thus
'P
<P = (ai b 1 ++ aj ba + Eg bg . .
.)
+ bj • Cj aj dj <Jf
bj • C2 aj dj + bg • Cg aj dg + • •
+ bg • Ci ag dj + bg . C2 ag da + bg • Cg aj dg + • • •
+ (23)"
Let S^0'¥.
To show J2.r= (p. (f .r),
or (<P . ?F) . r = (? . ( r . r). (24)
(a b • c d) •= b.c (a d
r • r) = (b c) (d r) a,
• •
ab . (^c d• r) = a b c (d • • r) = (b c) (d r) a. • •
Hence (a b c d) r =
• • ab • (c d • r).
(0.¥).S = 0.(¥.S).
a :
(ab • r) • c d = b r (a cd) = (b
• • • r) (a • c) d,
ab • (r • c d) = ab d (r c) = b
• • • d (r • c)
Hence (a b • r) • c d ^ ab • (r • e d).
(r r = r X ((? f) = r X
X (/>) . W, • >
s . r X <? and (? x r • s ^
s . = s (r X 0),
r X <2> .
^^^^^
0xi'S = (0xt)'a,
LINEAR VECTOR FUNCTIONS 281
0-(TXa) ?t (<2>.r) X s,
(ab) X (o d) =a (b X c) d.
<P = al + bm + cn.
If 1, m, n are coplanar one of the three may be expressed
in terms of the other two as
l = xra. + ya.
(Z>=:dp + eq
and that the consequents 1, m, n of <? were non-coplanar.
This supposition leads to a contradiction. For let 1', m', n'
The vectors 1', m', n' exist and are non-coplanar because
1, m, & have been assumed to be non-coplanar. Any vector r
may be expressed in terms of them as
T = xV + ym' + gn,'
<P = dp.
Let ?r=<^xp = dp xp = dO = 0,
S'"=al xp-f-bmxp-f-cnxp.
284 VECTOR ANALYSIS
From the second equation it is evident that ¥ used as a
postfactor for any vector
= 0. 1-0
Hence aslxp+ymxp + ^nxp
must be zero for every value of r, that is, for every value of se,
y, z. Hence
with the antecedent of and t any value collinear with the con-
286 VECTOR ANALYSIS
sequent of ^ ; hut no other values. The dyadic used as a
prefactor reduces any vector r to the line of the antecedent
of In particular any vectors perpendicular to the con-
0.
sequent of <P are reduced to zero. The dyadic used as a
postfactor reduces any vector r to the line of the consequent
and Cj
s' = I « (bj . Ci) +y (bi . C2)} &i + {x (ba . Ci) + y (ba • Cj)} a^.
LINEAR VECTOR FUNCTIONS 287
(bi X bj) .
(ci X C2) ^ 0.
The vector b^ x h^ is perpendicular to the plane of the con-
sequents of ; and c^ X Cj, to the plane of the antecedents of
¥. Their scalar product vanishes when and only when the
vectors are perpendicular — that is, when the planes are per-
case.
I = ii + jj + kk. (33)
If I.i = i,
0.1 = 0.
In like manner it may be shown that 1.0 = 0.
Theorem : If a', b', c' and a, b, c be two reciprocal systems
of vectors the expressions
= &\ + bm + on
is ^n idemfactor 1, m, n must be the reciprocal system of
a,^, c.
10
290 VECTOR ANALYSIS
In the first place since is the idemfactor, it is a complete
dyadic. Hence the antecedents a, b, c are non-coplanar and
possess a set of reciprocals a', b', c'. Let
By hypothesis = r. r . (/>
Let 'W = l.
To show W .0 = 1.
T.(0- ¥)=T'1 = T,
T'(0 ' T) .0 =1 .0,
T. (0 .¥).0=(t .0^.(W .0) =! .0.
¥.0 = 1.
If the product of two dyadics is an idemfactor, that product
may be taken in either order.
109.] Definition: When two dyadics are so related that
their product is equal to the idemfactor, they are said to be
1This necessitates both the dyadics * and T to be complete. For the product
of two incomplete dyadics is incomplete and hence could not be equal to the
idemfactor.
LINEAR VECTOR FUNCTIONS 291
0=¥,
0.0-^ = I,
and W'W--^ = l.
To show 0-^ = W-\
= al + bm + en,
(?
For (al + bm + cn) • (ra' + n' b' + n' c') = a a' + bb' + cc'.
That is,
if 0.W=O'S, then?r=-0,
and if (P r = (?
. • s, then r = s, (37)
Hence t x r =t x s.
0. W' jT-i .
<p-i = (P
.
( ?r . ?r-i) .0-^ = 0. 0-^=zi.
Hence (<P . f) . ( ?F-i .
(«>-i) = I.
t' = (P'T
?p- = ii-jj-kk.
Geometrically the transformation
r'= ¥'T
T .0 = 0^.1 (9)
{^.¥)a=¥a'0o, (40)
{0 . W)a- r = T '(_(!>-¥) = {t . 0) . W,
I .0 = 0C.T,
(t.0^ .W^¥o'(i-0)=¥c'0o-i.
Hence (0 -
¥)c=^o' ^c-
(^0-i-y^=(^0^)-i^0-\
(42)
dyadics.
T .0 = -T, = 00-
For anti-self-conjugate dyadics
T .0 = -0 -T, = -0o.
Theorem : Any dyadic may be divided in one and only one
way into two parts of which one is self-conjugate and the
other anti-self-conjugate.
But (0 + = 0^ + 0^c=<l>c+
0c-)^ <^,
and = 0n,
^(^0-0^-)
0=z0' + 0".
unique.
Suppose = al + bm-|-cn,
<?
But a 1 r — 1 a r = — (a X
• X • 1) r,
the vector of 0.
This theorem follows as a corollary from equations (44).
Theorem : Any dyadic may be broken up into two parts
of which one is self-conjugate and the other equivalent to
minus one half the vector of used in cross multiplication.
cx(r-)-s) = cxr-fcxs.
Hence it must be possible to represent the operator c x as a
dyadic. This dyadic wiU be uniplanar with plane of its
antecedents and consequents perpendicular to c, so that it
will reduce all vectors parallel to c to zero. The dyadic may
be found as follows
=I • (c X I) • r= (c X I) r. •
Hence c x r = (I x c) •r = (e x I) • r,
and r X c = r . (I X c) = r (c X I),
.
(46)
This may be stated in words.
LINEAR VECTOR FUNCTIONS 299
(I X c)^ = (c X I)^ = I X c = c X I.
X^=;:^^rX^ (48)
X* = I,
X« = X.
The dyadic X therefore appears as a fourth root of the
idemfactor. The quadrantal versor X is analogous to the
imaginary V— 1 of a scalar algebra. The dyadic X is com-
plete and consists of two parts of which I x c is anti-self-
conjugate ; and c c, self-conjugate.
I = il + jj + kk
into i, j, and k successively.
These expressions represent
quadrantal versors about the axis 1, j, k respectively combined
with annihilators along those axes. They are equivalent,
when used in direct multiplication, to i X, jx, kX respectively,
(a X b) X I =I X (a X b) = b a — a b. (50)
For
r'= (2>.r
By expressing V in nonion form, the equation r' Y r' = 1 is seen to be of the second
• •
degree. Hence i' describes a qnadric surface. The only closed qoadric surface
is the ellipsoid.
LINEAR VECTOR FUNCTIONS 303
<P = ai + bj + ck.
= (ai +bj +
r' ck).r,
= (bj + ck)T,
r'
of finally a sphere.
Theorem: Any self-conjugate dyadic may be expressed in
the form ^ ..,,..,
= an +bjj + ,,
cki. ,cc\
(55)
Since ^= ^0,
(<?2-a2i).i/ =o
(P2_a21=(j2_a2)jj + (c2_a2)]£t^
(<?2_a2I).i=0.
k and k' are parallel. The dyadic therefore takes the form
*
Double Multiplication
•' '
ab:cd = a«cb*d. (56)
ab:cd = cd:ab,
1 The researches of Professor Gibbs upon Double Multiplication are here
printed for the first time.
LINEAR VECTOR FUNCTIONS 307
and the distributive law both with regard to the dyads and
with regard to the vectors in the dyads. The double dot
product of two dyadics is obtained by multiplying the prod-
uct out formally according to the distributive law into the
sum of a number of double dot products of dyads.
+ e3d8 + ...)
ab H cd = cd^ ab,
308 VECTOR ANALYSIS
and the distributive law both with regard to the dyads and
with regard to the vectors of which the dyads are composed.
The double cross product of two dyadics is therefore defined
as the formal expansion of the product according to the
distributive law into a sum of double cross products of
dyads.
+ (67)"
(i<Pl¥)lii^0-Q¥^^Q-). -
ij:ij=i.i j.j =l
ij:ki = i.k j«i = 0.
0^ ¥:Q
is a scalar quantity. The multiplication with the double
cross must be performed first. This product is entirely in-
dependent of the order in which the factors are arranged or
the position of the dot and crosses. Let a b, c d, and ef be
three dyads,
ab^cd:ef = [ace] [bdf]. (59)
0^0:0 = ^J>2 :
0''0:0
0s = -^ = [ab c] [Imn]. (62)
(64)
Let = al + bm4-en
(P
^^^>
i0.W\=0,W,.
Choose any three non-coplanar vectors 1, m, n as consequents
of and let 1', m',n' be the antecedents of W.
(P= al + b m + en,
r = l'd + m'e + n'f,
?r=ad + be + cf,
= [a b c] [1 m n],
<?g
Ig — X.
Theorem: The product of the second and conjugate of
a dyadic is equal to the product of the third and the
idemfactor.
0^.0,= 0^1, (68)
0o = l& + mb + nc,
0^ 0g= [Imri] (h X c a + c x a b + a x b c).
.
= [a be] I.
(? = aii ii + + ^isik
ai2 ij (13)
+ a2iJi + a22JJ + Hzi'^
+ agi k i + agj k j + agg k k.
The conjugate of has the same scalar coefficients as 0, but
they are arranged symmetrically with respect to the main
diagonal. Thus
0o=a^^ii + + agj ik,
aji ij
"21 *28
This minor is taken with the negative sign. That is, the
coefficient of 1j in (Pg is what is termed the cofactor of the
coefficient of ij in the determinant. The cofactor is merely
the first minor taken with the positive or negative sign
according as the sum of the subscripts of the term whose
first is under consideration is even or odd.
minor The co-
any dyad in 0^ is easily seen to be the cofactor of
efficient of
^11=
LINEAR VECTOR FUNCTIONS 317
*u *ia *18
and is written
0^=\0\ (72)'
0..0,= 0A (68)
or 0„-^ = —- 0„
or 0-^ = —- 0„„
""
0s
+ ^12Ji+^22JJ+ A2Jl^
.+ ^igki + ^agkj + ^askk.
Hence 0-^ =
«11
318 VECTOR ANALYSIS
If the determinant be denoted by D
0.W= (ttji 611 + ai2 ^21 + «18 ^31) " + («ii hi + "'la ^22
"''
^31 "31 + *32 "32 + *33 "SS"
"11 *12 "IS "11 "12 "18 «11 ^11 + «:12 ''21 "f" *19 6,
18 "81
Hi *22 «2S "21 "22 "28 0521 &11 + a.'22 ''21 "" "^'O °81
28 "i
a 11 +
^12 ^22 + a,„6.
18 "82
''^12 *11 "is "'" '^12 ''28 + ais 6.
IS "38'
*
«21 ^12 + S2 *22 + 28 ''32 ttgi 613 + a.'22
23
'' '''23 "88' (76)
*81 ^12 + 32 "22 + ^83 "32
"'
''Sl "is + *S2 ^28 + «;33 °83'
If (? = aH-bm + en,
Hence ^11 12
•^10 -" 18 '11 "12 *18
<P. "^8'
the scalar of (P, the scalar of the second of (P, and the third
or determinant of (?. If be expressed in nonion form these
quantities are
320 VECTOR ANALYSIS
^s
LINEAR VECTOR FUNCTIONS 321
Summary of Chapter V
A vector r' is said to be a linear function of a vector r
when the components of r' are linear homogeneous functions
of the components of r. Or a function of r is said to be a
linear vector function of r when the function of the sum of
two vectors is the sum of the functions of those vectors.
<P • r = aj bj • r + ag bg • r + aj bj • r + • • • (8)
21
322 VECTOR ANALYSIS
Two dyadics are equal when they are equal as operators
upon all vectors or upon three non-coplanar vectors. That
is, when
,x = W • X f or aU values or for three non-
coplanar values of r, (10)
coplanar values of r,
(a + b + c+.--)(l-l-m-|-nH ) = al + am + aiiH
-I- b 1 -1- b m + bn -I- • • •
-I- cl + cm + en -I- • •
+
(11)'
uaique.
The symbolic product ab known as a dyad is the most
general product of two vectors in which multiplication by a
scalar is associative. It is called the indeterminate product.
The product imposes five conditions upon the vectors a and
b. Their directions and the product of their lengths are
determined by the product. The
and vector products
scalar
= ai bi + aa bj + ag bj +
<?>^ . • . • • .
(18)
^x = ai X bi + aj X bj + ag X bg + •
. •
(19)
<p^ = (j k-k
. <?> . j) 1+ (k
. (^ -1 - . <P - 1 . (?> • k) j
+ (i.0.j-j.0.i)]!. (21)
I = ii + ji + kk. (33)
Or I = aa' + bb' + cc'. (34)
I
r. (<?-(?,)=- -J
rx<P,. (44)
Also c X r = (I X e) • r = (c X I) • r, (46)
c X (2> = (I X c) . = (c X I)
(? . <^>.
(a X b) X r = (b a — a b) r •
r X (a X b) = r (b a - ab). . (51)
ab = a c b .d,
: cd •
(56)
ab|^cd = axc bxd. (57)
/ = al + bm + cn,
<P
0^=\0^0'.0=\a,\ic]l\mn\. (62)
i^c)s = K (63)
(<?H1)^=((Z»,)-1,
0^.0^=0^1 (68)
<^>3 9i 0, is complete
0^ = 0, 0^ = 0, 0^0, is linear.
satisfied by a dyadic 0.
0s -<Ps 0^+ 0^+03 1 = 0. (79)
0^s
cients 0g, 0^g, and (Pg are the three fundamental scalar in-
variants of 0.
Exercises on Chapter V
1. Show that the two definitions given in Art. 98 for
a linear vector function are equivalent.
abxc+bcxa+caxb = [abc]I
and b X c a + c-x a b + a x b c = [a b c] I.
x' ^0-x.
This equation therefore may be regarded as defining a trans-
formation of the points P of space situated at the terminus of
r into the point P', situated at the terminus of r'. The origin
r =b -1- a; a
to unity.
Let = al + bm-F en
<?>
s' = <^>o • s.
334 VECTOR ANALYSIS
It is important to notice that the vector s denoting a plane
[d'e'f']-^ = [def].
The volume [a b cj to [d e f] is as 0g is to
ratio of the unity.
But the vectors d, e, f were any three vectors which deter-
mine a parallelopiped. Hence all volumes are changed by
the action of in the same ratio and this ratio is as 0^ is to 1.
Let = xi + yj + zk
I
are changed into the vectors i', j', k' and any vector r is
changed from its position relative to i, j, k into the same posi-
tion relative to i',j',k'. Hence by the transformation no
change of shape is effected. The strain reduces to a rotation
which carries i, j, k into i', j', k'. Conversely suppose the
body suffers no change of shape — that is, suppose it subjected
to a rotation. The vectors i, j, k must be carried into another
right-handed rectangular system of unit vectors. Let these
be i', j', k'. The dyadic may therefore be reduced to the
form
<P = i'i + j'j+k'k.
Definition : A dyadic which is reducible to the form
<?>o
= ii' + Jj' + kk',
(?>.<?>^ = i'i' + j'j' + k'k' = I
0-^=0^.
Hence the first part of the theorem is proved. To prove the
second part let
<P= ai + b j + ck,
<p^ = ia-|-j b + kc,
If 0--^
=0c, 0'0c = I-
Hence aa + bb + cc = I.
336 VECTOR ANALYSIS
Hence (Art. 108) the antecedents a, b, c and the consequents
a, b, c must be reciprocal systems. Hence (page 87) they
must be either a right-handed or a left-handed rectangular
system of unit vectors. The left-handed system may be
changed to a right-handed one by prefixing the negatiTC
sign to each vector. Then
<?>-i'i + yj + k'k,
<?> .
(2>c
= I, 0^=\0\=1. (2)
(P is a versor.
4>. 00 = 1, I (? I = - 1. (3)
(P is a perversor.
i'= i'
22
338 VECTOR ANALYSIS
<P = ii + cos g- (jj + kk) + sin q (kj - jk). (4)
+kk = I-ii,
jj
kj — k = I X j i.
Hence *
<^ = i i + cos 2 (I — i) + sin g i I x i.
(5)
Hence the vector r has been rotated in its plane through the
angle q. If r were any vector in space its component parallel
to a suffers no change ; but its component perpendicular to a
is rotated about a through an angle of q degrees. The whole
vector is therefore rotated about a through that angle.
Let a be given in terms of i, j, k as
a = ttj i + aj j + as k,
+ + 0jj-aiJk,
a3Ji
— a2ki + «! kj + Okk.
Hence
<^ = {^i^ (1 — cos q) + cos 2} i i
(?5 = i • i + cos 2 (J
.
j + kk) + sin 2 (k o
j —j . k),
sm a •
ta,nq =- ^— = -\/<?x
\, ^ ^
"
<^x
•
(8)
-|^ 1 + cos 2 1+ 0a
The tangent of one-half the angle of version is therefore
determined when the values of 0^ and 0g are known. The
vector <?x and the scalar 0^, which are invariants of 0, deter-
mine completely the versor 0. Let ft be a vector drawn
in the direction of the axis of rotation. Let the magnitude
of Q, be equal to the tangent of one-half the angle q of
version,
a = a.a-tan^l^.
i^'
The vector Q, determines the versor completely. Q will be
called the vector semi-tangent of version.
By (6) a versor was expressed in terms of a unit vector
parallel to the axis of rotation.
c — ft X c into c + ft X c.
^
(C —a X C) • (C + ft X c) _ C-C — ftXC'ftXC
c2 (1 + tan2 I q) c2 (1 + tan2 1 2)
c2 (1 - tan2 I g)
= cos q
c2 (1 + tan2 1 q-)
(c — a X c) X (c + a X c) _ 2 c X (a X c)
~ + tan^ 1 q)
c2 (1 + tan2 Iq-)
2
ca (1
2 c^ tan
2f-== sin q. V v/
<?> = (I + 1 X a) .
(I -1 X a)"^ (10)'
(i-t-a.a)c.
TTa:a=('-'^*>
carries the vector c — a x c into the vector c, if c be perpen-
dicular to a as has been supposed. Consequently the dyaditj
(I + Ixa)''
i + a-a
^^gg + d + ixa)'^ ^
,
i + a-a ^
Hence substituting
ft = ai + 6j + ck.
i2 = (2bb-I).(2aa-I). (14)
^ = 4a.b ba - 2 aa — 2bb +
— 2aa I,
fx = 4 c b c X b, .
?F^ = 4(c.b)2-l
(f (p)^ = 4c-a c X a, .
...-—-^
_ (b X c) X (a X b) _ [a b c] b
^"~
^ a«bb'C a»bb»c'
T
fti >: ftg + fta X fti
"
a -
a • b b • c
b b b
= (a X b) X
• (b c) a . . c a • c b •
fti • fta
a.bb«c a«bb»c a>bb*c
— r— — =
9i * G
Hence ;
1 — ft»^ • ft,.
a.b b. c ^
ii. Hence
(?» = i i +' (cos 2'
<?i + sin g- 0^y
(P"= ii + cos" 2 <Pi" +™ cos"~^ q sin g #^"-1 • (?>2 H
n(n-l){n-'2)
sm.nq = ncos"~'^qsva.q
jn— — cos"~''2'sm''2'H- • •
Thus the ordinary expansions for cos wg' and sinwg are
obtained in a manner very similar to the manner in which
they are generally obtained.
The expression for a versor may be generalized as follows.
Let a, b, c be any three non-coplanar vectors ; and a', b', c', the
reciprocal system. Consider the dyadic
That is, the ellipse and the circle are cut from the same
cylinder. The two semi-diameters i and j of the circle pro-
ject into the conjugate semi-diameters a and b of the ellipse.
The radius vector r in the ellipse projects into the radius vector
f in the unit circle. The radius vector r' in the ellipse which
is equal to (?• r, projects into a radius vector r' in the circle
such that
f = cos Qj -f- 2') i -f- sin (p -f- q) j.
—=
27r
S
m,
(P'n =I
may then be regarded as the mth root of the idemfactor.
In like manner if q and 2 ir are commensurable it is possible
to raise to such a power that it becomes equal to the idem-
factor and even if q and 2 tt are incommensurable a power of
may be found which differs by as little as one pleases from
the idemfactor. Hence any cyclic dyadic may be regarded as
a root of the idemfactor.
1 It is evident that fixing the result of the application of * to all radii vectors
in an ellipse practically fixes it for all vectors in the plane of b and c. For any
vector in that plane may be regarded as a scalar multiple of a radius vector of
the ellipse.
RBTAWieifS-^ff:B STRAINS 351
case the strain has the lines i', j', k' for principal axes: in
the second, i, j, k. In both cases the ratios of elongation are
the same, — a to 1, & to 1, c to 1. If the negative sign occurs
before the product the version and pure strain must have
associated with them a reversal of directions of all vectors in
space — that is, a perversion. Hence
Theorem: Any dyadic is reducible to the product of a
versor and a right tensor taken in either order and a positive
or negative sign. Hence the most general transformation
representable by a dyadic consists of the product of a rota-
tion or version about a definite axis through a definite angle
accompanied by a pure strain either with or without perver-
sion. The rotation and strain may be performed in either
order. In the two cases the rotation and the ratios of elonga-
tion of the strain are the same ; but the principal axes of the
strain differ according as it is performed before or after the
ROTATIONS AND STRAINS 353
and a', b', c' are two reciprocal systems of vectors. Neces-
sarily a, b, c and a', b', c' are each three non-coplanar.
Definition : A dyadic that may be reduced to the form
h =p cos q
and c=^sin2'. (24)
Then 'V
0. W= ¥. = aia2aa' + i?ii>2Cos
(P + (bb' + cc')
(g-j g-j)
</>'a=aa (27)
or (<?> - a I) . a = 0.
The dyadic <P — aI is therefore planar since it reduces vectors
in the direction a to zero. In special cases, which are set
But 0:1= 0s
<P2 : I = <^>25.
•2i, = aa.
is a solution of a cubic equation. Let x replace a. The
cubic equation becomes
For let x = a, x = b, x =c
be the three roots of the equation. The dyadics
Then <? a = a a,
•
Hence m (a — c) a + m (& — c) b = 0,
and m (a — c) = 0, w (b — c) = 0.
Hence m = or a = c, = or 5 = ti c.
a' •(<?>- a I) =
and let the lengths of a and a' be so adjusted that a' • a = l.
This cannot be accomplished in the special case in which a
360 VECTOR ANALYSIS
and a' are mutually perpendicular. Let b be any vector in
the plane perpendicular to a'.
a' • ((2> - a I) . b = 0.
Hence (#— al)»b is perpendicular to a'. Hence (?«b is
\_0.a. (?2 . b . = [a
b] <?>3 0-la b]. (31)
But = aa,.
• a.
^2 = a-i 0y
ROTATIONS AND STRAINS 361
= cos
n q. (36)
bj = cos 2 b + sin q c.
0-aI
and the plane of the consequents are perpendicular the
and a' used in the reduction to cyclotonic form are
vectors a
perpendicular and it is impossible to determine a' such that
a • a' shall be unity. The reduction falls through.
If « =± 1, b_i + bi =± 2 b.
Let = 2b.b_i + bi
Choose c = bi — b = b — b_j.
W'e=pc=p})^— p'b= • 0.
order or arrangement
1 In these cases it will be seen that the cubic equation has three real roots.
In one case two of them are equal and in the other case three of them. Thus
these dyadics may be regarded as limiting cases lying between the cyclotonic in
which two of the roots are imaginary and the tonic in which all the roots are real
and distinct. The limit may be regarded as taking place either by the pure
imaginary part of the two imaginary roots of the cyclotonic becoming zero or by
two of the roots of the touic approaching each other.
ROTATIONS AND STRAINS 863
I + cV.
(I + cb') • a3a = iea,
(I + cb') • icb = ajb + xa,
(I + c b') • asc = c. a;
?'=(?-aI = ^ji-f-CkX+^kj.
This may be reduced as follows to the form
ab' + bc'
'^-^H^^Y^^
Then (j^ =ADk, y=21f^
b =^ j e' = i.
ROTATIONS AND STRAINS 365
dyadics.
366 VECTOR ANALYSIS
135.] A more systematic treatment of the various kinds
of dyadics which may arise may be given by means of the
Hamilton-Cayley equation
0^-0^0-^+0^^0-0^1 = (39)
(<?-aI).((^>-5I).((Z>-cI) = 0. (40)
If, however, the cubic has only one root the Hamilton-Cayley
equation takes the form
(i0-al).(^0^-2pcosq0+p^T)=O. (41)
I. (<P-aI).((2>-6I).(<?-cI)= 0.
11. {0-al)'{0^-'2,pQosq0 + pn) = (i.
Summary of Chapter VI
The transformation due to a dyadic is a linear homogeneous
strain. The dyadic itself gives the transformation of the
points in space. The second of the dyadic gives the trans-
formation of plane areas. The third of the dyadic gives the
ratio in which volumes are changed.
that
<?= -(i'i-i-j'j + k'k) (ly
or that (P <?o = I>
. <P3 = - 1
^= a.a = tan='^^
rT^ (9)
(10)
a-a ' "V a-a; ' VoTa
or (2> = (I-f-I xa)-(I-IxQ)-i (10)'
^ aa
^= + (I +1 X a)2
(!«)"
l + ftft
^^(i-a.ft)i + 2aft+2ixft
^
i + a-a ^
l-tti-tta
ttg fti
((P-aI).((P-6I).((?-cI) = tonic
CHAPTER VII
MISCELLANEOUS APPLICATIONS
Quadric Surfaces
r • • r = const. (1)
If x = xi + yi-\- »k,
r . <? . r = ±— ± —± — (4:')
^ ''
a^ h^ c^
r • <? • r = const.
This is evident from the equations of the central quadric
surfaces when reduced to the normal form. They are
x^ y^ z^
± -r ± — ± — = const.
a^ 0^ c^
i i 11 kk
The corresponding dyadic <P is (?=±—5±-r-± —r-
0^ a^ c^
r • r, (r • a) (b • )r', r • c, d • c,
But r • r =r I • • r,
and (r • a) (b r) = r
• • ab • r.
374 VECTOR ANALYSIS
Hence the most general quadratic expression may be reduced
to
T'0'T + r'A + C = O,
where <? is a constant dyadic, A a constant vector, and
a constant scalar. The dyadic may be regarded as self-
conjugate if desired.
by replacing r by r' — t.
(r' - t) - t) + (r' - t) A + C =
. (? • (r' .
equal. Hence
iA = (?).t
r • (P . r = 1.
dr ' (P 'T + i • dT — 0.
Since is self-conjugate these two terms are equal and
dx-0-T = O. (5)
H= 'T (6)
-T
n = .
i (P-r = N.
Hence r . (2> . r =W .
<2>-i . .
0-i . N = N 0'^
• • N.
N.<^>-i.W = l. (8)
a^ 0^ c^
Let r =s+ a; a.
r . (P . r = (s + a; a) • (P • (s + a; a) = 1.
Hence s + (?>.s + 2a;s.<P.a + a;2a.<?>-a = l.
b • (P • a = 0.
The symmetiy of this equation shows that a is a radius
vector in the plane conjugate with b. Let c be a third radius
vector in the ellipsoid and let it be chosen as the line of
intersection of the diametral planes conjugate respectively
with a and b. Then
a . (? . b = 0,
b . <P . c = 0, (11)
c . (2> . a = 0.
378 VECTOR ANALYSIS
The vectors a, b, c are changed into <P • a, <P • b, ^ • c by
the dyadic 0. Let
)-i — aa + bb + cc.
(pj-i = [abc]2.
:
Hence [a b c] = ahc.
0f^ = i .
<p-i . i -i- j .
(^>-i
. j -I- k .
(?*-i
. k.
Hence Os =
i.(?.i
:^ H
\
,
—s—
' '\
r H
k.<?>.k
;;
But the three terms in this expression are the squares of the
reciprocals of the radii vectors drawn respectively in the i, j,
k directions. Hence
380 VECTOR ANALYSIS
The sum of the squares of the reciprocals of three mutually
perpendicular radii vectors in an ellipsoid is constant. And
in a similar manner : the sum of the squares of the perpen-
diculars dropped from the origin upon three mutually perpen-
dicular tangent planes is constant.
139.] The equation of the polar plane of the point deter-
mined by the vector a is ^
s • <Z> . a = 1. (13)
(x + yy (x + yy (x + yy
If the terminus of s lies in the polar plane of a the two values
of the ratio x:y determined by this equation must be equal
in magnitude and opposite in sign. Hence the term vaxy
vanishes.
Hence s • (P • a =1
is the desired equation of the polar plane of the terminus
of a.
s • ' z^ — 1,
1
or 8 . ^ - a = - •
Q = T -0 'T-1 =
and P=r.c-C=0,
the equation (r • <P • r — 1) + A (r • c — C)^ =
represents a quadric surface which passes through the curve
of intersection of Q and F and is tangent to Q along that
curve. In like manner if two quadrics Q and Q' are given,
Q =i . (P-r-1 =
Q'=T' <P' 'T-1 = 0,
the equation (r • -1 — 1) + k(i • 0' •x-l) =
represents a quadric surface which passes through the curves
of intersection of Q and Q' and which cuts Q and Q' at no
other points. In case this equation is factorable into two
equations which are linear in and which consequently rep- r,
s . (2> . r = 1.
If this plane contains A, its equation is satisfied by a. Hence
the conditions which must be satisfied by r if its tangent
plane passes through A are
a . (P = 1,
• r
r . 'T = l.
Let c + a =p and c — a = q.
Then <P = 51 + ^(p q + q p). (14)
pq + q p.
r . <P . r = 1.
Substituting the value oi 0, x • = 1 becomes
> t
BI ' T + T • -p q«r = l.
Let r -
p = 71
B I ' T + n q.r — 1 = 0.
This is a sphere because the terms of the second order all
have the same coefficient B. If the equation of this sphere
be subtracted from that of the given quadric, the resulting
equation is that of a quadric which passes through the inter-
section of the sphere and the given quadric. The difference
q • r (r • p — ») = 0.
Hence the sphere and the quadric intersect in two plane
curves lying in the planes
q • r = and r • p = ».
:
5 r • r =l
and is of constant length. The section is therefore a circular
section. The radius of the section is equal in length to the
mean semi-axis of the quadric.
For convenience let the quadric be an ellipsoid. The con-
stants A, B, are then positive. The reciprocal dyadic 0~^
may be reduced in a similar manner.
^ , i i i i
Is k
A £ G
11.
B \B J \A bJ
r . (P . r =1 or N • (p-i . N= 1.
For -N N+N . . uV . W = 1.
N N = ^.
.
N-u =
is the equation of a cylinder of which the elements and tan-
gent planes are parallel to u. If then N N • is constant the
cylinder is a circular cylinder enveloping the ellipsoid. The
radius of the cylinder is equal in length to the mean semi-axis
of the ellipsoid.
Let r . <P . r =1
and r = 0.
. a
Differentiate : rf r <P
• r = 0, -
£?r a = 0.•
[a r (?> . r] = 0. (16)
Conversely if [a r <P • r] = 0,
dr may be chosen perpendicular to their common plane.
a r • r = 0.
QUADRIC SURFACES 387
r . ?r2 . r = 1, (17)
instead of r • <P • r = 1.
. ii jj kk
and ii_U_!ii.
a c
r . W ' r . r = 1,
or '(?r.r). (?P'-r) = l.
¥.e,
a' . a' = b' . b' = c' . c' = 1,
a' . V - b' . c' = c' . a' = 0.
Hence the three radii vectors a', b', c' of the unit sphere into
which three conjugate radii vectors in the eUipsoid are trans-
r . <P = 1,
. r
and f • ¥ 'i = 1.
By means of the strain <?* the radii vectors r of the first
ellipsoid are changed into the radii vectors r' of a unit sphere
T'=0i-T, t''t'=1.
By means ¥-i the radii vectors r' of this unit
of the strain
sphere are transformed in like manner into the radii vectors f
of the second ellipsoid. Hence by the product r is changed
into r.
r = ¥-i ^» . r. (19)
QUADRIC SURFACES 389
faces is
a^ — +ji
r^+
6^ —n=
— w 7^-7 c'^
1- (20)
= -^^— + "5
a^ — «!
1
&2 — «! c^— «!
i i i i k k
= (n^ - n{) I.
I . .T =l
and r . ?? • r = 1
and r • ^T r = 1. •
s- (p-i. = 1, 8
s . s' = 0.
But r = f-i • s' = (?-i • s = (f-i + a; I) . s
= ¥~^ . 8 + a; s,
In like manner
r _ (2)-i
. 8 = f-i . s' - = xl)'s'=0-^'s' -X
(<P-^ s'.
a; 8 . s' =s .
<p-i • - 8 <P-i 8 = 8 <p-i - 1.
s' • • . . s'
Hence s = 0, • s'
. i i ii kk
w= a^ —''dn + b^ —^idn + ^^
— dn c^
,
r . (/> . r = 1,
and f • (<P + Idn) f = 1, •
By (19) x = i<^
+ ldn)-^' 0i'T,
f = (I -h dn^-i • r,
_ dn dn x
2 2 a*
dn
= 1 — ——-
X •
The ratio of these components is -
Ja^ -dn a - 1
H —
dn
i «'* «
a 2 a^ X
a
T
In
vi
like manner -!^
by
\^l^—dn —
d^H
Pot
d ^ + ^D + VF=0, V.D = (1)
Operate by V X V X.
V X V X Pot —
d^B
-f-47rVxVx(2>.D = 0. (3)
The last term disappears owing to the fact that the curl of
PotVxVx^+47rVxVx(^>.D = 0. (3)'
But V x V X = VV - V . . V.
. „ dB
Remembering that V • D and consequently V -jj
• and
,
d^ D — 4 tt
V • vanish and that Pot V V • is equal to the
dt^
equation reduces at once to
'^^^
dt^
V V. <? • D- VV . (? . D, V . D = 0. (4)
^^"^
Moreover -r—s = — n^H.
a V
Hence if the harmonic vibration D is to satisfy the equa-
tions (4) of the medium
and m A = 0.•
(6)
THE PROPAGATION OF LIGHT IN CRYSTALS 395
m • m. ^ _
. D. (5)'
396 VECTOR ANALYSIS
Hence the wave slowness s due to a displacement in the
direction ais equal in magnitude (but not in direction) to the
But the first term contains <P - a twice and vanishes. Hence
(I -s = 0, . 8 C?> + s s . <^>) 3
Hence
((p-l_8.sI)3 + S.(<2>-l-8.8l)o"^.8((?^l-S.sI)3 = 0.
Dividing out the common factor and remembering that is
self-conjugate.
1 + 8 . (<P-1 - 8 • 8 I)-l . 8 = 0.
or 1 + 8 •
—B
-, '8 = 0,
1 • a
8.1.8
+s- -8=0 rt
8 . S
r
I —8 .
7n
8 (P
/ 8 • 8 <P \ .
g . ( I J 1.8 = 0.
V I -8 • 8 (^/
Hence s . := ^ -8 = 0. (14)
I —8 . 8 <P
^ ii jj kk
Let ^ = 7^ + P + ^
898 VECTOR ANALYSIS
^ ^ii+/'-^'lJJ+/^-.Vfe-
<k2 0.2 «2
a2 52 c2
X y X z
xy h'^ — s^ + y^ y z --
0. (13)'
X z y z c^ — s^ + z^
s^ — a^ s^ — h^ s^ — c^
a^ x^ b^ y^ c^ z^
_
s2 _ ^[2 s2 _ 52 s2 _ g2
*^ y^ ^^
or
"
mZ"*' l^Z^TITf!
62 a2 c2
The equation is
s
+
-t
:5+
a
:2
o
= (15)
1
wave's advance. And the velocity with which the ray travels
B • da — a, • ds S'(P«a = 0,
or d a • (s — a. s • - a,) = 0.
But since v' da = 0, v' and s — a s • <? • a have the same
direction.
v' = a; (s —as . <P • a),
s.<P«a — a'<P«as'(?'a
v' . ' a,= = 0.
s • 8
Hence the ray velocity v' is perpendicular to • a,, that is, the
ray velocity lies in the tangent plane to the eEipsoid at the
extremity of the radius vector a drawn in the direction of the
displacement. Equation (17) shows that v' is coplanar with
a and s. The vectorsa, s, <P • a, and v' therefore lie in one
The angle from s to v' is equal to the angle from a to > a,.
=
s»a s»s = a-a a»a' = l,
are called the primary optic axes. These are the axes of equal
wave velocities but unequal ray velocities.
In like manner v' being coplanar with a and • a
have been constants. For the elements of the theory and for
elementary applications these constant dyadics suffice. The
introduction of variable dyadics, however, leads to a simplifica-
tion and unification of thedifferential and integral calculus of
— xi + yj + zk,
I
dt = dx i + dy + dzk, i
.3W
XT
Hence «W = aJ r (.5W
^^ • J i
( dX
,
h ] -z
ay
,
f-
,
k
5W)
—— \
az
SW 5W + 3W
VW = i^
ax
+ ay k-;f^.
3z
j ^ (1)
3x dx Sx
By 9y * dy
.9X .SF 3Z ,
az az a z
if W = Xi + Fi + Zk.
z
dW = dVV^dT.VVV,
where VVr= ii ^^
9
+
x''^
ij- ^-^
3x5y
+ ik
dxdz
dy dx Sy^ 3y d
v.vr=(vvF).= ^ + .^ + 5^.
406 VECTOR ANALYSIS
If an attempt were made to apply the operator V symboli-
cally to a scalar function V three times, the result would be a
sum of twenty-seven terms Uke
iiiTp-^,
o'a!'*
ijk^i
c)
—X —y dz
a
T
^,etc.
„
V X
^.90.90
= ix^^ + jx
90
9.x dy 9z
9
V. (^ = i.-_ + j
9X
If
<j] ^^= u i + V j +w k,
Or if (P = i u +j V + k w,
VARIABLE DYADICS 407
„, ^ 5 9 3
,
(a.V)(? = a,^ + a,^ + a3^. (9)
^ 9^ 9^ 9^
(V . V) W = V . (VW) = V V W. .
V (v w) = V V w + V V w.
• •
[V (v X w)]w = V V X w,
[V (V X w)]v = [- V (w X v)], = - V w X V.
Hence V(vxw) = Vvxw — v^xv.
148.] It was seen (Art. 79) that if denote a curve of
which the initial point is r„ and the final point is r the line in-
tegral of the derivative of a scalar function taken along the
curve is equal to the difference between the values of that
function at r and r„.
In like manner
J
I dr - VW = W Qc') — Vf (r^),
and fdi-VW = 0.
Jo
It may be well to note that the integrals
j j
da.V X 0= jdT'0.
The line integrals are taken over the complete bounding curve
of the surface over which the surface integrals are taken. In
like manner the following relations exist between volume and
surface integrals.
C C Cdv V X W= Cfd&xW
f f f dv V' 0=ffd&'
f f fdvV X 0=ffdax 0.
410 VECTOR ANALYSIS
The surface integrals are taken over the complete bounding
surface of the region throughout which the volume integrals
are taken.
Numerous formulae of integration by parts like those upon
page 250 might be added. The reader will find no difficulty in
obtaining them The integrating operators may
for himself.
also be extended to other cases. To the potentials of scalar
and vector functions the potential, Pot (?, of a dyadic may be
added. The Newtonian of a vector function and the Lapla-
cian and Maxwellian of dyadics may be defined.
NewW= fff'-M^^-^q^l^L^dv,,
'JJ J ^12 ;i
= u + vj + w
i k.
a;
=/i (^, v), y=f^ (m, V), z=zf^ (m, v),
or r =f (m, v').
dF=dT' VF=0.
Hence the derivative Vi^is collinear with the normal to the
surface. Moreover, inasmuch as F and the negative of i^when
1 Much of what follows is practically free from the use of dyadics. This is
F (x, y, z) = const.
the normal V F lies upon that side upon which the constant
increases. Let V i^ be represented by N the magnitude of
which may be denoted by N, and let n be a unit normal drawn
in the direction of N. Then
(s-r). Vi?'=0.
and in hke manner the equation of the normal line is
(s-r) X Vi^'^O,
or s = r + A;V.F
where ^ is a variable parameter. These equations may be
translated into Cartesian form and give the familiar results.
150.] The variation (Zn of the unit normal to a suiface
plays an important part in the theory of curvature. c?n is
perpendicular to n because n is a unit vector.
d-a. =- —V jy2
Fa- -dVF,
iV
THE CURVATURE OF SURFACES 413
iV iV'^
dn = d(^VF)=zVFd^
\J!i N + ^dVF.
N
J
^^(I-nn).VV^.(I-nn).
Let
N
Then dn = (fr • (P. (4)
radius vector r.
where i' and j' are two perpendicular unit vectors lying in the
tangent plane and a and I are positive or negative scalars.
The vectors i', j' and the scalars a, b vary from point to point
of the surface. The dyadic <P is variable.
151.] The conic r • (? • r = 1 is called the indicatrix of the
surface at the point in question. If this conic is an ellipse,
that is, if a and b have the same sign, the surface is convex at
the point ; but if the conic is an hyperbola, that is, if a and b
have opposite signs the surface is concavo-convex. The curve
r . (p . r = 1 may be regarded as approximately equal to the
intersection of the surface with a plane drawn parallel to the
tangent plane and near to it. If r • <P • r be set equal to zero
the result is a pair of straight lines. These are the asymp-
totes of the conic. If they are real the conic is an hyperbola
if imaginary, an ellipse. Two directions on the surface which
are parallel to conjugate diameters of the conic are called con-
jugate directions. The directions on the surface which coin-
cide with the directions of the principal axes i', j' of the
indicatrix are known as the principal directions. They are a
special case of conjugate directions. The directions upon the
surface which coincide with the directions of the asymptotes
of the indicatrix are known as asymptotic directions. In case
the surface is convex, the indicatrix is an ellipse and the
asymptotic directions are imaginary.
In special cases the dyadic <?may be such that the coeffi-
cients a and 6 are equal. may then be reduQcd to the
form
<P = a(i'i'+j'j') (5)'
THE CURVATURE OF SURFACES 415
dnxdr = 0. (6)
[n dn dT^=:0. (7)
= dv dn •
dn dv = dr
' ' dr.
dn . (Z r = 0. (8)
153.} Let P
be a given point upon a surface and n the
normal to the surface at P. Pass a plane p through n. This
plane p is normal to the surface and cuts out a plan© section.
Consider the curvature of this plane section at the point P.
Let a' B© normal to the plane section in the plane of the
gectian. n! coincides with n at the point P. But unless the
plane p cuts the surface everywhere orthogonally, the normal
n' to the plan© section and the normal n to the surface wiU not
coincide, d n and d n' wfll also be difEereiit, The cuorvature
of ihe plane eeotiott lying in p is (Art. §7).
£t^^
~ ds"ds^'
As far as numerical value is concerned the increment of the
Unit iangent t and th© increment of the unit normal n' are
^dn' dt dn' • dt
~ ds d s ds"^
«? (n • d r) = d n . dr + n • (?2 r = 0,
d(v! • dt^ =dn' dx + . n' > d^i=0.
Hence dn dx +
- n ' d^i = dn' > dx + ji! * d^x.
Since n
418 VECTOR ANALYSIS
at • di dr ' ax
iV
VVjP— 2iiii.VVi^+nn.VVi^.nn
JSf
(yVF-)a = V .VF,
0^ = 1^^- VF-WF-VF
V'VF - \- " • (13)'
as follows:
(I - nn)2 ' (VV-?')2 • <J - °°)2
(I-nn)2 = nn.
Hence
V^.(VViOo-V^ = VJ^VJ^;(VVi02 ....
^^B = ^, '
T' ^ ^
tx dn = 0. (15)
Hence m. • d n =
d (m • n) = =m ' dn + n • d m.
Hence n • dm = 0.
Moreover
THE CURVATURE OF SURFACES 421
n • t X d t = 0,
n X t . dt= (IT)
or m • d t = 0.
[n d r d^x] = 0. (18)
=a i'i' + 6 j'j'
integral.
m- dt. (20)
L'
For any closed curve this integral may be evaluated in a
manner analogous to that employed (page 190) in the proof
of Stokes's theorem. Consider two curves G and C" near
THE CURVATURE OF SURFACES 423
together. The variation which the integral undergoes when
the curve of integration is changed from C to C" is
S / TO. » dt.
8Jm.dt=Js(m.dt)=rSm.dt+ fm-Sdt
d (m • 8 1) =dm ' St +m • d St
S
I
m-dt= I
Sm- dt— f dm- St
Sni«dt=:8m-I«c?t = Sm»nn-dt,
for t • dt and S m m vanish. A similar transformation may
•
S
I
m' dt= /(Sm«n n dt — dm n n«St). ' •
seen that
Sm»n:= — ni'Sn n»St^ — Sn-t
(Zin«n = — m-dn n ' dt = — dn • t.
Hence S j m dt=
> i (m • Sn t ' dn — m • dn t • 5n)
S fm dt= '
j
(m X t - Snx dn) = — I n- Snx dn.
424 VECTOR ANALYSIS
The differential 8n X f^n represents the element of area in
the hodogram upon the unit sphere. The integral
/ n«5nx(in=5 / n- da!
surface wtiich lies between the curves G and C". Let the
curve C start at a point upon the surface and spread out to
any desired size. The total amount of turning which is re-
/ S /m«dt=/ m»<£t-'2w.
JJn.da!=.E, ^21)
or S+ I
m ' dt=:2'jr.
•=~n^x.
df'
x^ A sin n t.
D =1 -4 sin n t.
r = Kcos nt + 'Bsinnt.
The discussion of waves may be carried through as pre-
= — » A sin (m
I
• r — » + B sin (m
^) • T —M f (5)
'dt
D
V — IdJi = A cos Cm • r — «. — B sin (m » t — nf)
^
d^
-
^
+ V—1 I
A sin (m • r — n ^) + B sin (m • r — ji }•
^"-"^I7^^^+^~^^^« (6)'
If r = ij + i Ta
and Tj = aj^ i + ^j j + Zj k,
Tj = a;a i + 2/2 j + «2 ^j
or r = a; i -t- y j + a k.
Two bivectora are equal when their real and their imaginary
parts are equal. Two bivectois are parallel when one is the
product of the other by a scalar (real or imaginary). If
a bivector is parallel to a real vector it is said to have a real
direction. In other cases it has a complex or imaginary
direction. The value of the sum, difference, direct, skew,
and indeterminate products of two bivectors is obvioTiS' with-
out special definition. These statements may be put into
analytic form as follows.
Then if r = s, = Sj and t^ — s^
rj
if r Hs r = g = («! + ix^^
j& s,
430 VECTOR ANALYSIS
r +s= (ti + Sj) + i (rg + Sj),
where a^ + b'^ — 1,
the conjugate is multiplied by a — i b, and hence the four
products
160.] A
closer examination of the effect of multiplying a
bivector by a cyclic factor yields interesting and important
geometric results. Let
If a . b = 0,
r • r = (cos 2 y + i sin 2 2') (a • a — b • b).
Let r • r =a+ « 5,
and tan 2 tf = -.
a
a 4- i b = (cos g — i sin 2) r.
r . r = 0, r circular.
If t— a? i + 2/ j 4 « k,
r . r = *!2 + 2/2 + s2 ^ 0.
The condition r • r = 0, which for real vectors implies r = 0,
is not sufficient to ensure the> vanishing of a biveetos. The
HARMONIC VIBRATIONS AND BIVECTORS 433
ellipses are the same. Hence two parallel vectors have their
directional ellipse similar and similarly placed the ratio of —
similitude being the modulus of the biscalar. It is evident
that any two circular bivectors whose planes coincide are
parallel. A circular vector and a non-circular vector cannot
be parallel.
The condition that two bivectors be perpendicular
is r • s = 0,
or '^i
' ^1 — Tg • Sg = Pj • S2 + Tj • Sj := 0.
Tg • 83 = and ij • S2 + Tg • Sj = 0.
28
434 VECTOR ANALYSIS
The first equation shows that t^ and s^ are perpendicular and
hence Sj and Sj are perpendicular. Moreover, the second
shows that the angular directions from r^ to t^ and from Sj to
Sg are the same, and that the axes of the directional ellipses
m = nij + 4 m^,
n= +
n-^ in2,
"-""•
the resultant is (A + B) e "" •
^ g— mj • r gi (— mj • r __
1 Such use of bivectors .is made by Professor Gibbs in his course of lectures on
Bivectors were not used in the second part of this chapter, because in the opinion
of the present author they possess no essential advantage over real vectors until
the more advanced parts of the theory, rotation of the plane of polarization by
magnets and crystals, total and metallic reflection, etc., are reached?