Course PDF

1
2A1 Vector Algebra and Calculus
8 Lectures MT 2013 Stephen Roberts

Rev Oct 2013 sjrob@robots.ox.ac.uk
Course Page www.robots.ox.ac.uk/∼sjrob/Teaching/Vectors
Overview
Many of you will know a good deal already about Vector Algebra — how to add and subtract
vectors, how to take scalar and vector products of vectors, and something of how to describe
geometric and physical entities using vectors. This course will remind you about that good stuff,
but goes on to introduce you to the subject of Vector Calculus which, like it says on the can,
combines vector algebra with calculus.
To give you a feeling for the issues, suppose you were interested in the temperature T of water in
a river. Temperature T is a scalar, and will certainly be a function of a position vector x = (x , y , z )
and may also be a function of time t: T = T (x, t). It is a scalar field.
Suppose now that you kept y , z , t constant, and asked what is the change in temperature as you
move a small amount in x ? No doubt you’d be interested in calculating ∂T /∂x. Similarly if you
kept the point fixed, and asked how does the temperature change of time, you would be interested
in ∂T /∂t.
But why restrict ourselves to movements up-down, left-right, etc? Suppose you wanted to know
what the change in temperature along an arbitrary direction. You would be interested in
∂T
,
∂x
but how would you calculate that? Is ∂T /∂x a vector or a scalar?

Now let’s dive into the flow. At each point x in the stream, at each time t, there will be a stream
velocity v(x, t). The local stream velocity can be viewed directly using modern techniques such
as laser Doppler anemometry, or traditional techniques such as throwing twigs in. The point now
is that v is a function that has the same four input variables as temperature did, but its output
result is a vector. We may be interested in places x where the stream suddenly accelerates, or
vortices where the stream curls around dangerously. That is, we will be interested in finding the
acceleration of the stream, the gradient of its velocity. We may be interested in the magnitude of
the acceleration (a scalar). Equally, we may be interested in the acceleration as a vector, so that
we can apply Newton’s law and figure out the force.
This is the stuff of vector calculus.
2
Grey book
Vector algebra: scalar and vector products; scalar and vector triple products; geometric appli-
cations. Differentiation of a vector function; scalar and vector fields. Gradient, divergence and
curl - definitions and physical interpretations; product formulae; curvilinear coordinates. Gauss’
and Stokes’ theorems and evaluation of integrals over lines, surfaces and volumes. Derivation
of continuity equations and Laplace’s equation in Cartesian, cylindrical and spherical coordinate
systems.
Course Content
• Introduction and revision of elementary concepts, scalar product, vector product.
• Triple products, multiple products, applications to geometry.
• Differentiation and integration of vector functions of a single variable.
• Curvilinear coordinate systems. Line, surface and volume integrals.
• Vector operators.
• Vector Identities.
• Gauss’ and Stokes’ Theorems.
• Engineering Applications.
Learning Outcomes
You should be comfortable with expressing systems (especially those in 2 and 3 dimensions) using
vector quantities and manipulating these vectors without necessarily going back to some underlying
coordinates.
You should have a sound grasp of the concept of a vector field, and be able to link this idea to
descriptions of various physical phenomena.
You should have a good intuition of the physical meaning of the various vector calculus operators
and the important related theorems. You should be able to interpret the formulae describing
physical systems in terms of this intuition.
References
Although these notes cover the material you need to know you should, wider reading is essen-
tial. Different explanations and different diagrams in books will give you the perspective to glue
everything together, and further worked examples give you the confidence to tackle the tute sheets.
• J Heading, ”Mathematical Methods in Science and Engineering”, 2nd ed., Ch.13, (Arnold).
• G Stephenson, ”Mathematical Methods for Science Students”, 2nd ed., Ch.19, (Longman).
• E Kreyszig, ”Advanced Engineering Mathematics”, 6th ed., Ch.6, (Wiley).
• K F Riley, M. P. Hobson and S. J. Bence, ”Mathematical Methods for the Physics and
Engineering” Chs.6, 8 and 9, (CUP).
• A J M Spencer, et. al. ”Engineering Mathematics”, Vol.1, Ch.6, (Van Nostrand Reinhold).
• H M Schey, “Div, Grad, Curl and all that”, Norton
3
Course WWW Pages

Pdf copies of these notes (including larger print versions), tutorial sheets, FAQs etc will be accessible
from
www.robots.ox.ac.uk/∼sjrob/Teaching/Vectors
4
Lecture 1
Vector Algebra
1.1 Vectors
Many physical quantities, such a mass, time, temperature are fully specified by one
number or magnitude. They are scalars. But other quantities require more than
one number to describe them. They are vectors. You have already met vectors in
their more pure mathematical sense in your course on linear algebra (matrices and
vectors), but often in the physical world, these numbers specify a magnitude and
a direction — a total of two numbers in a 2D planar world, and three numbers in
3D.
Obvious examples are velocity, acceleration, electric field, and force. Below, prob-
ably all our examples will be of these “magnitude and direction” vectors, but we
should not forget that many of the results extend to the wider realm of vectors.
There are three slightly different types of vectors:
• Free vectors: In many situtations only the magnitude and direction of a

vector are important, and we can translate them at will (with 3 degrees of
freedom for a vector in 3-dimensions).
• Sliding vectors: In mechanics the line of action of a force is often important
for deriving moments. The force vector can slide with 1 degree of freedom.
• Bound or position vectors: When describing lines, curves etc in space, it is
obviously important that the origin and head of the vector are not translated
about arbitrarily. The origins of position vectors all coincide at an overall
origin O.
One the advantages of using vectors is that it frees much of the analysis from
the restriction of arbitrarily imposed coordinate frames. For example, if two free
vectors are equal we need only say that their magnitudes and directions are equal,
and that can be done with a drawing that is independent of any coordinate system.
However, coordinate systems are ultimately useful, so it useful to introduce the
idea of vector components. Try to spot things in the notes that are independent
5
6 LECTURE 1. VECTOR ALGEBRA
r r
r1 2 3
O
Free vectors Sliding vectors Position vectors
Figure 1.1:
of coordinate system.
1.1.1 Vector elements or components in a coordinate frame

A method of representing a vector is to list the values of its elements or components
in a sufficient number of different (preferably mutually perpendicular) directions,
depending on the dimension of the vector. These specified directions define a
coordinate frame. In this course we will mostly restrict our attention to the
3-dimensional Cartesian coordinate frame O(x , y , z ). When we come to examine
vector fields later in the course you will use curvilinear coordinate frames, especially
3D spherical and cylindrical polars, and 2D plane polar, coordinate systems.
x
2
a3
a
x j
1
a2
a1
i
Figure 1.2: Vector components.
In a Cartesian coordinate frame we write

a = [a1 , a2, a3 ] = [x2 − x1 , y2 − y1 , z2 − z1 ] or a = [ax , ay , az ]
as sketched in Figure 1.2. Defining ı̂ı ,̂ , k̂k as unit vectors in the x , y , z directions
ı̂ı = [1, 0, 0] ̂ = [0, 1, 0] k̂k = [0, 0, 1]
1.1. VECTORS 7
a+b
a a−b c
k
c b+c
b
b
a+b
i j
a
(a) (b)
Figure 1.3: (a) Addition of two vectors is commutative, but subtraction isn’t. Note that the
coordinate frame is irrelevant. (b) Addition of three vectors is associative.
we could also write

a = a1ı̂ı + a2̂ + a3k̂k .
Although we will be most often dealing with vectors in 3-space, you should not
think that general vectors are limited to three components.
In these notes we will use bold font to represent vectors a, ω , In your written work,
underline the vector symbol a, ω and be meticulous about doing so. We shall use
the hat to denote a unit vector.
1.1.2 Vector equality

Two free vectors are said to be equal iff their lengths and directions are the same.
If we use a coordinate frame, we might say that corresponding components of
the two vectors must be equal. This definition of equality will also do for position
vectors, but for sliding vectors we must add that the line of action must be identical
too.
1.1.3 Vector magnitude and unit vectors

Provided we use an orthogonal coordinate system, the magnitude of a 3-vector is
q
a = |a| = a12 + a22 + a32
To find the unit vector in the direction of a, simply divide by its magnitude
a
â = .
|a|
1.1.4 Vector Addition and subtraction

Vectors are added/subtracted by adding/subtracting corresponding components,
exactly as for matrices. Thus
a + b = [a1 + b1, a2 + b2 , a3 + b3 ]
Addition follows the parallelogram construction of Figure 1.3(a). Subtraction (a −

b) is defined as the addition (a + (−b)). It is useful to remember that the vector
a − b goes from b to a.
The following results follow immediately from the above definition of vector addi-
tion:
(a) a + b = b + a (commutativity) (Figure 1.3(a))
(b) (a + b) + c = a + (b + c) = a + b + c (associativity) (Figure 1.3(b))
(c) a + 0 = 0 + a = a, where the zero vector is 0 = [0, 0, 0].
(d) a + (-a) = 0
1.1.5 Multiplication of a vector by a scalar. (NOT the scalar product!)

Just as for matrices, multiplication of a vector a by a scalar c is defined as multi-
plication of each component by c, so that
ca = [ca1, ca2, ca3].
It follows that:
p
|ca| = (ca1)2 + (ca2)2 + (ca3)2 = |c||a|.
The direction of the vector will reverse if c is negative, but otherwise is unaffected.
(By the way, a vector where the sign is uncertain is called a director.)
♣ Example
Q. Coulomb’s law states that the electrostatic force on charged particle Q due
to another charged particle q1 is
Qq1
F = K 2 êr
r
where r is the vector from q1 to Q and r̂ is the unit vector in that same
direction. (Note that the rule “unlike charges attract, like charges repel” is
built into this formula.) The force between two particles is not modified by
the presence of other charged particles.
Hence write down an expression for the force on Q at R due to N charges qi
at ri .
A. The vector from qi to Q is R − ri . The unit vector in that direction is
(R − ri )/|R − ri |, so the resultant force is
N
X Qqi
F(R) = K (R − ri ) .
|R − ri |3
i =1
Note that F(R) is a vector field.

1.2. SCALAR, DOT, OR INNER PRODUCT 9
1.2 Scalar, dot, or inner product

This is a product of two vectors results in a scalar quantity and is defined as follows
for 3-component vectors:
a · b = a1 b1 + a2 b2 + a3 b3 .
Note that
a · a = a12 + a22 + a32 = |a|2 = a2.
The following laws of multiplication follow immediately from the definition:
(a) a · b = b · a (commutativity)
(b) a · (b + c) = a · b + a · c (distributivity with respect to vector addition)
(c) (λa) · b = λ(a · b) = a · (λb) scalar multiple of a scalar product of two vectors
1.2.1 Geometrical interpretation of scalar product
B
θ
b a
a−b
θ
Projection of b onto
O A
a direction of a
(a) (b)
Figure 1.4: (a) Cosine rule. (b) Projection of b onto a.
Consider the square magnitude of the vector a − b. By the rules of the scalar
product, this is
|a − b|2 = (a − b) · (a − b)
= a · a + b · b − 2(a · b)
= a2 + b2 − 2(a · b)
But, by the cosine rule for the triangle OAB (Figure 1.4a), the length AB 2 is given
by
|a − b|2 = a2 + b2 − 2ab cos θ
where θ is the angle between the two vectors. It follows that
a · b = ab cos θ,
which is independent of the co-ordinate system used, and that |a · b| ≤ ab. Con-
versely, the cosine of the angle between vectors a and b is given by cos θ = a·b/ab.
1.2.2 Projection of one vector onto the other

Another way of describing the scalar product is as the product of the magnitude
of one vector and the component of the other in the direction of the first, since
b cos θ is the component of b in the direction of a and vice versa (Figure 1.4b).
Projection is particularly useful when the second vector is a unit vector — a · ı̂ı is
the component of a in the direction of ı̂ı .
Notice that if we wanted the vector component of b in the direction of a we
would write
(b · a)a
(b · â)â = .
a2
In the particular case a · b = 0, the angle between the two vectors is a right angle
and the vectors are said to be mutually orthogonal or perpendicular — neither
vector has any component in the direction of the other.
An orthonormal coordinate system is characterised by ı̂ı · ı̂ı = ̂ · ̂ = k̂k · k̂k = 1; and
ı̂ı · ̂ = ̂ · k̂k = k̂k · ı̂ı = 0.
1.2.3 A scalar product is an “inner product”

So far we have been writing our vectors as row vectors a = [a1 , a2, a3 ]. This is
convenient because it takes up less room than writing column vectors
 
a1
a =  a2  .
a3
In matrix algebra vectors are more usually defined as column vectors, as in
    
M11 M12 M13 a1 v1
 M21 M22 M23   a2  =  v2 
M31 M32 M33 a3 v3
1.2. SCALAR, DOT, OR INNER PRODUCT 11
and a row vector is written as a⊤ . Now for most of our work we can be quite
relaxed about this minor difference, but here let us be fussy.
Why? Simply to point out at that the scalar product is also the inner product
more commonly used in linear algebra. Defined as a⊤ b when vectors are column
vectors as
 
b1
a · b = a⊤ b = [a1 , a2, a3]  b2  = a1b1 + a2b2 + a3b3 .
b3
Here we treat a n-dimensional column vector as an n × 1 matrix.
(Remember that if you multiply two matrices Mm×n Nn×p then M must have the
same columns as N has rows (here denoted by n) and the result has size (rows ×
columns) of m × p. So for n-dimensional column vectors a and b, a⊤ is a 1 × n
matrix and b is n × 1 matrix, so the product a⊤ b is a 1 × 1 matrix, which is (at
last!) a scalar.)
♣ Examples
Q1. A force F is applied to an object as it moves by a small amount δr. What
work is done on the object by the force?
A1. The work done is equal to the component of force in the direction of the dis-
placement multiplied by the displacement itself. This is just a scalar product:
δW = F · δr .
Q2. A cube has four diagonals, connecting opposite vertices. What is the angle
between an adjacent pair?
A2. Well, you could plod through using

Pythagoras’ theorem to find the length
of the diagonal from cube vertex to cube [−1,1,1] [1,1,1]
centre, and perhaps you should to check
the following answer.
k
The directions of the diagonals are
[±1, ±1, ±1]. The ones shown in the
figure are [1, 1, 1] and [−1, 1, 1]. The
j
angle is thus
i
[1, 1, 1] · [−1, 1, 1] −1 1
θ = cos−1 √ 2 √ = cos
1 + 12 + 12 −12 + 12 + 12 3
Q3. A pinball moving in a plane with velocity s bounces (in a purely elastic impact)
from a baffle whose endpoints are p and q. What is the velocity vector after
the bounce?
A3. Best to refer everything to a coordi-

nate frame with principal directions
û along and v̂ perpendicular to the
baffle: q
q−p
û =
|q − p| ^u
v̂ = u⊥ = [−uy , ux ]
p
Thus the velocity before impact is
sbefore = (s.û)û + (s.v̂)v̂ ^v
After the impact, the component of s

velocity in the direction of the baf-
fle is unchanged and the component
normal to the baffle is negated:
safter = (s.û)û − (s.v̂)v̂
1.2.4 Direction cosines use projection

Direction cosines are commonly used in the field of crystallography. The quantities
a · ı̂ı a · ̂ a · k̂k

λ= , µ= , ν=
a a a
represent the cosines of the angles which the vector a makes with the co-ordinate
vectors ı̂ı ,̂ , k̂k and are known as the direction cosines of the vector a. Since
a ·ı̂ı = a1 etc, it follows immediately that a = a(λı̂ı + µ̂ + ν k̂k ) and λ2 + µ2 + ν 2 =
1 2 2 2
a2 [a1 + a2 + a3 ] = 1
1.3 Vector or cross product

The vector product of two vectors a and b is denoted by a × b and is defined as
follows
a × b = (a2b3 − a3b2 )ı̂ı + (a3b1 − a1b3 )̂ + (a1b2 − a2 b1)k̂k .

1.3. VECTOR OR CROSS PRODUCT 13
j
i
Figure 1.5: The direction cosines are cosines of the angles shown.
It is MUCH more easily remembered in terms of the (pseudo-)determinant

ı̂ı ̂ k̂k

a × b = a1 a2 a3
b1 b2 b3
where the top row consists of the vectors ı̂ı , ̂ , k̂k rather than scalars.
Since a determinant with two equal rows has value zero, it follows that a × a = 0.
It is also easily verified that (a × b) · a = (a × b) · b = 0, so that a × b is orthogonal
(perpendicular) to both a and b, as shown in Figure 1.6.
Note that ı̂ı × ̂ = k̂k , ̂ × k̂k = ı̂ı , and k̂k × ı̂ı = ̂ .
The magnitude of the vector product can be obtained by showing that
|a × b|2 + (a · b)2 = a2 b2
from which it follows that
|a × b| = ab sin θ ,
which is again independent of the co-ordinate system used. This is left as an

exercise.
Unlike the scalar product, the vector product does not satisfy commutativity but
is in fact anti-commutative, in that a × b = −b × a. Moreover the vector product
does not satisfy the associative law of multiplication either since, as we shall see
later a × (b × c) 6= (a × b) × c.
Since the vector product is known to be orthogonal to both the vectors which form
the product, it merely remains to specify its sense with respect to these vectors.
Assuming that the co-ordinate vectors form a right-handed set in the order ı̂ı ,̂ , k̂k
it can be seen that the sense of the the vector product is also right handed, i.e
the vector product has the same sense as the co-ordinate system used.

ı̂ı ̂ k̂k

ı̂ı × ̂ = 1 0 0 = k̂k .
0 1 0
In practice, figure out the direction from a right-handed screw twisted from the
first to second vector as shown in Figure 1.6(a).
ax b b
axb
in right−hand screw sense
b bsin θ
θ
a
Plane of vectors a and b a
Figure 1.6: (a)The vector product is orthogonal to both a and b. Twist from first to second and
move in the direction of a right-handed screw. (b) Area of parallelogram is ab sin θ.
1.3.1 Geometrical interpretation of vector product

The magnitude of the vector product (a × b) is equal to the area of the parallelo-
gram whose sides are parallel to, and have lengths equal to the magnitudes of, the
vectors a and b (Figure 1.6b). Its direction is perpendicular to the parallelogram.
♣ Example
Q. g is vector from A [1,2,3] to B [3,4,5].
ℓ̂ℓ is the unit vector in the direction from O to A.
Find m̂, a UNIT vector along g × ℓ̂ℓ
Verify that m̂ is perpendicular to ℓ̂ℓ.
Find n̂, the third member of a right-handed coordinate set ℓ̂ℓ, m̂, n̂.
A.
g = [3, 4, 5] − [1, 2, 3] = [2, 2, 2]
1
ℓ̂ℓ = √ [1, 2, 3]
14

ı̂ı ̂ k̂k
1 = √1 [2, −4, 2]

ℓ
g × ℓ̂ = √ 2 2 2
14 14
1 2 3
1.3. VECTOR OR CROSS PRODUCT 15
Hence
1
m̂ = · √ [2, −4, 2]
24
and
n̂ = ℓ̂ℓ × m̂
Revised Oct 2013

Lecture 2
Multiple Products. Geometry using Vectors
2.1 Triple and multiple products

Using mixtures of the pairwise scalar product and vector product, it is possible to
derive “triple products” between three vectors, and indeed n-products between n
vectors.
There is nothing about these that you cannot work out from the definitions of pair-
wise scalar and vector products already given, but some have interesting geometric
interpretations, so it is worth looking at these.
2.1.1 Scalar triple product

This is the scalar product of a vector product and a third vector, i.e. a · (b × c).
Using the pseudo-determinant expression for the vector product, we see that the
scalar triple product can be represented as the true determinant

a1 a2 a3

a · (b × c) = b1 b2 b3
c1 c2 c3
You will recall that if you swap a pair of rows of a determinant, its sign changes;
hence if you swap two pairs, its sign stays the same.

a1 a2 a3 c1 c2 c3 c1 c2 c3

b1 b2 b3 1st SWAP b1 b2 b3 2nd SWAP a1 a2 a3

c1 c2 c3 a1 a2 a3 b1 b2 b3
+ − +
This says that
(i) a · (b × c) = b · (c × a) = c · (a × b) (Called cyclic permutation.)
17
18 LECTURE 2. MULTIPLE PRODUCTS. GEOMETRY USING VECTORS
(ii) a · (b × c) = −b · (a × c) and so on. (Called anti-cyclic permutation.)

(iii) The fact that a · (b × c) = (a × b) · c allows the scalar triple product to be
written as [a, b, c]. This notation is not very helpful, and we will try to avoid
it below.
2.1.2 Geometrical interpretation of scalar triple product

The scalar triple product gives the volume of the parallelopiped whose sides are
represented by the vectors a, b, and c.
We saw earlier that the vector product (a × b) has magnitude equal to the area
of the base, and direction perpendicular to the base. The component of c in this
direction is equal to the height of the parallelopiped shown in Figure 2.1(a).
n
c cos β b
c
c
β
b
a
a
Figure 2.1: (a) Scalar triple product equals volume of parallelopiped. (b) Coplanarity yields zero
scalar triple product.
2.1.3 Linearly dependent vectors

If the scalar triple product of three vectors is zero
a · (b × c) = 0
then the vectors are linearly dependent. That is, one can be expressed as a linear
combination of the others. For example,
a = λb + µc
where λ and µ are scalar coefficients.
You can see this immediately in two ways:
• The determinant would have one row that was a linear combination of the
others. You’ll remember that by doing row operations, you could reach a row
of zeros, and so the determinant is zero.
• The parallelopiped would have zero volume if squashed flat. In this case all
three vectors lie in a plane, and so any one is a linear combination of the
other two. (Figure 2.1b.)
2.1. TRIPLE AND MULTIPLE PRODUCTS 19
2.1.4 Vector triple product

This is defined as the vector product of a vector with a vector product, a × (b × c).
Now, the vector triple product a × (b × c) must be perpendicular to (b × c), which
in turn is perpendicular to both b and c. Thus a × (b × c) can have no component
perpendicular to b and c, and hence must be coplanar with them. It follows that
the vector triple product must be expressible as a linear combination of b and c:
a × (b × c) = λb + µc .
The values of the coefficients can be obtained by multiplying out in component
form. Only the first component need be evaluated, the others then being obtained
by symmetry. That is
(a × (b × c))1 = a2 (b × c)3 − a3(b × c)2

= a2 (b1c2 − b2c1 ) + a3(b1 c3 − b3 c1)
= (a2 c2 + a3 c3)b1 − (a2b2 + a3b3 )c1
= (a1 c1 + a2 c2 + a3 c3 )b1 − (a1b1 + a2b2 + a3 b3)c1
= (a · c)b1 − (a · b)c1
The equivalents must be true for the 2nd and 3rd components, so we arrive at the
identity
a × (b × c) = (a · c)b − (a · b)c .
a
bxc
In arbitrary direction
a x (bx c)
Figure 2.2: Vector triple product.
2.1.5 Projection using vector triple product

An example of the application of this formula is as follows. Suppose v is a vector
and we want its projection into the x y -plane. The component of v in the z
direction is v · k̂k , so the projection we seek is v − (v · k̂k )k̂k . Writing k̂k ← a, v ← b,
k̂k ← c,
a × (b × c) = (a · c)b − (a · b)c
↓ ↓ ↓
k̂k × (v × k̂k ) = (k̂k · k̂k )v − (k̂k · v)k̂k
= v − (v · k̂k )k̂k
So v − (v · k̂k )k̂k = k̂k × (v × k̂k ).

(Hot stuff! But the expression v − (v · k̂k )k̂k is much easier to understand, and
cheaper to compute!)
2.1.6 Other repeated products

Many combinations of vector and scalar products are possible, but we consider only
one more, namely the vector quadruple product (a × b) × (c × d). By regarding
a × b as a single vector, we see that this vector must be representable as a linear
combination of c and d. On the other hand, regarding c × d as a single vector, we
see that it must also be a linear combination of a and b. This provides a means
of expressing one of the vectors, say d, as linear combination of the other three,
as follows:
(a × b) × (c × d) = [(a × b) · d]c − [(a × b) · c]d
= [(c × d) · a]b − [(c × d) · b]a
Hence
[(a × b) · c] d = [(b × c) · d] a + [(c × a) · d] b + [(a × b) · d] c
or
[(b × c) · d] a + [(c × a) · d] b + [(a × b) · d] c
d= = αa + βb + γc .
[(a × b) · c]
This is not something to remember off by heart, but it is worth remembering that
the projection of a vector on any arbitrary basis set is unique.
♣ Example
Q1 Use the quadruple vector product to express the vector d = [3, 2, 1] in terms
of the vectors a = [1, 2, 3], b = [2, 3, 1] and c = [3, 1, 2].
A1 Grinding away at the determinants, we find
[(a×b)·c] = −18; [(b×c)·d] = 6; [(c×a)·d] = −12; [(a×b)·d] = −12
So, d = (−a + 2b + 2c)/3.
2.2. GEOMETRY USING VECTORS: LINES, PLANES 21
a
d
Figure 2.3: The projection of a (3-)vector onto a set of (3) basis vectors is unique. Ie in d =
αa + βb + γc, the set {α, β, γ} is unique.
2.2 Geometry using vectors: lines, planes

2.2.1 The equation of a line
The equation of the line passing through the point whose position vector is a and
lying in the direction of vector b is
r = a + λb
where λ is a scalar parameter. If you make b a unit vector, r = a + λb̂ then λ will
represent metric length.
For a line defined by two points a1 and a2
r = a1 + λ(a2 − a1)
or for the unit version
(a2 − a1)
r = a1 + λ
|a2 − a1|
^
λb
Point r traces
a out line.
r
Figure 2.4: Equation of a line. With b̂ a unit vector, λ is in the length units established by the
definition of a.
2.2.2 The shortest distance from a point to a line

Referring to Figure 2.5(a) the vector p from c to any point on the line is p =
a + λb̂ − c = (a − c) + λb̂ which has length squared p 2 = (a − c)2 + λ2 + 2λ(a −
c) · b̂ . Rather than minimizing length, it is easier to minimize length-squared. The
minumum is found when d p 2/d λ = 0, ie when
λ = −(a − c) · b̂ .
So the minimum length vector is
p = (a − c) − ((a − c) · b̂)b̂.
You might spot that is the component of (a−c) perpendicular to b̂ (as expected!).
Furthermore, using the result of Section 2.1.5,
p = b̂ × [(a − c) × b̂] .
Because b̂ is a unit vector, and is orthogonal to [(a − c) × b̂], the modulus of the
vector can be written rather more simply as just
pmin = |(a − c) × b̂| .
c
λb λb
P
a
r-c a µd
r Q
(a) (b)
Figure 2.5: (a) Shortest distance point to line. (b) Shortest distance, line to line.
2.2.3 The shortest distance between two straight lines

If the shortest distance between a point and a line is along the perpendicular, then
the shortest distance between the two straight lines r = a + λb̂ and r = c + µd̂
must be found as the length of the vector which is mutually perpendicular to the
lines.
The unit vector along the mutual perpendicular is
p̂ = (b̂ × d̂)/|b̂ × d̂| .
(Yes, don’t forget that b̂ × d̂ is NOT a unit vector. b̂ and d̂ are not orthogonal,
so there is a sin θ lurking!)
The minimum length is therefore the component of a − c in this direction

pmin = (a − c) · (b̂ × d̂)/|b̂ × d̂| .
2.2. GEOMETRY USING VECTORS: LINES, PLANES 23
♣ Example
Q Two long straight pipes are specified using Cartesian co-ordinates as follows:
Pipe A has diameter 0.8 and its axis passes through points (2, 5, 3) and
(7, 10, 8).
Pipe B has diameter 1.0 and its axis passes through the points (0, 6, 3) and
(−12, 0, 9).
Determine whether the pipes need to be realigned to avoid intersection.
A Each pipe axis is defined using two points. The vector equation of the axis
of pipe A is
√
r = [2, 5, 3] + λ′ [5, 5, 5] = [2, 5, 3] + λ[1, 1, 1]/ 3
The equation of the axis of pipe B is
√
r = [0, 6, 3] + µ′ [12, 6, 6] = [0, 6, 3] + µ[−2, −1, 1]/ 6
The perpendicular to the two axes has direction

ı̂ı 
̂ k
k̂

[1, 1, 1] × [−2, −1, 1] = 1 1 1 = [2, −3, 1] = p

−2 −1 1
The length of the mutual perpendicular is
[2, −3, 1] [2, −3, 1]
(a − c) · √ = [2, −1, 0] · √ = 1.87 .
14 14
But the sum of the radii of the two pipes is 0.4 + 0.5 = 0.9. Hence the pipes
do not intersect.
2.2.4 The equation of a plane
There are a number of ways of specifying the equation of a plane.
1. If b and c are two non-parallel vectors (ie b × c 6= 0), then the equation of
the plane passing through the point a and parallel to the vectors b and c may
be written in the form
r = a + λb + µc
where λ, µ are scalar parameters. Note that b and c are free vectors, so don’t
have to lie in the plane (Figure 2.6(a).)
2. Figure 2.6(b) shows the plane defined by three non-collinear points a, b and
c in the plane (note that the vectors b and c are position vectors, not free
vectors as in the previous case). The equation might be written as
r = a + λ(b − a) + µ(c − a)
3. Figure 2.6(c) illustrates another description is in terms of the unit normal to
the plane n̂ and a point a in the plane
r · n̂ = a · n̂ .
r ^
n
r c
b r
a
c
a a
NB that these are

O b parallel to the plane, not
necessarily in the plane
O O
(a) (b) (c)
Figure 2.6: (a) Plane defined using point and two lines. (b) Plane defined using three points. (c)
Plane defined using point and normal. Vector r is the position vector of a general point in the
plane.
2.2.5 The shortest distance from a point to a plane

The shortest distance from a point d to the plane is along the perpendicular.
Depending on how the plane is defined, this can be written as
|(d − a) · (b × c)|
D = |(d − a) · n̂| or D= .
|b × c|
2.3 Solution of vector equations

It is sometimes required to obtain the most general vector which satisfies a given
vector relationship. This is very much like obtaining the locus of a point. The best
method of proceeding in such a case is as follows:
(i) Decide upon a system of three co-ordinate vectors using two non-parallel vectors
appearing in the vector relationship. These might be a, b and their vector product
(a × b).
(ii) Express the unknown vector x as a linear combination of these vectors
x = λa + µb + νa × b
where λ, µ, ν are scalars to be found.
(iii) Substitute the above expression for x into the vector relationship to determine
the constraints on λ, µ and ν for the relationship to be satisfied.
♣ Example
Q Solve the vector equation x = x × a + b.
A Step (i): Set up basis vectors a, b and their vector product a × b.
Step (ii): x = λa + µb + νa × b.
2.4. ROTATION, ANGULAR VELOCITY/ACCELERATION AND MOMENTS 25
Step (iii): Bung this expression for x into the equation!
λa + µb + νa × b = (λa + µb + νa × b) × a + b
= 0 + µ(b × a) + ν(a × b) × a + b
= −ν(a · b)a + (νa2 + 1)b − µ(a × b)
We have learned that any vector has a unique expression in terms of a basis
set, so that the coefficients of a, b and a × b on either side of the equation
much be equal.
⇒ λ = −ν(a · b)
µ = νa2 + 1
ν = −µ
so that
1 1 a·b
µ= ν=− λ= .
1 + a2 1 + a2 1 + a2
So finally the solution is the single point:
1
x= ((a · b)a + b − (a × b))
1 + a2
2.4 Rotation, angular velocity/acceleration and moments

A rotation can be represented by a vector whose direction is along the axis of
rotation in the sense of a r-h screw, and whose magnitude is proportional to the
size of the rotation (Fig. 2.7). The same idea can be extended to the derivatives,
ω.
that is, angular velocity ω and angular acceleration ω̇
Angular accelerations arise because of a moment (or torque) on a body. In me-
chanics, the moment of a force F about a point Q is defined to have magnitude
M = F d , where d is the perpendicular distance between Q and the line of action
L of F.
The vector equation for moment is
M=r×F
where r is the vector from Q to any point on the line of action L of force F.
The resulting angular acceleration vector is in the same direction as the moment
vector.
The instantaneous velocity of any point P on a rigid body undergoing pure rotation
can be defined by a vector product as follows. The angular velocity vector ω has
ω
in right−hand screw sense
dθ δt
dt
Figure 2.7: The angular velocity vector ω is along the axis of rotation and has magnitude equal to
the rate of rotation.
magnitude equal to the angular speed of rotation of the body and with direction
the same as that of the r-h screw. If r is the vector OP , where the origin O can
be taken to be any point on the axis of rotation, then the velocity v of P due to
the rotation is given, in both magnitude and direction, by the vector product
v = ω × r.
ω
M
r r
α v
F
Figure 2.8:
Revised Oct 2013

Lecture 3
Differentiating Vector Functions of a Single

Variable
Your experience of differentiation and integration has extended as far as scalar

functions of single and multiple variables — functions like f (x ) and f (x , y , t).
It should be no great surprise that we often wish to differentiate vector func-
tions. For example, suppose you were driving along a wiggly road with position
r(t) at time t. Differentiating r(t) wrt time should yield your velocity v(t), and
differentiating v(t) should yield your acceleration. Let’s see how to do this.
3.1 Differentiation of a vector

The derivative of a vector function a(p) of a single parameter p is
a(p + δp) − a(p)
a′ (p) = lim .
δp→0 δp
If we write a in terms of components relative to a FIXED coordinate system (ı̂ı ,̂ , k̂k
constant)
a(p) = a1 (p)ı̂ı + a2 (p)̂ + a3(p)k̂k
then
d a1 d a2 d a3
a′ (p) = ı̂ı + ̂ + k̂k .
dp dp dp
That is, in order to differentiate a vector function, one simply differentiates each
component separately. This means that all the familiar rules of differentiation
apply, and they don’t get altered by vector operations like scalar product and
vector products.
Thus, for example:
d da db d da db
(a × b) = ×b+a× (a · b) = ·b+a· .
dp dp dp dp dp dp
27
28 LECTURE 3. DIFFERENTIATING VECTOR FUNCTIONS OF A SINGLE VARIABLE
Note that d a/d p has a different direction and a different magnitude from a.
Likewise, as you might expect, the chain rule still applies. If a = a(u) and u = u(t),
say:
d da du
a=
dt du dt
♣ Examples
Q A 3D vector a of constant magnitude is varying over time. What can you say
about the direction of ȧ?
A Using intuition: if only the direction is changing, then the vector must be
tracing out points on the surface of a sphere. We would guess that the
derivative ȧ is orthogonal to a.
To prove this write
d da da da
(a · a) = a · + · a = 2a · .
dt dt dt dt
But (a · a) = a2 which we are told is constant. So
d da
(a · a) = 0 ⇒ 2a · =0
dt dt
and hence a and d a/d t must be perpendicular.
Q The position of a vehicle is r(u) where u is the amount of fuel consumed by
some time t. Write down an expression for the acceleration.
A The velocity is
dr dr du
v= =
dt du dt
2
d 2r d r d 2u

d dr du
a= = +
dt dt d u2 dt d u d t2
3.1.1 Geometrical interpretation of vector derivatives

Let r(p) be a position vector tracing a space curve as some parameter p varies.
The vector δr is a secant to the curve, and δr/δp lies in the same direction. (See
Fig. 3.1.) In the limit as δp tends to zero δr/δp = d r/d p becomes a tangent to
the space curve. If the magnitude of this vector is 1 (i.e. a unit tangent), then
3.1. DIFFERENTIATION OF A VECTOR 29
|d r| = d p so the parameter p is arc-length (metric distance). More generally,

however, p will not be arc-length and we will have:
dr dr ds
=
dp ds dp
So, the direction of the derivative is that of a tangent to the curve, and its
magnitude is |d s/d p|, the rate of change of arc length w.r.t the parameter.
Of course if that parameter p is time, the magnitude |d r/d t| is the speed.
♣ Example
Q Draw the curve
s s hs
r = a cos( √ 2 ı + a sin( √
)ı̂  + √
)̂ k̂k
a + h2 a2 + h 2 a2 + h 2
where s is arc length and h, a are constants. Show that the tangent d r/d s
to the curve has a constant elevation angle w.r.t the x y -plane, and determine
its magnitude.
A
dr a a h
= −√ sin ()ı̂ı + √ cos ()̂ + √ k̂k
ds a2 + h 2 a2 + h 2 a2 + h 2
√
The projection
√ on the x y plane has magnitude a/ a2 + h2 and in the z
direction h/ a2 + h2, so the elevation angle is a constant, tan−1(h/a).
We are expecting d r/d s = 1, and indeed
q p
2
a sin () + a cos () + h / a2 + h2 = 1.
2 2 2 2
3.1.2 Arc length is a special parameter!

It might seem that we can be completely relaxed about saying that any old pa-
rameter p is arc length, but this is not the case. Why not? The reason is that arc
length is special is that, whatever the parameter p,
Z p
dr
s= dp .
p0 d p

Perhaps another way to grasp the significance of this is using Pythagoras’ theorem
on a short piece of curve: in the limit as d x etc tend to zero,
ds2 = dx2 + dy 2 + dz2 .
δr δr
r (p) r (s)
r (p + δ p) r (s + δ s)
dr ds dr
dp dp ds
1
Figure 3.1: Left: δr is a secant to the curve but, in the limit as δp → 0, becomes a tangent.
Right: if the parameter is arc length s, then |dr| = ds.
So if a curve is parameterized in terms of p

s
ds dx 2 dy 2 dz 2
= + + .
dp dp dp dp
As an example, suppose in our earlier example we had parameterized our helix as
r = a cos pı̂ı + a sin p̂ + hp k̂k
It would be easy just to say that p was arclength, but it would not be correct
because
s
ds dx 2 dy 2 dz 2
= + +
dp dp dp dp
q p
2
= a sin p + a cos p + h = a2 + h2
2 2 2 2
√
If p really was arclength, d s/d p = 1. So p/ a2 + h2 is arclength, not p.
3.2 Integration of a vector function

The integration of a vector function of a single scalar variable can be regarded
simply as the reverse of differentiation. In other words
Z p2
d a(p)
dp
p1 dp
For example the integral of the acceleration vector of a point over an interval of
time is equal to the change in the velocity vector during the same time interval.
However, many other, more interesting and useful, types of integral are possible,
especially when the vector is a function of more than one variable. This requires
the introduction of the concepts of scalar and vector fields. See later!
3.3. CURVES IN 3 DIMENSIONS 31
3.3 Curves in 3 dimensions

In the examples above, parameter p has been either arc length s or time t. It
doesn’t have to be, but these are the main two of interest. Later we shall look
at some important results when differentiating w.r.t. time, but now let use look
more closely at 3D curves defined in terms of arc length, s.
Take a piece of wire, and bend it into some arbitrary non-planar curve. This is a
space curve. We can specify a point on the wire by specifying r(s) as a function
of distance or arc length s along the wire.
3.3.1 The Frénet-Serret relationships

We are now going to introduce a local orthogonal coordinate frame for each point
s along the curve, ie one with its origin at r(s). To specify a coordinate frame we
need three mutually perpendicular directions, and these should be intrinsic to the
curve, not fixed in an external reference frame. The ideas were first suggested by
two French mathematicians, F-J. Frénet and J. A. Serret.
1. Tangent t̂
There is an obvious choice for the first direction at the point r(s), namely the
unit tangent t̂. We already know that
d r(s)
t̂ =
ds
2. Principal Normal n̂
Recall that earlier we proved that if a was a vector of constant magnitude
that varies in direction over time then d a/d t was perpendicular to it. Because
t̂ has constant magnitude but varies over s, d t̂/d s must be perpendicular to
t̂.
Hence the principal normal n̂ is
d t̂
= κn̂ : where κ ≥ 0 .
ds
κ is the curvature, and κ = 0 for a straight line. The plane containing t̂ and
n̂ is called the osculating plane.
3. The Binormal b̂
The local coordinate frame is completed by defining the binormal
b̂(s) = t̂(s) × n̂(s) .

Since b̂ · t̂ = 0,
d b̂ d t̂ d b̂
· t̂ + b̂ · = · t̂ + b̂ · κn̂ = 0
ds ds ds
from which
d b̂
· t̂ = 0.
ds
But this means that d b̂/d s is along the direction of n̂, or
d b̂
= −τ(s)n̂(s)
ds
where τ is the torsion, and the negative sign is a matter of convention.
Differentiating n̂ · t̂ = 0 and n̂ · b̂ = 0, we find
d n̂
= −κ(s)t̂(s) + τ(s)b̂(s).
ds
The Frénet-Serret relationships:
d t̂/d s = κn̂
d n̂/d s = −κ(s)t̂(s) + τ(s)b̂(s)
d b̂/d s = −τ(s)n̂(s)
♣ Example
Q Derive κ(s) and τ(s) for the helix

s s s p
r(s) = a cos ı̂ı + a sin ̂ + h k̂k ; β= a2 + h 2
β β β
and comment on their values.
A We found the unit tangent earlier as

dr a s a s h
t̂ = = − sin , cos , .
ds β β β β β
Differentiation gives

d t̂ a s a s
κn̂ = = − 2 cos , − 2 sin , 0
ds β β β β
3.4. RADIAL AND TANGENTIAL COMPONENTS IN PLANE POLARS 33
Curvature is always positive, so

a s s
κ= 2 n̂ = − cos , − sin , 0 .
β β β
So the curvature is constant, and the normal is parallel to the x y -plane.
Now use

ı
ı̂ 
̂ k
k̂
h s h s a
b̂ = t̂×n̂ = (−a/β)S (a/β)C (h/β) = sin , − cos ,
β β β β β
−C −S 0
and differentiate b̂ to find an expression for the torsion

d b̂ h s h s −h
= cos , sin , 0 = n̂
ds β2 β β2 β β2
so the torsion is
h
τ=
β2
again a constant.
3.4 Radial and tangential components in plane polars

In plane polar coordinates, the radius vector
of any point P is given by
r = r cos θı̂ı + r sin θ̂
= rêr êθ
where we have introduced the unit radial vec- êr
tor
P
êr = cos θı̂ı + sin θ̂ . ̂ r
θ
The other “natural” (we’ll see why in a later
lecture) unit vector in plane polars is orthog- ı̂ı
onal to êr and is
êθ = − sin θı̂ı + cos θ̂
so that êr · êr = êθ · êθ = 1 and êr · êθ = 0.

Now suppose P is moving so that r is a function of time t. Its velocity is
d dr dêr
ṙ = (rêr ) = êr + r
dt dt dt
dr dθ
= êr + r (− sin θı̂ı + cos θ̂ )
dt dt
dr dθ
= êr + r êθ
dt dt
= radial + tangential
The radial and tangential components of velocity of P are therefore d r /d t and

r d θ/d t, respectively.
Differentiating a second time gives the acceleration of P
d 2r dr dθ dr dθ d 2θ dθ dθ
r̈ = 2
êr + êθ + êθ + r êθ − r êr
"d t dt dt # dt dt d t2 dt dt
2
d 2r d 2θ

dθ dr dθ
= −r êr + 2 + r 2 êθ
d t2 dt dt dt dt
3.5. ROTATING SYSTEMS 35
3.5 Rotating systems

Consider a body which is rotating with constant angular velocity ω about some
axis passing through the origin. Assume the origin is fixed, and that we are sitting
in a fixed coordinate system Ox y z .
If ρ is a vector of constant magnitude and constant direction in the rotating system,
then its representation r in the fixed system must be a function of t.
r(t) = R(t)ρρ
At any instant as observed in the fixed system

dr
= Ṙρρ + Rρ̇ρ
dt
but the second term is zero since we assumed ρ to be constant so we have
dr
= ṘR⊤r
dt
Note that:
• d r/d t will have fixed magnitude;

• d r/d t will always be perpendicular to the axis of rotation;
• d r/d t will vary in direction within those constraints;
• r(t) will move in a plane in the fixed system.
ρ
Now let’s consider the term ṘR⊤. First, note that RR⊤ = I (the identity), so
differentiating both sides yields
⊤
ṘR⊤ + RṘ = 0
⊤
ṘR⊤ = −RṘ
Thus ṘR⊤ is anti-symmetric:

 
0 −z y
ṘR⊤ =  z 0 −x 
−y x 0
Now you can verify for yourself that application of a matrix of this form to an
arbitrary vector has precisely the same effect as the cross product operator, ω ×,
where ω = [x y z ]⊤. Loh-and-behold, we then we have
ṙ = ω × r
matching the equation at the end of lecture 2, v = ω ×r, as we would hope/expect.
3.5.1 Rotation: Part 2

Now suppose ρ is the position vector of a point P which moves in the rotating
frame. There will be two contributions to motion with respect to the fixed frame,
one due to its motion within the rotating frame, and one due to the rotation itself.
So, returning to the equations we derived earlier:
r(t) = R(t)ρρ(t)
and the instantaenous differential with respect to time:
dr
= Ṙρρ + Rρ̇ρ = ṘR⊤r + Rρ̇ρ
dt
Now ρ is not constant, so its differential is not zero; hence rewriting this last
equations we have that
The instantaneous velocity of P in the fixed frame is
dr
= Rρ̇ρ + ω × r
dt
The second term of course, is the contribution from the rotating frame which we
saw previously. The first is the linear velocity measured in the rotating frame ρ̇ρ,
referred to the fixed frame (via the rotation matrix R which aligns the two frames)
P at t+ δ t
ω δρ δr
P at t
(ω r) δ t
r= ρ at t
3.5.2 Rotation 3: Instantaneous acceleration

Our previous result is a general one relating the time derivatives of any vector in
rotating and non-rotating frames. Let us now consider the second differential:
ω × r + ω × ṙ + Ṙρ̇ρ + Rρ̈ρ
r̈ = ω̇
We shall assume that the angular acceleration is zero, which kills off the first term,
and so now, substituting for ṙ we have
r̈ = ω × r + Rρ̇ρ) + Ṙρ̇ρ + Rρ̈ρ

ω × (ω
= ω × r) + ω × Rρ̇ρ + Ṙρ̇ρ + Rρ̈ρ
ω × (ω
= ω × r) + ω × Rρ̇ρ + Ṙ(R⊤ R)ρ̇ρ + Rρ̈ρ
ω × (ω
= ω × (ω
ω × r) + 2ω ω × (Rρ̇ρ) + Rρ̈ρ
The instantaneous acceleration is therefore
r̈ = Rρ̈ρ + 2ω
ω × (Rρ̇ρ) + ω × (ω
ω × r)
• The first term is the acceleration of the point P in the rotating frame mea-
sured in the rotating frame, but referred to the fixed frame by the rotation
R
• The last term is the centripetal acceleration to due to the rotation. (Yes! Its
magnitude is ω 2 r and its direction is that of −r. Check it out.)
γt
ω = ωm̂
r
m̂
ℓ̂ℓ
n̂ γ = γ ℓ̂ℓ
Figure 3.2: Coriolis example.
• The middle term is an extra term which arises because of the velocity of P
in the rotating frame. It is known as the Coriolis acceleration, named after
the French engineer who first identified it.
Because of the rotation of the earth, the Coriolis acceleration is of great im-
portance in meteorology and accounts for the occurrence of high pressure anti-
cyclones and low pressure cyclones in the northern hemisphere, in which the Coriolis
acceleration is produced by a pressure gradient. It is also a very important compo-
nent of the acceleration (hence the force exerted) by a rapidly moving robot arm,
whose links whirl rapidly about rotary joints.
♣ Example
Q Find the instantaneous acceleration of a projectile fired along a line of longi-
tude (with angular velocity of γ constant relative to the sphere) if the sphere
is rotating with angular velocity ω .
A Consider a coordinate frame defined by mutually orthogonal unit vectors,

ℓ̂ℓ, m̂ and n̂, as shown in Fig. 3.2. We shall assume, without loss of generality,
that the fixed and rotating frames are instantaneously aligned at the moment
shown in the diagram, so that R = I, the identity, and hence r = ρ .
In the rotating frame
ρ̇ρ = γ × ρ and ρ̈ρ = γ × ρ̇ρ = γ × (γ
γ × ρ)
So the in the fixed reference frame, because these two frames are instanta-
neously aligned
γ × ρ ) + 2ω
r̈ = γ × (γ ω × (γ
γ × ρ ) + ω × (ω
ω × r) .
The first term is the centripetal acceleration due to the projectile moving
around the sphere — which it does because of the gravitational force. The
last term is the centripetal acceleration resulting from the rotation of the
sphere. The middle term is the Coriolis acceleration.
Using Fig. 3.2, at some instant t
r(t) = ρ (t) = r cos(γt)m̂ + r sin(γt)n̂
and
γ = γ ℓ̂ℓ
Then
γ × ρ ) = (γ
γ × (γ γ − γ 2 ρ = −γ 2 ρ = −γ 2 r,
γ · ρ )γ
Check the direction — the negative sign means it points towards the centre
of the sphere, which is as expected.
Likewise the last term can be obtained as
ω × r) = −ω 2 r sin(γt)n̂
ω × (ω
Note that it is perpendicular to the axis of rotation m̂, and because of the
minus sign, directed towards the axis)
The Coriolis term is derived as:
ω × ρ̇ρ = 2ω
2ω ω × (γγ × ρ)
     
0 γ 0
= 2 ω  ×  0  × r cos γt 
0 0 r sin γt
= 2ωγr cos γt ℓ̂ℓ
Instead of a projectile, now consider a rocket on rails which stretch north

from the equator. As the rocket travels north it experiences the Coriolis force
(exerted by the rails):
2 γ ω R cos γt ℓ̂ℓ
+ve -ve +ve +ve
Hence the coriolis force is in the direction opposed to ℓ̂ℓ (i.e. in the opposite
direction to the earth’s rotation). In the absence of the rails (or atmosphere)
the rocket’s tangetial speed (relative to the surface of the earth) is greater
than the speed of the surface of the earth underneath it (since the radius
of successive lines of latitude decreases) so it would (to an observer on the
earth) appear to deflect to the east. The rails provide a coriolis force keeping
it on the same meridian.
Revised Oct 2013

Rocket’s velocity in direction of meridian

Tangential velocity of earth’s surface
Tangential component of velocity

(NB instantaneously common to earth’s surface and rocket)
Figure 3.3: Rocket example

Figure 3.4: Coriolis effect giving rise to weather systems

Lecture 4
Line, Surface and Volume Integrals.

Curvilinear coordinates.
We started off the course being concerned with individual vectors a, b, c, and so
on.
We went on to consider how single vectors vary over time or over some other
parameter such as arc length.
In much of the rest of the course, we will be concerned with scalars and vectors
which are defined over regions in space — scalar and vector fields
In this lecture we introduce line, surface and volume integrals, and consider how
these are defined in non-Cartesian, curvilinear coordinates
4.1 Scalar and vector fields

When a scalar function u(r) is determined or defined at each position r in some
region, we say that u is a scalar field in that region.
Similarly, if a vector function v(r) is defined at each point, then v is a vector field
in that region. As you will see, in field theory our aim is to derive statements about
the bulk properties of scalar and vector fields, rather than to deal with individual
scalars or vectors.
Familiar examples of each are shown in figure 4.1.
In Lecture 1 we worked out the force F(r) on a charge Q arising from a number
of charges qi . The electric field is F/Q, so
N
X qi 1
E(r) = K (r − ri ) . (K = )
|r − ri |3 4πǫr ǫ0
i =1
For example; you could work out the velocity field, in plane polars, at any point on
43
44LECTURE 4. LINE, SURFACE AND VOLUME INTEGRALS. CURVILINEAR COORDINATES.
(a) (b)
Figure 4.1: Examples of (a) a scalar field (pressure); (b) avector field (wind velocity)
a wheel spinning about its axis

v(r) = ω × r
or the fluid flow field around a wing.
If the fields are independent of time, they are said to be steady. Of course, most
vector fields of practical interest in engineering science are not steady, and some
are unpredictable.
Let us first consider how to perform a variety of types of integration in vector and
scalar fields.
4.2 Line integrals through fields

Line integrals are concerned with measuring the integrated interaction with a field
as you move through it on some defined path. Eg, given a map showing the
pollution density field in Oxford, you may wish to work out how much pollution
you breathe in when cycling from college to the Department via different routes.
First recall the definition of an integral for a scalar function f (x ) of a single scalar
variable x . One assumes a set of n samples fi = f (xi ) spaced by δxi . One forms
the limit of the sum of the products f (xi )δxi as the number of samples tends to
infinity
Z Xn
f (x )d x = lim fi δxi .
n → ∞ i =1
δxi → 0
For a smooth function, it is irrelevant how the function is subdivided.
4.2.1 Vector line integrals

In a vector line integral, the path L along which the integral is to be evaluated
is split into a large number of vector segments δri . Each line segment is then
4.2. LINE INTEGRALS THROUGH FIELDS 45
F(r)
r δr
Figure 4.2: Line integral. In the diagram F(r) is a vector field, but it could be replaced with scalar
field U(r).
multiplied by the quantity associated with that point in space, the products are
then summed and the limit taken as the lengths of the segments tend to zero.
There are three types of integral we have to think about, depending on the nature
of the product:
1. Integrand U(r) is a scalar field, hence the integral is a vector.

Z !
X
I = U(r)d r = lim Ui δ r i .
L δri →0
i
2. Integrand a(r) is a vector field dotted with d r hence the integral is a scalar:
Z !
X
I = a(r) · d r = lim ai · δ r i .
L δri →0
i
3. Integrand a(r) is a vector field crossed with d r hence vector result.

Z !
X
I = a(r) × d r = lim ai × δ r i .
L δri →0
i
Note immediately that unlike an integral in a single scalar variable, there are many
paths L from start point rA to end point rB , and the integral will in general depend
on the path taken.
Physical examples of line integrals

• The total work done by a force F as it moves a point from A to B along
a given path C is given by a line integral of type 2 above. If the force acts
at point r and the instantaneous displacement along curve C is d r then the

infinitessimal work done is d W = F.d r, and so the total work done traversing
the path is
Z
WC = F.d r
C
• Ampère’s law relating magnetic field B to linked current can be written as

I
B.d r = µ0 I
C
where I is the current enclosed by (closed) path C.
• The force on an element of wire carrying current I, placed in a magnetic field
of strength B, is d F = Id r × B. So if a loop this wire C is placed in the field
then the total force will be and integral of type 3 above:
I
F = I dr × B
C
Note that the expressions above are beautifully compact in vector notation, and are
all independent of coordinate system. Of course when evaluating them we need
to choose a coordinate system: often this is the standard Cartesian coordinate
system (as in the worked examples below), but need not be, as we shall see in
section 4.6.
♣ Examples
Q1 An example in the x y -plane. A force F = x 2yı̂ı + x y 2̂ acts on a body as it
moves between (0, 0) and (1, 1).
Determine the work done when the path is
1. along the line y = x .
2. along the curve y = x n , n > 0.
3. along the x axis to the point (1, 0) and then along the line x = 1.
A1 This is an example of the “type 2” line integral. In planar Cartesians, d r =
ı̂ı d x + ̂ d y . Then the work done is
Z Z
F · d r = (x 2y d x + x y 2d y ) .
L L
1. For the path y = x we find that d y = d x . So it is easiest to convert all

y references to x .
Z (1,1) Z x=1 Z x=1
2 2 2 2
x=1
2x 3d x = x 4/2x=0 = 1/2 .

(x y d x +x y d y ) = (x x d x +x x d x ) =
(0,0) x=0 x=0
4.2. LINE INTEGRALS THROUGH FIELDS 47
1,1
1
2 3
0,0 0,1
Figure 4.3: Line integral taken along three difference paths
2. For the path y = x n we find that d y = nx n−1d x , so again it is easiest to

convert all y references to x .
Z (1,1) Z x=1
2 2
(x y d x + x y d y ) = (x n+2d x + nx n−1.x .x 2n d x )
(0,0) x=0
Zx=1
= (x n+2d x + nx 3n d x )
x=0
1 n
= +
n + 3 3n + 1
3. This path is not smooth, so break it into two. Along the first section,
y = 0 and d y = 0, and on the second x = 1 and d x = 0, so
Z B Z x=1 Z y =1 y =1
2 2 2
1y 2d y = 0+ y 3/3y =0 = 1/3 .

(x y d x +x y d y ) = (x 0d x )+
A x=0 y =0
So in general the integral depends on the path taken. Notice that answer (1)
is the same as answer (2) when n = 1, and that answer (3) is the limiting
value of answer (2) as n → ∞.
Q2 Repeat part (2) using the Force F = x y 2ı̂ı + x 2y̂ .

A2 For the path y = x n we find that d y = nx n−1d x , so

Z (1,1) Z x=1
2 2
(y x d x + y x d y ) = (x 2n+1d x + nx n−1.x 2.x n d x )
(0,0) x=0
Zx=1
= (x 2n+1d x + nx 2n+1d x )
x=0
1 n
= +
2n + 2 2n + 2
1
= independent of n
2
4.3 Line integrals in Conservative fields

In the second example, the line integral has the same value for the whole range
of paths. In fact it is wholly independent of path. This is easy to see if we write
g(x , y ) = x 2y 2/2. Then using the definition of the perfect differential
∂g ∂g
dg = dx + dy
∂x ∂y
we find that
Z B Z B
2 2
(y x d x + y x d y ) = dg
A A
= gB − gA
which depends solely on the value of g at the start and end points, and not at all
on the path used to get from A to B. Such a vector field is called conservative.
One sort of line integral performs the integration around a complete loop and is
denoted with a ring. If E is a conservative field, determine the value of
I
E · dr .
In electrostatics, if E is the electric field then the potential function is

Z
φ = − E · dr .
Do you think E is conservative?

4.4. SURFACE INTEGRALS 49
4.3.1 A note on line integrals defined in terms of arc length

Line integrals are often defined in terms of scalar arc length. They don’t appear
to involve vectors (but actually they are another form of type 2 defined earlier).
The integrals usually appear as follows
Z
I = F (x , y , z )d s
L
and most often the path L is along a curve defined parametrically as x = x (p),
y = y (p), z = z (p) where p is some parameter. Convert the function to F (p),
writing
Z pend
ds
I= F (p) dp
pstart dp
where
" 2 2 2#1/2
ds dx dy dz
= + + .
dp dp dp dp
Note that the parameter p could be arc-length s itself, in which case d s/d p = 1
of course! Another possibility is that the parameter p is x — that is we are told
y = y (x ) and z = z (x ). Then
Z xend " 2 2 #1/2
dy dz
I= F (x ) 1 + + dx .
xstart dx dx
4.4 Surface integrals

These can be defined by analogy with line integrals.
The surface S over which the integral is to be evaluated is now divided into in-
finitesimal vector elements of area d S, the direction of the vector d S representing
the direction of the surface normal and its magnitude representing the area of the
element.
Again there are three possibilities:
R
• S Ud S — scalar field U; vector integral.
R
• S a · d S — vector field a; scalar integral.
R
• S a × d S — vector field a; vector integral.
R
(in addition, of course, to the purely scalar form, S Ud S).
Physical example of surface integral
• Physical examples of surface integrals with

R vectors often involve the idea of
flux of a vector field through a surface, S a.d S For example the mass of fluid
crossing a surface S in time d t is d M = ρv.d Sd t where ρ(r) is the fluid
density and v(r) is the fluid velocity. The total mass flux can be expressed as
a surface integral:
Z
ΦM = ρ(r )v(r).d S
S
Again, though this expression is coordinate free, we evaluate an example below

using Cartesians. Note, however, that in some problems, symmetry may lead us
to a different more natural coordinate system.
♣ Example
R
Evaluate F · d S over the x = 1 side of
the cube shown in the figure when F =
yı̂ı + z̂ + x k̂k .
d S is perpendicular to the surface. Its ±
direction actually depends on the nature 11111
00000
00000
11111
of the problem. More often than not,
z 00000
11111
the surface will enclose a volume, and the 00000
11111
00000
11111
surface direction is taken as everywhere 1 y 11111
00000
emanating from the interior.
1
00000
11111
00000
11111
Hence for the x = 1 face of the cube 00000
11111
00000
11111
dS = dydz i
d S = d y d zı̂ı 00000
11111
00000
11111
1 x
and
Z Z Z
F · dS = ydydz
1 21 1 1
= y 0 z |0 = .
2 2
4.5 Volume integrals

The definition of the volume integral is again taken as the limit of a sum of products
as the size of the volume element tends to zero. One obvious difference though is
that the element of volume is a scalar (how could you define a direction with an
infinitesimal volume element?). The possibilities are:
4.6. CHANGING VARIABLES: CURVILINEAR COORDINATES 51
R
• V U(r)d V — scalar field; scalar integral.
R
• V ad V — vector field; vector integral.
You have covered these (more or less) in your first year course, so not much more
to say here. The next section considers these again in the context of a change of
coordinates.
4.6 Changing variables: curvilinear coordinates

Up to now we have been concerned with Cartesian coordinates x , y , z with coor-
dinate axes ı̂ı ,̂ , k̂k . When performing a line integral in Cartesian coordinates, you
write
r = xı̂ı + y̂ + z k̂k and d r = d xı̂ı + d y̂ + d z k̂k
and can be sure that length scales are properly handled because – as we saw in
Lecture 3 –
p
|d r| = d s = d x 2 + d y 2 + d z 2 .
The reason for using the basis ı̂ı ,̂ , k̂k rather than any other orthonormal basis set is
that ı̂ı represents a direction in which x is increasing while the other two coordinates
remain constant (and likewise for ̂ and k̂k with y and z respectively), simplifying
the representation and resulting mathematics.
Often the symmetry of the problem strongly hints at using another coordinate
system:
• likely to be plane, cylindrical, or spherical polars,
• but can be something more exotic
The general name for any different “u, v , w ” coordinate system is a curvilinear
coordinate system. We will see that the idea hinted at above – of defining a
basis set by considering directions in which only one coordinate is (instantaneously)
increasing – provides the approriate generalisation.
We begin by discussing common special cases: cylindrical polars and spherical
polars, and conclude with a more general formulation.
4.6.1 Cylindrical polar coordinates

As shown in figure 4.4 a point in space P having cartesian coordinates x , y , z can
be expressed in terms of cylindrical polar coordinates, r, φ, z as follows:
r = xı̂ı + y̂ + z k̂k

= r cos φı̂ı + r sin φ̂ + z k̂k
z
Lines of
constant r
Lines of
constant z
^
k
^
j y
^i r
r
θ
Lines of
x constant φ
(a) (b)
Figure 4.4: Cylindrical polars: (a) coordinate definition; (b) “iso” lines in r , φ and z .
Note that, by definition, ∂∂rr represents a direction in which (instantaneously) r is

changing while the other two coordinates stay constant. That is, it is tangent to
∂r ∂r
lines of constant φ and z . Likewise for ∂φ and ∂z , Thus the vectors:
Aside on notation: some texts

use the notation r̂r , φ̂
φ, . . . to rep-
resent the unit vectors that form
∂r
er = = cos φı̂ı + sin φ̂ the local basis set. Though I pre-
∂r
∂r fer the notation used here, where
eφ = = −r sin φı̂ı + r cos φ̂ the basis vectors are written as
∂φ
∂r êe with appropriate subscripts (as
ez = = k̂k used in Riley et al), you should be
∂z
aware of, and comfortable with,
either possibility.
form a basis set in which we may describe infinitessimal vector displacements in

the position of P , d r. It is more usual, however, first to normalise the vectors to
obtain their corresponding unit vectors, êr , êφ , êz . Following the usual rules of
calculus we may write:
∂r ∂r ∂r
dr = dr + dφ + dz
∂r ∂φ ∂z
= d r er + d φeφ + d z ez
= d rêr + r d φêφ + d zêz
Now here is the important thing to note. In cartesian coordinates, a small change
in (eg) x while keeping y and z constant would result in a displacement of

√ p
d s = |d r| = d r.d r = d x 2 + 0 + 0 = d x
But in cylindrical polars, a small change in φ of d φ while keeping r and z constant
results in a displacement of
p
d s = |d r| = r 2 (d φ)2 = r d φ
Thus the size of the (infinitessimal) displacement is dependent on the value of r .
Factors such as this r are known as scale factors or metric coefficients, and we
must be careful to take them into account when, eg, performing line, surface or
volume integrals, as you will below. For cylindrical polars the metric coefficients
are clearly 1, r and 1.
Example: line integral in cylindrical coordinates
Q Evaluate C a · d l , where a = x 3 ̂ − y 3ı̂ı + x 2y k̂k and C is the circle of radius r

H
in the z = 0 plane, centred on the origin.
A Consider figure 4.5. In this case our cylindrical coordinates effectively reduce
to plane polars since the path of integration is a circle in the z = 0 plane, but
let’s persist with the full set of coordinates anyway; the k̂k component of a
will play no role (it is normal to the path of integration and therefore cancels
as seen below).
On the circle of interest
a = r 3 (− sin3 φı̂ı + cos3 φ̂ + cos2 φ sin φk̂k )
and (since dz = dr = 0 on the path)
d r = r d φ êφ
= r d φ(− sin φı̂ı + cos φ̂ )
so that
2π
3π 4
I Z
a · dr = r 4 (sin4 φ + cos4 φ)d φ = r
C 0 2
since
2π 2π
3π
Z Z
4
sin φd φ = cos4 φd φ =
0 0 4
z
From above
y
dR = r dφêφ
dφ
y
φ x
φ r
dR = r dφêφ
Figure 4.5: Line integral example in cylindrical coordinates
Volume integrals in cylindrical polars
In Cartesian coordinates a volume element is given by (see figure 4.6a):
dV = dxdydz
Recall that the volume of a parallelopiped is given by the scalar triple product of
the vectors which define it (see section 2.1.2). Thus the formula above can be
derived (even though it is “obvious”) as:
d V = d xı̂ı .(d y̂ × d z k̂k ) = d x d y d z
since the basis set is orthonormal.

In cylindrical polars a volume element is given by (see figure 4.6b):
d V = d rêr .(r d φêφ × d zêz ) = r d φd r d z
Note also that this volume, because it is a scalar triple product, can be written as
a determinant:
∂y ∂z

êr d r er d r ∂x

∂r ∂r ∂r
∂y ∂z
∂x
d V = êφ r d φ = eφ d φ = ∂φ ∂φ

∂φ d r d φd z
êz d z ez d z ∂x ∂y ∂z
∂z ∂z ∂z
where the equality on the right-hand side follows from the definitions of êr = ∂∂rr =
∂x
ı ∂y  ∂z k
∂r ı̂ + ∂r ̂ + ∂r k̂ , etc. This is the explanation for the “magical” appearance of the
determinant in change-of-variables integration that you encountered in your first
year maths!
y dV = dx dy dz dV = r dr dφdz
dzêz
r dφêφ
dy̂
drêr
x dz k̂k
y
dxı̂ı
x dφ
z φ
r dφ
(a) (b)
Figure 4.6: Volume elements dV in (a) Cartesian coordinates; (b) Cylindrical polar coordinates
Surface integrals in cylindrical polars

Recall from section 4.4 that for a surface element with normal along ı̂ı we have:
d S = d y d zı̂ı
More explicitly this comes from finding normal to the plane that is tangent to the
surface of constant x and from finding the area of an infinitessimal area element
on the plane. In this case the plane is spanned by the vectors ̂ and k̂k and the area
of the element given by (see section 1.3):

d S = d y̂ × d z k̂k

Thus
d S = d y̂ × d z k̂k = ı̂ı d S = d y d zı̂ı
In cylindrical polars, surface area elements (see figure 4.7) are given by:
d S = d rêr × r d φêφ = r d r d φêz (for surfaces of constant z )
d S = r d φêφ × d zêz = r d φd zêr (for surfaces of constant r )
Similarly we can find d S for surfaces of constant φ, though since these aren’t as
common this is left as a (relatively easy) exercise.
dSz = r dr dφêz
r dφêφ drêr
dzêz
dSr = r dφdzêr
dzêz
drêr r dφêφ
dSφ = dr dzêφ
y
Figure 4.7: Surface elements in cylindrical polar coordinates
4.6.2 Spherical polars

Much of the development for spherical polars is similar to that for cylindrical polars.
As shown in figure 4.6.2 a point in space P having cartesian coordinates x , y , z
can be expressed in terms of spherical polar coordinates, r, θ, φ as follows:
r = xı̂ı + y̂ + z k̂k

= r sin θ cos φı̂ı + r sin θ sin φ̂ + r cos θk̂k
The basis set in spherical polars is obtained in an analogous fashion: we find unit
z
z Lines of
constant φ
θ (longitude)
êr Lines of
constant r
P êφ
k̂k êθ
r
̂
y
ı̂ı y
φ
Lines of
x x constant θ
(latitude)
vectors which are in the direction of increase of each coordinate:
∂r
er = = sin θ cos φı̂ı + sin θ sin φ̂ + cos θk̂k = êr
∂r
∂r
eθ = = r cos θ cos φı̂ı + r cos θ sin φ̂ − r sin θk̂k = rêθ
∂θ
∂r
eφ = = −r sin θ sin φı̂ı + r sin θ cos φ̂ = r sin θêφ
∂φ
As with cylindrical polars, it is easily verified that the vectors êr , êθ , êφ form an
orthonormal basis.
A small displacement d r is given by:
∂r ∂r ∂r
dr = dr + dθ + dφ
∂r ∂θ ∂φ
= d r er + d θeθ + d φeφ
= d rêr + r d θêθ + r sin θd φêφ
Thus the metric coefficients are 1, r, r sin θ.
Volume integrals in spherical polars
In spherical polars a volume element is given by (see figure 4.8):
d V = d rêr .(r d θêθ × r sin θd φêφ) = r 2 sin θd r d θd φ
Note again that this volume could be written as a determinant, but this is left as
an exercise.
Surface integrals in spherical polars
The most (the only?) useful surface elements in spherical polars are those tangent
to surfaces of constant r (see figure 4.9). The surface direction (unnormalised) is
given by êθ × êφ = êr and the area of an infinitessimal surface element is given by
|r d θêθ × r sin θd φêφ | = r 2 sin θd θd φ.
Thus a surface element d S in spherical polars is given by
d S = r d θêθ × r sin θd φêφ = r 2 sin θêr

drêr
r sin θ
r sin θdφêφ
dφ dV = r 2 sin θdr dθdφ
r dθêθ
r
θ
dθ
x φ
dφ r sin θdφ
Figure 4.8: Volume element dV in spherical polar coordinates
♣ Example: surface integral in spherical polars
Q Evaluate S a · d S, where a = z 3k̂k

R
and S is the sphere of radius A cen-

tred on the origin.
A On the surface of the sphere:
a = A3 cos 3θk̂k d S = A2 sin θ d θ d φêr
Hence
Z Z 2π Z π
a · dS = A3cos 3θ A2 sin θ [êr · k̂k ] d θd φ
S φ=0 θ=0
Z 2π Z π
5
= A dφ cos 3θ sin θ[cos θ] d θ
0 0
1 π
= 2πA5 − cos5 θ 0
5
5
4πA
=
5
z
r sin θdφêφ
dSr = r 2 sin θdθdφêr
r dθêθ
Figure 4.9: Surface element dS in spherical polar coordinates
4.6.3 General curvilinear coordinates

Cylindrical and spherical polar coordinates are two (useful) examples of general
curvilinear coordinates. In general a point P with Cartesian coordinates x , y , z can
be expressed in terms of the curvilinear coordinates u, v , w where
x = x (u, v , w ), y = y (u, v , w ), z = z (u, v , w )
Thus
r = x (u, v , w )ı̂ı + y (u, v , w )̂ + z (u, v , w )k̂k
and
∂r ∂x ∂y ∂z
= ı̂ı + ̂ + k̂k
∂u ∂u ∂u ∂u
and similarly for partials with respect to v and w , so
∂r ∂r ∂r
dr = du + dv + dw
∂u ∂v ∂w
We now define the local coordinate system as before by considering the directions
in which each coordinate “unilaterally” (and instantaneously) increases:

∂r ∂r
eu = = êu = hu êu
∂u ∂u

∂r ∂r
ev = = êv = hv êv
∂v ∂v

∂r ∂r
ew = = êw = hw êw
∂w ∂w
∂r ∂r ∂r
The metric coefficients are therefore hu = | ∂u |, hv = | ∂v | and hw = | ∂w |.
A volume element is in general given by
d V = hu d uêu .(hv d vêv × hw d wêw )
and simplifies if the coordinate system is orthonormal (since êu .(êv × êw ) = 1) to
d V = hu hv hw d ud v d w
A surface element (normal to constant w , say) is in general
d S = hu d uêu × hv d vêv
and simplifies if the coordinate system is orthogonal to
d S = hu hv d ud vêw
4.6.4 Summary
To summarise:
General curvilinear coordinates
x = x (u, v , w ), y = y (u, v , w ), z = z (u, v , w )

r = x (u, v , w )ı̂ı + y (u, v , w )̂ + z (u, v , w )k̂k

∂r ∂r ∂r
hu = , hv = , hw =
∂u ∂v ∂w
1 ∂r 1 ∂r 1 ∂r
û = êu = , v̂ = êv = , ŵ = êw =
hu ∂u hv ∂v hw ∂w
d r = hu d uû + hv d v v̂ + hw d w ŵ
d V = hu hv hw d ud v d w û.(v̂ × ŵ)
d S = hu hv d ud v û × v̂ (for surface element tangent to constant w )
Plane polar coordinates
x = r cos θ, y = r sin θ
r = r cos θı̂ı + r sin θ̂
hr = 1, hθ = r
êr = cos θı̂ı + sin θ̂ , êθ = − sin θı̂ı + cos θ̂
dr = d rêr + r d θêθ
dS = r d r d θk̂k
Cylindrical polar coordinates
x = r cos φ, y = r sin φ, z =z
r = r cos φı̂ı + r sin φ̂ + z k̂k
hr = 1, hφ = r, hz = 1
êr = cos φı̂ı + sin φ̂ , êφ = − sin φı̂ı + cos φ̂ , êz = k̂k
dr = d rêr + r d φêφ + d zêz
dS = r d r d φk̂k (on the flat ends)
dS = r d φd zêr (on the curved sides)
dV = r d r d φd z
Spherical polar coordinates
x = r sin θ cos φ, y = r sin θ sin φ, z = r cos θ

r = r sin θ cos φı̂ı + r sin θ sin φ̂ + r cos θk̂k
hr = 1, hθ = r, hφ = r sin θ
êr = sin θ cos φı̂ı + sin θ sin φ̂ + cos θk̂k
êθ = cos θ cos φı̂ı + cos θ sin φ̂ + sin θk̂k
êφ = − sin φı̂ı + cos φ̂
dr = d rêr + r d θêθ + r sin θd φêφ
dS = r 2 sin θd r d θd φêr (on a spherical surface)
dV = r 2 sin θd r d θd φ
IDR October 10, 2013

Lecture 5
Vector Operators: Grad, Div and Curl
In the first lecture of the second part of this course we move more to consider
properties of fields. We introduce three field operators which reveal interesting
collective field properties, viz.
• the gradient of a scalar field,
• the divergence of a vector field, and
• the curl of a vector field.
There are two points to get over about each:
• The mechanics of taking the grad, div or curl, for which you will need to brush
up your multivariate calculus.
• The underlying physical meaning — that is, why they are worth bothering
about.
In Lecture 6 we will look at combining these vector operators.
5.1 The gradient of a scalar field

Recall the discussion of temperature distribution throughout a room in the overview,
where we wondered how a scalar would vary as we moved off in an arbitrary direc-
tion. Here we find out how.
If U(x , y , z ) is a scalar field, ie a scalar function of position r = [x , y , z ] in 3
dimensions, then its gradient at any point is defined in Cartesian co-ordinates by
∂U ∂U ∂U
gradU = ı̂ı + ̂ + k̂k .
∂x ∂y ∂z
It is usual to define the vector operator which is called “del” or “nabla”
∂ ∂ ∂
∇ = ı̂ı + ̂ + k̂k .
∂x ∂y ∂z
63
64 LECTURE 5. VECTOR OPERATORS: GRAD, DIV AND CURL
Then
gradU ≡ ∇ U .
Note immediately that ∇ U is a vector field!
Without thinking too carefully about it, we can see that the gradient of a scalar
field tends to point in the direction of greatest change of the field. Later we will
be more precise.
♣ Worked examples of gradient evaluation

1. U = x 2

∂ ∂ ∂
⇒ ∇U = ı̂ı + ̂ + k̂k x 2 = 2xı̂ı .
∂x ∂y ∂z
2. U = r 2
r 2 = x2 + y 2 + z 2
∂ ∂ ∂
⇒ ∇U = ı̂ı + ̂ + k̂k (x 2 + y 2 + z 2 )
∂x ∂y ∂z
= 2xı̂ı + 2y̂ + 2z k̂k = 2 r .
3. U = c · r, where c is constant.

∂ ∂ ∂
⇒ ∇ U = ı̂ı + ̂ + k̂k (c1x + c2y + c3z ) = c1 ı̂ı +c2̂ +c3 k̂k = c .
∂x ∂y ∂z
p
4. U = f (r ), where r = (x 2 + y 2 + z 2 )
U is a function of r alone so df /d r exists. As U = f (x , y , z ) also,
∂f df ∂r ∂f df ∂r ∂f df ∂r
= = = .
∂x d r ∂x ∂y d r ∂y ∂z d r ∂z

∂f ∂f ∂f df ∂r ∂r ∂r
⇒ ∇U = ı̂ı + ̂ + k̂k = ı̂ı + ̂ + k̂k
∂x ∂y ∂z d r ∂x ∂y ∂z
p
But r = x 2 + y 2 + z 2 , so ∂r /∂x = x /r and similarly for y , z .
!
df xı̂ı + y̂ + z k̂k df r
⇒ ∇U = = .
dr r dr r
5.2. THE SIGNIFICANCE OF GRAD 65
gradU r
U(r)
dr
r + dr
U(r + dr)
Figure 5.1: The directional derivative
5.2 The significance of grad

If our current position is r in some scalar field U (Fig. 5.1), and we move an
infinitesimal distance d r, we know that the change in U is
∂U ∂U ∂U
dU = dx + dy + dz .
∂x ∂y ∂z
But we know that d r = (ı̂ı d x + ̂ d y + k̂k d z ) and ∇ U = (ı̂ı ∂U/∂x + ̂ ∂U/∂y +
k̂k ∂U/∂z), so that the change in U is also given by the scalar product
d U = ∇U · d r .
Now divide both sides by d s
dU dr
= ∇U · .
ds ds
But remember that |d r| = d s, so d r/d s is a unit vector in the direction of d r.
This result can be paraphrased as:
• gradU has the property that the rate of change of U wrt distance in a
particular direction (d̂) is the projection of gradU onto that direction
(or the component of gradU in that direction).
The quantity d U/d s is called a directional derivative, but note that in general it
has a different value for each direction, and so has no meaning until you specify
the direction.
We could also say that
• At any point P, gradU points in the direction of greatest change of

U at P, and has magnitude equal to the rate of change of U wrt
distance in that direction.
0.1
0.08
0.06
0.04
0.02
0
4
2 4
2
0
0
−2
−2
−4 −4
Another nice property emerges if we think of a surface of constant U – that is the

locus (x , y , z ) for
U(x , y , z ) = constant .
If we move a tiny amount within that iso-U surface, there is no change in U, so
d U/d s = 0. So for any d r/d s in the surface
dr
∇U · = 0.
ds
But d r/d s is a tangent to the surface, so this result shows that
• gradU is everywhere NORMAL to a surface of constant U.
gradU
Surface of constant U Surface of constant U

These are called Level Surfaces
5.3. THE DIVERGENCE OF A VECTOR FIELD 67
5.3 The divergence of a vector field

The divergence computes a scalar quantity from a vector field by differentiation.
If a(x , y , z ) is a vector function of position in 3 dimensions, that is a = a1ı̂ı + a2 ̂ +
a3k̂k , then its divergence at any point is defined in Cartesian co-ordinates by
∂a1 ∂a2 ∂a3
diva = + +
∂x ∂y ∂z
We can write this in a simplified notation using a scalar product with the ∇ vector
differential operator:

∂ ∂ ∂
diva = ı̂ı + ̂ + k̂k ·a= ∇·a
∂x ∂y ∂z
Notice that the divergence of a vector field is a scalar field.
♣ Examples of divergence evaluation

a diva
1) xı̂ı 1
2) r(= xı̂ı + y̂ + z k̂k ) 3
3) r/r 3 0
4) r c, for c constant (r · c)/r
We work through example 3).
The x component of r/r 3 is x .(x 2 + y 2 + z 2)−3/2, and we need to find ∂/∂x of it.
∂ −3 2
x .(x 2 + y 2 + z 2)−3/2 = 1.(x 2 + y 2 + z 2)−3/2 + x (x + y 2 + z 2)−5/2.2x
∂x 2
−3 2 −2

= r 1 − 3x r .
The terms in y and z are similar, so that
div(r/r 3 ) = r −3 3 − 3(x 2 + y 2 + z 2)r −2 = r −3 (3 − 3)

= 0
5.4 The significance of div

Consider a typical vector field, water flow, and denote it by a(r). This vector has
magnitude equal to the mass of water crossing a unit area perpendicular to the
direction of a per unit time.
Now take an infinitesimal volume element d V and figure out the balance of the
flow of a in and out of d V .
To be specific, consider the volume element d V = d x d y d z in Cartesian co-

ordinates, and think first about the face of area d x d z perpendicular to the y axis
and facing outwards in the negative y direction. (That is, the one with surface
area d S = −d x d z̂ .)
dz
dS = -dxdz j dS = +dxdz j
y
dx
dy
x
Figure 5.2: Elemental volume for calculating divergence.
The component of the vector a normal to this face is a · ̂ = ay , and is pointing
inwards, and so its contribution to the OUTWARD flux from this surface is
a · dS = − ay (y )d z d x ,
where ay (y ) means that ay is a function of y . (By the way, flux here denotes mass
per unit time.)
A similar contribution, but of opposite sign, will arise from the opposite face, but
we must remember that we have moved along y by an amount d y , so that this
OUTWARD amount is

∂ay
ay (y + d y )d z d x = ay + dy dxdz
∂y
The total outward amount from these two faces is
∂ay ∂ay
dydxdz = dV
∂y ∂y
Summing the other faces gives a total outward flux of

∂ax ∂ay ∂az
+ + dV = ∇ · a dV
∂x ∂y ∂z
So we see that
5.5. THE LAPLACIAN: DIV(GRADU) OF A SCALAR FIELD 69
The divergence of a vector field represents the flux generation per unit
volume at each point of the field. (Divergence because it is an efflux not
an influx.)
Interestingly we also saw that the total efflux from the infinitesimal volume was
equal to the flux integrated over the surface of the volume.
(NB: The above does not constitute a rigorous proof of the assertion because we
have not proved that the quantity calculated is independent of the co-ordinate
system used, but it will suffice for our purposes.)
5.5 The Laplacian: div(gradU) of a scalar field

Recall that gradU of any scalar field U is a vector field. Recall also that we
can compute the divergence of any vector field. So we can certainly compute
div(gradU), even if we don’t know what it means yet.
Here is where the ∇ operator starts to be really handy.

∂ ∂ ∂ ∂ ∂ ∂
∇U) = ı̂ı
∇ · (∇ + ̂ + k̂k · ı̂ı + ̂ + k̂k U
∂x ∂y ∂z ∂x ∂y ∂z

∂ ∂ ∂ ∂ ∂ ∂
= ı̂ı + ̂ + k̂k · ı̂ı + ̂ + k̂k U
∂x ∂y ∂z ∂x ∂y ∂z
2
∂2 ∂2

∂
= + + U
∂x 2 ∂y 2 ∂z 2
2
∂ U ∂ 2U ∂ 2U

= + +
∂x 2 ∂y 2 ∂z 2
This last expression occurs frequently in engineering science (you will meet it next
in solving Laplace’s Equation in partial differential equations). For this reason, the
operator ∇2 is called the “Laplacian”
2 2 2

∂ ∂ ∂
∇2 U = + + U
∂x 2 ∂y 2 ∂z 2
Laplace’s equation itself is
∇2 U = 0
♣ Examples of ∇2 U evaluation
U ∇2 U
1) r 2 (= x 2 + y 2 + z 2 ) 6
2) x y 2z 3 2x z 3 + 6x y 2z
3) 1/r 0
Let’s prove example (3) (which is particularly significant – can you guess why?).
1/r = (x 2 + y 2 + z 2)−1/2
∂ ∂ 2 ∂
(x + y 2 + z 2)−1/2 = − x .(x 2 + y 2 + z 2)−3/2
∂x ∂x ∂x
= −(x 2 + y 2 + z 2 )−3/2 + 3x .x .(x 2 + y 2 + z 2)−5/2
= (1/r 3)(−1 + 3x 2/r 2 )
Adding up similar terms for y and z
2 2 2

1 1 (x + y + x )
∇2 = 3 −3 + 3 =0
r r r2
5.6 The curl of a vector field

So far we have seen the operator ∇ applied to a scalar field ∇ U; and dotted with
a vector field ∇ · a.
We are now overwhelmed by an irrestible temptation to
• cross it with a vector field ∇ × a
This gives the curl of a vector field
∇ × a ≡ curl(a)
We can follow the pseudo-determinant recipe for vector products, so that

ı̂ı ̂ k̂k
∂ ∂ ∂
∇ × a = ∂x ∂y ∂z
(remember it this way)
a a a
x y z

∂az ∂ay ∂ax ∂az ∂ay ∂ax
= − ı̂ı + − ̂ + − k̂k
∂y ∂z ∂z ∂x ∂x ∂y
5.7. THE SIGNIFICANCE OF CURL 71
♣ Examples of curl evaluation

a ∇×a
1) −yı̂ı + x̂ 2k̂k

2) x 2y 2k̂k 2x 2yı̂ı − 2x y 2̂
5.7 The significance of curl

Perhaps the first example gives a clue. The field a = −yı̂ı + x̂ is sketched in
Figure 5.3(a). (It is the field you would calculate as the velocity field of an object
rotating with ω = [0, 0, 1].) This field has a curl of 2k̂k , which is in the r-h screw
sense out of the page. You can also see that a field H like this must give a finite
value to the line integral around the complete loop C a · d r.
y
y ax (y+dy)
y+dy
ay (x+dx)
dy
x
ay (x)
dx
y x+dx
x
ax (y)
(a) (b)
Figure 5.3: (a) A rough sketch of the vector field −yı̂ı + x̂ . (b) An element in which to calculate
curl.
In fact curl is closely related to the line integral around a loop.
The
H circulation of a vector a round any closed curve C is defined to be
C a · dr
and the curl of the vector field a represents the vorticity, or circulation
per unit area, of the field.
Our proof uses the small rectangular element d x by d y shown in Figure 5.3(b).
Consider the circulation round the perimeter of a rectangular element.
The fields in the x direction at the bottom and top are
∂ax
ax (y ) and ax (y + d y ) = ax (y ) + dy,
∂y
where ax (y ) denotes ax is a function of y , and the fields in the y direction at the
left and right are
∂ay
ay (x ) and ay (x + d x ) = ay (x ) + dx
∂x
Starting at the bottom and working round in the anticlockwise sense, the four
contributions to the circulation d C are therefore as follows, where the minus signs
take account of the path being opposed to the field:
d C = + [ax (y ) d x ] + [ay (x + d x ) d y ] − [ax (y + d y ) d x ] − [ay (x ) d y ]

∂ay ∂ax
= + [ax (y ) d x ] + ay (x ) + d x d y − ax (y ) + d y d x − [ay (x ) d y ]
∂x ∂y

∂ay ∂ax
= − dx dy
∂x ∂y
∇ × a) · d S
= (∇
where d S = d x d y k̂k .
NB: Again, this is not a completely rigorous proof as we have not shown that the
result is independent of the co-ordinate system used.
5.8 Some definitions involving div, curl and grad

• A vector field with zero divergence is said to be solenoidal.
• A vector field with zero curl is said to be irrotational.
• A scalar field with zero gradient is said to be, er, constant.
Revised Oct 2013

Lecture 6
Vector Operator Identities
In this lecture we look at more complicated identities involving vector operators.

The main thing to appreciate it that the operators behave both as vectors and
as differential operators, so that the usual rules of taking the derivative of, say, a
product must be observed.
There could be a cottage industry inventing vector identities. HLT contains a lot
of them. So why not leave it at that?
First, since grad, div and curl describe key aspects of vectors fields, they arise often
in practice, and so the identities can save you a lot of time and hacking of partial
derivatives, as we will see when we consider Maxwell’s equation as an example
later.
Secondly, they help to identify other practically important vector operators. So,
although this material is a bit dry, the relevance of the identities should become
clear later in other Engineering courses.
6.1 Identity 1: curl grad U = 0

ı
ı̂ 
̂ k
k̂

∇ × ∇U = ∂/∂x ∂/∂y ∂/∂z
∂U/∂x ∂U/∂y ∂U/∂z
2
∂ 2U

∂ U
= ı̂ı − + ̂ () + k̂k ()
∂y ∂z ∂z ∂y
= 0 ,
as ∂ 2/∂y ∂z = ∂ 2/∂z ∂y .
Note that the output is a null vector.
73
74 LECTURE 6. VECTOR OPERATOR IDENTITIES
6.2 Identity 2: div curl a = 0

∂/∂x ∂/∂y ∂/∂z

∇ · ∇ × a = ∂/∂x ∂/∂y ∂/∂z
ax ay az
∂ 2 az ∂ 2 ay ∂ 2 az ∂ 2 ax ∂ 2 ay ∂ 2 ax
= − − + + −
∂x ∂y ∂x ∂z ∂y ∂x ∂y ∂z ∂z ∂x ∂z ∂y
= 0
6.3 Identity 3: div and curl of Ua

Suppose that U(r) is a scalar field and that a(r) is a vector field and we are inter-
ested in the product Ua. This is a vector field, so we can compute its divergence
and curl. For example the density ρ(r) of a fluid is a scalar field, and the instan-
taneous velocity of the fluid v(r) is a vector field, and we are probably interested
in mass flow rates for which we will be interested in ρ(r)v(r).
The divergence (a scalar) of the product Ua is given by:
∇ · (Ua) = U(∇
∇ · a) + (∇
∇U) · a
= Udiva + (gradU) · a
In a similar way, we can take the curl of the vector field Ua, and the result should
be a vector field:
∇ × (Ua) = U∇
∇ × a + (∇
∇U) × a .
6.4 Identity 4: div of a × b

Life quickly gets trickier when vector or scalar products are involved: For example,
it is not that obvious that
div(a × b) = curla · b − a · curlb
To show this, use the determinant:

∂/∂xi ∂/∂xj ∂/∂xk
= ∂ [ay bz − az by ] + ∂ [az bx − ax bz ] + ∂ [ax by − ay bx ]

ax a y a z

bx
∂x ∂y ∂z
by bz
= . . . bash out the products . . .
= curla · b − a · (curl b)
6.5. IDENTITY 5: CURL(A × B) 75
6.5 Identity 5: curl(a × b)

ı̂ı ̂ k̂k

curl(a × b) = ∂/∂x ∂/∂y ∂/∂z

ay bz − az by az bx − ax bz ax by − ay bx
so the ı̂ı component is

∂ ∂
(ax by − ay bx ) − (az bx − ax bz )
∂y ∂z
which can be written as the sum of four terms:

∂by ∂bz ∂ay ∂az ∂ ∂ ∂ ∂
ax + − bx + + by + bz ax − ay + az bx
∂y ∂z ∂y ∂z ∂y ∂z ∂y ∂z
Adding ax (∂bx /∂x) to the first of these, and subtracting it from the last, and
doing the same with bx (∂ax /∂x) to the other two terms, we find that (you should
of course check this):
∇ × (a × b) = (∇
∇ · b)a − (∇
∇ · a)b + [b · ∇]a − [a · ∇]b
where [a · ∇ ] can be regarded as new, and very useful, scalar differential operator.
6.6 Definition of the operator [a · ∇ ]

This is a scalar operator, but it can obviously can be applied to a scalar field,
resulting in a scalar field, or to a vector field resulting in a vector field:

∂ ∂ ∂
[a · ∇ ] ≡ ax + ay + az .
∂x ∂y ∂z
6.7 Identity 6: curl(curla) for you to derive

The following important identity is stated, and left as an exercise:
curl(curla) = graddiva − ∇2a
where
∇2a = ∇2 ax ı̂ı + ∇2ay ̂ + ∇2 az k̂k

♣ Example of Identity 6: electromagnetic waves

Q: James Clerk Maxwell established a set of four vector equations which are fun-
damental to working out how electromagnetic waves propagate. The entire
telecommunications industry is built on these.
divD = ρ
divB = 0
∂
curlE = − B
∂t
∂
curlH = J + D
∂t
In addition, we can assume the following, which should all be familiar to you:
B = µr µ0 H, J = σE, D = ǫr ǫ0E,
where all the scalars are constants.
Now show that in a material with zero free charge density, ρ = 0, and with
zero conductivity, σ = 0, the electric field E must be a solution of the wave
equation
∇2 E = µr µ0 ǫr ǫ0 (∂ 2E/∂t 2 ) .
A: First, a bit of respect. Imagine you are the first to do this — this is a tingle
moment.
divD = div(ǫr ǫ0E) = ǫr ǫ0divE = ρ = 0 ⇒ divE = 0. (a)
divB = div(µr µ0 H) = µr µ0 divH = 0 ⇒ divB = 0 (b)
curlE = −∂B/∂t = −µr µ0 (∂H/∂t) (c)
curlH = J + ∂D/∂t = 0 + ǫr ǫ0 (∂E/∂t) (d )
But we know (or rather you worked out in Identity 6) that curlcurl = graddiv−
∇2 , and using (c)
curlcurlE = graddivE − ∇2E = curl (−µr µ0 (∂H/∂t))

so interchanging the order of partial differentation, and using (a) divE = 0:
∂
−∇2E = −µr µ0 (curlH)
∂t
∂ ∂E
= −µr µ0 ǫr ǫ0
∂t ∂t
2 ∂ 2E
⇒ ∇ E = µr µ0ǫr ǫ0 2
∂t
6.8. GRAD, DIV, CURL AND ∇2 IN CURVILINEAR CO-ORDINATE SYSTEMS 77
This equation is actually three equations, one for each component:
2 ∂ 2 Ex
∇ Ex = µr µ0 ǫr ǫ0 2
∂t
and so on for Ey and Ez .
6.8 Grad, div, curl and ∇2 in curvilinear co-ordinate systems

It is possible to obtain general expressions for grad, div and curl in any orthogonal
curvilinear co-ordinate system by making use of the h factors which were introduced
in Lecture 4.
We recall that the unit vector in the direction of increasing u, with v and w being
kept constant, is
1 ∂r
û =
hu ∂u
where r is the position vector, and

∂r
hu =
∂u
is the metric coefficient. Similar expressions apply for the other co-ordinate direc-
tions. Then
d r = hu ûd u + hv v̂d v + hw ŵd w .
6.9 Grad in curvilinear coordinates

Noting that U = U(r) and U = U(u, v , w ), and using the properties of the gradient
of a scalar field obtained previously
∂U ∂U ∂U
∇U · d r = d U = du + dv + dw
∂u ∂v ∂w
It follows that
∂U ∂U ∂U
∇ U · (hu ûd u + hv v̂d v + hw ŵd w ) = du + dv + dw
∂u ∂v ∂w
The only way this can be satisfied for independent d u, d v , d w is when
1 ∂U 1 ∂U 1 ∂U
∇U = û + v̂ + ŵ
hu ∂u hv ∂v hw ∂w
6.10 Divergence in curvilinear coordinates

Expressions can be obtained for the divergence of a vector field in orthogonal
curvilinear co-ordinates by making use of the flux property.
We consider an element of volume d V . If the curvilinear coordinates are orthogonal
then the little volume is a cuboid (to first order in small quantities) and
d V = hu hv hw d u d v d w .
However, it is not quite a cuboid: the area of two opposite faces will differ as the
scale parameters are functions of u, v and w in general.
h (v) dw
w
h (v+dv) dw
w
y
h (v) du
u h (v+dv) du
u
h v dv
u The scale params are
functions of u,v,w
Figure 6.1: Elemental volume for calculating divergence in orthogonal curvilinear coordinates
So the net efflux from the two faces in the v̂ direction shown in Figure 6.1 is

∂av ∂hu ∂hw
= av + d v hu + d v hw + d v d ud w − av hu hw d ud w
∂v ∂v ∂v
∂(av hu hw )
= d ud v d w
∂v
which is easily shown by multiplying the first line out and dropping second order
terms (i.e. (d v )2).
By definition div is the net efflux per unit volume, so summing up the other faces:

∂(au hv hw ) ∂(av hu hw ) ∂(aw hu hv )
diva d V = + + d ud v d w
∂u ∂v ∂w

∂(au hv hw ) ∂(av hu hw ) ∂(aw hu hv )
⇒ diva hu hv hw d ud v d w = + + d ud v d w
∂u ∂v ∂w
6.11. CURL IN CURVILINEAR COORDINATES 79
So, finally,

1 ∂(au hv hw ) ∂(av hu hw ) ∂(aw hu hv )
diva = + +
hu hv hw ∂u ∂v ∂w
6.11 Curl in curvilinear coordinates

Recall from Lecture 5 that we computed the z component of curl as the circulation
per unit area from

∂ay ∂ax
dC = − dx dy
∂x ∂y
By analogy with our derivation of divergence, you will realize that for an orthogonal
curvilinear coordinate system we can write the area as hu hv d ud w . But the opposite
sides are no longer quite of the same length. The lower of the pair in Figure 6.2
is length hu (v )d u, but the upper is of length hu (v + d v )d u
y
au (v+dv)
v+dv
h u(v+dv) du
dv
hu (v) du
y u+du
u
au (v)
Figure 6.2: Elemental loop for calculating curl in orthogonal curvilinear coordinates
Summing this pair gives a contribution to the circulation

∂(hu au )
au (v )hu (v )d u − au (v + d v )hu (v + d v )d u = − dvdu
∂v
and together with the other pair:

∂(hu au ) ∂(hv av )
dC = − + d ud v
∂v ∂u
So the circulation per unit area is

dC 1 ∂(hv av ) ∂(hu au )
= −
hu hv d ud v hu hv ∂u ∂v
and hence curl is

1 ∂(hw aw ) ∂(hv av )
curla(u, v , w ) = − û +
hv hw ∂v ∂w

1 ∂(hu au ) ∂(hw aw )
− v̂ +
hw hu ∂w ∂u

1 ∂(hv av ) ∂(hu au )
− ŵ
hu hv ∂u ∂v
You should check that this can be written as

Curl in curvilinear coords:

h û hv v̂ hw ŵ
1 u∂

∂ ∂

curla(u, v , w ) = ∂u ∂v ∂w

hu hv hw

hu au hv av hw aw
6.12 The Laplacian in curvilinear coordinates

Substitution of the components of gradU into the expression for diva immediately
(!*?) gives the following expression for the Laplacian in general orthogonal co-
ordinates:

2 1 ∂ hv hw ∂U ∂ hw hu ∂U ∂ hu hv ∂U
∇U= + + .
hu hv hw ∂u hu ∂u ∂v hv ∂v ∂w hw ∂w
6.13 Grad Div, Curl, ∇2 in cylindrical polars

Here (u, v , w ) → (r, φ, z ). The position vector is r = r cos φı̂ı + r sin φ̂ + z k̂k , and
hr = |∂r/∂r |, etc.
p
⇒ hr = (cos2 φ + si n2φ) = 1,
q
hφ = (r 2 sin2 φ + r 2 cos2 φ) = r,
hz = 1
6.14. GRAD DIV, CURL, ∇2 IN SPHERICAL POLARS 81
∂U 1 ∂U ∂U
⇒ gradU = êr + êφ + k̂k
∂r r ∂φ ∂z
1 ∂(r ar ) ∂aφ ∂az
diva = + +
r ∂r ∂φ ∂z

1 ∂az ∂aφ ∂ar ∂az 1 ∂(r aφ) ∂ar
curla = − êr + − êφ + − k̂k
r ∂φ ∂z ∂z ∂r r ∂r ∂φ
∇2U = Tutorial Exercise
6.14 Grad Div, Curl, ∇2 in spherical polars

Here (u, v , w ) → (r, θ, φ). The position vector is r = r sin θ cos φı̂ı + r sin θ sin φ̂ +
r cos θk̂k .
q
⇒ hr = (sin2 θ(cos2 φ + sin2 φ) + cos2 θ) = 1
q
hθ = (r 2 cos2 θ(cos2 φ + sin2 φ) + r 2 sin2 θ) = r
q
hφ = (r 2 sin2 θ(sin2 φ + cos2 φ) = r sin θ
∂U 1 ∂U 1 ∂U
⇒ gradU = êr + êθ + êφ
∂r r ∂θ r sin θ ∂φ
1 ∂(r 2ar ) 1 ∂(aθ sin θ) 1 ∂aφ
diva = 2 + +
r ∂r r sin θ ∂θ r sin θ ∂φ
êr ∂ ∂ êθ ∂ ∂
curla = (aφ sin θ) − (aθ ) + (ar ) − (aφr sin θ) +
r sin θ ∂θ ∂φ r sin θ ∂φ ∂r

êφ ∂ ∂
(aθ r ) − (ar )
r ∂r ∂θ
2
∇ U = Tutorial Exercise
♣ Examples
Q1 Find curla in (i) Cartesians and (ii) Spherical polars when a = x (xı̂ı + y̂ + z k̂k ).
A1 (i) In Cartesians

ı̂ı 
̂ k
k̂

curla = ∂/∂x ∂/∂y ∂/∂z = −z̂ + y k̂k .
x2 xy xz
(ii) In spherical polars, x = r sin θ cos φ and r = (xı̂ı + y̂ + z k̂k ). So
a = r 2 sin θ cos φêr

⇒ ar = r 2 sin θ cos φ; aθ = 0; aφ = 0 .
Hence as

êr ∂ ∂ êθ ∂ ∂ êφ ∂ ∂
curla = (aφ sin θ) − (aθ ) + (ar ) − (aφ r sin θ) + (aθ r ) − (ar )
r sin θ ∂θ ∂φ r sin θ ∂φ ∂r r ∂r ∂θ

êθ ∂ 2 êφ ∂ 2
curla = (r sin θ cos φ) + − (r sin θ cos φ)
r sin θ ∂φ r ∂θ
êθ êφ
(−r 2 sin θ sin φ) + −r 2 cos θ cos φ)

=
r sin θ r
= êθ (−r sin φ) + êφ (−r cos θ cos φ)
Checking: these two results should be the same, but to check we need ex-
pressions for êθ , êφ in terms of ı̂ı etc.
Remember that we can work out the unit vectors êr and so on in terms of ı̂ı
etc using
1 ∂r 1 ∂r 1 ∂r
êr = ; êθ = ; êφ = where r = xı̂ı +y̂ +z k̂k .
h1 d r h2 d θ h3 d φ
Grinding through we find
ı̂ı ı̂ı
      
êr sin θ cos φ sin θ sin φ cos θ
 êθ  =  cos θ cos φ cos θ sin φ − sin θ   ̂  = [R]  ̂ 
êφ − sin φ cos φ 0 k̂k k̂k
Don’t be shocked to see a rotation matrix [R]: we are after all rotating one
right-handed orthogonal coord system into another.
So the result in spherical polars is
curla = (cos θ cos φı̂ı + cos θ sin φ̂ − sin θk̂k )(−r sin φ) + (− sin φı̂ı + cos φ̂ )(−r cos θ cos φ
= −r cos θ̂ + r sin θ sin φk̂k
= −z̂ + y k̂k
which is exactly the result in Cartesians.
Q2 Find the divergence of the vector field a = r c where c is a constant vector

(i) using Cartesian coordinates and (ii) using Spherical Polar coordinates.
6.14. GRAD DIV, CURL, ∇2 IN SPHERICAL POLARS 83
A2 (i) Using Cartesian coords:

∂ 2
diva = (x + y 2 + z 2 )1/2cx + . . .
∂x
= x .(x 2 + y 2 + z 2 )−1/2cx + . . .
1
= r·c .
r
(ii) Using Spherical polars

a = ar êr + aθ êθ + aφ êφ
and our first task is to find ar and so on. We can’t do this by inspection, and
finding their values requires more work than you might think! Recall
ı̂ı ı̂ı
      
êr sin θ cos φ sin θ sin φ cos θ
 êθ  =  cos θ cos φ cos θ sin φ − sin θ   ̂  = [R]  ̂ 
êφ − sin φ cos φ 0 k̂k k̂k
Now the point is the same point in space whatever the coordinate system, so
ar êr + aθ êθ + aφêφ = ax ı̂ı + ay ̂ + az k̂k
and using the inner product
⊤  ⊤  
ı̂ı
  
ar êr ax
 aθ   êθ  =  ay   ̂ 
aφ êφ az k̂k
⊤ ⊤  
ı̂ı ı̂ı
   
ar ax
 aθ  [R]  ̂  =  ay   ̂ 
aφ k̂k az k̂k
 ⊤  ⊤
ar ax
⇒  aθ  [R] =  ay 
aφ az
 ⊤  ⊤
ar ax
⇒  aθ  =  ay  [R]⊤
aφ az
   
ar ax
⇒  aθ  = [R]  ay 
aφ az
For our particular problem, ax = r cx , etc, where cx is a constant, so now we

can write down
ar = r (sin θ cos φcx + sin θ sin φcy + cos θcz )

aθ = r (cos θ cos φcx + cos θ sin φcy − sin θcz )
aφ = r (− sin φcx + cos φcy )
Now all we need to do is to bash out

1 ∂(r 2 ar ) 1 ∂(aθ sin θ) 1 ∂aφ
diva = + +
r 2 ∂r r sin θ ∂θ r sin θ ∂φ
In glorious detail this is
diva = 3 (sin θ cos φcx + sin θ sin φcy + cos θcz ) +

1
cos2 θ − sin2 θ)(cos φcx + sin φcy ) − 2 sin θ cos θcz +

sin θ
1
(− cos φcx − sin φcy )
sin θ
A bit more bashing and you’ll find
diva = sin θ cos φcx + sin θ sin φcy + cos θcz

= êr · c
This is EXACTLY what you worked out before of course.
Take home messages from these examples:
• Just as physical vectors are independent of their coordinate systems, so are

differential operators.
• Don’t forget about the vector geometry you did in the 1st year. Rotation
matrices are useful!
• Spherical polars were NOT a good coordinate system in which to think about
this problem. Let the symmetry guide you.
Revised Oct 2013

Lecture 7
Gauss’ and Stokes’ Theorems
This section finally begins to deliver on why we introduced div grad and curl. Two
theorems, both of them over two hundred years old, are explained:
• Gauss’ Theorem enables an integral taken over a volume to be replaced by

one taken over the surface bounding that volume, and vice versa. Why would
we want to do that? Computational efficiency and/or numerical accuracy!
• Stokes’ Law enables an integral taken around a closed curve to be replaced
by one taken over any surface bounded by that curve.
7.1 Gauss’ Theorem

Suppose that a(r) is a vector field and we want to compute the total flux of the
field across the surface S that bounds a volume V . That is, we are interested in
calculating:
Z
a · dS
S
dS
dS
dS dS
dS
Figure 7.1: The surface element dS must stick out of the surface.
85
86 LECTURE 7. GAUSS’ AND STOKES’ THEOREMS
where recall that d S is normal to the locally planar surface element and must
everywhere point out of the volume as shown in Figure 7.1.
Gauss’ Theorem tells us that we can do this by considering the total flux generated
inside the volume V :
Gauss’ Theorem
Z Z
a · dS = div a d V
S V
obtained by integrating the divergence over the entire volume.
7.2 Informal proof

An non-rigorous proof can be realized by recalling that we defined div by considering
the efflux d E from the surfaces of an infinitesimal volume element
dE = a · dS
and defining it as
div a d V = dE = a · dS .
If we sum over the volume elements, this results in a sum over the surface elements.
But if two elemental surface touch, their d S vectors are in opposing direction and
cancel as shown in Figure 7.2. Thus the sum over surface elements gives the
overall bounding surface.
Figure 7.2: When two elements touch, the dS vectors at the common surface cancel out. One
can imagine building the entire volume up from the infinitesimal units.
7.2. INFORMAL PROOF 87
♣ Example of Gauss’ Theorem

This is a typical example, in which the surface integral is rather tedious, whereas
the volume integral is straightforward.
Q Derive S a · d S where a = z 3k̂k and S is the surface of a sphere of radius R

R
centred on the origin:
1. directly;
2. by applying Gauss’ Theorem
3
z k
2
R sin θ d θ d φ r
dz
z
R
R
Figure 7.3:
A (1) On the surface of the sphere, a = R3 cos3 θk̂k and d S = R2 sin θd θd φr̂.
Everywhere r̂ · k̂k = cos θ.
Z Z 2π Z π
⇒ a · dS = R3 cos3 θ . R2 sin θd θd φêr · k̂k
S φ=0 θ=0
Z 2π Z π
= R3 cos3 θ . R2 sin θd θd φ . cos θ
φ=0 θ=0
Z π
= 2πR5 cos4 θ sin θd θ
0
2πR5 5
π 4πR5
= − cos θ 0 =
5 5
(2) To apply Gauss’ Theorem, we need to figure out div a and decide how to
compute the volume integral. The first is easy:
diva = 3z 2
For the second, because diva involves just z , we can divide the sphere into
discs of constant z and thickness d z , as shown in Fig. 7.3. Then
d V = π(R2 − z 2)d z
and
Z Z R
div ad V = 3π z 2(R2 − z 2)d z
V −R
2 3
R
z5

R z
= 3π −
3 5 −R
4πR5
=
5
7.3 Surface versus volume integrals

At first sight, it might seem that with a computer performing surface integrals
might be better than a volume integral, perhaps because there are, somehow,
“fewer elements”. However, this is not the case. Imagine doing a surface integral
over a wrinkly surface, say that of the moon. All the elements involved in the
integration are “difficult” and must be modelled correctly. With a volume integral,
most of the elements are not at the surface, and so the bulk of the integral is
done without accurate modelling. The computation is easier, faster, and better
conditioned numerically.
7.4 Extension to Gauss’ Theorem

Suppose the vector field a(r) is of the form a = U(r)c, where U(r) as scalar field
and c is a constant vector. Then, as we showed in the previous lecture,
div a = gradU · c + Udiv c

= gradU · c
since divc = 0 because c is constant.

Gauss’ Theorem becomes
Z Z
Uc · d S = grad U · cd V
S V
or, alternatively, taking the constant c out of the integrals

Z Z
c· Ud S = c · grad Ud V
S V
7.4. EXTENSION TO GAUSS’ THEOREM 89
This is still a scalar equation but we now note that the vector c is arbitrary so
that the result must be true for any vector c. This can be true only if the vector
equation
Z Z
Ud S = grad Ud V
S V
is satisfied.
If you think this is fishy, just write c =R ı̂ı , then c = ̂ , and c = k̂k in turn, and you
must obtain the three components of S Ud S in turn.
Further “extensions” can be obtained of course. For example one might be able
to write the vector field of interest as
a(r) = b(r) × c
where c is a constant vector.
♣ Example of extension to Gauss’ Theorem
Q U = x 2 + y 2 + z 2 is a scalar field,
and volume V is the cylinder x 2 +
y 2 ≤ a2, 0 ≤ z ≤ h. Compute the
surface integral
Z z dS
Ud S
S
over the surface of the cylinder. z=h

A It is immediately clear from sym-
metry that there is no contribution
from the curved surface of the cylin- dS
der since for every vector surface el- z=0
ement there exists an equal and op-
posite element with the same value
of U. We therefore need consider
only the top and bottom faces.
Top face:
U = x 2 + y 2 + z 2 = r 2 + h2 and d S = r d r d φk̂k
so
a 2π a
1 1
Z Z Z
Ud S = (h2 + r 2 )2πr d r d φk̂k = k̂k π h2 r 2 + r 4 = π[h2 a2 + a4]k̂k
r =0 φ=0 2 0 2
Bottom face:
U = r 2 and d S = −r d r d φk̂k
4
The contribution from this face is thus − πa2 k̂k , and the total integral is πh2 a2 k̂k .
On the other hand, using Gauss’ Theorem we have to compute
Z
grad Ud V
V
In this case, grad U = 2r,

Z
2 (xı̂ı + y k̂k + z k̂k )r d r d z d φ
V
The integrations over x and y are zero by symmetry, so that the only remaining
part is
Z h Z a Z 2π
2 zdz r dr d φk̂k = πa2 h2 k̂k
z=0 r =0 φ=0
7.5 Stokes’ Theorem

Stokes’ Theorem relates a line integral around a closed path to a surface integral
over what is called a capping surface of the path.
Stokes’ Theorem states:

I Z
a · dl = curl a · d S
C S
where S is any surface capping the curve C.
Why have we used d l rather than d r, where r is the position vector?

There is no good reason for this, as d l = d r. It just seems to be common usage
in line integrals!
7.6. INFORMAL PROOF 91
7.6 Informal proof

You will recall that in Lecture 5 that we defined curl as the circulation per unit
area, and showed that
X
∇ × a) · d S .
a · d l = d C = (∇
around elemental loop
Now if we add these little loops together, the internal line sections cancel out
because the d l’s are in opposite direction but the field a is not. This gives the
larger surface and the larger bounding contour as shown in Fig. 7.4.
y
ax (y+dy)
y+dy
ay (x+dx)
dy
ay (x)
dx
y x+dx
x
ax (y)
Figure 7.4: An example of an elementary loop, and how they combine together.
For a given contour, the capping surface can be ANY surface bound by
the contour. The only requirement is that the surface element vectors point in
the “general direction” of a right-handed screw with respect to the sense of the
contour integral. See Fig. 7.5.
Back
Front
Back
Front
Figure 7.5: For a given contour, the bounding surface can be any shape. dS’s must have a positive
component in the sense of a r-h screw wrt the contour sense.
♣ Example of Stokes’ Theorem

In practice, (and especially in exam questions!) the bounding contour is often
planar, and the capping surface flat or hemispherical or cylindrical.
Q Vector field a = x 3̂ − y 3ı̂ı and C is the circle of radius R centred on the origin.
Derive
I
a · dl
C
directly and (ii) using Stokes’ theorem where the surface is the planar surface
bounded by the contour.
A(i) Directly. On the circle of radius R
a = R3 (− sin3 θı̂ı + cos3 θ̂ )

and
d l = Rd θ(− sin θı̂ı + cos θ̂ )
so that:
2π
3π 4
I Z
a · dl = R4 (sin4 θ + cos4 θ)d θ = R ,
C 0 2
since
2π 2π
3π
Z Z
4
sin θd θ = cos4 θd θ =
0 0 4
A(ii) Using Stokes’ theorem ...

ı̂ı ̂ k̂k
∂ ∂ ∂ = 3(x 2 + y 2 )k̂k = 3r 2 k̂k

curl a = ∂x ∂y ∂z
−y 3 x 3 0
We choose area elements to be circular strips of radius r thickness d r . Then

Z R
3π 4
Z
d S = 2πr d r k̂k and curl a · d S = 6π r 3d r = R
S 0 2
7.7. AN EXTENSION TO STOKES’ THEOREM 93
7.7 An Extension to Stokes’ Theorem

Just as we considered one extension to Gauss’ theorem (not really an extension,
more of a re-expression), so we will try something similar with Stoke’s Theorem.
Again let a(r) = U(r)c, where c is a constant vector. Then
curl a = Ucurl c + grad U × c)
Again, curl c is zero. Stokes’ Theorem becomes in this case:
I Z Z
U(c · d l) = (grad U × c · d S = c · (d S × grad U)
C S S
or, rearranging the triple scalar products and taking the constant c out of the
integrals gives
I Z
c · Ud l = −c · gradU × d S .
C S
But c is arbitrary and so

I Z
Ud l = − grad U × d S
C S
7.8 ♣ Example of extension to Stokes’ Theorem

H y
Q Derive C Ud r (i) directly and (ii) using Stokes’, d r = dρ
where U = x 2 + y 2 + z 2 and the line integral is taken
around C the circle (x − a)2 + y 2 = a2 and z = 0.
r ρ
α x
Note that, for no special reason, we have used d r a
here not d l.
A(i) First some preamble.

If the circle were centred at the origin, we would write d r = ad θêθ =
ad θ(− sin θı̂ı +cos θ̂ ). For such a circle the magnitude r = |r| = a, a constant
and so d r = 0.
However, in this example d r is not always in the direction of êθ , and d r 6= 0.
Could you write down d r? If not, revise Lecture 3, where we saw that in plane
polars x = r cos θ, y = r sin θ and the general expression is
d r = d xı̂ı + d y̂ = (cos θd r − r sin θd θ)ı̂ı + (sin θd r + r cos θd θ)̂
To avoid having to find an expression for r in terms of θ, we will perform a

coordinate transformation by writing r = [a, 0]⊤ + ρ . So, x = (a + ρ cos α)
and y = ρ sin α, and on the circle itself where ρ = a
r = a(1 + cos α)ı̂ı + a sin α̂ ,
d r = ad α(− sin αı̂ı + cos α̂ ) ,
and, as z = 0 on the circle,
U = a2 (1 + cos α)2 + a2 sin2 α = 2a2(1 + cos α) .
The line integral becomes
I Z 2π
3
Ud r = 2a (1 + cos α)(− sin αı̂ı + cos α̂ )d α = 2πa3 ̂
α=0
A(ii) Now using Stokes’ ...

For a planar surface covering the disc, the surface element can be written
using the new parametrization as
d S = ρ d ρ d αk̂k
Remember that U = x 2 + y 2 + z 2 = r 2, and as z = 0 in the plane
grad U = 2(xı̂ı + y̂ + z k̂k ) = 2(a + ρ cos α)ı̂ı + 2ρ sin α̂ .

Be careful to note that x , y are specified for any point on the disc, not on its
circular boundary!
So

ı
ı̂ 
̂ k
k̂

d S×gradU = 2ρ d ρ d α 0 0 1 = 2ρ[−ρ sin αı̂ı +(a+ρ cos α)̂ ] d ρ d α
(a + ρ cos α) ρ sin α 0
R 2π R 2π
Both 0 sin αd α = 0 and 0 cos αd α = 0, so we are left with
Z Z a Z 2π
d S × gradU = 2ρa̂ d ρ d α = 2πa3 ̂
S ρ=0 α=0
Revised Oct 2013

Lecture 8
Engineering applications
In Lecture 6 we saw one classic example of the application of vector calculus to

Maxwell’s equation.
In this lecture we explore a few more examples from fluid mechanics and heat
transfer. As with Maxwell’s eqations, the examples show how vector calculus
provides a powerful way of representing underlying physics.
The power come from the fact that div, grad and curl have a significance or
meaning which is more immediate than a collection of partial derivatives. Vector
calculus will, with practice, become a convenient shorthand for you.
• Electricity – Ampère’s Law

• Fluid Mechanics - The Continuity Equation
• Thermo: The Heat Conduction Equation
• Mechanics/Electrostatics - Conservative fields
• The Inverse Square Law of force
• Gravitational field due to distributed mass
• Gravitational field inside body
• Pressure forces in non-uniform flows
95
96 LECTURE 8. ENGINEERING APPLICATIONS
8.1 Electricity – Ampère’s Law

If the frequency is low, the displacement current in Maxwell’s equation curlH =
J + ∂D/∂t is negligible, and we find
curlH = J
Hence
Z Z
curlH · d S = J · dS
S S
or
I Z
H · dl = J · dS
S
R
where S J · d S is total current through the surface.
Now consider the H around a straight wire carrying current I. Symmetry tells us
the H is in the êθ direction, in a rhs screw sense with respect to the current. (You
might check this against Biot-Savart’s law.)
Suppose we asked what is the magnitude of H?
I
Current
r H
1111
0000
Back r H
H a r
C Top view
H Front In wire Outside wire
for r<a
Inside the wire, the bounding contour only encloses a fraction (πr 2 )/(πa2 ) of the
current, and so
Z
H2πr = J · d S = I(r 2 /A2)
⇒ H = Ir /2πA2
whereas outside we enclose all the current, and so

Z
H2πr = J · dS = I
⇒ H = I/2πr
A plot is shown in the Figure.

8.2. FLUID MECHANICS - THE CONTINUITY EQUATION 97
8.2 Fluid Mechanics - The Continuity Equation

The Continuity Equation expresses the condition of conservation of mass in a
fluid flow. The continuity principle applied to any volume (called a control volume)
may be expressed in words as follows:
“The net rate of mass flow of fluid out of the control volume must equal
the rate of decrease of the mass of fluid within the control volume”
dS q
Control Volume V
Figure 8.1:
To express the above as a mathematical equation, we denote the velocity of the

fluid at each point of the flow by q(r) (a vector field) and the density by ρ(r) (a
scalar field). The element of rate-of-volume-loss through surface d S is d V̇ = q·d S,
so the rate of mass loss is
d Ṁ = ρq · d S,
so that the total rate of mass loss from the volume is
∂
Z Z
− ρ(r)d V = ρq · d S.
∂t V S
Assuming that the volume of interest is fixed, this is the same as
∂ρ
Z Z
− dV = ρq · d S .
V ∂t S
Now we use Gauss’ Theorem to transform the RHS into a volume integral
∂ρ
Z Z
− dV = div (ρq)d V .
V ∂t V
The two volume integrals can be equal for any control volume V only if the two
integrands are equal at each point of the flow. This leads to the mathematical
formulation of
The Continuity Equation:

∂ρ
div (ρq) = −
∂t
Notice that if the density doesn’t vary with time, div (ρq) = 0, and if the density
doesn’t vary with position then
The Continuity Equation for uniform, time-invariant density:
div (q) = 0 .
In this last case, we can say that the flow q is solenoidal.
8.3 Thermodynamics - The Heat Conduction Equation

Flow of heat is very similar to flow of fluid, and heat flow satisfies a similar con-
tinuity equation. The flow is characterized by the heat current density q(r) (heat
flow per unit area and time), sometimes misleadingly called heat flux.
Assuming that there is no mass flow across the boundary of the control volume and
no source of heat inside it, the rate of flow of heat out of the control volume by
conduction must equal the rate of decrease of internal energy (constant volume)
or enthalpy (constant pressure) within it. This leads to the equation
∂T
div q = −ρc ,
∂t
where ρ is the density of the conducting medium, c its specific heat (both are
assumed constant) and T is the temperature.
In order to solve for the temperature field another equation is required, linking q
to the temperature gradient. This is
q = −κgrad T,
where κ is the thermal conductivity of the medium. Combining the two equations
gives the heat conduction equation:
∂T
−div q = κdiv grad T = κ∇2 T = ρc
∂t
where it has been assumed that κ is a constant. In steady flow the temperature
field satisfies Laplace’s Equation ∇2T = 0.
8.4. MECHANICS - CONSERVATIVE FIELDS OF FORCE 99
8.4 Mechanics - Conservative fields of force

A conservative field of force is one for which the work done
Z B
F · d r,
A
moving from A to B is indep. of path taken. As we saw in Lecture 4, conservative

fields must satisfy the condition
I
F · d r = 0,
C
Stokes’ tells us that this is

Z
curl F · d S = 0,
S
where S is any surface bounded by C.

But if true for any C containing A and B, it must be that
curl F = 0
Conservative fields are irrotational
All radial fields are irrotational
One way (actually the only way) of satisfying this condition is for
F=∇ U
The scalar field U(r) is the Potential Function
8.5 The Inverse Square Law of force

Radial forces are found in electrostatics and gravitation — so they are certainly
irrotational and conservative.
But in nature these radial forces are also inverse square laws. One reason why this
may be so is that it turns out to be the only central force field which is solenoidal,
i.e. has zero divergence.
If F = f (r )r,
div F = 3f (r ) + r f ′ (r ).
For div F = 0 we conclude

df
r + 3f = 0
dr
or
df dr
+3 = 0.
f r
Integrating with respect to r gives f r 3 = const = A, so that
Ar A
F= 3
, |F| = 2 .
r r
The condition of zero divergence of the inverse square force field applies everywhere
except at r = 0, where the divergence is infinite.
To show this, calculate the outward normal flux out of a sphere of radius R centered
on the origin when F = Fr̂. This is
Z Z Z
F · dS = F r̂ · d S = F d = F 4πR2 = 4πA = Constant.
S Sphere Sphere
Gauss tells us that this flux must be equal to

Z Z R
div Fd V = div F4πr 2 d r
V 0
where we have done the volume integral as a summation over thin shells of surface
area 4πr 2 and thickness d r .
But for all finite r , divF = 0, so divF must be infinite at the origin.
The flux integral is thus
• zero — for any volume which does not contain the origin
• 4πA for any volume which does contain it.
8.6. GRAVITATIONAL FIELD DUE TO DISTRIBUTED MASS: POISSON’S EQUATION 101
8.6 Gravitational field due to distributed mass: Poisson’s Equa-

tion
If one tried the same approach as §8.4 for the gravitational field, A = Gm, where
m is the mass at the origin and G the universal gravitational constant, one would
run into the problem that there is no such thing as point mass.
We can make progress though by considering distributed mass.
The mass contained in each small volume element d V is ρd V and this will make a
contribution −4πρGd V to the flux integral from the control volume. Mass outside
the control volume makes no contribution, so that we obtain the equation
Z Z
F · d S = −4πG ρd V.
S V
Transforming the left hand integral by Gauss’ Theorem gives

Z Z
div Fd V = −4πG ρd V
V V
which, since it is true for any V , implies that

−div F = 4πρG.
Since the gravitational field is also conservative (i.e. irrotational) it must have
an associated potential function U, so that F = grad U. It follows that the
gravitational potential U satisfies
Poisson’s Equation
∇2U = 4πρG .
Using the integral form of Poisson’s equation, it is possible to calculate the gravi-
tational field inside a spherical body whose density is a function of radius only. We
have
Z R
4πR2 F = 4πG 4πr 2 ρd r,
0
where F = |F|, or
Z R
G MG
|F | = 2 4πr 2 ρd r = 2 ,
R 0 R
where M is the total mass inside radius R. For the case of uniform density, this is
equal to M = 34 πρR3 and |F | = 43 πρGR.
8.7 Pressure forces in non-uniform flows

When a body is immersed in a flow it experiences a net pressure force
Z
Fp = − pd S,
S
where S is the surface of the body. If the pressure p is non-uniform, this integral
is not zero. The integral can be transformed using Gauss’ Theorem to give the
alternative expression
Z
Fp = − grad p d V,
V
where V is the volume of the body. In the simple hydrostatic case p + ρgz =
constant, so that
grad p = −ρgk
and the net pressure force is simply
Z
Fp = g k̂ ρd V
V
which, in agreement with Archimedes’ principle, is equal to the weight of fluid

displaced.
Revised Oct 2013

Course PDF

Uploaded by

Copyright:

Available Formats

Course PDF

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Course PDF

Uploaded by

Copyright:

Available Formats

What is the course about?

What is the course about?

What are some of the main topics covered in vector calculus?

What are some of the main topics covered in vector calculus?

1

2A1 Vector Algebra and Calculus

8 Lectures MT 2013 Stephen Roberts

but how would you calculate that? Is ∂T /∂x a vector or a scalar?

Course WWW Pages

• Free vectors: In many situtations only the magnitude and direction of a

1.1.1 Vector elements or components in a coordinate frame

Figure 1.2: Vector components.

In a Cartesian coordinate frame we write

we could also write

1.1.2 Vector equality

1.1.3 Vector magnitude and unit vectors

1.1.4 Vector Addition and subtraction

Addition follows the parallelogram construction of Figure 1.3(a). Subtraction (a −

1.1.5 Multiplication of a vector by a scalar. (NOT the scalar product!)

Note that F(R) is a vector field.

1.2 Scalar, dot, or inner product

1.2.1 Geometrical interpretation of scalar product

Figure 1.4: (a) Cosine rule. (b) Projection of b onto a.

1.2.2 Projection of one vector onto the other

1.2.3 A scalar product is an “inner product”

A2. Well, you could plod through using

A3. Best to refer everything to a coordi-

sbefore = (s.û)û + (s.v̂)v̂ ^v

After the impact, the component of s

safter = (s.û)û − (s.v̂)v̂

1.2.4 Direction cosines use projection

a · ı̂ı a · ̂ a · k̂k

1.3 Vector or cross product

a × b = (a2b3 − a3b2 )ı̂ı + (a3b1 − a1b3 )̂ + (a1b2 − a2 b1)k̂k .

It is MUCH more easily remembered in terms of the (pseudo-)determinant

from which it follows that

which is again independent of the co-ordinate system used. This is left as an

1.3.1 Geometrical interpretation of vector product

Revised Oct 2013

Multiple Products. Geometry using Vectors

2.1 Triple and multiple products

2.1.1 Scalar triple product

(i) a · (b × c) = b · (c × a) = c · (a × b) (Called cyclic permutation.)

(ii) a · (b × c) = −b · (a × c) and so on. (Called anti-cyclic permutation.)

2.1.2 Geometrical interpretation of scalar triple product

2.1.3 Linearly dependent vectors

2.1.4 Vector triple product

(a × (b × c))1 = a2 (b × c)3 − a3(b × c)2

Figure 2.2: Vector triple product.

2.1.5 Projection using vector triple product

So v − (v · k̂k )k̂k = k̂k × (v × k̂k ).

2.1.6 Other repeated products

2.2 Geometry using vectors: lines, planes

2.2.2 The shortest distance from a point to a line

2.2.3 The shortest distance between two straight lines

NB that these are

2.2.5 The shortest distance from a point to a plane

2.3 Solution of vector equations

Step (iii): Bung this expression for x into the equation!