Chapter 1. Vector Analysis
1.1. Introduction
This chapter will give an introduction to vector analysis. A
vector, or displacement, has a
direction and a
magnitude.
The following vector notations will be used in this
course:
1.

:
the vector
A 2.

:
the magnitude of vector
AThe
opposite of a vector

is a vector with the same magnitude as

but pointing in a direction opposite to that of

.
The opposite of vector

is written as

.
There
are four different vector operations that we will be using in this course:
vector addition, vector multiplication by a scalar, the dot product of two
vectors, and the cross product of two vectors. We will start Chapter 1 by
discussing these four vector operations in some detail.
1.1.1. Vector Addition
Two vectors

and

can be added. The result of this operation is a new vector

(see Figure 1.1a). Using vector notation we can write vector addition as
follows:
Vector addition is commutative. This means that the order in which
two vectors are added does not affect the result (see Figure 1.1). The
commutative properties of vector addition can be written as
To subtract vector

from vector

is equivalent to adding the opposite of vector

to

.
In other words:

Figure
1.1. Vector Addition.
1.1.2. Multiplication by a scalar
A vector can be multiplied by a scalar. If the scalar
a is a
positive number (see Figure 1.2a) than the result of the multiplication of the
vector

by
a is a new vector with a magnitude equal to

and a direction equal to the direction of

.
If the scalar
a is a negative number (see Figure 1.2b) than the result of
the multiplication of the vector

by
a is a new vector with a magnitude equal to

and a direction opposite to the direction of

.
Figure
1.2. Vector multiplication.
Scalar multiplication is
distributive. This means that the result
of the multiplication of the vector sum of vector

and

by a scalar
a is equal to the vector sum of

and

:
1.1.3. Dot product of two vectors
Figure
1.3. The scalar product of two vectors.
The dot product, also called the scalar product, is defined
as
where

is the angle between vectors

and

(see Figure 1.3). The dot product is
commutative which means that the
order of the vectors

and

does not effect the result of the dot product. In other words:
The dot product is also distributive:
The dot product will be frequently used to determine whether two vectors

and

are perpendicular. When

and

are perpendicular the angle

is equal to 90°. The definition of the dot product shows that in this case
the dot product between

and

is equal to zero.
1.1.4. Cross products of two vectors
The cross product of two vectors

and

is a third vector

.
The vector

is perpendicular to both

and

,
and has a length equal to
where

is the smallest angle between

and

(see Figure 1.4). The direction of the vector

can be determined using the
right-hand rule. By applying the right-hand
rule it can be shown easily that
the vector product is not commutative.
This implies that the order of the vectors

and

is important, and reversing the order will change the result of the vector
product:
The vector product is distributive which requires that

Figure
1.4. The vector product of
and
.
Example: Problem 1.2
Is the cross product associative? If
so, prove it; if not, provide a counter example.
If the cross product of
two vectors is associative then the following relation must hold:
Consider the special case in which

=

and

is perpendicular to

and

.
In this case the cross product between

and

is equal to the null vector. The left-hand side of the equation is thus equal
to the null vector. The cross product of

and

is a new vector, perpendicular to both

and

,
and with a length equal to

.
The cross product between

and

is a new vector. Since

and

are perpendicular, the magnitude of the cross product of

and

is equal to

which is non-zero. We thus conclude that

is not equal to

which shows that the cross product is not associative.
1.2. Vector components
A vector in three dimensions can be identified uniquely by specifying
three coordinates. In a
Cartesian coordinate frame three base vectors

,

,
and

are defined, parallel to the
x,
y, and
z axes, respectively
(see Figure 1.5). Each of these
base vectors has unit length and is
perpendicular to the other two base vectors. Any vector is uniquely defined by
specifying its components along the
x,
y, and
z axes. A
vector

is defined in terms of the three Cartesian coordinates
Ax,
Ay, and
Az. Using the three base vectors,
one can reconstruct the vector

:

Figure
1.5. The Cartesian coordinate frame.
The vector operations discussed in Section 1.1 can be easily expressed in
terms of vector components. We will now discuss each of these four vector
operations.
1.2.1. Vector addition
The components of the vector sum of vectors

and

is equal to the sum of their like components. If a vector

is the vector sum of the vectors

and

then the components of

are equal to
In these equations,
Ax,
Ay, and
Az are the vector components of vector

,
and
Bx,
By, and
Bz are the
vector components of vector

.
1.2.2. Multiplication by a scalar
The components of a vector

,
which is the result of a scalar multiplication of vector

with a scalar
a, are equal to the components of

multiplied by
a:
1.2.3. Scalar product
The scalar product of vectors

and

is equal to the sum of the products of their like components:
This relation can be derived using the commutative properties of the scalar
product:
In the last step of this derivation we have used the fact that the vectors

,

,
and

are perpendicular to each other and have unit length. The scalar products of
these unit vectors are easy to evaluate and are equal to
1.2.4. Vector Product
The components of the vector

,
which is the vector product of vectors

and

,
can be calculated using the distributive properties of the vector
product:
In this derivation we have used the fact that the vectors

,

,
and

are perpendicular to each other and have unit length. The vector products of
these unit vectors are easy to evaluate and are equal to
The expression of the vector product of

and

in terms of the components of these vectors can be neatly rewritten as a
determinant:
Example: Problem 1.4 Use the vector product to find the
components of the unit vector

perpendicular to the plane shown in Figure 1.6.
Figure
1.6. Problem 1.4
The orientation of a plane in three dimensions is completely determined by
specifying two vectors that are parallel to this plane. Two possible vectors
are the vector

,
connecting (1,0,0) and (0,2,0), and the vector

,
connecting (1,0,0) and (0,0,3). The vector product of

and

is a vector

which is perpendicular to both

and

.
The vector

will therefore be pointing in the same direction as

.
Using the definition of the vector product in terms of the components of vectors

and

we can calculate

:
The length of the vector

is equal to
The unit vector

is parallel to

and has a length equal to 1. Therefore
1.3. Vector Transformation
The three coordinates required to specify the direction and length of a
vector are not uniquely defined. The same vector will in general have different
coordinates in different coordinate systems. There is an infinite number of
ways to choose a coordinate system, although the choice of the coordinate axes
is usually influenced by the symmetry of the problem. Figure 1.7 illustrates
two possible choices of the
y and
z axes. Coordinate system
S is defined by the
x,
y, and
z axes. Coordinate
system
S' is defined by the
x',
y', and
z' axes.
Coordinate system
S' can be obtained by rotating coordinate system
S through an angle

around the
x axis.
Figure
1.7. Coordinate Transformation.
The angle between the vector

and the
y axis is equal to

.
The
y and
z coordinates of the vector

are therefore equal to
The angle between the vector

and the
y' axis is equal to

.
The
y' and
z' coordinates of the vector

are therefore equal to
Using simple trigonometric relations we can obtain the component of the
vector

in
S' in terms of the components of the vector

in
S:
These relations show, not unexpectedly, that the coordinates of the vector

in coordinate system
S are related to the coordinates of vector

in coordinate system
S'. This relation can be rewritten in matrix
notation as
The rotation of the coordinate system S around the x axis
leaves the x axis unchanged. The x coordinate of any vector in
S will therefore be equal to the x' coordinate of this vector in
S'. The three dimensional version of this coordinate transformation is
therefore given by
The rotation of the coordinate system
S around the
x axis
will not change the length of the vector

.
This can be verified by calculating the length of the vector

in coordinate system
S. This length is equal to
The vector

in coordinate system
S' is specified by the following
coordinates:
The length of

in coordinate system
S' is defined in terms of its coordinates in
S' as
Using the relation between the coordinates of

in
S' and the coordinates of

in
S we obtain
In general the coordinate transformation describing a rotation around an
arbitrary axis can be written as
or, more compactly,
where A1 = Ax, A2 =
Ay, and A3 =
Az.
1.4. Differential calculus
The derivative df/dx of a function f(x) tells
us how rapidly the function varies when the argument x is changed by a
tiny amount dx:
A position dependent scalar function in three dimensions f(x,
y, z) will in general be a function of three variables. The
variation of f(x, y, z) between two closely spaced
points depends not only on the distance between the two points, but also on
their orientation. In this case, the theorem of partial derivatives can
be used to calculate the change in f(x, y,
z):
This equation can be rewritten as:
The first term in this expression is called the
gradient of
f, written as

,
and is defined as
The second term is called the
infinitesimal displacement vector

which is equal to
The change in the scalar function f can thus be written as
From the definition of the gradient

it is clear that it is a vector. The
direction of

points in the direction of maximum increase of the function
f. This
follows immediately from the expression of
df in terms of

:
The change in
f,
df, will be maximum when

= 0°. In this case,

and

are parallel. The
magnitude of

gives the slope of the function
f along the direction of maximum change.
If the gradient of
f vanishes at a certain point then the function
f has a maximum, a minimum, a saddle point, or a shoulder at that
point.
Example: Problem 1.13 Let

be the vector from some fixed point (
x0,
y0,
z0) to the point (
x,
y,
z), and let
r be its length. Show that
a)

b)

c)
What is he
general formula for

?
a) The vector

is equal to
The scalar function
r2 is equal to the square of the
length of the vector

and will be a function of
x,
y, and
z:
The gradient of r2 can be obtained using the following
partial derivatives:
The gradient of r2 is equal to
where

is the unit vector in the direction of

.
b)
The scalar function 1/
r can be written as
The gradient of 1/r can be obtained using the following partial
derivatives:
The gradient of 1/r is thus equal to
c) The most general form of

can be guessed by comparing the results obtained in part a) and b). These
answers suggest that the general form of

is given by
The operator

has the formal appearance of a vector and can be written as
The operator

is also called the
vector operator. If

behaves like a vector than we expect that

behaves like a vector, and consequently rotates like a
vector.
Example: Problem 1.14 Suppose that
f is a
function of two variables (
y and
z) only. Show that the gradient

transforms as a vector under a rotation about the
x
axis.
Consider a coordinate system
S with
x,
y, and
z axes. Coordinate system
S is related to a coordinate system
S' via a rotation by an angle

about the
x axis. The
y and
z coordinates of a vector in
S are related to the
y' and
z' coordinates in
S':
The gradient of f in the (y, z) frame is equal to
The gradient of f in the (y', z') frame is equal
to
If

transforms like a vector than the components of

must be related to the components of

in the following manner:
Using standard differential algebra we obtain the following relations for
∂f/∂y' and ∂f/∂z':
To evaluate ∂y/∂y',
∂y/∂z', ∂z/∂y', and
∂z/∂z' we must express y and z in terms
of y' and z'. This can be achieved by manipulating the following
two relations:
After a little algebra we obtain for y and z
Using these two relations we now can calculate the various partial
derivatives of y and z:
Using these partial derivatives we can now express

in terms of

:
These last two equations clearly show that

rotates like a vector.
If

mimics a vector it may be used in vector multiplication. Three different types
of vector multiplication can be carried out using

:
1. scalar
multiplication with scalar function
f,

,
also called the
gradient of
f.
2. dot (scalar) product with
vector function

,

,
also called the
divergence of

.
3. cross
(vector) product with vector function

,

,
also called the
curl of

.
We
will now discuss in detail the divergence and curl of

.
1.4.1. The Divergence
The divergence of a vector

is defined as
and is a scalar. The geometrical interpretation of the divergence of a
vector function

can be obtained by considering three special cases. First consider the vector
function defined as
This vector function

is shown at a few points in the
x-
y plane in Figure 1.8a. The
divergence of

is equal to
Now consider a vector function

defined such that the vectors point toward the origin (see Figure
1.8b):
For this vector function the divergence is equal to

Figure
1.8. Various vector functions discussed in the text.
Finally consider the case in which the vector function

is a constant vector of unit length along the
y axis (see Figure 1.8c),
independent of position:
For this vector function the divergence is equal to
We conclude that the divergence of a vector function

evaluated at a particular point
P is a measure of how much the vector
function

spreads
out:
1.

> 0: vector function

is spreading
out.
2.

< 0: vector function

is narrowing
in.
3.

= 0: vector function

is not spreading out or narrowing in.
One vector field that we will be
using frequently is the electric field

.
The three cases just discussed are relevant in electrostatics:
1. the
divergence of

generated by a positive point charge is positive (see Figure
1.9a).
2. the divergence of

generated by a negative point charge is negative (see Figure
1.9b).
3. the divergence of

generated by an infinitely large parallel-plate capacitor is zero (see Figure
1.9c).
Figure
1.9. a) Electric field generated by a positive point charge. b) Electric field
generated by a negative point charge. c) Electric field generated by an
infinitely large parallel-plate capacitor.
Example: Problem 1.17 In two dimensions, show that the
divergence transforms as a scalar under rotations.
Hint: use the method
of Problem 1.14 to calculate the derivatives.
If the divergence of a
vector function

transforms as a scalar under rotation then the divergence of

in a coordinate system
S(
x,
y,
z) must be identical
to the divergence of

in a coordinate system
S'(
x',
y',
z'). This
requires that
Consider a rotation about the x axis. In this case, x =
x' and vx = vx', and consequently
The
y' and
z' components of vector

are related to the
y and
z components of vector

in the following manner:
The partial derivative of vy' with respect to y'
is equal to
Using the expressions for ∂y/∂y' and
∂z/∂y' obtained in problem 1.14 we can rewrite this
expression as
In a similar manner the partial derivative of vz' with
respect to z' can be obtained:
Combining these last two equations, we obtain
and therefore, the divergence of a vector function

translates like a scalar under rotation.
1.4.2. The curl
The curl of vector function

is equal to
The curl of a vector function

evaluated at a certain point
P is a measure of how much the vector
function

curls around this point. For the vector functions

shown in Figure 1.8 the curl is zero.
Example: The magnetic
field Consider a wire of radius
R carrying a current
I
along the
z axis. The magnetic field produced by this wire depends only
on the distance to the center of the wire. Its strength is equal to
and its direction is tangent to the circle centered on the wire (see Figure
1.10).
Figure
1.10. Magnetic vector field.
The magnetic field is a vector field, and the vector function is equal
to
In order to calculate the curl of

at
r >
R we have to evaluate
∂
By/∂
x and
∂
Bx/∂
y at
r >
R:
The curl of

at
r >
R is thus equal to
In order to calculate the curl of

at
r <
R we have to evaluate
∂
By/∂
x and
∂
Bx/∂
y at
r <
R:
The curl of

at
r <
R is thus equal to
These calculations suggest that the curl of

evaluated at a certain point is proportional to the current density at that
point: outside the wire (
r >
R) the current density is 0, while
inside the wire (
r <
R) the current density is
I/(
πR2).
1.4.3. Product Rules
Without proof I mention the following product rules involving

that will be used frequently in this
course:
1.

:
f and
g are scalar
functions
2.

:

and

are vector
functions
3.

:
f is a scalar function and

is a vector
function
4.

:

and

are vector
functions
5.

:
f is a scalar function and

is a vector
function
6.

:

and

are vector functions
1.4.4. Second derivatives
By applying

twice we can construct five species of second derivatives:
1. The
divergence of a gradient:
where
T is a scalar function of
x,
y, and
z.

is also called the
Laplacian of
T.
2. The
curl of a
gradient (is always zero):
This can be shown easily:
3. The gradient of a divergence:
4. The divergence of a curl (is always zero):
This can be shown easily:
5. The
curl of a curl of a vector function

can be expressed in terms of the Laplacian and the gradient of the divergence of

:
1.5. Integral calculus
Consider a one-dimensional function df(x)/dx to be
integrated between x = a and x = b. This integral
can be obtained using the fundamental theorem of calculus, which states
that
In three dimensions this fundamental theorem is replaced by the
fundamental theorem for gradients, which states that
In a one-dimensional world there is just a single path from
a to
b. In a three-dimensional world there are an infinite number of ways to
move from
a to
b. Nevertheless, the
line integral of

depends only on the function values at the end points and not on the path
chosen. Thus we conclude:
1. The line integral

is independent of the path chosen between
a and
b.
2. The
line integral

around any closed loop is zero.
Example 1.6 Although the line
integral of the gradient of a vector function is independent of the path, the
same is not true for arbitrary vector functions. Calculate

from
a = (1, 1, 0) to
b = (2, 2, 0)
for the vector function

.
Do it first by path (1), shown in Figure 1.11, and then by path (2).
Figure
1.11. Example 1.6
Consider first path (1). The first part of this path is parallel to the
x axis (
y = 1) and along this segment the scalar product between

and

is equal to
The line integral along this segment is equal to
The second part of path (1) is parallel to the
y axis (
x = 2)
and along this segment the scalar product between

and

is equal to
The line integral along this segment is equal to
The line integral of

along path (1) is equal to the sum of the line integrals along the two segments,
and is thus equal to 1 + 10 = 11.
Consider now path (2). Along this path
the
x and
y coordinates are equal (
x =
y and
dx =
dy). The displacement vector

along this path is equal to
The scalar product between

and

along path (2) is equal to
The line integral along segment (2) is equal to
Comparing the line integral of

for path (1) and the line integral of

for path (2), we conclude that for the vector function

the line integral is path dependent.
In a three-dimensional world we can
have, besides line integrals, surface integrals and volume integrals. Various
fundamental theorems can simplify the calculation of surface and volume
integrals. For example, the
fundamental theorem of divergences states
that the integral of a divergence of a vector function

over a volume is equal to the surface integral of the vector function

over the surface that bounds the volume:

The right-hand side of this equation is called the
flux of

through the surface that bounds the volume. The vector

is a vector whose magnitude is equal to the area of an infinitesimal surface
element and whose direction is perpendicular to the surface, pointing outwards.
The left-hand side of this equation represents a
source term (for
example,

integrated over a volume is proportional to the total charge present in that
volume).
Example: Problem 1.32 Test the divergence theorem
for the function

.
Take as your volume the cube shown in Figure 1.12, with sides of length
2.
The divergence of

is equal to
The volume integral of

over the volume of the cube is equal to

Figure
1.12. Problem 1.32.
To evaluate the surface integral of

across the surface of the cube, we have to consider each of the six faces
separately. Start with considering the face of the cube in the
x-
y plane (with
z = 0). The vector

of this surface is equal to

.
The scalar product of

and

is equal to
In the last step we have used the fact that
z = 0 on this face. For
the same reason the scalar product of

and

is equal to zero on the face in the
x-
z plane and on the face in
the
y-
z plane. On the face of the cube parallel to the
x-
y plane and with
z = 2 the scalar product of

and

is equal to
The surface integral of

over this face is equal to
On the face the cube parallel to the
x-
z plane and with
y = 2 the scalar product of

and

is equal to
The surface integral of

over this face is equal to
On the face of the cube parallel to the
y-
z plane and with
x = 2 the scalar product of

and

is equal to
The surface integral of

over this face is equal to
The surface integral across the surface of the cube is equal to the sum of
the surface integrals across each of the 6 faces of the cube, and is thus equal
to 24 + 16 + 8 = 48. This is equal to the volume integral of

over the volume of the cube.
The
fundamental theorem for curls,
also known as Stokes' theorem, states that
Here,

is an infinitesimal surface element. It is a vector whose magnitude is equal to
the area of the surface element and whose direction is perpendicular to the
surface. The vector

is tangential to the boundary of the surface. The orientation of the surface
vector

and the direction of integration of the boundary should be consistent with he
right-hand rule. The following two corollaries follow from Stokes'
theorem:
1.

depends only on the boundary line, not on the particular surface
used.
2.

around any closed surface.
Example: Problem 1.33 Test
Stokes' theorem for the vector function

,
using the triangular shaded area shown in Figure 1.13.
The curl of

is equal to

Fig.
1.13. Problem 1.33.
The surface vector

is perpendicular to the surface and is equal to

.
The scalar product between

and

is equal to
The surface integral of

is equal to
To evaluate the line integral around the boundary of the surface we have to
evaluate the line integral along each of the three sides of the triangle. The
direction of evaluation of the integral must be consistent with the chosen
direction of

.
The right-hand rule requires that the line integral is evaluated counter
clockwise. First consider the segment between (0, 0, 0) and (0, 2, 0). Along
this segment

.
The vector product between

and

is equal to
since
z = 0 along this segment. The line integral along this
segment is therefore equal to zero. Now consider the line segment between (0,
2, 0) and (0, 0, 2). Along this segment

and

.
The vector product between

and

is equal to
since x = 0 along this segment. Along this segment z = 2 -
y and thus
The line integral of

along this segment is thus equal to
Note: the limits of the integral are chosen such that the line integral is
evaluated in a counter clockwise direction. The third segment of the boundary
to be considered connects (0, 0, 2) and (0, 0, 0). The vector

along this segment is equal to

.
The vector product between

and

is equal to
since
x = 0 along this segment. The line integral of

is therefore equal to zero. The line integral along the boundary of the surface
is equal to the sum of the line integral of along each of the three segments.
Thus
which is equal to the surface integral of

.
1.6. Curvilinear coordinates
The Cartesian coordinate system is a coordinate system that is often used
in calculations involving systems with no apparent symmetry. To describe
systems that have spherical or cylindrical symmetry it is often more convenient
to use spherical coordinates or cylindrical coordinates,
respectively. These two coordinate systems will be discussed in this
section.
1.6.1. Spherical coordinates
Spherical coordinates are always used when the system under consideration
has spherical symmetry. The location of a point
P (see Figure 1.14) is
completely determined by specifying the following three
coordinates:
1.
r : the distance between the origin and
P.
2.
: the
polar angle which is the angle between the vector
P and
the
z-axis (see Figure
1.14).
3.
: the
azimuthal angle which is the angle between the projection of
the vector to
P in the
x-
y plane and the
x
axis.
Figure
1.14. Spherical coordinates.
In general, any vector

can be expressed in terms of these three coordinates:
In contrast to the unit vectors

,

,
and

in a Cartesian coordinate system, the unit vectors

,

,
and

in a spherical coordinate system are not constant; they change direction as
P moves around.
Consider a point
P (see Figure 1.14) which is
defined by the three spherical coordinates (
r,

,
and
).
The corresponding Cartesian coordinates are
The unit vectors in the spherical coordinate system can be calculated as
follows:
1.

: The
direction of this unit vector points from the origin to point
P:
2.

: A
change
d
in the polar angle

will cause a change in the direction of the vector

,
without changing its length:
The unit vector

is defined as
3.

: A
change
d
in the polar angle

will cause a change in the direction of the vector

without changing its length:
The unit vector

is defined as
The expressions for

,

,
and

clearly show that the direction of these unit vectors depends on the point being
described. In contrast, in a Cartesian coordinate system the direction of the
three unit vectors is fixed, and independent of the point being
described.
An infinitesimal element of length in the

direction is simply
dr:
The length of an infinitesimal element of length in the

direction (as a result of a change in the polar angle of
d
)
is equal to
The length of an infinitesimal element of length in the

direction (as a result of a change in the azimuthal angle
d
)
is equal to
The most general infinitesimal displacement is thus equal to
and this expression is used to evaluate line integrals in spherical
coordinates. The infinitesimal volume element

,
in spherical coordinates, is the product of

,

,
and

:
There is no general expression for an infinitesimal surface element

since it depends entirely on the orientation of the surface. For example,
consider the surface of a sphere of radius
r. On this surface
r
is constant, and the orientation of the surface is perpendicular to

.
In this case

is equal to
Most calculations in spherical coordinates involve the application of
vector derivatives. I will therefore summarize here, without derivation, the
vector derivatives in spherical coordinates:
1. the gradient of a
scalar function T:
2. the
divergence of a vector function

:
3. the
curl of a vector function

:
4. the Laplacian of a scalar function T:
1.6.2. Cylindrical coordinates
Cylindrical coordinates are always used when the system under
consideration has cylindrical symmetry. The
z axis is defined such that
it coincides with the center axis of the cylinder. The location of a point
P (see Figure 1.15) is completely determined by specifying the following
three coordinates:
1.
r : the distance between
P and the
z axis (see Figure
1.15).
2.
: the azimuthal angle which is the angle between the projection of the
vector to
P in the
x-
y plane and the
x axis (see
Figure 1.15).
3.
z : the
z coordinate of point
P
(see Figure 1.15).
In general, any vector

can be expressed in terms of these three coordinates:
However, the unit vectors

and

are not constant; they change direction as
P moves around.
Figure
1.15. Cylindrical coordinates.
Consider a point
P (see Figure 1.15) which is defined by three
spherical coordinates (
r,

,
and
z). The corresponding Cartesian coordinates are
An infinitesimal change in
r,

,
or
z will produce the following infinitesimal displacements:
The infinitesimal displacement vector, to be used when evaluating line
integrals in spherical coordinates, is equal to
The infinitesimal volume element dτ, to be used in volume
integration, is equal to
The infinitesimal area element

,
to be used in surface integration, is not uniquely determined and depends on the
orientation of the surface. For example, if we carry out a surface integration
over the surface of a cylinder of radius
r (fixed), the infinitesimal
surface element to be used is equal to
Most calculations in cylindrical coordinates involve the application of
vector derivatives. I will therefore summarize here, without derivation, the
vector derivatives in cylindrical coordinates:
1. The gradient of
a scalar function T:
2. The
divergence of a vector function

:
3. The
curl of a vector function

:
4. The Laplacian of a scalar function T:
Example: Problem 1.42
a) Find the divergence of the
function
b) Test the divergence theorem for this function, using the
quarter-cylinder (radius 2, height 5) shown in Figure 1.16.
c) Find the curl
of

.
Figure
1.16. Problem 1.42.
a) To determine the divergence of

we have to calculate the following partial derivatives:
Using the definition of the divergence of

in terms of cylindrical coordinates we find that
b) The volume integral of the divergence of

is equal to
To calculate the surface integral of

we have to consider 5 different surfaces:
1. The surface parallel to
the x-y plane at z = 0. The infinitesimal surface
element for this surface is equal to
(Note: the direction of

is perpendicular to the surface and pointing outwards). The scalar product
between

and

is equal to
since
z = 0 on this surface. The surface integral of

across this surface is therefore equal to zero.
2. The surface
parallel to the x-y plane at z = 5. The infinitesimal
surface element for this surface is equal to
The scalar product between

and

is equal to
since
z = 5 on this surface. The surface integral of

across this surface is equal to
3. The surface in the x-z plane with
= 0. The infinitesimal surface element for this surface is equal
to
The scalar product between

and

is equal to
since

= 0 on this surface. The surface integral of

across this surface is therefore equal to zero.
4. The surface in the
y-z plane with
= π/2. The infinitesimal surface element for this surface is equal
to
The scalar product between

and

is equal to
since

= π/2 on this surface. The surface integral of

across this surface is therefore equal to zero.
5. The cylinder wall
(r = 2). The infinitesimal surface element for this surface is equal
to
The scalar product between

and

is equal to
The surface integral of

across this surface is equal to
The surface integral of

across the whole surface is equal to the sum of the surface integral of

across each of these five surfaces:
which is equal to the volume integral of

.
1.7. The Dirac delta function
Consider the vector function

:
Consider the surface integral of

across the surface of a sphere of radius
R. The surface element vector

for this surface is equal to
The scalar product of

and

is equal to
The surface integral of

across the surface of this sphere is equal to
and is independent of R. Applying the divergence theorem we
conclude that
for every sphere, centered on
r = 0, and independent of the radius
R. This suggests that the only contribution to the volume integral of

comes from a single point at
r = 0. The divergence of

is zero at every point except at
r = 0:
At
r = 0, the divergence of

is undefined since 0/0 is undefined.
The function

is called the
Dirac delta function

and will be used frequently in this course. The Dirac delta function has the
following
properties:

for
any volume that includes
r =
0

for
any volume that does not include
r = 0
The Dirac delta function

is thus defined as
In Problem 1.13 we showed that the gradient of 1/
r is equal to

.
This relation can be used to rewrite the expression for the Dirac delta function
as
In this section we will discuss the use of the one- and three-dimensional
Dirac delta functions.
1.7.1. The one-dimensional Dirac delta function
The one-dimensional Dirac delta function

is defined as
and its integral is equal to
Any integral of

between
x =
a and
x =
b will be equal to 1 if
a < 0 <
b.
Every time we will be using the Dirac delta
function it will be used as part of the integrand of an integral. For
example:
Example: Problem 1.43 Evaluate the following
integrals:
a)

b)

c)

d)

a) The
Dirac delta function

will be 0 for all
x except
x = 3. The limits of the integral
include
x = 3. Thus
b) The Dirac delta function

will be 0 for all
x except
x = π. The limits of the integral
include
x = π. Thus
c) The Dirac delta function

will be 0 for all
x except
x = - 1. The limits of the integral do
not include
x = - 1. Thus
d) The Dirac delta function

will be 0 for all
x except
x = - 2. The limits of the integral
include
x = - 2. Thus
Note: Most mistakes made in the evaluation of integrals
containing the Dirac delta function occur when the argument of the Dirac delta
function has a form other than (x + a). This is best illustrated
by looking at the following example:
1.7.2. The three-dimensional Dirac delta function
The three-dimensional Dirac delta function

is zero everywhere except at

= 0 (
x = 0,
y = 0,
z = 0). Every time we will be using the
three-dimensional Dirac delta function it will be used as part of the integrand
of an integral. For example:
if the integration volume
V includes

= 0.
Example: Problem 1.47a Evaluate the following
integral:
The Dirac delta function is zero everywhere except at

.
The volume integral is thus equal to
1.8. Helmholtz Theorem
Consider being told that the divergence of a vector function

is the scalar function
D and that the curl of the vector function

is the vector function

.
Is this sufficient information to determine the vector function

uniquely?
The vector function

satisfies the following relations:
The divergence of

must be zero since the divergence of the curl of any vector function is zero
(

).
Consider the following solution:
where
If

is a correct solution than the divergence of

must be equal to
D:
Here we have used one of the properties of second derivatives, which states
that

for any vector function

.
If

is a correct solution than the curl of

must be equal to

:
Here we have used one of the properties of second derivatives which states
that

for any scalar function
U. The Laplacian of

is equal to
The divergence of

is equal to
In this derivation we have used the product rule for divergences (Griffiths
page 21), the divergence theorem and the fact that the divergence of

is equal to zero. The surface integral of

/(
r
-
r') will be equal to zero if

goes to zero faster than 1/
r2 when
r → ∞.
If this is the case, the divergence of

is equal to zero, and consequently
So far we have verified that we can obtain a vector function

if its divergence and curl are given. However, the vector function

found is not unique. Consider a vector function

whose curl and divergence vanish. Since the curl is distributive we can express
the curl of the new vector function in terms of the curl of

and the curl of

:
Since the divergence is distributive we can express the divergence of the
new vector function in terms of the divergence of

and the divergence of

:
However, there is no function that has zero divergence and zero curl and
goes to zero at infinity. Since it is expected that the electric and magnetic
fields go to zero at large distances we will exclude this type of contributions
to

by requiring that

goes to zero at large distances. With this requirement,

is uniquely defined if its curl and divergence are given. This conclusion is
know as the
Helmholtz theorem:
If the divergence

and the curl

of a vector function

are specified, and if they both go to zero faster than 1/
r2 as
r → ∞ and if

goes to zero as
r → ∞ then

is given uniquely by
where
1.9. Scalar and vector potentials
Consider a vector function

which is the gradient of a scalar function

:

.
The scalar function

is also called the
scalar potential of the field

.
The vector field

,
defined by the scalar function

,
has the following
properties:
1.

for any closed
loop.
2.

is independent of path, for any given set of end
points.
3.

everywhere.
The vector field

generated by the scalar function

is called a
curl-less field. In Example 1.6 we showed that for the
vector function

the line integral is path dependent. The curl of this vector function is equal
to

and consequently this vector function can not be written as the gradient of a
scalar function
U. Therefore, the line integral of this vector function
is path dependent, as was demonstrated in Example 1.6.
Now consider a vector
field

which is the curl of a vector function

:

.
The vector function

is called the
vector potential for the field

.
The vector field

defined by the vector potential
W has the following
properties:
1.

for any closed
surface.
2.

is independent of the surface chosen, and depends only on the boundary
line.
3.

everywhere.
The vector field

generated by the vector potential

is called a
divergence-less field.
The scalar potential
U and
the vector potential

are not unique. Any constant can be added to
U without effecting its
gradient, since the gradient of a constant is equal to zero. Any gradient of a
scalar function can be added to

without effecting its curl, since the curl of a gradient is equal to
zero.