Vector Analysis (PART 1)
Vector Analysis (PART 1)
Introduction
Summary of prerequisites
1. A scalar quatity has magnitude only; a vector quatity has both magnitude and direction.
2. The axes of reference, OX,OY,OZ, form a right-handed set. The symbols I,j,k
denote unit vectors in the direction OX,OY,OZ, respectively.
And l 2 + m2 + n 2 = 1
A.B = a x bx +a y by + a z bz
1
i j k
Also A × B = ax ay az
bx by bz
One or two examples will no doubt help to recall the main points.
If i, j , k are unit vectors in the directions OX,OY,OZ, respectively, then any position
vector OP ( = r ) can be represent in the form
OP = r = a xi +a y j +a z k
Then r =
_______________________________________________________________________
_
r = a x2 + a y2 + a z2
2
The direction of OP is denoted by stating the direction cosines of the angles made by OP
and the three coordinate axes.
OL a x
l = cos α = =
OP r
cos α, a
∴l , m, n, =OM cos β, cos γ
m = cos β =
y
=
OP r
So, if P is the point (3,2,6), then
ON a z
=
r n= cos
γ; = =
OP r
l = ; m = ; n =
_______________________________________________________________________
_
r =7;
l = 0.429 ; m = 0.286 ; n = 0.857
Since ( r )2 = 9 + 4 + 36 = 49 ∴r =7
3
l = cos α = = 0.4286
7
2
m = cos β = = 0.2857
7
6
n = cos γ = = 0.8571
7
If the direction cosines of A are l1 , m1 , n1 and those B are l2 , m2 , n2 , then the angle
between the vectors is given by
3
0 = 66 36 '
For A : r1 = 4 +9 +16 = 29
2 3 4
∴l1 = ; m1 = ;n =
29 29 29
For B : r2 = 1 + 4 +9 = 14
1 −2 3
∴l2 = ; m2 = ; n2 =
14 14 14
1
Then cos 0 = {2 − 6 +12 } = 0.3970
14 × 29
∴0 = 66 36 '
A.B = AB cos 0
If we consider the scalar products of the unit vectors i,j,k which are mutually
perpendicular, then
A.B =
_______________________________________________________________________
_
4
A.B = 2 − 6 + 20 = 16
Also, since A.B = AB cos 0 ,we can determine the angle between the vectors. In this
case 0 =.......... .......... .
_______________________________________________________________________
_
0 = 57 9'
A = 2i − 3 j + 4k ∴A = A = 4 +9 +16 = 29
B = i + 2 j + 5k ∴B = B = 1 +4 +25 = 30
We have already found that A.B =16 and A.B = AB cos 0
i j k
Also A × B = ax ay az
bx by bz
i × j = (1)(1) sin 90 = 1 ∴i × j = j ×k = k ×i =1
5
Note that j ×i = −( i × j ) = −1 ∴ j ×i = k × j = i ×k = −1
Also i × i = (1)(1) sin 0 = 0 ∴i ×i = j × j = k × k = 0
It can also be shown that
( i ) A ×( B + C ) A × B + A ×C
And ( ii ) A × B = −( B × A)
Make a note of these result (3) , ( 4 ) and (5)
Then, if A = 3i − 2 j + 4k and B = 2i − 3 j − 2k
A × B =
_______________________________________________________________________
_
i j k
A ×B = 3 −2 4
2 −3 −2
0 =
_______________________________________________________________________
_
0 = 79 40 '
6
B = B = 2 2 +32 + 2 2 = 17 = 4.123
∴21 .84 =(5.385 )( 4.123 ) sin 0
∴sin 0 = 0.9838 ∴0 = 79 40 '
So, to recapitulate:
Triple products
We now deal with the various products that we form with three vectors.
Scalar triple product of three vectors.
If A,B,C are three vectors, the scalar formed by the product A.( B ×C ) is called the
scalar triple product.
i j k
∴A.( B ×C ) = (a i + a j + a k ). bx by bz
cx cy cz
Multiplying the top row by the external bracket and remembering that
Example If A = 2i − 3 j + 4k ; B = i − 2 j − 3k ; C = 2i + j + 2k
7
2 −3 4
Then A.( B ×C ) = 1 −2 −3
2 1 2
=
_______________________________________________________________________
_
A.( B ×C ) = 42
2 −3 4
For A.( B ×C ) = 1 −2 −3
2 1 2
= 2( − 4 + 3) + 3( 2 + 6 ) + 4(1 + 4 ) = 42
As simple as that.
bx by bz ax ay az
(a) B.( C × A) = cx cy c z = − cx cy cz
ax ay az bx by bz
Since interchanging two rows in a determinant reverses the sign. If we now interchange
rows 2 and 3 and again change the sign, we have
ax ay az
B.( C × A) = bx by bz = A.( B ×C )
cx cy cz
i.e. the scalar triple product is unchanged by a cyclic change of the vectors involved
bx by bz ax ay az
(b) B.( A ×C ) = a x ay a z = − bx by bz
cx cy cz cx cy cz
i.e. a change of vectors not in cyclic order, changes the sign of the scalar triple product.
8
ax ay az
(c) A.( B × A) = bx by bz = 0 since two rows are identical
ax ay az
Example If A = i + 2 j + 3k ; B = 2i − 3 j + k ; C = 3i + j − 2k
1 2 3
For A.( B ×C ) = 2 −3 1 =1( 6 −1) −2( −4 −3) +3( 2 +9 ) = 52
3 1 −2
C.( B × A) = −A.( B ×C ) = − 52
Coplanar vectors
The scalar triple products provides a test to show whether three given vectors lie in the
same plane.
9
The scalar product of two vectors A and D is A.D where A.D =A D cos 0,0
being the angle between A and D
π π
If 0 = , A.D = A D cos =0
2 2
∴A.D = 0
10
_______________________________________________________________________
_
A.( B ×C ) = 0
1 2 −3
For A.( B ×C ) = 2 −1 2 =1(1 −2 ) −2( −2 −6 ) −3( 2 +3) = 0
3 1 −1
Example 2
If A = 2i − j + 3k ; B = 3i + 2 j + k ; C = i + pj + 4k
are coplanar, find the value of p .
The method is clear enough. We merely set up and evaluate the determinant and solve the
equation A.( B ×C ) = 0
_______________________________________________________________________
_
p =
p = −3
3 2 −1
In this case A.( B ×C ) = 2 −1 3
1 −2 2
∴ A.( B × C ) =/ 0
∴A, B, C are not coplanar
11
_______________________________________________________________________
_
A × ( B ×C )
And are called the vector triple products (10 )
( A × B ) ×C
Now
i j k
by bz b bz bx by
( B ×C ) = bx bx bx = i −j x +k
cy cz cx cz cx cy
cx cy cz
i j k
Then A ×( B ×C ) = ax ay az
by bz b bz bx by
− x
cy cz cx cz cx cy
i j k
= ax ay az
by bz bz bx bx by
cy cz cz cx cx cy
Example 1 If A = 2i − 3 j + k ; B = i + 2 j − k ; C = 3i + j + 3k ;
Determine the vector triple product A ×( B ×C )
12
B × C = 7i − 6 j − 5k
i j k
=i ( 6 +1) − j (3 +3) +k (1 −6 )
For B ×C = 1 2 −1
= 7i −6 j −5k
3 1 3
Then A × ( B × C ) =
_______________________________________________________________________
_
A × ( B × C ) = 2li +17 j + 9k
i j k
For A ×( B ×C ) = 2 −3 1
7 −6 5
That is fundamental enough. There is, however, an even easier way of determining a
vector triple product. It can be proved that
A ×( B ×C ) = ( A.C ) B − A.B )C 11
And ( A × B ) ×C = ( C. A) B − ( C.B ) A
The proof of this is given in the Appendix. For the moment, make a careful note of the
expressions: then we will apply the method to the example we have just completed.
_______________________________________________________________________
_
A = 2i − 3 j + k ; B = i + 2 j − k ; C = 3i + j + 3k and we have
A ×( BXC ) = ( A.C ) B − ( A.B )C
= ( 6 − 3 + 3)( i + 2 j − k ) − ( 2 − 6 −1)( 3i + j + 3k )
= 6( i + 2 j − k ) + 5( 3i + j + 3k ) = 2li +17 j + 9k
13
Here is another.
Example 2 If A = 3i + 2 j − 2k ; B = 4i − j + 3k C = 2i − 3 j + k determine ( A × B ) ×C
using the relationship.
( A × B ) ×C = ( C. A) B − ( C.BA )
( A × B ) ×C =
_______________________________________________________________________
_
−50 i − 26 j + 22 k
For ( A × B ) ×C = ( C. A) B − ( C.B ) A
= ( 6 − 6 − 2 )( 4i − j + 3k ) −8 + 3 + 3( 3i + 2 j − 2k )
= −2( 4i − j + 3k ) −14 ( 3i + 2 j − 2k )
= − 50 i − 26 j + 22 k
Example 3 A = i + 3 j + 2k ; B = 2i + 5 j − k ; C = i + 2 j + 3k
A × ( B × C ) =
( A × B ) × C =
And ( A × B ) ×C = ( C. A) B − ( C.B ) A
= (1 + 6 + 6 )( 2i + 5 j − k ) − ( 2 +10 − 3)( i + 3 j + 2k )
=13 ( 2i + 5 j − k ) − 9( i + 3 j + 2k ) =17 i + 38 j − 31k
14
(i )(i × j ) = k
∴i ×( i × j ) = i ×k = − j
∴i ×( i × j ) − j
(ii )( i ×i ) × j = ( 0 ) × j = 0
∴( i ×i ) × j = 0
_______________________________________________________________________
_
( a ) A.( B ×C ) = −12
(b ) A ×( B ×C ) = 62 i + 44 j − 74 k
( A × B ) × C =109 i + 7 j − 22 k
Here is the working.
5 −2 3
(a) A.( B ×C ) = 3 1 −2
1 −3 4
15
= 5( 4 − 6 ) + 2(12 + 2 ) + 3( − 9 −1) = −12
(ii) ( A × B ) ×C = ( C. A) B − ( C.B ) A
= 23( 3i + j − 2k ) − ( − 8)( 5i − 2 j + 3k )
=109 i +7 j −22 k
Differential of vectors
In many practical problems, we often deal with vectors that change with time, e.g.
velocity, acceleration, etc. If a vector A depends on a scalar variable t, then A can be
represented as A(t ) and A is then said to be a fuction of t .
If A = a xi + a y j +a z k then ax , a y , a z will also be dependent on the parameter
t.
i.e. A( t ) = a x ( t )i + a y ( t ) j + a z ( t ) k
_______________________________________________________________________
_
d
{ A( t )} = i d { ax ( t )} + j d d
a y ( t ) + k { az ( t )}
dt dt dt dt
dA da da da
In short =i x + j y +k z .
dt dt dt dt
dA da da da
=i x + j y +k z
du du du du
16
Example 1 If A = (3u 2 + 4)i + ( 2u − 5) j + 4u 3k , then
dA
=.......... ......
du
_______________________________________________________________________
_
dA
= 6ui + 2 j + 12 u 2 k
du
d2A
If we differentiate this again, we get = 6i + 24uk
du 2
dA d2A
When u = 2 = 12 i + 2 j + 48 k and = 6i + 48uk
du du 2
Then
dA d2A
=.......... .. and =.......
du du 2
dA d2A
=49 .52 ; and = 48 .37
du du 2
For
dA
du
{
= 12 2 + 2 2 + 48 2 }
12
= {2452 }
12
= 49 .52
17
And
d2A
du 2
{
= 62 + 48 2 }12
= {2340 }
12
= 48 .37
_______________________________________________________________________
_
dF
= 2 cos 2i + 3e3 j − k
dt
d 2F
2
= −4 sin 2i + 9e3 j + 6k
dt
dF d 2F
From these, we could if required find the magnitudes of and
dt dt 2
dF d 2F
=...... =.......... ...
dt dt 2
_______________________________________________________________________
_
dF d 2F
= 60 .27 ; =180 .9
dt dt 2
And
d 2F
dt 2
{
= ( − 4 sin 2 ) + 81 e 6 + 36
2
} 12
dA d2A dA d2A
(a) (b) (c) (d) at u =3
du du 2 du du 2
18
_______________________________________________________________________
_
(a) dA = i − 2uj + 6u 2 k at u = 3, dA = i − 6 j + 54 k
du du
2 2
(b) d A = −2 j + 12uk at u = 3, d A2
= −2 j + 36 k
du 2 du
2
(d) d A2
= {4 +1296 } = (1300 ) = 36 .06
12 12
du
_______________________________________________________________________
_
r = ( 2t + 3) i + ( t 2 + 3t ) j + (t 3 + 2t 2 )k
Then, at t =1
D = ( B × C) d 2r
dt 2
= ........
_______________________________________________________________________
_
dr d 2r
= 2i + 5 j + 7 k ; = 2 j + 10 k
dt dt 2
For
dr
dt
(
= 2i + ( 2t + 3) j + 3t 2 + 4t k)
dr
∴at t =1, dt
= 2i + 5 j + 7k
19
d 2r
And 2 = 2 j + ( 6t + 4) k
dt
d 2r
∴At t =1,
dt 2
= 2 j + 10 k
2i + 3 j + 4k 1
= ( 2i + 3 j + 4k )
4 + 9 + 16 29
dr
Denote this unit vector by I. Then the component of in the direction of I
dt
d
r d
r
= c
os 0 = .I
dt dt
1
= (2i +5 j +7 k )(
. 2i +3 j +4k )
29
=.......... .......... ..
_______________________________________________________________________
_
8.73
1
Since ( 2i + 5 j + 7k ).( 2i + 3 j + 4k ) = 1 ( 4 + 15 + 28 )
29 29
47
= = 8.73
29
d 2r
Similarly, the component of in the direction of I is
..........
dt 2
..........
_______________________________________________________________________
_
8.54
d 2r d 2r 1
For 2
cos 0 = 2
.I = ( 2 j +10 k ).( 2i + 3 j + 4k )
dt dt 29
20
1 46
= (6 + 40 ) = = 8.54
29 29
d
(b) { A + B} = dA + dB
du du du
d
(c) { A.B} = A. dB + dA .B
du du du
d
(d) { A × B} = A × dB + dA × B.
du du du
A(u ). A( u ) = a x2 + a y2 + a z2 = A = A2
2
i.e. constant
d
Also { A( u ). A( u )} = A( u ). d { A( u )} + A( u ). d { A( u )}
du du du
= 2 A( u ).
d
{ A( u )} = d { A2 } = 0
du du
d d
A(u ). { A(u )} = A(u ) { A(u )} cos 0 = 0
du du
∴cos 0 =0
π
∴0 =
2
Now let us deal with unit tangent vectors.
21
Unit tangent vectors
_______________________________________________________________________
_
d
{ A( u )}
T = du
d
{ A( u )}
du
Example 1 Determine the unit tangent at the point (2,4,7) for the curve with parametric
equations x = 2u; y =u 2 + 3; z = 2u 2 + 5
( ) ( )
r = a x i + a y j + a z k = 2ui + u 2 + 3 j + 2u 2 + 5 k
dr
∴ =.......... .
du
_______________________________________________________________________
_
22
dr
= 2i + 2uj + 4uk
du
dr
And at u =1, = 2i + 2 j + 4k
du
dr
Hence du = T =
_______________________________________________________________________
_
dr 1
=2 6 ; T= { i + j + 2k }
du 6
dr
For du = {4 + 4 +16 } = 24 = 2 6
12 12
dr
2i + 2 j + 4k 1
T= du = = { i + j + 2k }
dr 2 6 6
du
Let us do another.
Example 2 Find the unit tangent vector at the point ( 2,0, π ) for the curve with parametric
equations x =2 sin 0; y =3 cos 0; z = 20.
Then, at dr
0 =π 2, =.......... ...
d0
dr
=
d0
23
dr dr
= −3 j + 2k ; = 13
d0 d0
1
T = ( − 3 j + 2k )
13
And now
x = 3t ; y = 2t 2 ; z = t 2 + t
1
T= ( 3i + 8 j + 5k )
98
At t = 2, r = 6i + 8 j + 6k and dr = 3i + 8 j + 5k
dt
dr
= ( 9 + 64 + 25 ) = 98
12
∴
dt
dr dt 1
∴T = = ( 3i + 8 j + 5k )
dr dt 98
_______________________________________________________________________
_
24
If F = xi + yj + zk then x,y,z will also be function of u and v
∂F ∂x ∂y ∂z
Then = i+ j+ k
∂u ∂u ∂u ∂u
∂F ∂x ∂y ∂z
= i+ j+ k
∂v ∂v ∂v ∂v
∂2 F ∂2 x ∂2 y ∂2 z
= 2i+ 2 j+ 2k
∂u 2 ∂u ∂u ∂u
∂ 2 F ∂2 x ∂2 y ∂2 z
= i + j + k
∂v 2 ∂v 2 ∂v 2 ∂v 2
∂2 F ∂2 x ∂2 y ∂2 z
= i+ j+ k
∂u∂v ∂u∂v ∂u∂v ∂u∂v
Example If F = 2uvi + (u 2 − 2v ) j + (u + v 2 )k
∂F ∂F ∂2 F ∂2 F
= ; = ; 2 = ; =
∂u ∂v ∂u ∂u∂v
∂F ∂F
= 2vi + 2uj + k ; = 2ui − 2 j + 2vk
∂u ∂v
∂2 F ∂2 F
= 2 j ; = 2i
∂u 2 ∂u∂v
25
b b b b
∫ Fdu
a
= i ∫ xdu + j ∫ ydu + k ∫ zdu
a a a
42 i − 2 j + 80 k
[
Fdt = i (t 3 + 2t 2 ) + j (t 2 −5t ) + kt 4 ]
3 3
For ∫1 1
= ( 45 i − 6 j + 81 k ) − ( 3i − 4 j + k ) = 42 i − 2 j + 80 k
Here is slightly different one.
Example 2 If F = 3ui + u 2 j + ( u + 2) k
∫ ( F ×V )du
2
Evaluate 0
F ×V = (u 3 + u 2 + 6u )i − (u 2 − 10u ) j − ( 2u 3 + 9u 2 )k
i j k
For F ×V = 3u u 2
(u + 2)
2u −3u (u − 2)
∫ ( F ×V )du
2
Then =.......... .......
0
4
{14 i +13 j − 24 k }
3
u 4
u3 u3 u4
∫ ( F ×V )du = 4
+
3
+ 3u 2 i − − 5u 2 j −
3 2
+ 3u 3 k
26
8 8
∴∫ ( F ×V )du = 4 + + 12 i − − 20 j − ( 8 + 24 ) k
2
0
3 3
4
= {14 i +13 j − 24 k }
3
Example 3 If F = A × ( B × C ) where
A = 3t 2i + ( 2t − 3) j + 4tk
B = 2i + 4tj + 3(1 − t ) k
C = 2ti −3t 2 j − 2tk
1
Determine ∫ Fdt
0
A ×( B ×C ) − ( A.C ) B − ( A.B )C
So A.C =
And A.B =
1
{− 20 i +132 j + 75 k }
60
For F = A × B × C = ( 2t 2 − 4t 3 )i + ( 4t 3 + 6t 4 ) j + (3t 2 + t 3 )k
Intergration with respect to t then gives the result stated above.
27
Scalar and vector fields
28
∂o
/ ∂o
/ ∂o
/
/=
grad o i+ j+ k
∂x ∂y ∂z
/ = 4 xyz i +2 x z j +6 x yz k
3 2 3 2 2
grad o
∂ ∂ ∂
/ = i
grad o +j +k o
/
∂x ∂y ∂z
∂ ∂ ∂
Where is called a vector differential operator and is denoted by
i ∂x + j ∂y + k ∂z
the symbol ∇ (pronounced ‘del’or sometimes ‘nabla’)
∂ ∂ ∂
i.e. ∇=
i ∂x + j ∂y + k ∂z
o/
Beware! ∇cannot exist alone:it is an operator and must operate on a stated scalar
function o/( x, y, z )
If F is a vector function, ∇F has no meaning.
So we have:
∂ ∂ ∂
∇o / =
/ = grad o i ∂x + j ∂y +k ∂z /
o
∂o
/ ∂o
/ ∂o
/
=i +j +k
∂x ∂y ∂z
4( 2li + 8 j +18 k )
∂o
Since o/ = x 2 yz 3 + xy 2 z 2 ∴ / = 2 xyz 3 + y 2 z 2
∂x
∂o/ ∂o/
= x 2 z 3 + 2xyz 2
= 3x 2 yz 2 + 2 xy 2 z
∂y ∂z
29
∂o/ ∂o
Then at (1,3,2 ) = 48 + 36 ∴ / = 84
∂x ∂x
∂o
∴ / = 32
∂o
/
= 8 + 24
∂y ∂y
∂o/ ∂o
= 36 + 36 ∴ / = 72
∂z ∂z
A.B = x 2 yz 3 + x 2 yz + x 2 y 2 z 2
2 xyz ( z 2 +1 + yz )i + x 2 z ( z 2 +1 + 2 yz ) j + x 2 y (3z 2 +1 + 2 yz )k
( ) ( ) (
∴∇( A.B ) = 2 xyz z 2 +1 + yz i + x 2 z z 2 +1 + 2 yz j + x 2 y 3 z 2 +1 + 2 yz k )
Now let us obtain another useful relationship.
30
If o/( x, y , z ) is a scalar function at P , we know
∂o
/ ∂o
/ ∂o
/
/ =∇o
grad o /= i+ j+ k
∂x ∂y ∂z
Then / dr =
grad o.
∂o
/ ∂x + ∂o
/ ∂y + ∂o
/ ∂z
/.dr =
grad o
∂x ∂y ∂z
∂o/ ∂o ∂o
/ j + / k .( dxi + dyj + dzk )
For grad o/.dr = i +
∂x ∂y ∂z
∂o ∂o ∂o
= / dx + / ∂y + / ∂z
∂x ∂y ∂z
= the total different d o/ /
of o
∴ /=
do dr . grad o
/
Directional derivatives
dr dr
=
ds dr
dr
and is thus a unit vector in the direction
ds
of dr .
If we denote this unti vector by â, i.e
dr
= â, the result becomes
ds
do/ =a
ˆ.grad o
/
ds
do/
/ on the unit vector â and is called the directional
is thus the projection of grad o
ds
derivative of o
/ in the direction of â . It gives the rate of change of o
/ with distance
do/
measured in the direction of â and =a ˆ.grad o / be a maximum when â and
ds
/ have the same direction of ,since then
grad o
31
aˆ.grado/ = aˆ grado/ cos0 and 0 will be zero
/ gives the direction in which the maximum rate of changes
Thus the direction of grad o
of / occurs
o
∇o/ = ( 2 xz + 2 y 2 )i + ( 4 xy + z 2 ) j + ( x 2 + 2 yz )k
∂o/ ∂o
/ = 4 xy + z 2 ; ∂o/
Since = 2 xz + 2 y 2 ; = x 2 + 2 yz
∂x ∂y ∂z
Then at (1,2,−1)
∇o / = ( − 2 +8)i + (8 +1) j + (1 − 4 )k = 6i +9 j −3k
1
aˆ = ( 2i + 3 j − 4k )
29
For A = 2i + 3 j − 4k ∴A = 4 +9 +16 = 29
A 1
aˆ = = ( 2i + 3 j − 4k )
A 29
1
So we have ∇o/ = 6i + 9 j − 3k and aˆ = ( 2i + 3 j − 4k )
29
do/
∴ = aˆ.∇o
/ =.......... ..
ds
do/ 51
= = 9.47
ds 29
do/ 1
Since = aˆ.∇o/ = ( 2i + 3 j − 4k ).( 6i + 9 j − 3k ) )
ds 29
1
= (12 + 27 +12 ) = 51 = 9.47
29 29
32
That is all there is to it
(i) from the given scalar function o/ , determine ∇ /
o
(ii)Find the unit vector â in the direction of the given vector A
/ = aˆ.∇o
do
(iii) Then /
ds
do/ − 23
= = −3.43
ds 3 5
For o/ = x 2 y + y 2 z + z 2 x
( ) ( ) (
∴∇o/ = 2 xy + z 2 i + x 2 + 2 yz j + y 2 + 2 zx k )
∴ at (1−, 1,2) , ∇o/ = 2i − 3 j + 5k
A = 4i + 2 j − 5k ∴A = 16 +4 +25 = 45 =3 5
1
∴aˆ = ( 4i + 2 j − 5k )
3 5
do/ 1
∴ = aˆ.∇o/ = ( 4i + 2 j − 5k ).( 2i − 3 j + 5k )
ds 3 5
=
1
( 8 − 6 − 25 ) = − 23 = − 3.43
3 5 3 5
Example 3. Find the direction from the point 1,1,0 which gives the greatest rate of
increase rate of increase of the function o/ = ( x + 3 y ) 2 + ( 2 y − z ) 2
/
This appears to be different, but it rests on the fact that the greatest rate of increase of o
with respect to distance is in…..
The direction of
∇ /
o
All we need then is to find the vector ∇/ , which is
o
∇o/ = 4( 2i + 8 j − k )
For o/ = ( x + 3 y ) 2 + ( 2 y − z ) 2
∂o ∂o/ ∂o/
∴ / = 2( x + 3 y ); = 6( x + 3 y ) + 4( 2 y − z ); = −2( 2 y − z )
∂x ∂y ∂z
∂o/ ∂o/ =32 ; ∂o /
∴At (1,1,0), ∂x
=8;
∂y ∂z
= −4
∴∇o/ = 8i + 32 j − 4k = 4( 2i + 8 j − k )
33
Unit normal vectors
If o/( x, y, z ) = constant, this relationship represents a surface in space, depending on
the value ascribed to the constant.
( ) ( ) (
∇o/ = 3 x 2 y + 4 z 2 + y 2 z i + x 3 + 2 xyz j + 8 xz + xy 2 )
Then, at (1,3,−1) ∇o/ = 4i − 5 j + k
And the unit normal at (1,3,−1) is
1
( 4i − 5 j + k )
42
Since ∇o
/ = 16 +25 +1 = 42
/ ∇o 1
And N = ∇o/ = 42 ( 4i − 5 j + k )
34
One more
Example 2 Determine the unit normal to the surface
xyz + x 2 y − 5 yz − 5 = 0 at the point ( 3,1,2 )
All very straightforward. Complete it.
1
Unit normal = N = ( 8i + 5 j − 2k )
93
For o/ = xyz + x 2 y − 5 yz − 5
( )
∴∇o/ = ( yz + 2 xy )i + xz + x 2 − 5 z j + ( xy − 5 y ) k
At (3,1,2), ∇o/ =8i +5 j −2k ; ∇o/ = 64 +25 +4 = 93
∇o 1
∴unit normal = N = ∇o// =
93
(8i + 5 j − 2k )
Collecting our results so far, we have, for o/( x, y , z ) a scalar function
∂o
/ / ∂o / ∂o
(a) grad o/ = ∇o/ = ∂x i + ∂y j + ∂z k
/ / ∂o
/ ∂o ∂o
/ =dr .grad o
(b) do /=
/ where do ∂x + ∂y + ∂z
∂x ∂y ∂z
do/
(c) directional derivative =a
ˆ.grad o
/
ds
/ ∇o
(d)unit normal vector N = ∇o/
(a) ∇( A + B ) = i ∂ ∂ ∂
( A + B) + j ( A + B) + k + ( A + B )
∂x ∂y ∂z
∂A ∂A ∂A ∂B ∂B ∂B
= i + j+ k + i+ j+ k
∂x ∂y ∂z ∂x ∂y ∂z
(b) ∇( AB ) = i ∂ ∂ ∂
( AB ) + j ( AB ) + k ( AB )
∂x ∂y ∂z
∂B ∂A ∂B ∂A ∂B ∂A
= i A +B + j A +B + k A +B
∂x ∂x ∂y ∂y ∂z ∂z
∂B ∂B ∂B ∂A ∂A ∂A
= A i +A j +A k + B i +B j +B k
∂x ∂y ∂z ∂x ∂y ∂z
∂B ∂B ∂B ∂A ∂A ∂A
= A i+ j+ k + B i + j+ k
∂x ∂y ∂z ∂x ∂y ∂z
35
∴∇( AB ) = A(∇B ) +B (∇A)
Remember that in these results A and B are scalars. The operator ∇acting on a vector
Has no meaning
∂B ∂B ∂B
∇B =
∂x
i+
∂y
j+
∂z
(
k = y 2 zi + 2 xyzj + xy 2 − 3 z 2 k )
=3i +12 j −25 k at ( 2,1,3)
∇( AB ) = 3( −117 i + 36 j − 362 k )
∇( AB ) = A(∇B ) + B ( ∇A)
A = x yz + xz
2 2
∴at ( 2,1,3), A = 12 + 18 = 30
B = xy z − z 3
2
∴at ( 2,1,3), B = 6 − 27 = −21
∴∇( AB ) = 30 ( 3i + 12 j − 25 k − 21( 21i +12 j + 16 k ) )
= −351 i +108 j −1086 k =3( −117 i +36 j −362 k )
The operator ∇. (notice the ‘dot’; it makes all the difference) can be applied to a vector
function A( x, y, z ) to give the divergence of A, written in short as div A
If A = axi + a y j + az k
36
∂ ∂ ∂
A = ∇. A = .(a xi + a y j + a z k )
i ∂x + j ∂y + k ∂z
Div
∂a x ∂a ∂a
∴divA = ∇. A = + y + z
∂x ∂y ∂z
Note that
(a) The grad operator ∇acts on a scalar and gives a vector
(b) The div operation ∇. acts on a vector and gives a scalar
divA = ∇. A = 2 xy − xz + 2 yz
We simply take the appropraite partial derivatives of the coefficients of i,j and k. It
couldhardly be easier.
∇A
. =0
For ( )
A = 2 x 2 yi − 2 xy 2 + y 3 z j + 3 y 2 z 2 k
∂a x ∂a y ∂a z
∇. A = + +
∂x ∂y ∂z
( )
= 4 xy − 2 2 xy + 3 y 2 z + 6 y 2 z
= 4 xy − 4 xy − 6 y z + 6 y 2 z = 0
2
Such a vector A for which ∇A . = 0 at all points, i.e for all values of x, y , z is called a
solenoid vector. It is rather a special case.
The curl operator denoted by ∇x ,acts on a vector and gives another vector as a result.
∂ ∂ ∂
i.e. curl A = ∇× A =
i +j +k ×( axi + a y j + az k )
∂x ∂y ∂z
i j k
∂ ∂ ∂
=
∂x ∂y ∂z
ax ay az
37
∂a z ∂a y ∂a x ∂a z ∂a y ∂a x
∴∇× A = i
∂y − ∂z
+ j ∂z − ∂x + k
∂x − ∂y
Curl A is thus a vector function. It is best remembered in its determinant form, so make a
note of it.
_______________________________________________________________________
_
∂ ∂ 2 ∂ ∂ 4
∇× A = i (
− x 2 yz − ) (
x + y2 − j )
− x 2 yz − (
y − x2 z 2 ) ( )
∂y ∂z ∂x ∂z
All that now remains is to obtain the partial derivatives and substitute the values of
x, y , z
2i − 8 j −106 k
{ } { } {
∇ × A = i − x 2 z − j − 2 xyz + 2 x 2 z + k 2 x − 4 y 3 }
∴At (1,3,−2), ∇× A = i ( 2 ) − j{12 − 4} + k ( 2 −108 )
=2i −8 j −106 k
_______________________________________________________________________
_
38
i j k
∂ ∂ ∂
∂x ∂y ∂z
ze 2 xy 2 xz cos y x +2 y
Now expand the determinant and substitutethe values of x,y, and z, finally obtaining curl
F =.......... .......... ...
_______________________________________________________________________
_
Curl F = ∇× F = −2( i + 3k )
{ } {
∇ × F = i{ 2 − 2 x cos y} − j 1 − e 2 xy + k 2 z cos y − 2 xze 2 xy }
∴at ( 2,0,3) ∇× F = i ( 2 − 4 ) − j (1 −1) + k ( 6 −12 )
= −2i − 6k = −2( i + 3k )
i j k
∂ ∂ ∂
(ii) curl A = ∇× A =
∂x ∂y ∂z
ax ay az
_______________________________________________________________________
_
Example 3 If o/ = x 2 y 2 + x 3 yz − yz 2
And F = xy 2i − 2 yzj + xyzk
39
Determine for the point P (1,−1,2 )
(a) ∇ /,
o
(b)unit normal,
(c) ∇ .F
(d) ∇×F .
_______________________________________________________________________
_
/∂o / ∂o / ∂o
(a) ∇o/ = ∂x i + ∂y j + ∂z k
( ) ( ) ( )
= 2 xy 2 + 3 x 2 yz i + 2 x 2 y + x 3 z − z 2 j + x 3 y − 2 yz k
∴at (1,−1,2) ∇o/ = −4i − 4 j + 3k
∇o/
(b) N = ∇o
/ = 16 +16 +9 = 41
∇o/
−1
∴N = ( 4i + 4 j − 3k )
41
∂a x ∂a ∂a
(c) F = xy 2i − 2 yzj + xyzk ∇.F = + y + z
∂x ∂y ∂z
∴∇.F = y 2 −2 z + xy
∴at (1,−1,2 ) ∇.F = 1 − 4 −1 = −4 ∴∇.F = −4
i j k
∂ ∂ ∂
(d) ∇× F =
∂x ∂y ∂z
xy 2 − 2 yz xyz
∴∇× F = i ( xz + 2 y ) − j ( yz − 0) + k ( 0 − 2 xy )
= ( xz + 2 y )i − yzj − 2 xyk
Multiple operations
We can combine the operators grad, div and curl in multiple operations, as in the
examples that follow.
40
Example 1 If A = x 2 yi + yz 3 j − zx 3k
∂ ∂ ∂ 2
Then div A = ∇. A =
i+ j+
∂z
(
k . x yi + yz 3 j − zx 3k )
∂x ∂y
= 2 xy + z + x = o/ say
3 3
( )
i.e. grad div A = ∇( ∇. A) = 2 y + 3 x 2 i + 2 xj + 3 z 2 k
_______________________________________________________________________
_
div grad o/ = 6
We have o/ − xyz − 2 y 2 z + x 2 z 2
∂o/ ∂o
/ ∂o
/
Grad / = ∇o/ =
o i+ j+ k
∂x ∂y ∂z
= ( yz + 2 xz 2 )i + ( xz − 4 yz ) j + ( xy − 2 y 2 + 2 x 2 z )k
Determine an expression for curl F in the usual way, which will be a vector, and then the
curl of the result. Finally substitute values
_______________________________________________________________________
_
Curl curl F = ∇× ( ∇× F ) = i + 2 j + 6k
41
i j k
∂ ∂ ∂
For curl F =
∂x ∂y ∂z
x 2 yz xyz 2 y2z
(
= ( 2 yz − 2 xyz )i + x 2 yj + yz 2 − x 2 z k)
i j k
∂ ∂ ∂
Then curl curl F =
∂x ∂y ∂z
2 yz − 2 xyz x2 y yz 2 − x 2 z
= z 2i − ( − 2 xz − 2 y + 2 xy ) j + ( 2 xy − 2 z + 2 xz ) k
_______________________________________________________________________
_
Remember that grad, div and curl are operators and that they must act on a scalar or
vector as appropriate. They cannot exist alone and must be followed by a function.
One or two interesting general results appear.
∂o
/ /∂o / ∂o
grad o/ = ∂x i + ∂y j + ∂z k
i j k
∂ ∂ ∂
∴curl grad /=
o
∂x ∂y ∂z
∂o
/ ∂o/ ∂o
/
∂x ∂y ∂z
∂2 o / − j ∂ o
∂2 o 2
/ − ∂o
2
/ + k ∂ o
2
/ − ∂o
2
/ =0
=i / −
∂y∂ z ∂z∂ y ∂z∂ x ∂x∂ z ∂x∂ y ∂y∂ x
∴curl grad / =∇×(∇o
o / ) =0
i j k
∂ ∂ ∂
curl A = ∇× A =
∂x ∂y ∂z
ax ay az
42
∂a ∂a ∂a ∂a ∂a ∂a
= i z − y − j z − x + k y − x
∂y ∂z ∂x ∂z ∂x ∂y
∂ ∂ ∂
Then div curl A = ∇.( ∇× A) =
i +j +k .( ∇× A)
∂x ∂y ∂z
∂2 a z ∂2 a y ∂2 a z ∂2 a x ∂2 a y ∂2 a x
= − − + + − =0
∂x∂y ∂z∂x ∂x∂y ∂y∂z ∂z∂x ∂y∂z
/ where o
(c) Div grad o / is a scalar
/ ∂o / ∂o / ∂o
grad o/ = ∂x i + ∂y j + ∂z k
43
That brings us to the end of this particular programme. We have covered quite a lot of
new material, so check carefully through the Revision Summary that follows; then you
can deal with the test exercise.
44