Skip to main content
Mathematics LibreTexts

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:

    Divergence 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 \]

    Moving to three dimensions, the divergence theorem provides us with a relationship between a triple integral over a solid and the surface integral over the surface that encloses the solid.
    Example \(\PageIndex{1}\)

    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\).

    Partial Proof

    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) \).

    alt

    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.

    alt

    Contributors and Attributions


    16.8: The Divergence Theorem and a Unified Theory is shared under a CC BY-NC-SA license and was authored, remixed, and/or curated by LibreTexts.

    • Was this article helpful?