Calculus Iii: Surface Integrals
Calculus Iii: Surface Integrals
Calculus Iii: Surface Integrals
Surface Integrals
Paul Dawkins
Calculus III
Table of Contents
Preface............................................................................................................................................. ii Surface Integrals ............................................................................................................................ 3
Introduction ................................................................................................................................................ 3 Parametric Surfaces .................................................................................................................................... 4 Surface Integrals....................................................................................................................................... 10 Surface Integrals of Vector Fields ............................................................................................................ 19 Stokes Theorem ...................................................................................................................................... 29 Divergence Theorem ................................................................................................................................ 34
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
Preface
Here are my online notes for my Calculus III course that I teach here at Lamar University. Despite the fact that these are my class notes, they should be accessible to anyone wanting to learn Calculus III or needing a refresher in some of the topics from the class. These notes do assume that the reader has a good working knowledge of Calculus I topics including limits, derivatives and integration. It also assumes that the reader has a good knowledge of several Calculus II topics including some integration techniques, parametric equations, vectors, and knowledge of three dimensional space. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed. 1. Because I wanted to make this a fairly complete set of notes for anyone wanting to learn calculus I have included some material that I do not usually have time to cover in class and because this changes from semester to semester it is not noted here. You will need to find one of your fellow class mates to see if there is something in these notes that wasnt covered in class. 2. In general I try to work problems in class that are different from my notes. However, with Calculus III many of the problems are difficult to make up on the spur of the moment and so in this class my class work will follow these notes fairly close as far as worked problems go. With that being said I will, on occasion, work problems off the top of my head when I can to provide more examples than just those in my notes. Also, I often dont have time in class to work all of the problems in the notes and so you will find that some sections contain problems that werent worked in class due to time restrictions. 3. Sometimes questions in class will lead down paths that are not covered here. I try to anticipate as many of the questions as possible in writing these up, but the reality is that I cant anticipate all the questions. Sometimes a very good question gets asked in class that leads to insights that Ive not included here. You should always talk to someone who was in class on the day you missed and compare these notes to their notes and see what the differences are. 4. This is somewhat related to the previous three items, but is important enough to merit its own item. THESE NOTES ARE NOT A SUBSTITUTE FOR ATTENDING CLASS!! Using these notes as a substitute for class is liable to get you in trouble. As already noted not everything in these notes is covered in class and often material or insights not in these notes is covered in class.
ii
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
Surface Integrals
Introduction
In the previous chapter we looked at evaluating integrals of functions or vector fields where the points came from a curve in two- or three-dimensional space. We now want to extend this idea and integrate functions and vector fields where the points come from a surface in threedimensional space. These integrals are called surface integrals. Here is a list of the topics covered in this chapter. Parametric Surfaces In this section we will take a look at the basics of representing a surface with parametric equations. We will also take a look at a couple of applications. Surface Integrals Here we will introduce the topic of surface integrals. We will be working with surface integrals of functions in this section. Surface Integrals of Vector Fields We will look at surface integrals of vector fields in this section. Stokes Theorem We will look at Stokes Theorem in this section. Divergence Theorem Here we will take a look at the Divergence Theorem.
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
Parametric Surfaces
Before we get into surface integrals we first need to talk about how to parameterize a surface. When we parameterized a curve we took values of t from some interval [ a, b ] and plugged them into
r r r r r (t ) = x (t ) i + y (t ) j + z (t ) k
and the resulting set of vectors will be the position vectors for the points on the curve. With surfaces well do something similar. We will take points, ( u , v ) , out of some twodimensional space D and plug them into
r r r r r ( u, v ) = x ( u, v ) i + y ( u, v ) j + z ( u, v ) k
and the resulting set of vectors will be the position vectors for the points on the surface S that we are trying to parameterize. This is often called the parametric representation of the parametric surface S. We will sometimes need to write the parametric equations for a surface. There are really nothing more than the components of the parametric representation explicitly written down.
x = x ( u, v )
y = y ( u, v )
z = z ( u, v )
Example 1 Determine the surface given by the parametric representation r r r r r ( u, v ) = u i + u cos v j + u sin v k
Solution Lets first write down the parametric equations.
x=u
y = u cos v
z = u sin v
So, we were able to eliminate the parameters and the equation in x, y, and z is given by,
From the Quadric Surfaces section notes we can see that this is a cone that opens along the x-axis. We are much more likely to need to be able to write down the parametric equations of a surface than identify the surface from the parametric representation so lets take a look at some examples of this.
Example 2 Give parametric representations for each of the following surfaces. (a) The elliptic paraboloid x = 5 y 2 + 2 z 2 - 10 . [Solution] (b) The elliptic paraboloid x = 5 y 2 + 2 z 2 - 10 that is in front of the yz-plane.
[Solution]
Calculus III
Solution (a) The elliptic paraboloid x = 5 y 2 + 2 z 2 - 10 . This one is probably the easiest one of the four to see how to do. Since the surface is in the form x = f ( y, z ) we can quickly write down a set of parametric equations as follows,
x = 5 y 2 + 2 z 2 - 10
y=y
z=z
The last two equations are just there to acknowledge that we can choose y and z to be anything we want them to be. The parametric representation is then,
r r r r r ( y, z ) = ( 5 y 2 + 2 z 2 - 10 ) i + y j + z k
[Return to Problems]
(b) The elliptic paraboloid x = 5 y 2 + 2 z 2 - 10 that is in front of the yz-plane. This is really a restriction on the previous parametric representation. The parametric representation stays the same.
r r r r r ( y, z ) = ( 5 y 2 + 2 z 2 - 10 ) i + y j + z k
However, since we only want the surface that lies in front of the yz-plane we also need to require that x 0 . This is equivalent to requiring,
5 y 2 + 2 z 2 - 10 0
(c) The sphere x 2 + y 2 + z 2 = 30 .
or
5 y 2 + 2 z 2 10
[Return to Problems]
This one can be a little tricky until you see how to do it. In spherical coordinates we know that the equation of a sphere of radius a is given by,
r =a
and so the equation of this sphere (in spherical coordinates) is r = 30 . Now, we also have the following conversion formulas for converting Cartesian coordinates into spherical coordinates.
x = r sin j cos q
y = r sin j sin q
z = r cos j
However, we know what r is for our sphere and so if we plug this into these conversion formulas we will arrive at a parametric representation for the sphere. Therefore, the parametric representation is,
All we need to do now is come up with some restriction on the variables. First we know that we have the following restriction.
0 j p
This is enforced upon us by choosing to use spherical coordinates. Also, to make sure that we only trace out the sphere once we will also have the following restriction.
0 q 2p
5
Calculus III
(d) The cylinder y 2 + z 2 = 25 . As with the last one this can be tricky until you see how to do it. In this case it makes some sense to use cylindrical coordinates since they can be easily used to write down the equation of a cylinder. In cylindrical coordinates the equation of a cylinder of radius a is given by and so the equation of the cylinder in this problem is r = 5 . Next, we have the following conversion formulas.
r=a
x=x
y = r sin q
z = r cos q
Notice that they are slightly different from those that we are used to seeing. We needed to change them up here since the cylinder was centered upon the x-axis. Finally, we know what r is so we can easily write down a parametric representation for this cylinder.
r r r r r ( x, q ) = x i + 5sin q j + 5cos q k
We will also need the restriction 0 q 2p to make sure that we dont retrace any portion of the cylinder. Since we havent put any restrictions on the height of the cylinder there wont be any restriction on x.
[Return to Problems]
In the first part of this example we used the fact that the function was in the form x = f ( y, z ) to quickly write down a parametric representation. This can always be done for functions that are in this basic form.
z = f ( x, y ) x = f ( y, z ) y = f ( x, z )
r r r r r ( x, y ) = x i + y j + f ( x, y ) k r r r r r ( x, y ) = f ( y , z ) i + y j + z k r r r r r ( x, y ) = x i + f ( x, z ) j + z k
Okay, now that we have practice writing down some parametric representations for some surfaces lets take a quick look at a couple of applications. Lets take a look at finding the tangent plane to the parametric surface S given by,
r r r r r ( u, v ) = x ( u, v ) i + y ( u, v ) j + z ( u, v ) k
First, define
r r y r z r x ru ( u , v ) = ( u , v ) i + ( u , v ) j + ( u, v ) k u u u r r y r z x r rv ( u , v ) = ( u , v ) i + ( u , v ) j + ( u , v ) k v v v
6 http://tutorial.math.lamar.edu/terms.aspx
Calculus III
Now, provided ru rv 0 it can be shown that the vector ru rv will be orthogonal to the surface S. This means that it can be used for the normal vector that we need in order to write down the equation of a tangent plane. This is an important idea that will be used many times throughout the next couple of sections. Lets take a look at an example.
Example 3 Find the equation of the tangent plane to the surface given by r r r r r ( u, v ) = u i + 2v 2 j + ( u 2 + v ) k
at the point ( 2, 2,3) . Solution r r Lets first compute ru rv . Here are the two individual vectors.
r r r ru ( u , v ) = i + 2u k
r r r rv ( u, v ) = 4v j + k
Now the cross product (which will give us the normal vector n ) is,
r r r r i j k r r r r r r n = ru rv = 1 0 2u = -8uv i - j + 4v k 0 4v 1
Now, this is all fine, but in order to use it we will need to determine the value of u and v that will give us the point in question. We can easily do this by setting the individual components of the parametric representation equal to the coordinates of the point in question. Doing this gives,
2=u
u=2
2 = 2v
v = 1
3 = u2 + v
Now, as shown, we have the value of u, but there are two possible values of v. To determine the correct value of v lets plug u into the third equation and solve for v. This should tell us what the correct value is.
3= 4+v
v = -1
Okay so we now know that well be at the point in question when u = 2 and v = -1 . At this point the normal vector is,
r r r r n = 16 i - j - 4 k
16 ( x - 2 ) - ( y - 2 ) - 4 ( z - 3) = 0 16 x - y - 4 z = 18
You do remember how to write down the equation of a plane, right? The second application that we want to take a quick look at is the surface area of the parametric surface S given by,
2007 Paul Dawkins 7 http://tutorial.math.lamar.edu/terms.aspx
Calculus III
r r A = ru rv dA
D
Example 4 Find the surface area of the portion of the sphere of radius 4 that lies inside the cylinder x 2 + y 2 = 12 and above the xy-plane.
Solution Okay weve got a couple of things to do here. First we need the parameterization of the sphere. We parameterized a sphere earlier in this section so there isnt too much to do at this point. Here is the parameterization.
Next we need to determine D. Since we are not restricting how far around the z-axis we are rotating with the sphere we can take the following range for q .
0 q 2p
Now, we need to determine a range for j . This will take a little work, although its not too bad. First, lets start with the equation of the sphere.
x 2 + y 2 + z 2 = 16
Now, if we substitute the equation for the cylinder into this equation we can find the value of z where the sphere and the cylinder intersect.
x 2 + y 2 + z 2 = 16 12 + z 2 = 16 z2 = 4 z = 2
Now, since we also specified that we only want the portion of the sphere that lies above the xyplane we know that we need z = 2 . We also know that r = 4 . Plugging this into the following conversion formula we get,
j=
p 3
p 3
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
r i
4sin j cos q 0 4 cos j sin q -4sin j r r r r = -16sin 2 j cos q i - 16sin j cos j sin 2 q k - 16sin 2 j sin q j - 16sin j cos j cos 2 q k r r r = -16sin 2 j cos q i - 16sin 2 j sin q j - 16sin j cos j ( sin 2 q + cos 2 q ) k r r r = -16sin 2 j cos q i - 16sin 2 j sin q j - 16sin j cos j k r r rq rj = 256sin 4 j cos 2 q + 256sin 4 j sin 2 q + 256sin 2 j cos 2 j = 256sin 4 j ( cos 2 q + sin 2 q ) + 256sin 2 j cos 2 j = 256sin 2 j ( sin 2 j + cos 2 j ) = 16 sin 2 j = 16 sin j = 16sin j
We can drop the absolute value bars in the sine because sine is positive in the range of j that we are working with. We can finally get the surface area.
A = 16sin j dA
D
= =
2p 0 2p 0 2p
p 3 0
16sin j dj dq
p 3
-16 cos j 0 dq
= 8 dq
0
= 16p
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
Surface Integrals
It is now time to think about integrating functions over some surface, S, in three-dimensional space. Lets start off with a sketch of the surface S since the notation can get a little confusing once we get into it. Here is a sketch of some surface S.
The region S will lie above (in this case) some region D that lies in the xy-plane. We used a rectangle here, but it doesnt have to be of course. Also note that we could just as easily looked at a surface S that was in front of some region D in the yz-plane or the xz-plane. Do not get so locked into the xy-plane that you cant do problems that have regions in the other two planes. Now, how we evaluate the surface integral will depend upon how the surface is given to us. There are essentially two separate methods here, although as we will see they are really the same. First, lets look at the surface integral in which the surface S is given by z = g ( x, y ) . In this case the surface integral is,
2 2 g g f ( x, y, z ) dS = f ( x, y, g ( x, y ) ) + + 1 dA x y D
Now, we need to be careful here as both of these look like standard double integrals. In fact the integral on the right is a standard double integral. The integral on the left however is a surface integral. The way to tell them apart is by looking at the differentials. The surface integral will have a dS while the standard double integral will have a dA. In order to evaluate a surface integral we will substitute the equation of the surface in for z in the integrand and then add on the often messy square root. After that the integral is a standard double integral and by this point we should be able to deal with that.
10
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
Note as well that there are similar formulas for surfaces given by y = g ( x, z ) (with D in the xzplane) and x = g ( y, z ) (with D in the yz-plane). We will see one of these formulas in the examples and well leave the other to you to write down. The second method for evaluating a surface integral is for those surfaces that are given by the parameterization,
r r r r r ( u, v ) = x ( u, v ) i + y ( u, v ) j + z ( u, v ) k r r r
f ( x, y, z ) dS = f ( r ( u, v ) ) r r
u S D
dA
where D is the range of the parameters that trace out the surface S. Before we work some examples lets notice that since we can parameterize a surface given by z = g ( x, y ) as,
r r r r r ( x, y ) = xi + yj + g ( x, y ) k
2 2
we can always use this form for these kinds of surfaces as well. In fact it can be shown that,
r r g g rx ry = + + 1 x y
for these kinds of surfaces. You might want to verify this for the practice of computing these cross products. Lets work some examples.
Example 1 Evaluate
front of the yz-plane.
6 xy dS
S
Solution Okay, since we are looking for the portion of the plane that lies in front of the yz-plane we are going to need to write the equation of the surface in the form x = g ( y, z ) . This is easy enough to do.
x = 1- y - z
Next we need to determine just what D is. Here is a sketch of the surface S.
11
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
Notice that the axes are labeled differently than we are used to seeing in the sketch of D. This was to keep the sketch consistent with the sketch of the surface. We arrived at the equation of the hypotenuse by setting x equal to zero in the equation of the plane and solving for z. Here are the ranges for y and z.
0 y 1
0 z 1- y
Now, because the surface is not in the form z = g ( x, y ) we cant use the formula above. However, as noted above we can modify this formula to get one that will work for us. Here it is,
2 2 g g f ( x, y, z ) dS = f ( g ( y, z ) , y, z ) 1 + + dA y z D
The changes made to the formula should be the somewhat obvious changes. So, lets do the integral.
6 xy dS = 6 (1 - y - z ) y
S D
1 + ( -1) + ( -1) dA
2 2
Notice that we plugged in the equation of the plane for the x in the integrand. At this point weve got a fairly simple double integral to do. Here is that work.
2007 Paul Dawkins 12 http://tutorial.math.lamar.edu/terms.aspx
Calculus III
6 xy dS = 3 6 ( y - y
S D 1 1- y 1
- zy ) dA
= 6 3 y - y 2 - zy dz dy 0 0 1 = 6 3 yz - zy 2 - z 2 y dy 2 0 0
1 1 1 = 6 3 y - y 2 + y 3 dy 0 2 2 1- y
1 1 3 1 = 6 3 y 2 - y3 + y 4 = 3 8 0 4 4 Example 2 Evaluate
z dS
S
Solution We gave the parameterization of a sphere in the previous section. Here is the parameterization for this sphere.
Since we are working on the upper half of the sphere here are the limits on the parameters.
p 2
r r j k 2sin j cos q 0 2 cos j sin q -2sin j r r r r = -4sin 2 j cos q i - 4sin j cos j sin 2 q k - 4sin 2 j sin q j - 4sin j cos j cos 2 q k r r r = -4sin 2 j cos q i - 4sin 2 j sin q j - 4sin j cos j ( sin 2 q + cos 2 q ) k r r r = -4sin 2 j cos q i - 4sin 2 j sin q j - 4sin j cos j k
13
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
r r rq rj = 16sin 4 j cos 2 q + 16sin 4 j sin 2 q + 16sin 2 j cos 2 j = 16sin 4 j ( cos 2 q + sin 2 q ) + 16sin 2 j cos 2 j = 16sin 2 j ( sin 2 j + cos 2 j ) = 4 sin 2 j = 4 sin j = 4sin j
We can drop the absolute value bars in the sine because sine is positive in the range of j that we are working with. The surface integral is then,
z dS = 2 cos j ( 4sin j ) dA
S D
Dont forget that we need to plug in for x, y and/or z in these as well, although in this case we just needed to plug in z. Here is the evaluation for the double integral.
z dS =
S 0
2p
p 2 0
4sin ( 2j ) dj dq
p 2 0
= =
2p 0 2p 0
( -2 cos ( 2j ) )
4 dq
dq
= 8p Example 3 Evaluate
y dS
S
between z = 0 and z = 6 . Solution We parameterized up a cylinder in the previous section. Here is the parameterization of this cylinder.
r r r r r ( z ,q ) = 3 cos q i + 3 sin q j + z k
0 z6 0 q 2p
r r rz ( z , q ) = k r r r rq ( z ,q ) = - 3 sin q i + 3 cos q j
14
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
r r rz rq =
r i 0 - 3 sin q
r j 0
r k 1
r r rz rq = 3cos 2 q + 3sin 2 q = 3
y dS =
S D
3 sin q
2p 0 2p 0
( 3 ) dA
= 3 = 3 =0 Example 4 Evaluate
6 0
sin q dz dq
6sin q dq
2p
= ( -18cos q ) 0
y + z dS
S
whose bottom is the disk x 2 + y 2 3 in the xy-plane and whose top is the plane z = 4 - y . Solution There is a lot of information that we need to keep track of here. First, we are using pretty much the same surface (the integrand is different however) as the previous example. However, unlike the previous example we are putting a top and bottom on the surface this time. Lets first start out with a sketch of the surface.
15
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
Actually we need to be careful here. There is more to this sketch than the actual surface itself. Were going to let S1 be the portion of the cylinder that goes from the xy-plane to the plane. In other words, the top of the cylinder will be at an angle. Well call the portion of the plane that lies inside (i.e. the cap on the cylinder) S2 . Finally, the bottom of the cylinder (not shown here) is the disk of radius
In order to do this integral well need to note that just like the standard double integral, if the surface is split up into pieces we can also split up the surface integral. So, for our example we will have,
y + z dS = y + z dS + y + z dS + y + z dS
S S1 S2 S3
Were going to need to do three integrals here. However, weve done most of the work for the first one in the previous example so lets start with that.
S1 : The Cylinder
The parameterization of the cylinder and rz rq is,
r r r r r ( z ,q ) = 3 cos q i + 3 sin q j + z k
r r rz rq = 3
The difference between this problem and the previous one is the limits on the parameters. Here they are.
0 q 2p
0 z 4 - y = 4 - 3 sin q
16
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
The upper limit for the zs is the plane so we can just plug that in. However, since we are on the cylinder we know what y is from the parameterization so we will also need to plug that in. Here is the integral for the cylinder.
y + z dS = (
S1 D
3 sin q + z
2p 0 2p 0
)( 3 ) dA
3 sin q + z dz dq
= 3 = 3
4 - 3 sin q 0
3 sin q 4 - 3 sin q +
) 1 (4 2
3 sin q
dq
2p 3 = 3 8 - sin 2 q dq 0 2 2p 3 = 3 8 - (1 - cos ( 2q ) ) dq 0 4
3 29 = 3 q + sin ( 2q ) 4 8 0 = 29 3 p 2
2p
y + z dS = ( y + 4 - y ) ( 0 ) + ( -1)
2 S2 D
+ 1 dA
= 2 4 dA
D
Dont forget that we need to plug in for z! Now at this point we can proceed in one of two ways. Either we can proceed with the integral or we can recall that of D and we know that D is the disk of radius
17
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
y + z dS = 4 2 dA
S2
( 3) )
2
Again, this is set up to use the initial formula we gave in this section once we realize that the equation for the bottom is given by g ( x, y ) = 0 and D is the disk of radius origin. Also, dont forget to plug in for z. Here is the work for this integral.
3 centered at the
y + z dS = ( y + 0 ) ( 0 ) + ( 0 ) + (1)
2 2 S3 D
dA
= y dA
D
= 0
2p
r 2 sin q dr dq
3 1 3 = r sin q dq 0 0 3
2p
2p
3 sin q dq
2p 0
= - 3 cos q =0
We can now get the value of the integral that we are after.
y + z dS = y + z dS + y + z dS + y + z dS
S S1 S2 S3
29 3 p + 12 2 p + 0 2 p = 29 3 + 24 2 2 =
18
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
f ( x, y , z ) = z - g ( x, y )
In terms of our new function the surface is then given by the equation f ( x, y, z ) = 0 . Now, recall that f will be orthogonal (or normal) to the surface given by f ( x, y, z ) = 0 . This means that we have a normal vector to the surface. The only potential problem is that it might not be a unit normal vector. That isnt a problem since we also know that we can turn any vector into a unit vector by dividing the vector by its length. In or case this is,
r f n= f
In this case it will be convenient to actually compute the gradient vector and plug this into the formula for the normal vector. Doing this gives,
r r r - gx i - g y j + k r f n= = 2 2 f ( gx ) + ( g y ) +1
Now, from a notational standpoint this might not have been so convenient, but it does allow us to make a couple of additional comments.
19
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
First, notice that the component of the normal vector in the z-direction (identified by the k in the normal vector) is always positive and so this normal vector will generally point upwards. It may not point directly up, but it will have an upwards component to it. This will be important when we are working with a closed surface and we want the positive orientation. If we know that we can then look at the normal vector and determine if the positive orientation should point upwards or downwards. Remember that the positive orientation must point out of the region and this may mean downwards in places. Of course if it turns out that we need the downward orientation we can always take the negative of this unit vector and well get the one that we need. Again, remember that we always have that option when choosing the unit normal vector. Before we move onto the second method of giving the surface we should point out that we only did this for surfaces in the form z = g ( x, y ) . We could just as easily done the above work for surfaces in the form y = g ( x, z ) (so f ( x, y, z ) = y - g ( x, z ) ) or for surfaces in the form
x = g ( y, z ) (so f ( x, y, z ) = x - g ( y, z ) ).
Now, we need to discuss how to find the unit normal vector if the surface is given parametrically as,
r r r r r ( u, v ) = x ( u, v ) i + y ( u, v ) j + z ( u, v ) k r r In this case recall that the vector ru rv will be normal to the tangent plane at a particular point.
But if the vector is normal to the tangent plane at a point then it will also be normal to the surface at that point. So, this is a normal vector. In order to guarantee that it is a unit normal vector we will also need to divide it by its magnitude. So, in the case of parametric surfaces one of the unit normal vectors will be,
r r r ru r n = r rv ru rv
As with the first case we will need to look at this once its computed and determine if it points in the correct direction or not. If it doesnt then we can always take the negative of this vector and that will point in the correct direction. Finally, remember that we can always parameterize any surface given by z = g ( x, y ) (or
F gdS = F gn dS
S S
r r
20
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
where the right hand integral is a standard surface integral. This is sometimes called the flux of r F across S. Before we work any examples lets notice that we can substitute in for the unit normal vector to get a somewhat easier formula to use. We will need to be careful with each of the following formulas however as each will assume a certain orientation and we may have to change the normal vector to match the given orientation. Lets first start by assuming that the surface is given by z = g ( x, y ) . In this case lets also assume that the vector field is given by F = P i + Q j + R k and that the orientation that we are after is the upwards orientation. Under all of these assumptions the surface integral of F over S is,
r r r r F gdS = F gn dS
S
r r r r r P i + Q r + R k g - g x i - g y j + k = j 2 2 ( gx ) + ( g y ) + 1 D r r r r r r = P i + Q j + R k g - g x i - g y j + k dA
( gx ) + ( g y )
2
+ 1 dA
)(
= - Pg x - Qg y + R dA
D
Now, remember that this assumed the upward orientation. If wed needed the downward orientation then we would need to change the signs on the normal vector. This would in turn change the signs on the integrand as well. So, we really need to be careful here when using this formula. In general it is best to rederive this formula as you need it. When weve been given a surface that is not in parametric form there are in fact 6 possible integrals here. Two for each form of the surface z = g ( x, y ) , y = g ( x, z ) and x = g ( y, z ) . Given each form of the surface there will be two possible unit normal vectors and well need to choose the correct one to match the given orientation of the surface. However, the derivation of each formula is similar to that given here and so shouldnt be too bad to do as you need to. Notice as well that because we are using the unit normal vector the messy square root will always drop out. This means that when we do need to derive the formula we wont really need to put this in. All well need to work with is the numerator of the unit vector. We will see at least one more of these derived in the examples below. It should also be noted that the square root is nothing more than,
( gx ) + ( g y )
2
+ 1 = f
so in the following work we will probably just use this notation in place of the square root when we can to make things a little simpler.
2007 Paul Dawkins 21 http://tutorial.math.lamar.edu/terms.aspx
Calculus III
Lets now take a quick look at the formula for the surface integral when the surface is given r parametrically by r ( u , v ) . In this case the surface integral is,
F gdS = F gn dS
S
r r
r r r ru rv r r = F g r r ru rv dA ru rv D r r r = F g( ru rv ) dA
D
Again note that we may have to change the sign on ru rv to match the orientation of the surface and so there is once again really two formulas here. Also note that again the magnitude cancels in this case and so we wont need to worry that in these problems either. Note as well that there are even times when we will used the definition, directly. We will see an example of this below. Lets now work a couple of examples.
r r r r F gdS = F gn dS ,
S S
Example 1 Evaluate
2 2
paraboloid y = x + z , 0 y 1 and the disk x 2 + z 2 1 at y = 1 . Assume that S has positive orientation. Solution Okay, first lets notice that the disk is really nothing more than the cap on the paraboloid. This means that we have a closed surface. This is important because weve been told that the surface has a positive orientation and by convention this means that all the unit normal vectors will need to point outwards from the region enclosed by S. Lets first get a sketch of S so we can get a feel for what is going on and in which direction we will need to unit normal vectors to point.
22
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
As noted in the sketch we will denote the paraboloid by S1 and the disk by S2 . Also note that in order for unit normal vectors on the paraboloid to point away from the region they will all need to point generally in the negative y direction. On the other hand, unit normal vectors on the disk will need to point in the positive y direction in order to point away from the region. Since S is composed of the two surfaces well need to do the surface integral on each and then add the results to get the overall surface integral. Lets start with the paraboloid. In this case we have the surface in the form y = g ( x, z ) so we will need to derive the correct formula since the one given initially wasnt for this kind of function. This is easy enough to do however. First define,
f ( x, y , z ) = y - g ( x , z ) = y - x 2 - z 2 f = -2 x,1, -2 z
Now, the y component of the gradient is positive and so this vector will generally point in the positive y direction. However, as noted above we need the normal vector point in the negative y direction to make sure that it will be pointing away from the enclosed region. This means that we will need to use
2 x, -1, 2 z r -f n= = -f f
Lets note a couple of things here before we proceed. We dont really need to divide this by the magnitude of the gradient since this will just cancel out once we actually do the integral. So, because of this we didnt bother computing it. Also, the dropping of the minus sign is not a typo. When we compute the magnitude we are going to square each of the components and so the minus sign will drop out.
S1 : The Paraboloid
Okay, here is the surface integral in this case.
2007 Paul Dawkins 23 http://tutorial.math.lamar.edu/terms.aspx
Calculus III
F gdS = ( y j - z k )g
S1 D
r 2 x, -1, 2 z f
f dA
= - y - 2 z 2 dA = - ( x 2 + z 2 ) - 2 z 2 dA
D D
= - x 2 + 3z 2 dA
D
Dont forget that we need to plug in the equation of the surface for y before we actually compute the integral. In this case D is the disk of radius 1 in the xz-plane and so it makes sense to use polar coordinates to complete this integral. Here are polar coordinates for this region.
x = r cos q 0 q 2p
z = r sin q 0 r 1
Note that we kept the x conversion formula the same as the one we are used to using for x and let z be the formula that used the sine. We could have done it any order, however in this way we are at least working with one of them as we are used to working with. Here is the evaluation of this integral.
S1 D
r r F gdS = - x 2 + 3z 2 dA = - 0 = - 0
2p
(r
1 0 1 0
cos 2 q + 3r 2 sin 2 q ) r dr dq
2
2p
( cos
q + 3sin 2 q ) r 3 dr dq
1
2p
We can now do the surface integral on the disk (cap on the paraboloid). This one is actually fairly easy to do and in fact we can use the definition of the surface integral directly. First lets notice that the disk is really just the portion of the plane y = 1 that is in front of the disk of radius 1 in the xz-plane. Now we want the unit normal vector to point away from the enclosed region and since it must
2007 Paul Dawkins 24 http://tutorial.math.lamar.edu/terms.aspx
Calculus III
also be orthogonal to the plane y = 1 then it must point in a direction that is parallel to the y-axis, but we already have a unit vector that does this. Namely,
r r n= j
the standard unit basis vector. It also points in the correct direction for us to use. Because we have the vector field and the normal vector we can plug directly into the definition of the surface integral to get,
r r r r r F gdS = y j - z k g( j ) dS = y dS
S2
S2
S2
At this point we need to plug in for y (since S2 is a portion of the plane y = 1 we do know what it is) and well also need the square root this time when we convert the surface integral over to a double integral. In this case since we are using the definition directly we wont get the canceling of the square root that we saw with the first portion. To get the square root well need to acknowledge that
y = 1 = g ( x, z )
2
( gx )
The surface integral is then,
S2 S2
+1+ ( gz )
r r F gdS = y dS = 1 0 + 1 + 0 dA = dA
D D
At this point we can acknowledge that D is a disk of radius 1 and this double integral is nothing more than the double integral that will give the area of the region D so there is no reason to compute the integral. Here is the value of the surface integral.
r r F gdS = p
S2
Finally, to finish this off we just need to add the two parts up. Here is the surface integral that we were actually asked to compute.
Example 2 Evaluate
2 2 2
x + y + z = 9 and the disk x 2 + y 2 9 in the plane z = 0 . Assume that S has the positive
orientation. Solution So, as with the previous problem we have a closed surface and since we are also told that the surface has a positive orientation all the unit normal vectors must point away from the enclosed region. To help us visualize this here is a sketch of the surface.
25
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
We will call S1 the hemisphere and S2 will be the bottom of the hemisphere (which isnt shown on the sketch). Now, in order for the unit normal vectors on the sphere to point away from enclosed region they will all need to have a positive z component. Remember that the vector must be normal to the surface and if there is a positive z component and the vector is normal it will have to be pointing away from the enclosed region. On the other hand, the unit normal on the bottom of the disk must point in the negative z direction in order to point away from the enclosed region.
S1 : The Sphere
Lets do the surface integral on S1 first. In this case since the surface is a sphere we will need to use the parametric representation of the surface. This is,
Since we are working on the hemisphere here are the limits on the parameters that well need to use.
p 0 q 2p 0 j 2 r r Next, we need to determine rq rj . Here are the two individual vectors and the cross product. r r r rq (q , j ) = -3sin j sin q i + 3sin j cos q j r r r r rj (q , j ) = 3cos j cos q i + 3cos j sin q j - 3sin j k
r i r r rq rj = -3sin j sin q 3cos j cos q
r r j k 3sin j cos q 0 3cos j sin q -3sin j r r r r = -9sin 2 j cos q i - 9sin j cos j sin 2 q k - 9sin 2 j sin q j - 9sin j cos j cos 2 q k r r r = -9sin 2 j cos q i - 9sin 2 j sin q j - 9sin j cos j ( sin 2 q + cos 2 q ) k r r r = -9sin 2 j cos q i - 9sin 2 j sin q j - 9sin j cos j k
26
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
Note that we wont need the magnitude of the cross product since that will cancel out once we start doing the integral. Notice that for the range of j that weve got both sine and cosine are positive and so this vector will have a negative z component and as we noted above in order for this to point away from the enclosed area we will need the z component to be positive. Therefore we will need to use the following vector for the unit normal vector.
Again, we will drop the magnitude once we get to actually doing the integral since it will just cancel in the integral.
r r r r r F ( r (q , j ) ) = 3sin j cos q i + 3sin j sin q j + 81cos 4 j k r Remember that in this evaluation we are just plugging in the x component of r (q , j ) into the
vector field etc. We also may as well get the dot product out of the way that we know we are going to need.
r r r r F ( r (q , j ) )g( rq rj ) = 27 sin 3 j cos 2 q + 27sin 3 j sin 2 q + 729sin j cos5 j = 27 sin 3 j + 729sin j cos5 j
r r r F gdS = F g S1
D
r r ru rv r r rq rj
r r rq rj dA
= 0 = 0
2p
p 2 0 p 2 0
2p
2p
Calculus III
This also means that we can use the definition of the surface integral here with
r r n = -k
We need the negative since it must point away from the enclosed region. The surface integral in this case is,
F gdS = ( x i + y j + z k )g( -k ) dS
4 S2 S2
= - z 4 dS
S2
Remember, however, that we are in the plane given by z = 0 and so the surface integral becomes,
r r F gdS = - z 4 dS = 0 dS = 0
S2 S2 S2
The last step is to then add the two pieces up. Here is surface integral that we were asked to look at.
We will leave this section with a quick interpretation of a surface integral over a vector field. If r v is the velocity field of a fluid then the surface integral
r r v gdS
S
represents the volume of fluid flowing through S per time unit (i.e. per second, per minute, or whatever time unit you are using).
28
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
Stokes Theorem
In this section we are going to take a look at a theorem that is a higher dimensional version of Greens Theorem. In Greens Theorem we related a line integral to a double integral over some region. In this section we are going to relate a line integral to a surface integral. However, before we give the theorem we first need to define the curve that were going to use in the line integral. Lets start off with the following surface with the indicated orientation.
Around the edge of this surface we have a curve C. This curve is called the boundary curve. The orientation of the surface S will induce the positive orientation of C. To get the positive orientation of C think of yourself as walking along the curve. While you are walking along the curve if your head is pointing in the same direction as the unit normal vectors while the surface is on the left then you are walking in the positive direction on C. Now that we have this curve definition out of the way we can give Stokes Theorem. Stokes Theorem Let S be an oriented smooth surface that is bounded by a simple, closed, smooth boundary curve r C with positive orientation. Also let F be a vector field then,
C
F g d r = curl F g dS
S
In this theorem note that the surface S can actually be any surface so long as its boundary curve is given by C. This is something that can be used to our advantage to simplify the surface integral on occasion. Lets take a look at a couple of examples.
29
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
curl F g dS
S
where F = z 2 i - 3 xy j + x 3 y 3 k
and S is the part of z = 5 - x - y above the plane z = 1 . Assume that S is oriented upwards. Solution Lets start this off with a sketch of the surface.
In this case the boundary curve C will be where the surface intersects the plane z = 1 and so will be the curve
1 = 5 - x2 - y 2 x2 + y2 = 4 at z = 1
So, the boundary curve will be the circle of radius 2 that is in the plane z = 1 . The parameterization of this curve is,
r r r r r ( t ) = 2 cos t i + 2sin t j + k , 0 t 2p
The first two components give the circle and the third component makes sure that it is in the plane z = 1. Using Stokes Theorem we can write the surface integral as the following line integral.
r r r r 2p r r r curl F g dS = F g d r = F ( r ( t ) ) g r ( t ) dt
S C 0
So, it looks like we need a couple of quantities before we do this integral. Lets first get the vector field evaluated on the curve. Remember that this is simply plugging the components of the parameterization into the vector field.
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
Next, we need the derivative of the parameterization and the dot product of this and the vector field.
curl F g dS =
= ( 2cos t + 8cos3 t ) =0
r i r curl F = x z2
r j y y2
r k r r r = 2 z j - j = ( 2 z - 1) j z x
Now, all we have is the boundary curve for the surface that well need to use in the surface integral. However, as noted above all we need is any surface that has this as its boundary curve. So, lets use the following plane with upwards orientation for the surface.
31
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
Since the plane is oriented upwards this induces the positive direction on C as shown. The equation of this plane is,
x + y + z =1
z = g ( x, y ) = 1 - x - y
Now, lets use Stokes Theorem and get the surface integral set up.
C
r r r r F g d r = curl F g dS
S S
r r = ( 2 z - 1) j gdS r f = ( 2 z - 1) j g f D f dA
Okay, we now need to find a couple of quantities. First lets get the gradient. Recall that this comes from the function of the surface.
f ( x, y , z ) = z - g ( x, y ) = z - 1 + x + y r r r f = i + j + k
Note as well that this also points upwards and so we have the correct direction. Now, D is the region in the xy-plane shown below,
We get the equation of the line by plugging in z = 0 into the equation of the plane. So based on this the ranges that define D are,
0 x 1
0 y -x +1
r r r r r r F g d r = ( 2 z - 1) j g i + j + k dA
1 0
- x +1 0
2 (1 - x - y ) - 1 dy dx
Dont forget to plug in for z since we are doing the surface integral on the plane. Finishing this out gives,
32
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
r r 1 - x +1 F g d r = 1 - 2 x - 2 y dy dx
0 0
= ( y - 2 xy - y 2 )
1 0
- x +1 0
dx
= x 2 - x dx
0
1 1 = x3 - x 2 2 0 3 1 =6
In both of these examples we were able to take an integral that would have been somewhat unpleasant to deal with and by the use of Stokes Theorem we were able to convert it into an integral that wasnt too bad.
33
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
Divergence Theorem
In this section we are going to relate surface integrals to triple integrals. We will do this with the Divergence Theorem. Divergence Theorem Let E be a simple solid region and S is the boundary surface of E with positive orientation. Let r F be a vector field whose components have continuous first order partial derivatives. Then,
r r r F gdS = div F dV
S E
r r r r r r F gdS where F = xy i - 1 y 2 j + z k 2
S
The region E for the triple integral is then the region enclosed by these surfaces. Note that cylindrical coordinates would be a perfect coordinate system for this region. If we do that here are the limits for the ranges.
0 z 4 - 3r 2 0 r 1 0 q 2p
Well also need the divergence of the vector field so lets get that.
r div F = y - y + 1 = 1
34
http://tutorial.math.lamar.edu/terms.aspx
Calculus III
r r r F gdS = div F dV
S E 1 4-3r 2 = r dz dr dq 0 0 0 2p
= 0
2p
1 0
4r - 3r 3 dr dq
1
3 = 2r 2 - r 4 dq 4 0 0
2p 5 = dq 0 4 5 = p 2
2p
35
http://tutorial.math.lamar.edu/terms.aspx