# 15.7 The Divergence Theorem and Stokes’ Theorem

## The Divergence Theorem

Theorem 15.4.2 gives the Divergence Theorem in the plane, which states that the flux of a vector field across a closed curve equals the sum of the divergences over the region enclosed by the curve. Recall that the flux was measured via a line integral, and the sum of the divergences was measured through a double integral.

We now consider the three-dimensional version of the Divergence Theorem. It states, in words, that the flux across a closed surface equals the sum of the divergences over the domain enclosed by the surface. Since we are in space (versus the plane), we measure flux via a surface integral, and the sums of divergences will be measured through a triple integral.

margin: Note: the term “outer unit normal vector” used in Theorem 15.7.1 means $\vec{n}$ points to the outside of $\mathcal{S}$. Λ
###### Theorem 15.7.1 The Divergence Theorem (in space)

Let $D$ be a closed domain in space whose boundary is an orientable, piecewise smooth surface $\mathcal{S}$ with outer unit normal vector $\vec{n}$, and let $\vec{F}$ be a vector field whose components are differentiable on $D$. Then

 $\iint_{\mathcal{S}}\vec{F}\cdot\vec{n}\operatorname{d}\!S=\iiint_{D}% \operatorname{div}\vec{F}\operatorname{d}\!V.$
margin: Figure 15.7.1: The surfaces used in Example 15.7.1. Λ
###### Example 15.7.1 Using the Divergence Theorem in space

Let $D$ be the domain in space bounded by the planes $z=0$ and $z=2x$, along with the cylinder $x=1-y^{2}$, as graphed in Figure 15.7.1, let $\mathcal{S}$ be the boundary of $D$, and let $\vec{F}=\left\langle x+y,y^{2},2z\right\rangle$.

Verify the Divergence Theorem by finding the total outward flux of $\vec{F}$ across $\mathcal{S}$, and show this is equal to $\iiint_{D}\operatorname{div}\vec{F}\operatorname{d}\!V$.

SolutionThe surface $\mathcal{S}$ is piecewise smooth, comprising surfaces $\mathcal{S}_{1}$, which is part of the plane $z=2x$, surface $\mathcal{S}_{2}$, which is part of the cylinder $x=1-y^{2}$, and surface $\mathcal{S}_{3}$, which is part of the plane $z=0$. To find the total outward flux across $\mathcal{S}$, we need to compute the outward flux across each of these three surfaces.

We leave it to the reader to confirm that surfaces $\mathcal{S}_{1}$, $\mathcal{S}_{2}$ and $\mathcal{S}_{3}$ can be parameterized by $\vec{r}_{1}$, $\vec{r}_{2}$ and $\vec{r}_{3}$ respectively as

 $\displaystyle\vec{r}_{1}(u,v)$ $\displaystyle=\left\langle v(1-u^{2}),u,2v(1-u^{2})\right\rangle,$ $\displaystyle\vec{r}_{2}(u,v)$ $\displaystyle=\left\langle(1-u^{2}),u,2v(1-u^{2})\right\rangle,$ $\displaystyle\vec{r}_{3}(u,v)$ $\displaystyle=\left\langle v(1-u^{2}),u,0\right\rangle,$

where $-1\leq u\leq 1$ and $0\leq v\leq 1$ for all three functions.

We compute a unit normal vector $\vec{n}$ for each as $\frac{\vec{r}_{u}\times\vec{r}_{v}}{\left\lVert\vec{r}_{u}\times\vec{r}_{v}% \right\rVert}$, though recall that as we are integrating $\vec{F}\cdot\vec{n}\operatorname{d}\!S$, we actually only use $\vec{r}_{u}\times\vec{r}_{v}$. Finally, in previous flux computations, it did not matter which direction $\vec{n}$ pointed as long as we made note of its direction. When using the Divergence Theorem, we need $\vec{n}$ to point to the outside of the closed surface, so in practice this means we’ll either use $\vec{r}_{u}\times\vec{r}_{v}$ or $\vec{r}_{v}\times\vec{r}_{u}$, depending on which points outside of the closed surface $\mathcal{S}$.

We leave it to the reader to confirm the following cross products and integrations are correct.

For $\mathcal{S}_{1}$, we need to use $\vec{r}_{1v}\times\vec{r}_{1u}=\left\langle 2(u^{2}-1),0,1-u^{2}\right\rangle$. (Note the $z$-component is nonnegative as $u\leq 1$, therefore this vector always points up, meaning to the outside, of $\mathcal{S}$.) The flux across $\mathcal{S}_{1}$ is:

 $\displaystyle\iint_{\mathcal{S}_{1}}$ $\displaystyle\vec{F}\cdot\vec{n}_{1}\operatorname{d}\!S$ $\displaystyle=\int_{0}^{1}\int_{-1}^{1}\vec{F}\bigl{(}\vec{r}_{1}(u,v)\bigr{)}% \cdot\bigl{(}\vec{r}_{1v}\times\vec{r}_{1u}\bigr{)}\operatorname{d}\!u% \operatorname{d}\!v$ $\displaystyle=\int_{0}^{1}\int_{-1}^{1}\left\langle v(1-u^{2})+u,u^{2},4v(1-u^% {2})\right\rangle\cdot\left\langle 2(u^{2}-1),0,1-u^{2}\right\rangle% \operatorname{d}\!u\operatorname{d}\!v$ $\displaystyle=\int_{0}^{1}\int_{-1}^{1}\bigl{(}2u^{4}v+2u^{3}-4u^{2}v-2u+2v% \bigr{)}\operatorname{d}\!u\operatorname{d}\!v$ $\displaystyle=\frac{16}{15}.$

For $\mathcal{S}_{2}$, we use $\vec{r}_{2u}\times\vec{r}_{2v}=\left\langle 2(1-u^{2}),4u(1-u^{2}),0\right\rangle$. (Note that the first is always nonnegative, meaning this vector points outside $\mathcal{S}$.) The flux across $\mathcal{S}_{2}$ is:

 $\displaystyle\iint_{\mathcal{S}_{2}}\vec{F}\cdot\vec{n}_{2}\operatorname{d}\!S$ $\displaystyle=\int_{0}^{1}\int_{-1}^{1}\vec{F}\bigl{(}\vec{r}_{2}(u,v)\bigr{)}% \cdot\bigl{(}\vec{r}_{2u}\times\vec{r}_{2v}\bigr{)}\operatorname{d}\!u% \operatorname{d}\!v$ $\displaystyle=\int_{0}^{1}\int_{-1}^{1}\left\langle 1-u^{2}+u,u^{2},4v(1-u^{2}% )\right\rangle$ $\displaystyle\hphantom{=\int_{0}^{1}\int_{-1}^{1}}\cdot\left\langle 2(1-u^{2})% ,4u(1-u^{2}),0\right\rangle\operatorname{d}\!u\operatorname{d}\!v$ $\displaystyle=\int_{0}^{1}\int_{-1}^{1}\bigl{(}4u^{5}-2u^{4}-2u^{3}+4u^{2}-2u-% 2\bigr{)}\operatorname{d}\!u\operatorname{d}\!v$ $\displaystyle=\frac{32}{15}.$

For $\mathcal{S}_{3}$, we use $\vec{r}_{3u}\times\vec{r}_{3v}=\left\langle 0,0,u^{2}-1\right\rangle$. (Note the $z$-component is never positive, meaning this vector points down, outside of $\mathcal{S}$.) The flux across $\mathcal{S}_{3}$ is:

 $\displaystyle\iint_{\mathcal{S}_{3}}\vec{F}\cdot\vec{n}_{3}\operatorname{d}\!S$ $\displaystyle=\int_{0}^{1}\int_{-1}^{1}\vec{F}\bigl{(}\vec{r}_{3}(u,v)\bigr{)}% \cdot\bigl{(}\vec{r}_{3u}\times\vec{r}_{3v}\bigr{)}\operatorname{d}\!u% \operatorname{d}\!v$ $\displaystyle=\int_{0}^{1}\int_{-1}^{1}\left\langle v(1-u^{2})+u,u^{2},0\right% \rangle\cdot\left\langle 0,0,u^{2}-1\right\rangle\operatorname{d}\!u% \operatorname{d}\!v$ $\displaystyle=\int_{0}^{1}\int_{-1}^{1}0\operatorname{d}\!u\operatorname{d}\!v$ $\displaystyle=0.$

Thus the total outward flux, measured by surface integrals across all three component surfaces of $\mathcal{S}$, is $16/15+32/15+0=48/15=16/5=3.2$. We now find the total outward flux by integrating $\operatorname{div}\vec{F}$ over $D$.

Following the steps outlined in Section 14.6, we see the bounds of $x$, $y$ and $z$ can be set as (thinking “surface to surface, curve to curve, point to point”):

 $0\leq z\leq 2x;\quad 0\leq x\leq 1-y^{2};\quad-1\leq y\leq 1.$

With $\operatorname{div}\vec{F}=1+2y+2=2y+3$, we find the total outward flux of $\vec{F}$ over $\mathcal{S}$ as:

 $\text{Flux}=\iiint_{D}\operatorname{div}\vec{F}\operatorname{d}\!V=\int_{-1}^{% 1}\int_{0}^{1-y^{2}}\int_{0}^{2x}\bigl{(}2y+3\bigr{)}\operatorname{d}\!z% \operatorname{d}\!x\operatorname{d}\!y=16/5,$

the same result we obtained previously.

In Example 15.7.1 we see that the total outward flux of a vector field across a closed surface can be found two different ways because of the Divergence Theorem. One computation took far less work to obtain. In that particular case, since $\mathcal{S}$ was comprised of three separate surfaces, it was far simpler to compute one triple integral than three surface integrals (each of which required partial derivatives and a cross product). In practice, if outward flux needs to be measured, one would choose only one method. We will use both methods in this section simply to reinforce the truth of the Divergence Theorem.

We practice again in the following example.

margin: Figure 15.7.2: The surfaces used in Example 15.7.2. Λ
###### Example 15.7.2 Using the Divergence Theorem in space

Let $\mathcal{S}$ be the surface formed by the paraboloid $z=1-x^{2}-y^{2}$, $z\geq 0$, and the unit disk centered at the origin in the $x$-$y$ plane, graphed in Figure 15.7.2, and let $\vec{F}=\left\langle 0,0,z\right\rangle$. (This surface and vector field were used in Example 15.6.3.)

Verify the Divergence Theorem; find the total outward flux across $\mathcal{S}$ and evaluate the triple integral of $\operatorname{div}\vec{F}$, showing that these two quantities are equal.

SolutionWe find the flux across $\mathcal{S}$ first. As $\mathcal{S}$ is piecewise-smooth, we decompose it into smooth components $\mathcal{S}_{1}$, the disk, and $\mathcal{S}_{2}$, the paraboloid, and find the flux across each.

In Example 15.6.3, we found the flux across $\mathcal{S}_{1}$ is 0. We also found that the flux across $\mathcal{S}_{2}$ is $\pi/2$. (In that example, the normal vector had a positive $z$ component hence was an outer normal.) Thus the total outward flux is $0+\pi/2=\pi/2$.

We now compute $\iiint_{D}\operatorname{div}\vec{F}\operatorname{d}\!V$. We can describe $D$ as the domain bounded by (think “surface to surface, curve to curve, point to point”):

 $0\leq z\leq 1-x^{2}-y^{2},\quad-\sqrt{1-x^{2}}\leq y\leq\sqrt{1-x^{2}},\quad-1% \leq x\leq 1.$

This description of $D$ is not very easy to integrate. With polar, we can do better. Let $R$ represent the unit disk, which can be described in polar simply as $r$, where $0\leq r\leq 1$ and $0\leq\theta\leq 2\pi$. With $x=r\cos\theta$ and $y=r\sin\theta$, the surface $\mathcal{S}_{2}$ becomes

 $z=1-x^{2}-y^{2}\Rightarrow 1-(r\cos\theta)^{2}-(r\sin\theta)^{2}\Rightarrow 1-% r^{2}.$

Thus $D$ can be described as the domain bounded by:

 $0\leq z\leq 1-r^{2},\quad 0\leq r\leq 1,\quad 0\leq\theta\leq 2\pi.$

With $\operatorname{div}\vec{F}=1$, we can integrate, recalling that $\operatorname{d}\!V=r\operatorname{d}\!z\operatorname{d}\!r\operatorname{d}\!\theta$:

 $\iiint_{D}\operatorname{div}\vec{F}\operatorname{d}\!V=\int_{0}^{2\pi}\int_{0}% ^{1}\int_{0}^{1-r^{2}}r\operatorname{d}\!z\operatorname{d}\!r\operatorname{d}% \!\theta=\frac{\pi}{2},$

which matches our flux computation above.

###### Example 15.7.3 A “paradox” of the Divergence Theorem and Gauss’s Law

The magnitude of many physical quantities (such as light intensity or electromagnetic and gravitational forces) follow an “inverse square law”: the magnitude of the quantity at a point is inversely proportional to the square of the distance to the source of the quantity.

Let a point light source be placed at the origin and let $\vec{F}$ be the vector field which describes the intensity and direction of the emanating light. At a point $(x,y,z)$, the unit vector describing the direction of the light passing through that point is $\left\langle x,y,z\right\rangle/\sqrt{x^{2}+y^{2}+z^{2}}$. As the intensity of light follows the inverse square law, the magnitude of $\vec{F}$ at $(x,y,z)$ is $k/(x^{2}+y^{2}+z^{2})$ for some constant $k$. Taken together,

 $\vec{F}(x,y,z)=\frac{k}{(x^{2}+y^{2}+z^{2})^{3/2}}\left\langle x,y,z\right\rangle.$
margin: Figure 15.7.3: The cube used in Example 15.7.3. Λ

Consider the cube, centered at the origin, with sides of length $2a$ for some $a>0$ (hence corners of the cube lie at $(a,a,a)$, $(-a,-a,-a)$, etc., as shown in Figure 15.7.3). Find the flux across the six faces of the cube and compare this to $\iint_{D}\operatorname{div}\vec{F}\operatorname{d}\!V$.

SolutionLet $\mathcal{S}_{1}$ be the “top” face of the cube, which can be parameterized by $\vec{r}(u,v)=\left\langle u,v,a\right\rangle$ for $-a\leq u\leq a$, $-a\leq v\leq a$. We leave it to the reader to confirm that $\vec{r}_{u}\times\vec{r}_{v}=\left\langle 0,0,1\right\rangle$, which points outside of the cube.

The flux across this face is:

 Flux $\displaystyle=\iint_{\mathcal{S}_{1}}\vec{F}\cdot\vec{n}\operatorname{d}\!S$ $\displaystyle=\int_{-a}^{a}\int_{-a}^{a}\vec{F}\bigl{(}\vec{r}(u,v)\bigr{)}% \cdot\bigl{(}\vec{r}_{u}\times\vec{r}_{v}\bigr{)}\operatorname{d}\!u% \operatorname{d}\!v$ $\displaystyle=\int_{-a}^{a}\int_{-a}^{a}\frac{k\ a}{(u^{2}+v^{2}+a^{2})^{3/2}}% \operatorname{d}\!u\operatorname{d}\!v.$ This double integral is not trivial to compute, requiring multiple trigonometric substitutions. This example is not meant to stress integration techniques, so we leave it to the reader to confirm the result is $\displaystyle=\frac{2k\pi}{3}.$

Note how the result is independent of $a$; no matter the size of the cube, the flux through the top surface is always $2k\pi/3$.

An argument of symmetry shows that the flux through each of the six faces is $2k\pi/3$, thus the total flux through the faces of the cube is $6\times 2k\pi/3=4k\pi$.

It takes a bit of algebra, but we can show that $\operatorname{div}\vec{F}=0$. Thus the Divergence Theorem would seem to imply that the total flux through the faces of the cube should be

 $\text{Flux }=\iiint_{D}\operatorname{div}\vec{F}\operatorname{d}\!V=\iiint_{D}% 0\operatorname{d}\!V=0,$

but clearly this does not match the result from above. What went wrong?

Revisit the statement of the Divergence Theorem. One of the conditions is that the components of $\vec{F}$ must be differentiable on the domain enclosed by the surface. In our case, $\vec{F}$ is not differentiable at the origin — it is not even defined! As $\vec{F}$ does not satisfy the conditions of the Divergence Theorem, it does not apply, and we cannot expect $\iint_{\mathcal{S}}\vec{F}\cdot\vec{n}\operatorname{d}\!A=\iiint_{D}% \operatorname{div}\vec{F}\operatorname{d}\!V$.

Since $\vec{F}$ is differentiable everywhere except the origin, the Divergence Theorem does apply over any domain that does not include the origin. Let $\mathcal{S}_{2}$ be any surface that encloses the cube used before, and let $\hat{D}$ be the domain between the cube and $\mathcal{S}_{2}$; note how $\hat{D}$ does not include the origin and so the Divergence Theorem does apply over this domain. The total outward flux over $\hat{D}$ is thus $\iint_{\hat{D}}\operatorname{div}\vec{F}\operatorname{d}\!V=0$, which means the amount of flux coming out of $\mathcal{S}_{2}$ is the same as the amount of flux coming out of the cube. The conclusion: the flux across any surface enclosing the origin will be $4k\pi$.

This has an important consequence in electrodynamics. Let $q$ be a point charge at the origin. The electric field generated by this point charge is

 $\vec{E}=\frac{q}{4\pi\epsilon_{0}}\frac{\left\langle x,y,z\right\rangle}{(x^{2% }+y^{2}+z^{2})^{3/2}},$

i.e., it is $\vec{F}$ with $k=q/(4\pi\epsilon_{0})$, where $\epsilon_{0}$ is a physical constant (the “permittivity of free space”). Gauss’s Law states that the outward flux of $\vec{E}$ across any surface enclosing the origin is $q/\epsilon_{0}$.

Our interest in the Divergence Theorem is twofold. First, it’s truth alone is interesting: to study the behavior of a vector field across a closed surface, one can examine properties of that field within the surface. Secondly, it offers an alternative way of computing flux. When there are multiple methods of computing a desired quantity, one has power to select the easiest computation as illustrated next.

###### Example 15.7.4 Using the Divergence Theorem to compute flux

Let $\mathcal{S}$ be the cube bounded by the planes $x=\pm 1$, $y=\pm 1$, $z=\pm 1$, and let $\vec{F}=\left\langle x^{2}y,2yz,x^{2}z^{3}\right\rangle$. Compute the outward flux of $\vec{F}$ over $\mathcal{S}$.

SolutionWe compute $\operatorname{div}\vec{F}=2xy+2z+3x^{2}z^{2}$. By the Divergence Theorem, the outward flux is the triple integral over the domain $D$ enclosed by $\mathcal{S}$:

 $\text{Outward flux:}\quad\int_{-1}^{1}\int_{-1}^{1}\int_{-1}^{1}(2xy+2z+3x^{2}% z^{2})\operatorname{d}\!z\operatorname{d}\!y\operatorname{d}\!x=\frac{8}{3}.$

The direct flux computation requires six surface integrals, one for each face of the cube. The Divergence Theorem offers a much more simple computation.

## Stokes’ Theorem

Just as the spatial Divergence Theorem of this section is an extension of the planar Divergence Theorem, Stokes’ Theorem is the spatial extension of Green’s Theorem. Recall that Green’s Theorem states that the circulation of a vector field around a closed curve in the plane is equal to the sum of the curl of the field over the region enclosed by the curve. Stokes’ Theorem effectively makes the same statement: given a closed curve that lies on a surface $\mathcal{S}$, the circulation of a vector field around that curve is the same as the sum of “the curl of the field” across the enclosed surface. We use quotes around “the curl of the field” to signify that this statement is not quite correct, as we do not sum $\operatorname{curl}\vec{F}$, but $\operatorname{curl}\vec{F}\cdot\vec{n}$, where $\vec{n}$ is a unit vector normal to $\mathcal{S}$. That is, we sum the portion of $\operatorname{curl}\vec{F}$ that is orthogonal to $\mathcal{S}$ at a point.

Green’s Theorem dictated that the curve was to be traversed counterclockwise when measuring circulation. Stokes’ Theorem will follow a right hand rule: when the thumb of one’s right hand points in the direction of $\vec{n}$, the path $C$ will be traversed in the direction of the curling fingers of the hand (this is equivalent to traversing counterclockwise in the plane).

###### Theorem 15.7.2 Stokes’ Theorem

Let $\mathcal{S}$ be a piecewise smooth, orientable surface whose boundary is a piecewise smooth curve $C$, let $\vec{n}$ be a unit vector normal to $\mathcal{S}$, let $C$ be traversed with respect to $\vec{n}$ according to the right hand rule, and let the components of $\vec{F}$ have continuous first partial derivatives over $\mathcal{S}$. Then

 $\oint_{C}\vec{F}\cdot\operatorname{d}\!\vec{r}=\iint_{\mathcal{S}}(% \operatorname{curl}\vec{F})\cdot\vec{n}\operatorname{d}\!S.$

In general, the best approach to evaluating the surface integral in Stokes’ Theorem is to parameterize the surface $\mathcal{S}$ with a function $\vec{r}(u,v)$. We can find a unit normal vector $\vec{n}$ as

 $\vec{n}=\frac{\vec{r}_{u}\times\vec{r}_{v}}{\left\lVert\vec{r}_{u}\times\vec{r% }_{v}\right\rVert}.$

Since $\operatorname{d}\!S=\left\lVert\vec{r}_{u}\times\vec{r}_{v}\right\rVert% \operatorname{d}\!A$, the surface integral in practice is evaluated as

 $\iint_{\mathcal{S}}(\operatorname{curl}\vec{F})\cdot(\vec{r}_{u}\times\vec{r}_% {v})\operatorname{d}\!A,$

where $\vec{r}_{u}\times\vec{r}_{v}$ may be replaced by $\vec{r}_{v}\times\vec{r}_{u}$ to properly match the direction of this vector with the orientation of the parameterization of $C$.

margin: Figure 15.7.4: As given in Example 15.7.5, the surface $\mathcal{S}$ is the portion of the plane bounded by the curve. Λ
###### Example 15.7.5 Verifying Stokes’ Theorem

Considering the planar surface $f(x,y)=7-2x-2y$, let $C$ be the curve in space that lies on this surface above the circle of radius 1 and centered at $(1,1)$ in the $x$-$y$ plane, let $\mathcal{S}$ be the planar region enclosed by $C$, as illustrated in Figure 15.7.4, and let $\vec{F}=\left\langle x+y,2y,y^{2}\right\rangle$. Verify Stokes’ Theorem by showing $\oint_{C}\vec{F}\cdot\operatorname{d}\!\vec{r}=\iint_{\mathcal{S}}(% \operatorname{curl}\vec{F})\cdot\vec{n}\operatorname{d}\!S$.

SolutionWe begin by parameterizing $C$ and then find the circulation. A unit circle centered at $(1,1)$ can be parameterized with $x=\cos t+1$, $y=\sin t+1$ on $0\leq t\leq 2\pi$; to put this curve on the surface $f$, make the $z$ component equal $f(x,y)$: $z=7-2(\cos t+1)-2(\sin t+1)=3-2\cos t-2\sin t$. All together, we parameterize $C$ with $\vec{r}(t)=\left\langle\cos t+1,\sin t+1,3-2\cos t-2\sin t\right\rangle$.

The circulation of $\vec{F}$ around $C$ is

 $\displaystyle\oint_{C}\vec{F}\cdot\operatorname{d}\!\vec{r}$ $\displaystyle=\int_{0}^{2\pi}\vec{F}\bigl{(}\vec{r}(t)\bigr{)}\cdot\vec{r}% \hskip 1.25pt^{\prime}(t)\operatorname{d}\!t$ $\displaystyle=\int_{0}^{2\pi}\bigl{(}2\sin^{3}t-2\cos t\sin^{2}t+3\sin^{2}t-3% \cos t\sin t\bigr{)}\operatorname{d}\!t$ $\displaystyle=3\pi.$

We now parameterize $\mathcal{S}$. (We reuse the letter “r” for our surface as this is our custom.) Based on the parameterization of $C$ above, we describe $\mathcal{S}$ with $\vec{r}(u,v)=\left\langle v\cos u+1,v\sin u+1,3-2v\cos u-2v\sin u\right\rangle$, where $0\leq u\leq 2\pi$ and $0\leq v\leq 1$.

We leave it to the reader to confirm that $\vec{r}_{u}\times\vec{r}_{v}=\left\langle 2v,2v,v\right\rangle$. As $0\leq v\leq 1$, this vector always has a non-negative $z$-component, which the right-hand rule requires given the orientation of $C$ used above. We also leave it to the reader to confirm $\operatorname{curl}\vec{F}=\left\langle 2y,0,-1\right\rangle$.

The surface integral of Stokes’ Theorem is thus

 $\displaystyle\iint_{\mathcal{S}}(\operatorname{curl}\vec{F})\cdot\vec{n}% \operatorname{d}\!S$ $\displaystyle=\iint_{\mathcal{S}}(\operatorname{curl}\vec{F})\cdot(\vec{r}_{u}% \times\vec{r}_{v})\operatorname{d}\!A$ $\displaystyle=\int_{0}^{1}\int_{0}^{2\pi}\left\langle 2v\sin u+2,0,-1\right% \rangle\cdot\left\langle 2v,2v,v\right\rangle\operatorname{d}\!u\operatorname{% d}\!v$ $\displaystyle=3\pi,$

which matches our previous result.

One of the interesting results of Stokes’ Theorem is that if two surfaces $\mathcal{S}_{1}$ and $\mathcal{S}_{2}$ share the same boundary, then $\iint_{\mathcal{S}_{1}}(\operatorname{curl}\vec{F})\cdot\vec{n}\operatorname{d% }\!S=\iint_{\mathcal{S}_{2}}(\operatorname{curl}\vec{F})\cdot\vec{n}% \operatorname{d}\!S$. That is, the value of these two surface integrals is somehow independent of the interior of the surface. We demonstrate this principle in the next example.

margin: (a) (b) Figure 15.7.5: As given in Example 15.7.6, the surface $\mathcal{S}$ is the portion of the plane bounded by the curve. Λ
###### Example 15.7.6 Stokes’ Theorem and surfaces that share a boundary

Let $C$ be the curve given in Example 15.7.5 and note that it lies on the surface $z=6-x^{2}-y^{2}$. Let $\mathcal{S}$ be the region of this surface bounded by $C$, and let $\vec{F}=\left\langle x+y,2y,y^{2}\right\rangle$ as in the previous example. Compute $\iint_{\mathcal{S}}(\operatorname{curl}\vec{F})\cdot\vec{n}\operatorname{d}\!S$ to show it equals the result found in the previous example.

SolutionWe begin by demonstrating that $C$ lies on the surface $z=6-x^{2}-y^{2}$. We can parameterize the $x$ and $y$ components of $C$ with $x=\cos t+1$, $y=\sin t+1$ as before. Lifting these components to the surface $f$ gives the $z$ component as $z=6-x^{2}-y^{2}=6-(\cos t+1)^{2}-(\sin t+1)^{2}=3-2\cos t-2\sin t$, which is the same $z$ component as found in Example 15.7.5. Thus the curve $C$ lies on the surface $z=6-x^{2}-y^{2}$, as illustrated in Figure 15.7.5.

Since $C$ and $\vec{F}$ are the same as in the previous example, we already know that $\oint_{C}\vec{F}\cdot\ d\vec{r}=3\pi$. We confirm that this is also the value of $\iint_{\mathcal{S}}(\operatorname{curl}\vec{F})\cdot\vec{n}\operatorname{d}\!S$.

We parameterize $\mathcal{S}$ with

 $\vec{r}(u,v)=\left\langle v\cos u+1,v\sin u+1,6-(v\cos u+1)^{2}-(v\sin u+1)^{2% }\right\rangle,$

where $0\leq u\leq 2\pi$ and $0\leq v\leq 1$, and leave it to the reader to confirm that

 $\vec{r}_{u}\times\vec{r}_{v}=\left\langle 2v\bigl{(}v\cos u+1\bigr{)},2v\bigl{% (}v\sin u+1\bigr{)},v\right\rangle,$

which also conforms to the right-hand rule with regard to the orientation of $C$. With $\operatorname{curl}\vec{F}=\left\langle 2y,0,-1\right\rangle$ as before, we have

 $\iint_{\mathcal{S}}(\operatorname{curl}\vec{F})\cdot\vec{n}\operatorname{d}\!S% \\ =\int_{0}^{1}\int_{0}^{2\pi}\left\langle 2v\sin u+2,0,-1\right\rangle\cdot% \left\langle 2v\bigl{(}v\cos u+1\bigr{)},2v\bigl{(}v\sin u+1\bigr{)},v\right% \rangle\operatorname{d}\!u\operatorname{d}\!v\\ =3\pi.$

Even though the surfaces used in this example and in Example 15.7.5 are very different, because they share the same boundary, Stokes’ Theorem guarantees they have equal “sum of curls” across their respective surfaces.

## A Common Thread of Calculus

We have threefold interest in each of the major theorems of this chapter: the Fundamental Theorem of Line Integrals, Green’s, Stokes’ and the Divergence Theorems. First, we find the beauty of their truth interesting. Second, each provides two methods of computing a desired quantity, sometimes offering a simpler method of computation.

There is yet one more reason of interest in the major theorems of this chapter. These important theorems also all share an important principle with the Fundamental Theorem of Calculus, introduced in Chapter 5.

Revisit this fundamental theorem, adopting the notation used heavily in this chapter. Let $I$ be the interval $[a,b]$ and let $y=F(x)$ be differentiable on $I$, with $F\,^{\prime}(x)=f(x)$. The Fundamental Theorem of Calculus states that

 $\int_{I}f(x)\operatorname{d}\!x=F(b)-F(a).$

That is, the sum of the rates of change of a function $F$ over an interval $I$ can also be calculated with a certain sum of $F$ itself on the boundary of $I$ (in this case, at the points $x=a$ and $x=b$).

Each of the named theorems above can be expressed in similar terms. Consider the Fundamental Theorem of Line Integrals: given a function $z=f(x,y)$, the gradient $\nabla f$ is a type of rate of change of $f$. Given a curve $C$ with initial and terminal points $A$ and $B$, respectively, this fundamental theorem states that

 $\int_{C}\nabla f\operatorname{d}\!s=f(B)-f(A),$

where again the sum of a rate of change of $f$ along a curve $C$ can also be evaluated by a certain sum of $f$ at the boundary of $C$ (i.e., the points $A$ and $B$).

Green’s Theorem is essentially a special case of Stokes’ Theorem, so we consider just Stokes’ Theorem here. Recalling that the curl of a vector field $\vec{F}$ is a measure of a rate of change of $\vec{F}$, Stokes’ Theorem states that over a surface $\mathcal{S}$ bounded by a closed curve $C$,

 $\iint_{\mathcal{S}}\bigl{(}\operatorname{curl}\vec{F}\bigr{)}\cdot\vec{n}% \operatorname{d}\!S=\oint_{C}\vec{F}\cdot\operatorname{d}\!\vec{r},$

i.e., the sum of a rate of change of $\vec{F}$ can be calculated with a certain sum of $\vec{F}$ itself over the boundary of $\mathcal{S}$. In this case, the latter sum is also an infinite sum, requiring an integral.

Finally, the Divergence Theorems state that the sum of divergences of a vector field (another measure of a rate of change of $\vec{F}$) over a region can also be computed with a certain sum of $\vec{F}$ over the boundary of that region. When the region is planar, the latter sum of $\vec{F}$ is an integral; when the region is spatial, the latter sum of $\vec{F}$ is a double integral.

The common thread among these theorems: the sum of a rate of change of a function over a region can be computed as another sum of the function itself on the boundary of the region. While very general, this is a very powerful and important statement.

In fact, mathematicians occasionally use the symbol $\partial$ to indicate the boundary of something (so the boundary of $D$ would be $\partial D$). This is the same symbol as a partial derivative, so you’ll have to look at the context to figure out which definition is being used. Having done this, we can express these fundamental theorems in a similar language in Key Idea 15.7.1.

###### Key Idea 15.7.1 Fundamental Theorems Relating Integrals and Domains
 Theorem Equation Orientation Fundamental Theorem of Calculus $\displaystyle\left.f\right|_{a}^{b}=\int_{a}^{b}f\,^{\prime}(x)\operatorname{d% }\!x$ Fundamental Theorem of Gradient Fields $\displaystyle\left.\phi\right|_{P}^{Q}=\int_{C}\nabla\phi\cdot\operatorname{d}\!s$ $C$ goes from $P$ to $Q$ Green’s Theorem $\displaystyle\oint_{\partial R}\vec{F}\cdot\operatorname{d}\!\vec{r}=\iint_{R}% \operatorname{curl}F\operatorname{d}\!A$ $\partial D$ oriented counterclockwise Stokes’ Theorem $\displaystyle\oint_{\partial S}\vec{F}\cdot\operatorname{d}\!\vec{r}=\iint_{S}% \operatorname{curl}\vec{F}\cdot\vec{n}\operatorname{d}\!S$ $\partial S$ oriented with $S$ to the left $\mathbb{R}^{2}$ Divergence Theorem $\displaystyle\iint_{\partial R}\vec{F}\cdot\vec{n}\operatorname{d}\!s=\iint_{R% }\operatorname{div}\vec{F}\operatorname{d}\!A$ $\partial R$ oriented outwards $\mathbb{R}^{3}$ Divergence Theorem $\displaystyle\iint_{\partial D}\vec{F}\cdot\operatorname{d}\!\vec{\sigma}=% \iiint_{D}\operatorname{div}\vec{F}\operatorname{d}\!V$ $\partial W$ oriented outwards

## Exercises 15.7

### Terms and Concepts

1. 1.

What are the differences between the Divergence Theorems of Section 15.4 and this section?

2. 2.

What property of a vector field does the Divergence Theorem relate to flux?

3. 3.

What property of a vector field does Stokes’ Theorem relate to circulation?

4. 4.

Stokes’ Theorem is the spatial version of what other theorem?

### Problems

In Exercises 5–8., a closed surface $\mathcal{S}$ enclosing a domain $D$ and a vector field $\vec{F}$ are given. Verify the Divergence Theorem on $\mathcal{S}$; that is, show $\iint_{\mathcal{S}}\vec{F}\cdot\vec{n}\operatorname{d}\!S=\iiint_{D}% \operatorname{div}\vec{F}\operatorname{d}\!V$.

1. 5.
$\mathcal{S}$ is the surface bounding the domain $D$ enclosed by the plane $z=2-x/2-2y/3$ and the coordinate planes in the first octant; $\vec{F}=\left\langle x^{2},y^{2},x\right\rangle$.
2. 6.
$\mathcal{S}$ is the surface bounding the domain $D$ enclosed by the cylinder $x^{2}+y^{2}=1$ and the planes $z=-3$ and $z=3$; $\vec{F}=\left\langle-x,y,z\right\rangle$.
3. 7.
$\mathcal{S}$ is the surface bounding the domain $D$ enclosed by $z=xy(3-x)(3-y)$ and the plane $z=0$; $\vec{F}=\left\langle 3x,4y,5z+1\right\rangle$.
4. 8.
$\mathcal{S}$ is the surface composed of $\mathcal{S}_{1}$, the paraboloid $z=4-x^{2}-y^{2}$ for $z\geq 0$, and $\mathcal{S}_{2}$, the disk of radius 2 centered at the origin; $\vec{F}=\left\langle x,y,z^{2}\right\rangle$.

In Exercises 9–14., a closed curve $C$ that is the boundary of a surface $\mathcal{S}$ is given along with a vector field $\vec{F}$. Verify Stokes’ Theorem on $C$; that is, show $\oint_{C}\vec{F}\cdot\operatorname{d}\!\vec{r}=\iint_{\mathcal{S}}\bigl{(}% \operatorname{curl}\vec{F}\,\bigr{)}\cdot\vec{n}\operatorname{d}\!S$.

1. 9.

$C$ is the unit circle in the $xy$ plane and $\mathcal{S}$ upper unit hemisphere; $\vec{F}=2y\,\vec{\imath}-x\,\vec{\jmath}+z\,\vec{k}$

2. 10.

$C$ is the curve parameterized by $\vec{r}(t)=\left\langle\cos t,\sin t,1\right\rangle$ and $\mathcal{S}$ is the portion of $z=x^{2}+y^{2}$ enclosed by $C$; $\vec{F}=xy\,\vec{\imath}+xz\,\vec{\jmath}+yz\,\vec{k}$

3. 11.

$C$ is the curve parameterized by $\vec{r}(t)=\left\langle\cos t,\sin t,1\right\rangle$ and $\mathcal{S}$ is the portion of $z=x^{2}+y^{2}$ enclosed by $C$; $\vec{F}=\left\langle z,-x,y\right\rangle$.

4. 12.

$C$ is the curve parameterized by $\vec{r}(t)=\left\langle\cos t,\sin t,e^{-1}\right\rangle$ and $\mathcal{S}$ is the portion of $z=e^{-x^{2}-y^{2}}$ enclosed by $C$; $\vec{F}=\left\langle-y,x,1\right\rangle$.

5. 13.

$C$ is the curve that follows the triangle with vertices at $(0,0,2)$, $(4,0,0)$ and $(0,3,0)$, traversing the the vertices in that order and returning to $(0,0,2)$, and $\mathcal{S}$ is the portion of the plane $z=2-x/2-2y/3$ enclosed by $C$; $\vec{F}=\left\langle y,-z,y\right\rangle$.

6. 14.

$C$ is the curve whose $x$ and $y$ coordinates follow the parabola $y=1-x^{2}$ from $x=1$ to $x=-1$, then follow the line from $(-1,0)$ back to $(1,0)$, where the $z$ coordinates of $C$ are determined by $f(x,y)=2x^{2}+y^{2}$, and $\mathcal{S}$ is the portion of $z=2x^{2}+y^{2}$ enclosed by $C$; $\vec{F}=\left\langle y^{2}+z,x,x^{2}-y\right\rangle$.

In Exercises 15–22., a closed surface $\mathcal{S}$ and a vector field $\vec{F}$ are given. Find the outward flux of $\vec{F}$ over $\mathcal{S}$ either through direct computation or through the Divergence Theorem.

1. 15.

$\mathcal{S}$ is the sphere of radius 3 centered at the origin; $\vec{F}=x\vec{\imath}+2y\vec{\jmath}+3z\vec{k}$

2. 16.

$\mathcal{S}$ is the boundary of the solid cube $S=\{\,(x,y,z):0\leq x,y,z\leq 1\,\}$; $\vec{F}=x\vec{\imath}+y\vec{\jmath}+z\vec{k}$

3. 17.

$\mathcal{S}$ is the unit sphere centered at the origin; $\vec{F}=x^{3}\vec{\imath}+y^{3}\vec{\jmath}+z^{3}\vec{k}$

4. 18.

$\mathcal{S}$ is the unit sphere centered at the origin; $\vec{F}=2\vec{\imath}+3\vec{\jmath}+5\vec{k}$

5. 19.

$\mathcal{S}$ is the surface formed by the intersections of $z=0$ and $z=(x^{2}-1)(y^{2}-1)$; $\vec{F}=\left\langle x^{2}+1,yz,xz^{2}\right\rangle$.

6. 20.

$\mathcal{S}$ is the surface formed by the intersections of the planes $z=\frac{1}{2}(3-x)$, $x=1$, $y=0$, $y=2$ and $z=0$; $\vec{F}=\left\langle x,y^{2},z\right\rangle$.

7. 21.

$\mathcal{S}$ is the surface formed by the intersections of the planes $z=2y$, $y=4-x^{2}$ and $z=0$; $\vec{F}=\left\langle xz,0,xz\right\rangle$.

8. 22.

$\mathcal{S}$ is the surface formed by the intersections of the cylinder $z=1-x^{2}$ and the planes $y=-2$, $y=2$ and $z=0$; $\vec{F}=\left\langle 0,y^{3},0\right\rangle$.

In Exercises 23–26., a closed curve $C$ that is the boundary of a surface $\mathcal{S}$ is given along with a vector field $\vec{F}$. Find the circulation of $\vec{F}$ around $C$ either through direct computation or through Stokes’ Theorem.

1. 23.

$C$ is the curve whose $x$- and $y$-values are determined by the three sides of a triangle with vertices at $(-1,0)$, $(1,0)$ and $(0,1)$, traversed in that order, and the $z$-values are determined by the function $z=xy$; $\vec{F}=\left\langle z-y^{2},x,z\right\rangle$.

2. 24.

$C$ is the curve whose $x$- and $y$-values are given by $\vec{r}(t)=\left\langle 2\cos t,2\sin t\right\rangle$ and the $z$-values are determined by the function $z=x^{2}+y^{3}-3y+1$; $\vec{F}=\left\langle-y,x,z\right\rangle$.

3. 25.

$C$ is the curve whose $x$- and $y$-values are given by $\vec{r}(t)=\left\langle\cos t,3\sin t\right\rangle$ and the $z$-values are determined by the function $z=5-2x-y$; $\vec{F}=\left\langle-\frac{1}{3}y,3x,\frac{2}{3}y-3x\right\rangle$.

4. 26.

$C$ is the curve whose $x$- and $y$-values are sides of the square with vertices at $(1,1)$, $(-1,1)$, $(-1,-1)$ and $(1,-1)$, traversed in that order, and the $z$-values are determined by the function $z=10-5x-2y$; $\vec{F}=\left\langle 5y^{2},2y^{2},y^{2}\right\rangle$.

Exercises 27–30. are designed to challenge your understanding and require no computation.

1. 27.

Let $\mathcal{S}$ be any closed surface enclosing a domain $D$. Consider $\vec{F}_{1}=\left\langle x,0,0\right\rangle$ and $\vec{F}_{2}=\left\langle y,y^{2},z-2yz\right\rangle$.

These fields are clearly very different. Why is it that the total outward flux of each field across $\mathcal{S}$ is the same?

2. 28.
(a) Green’s Theorem can be used to find the area of a region enclosed by a curve by evaluating a line integral with the appropriate choice of vector field $\vec{F}$. What condition on $\vec{F}$ makes this possible? (b) Likewise, Stokes’ Theorem can be used to find the surface area of a region enclosed by a curve in space by evaluating a line integral with the appropriate choice of vector field $\vec{F}$. What condition on $\vec{F}$ makes this possible?
3. 29.

The Divergence Theorem establishes equality between a particular double integral and a particular triple integral. What types of circumstances would lead one to choose to evaluate the triple integral over the double integral?

4. 30.

Stokes’ Theorem establishes equality between a particular line integral and a particular double integral. What types of circumstances would lead one to choose to evaluate the double integral over the line integral?

1. 31.

Construct a Möbius strip from a piece of paper, then draw a line down its center (like the dotted line in Figure 15.5.1(b)). Cut the Möbius strip along that center line completely around the strip. How many surfaces does this result in? How would you describe them? Are they orientable?

2. 32.
Use a computer algebra system to plot the Möbius strip parameterized as $\vec{r}(u,v)=$ $\left\langle\cos u\,(1+v\cos\tfrac{u}{2}),\sin u\,(1+v\cos\tfrac{u}{2}),v\sin% \tfrac{u}{2}\right\rangle,$ where $0\leq u\leq 2\pi$, $-\tfrac{1}{2}\leq v\leq\tfrac{1}{2}$
3. 33.

Let $\mathcal{S}$ be a closed surface and $\vec{F}$ a smooth vector field. Show that $\displaystyle\iint_{\mathcal{S}}(\operatorname{curl}{\vec{F}}\,)\cdot\vec{n}% \operatorname{d}\!S=0$. (Hint: Split $\mathcal{S}$ in half.)

4. 34.

Show that Green’s Theorem is a special case of Stokes’ Theorem. (Hint: Think of how a vector field $f(x,y)=P(x,y)i+Q(x,y)j$ in $\mathbb{R}^{2}$ can be extended in a natural way to be a vector field in $\mathbb{R}^{3}$.)