9.13 Surface Integrals
INTRODUCTION
In the xy-plane, the length of an arc of the graph of y = f(x) from x = a to x = b is given by the definite integral
(1)
The problem in three dimensions, which is the counterpart of the arc length problem, is to find the area A(s) of that portion of the surface S given by a function z = f(x, y) having continuous first partial derivatives on a closed region R in the xy-plane. Such a surface is said to be smooth.
Surface Area
Suppose, as shown in FIGURE 9.13.1(a), that an inner partition P of R is formed using lines parallel to the x- and y-axes. P then consists of n rectangular elements Rk of area ΔAk = Δxk Δyk that lie entirely within R. Let (xk, yk, 0) denote any point in an Rk. As we see in Figure 9.13.1(a), by projecting the sides of Rk upward, we determine two quantities: a portion or patch Sk of the surface and a portion Tk of a tangent plane at (xk, yk, f(xk, yk)). It seems reasonable to assume that when Rk is small, the area ΔTk of Tk is approximately the same as the area ΔSk of Sk.
FIGURE 9.13.1 What is the area of the surface above R?
To find the area of Tk let us choose (xk, yk, 0) at a corner of Rk as shown in Figure 9.13.1(b). The indicated vectors u and v, which form two sides of Tk, are given by
where fx(xk, yk) and fy(xk, yk) are the slopes of the lines containing u and v, respectively. Now from (16) of Section 7.4 we know that ΔTk = u ×v
where
In other words,
Consequently, the area A is approximately
Taking the limit of the foregoing sum as → 0 leads us to the next definition.
DEFINITION 9.13.1 Surface Area
Let f be a function for which the first partial derivatives fx and fy are continuous on a closed region R. Then the area of the surface over R is given by
.(2)
One could have almost guessed the form of (2) by naturally extending the one-variable structure of (1) to two variables.
EXAMPLE 1 Surface Area
Find the surface area of that portion of the sphere x2 + y2 + z2 = a2 that is above the xy-plane and within the cylinder x2 + y2 = b2, 0 < b < a.
SOLUTION
If we define z = f(x, y) by f(x, y) = then
and so 1 + [fx(x, y)]2 + [fy(x, y)]2 =
Hence, (2) is A(S) =
where R is indicated in FIGURE 9.13.2. To evaluate this double integral, we change to polar coordinates:
≡
FIGURE 9.13.2 Portion of a sphere in Example 1
Differential of Surface Area
The function
(3)
is called the differential of the surface area. We will use this function in the discussion that follows.
Surface Integral
As we have seen, a double integral
is a generalization of the definite integral f(x) dx. The surface area integral (2) is a generalization of the arc length integral (1). We are now going to consider a generalization of the line integral
G(x, y) ds. This generalization is called a surface integral.
- Let w = G(x, y, z) be defined in a region of 3-space that contains a surface S, which is the graph of a function z = f(x, y). Let the projection R of the surface onto the xy-plane be either a Type I or a Type II region.
- Divide the surface S into n patches Sk with areas ΔSk that correspond to a partition P of R into n rectangles Rk with areas ΔAk.
- Let
be the norm of the partition or the length of the longest diagonal of the Rk.
- Choose a sample point (
,
,
) on each patch Sk of surface area. See FIGURE 9.13.3.
- Form the sum
G(
,
,
) ΔSk.
FIGURE 9.13.3 Sample point on kth patch
DEFINITION 9.13.2 Surface Integral
Let G be a function of three variables defined over a region of 3-space containing the surface S. Then the surface integral of G over S is given by
.(4)
Method of Evaluation
If G, f, fx, and fy are continuous throughout a region containing S, we can evaluate (4) by means of a double integral. From (3) the left side of (4) becomes
(5)
Note that when G = 1, (5) reduces to formula (2) for surface area.
Projection of S into Other Planes
If y = g(x, z) is the equation of a surface S that projects onto a region R of the xz-plane, then
(6)
Similarly, if x = h(y, z) is the equation of a surface that projects onto the yz-plane, then the analogue of (5) is
(7)
Mass of a Surface
Suppose ρ(x, y, z) represents the density of a surface at any point, or mass per unit surface area; then the mass m of the surface is
(8)
EXAMPLE 2 Mass of a Surface
Find the mass of the surface of the paraboloid z = 1 + x2 + y2 in the first octant for 1 ≤ z ≤ 5 if the density at a point P on the surface is directly proportional to its distance from the xy-plane.
SOLUTION
The surface in question and its projection onto the xy-plane are shown in FIGURE 9.13.4. Now, since ρ(x, y, z) = kz and z = 1 + x2 + y2, (8) and (5) give
FIGURE 9.13.4 Surface in Example 2
By changing to polar coordinates, we obtain
≡
EXAMPLE 3 Evaluating a Surface Integral
Evaluate xz2 dS, where S is that portion of the cylinder y = 2x2 + 1 in the first octant bounded by x = 0, x = 2, z = 4, and z = 8.
SOLUTION
We shall use (6) with g(x, z) = 2x2 + 1 and R the rectangular region in the xz-plane shown in FIGURE 9.13.5. Since gx(x, z) = 4x and gz(x, z) = 0, it follows that
≡
FIGURE 9.13.5 Surface in Example 3
Orientable Surfaces
In Example 5, we are going to evaluate a surface integral of a vector field. In order to do this we need to examine the concept of an orientable surface. Roughly, an orientable surface S, such as that given in FIGURE 9.13.6(a), has two sides that could be painted different colors. The Möbius strip* shown in Figure 9.13.6(b) is not an orientable surface and is one-sided. A person who starts to paint the surface of a Möbius strip at a point will paint the entire surface and return to the starting point.
FIGURE 9.13.6 Oriented surface in (a); non-oriented surface in (b)
Specifically, we say a smooth surface S is orientable or is an oriented surface if there exists a continuous unit normal vector function n defined at each point (x, y, z) on the surface. The vector field n(x, y, z) is called the orientation of S. But since a unit normal to the surface S at (x, y, z) can be either n(x, y, z) or –n(x, y, z), an orientable surface has two orientations. See FIGURE 9.13.7(a)–(c). The Möbius strip shown again in Figure 9.13.7(d) is not an oriented surface, since if a unit normal n starts at P on the surface and moves once around the strip on the curve C, it ends up on the “opposite side” of the strip at P and so points in the opposite direction. A surface S defined by z = f(x, y) has an upward orientation (Figure 9.13.7(b)) when the unit normals are directed upward—that is, have positive k components, and it has a downward orientation (Figure 9.13.7(c)) when the unit normals are directed downward—that is, have negative k components.
FIGURE 9.13.7 Upward orientation in (b); downward orientation in (c)
If a smooth surface S is defined by g(x, y, z) = 0, then recall that a unit normal is
(9)
where ∇g = is the gradient of g. If S is defined by z = f(x, y), then we can use g(x, y, z) = z – f(x, y) = 0 or g(x, y, z) = f(x, y) – z = 0 depending on the orientation of S.
As we shall see in the next example, the two orientations of an orientable closed surface are outward and inward.
EXAMPLE 4 Orientations of a Surface
Consider the sphere of radius a > 0: x2 + y2 + z2 = a2. If we define g(x, y, z) = x2 + y2 + z2 – a2, then
Then the two orientations of the surface are
The vector field n defines an outward orientation, whereas n1 = –n defines an inward orientation. See FIGURE 9.13.8. ≡
FIGURE 9.13.8 Sphere in Example 4
Integrals of Vector Fields
If F(x, y, z) = P(x, y, z)i + Q(x, y, z)j + R(x, y, z)k is the velocity field of a fluid, then, as we saw in Figure 9.7.3, the volume of the fluid flowing through an element of surface area ΔS per unit time is approximated by
(height)(area of base) = (compn F) ΔS = (F · n) ΔS,
where n is a unit normal to the surface. See FIGURE 9.13.9. The total volume of a fluid passing through S per unit time is called the flux of F through S and is given by
(10)
In the case of a closed surface S, if n is the outer (inner) normal, then (10) gives the volume of fluid flowing out (in) through S per unit time.
FIGURE 9.13.9 Surface S in (10)
EXAMPLE 5 Flux Through a Surface
Let F(x, y, z) = zj + zk represent the flow of a liquid. Find the flux of F through the surface S given by that portion of the plane z = 6 – 3x – 2y in the first octant oriented upward.
SOLUTION
The vector field and the surface are illustrated in FIGURE 9.13.10. By defining the plane by g(x, y, z) = 3x + 2y + z – 6 = 0, we see that a unit normal with a positive k component is
Hence,
With R the projection of the surface onto the xy-plane, we find from (10) that
≡
FIGURE 9.13.10 Surface in Example 5
Depending on the nature of the vector field, the integral in (10) can represent other kinds of flux. For example, (10) could also give electric flux, magnetic flux, flux of heat, and so on.
REMARKS
If the surface S is piecewise defined, we express a surface integral over S as the sum of the surface integrals over the various pieces of the surface. For example, suppose S is the orientable piecewise-smooth closed surface bounded by the paraboloid z = x2 + y2 (S1) and the plane z = 1 (S2). Then the flux of a vector field F out of the surface S is
where we take S1 oriented downward and S2 oriented upward. See FIGURE 9.13.11 and Problem 35 in Exercises 9.13.
FIGURE 9.13.11 Piecewise-defined surface
9.13 Exercises Answers to selected odd-numbered problems begin on page ANS-26.
- Find the surface area of that portion of the plane 2x + 3y + 4z = 12 that is bounded by the coordinate planes in the first octant.
- Find the surface area of that portion of the plane 2x + 3y + 4z = 12 that is above the region in the first quadrant bounded by the graph r = sin 2θ.
- Find the surface area of that portion of the cylinder x2 + z2 = 16 that is above the region in the first quadrant bounded on the graphs of x = 0, x = 2, y = 0, y = 5.
- Find the surface area of that portion of the paraboloid z = x2 + y2 that is below the plane z = 2.
- Find the surface area of that portion of the paraboloid z = 4 – x2 – y2 that is above the xy-plane.
- Find the surface area of those portions of the sphere x2 + y2 + z2 = 2 that are within the cone z2 = x2 + y2.
- Find the surface area of the portion of the sphere x2 + y2 + z2 = 25 that is above the region in the first quadrant bounded by the graphs of x = 0, y = 0, 4x2 + y2 = 25. [Hint: Integrate first with respect to x.]
- Find the surface area of that portion of the graph of z = x2 – y2 that is in the first octant within the cylinder x2 + y2 = 4.
- Find the surface area of the portions of the sphere x2 + y2 + z2 = a2 that are within the cylinder x2 + y2 = ay.
- Find the surface area of the portions of the cone z2 =
(x2 + y2) that are within the cylinder (x – 1)2 + y2 = 1.
- Find the surface area of the portions of the cylinder y2 + z2 = a2 that are within the cylinder x2 + y2 = a2. [Hint: See Figure 9.10.12.]
- Use the result given in Example 1 to prove that the surface area of a sphere of radius a is 4πa2. [Hint: Consider a limit as b → a.]
- Find the surface area of that portion of the sphere x2 + y2 + z2 = a2 that is bounded between y = c1 and y = c2, 0 < c1 < c2 < a. [Hint: Use polar coordinates in the xz-plane.]
- Show that the area found in Problem 13 is the same as the surface area of the cylinder x2 + z2 = a2 between y = c1 and y = c2.
In Problems 15–24, evaluate the surface integral G(x, y, z) dS.
- G(x, y, z) = x; S the portion of the cylinder z = 2 – x2 in the first octant bounded by x = 0, y = 0, y = 4, z = 0
- G(x, y, z) = xy(9 – 4z); same surface as in Problem 15
- G(x, y, z) = xz3; S the cone z =
inside the cylinder x2 + y2 = 1
- G(x, y, z) = x + y + z; S the cone z =
between z = 1 and z = 4
- G(x, y, z) = (x2 + y2)z; S that portion of the sphere x2 + y2 + z2 = 36 in the first octant
- G(x, y, z) = z2; S that portion of the plane z = x + 1 within the cylinder y = 1 – x2, 0 ≤ y ≤ 1
- G(x, y, z) = xy; S that portion of the paraboloid 2z = 4 – x2 – y2 within 0 ≤ x ≤ 1, 0 ≤ y ≤ 1
- G(x, y, z) = 2z; S that portion of the paraboloid 2z = 1 + x2 + y2 in the first octant bounded by x = 0, y =
, z = 1
- G(x, y, z) = 24
z; S that portion of the cylinder y = x2 in the first octant bounded by y = 0, y = 4, z = 0, z = 3
- G(x, y, z) = (1 + 4y2 + 4z2)1/2; S that portion of the paraboloid x = 4 – y2 – z2 in the first octant outside the cylinder y2 + z2 = 1
In Problems 25 and 26, evaluate (3z2 + 4yz) dS, where S is that portion of the plane x + 2y + 3z = 6 in the first octant. Use the projection of S onto the coordinate plane indicated in the given figure.
FIGURE 9.13.12 Region R for Problem 25
FIGURE 9.13.13 Region R for Problem 26
In Problems 27 and 28, find the mass of the given surface with the indicated density function.
- S that portion of the plane x + y + z = 1 in the first octant; density at a point P directly proportional to the square of the distance from the yz-plane
- S the hemisphere z =
In Problems 29–34, let F be a vector field. Find the flux of F through the given surface. Assume the surface S is oriented upward.
- F = xi + 2zj + yk; S that portion of the cylinder y2 + z2 = 4 in the first octant bounded by x = 0, x = 3, y = 0, z = 0
- F = zk; S that part of the paraboloid z = 5 – x2 – y2 inside the cylinder x2 + y2 = 4
- F = xi + yj + zk; same surface S as in Problem 30
- F = –x3yi + yz3j + xy3k; S that portion of the plane z = x + 3 in the first octant within the cylinder x2 + y2 = 2x
- F =
x2i +
y2j + zk; S that portion of the paraboloid z = 4 – x2 – y2 for 0 ≤ z ≤ 4
- F = eyi + exj + 18yk; S that portion of the plane x + y + z = 6 in the first octant
- Find the flux of F = y2i + x2j + 5zk out of the closed surface S given in Figure 9.13.11.
- Find the flux of F = –yi + xj + 6z2k out of the closed surface S bounded by the paraboloids z = 4 – x2 – y2 and z = x2 + y2.
- Let T(x, y, z) = x2 + y2 + z2 represent temperature and let the “flow” of heat be given by the vector field F = −∇T. Find the flux of heat out of the sphere x2 + y2 + z2 = a2. [Hint: The surface area of a sphere of radius a is 4πa2.]
- Find the flux of F = xi + yj + zk out of the unit cube 0 ≤ x ≤ 1, 0 ≤ y ≤ 1, 0 ≤ z ≤ 1. See FIGURE 9.13.14. Use the fact that the flux out of the cube is the sum of the fluxes out of the sides.
FIGURE 9.13.14 Cube in Problem 38
- Coulomb’s law states that the electric field E due to a point charge q at the origin is given by E = kqr/
r
3, where k is a constant and r = xi + yj + zk. Determine the flux out of a sphere x2 + y2 + z2 = a2.
- If σ(x, y, z) is charge density in an electrostatic field, then the total charge on a surface S is
Find the total charge on that part of the hemisphere z =
that is inside the cylinder x2 + y2 = 9 if the charge density at a point P on the surface is directly proportional to distance from the xy-plane.
- The coordinates of the centroid of a surface are given by
where A(S) is the area of the surface. Find the centroid of that portion of the plane 2x + 3y + z = 6 in the first octant. - Use the information in Problem 41 to find the centroid of the hemisphere z =
.
- Let z = f(x, y) be the equation of a surface S and F be the vector field F(x, y, z) = P(x, y, z)i + Q(x, y, z)j + R(x, y, z)k. Show that
(F · n) dS equals
* To construct a Möbius strip cut out a long strip of paper, give one end a half-twist, and then attach the ends with tape.