MA20010: Vector Calculus and Partial Differential Equations
MA20010: Vector Calculus and Partial Differential Equations
MA20010: Vector Calculus and Partial Differential Equations
Differential Equations
Robert Scheichl
Department of Mathematical Sciences
University of Bath
October 2006 – January 2007
basic concepts and results in vector integration and vector calculus, Fourier series, and
the solution of partial differential equations by separation of variables.
It is
fundamental to most areas of applied Maths and many areas of pure Maths,
a prerequisite for a large number (11!!) of courses in Semester 2 and in Year 3/4.
Please take this course seriously and do not fall behind with the problem
sheets. The syllabus for this course is very dense and only through constant practice
will you be able to grasp the multitude of methods and concepts.
Please revise
The first AIM quiz is designed to help you focus on the relevant material in those courses.
Contents
1 Vector Integration 1
1.1 Line Integrals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1.1 Parametric Representation – Arclength . . . . . . . . . . . . . . . . . . . 3
1.1.2 Line Integrals of Scalar Fields . . . . . . . . . . . . . . . . . . . . . . . . 5
1.1.3 Line Integrals of Vector Fields . . . . . . . . . . . . . . . . . . . . . . . . 6
1.1.4 Application to Particle Motion . . . . . . . . . . . . . . . . . . . . . . . . 8
1.2 Surface Integrals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.2.1 The Surface Elements dS and dS . . . . . . . . . . . . . . . . . . . . . . 10
1.2.2 Surface Integrals of Scalar Fields . . . . . . . . . . . . . . . . . . . . . . 12
1.2.3 Flux – Surface Integrals of Vector Fields . . . . . . . . . . . . . . . . . . 13
1.3 Volume Integrals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.3.1 Change of Variables – Reparametrisation . . . . . . . . . . . . . . . . . . 16
2 Vector Calculus 20
2.1 Directional Derivatives and Gradients . . . . . . . . . . . . . . . . . . . . . . . . 20
2.1.1 Application: Level Surfaces and Grad . . . . . . . . . . . . . . . . . . . . 22
2.2 Divergence and Curl; the ∇–Operator . . . . . . . . . . . . . . . . . . . . . . . . 22
2.2.1 Second order derivatives – the Laplace Operator . . . . . . . . . . . . . . 23
2.2.2 Application: Potential Theory [Bourne, pp. 225–243] . . . . . . . . . . . 24
2.3 Differentiation in Curvilinear Coordinates . . . . . . . . . . . . . . . . . . . . . 25
2.3.1 Orthogonal Curvilinear Coordinates . . . . . . . . . . . . . . . . . . . . . 25
2.3.2 Grad, div, curl, ∇2 in Orthogonal Curvilinear Coordinates . . . . . . . . 27
3 Integral Theorems 30
3.1 The Divergence Theorem of Gauss . . . . . . . . . . . . . . . . . . . . . . . . . 31
3.2 Green’s Theorem in the Plane . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
3.3 Stokes’ Theorem . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
4 Fourier Series 36
4.1 Periodicity of Fourier Series . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
4.2 Fourier Convergence . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
4.3 Gibbs’ Phenomenon . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
4.4 Half-range Series . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
4.5 Application: Eigenproblems for 2nd-order ODEs . . . . . . . . . . . . . . . . . . 43
2
5.2.2 Circular Geometry . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
5.2.3 Hyperbolic PDEs – Wave Equation . . . . . . . . . . . . . . . . . . . . . 53
5.2.4 Parabolic PDEs – Diffusion Equation . . . . . . . . . . . . . . . . . . . . 56
5.3 Other Solution Methods for PDEs (not examinable!) . . . . . . . . . . . . . 56
5.3.1 The Laplace Transform Method . . . . . . . . . . . . . . . . . . . . . . . 56
5.3.2 D’Alembert’s Solution of the Wave Equation . . . . . . . . . . . . . . . . 58
Chapter 1
Vector Integration
In this course we will be dealing with functions in more than one unknown whose function
value might be vector–valued as well.
Definition 1.1. A function F : Rn → Rn is called a vector field on Rn . A function
f : Rn → R is called a scalar field on Rn . (Usually n will be 2 or 3.)
Example 1.2. Typical examples of vector fields are
• gravitational field G(x)
• electrical field E(x)
• magnetic field B(x)
• velocity field V (x)
Example 1.3. Typical examples of scalar fields are
• speed |V (x)|
1
• kinetic energy 2
m |V (x)|2
• electrical potential ϕE (x)
We will now learn how to integrate such fields in Rn .
Example 1.4. (Motivating Problem). What is the perimeter of the unit circle?
y
B
replacements 1
x
Perimeter = 4 ∗ length s of arc AB
A
Let us divide the interval [−1, 0] into N elements of length ∆xi = 1/N , i = 0, . . . , N − 1:
1
PSfrag replacements y
B
∆si ∆si
∆yi ∆yi
∆xi
A x
−1 ∆xi
meter = 4 ∗ length s of arc AB
Moreover
2
0
1
Z
s = √ dx (1.1)
−1 1 − x2
C = {r(t) : t0 ≤ t ≤ te } (1.2)
replacements
C
As t increases r(t) traces out the curve C.
r(t)
0
Note. In Ex. 1.4 t was chosen to be x. In general, t can be any parameter (e.g. angle θ).
3
x2 y 2
Example 1.5. Give a parametric representation of the ellipse 2 + 2 = 1.
a b
As in Example 1.4 we can now calculate the arclength s(t) of C from r(t0 ) to r(t) for every
t − t0
t ∈ [t0 , te ]. Let ∆t := for some N ∈ N and let ti := t0 + i∆t, for i = 0, . . . , N − 1.
N
PSfrag replacements
C
∆si
r(ti + ∆t)
r(ti )
Therefore
(1.3)
and
(1.4)
4
1.1.2 Line Integrals of Scalar Fields
Now assume we want to calculate the integral of a scalar field f : Rn → R along the curve C.
Geometrically this means, to calculate the following area A:
PSfrag replacements
z = f (x, y)
z
f (r(ti )) A
(in R2 )
y
r(te )
x C
∆si
r(t0 )
te − t 0
As in (1.4) with ∆t = we get
N
Hence for N → ∞
te
dr
Z
A = f (r(t)) dt
t0 dt
(an ordinary integral of a scalar valued function of t).
Definition 1.6. The line integral of a scalar field f : Rn → R along the curve C in (1.2)
is defined as Z te
dr
Z
f ds := f (r(t)) dt . (1.8)
C t0 dt
Remark 1.7.
R
(a) The value of C f ds does not depend on the choice of parametrisation, even if the orien-
tation of C is reversed [Bourne, p. 148], [Anton, pp. 1108–1109].
R
(b) The length of the curve C in (1.6) is a special line integral, i.e. LC = C 1 ds.
Example 1.8. Integrate f (x, y) = 2xy around the first quadrant of a circle with radius a as
shown:
y
PSfrag replacements
x
a
5
Step 1. Parametrise circle:
(dot product!)
Theorem 1.10. If T̂ is the unit tangent vector to C in (1.2) that points in the direction in
which t is increasing, then Z Z
F · dr = (F · T̂ ) ds , (1.10)
C C
i.e. the work integral of F = the line integral of the component of F parallel to C.
Proof.
6
Remark 1.11.
(a) It follows directly from Theorem 1.10 that a reversal of the orientation of C does change
the sign of the work integral (in contrast to the line integral of a scalar field, cf. Remark
1.7(a)) [Bourne, p. 153], [Anton, p. 1109].
R
(b) If F is a force field then C F · dr is the work done by moving a particle from r(t0 ) to
r(te ) along the curve C; hence the term work integral.
R
Example 1.12. (in tutorial). Evaluate the integral C F · dr where F (x, y) = (2y, x2 )T and
x2 y 2
C is the ellipse 2 + 2 = 1 in the upper half plane as shown:
a b
PSfrag replacements
y
C
b x
a
dr
Step 2. Calculate the vector line element dr = dt :
dt
7
Definition 1.13. (in tutorial).
R LetR F be a vector field such that F = F1 i + F2 j + F3 k. The
(vector valued) integrals C F ds and C F ∧ dr are defined as
Z Z Z Z T
F ds := F1 ds, F2 ds, F3 ds
C C C C
te
dr
Z Z
F ∧ dr := F (r(t)) ∧ dt .
C t0 dt
C1
x0
Recall:
dr
• particle velocity V (t) =
dt
1
• kinetic energy K(t) = 2
m |V (t)|2
8
Work done = Change in kinetic energy from t0 to te
1.2 Surface Integrals [Bourne, pp. 172–189] & [Anton, pp. 1130–1145]
You have already done the case n = 2 in MA10005, i.e. (flat) surfaces in R2 . Now we will
consider surfaces in R3 , i.e. n = 3.
PSfrag replacements
Recall (for curves): parametric representation
r(t)
C
t0 te
Now, let S be a (two-dimensional) surface in R3 with parametric representation r(u, v), i.e.
z
Sfrag replacements
r(u, v)
a v
b
b
D
v
y
u
u 2π
x
tric representation
can be parametrised by
9
1.2.1 The Surface Elements dS and dS
Recall: the equation of a plane in R3 is
x = x0 + λl + µm , λ, µ ∈ R
0
Now, let S = {r(u, v) : (u, v) ∈ D} as in (1.13) and let (u0 , v0 ) ∈ D. Fixing v0 we get
∂r
PSfrag replacements ∂v S
∂r
∂u
C u0
r(u0 , v0 ) C v0
0
Similarly,
Definition 1.16. (a) The tangent plane of S in (1.13) at the point r(u0 , v0 ) ∈ S is
3 ∂r ∂r
T := x ∈ R : x = r(u0 , v0 ) + λ (u0 , v0 ) + µ (u0 , v0 ) , λ, µ ∈ R . (1.14)
∂u ∂v
Note (on orientation). Reversing the rôles of u and v, changes the sign of n̂. This is
a matter of convention. For closed surfaces (e.g. sphere, cylinder, torus) it is common to
define n̂ to be the normal that points outward. (If it is not clear we will always specify the
orientation of the normal/surface.)
10
We will now define a surface element
dS for the surface S = {r(u, v) : (u, v) ∈ D} in (1.13):
dr
Recall the line element ds := dt which we defined in (1.6) for a
dt
PSfrag replacements
(one-dimensional) curve C.
ds C
P2 S
v + ∆v
P0 P3 D
v ∆S
r(u, v)
P1 v + ∆v
v
0 u u u + ∆u
u + ∆u
∂r ∂r
dS = ∧ du dv . (1.16)
∂u ∂v
Example 1.17. (Step 2). Calculate dS and find the outward unit normal of the cylindrical
shell of height b and radius a.
11
DIY (Recall Step 1 (i.e. Example 1.15): r(u, v) := (a cos u, a sin u, v)T .)
dS := dS n̂ , (1.17)
Example 1.19. (Step 3). Find the surface area of a cylindrical shell of height b, radius a.
[Recall: Step 1 (i.e. Example 1.15: r(u, v) := (a cos u, a sin u, v)T , 0 ≤ u ≤ 2π, 0 ≤ v ≤ b,
Step 2 (i.e. Example 1.17: dS = a du dv . ]
12
PSfrag replacements v ve (u)
v0 (u)
u
u0 ue
Remark 1.20.
(a) Another way to define the surface integral would be to
use a subdivision of S into (little) surface elements ∆Si S
with centre xi (as for line integrals), i.e.
ZZ N −1 ∆Si
f (xi ) ∆SiPSfrag replacements
X
f dS = lim . (1.20)
S N →∞
i=0 xi
See [Anton, p. 1130].
ZZ
(b) The area of a surface S is a special surface integral, i.e. AS = dS (cf. Remark 1.7(b)).
S
(b) (physical interpretation) If F is for example the velocity field of some fluid then the flux
across S is the net volume of fluid that passes through the surface per unit of time. For
more details:
DIY Please study [Anton, pp. 1137–1140] until the next lecture.
ZZ
Example 1.23. Find x · dS , where S is the sphere with radius a, i.e.
S
S := {x ∈ R3 : x2 + y 2 + z 2 = a2 }.
13
Step 1. Parametrise S: spherical polar coordinates
∂r ∂r
Step 2. Calculate ∧ :
∂θ ∂ϕ
Remark 1.24. (a) So far we have only considered smooth surfaces, i.e. surfaces S with
continuously differentiable parametric representation r. However, the notion of surface
integrals can be extended to surface integrals on piecewise smooth surfaces:
14
For example: Let S be the surface of the unit cube
and let S1 , . . . , S6 be the six faces (obviously each one
PSfrag replacements S6
of the Si , i = 1, . . . , 6 is smooth), and so
ZZ ZZ ZZ ZZ S4
= + +...+ . S2
S5 S3
S S1 S2 S6
S1
1.3 Volume Integrals [Bourne, pp. 189–194] & [Anton, pp. 1048–1090]
Ω := {x ∈ R3 : x2 + y 2 + z 2 ≤ a2 , x, y, z ≥ 0}.
y
lacements x
z
a
15
Things may become even nastier if we want to integrate over the entire ball.
(u, v, w)
r(u, v, w) = (x, y, z)T
Example 1.26. The parametrisation
D = {(ρ, θ, φ) ∈ R3 : 0 ≤ ρ ≤ a , 0 ≤ θ ≤ π , 0 ≤ φ < 2π }
to the ball with radius a. (ρ, θ, φ) are called spherical polar coordinates (ρ . . . radial dis-
tance, θ . . . polar angle, φ . . . azimuthal angle). See handout!
DIY Please study the handouts on spherical and cylindrical polar coordinates until the
replacements next lecture.
~ 0 = r(u, v, w), OP
Motivation. Let OP ~ 1 = r(u + ∆u, v, w), OP
~ 2 = r(u, v + ∆v, w), OP
~3 =
r(u, v, w + ∆w).
P3
r
w+∆w
P2 ∆V
v+∆v
w v P0 P1
u u+∆u
D Ω
16
For ∆u, ∆v, ∆w sufficiently small ∆V will be approximately a parallelepiped.
Recall (MA10006): PSfrag replacements
V
n̂ c
b A
a
17
Note. Definition 1.27 is no contradiction to the change of variables formula which you
have learnt in MA10005, since
∂r1 ∂r2 ∂r3
∂u ∂u ∂u
∂r ∂r ∂r ∂r1 ∂r2 ∂r3 ∂(x, y, z)
∂u ∧ ∂v · ∂w = det ∂v ∂v ∂v =: ∂(u, v, w) (1.26)
∂r ∂r ∂r
1 2 3
∂w ∂w ∂w
i.e. the Jacobian determinant.
[ Proof. (recall MA10006)
Remark 1.28.
(a) u, v, w can be regarded as curvilinear coordinates on Ω. (See the figure at the begin-
ning of Section 1.3.1.)
18
Step 3. Evaluate Formula (1.25) for the volume integral:
19