V9. Surface Integrals: 1. The Surface Integral For Flux
V9. Surface Integrals: 1. The Surface Integral For Flux
V9. Surface Integrals: 1. The Surface Integral For Flux
Surface Integrals
Surface integrals are a natural generalization of line integrals: instead of integrating over
a curve, we integrate over a surface in 3-space. Such integrals are important in any of the
subjects that deal with continuous media (solids, fluids, gases), as well as subjects that deal
with force fields, like electromagnetic or gravitational fields.
Though most of our work will be spent seeing how surface integrals can be calculated and
what they are used for, we first want to indicate briefly how they are defined. The surface
integral of the (continuous) function f (x, y, z) over the surface S is denoted by
You can think of d S as the area of an infinitesimal piece of the surface S . To define the
integral ( I ) , we subdivide the surface S into small pieces having area ASi, pick a point
(xi, yi, zi) in the i-th piece, and form the Riemann sum
As the subdivision of S gets finer and finer, the corresponding sums (2) approach a limit
which does not depend on the choice of the points or how the surface was subdivided. The
surface integral (1) is defined to be this limit. (The surface has to be smooth and not infinite
in extent, and the subdivisions have to be made reasonably, otherwise the limit may not
exist, or it may not be unique.)
the two integrals are the same, but the second is written using the common
and suggestive abbreviation dS = n d S .
If F represents the velocity field for the flow of an incompressible fluid of density 1, then
F . n represents the component of the velocity in the positive perpendicular direction to the
surface, and F . n d S represents the flow rate across the little infinitesimal piece of surface
2 V. VECTOR INTEGRAL CALCLUS
having area dS. The integral in (3) adds up these flows across the pieces of surface, so that
we may interpret (3) as saying
(4) flux of F through S = net flow rate across S,
where we count flow in the direction of n as positive, flow in the opposite direction as
negative. More generally, if the fluid has varying density, then the right side of (4) is the
net mass transport rate of fluid across S (per unit area, per time unit).
If F is a force field, then nothing is physically flowing, and one just uses the term "flux"
to denote the surface integral, as in (3).
We now show how to calculate the flux integral. It takes a few steps, and is best done by
examples. The basic things are finding n and dS. We will begin with two surfaces where
these are easy to calculate - the cylinder and the sphere. Then we will consider a general
surface.
As the parameters 8 and z vary, the whole cylinder is traced out ; the piece we want satisfies
0 5 8 5 ~ / 2 0, 5 z 5 h . The natural way to subdivide the cylinder is to use little pieces
of curved rectangle like the one shown, bounded by two horizontal circles and two vertical
lines on the surface. Its area dS is the product of its height and width:
Having obtained n and dS, the rest of the work is routine. We express the integrand of
our surface integral (3) in terms of z and 8:
ah2
inner integral = -cos 8
2
+ a2h sin 8 cos 8
outer integral = sin 8 + a2h-s i8 1 ~
= -(a + h) .
V9. SURFACE INTEGRALS 3
+
Example 2. Find the flux of F = xz i yz j + z2 k outward through that part of the
+ +
sphere x2 y2 z2 = a2 lying in the first octant (x, y, z, 2 0).
Solution. Once again, we begin by finding n and d S for the sphere. We take the
outside of the sphere as the positive side, so n points radially outward from the origin; we
see by inspection therefore that
(9) d S = a2 sinq5dq5d8.
Finally since the area element dS is expressed in terms of q5 and 8, the integration will
be done using these variables, which means we need to express x, y, z in terms of q5 and 8.
We use the formulas expressing Cartesian in terms of spherical coordinates (setting p = a
since (x, y, z) is on the sphere):
We can now calculate the flux integral (3). By (8) and (9), the integrand is
VF
n=f- (choose the right sign);
l VFl
n=
+
-fzi - f g j k
(if n points "up")
(f:+ +
fy2 1)lI2
The expression (11) for n uses the fact that the gradient V F is normal to the surface F = c;
we divided by lVFl to make n a unit vector. To get the second form (ll'), write the surface
in the form z - f (x, y) = 0, and then use ( l l ) , taking F ( x , y, z) to be z - f (x, y).
As an example of the use of ( l l ) , for the sphere x2 + y 2 + z2 = a2, it gives immediately
the unit normal vector
1
n = -(xi
a
+ yj +zk) .
Most often however it is (11') that is used. We will need therefore the
expression for d S in xy-coordinates. A natural choice for the infinitesimal
element of surface is the infinitesimal piece of S lying over the area element
dx dy in the xy-plane. The area of these two pieces are related by the basic
formula
dx dy
(12) d S = -.
In. k l
rJk.
The argument for (12) is as follows. Since d S is small, we can suppose to a
first-order approximation that it lies in a plane P . Then the formula (12) is just
a special case of a useful result relating the area of region R in the xy-plane with
the area of the region R' lying above it in the plane P, namely:
area R
(13) area R' = - ( n = unit normal to P ) .
In. k l
We can see that (13) is true in two steps:
Since (13) is true for each of the little pieces above and below, it is true for
their corresponding sums as well, and these two sums approximate the area of R
and R' arbitrarily closely. This completes the argument for (13), and therefore
also for (12).
If we apply (11') to (12), we get the explicit expression
One can think of this as the two-dimensional analog of the one-dimensional formula for the
infinitesimal arclength along a curve y = f (x):
Actually, for calculating the surface integral representing flux, what is needed is not d S
alone, but rather the combination n d S . If we combine (11') and (129, the square roots
miraculously cancel each other out, and we end up with the much simpler-looking vector
formula
and this is the formula you will use most frequently in calculating surface integrals for the
flux through a general surface given by z = f (x, y).
+
E x a m p l e 3. The portion of the plane 22 - 2y z = 1 lying in the first octant forms
+ +
a triangle S . Find the flux of F = x i y j z k through S; take the positive side of S as
the one where the normal points "up".
Solution. Writing the plane in the form z = 1- 22 + 2y, we get by (13),
where R is the region in the xy-plane over which S lies. (Note that since the integration
is to be in terms of x and y, we had to express z in terms of x and y for this last step.)
To see what R is explicitly, the plane intersects the three coordinate axes respectively at
x = 112, y = -112, z = 1 . So R is the region pictured; our integral has integrand 1, so its
value is the area of R, which is 118.
R e m a r k . When we write z = f (x, y), we are agreeing to parametrize our surface using
x and y as parameters. Thus the flux integral will be reduced to a double integral over
a region R in the xy-plane, involving only x and y. Therefore you must get rid of z by
using the relation z = f (x, y) after you have calculated the flux integral using (13). Then
determine R (the projection of S onto the xy-plane), and supply the limits for the iterated
integral over R.
6 V. VECTOR INTEGRAL CALCLUS
that we introduced at the beginning of this section can be used for many things other than
calculating flux. Here are some examples.
average value of f on s = -
area S
JJs f ( x , Y,z) d s .
In evaluating such integrals as (14), if the surface is a sphere or cylinder, then d S should
be written down using the formulas in section 2 of V9, while if the surface is a more general
one, given by z = f (x, y), then one uses the results of section 3 instead - unfortunately,
the radical which appears in d S (see (12') will not disappear, since there is no n to cancel
it.
Example 4. Find the average distance along the earth of the points in the northern
hemisphere from the North Pole. (Assume the earth is a sphere of radius a.)
Solution. - We use (15) and spherical coordinates, choosing the co-
ordinates so the North Pole is at z = a on the z-axis. The distance of
the point (a, ), 19) from (a, 0,O) is a), measured along the great circle, i.e.,
the longitude line - see the picture). We want to find the average of this
function over the upper hemisphere S. Integrating, and using (9), we get
JL a) d S = 12=
lTi2 a)a2 sin )d) dI9 = 2ra3
lTi2 ) sin ) d) = 2 r a 3 .
(The last integral used integration by parts.) Since the area of S = 27ra2, we get using (15)
the striking answer: average distance = a .
Exercises: Section 6B