16.8: Stokes's Theorem
( \newcommand{\kernel}{\mathrm{null}\,}\)
Recall that one version of Green's Theorem (see equation 16.5.1) is
∫∂DF⋅dr=∬D(∇×F)⋅kdA.
Here D is a region in the x-y plane and k is a unit normal to D at every point. If D is instead an orientable surface in space, there is an obvious way to alter this equation, and it turns out still to be true:
Provided that the quantities involved are sufficiently nice, and in particular if D is orientable,
∫∂DF⋅dr=∬D(∇×F)⋅NdS,
if ∂D is oriented counter-clockwise relative to N.
Note how little has changed: k becomes N, a unit normal to the surface, and dA becomes dS, since this is now a general surface integral. The phrase "counter-clockwise relative to N" means roughly that if we take the direction of N to be "up", then we go around the boundary counter-clockwise when viewed from "above''. In many cases, this description is inadequate. A slightly more complicated but general description is this: imagine standing on the side of the surface considered positive; walk to the boundary and turn left. You are now following the boundary in the correct direction.
Let F=⟨exycosz,x2z,xy⟩ and the surface D be x=√1−y2−z2, oriented in the positive x direction. It quickly becomes apparent that the surface integral in Stokes's Theorem is intractable, so we try the line integral. The boundary of D is the unit circle in the y-z plane, r=⟨0,cosu,sinu⟩, 0≤u≤2π. The integral is
∫2π0⟨exycosz,x2z,xy⟩⋅⟨0,−sinu,cosu⟩du=∫2π00du=0,
because x=0.
Consider the cylinder r=⟨cosu,sinu,v⟩, 0≤u≤2π, 0≤v≤2, oriented outward, and F=⟨y,zx,xy⟩. We compute
∬D∇×F⋅NdS=∫∂DF⋅dr
in two ways.
First, the double integral is
∫2π0∫20⟨0,−sinu,v−1⟩⋅⟨cosu,sinu,0⟩dvdu=∫2π0∫20−sin2udvdu=−2π.
The boundary consists of two parts, the bottom circle ⟨cost,sint,0⟩, with t ranging from 0 to 2π, and ⟨cost,sint,2⟩, with t ranging from 2π to 0. We compute the corresponding integrals and add the results:
∫2π0−sin2tdt+∫02π−sin2t+2cos2t=−π−π=−2π,
as before.
An interesting consequence of Stokes's Theorem is that if D and E are two orientable surfaces with the same boundary, then
∬D(∇×F)⋅NdS=∫∂DF⋅dr=∫∂EF⋅dr=∬E(∇×F)⋅NdS.
Sometimes both of the integrals
∬D(∇×F)⋅NdSand∫∂DF⋅dr
are difficult, but you may be able to find a second surface E so that
∬E(∇×F)⋅NdS
has the same value but is easier to compute.
In example 16.8.2 the line integral was easy to compute. But we might also notice that another surface E with the same boundary is the flat disk y2+z2≤1. The unit normal N for this surface is simply i=⟨1,0,0⟩. We compute the curl:
∇×F=⟨x−x2,−exysinz−y,2xz−xexycosz⟩.
Since x=0 everywhere on the surface,
(∇×F)⋅N=⟨0,−exysinz−y,2xz−xexycosz⟩⋅⟨1,0,0⟩=0,
so the surface integral is
∬E0dS=0,
as before. In this case, of course, it is still somewhat easier to compute the line integral, avoiding $nabla×F entirely.
Let F=⟨−y2,x,z2⟩, and let the curve C be the intersection of the cylinder x2+y2=1 with the plane y+z=2, oriented counter-clockwise when viewed from above. We compute ∫CF⋅dr in two ways.
First we do it directly: a vector function for C is({\bf r}=\langle \cos u,\sin u, 2-\sin u\rangle\), so r′=⟨−sinu,cosu,−cosu⟩, and the integral is then
∫2π0y2sinu+xcosu−z2cosudu=∫2π0sin3u+cos2u−(2−sinu)2cosudu=π.
To use Stokes's Theorem, we pick a surface with C as the boundary; the simplest such surface is that portion of the plane y+z=2 inside the cylinder. This has vector equation r=⟨vcosu,vsinu,2−vsinu⟩. We compute ru=⟨−vsinu,vcosu,−vcosu⟩, rv=⟨cosu,sinu,−sinu⟩, and ru×rv=⟨0,−v,−v⟩. To match the orientation of C we need to use the normal ⟨0,v,v⟩. The curl of F is ⟨0,0,1+2y⟩=⟨0,0,1+2vsinu⟩, and the surface integral from Stokes's Theorem is
∫2π0∫10(1+2vsinu)vdvdu=π.
In this case the surface integral was more work to set up, but the resulting integral is somewhat easier.
We can prove here a special case of Stokes's Theorem, which perhaps not too surprisingly uses Green's Theorem.
Suppose the surface D of interest can be expressed in the form z=g(x,y), and let F=⟨P,Q,R⟩. Using the vector function r=⟨x,y,g(x,y)⟩ for the surface we get the surface integral
∬D∇×F⋅dS=∬E⟨Ry−Qz,Pz−Rx,Qx−Py⟩⋅⟨−gx,−gy,1⟩dA=∬E−Rygx+Qzgx−Pzgy+Rxgy+Qx−PydA.
Here E is the region in the x-y plane directly below the surfaceD.
For the line integral, we need a vector function for ∂D. If ⟨x(t),y(t)⟩ is a vector function for ∂E then we may use r(t)=⟨x(t),y(t),g(x(t),y(t))⟩ to represent ∂D. Then
∫∂DF⋅dr=∫baPdxdt+Qdydt+Rdzdtdt=∫baPdxdt+Qdydt+R(∂z∂xdxdt+∂z∂ydydt)dt.
using the chain rule for dz/dt. Now we continue to manipulate this:
∫baPdxdt+Qdydt+R(∂z∂xdxdt+∂z∂ydydt)dt=∫ba[(P+R∂z∂x)dxdt+(Q+R∂z∂y)dydt]dt=∫∂E(P+R∂z∂x)dx+(Q+R∂z∂y)dy,
which now looks just like the line integral of Green's Theorem, except that the functions P and Q of Green's Theorem have been replaced by the more complicated P+R(∂z/∂x) and Q+R(∂z/∂y). We can apply Green's Theorem to get
∫∂E(P+R∂z∂x)dx+(Q+R∂z∂y)dy=∬E∂∂x(Q+R∂z∂y)−∂∂y(P+R∂z∂x)dA.
Now we can use the chain rule again to evaluate the derivatives inside this integral, and it becomes
∬EQx+Qzgx+Rxgy+Rzgxgy+Rgyx−(Py+Pzgy+Rygx+Rzgygx+Rgxy)dA=∬EQx+Qzgx+Rxgy−Py−Pzgy−RygxdA,
which is the same as the expression we obtained for the surface integral.
◻
Contributors
Integrated by Justin Marshall.