16.8: The Divergence Theorem and a Unified Theory
- Page ID
- 2642
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)When we looked at Green's Theorem, we saw that there was a relationship between a region and the curve that encloses it. This gave us the relationship between the line integral and the double integral. Now consider the following theorem:
Let \(Q\) be a solid region bounded by a closed surface oriented with outward pointing unit normal vector \(\text{n}\), and let \(\textbf{F}\) be a differentiable vector field (i.e., components have continuous partial derivatives). Then
\[ \iint \limits_{S} \textbf{F} \cdot \text{n} \, dS = \iiint \limits_{Q} \nabla \cdot \textbf{F}\,dv.\nonumber \]
Find
\[\iint \limits_{S} \textbf{F} \cdot \textbf{N} \,ds \nonumber \nonumber \]
where
\[ \textbf{F} (x,y,z) = y^2 \hat{\textbf{i}} + e^x(1-\cos{(x^2 + z^2)}\hat{\textbf{j}} + (x + z) \hat{\textbf{k}} \nonumber \nonumber \]
and \(S\) is the unit sphere centered at the point \((1,4,6)\) with outwardly pointing normal vector.
Solution
This seemingly difficult problem turns out to be quite easy once we have the divergence theorem. We have
\[ \nabla \cdot \textbf{F} = 0 + 0 + 1 = 1 \nonumber \nonumber \]
Now recall that a triple integral of the function \(1\) is the volume of the solid. Since the solid is a sphere of radius \(1\) we get \(\dfrac{4}{3}\pi\).
As usual, we will make some simplifying remarks and then prove part of the divergence theorem. We assume that the solid is bounded below by \( z = g_1(x,y) \) and above by \( z = g_2(x,y) \).
Notice that the outward pointing normal vector is upward on the top surface and downward for the bottom region. We also note that the divergence theorem can be written as
\[\begin{align*} \iint\limits_{S} \textbf{F} \cdot \textbf{N} dS &= \iint\limits_{S} (M \hat{\textbf{i}} \cdot \textbf{N} + N\hat{\textbf{j}} \cdot \textbf{N} + P \hat{\textbf{k}} \cdot \textbf{N} ) dS \\ &= \iiint\limits_{Q} \nabla \cdot \textbf{F} dv \\ &= \iiint\limits_{Q} (M_x + N_y + P_z) dv. \end{align*}\nonumber \]
We will show that
\[\iint\limits_{S} P \hat{\textbf{k}} \cdot \textbf{N} dS = \iiint \limits_{Q} P_x dv. \nonumber \]
We have on the top surface
\[ P\hat{\text{k}} \cdot \text{n}\, dS = P\hat{\text{k}} \cdot \left( (-g_2)_x \hat{\text{i}} - (g_2)_y \hat{\text{j}} + \hat{\text{k}} \right) = P(x,y,g_2(x,y)).\nonumber \]
On the bottom surface, we get
\[P\hat{\text{k}} \cdot \text{n} dS = P\hat{\text{k}} \cdot \left( (g_1)_x \hat{\text{i}} + (g_1)_y \hat{\text{j}} - \hat{\text{k}} \right) = -P(x,y,g_1(x,y)).\nonumber \]
Putting these together we get
\[\iint\limits_S P\hat{\text{k}} \cdot \textbf{N} dS = \iint\limits_R [P(x,y,g_2 (x,y))-P(x,y,g_1 (x,y)) ] dydx . \nonumber \]
For the triple integral, the Fundamental Theorem of Calculus tell us that
\[\begin{align*} \iiint\limits_Q P_z dzdydx &= \iint\limits_R [P(x,y,z) ]_{g_1(x,y)}^{g_2(x,y)} dydx \\ &= \iint\limits_R [P(x,y,g_2(x,y))-P(x,y,g_1(x,y))]dy\,dx. \end{align*}\nonumber \]
\(\square\)
An Interpretation of Divergence
We have seen that the flux is the amount fluid flow per unit time through a surface. If the surface is closed, then the total flux will equal the flow out of the solid minus the flow in. Often in the solid there is a source (such as a star when the flow is electromagnetic radiation) or a sink (such as the earth collecting solar radiation). If we have a small solid \(S(P)\) containing a point \(P\), then the divergence of the vector field is approximately constant, which leads to the approximation
\[ \iiint\limits_Q \nabla \cdot \textbf{F} \, dv \approx \nabla \cdot \textbf{F}(P) \text{Volume.} \nonumber \]
The divergence theorem expresses the approximation
Flux through \(S(P) \approx \nabla \cdot \textbf{F}(P) \)(Volume).
Dividing by the volume, we get that the divergence of \(\textbf{F}\) at \(P\) is the Flux per unit volume.
- If the divergence is positive, then the \(P\) is a source.
- If the divergence is negative, then \(P\) is a sink.
Contributors and Attributions
- Larry Green (Lake Tahoe Community College)
Integrated by Justin Marshall.