Math S21a: Multivariable calculus Oliver Knill, Summer 2011 The theorem explains what divergence means.

The theorem explains what divergence means. If we average the divergence over a small cube is
equal the flux of the field through the boundary of the cube. If this is positive, then more field
exists the cube than entering the cube. There is field ”generated” inside. The divergence measures
Lecture 24: Divergence theorem the expansion of the field.

1 Let F~ (x, y, z) = hx, y, zi and let S be sphere. The divergence of F~ is the constant function
RRR
There are three integral theorems in three dimensions. We have seen already the fundamental div(F~ ) = 3 and div(F~ ) dV = 3 · 4π/3 = 4π. The flux through the boundary is
RR RR
G R R
theorem of line integrals and Stokes theorem. Here is the divergence theorem, which completes r · ~ru × ~rv dudv = S |~r(u, v)|2 sin(v) dudv = 0π 02π sin(v) dudv = 4π also. We see that
S~
the list of integral theorems in three dimensions: the divergence theorem allows us to compute the area of the sphere from the volume of the
enclosed ball or compute the volume from the surface area.
Divergence Theorem. Let E be a solid with boundary surface S oriented so that
the normal vector points outside. Let F~ be a vector field. Then 2 What is the flux of the vector field F~ (x, y, z) = h2x, 3z 2 + y, sin(x)i through the solid
Z Z Z Z Z G = [0, 3] × [0, 3] × [0, 3] \ ([0, 3] × [1, 2] × [1, 2] ∪ [1, 2] × [0, 3] × [1, 2] ∪ [0, 3] × [0, 3] × [1, 2])
div(F~ ) dV = F~ · dS . which is a cube where three perpendicular cubic holes have been removed? Solution: Use
RRR RRR
E S
the divergence theorem: div(F~ ) = 2 and so ~
G div(F ) dV = 2 G dV = 2Vol(G) =
2(27 − 7) = 40. Note that the flux integral here would be over a complicated surface over
dozens of rectangular planar regions.
To prove this, one can look at a small box [x, x + dx] × [y, y + dy] × [z, z + dz]. The flux
of F~ = hP, Q, Ri through the faces perpendicular to the x-axes is [F~ (x + dx, y, z) · h1, 0, 0i + 3 Find the flux of curl(F ) through a torus if F~ = hyz 2 , z + sin(x) + y, cos(x)i and the torus
F~ (x, y, z) · h−1, 0, 0i]dydz = P (x + dx, y, z) − P (x, y, z) = Px dxdydz. Similarly, the flux through has the parametrization
the y-boundaries is Py dydxdz and the flux through the two z-boundaries is Pz dzdxdy. The
total flux through the faces of the cube is (Px + Py + Pz ) dxdydz = div(F~ ) dxdydz. A general ~r(θ, φ) = h(2 + cos(φ)) cos(θ), (2 + cos(φ)) sin(θ), sin(φ)i .
solid can be approximated as a union of small cubes. The sum of the fluxes through all the
cubes consists now of the flux through all faces without neighboring faces. and fluxes through Solution: The answer is 0 because the divergence of curl(F ) is zero. By the divergence
adjacent sides cancel. The sum of all the fluxes of the cubes is the flux through the boundary of theorem, the flux is zero.
the union. The sum of all the div(F~ ) dxdydz is a Riemann sum approximation for the integral
RRR
~
G div(F ) dxdydz. In the limit, where dx, dy, dz goes to zero, we obtain the divergence theorem.
4 Similarly as Green’s theorem allowed to calculate the area of a region by passing along the
boundary, the volume of a region can be computed as a flux integral: Take for example the
vector field F~ (x, y, z) = hx, 0, 0i which has divergence 1. The flux of this vector field through
RR
the boundary of a solid region is equal to the volume of the solid: ~
δG hx, 0, 0i· dS = Vol(G).

5 How heavy are we, at distance r from the center of the earth?
Solution: The law of gravity can be formulated as div(F~ ) = 4πρ, where ρ is the mass
density. We assume that the earth is a ball of radius R. By rotational symmetry, the
 gravitational force is normal to the surface: F~ (~x) = F~ (r)~x/||~x||. The flux of F~ through Remarks.
RR
a ball of radius r is ~ ~ 2~
RRR Sr F (x) · dS = 4πr F (r). By the divergence theorem, this is
4πMr = 4π 1) Greens theorem allows to switch from double integrals to one dimensional integrals.
Br ρ(x) dV , where Mr is the mass of the material inside Sr . We have
2) The curve is oriented in such a way that the region is to the left.
(4π)2 ρr 3 /3 = 4πr 2F~ (r) for r < R and (4π)2 ρR3 /3 = 4πr 2 F~ (r) for r ≥ R. Inside the earth,
3) The boundary of the curve can consist of piecewise smooth pieces.
the gravitational force F~ (r) = 4πρr/3. Outside the earth, it satisfies F~ (r) = M/r 2 with R
4) If C : t 7→ ~r(t) = hx(t), y(t)i, the line integral is ab hP (x(t), y(t)), Q(x(t), y(t))i ·
M = 4πR3 ρ/3.
hx (t), y (t)i dt.
′ ′

1.0
5) Green’s theorem was found by George Green (1793-1841) in 1827 and by Mikhail Ostro-
gradski (1801-1862).
6) If curl(F~ ) = 0 in a simply connected region, then the line integral along a closed curve is
0.8

0.6 zero. If two curves connect two points then the line integral along those curves agrees.
7) Taking F~ (x, y) = h−y, 0i or F~ (x, y) = h0, xi gives area formulas.
0.4

Example. Find the line integral of the vector field F~ (x, y) = hx4 + sin(x) + y, x + y 3 i along
0.2

the path ~r(t) = hcos(t), 5 sin(t) + log(1 + sin(t))i, where t runs from t = 0 to t = π.
0.5 1.0 1.5 2.0

Solution. curl(F~ ) = 0 implies that the line integral depends only on the end points
We are now at the end of the course. Lets have an overview over the integral theorems and give (0, 1), (0, −1) of the path. Take the simpler path ~r(t) = h−t, 0i, −1 ≤ t ≤ 1, which has velocity
R1 4
an other typical example in each case. ~r ′ (t) = h−1, 0i. The line integral is −1 ht − sin(t), −ti · h−1, 0i dt = −t5 /5|1−1 = −2/5.

The fundamental theorem for line integrals, Green’s theorem, Stokes theorem and divergence theo- Remark We could also find a potential f (x, y) = x5 /5 − cos(x) + xy + y 5/4. It has the
R R
rem are all incarnation of one single theorem A dF = δA F , where dF is a exterior derivative property that grad(f ) = F . Again, we get f (0, −1) − f (0, 1) = −1/5 − 1/5 = −2/5.
of F and where δA is the boundary of A. They all generalize the fundamental theorem of calculus.

Stokes theorem. If S is a surface with boundary C and F~ is a vector field, then

Fundamental theorem of line integrals: If C is a curve with boundary {A, B}
Z Z Z
and f is a function, then ~ .
curl(F~ ) · dS = F~ · dr
Z S C
~ = f (B) − f (A)
∇f · dr
C

Remarks.
Remarks. 1) Stokes theorem allows to derive Greens theorem: if F~ is z-independent and the surface S is
contained in the xy-plane, one obtains the result of Green.
R
1) For closed curves, the line integral C ∇f · dr ~ is zero. 2) The orientation of C is such that if you walk along C and have your head in the direction
2) Gradient fields are path independent: if F~ = ∇f , then the line integral between two of the normal vector ~ru × ~rv , then the surface to your left.
points P and Q does not depend on the path connecting the two points. 3) Stokes theorem was found by André Ampère (1775-1836) in 1825 and rediscovered by
3) The theorem holds Rin any dimension. In one dimension, it reduces to the fundamental George Stokes (1819-1903).
theorem of calculus ab f ′ (x) dx = f (b) − f (a) 4) The flux of the curl of a vector field does not depend on the surface S, only on the boundary
4) The theorem justifies the name conservative for gradient vector fields. of S.
5) The term ”potential” was coined by George Green who lived from 1783-1841. 5) The flux of the curl through a closed surface like the sphere is zero: the boundary of such a
Example. Let f (x, y, z) = x2 + y 4 + z. Find the line integral of the vector field surface is empty.
F~ (x, y, z) = ∇f (x, y, z) along the path ~r(t) = hcos(5t), sin(2t), t2 i from t = 0 to t = 2π.
Example. Compute the line integral of F~ (x, y, z) = hx3 + xy, y, zi along the polygonal path C
2 2
Solution. ~r(0) = h1, 0, 0i and ~r(2π) = h1, 0, 4π i and f (~r(0)) = 1 and f (~r(2π)) = 1 + 4π . The connecting the points (0, 0, 0), (2, 0, 0), (2, 1, 0), (0, 1, 0).
R Solution. The path C bounds a surface S : ~r(u, v) = hu, v, 0i parameterized by R = [0, 2] × [0, 1].
fundamental theorem of line integral gives C ∇f dr ~ = f (r(2π)) − f (r(0)) = 4π 2 .
By Stokes theorem, the line integral is equal to the flux of curl(F~ )(x, y, z) = h0, 0, −xi through
RR
S. The normal vector of S is ~rRu × ~ =
~rv = h1, 0, 0i × h0, 1, 0i = h0, 0, 1i so that S curl(F~ ) dS
Green’s theorem. If R is a region with boundary C and F~ is a vector field, then R2 R1 2 R1
Z Z Z 0 0 h0, 0, −ui · h0, 0, 1i dudv = 0 0 −u dudv = −2.

curl(F~ ) dxdy = ~ .
F~ · dr
R C
 Z Z
Divergence theorem: If S is the boundary of a region E in space and F~ is a
vector field, then Z Z Z Z Z δG
F = G
dF
div(F~ ) dV = ~ .
F~ · dS
B S

which becomes a single theorem called fundamental theorem of multivariable calculus. The
Remarks. theorem becomes much simpler in quantum calculus, where geometric objects and differential
1) The divergence theorem is also called Gauss theorem. forms are on the same footing. It tuns out that the theorem becomes then
2) It can be helpful to determine the flux of vector fields through surfaces.
3) It was discovered in 1764 by Joseph Louis Lagrange (1736-1813), later it was rediscovered by
Carl Friedrich Gauss (1777-1855) and by George Green. < δG, F >=< G, dF >
4) For divergence free vector fields F~ , the flux through a closed surface is zero. Such fields F~ are
also called incompressible or source free.
which you might see in linear algebra in the form < AT v, w >=< v, Aw >, where A is a matrix
Example. Compute the flux of the vector field F~ (x, y, z) = h−x, y, z 2 i through the boundary S and < v, w > is the dot product. If we deal with ”smooth” functions and fields that we have
of the rectangular box [0, 3] × [−1, 2] × [1, 2]. to pay a prize and consider in turn ”singular” objects like points or curves and surfaces. These
Solution. By Gauss theorem, the flux is equal to the triple integral of div(F ) = 2z over the box: are idealized objects which have zero diameter, radius or thickness. Nature likes simplicity and
R3 R2 R2
0 −1 1 2z dxdydz = (3 − 0)(2 − (−1))(4 − 1) = 27. elegance 1 : and has chosen quantum mathematics to be more fundamental but it manifests only in
the very small. While it is well understood mathematically, it will take a while until this formalism
How do these theorems fit together? In n-dimensions, there are n theorems. We have here seen will enter calculus courses.
the situation in dimension n=2 and n=3, but one could continue. The fundamental theorem of
line integrals generalizes directly to higher dimensions. Also the divergence theorem generalizes
directly since an n-dimensional integral in n dimensions. The generalization of curl and flux is
more subtle, since in 4 dimensions already, the curl of a vector field is a 6 dimensional object. It
Homework
is a n(n − 1)/2 dimensional object in general.
1 Compute using the divergence theorem the flux of the vector field F~ (x, y, z) = h3y, xy, 2yzi
1 through the unit cube [0, 1] × [0, 1] × [0, 1].

FTC 2 Find the flux of the vector field F~ (x, y, z) = hxy, yz, zxi through the solid cylinder x2 +y 2 ≤ 1,
1 −→ 1
0 ≤ z ≤ 1.
FTL Green
1 −→ 2 −→ 1 3 Use the divergence theorem to calculate the flux of F~ (x, y, z) = hx3 , y 3, z 3 i through the
FTL Stokes Gauss sphere S : x2 + y 2 + z 2 = 1 where the sphere is oriented so that the normal vector points
1 −→ 3 −→ 3 −→ 1 outwards.
In one dimensions, there is one derivative f (x) → f (x) from scalar to scalar functions. It cor-
′ Assume the vector field
responds to the entry 1 − 1 in the Pascal triangle. The next entry 1 − 2 − 1 corresponds to
differentiation in two dimensions, where we have the gradient f → ∇f mapping a scalar function F~ (x, y, z) = h5x3 + 12xy 2 , y 3 + ey sin(z), 5z 3 + ey cos(z)i
to a vector field with 2 components as well as the curl, F → curl(F ) which corresponds to the 4
is the magnetic field of the sun whose surface is a sphere
transition 2 − 1. The situation in three dimensions is captured by the entry 1 − 3 − 3 − 1 in the
Pascal triangle. The first derivative 1 − 3 is the gradient. The second derivative 3 − 3 is the curl of radius 3 oriented with the outward orientation. Com-
RR
~
pute the magnetic flux S F~ · dS.
and the third derivative 3 − 1 is the divergence. In n = 4 dimensions, we would have to look
at 1 − 4 − 6 − 4 − 1. The first derivative 1 − 4 is still the gradient. Then we have a first curl,
which maps a vector field with 4 components into an object with 6 components. Then there is a
second curl, which maps an object with 6 components back to a vector field, we would have to RR
~ where F~ (x, y, z) = hx, y, zi and S is the
Find S F~ · dS,
look at 1 − 4 − 6 − 4 − 1. When setting up calculus in dimension n, one talks about differential
forms instead of scalar fields or vector fields. Functions are 0 forms or n-forms. Vector fields can 5 boundary of the solid built with 9 unit cubes shown in
be described by 1 or n − 1 forms. The general formalism defines a derivative d called exterior the picture.
derivative on differential forms as well as integration of such k forms on k dimensional objects.
There is a boundary operation δ which maps a k-dimensional object into a k − 1 dimensional
object. This boundary operation is dual to differentiation. They both satisfy the same relation
dd(F ) = 0 and δδG = 0. Differentiation and integration are linked by the general Stokes theorem: 1
Leibniz: 1646-1716