$$\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 \|}$$ $$\newcommand{\inner}{\langle #1, #2 \rangle}$$ $$\newcommand{\Span}{\mathrm{span}}$$

# The Divergence Theorem

$$\newcommand{\vecs}{\overset { \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vecd}{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}}$$

We have examined several versions of the Fundamental Theorem of Calculus in higher dimensions that relate the integral around an oriented boundary of a domain to a “derivative” of that entity on the oriented domain. In this section, we state the divergence theorem, which is the final theorem of this type that we will study. The divergence theorem has many uses in physics; in particular, the divergence theorem is used in the field of partial differential equations to derive equations modeling heat flow and conservation of mass. We use the theorem to calculate flux integrals and apply it to electrostatic fields.

Overview of Theorems

Before examining the divergence theorem, it is helpful to begin with an overview of the versions of the Fundamental Theorem of Calculus we have discussed:

1. The Fundamental Theorem of Calculus: $\int_a^b f' (x) \, dx = f(b) - f(a).$ This theorem relates the integral of derivative $$f'$$ over line segment $$[a,b]$$ along the x-axis to a difference of $$f$$ evaluated on the boundary.
2. The Fundamental Theorem for Line Integrals: $\int_C \nabla f \cdot dr = f(P_1) - f(P_0),$ where $$P_0$$ is the initial point of C and $$P_1$$ is the terminal point of C. The Fundamental Theorem for Line Integrals allows path C to be a path in a plane or in space, not just a line segment on the x-axis. If we think of the gradient as a derivative, then this theorem relates an integral of derivative $$\nabla f$$ over path C to a difference of $$f$$ evaluated on the boundary of C.
3. Green’s theorem, circulation form: $\iint_D (Q_x - P_y)dA = \int_C F \cdot dr.$ Since $$Q_x - P_y = curl \, F \cdot k$$ and curl is a derivative of sorts, Green’s theorem relates the integral of derivative curl F over planar region D to an integral of F over the boundary of D.
4. Green’s theorem, flux form: $\iint_D (P_x + Q_y)dA = \int_C F \cdot N \, dS.$ Since $$P_x + Q_y = div \, F$$ and divergence is a derivative of sorts, the flux form of Green’s theorem relates the integral of derivative div F over planar region D to an integral of F over the boundary of D.
5. Stokes’ theorem: $\iint_S curl \, F \cdot dS = \int_C F \cdot dr.$ If we think of the curl as a derivative of sorts, then Stokes’ theorem relates the integral of derivative curlF over surface S (not necessarily planar) to an integral of F over the boundary of S.

## Stating the Divergence Theorem

The divergence theorem follows the general pattern of these other theorems. If we think of divergence as a derivative of sorts, then the divergence theorem relates a triple integral of derivative div F over a solid to a flux integral of F over the boundary of the solid. More specifically, the divergence theorem relates a flux integral of vector field F over a closed surface S to a triple integral of the divergence of F over the solid enclosed by S.

The Divergence Theorem

Let S be a piecewise, smooth closed surface that encloses solid E in space. Assume that S is oriented outward, and let F be a vector field with continuous partial derivatives on an open region containing E (Figure $$\PageIndex{1}$$). Then

$\iiint_E div \, F \, dV = \iint_S F \cdot dS. \label{divtheorem}$ Figure $$\PageIndex{1}$$: The divergence theorem relates a flux integral across a closed surface S to a triple integral over solid E enclosed by the surface.

Recall that the flux form of Green’s theorem states that

$\iint_D div \, dA = \int_C F \cdot N \, dS.$

Therefore, the divergence theorem is a version of Green’s theorem in one higher dimension.

The proof of the divergence theorem is beyond the scope of this text. However, we look at an informal proof that gives a general feel for why the theorem is true, but does not prove the theorem with full rigor. This explanation follows the informal explanation given for why Stokes’ theorem is true.

Proof

Let B be a small box with sides parallel to the coordinate planes inside E (Figure $$\PageIndex{2a}$$). Let the center of B have coordinates $$(x,y,z)$$ and suppose the edge lengths are $$\Delta x, \, \Delta y$$, and $$\Delta z$$. (Figure $$\PageIndex{1b}$$). The normal vector out of the top of the box is $$\mathbf{\hat k}$$ and the normal vector out of the bottom of the box is $$-k$$. The dot product of $$\vecs F = \langle P, Q, R \rangle$$ with $$\mathbf{\hat k}$$ is R and the dot product with $$-k$$ is $$-R$$. The area of the top of the box (and the bottom of the box) Figure $$\PageIndex{2}$$: (a) A small box B inside surface E has sides parallel to the coordinate planes. (b) Box B has side lengths$$\Delta x, \, \Delta y$$, and $$\Delta z$$ (c) If we look at the side view of B, we see that, since $$(x,y,z)$$ is the center of the box, to get to the top of the box we must travel a vertical distance of $$\Delta /z$$ up from $$(x,y,z)$$. Similarly, to get to the bottom of the box we must travel a distance $$\Delta /z$$ down from $$(x,y,z)$$.

The flux out of the top of the box can be approximated by $$R \left(x,\, y,\, z + \frac{\Delta z}{2}\right) \Delta x \Delta y$$ (Figure $$\PageIndex{2c}$$) and the flux out of the bottom of the box is $$- R \left(x,\, y,\, z - \frac{\Delta z}{2}\right) \Delta x \Delta y$$. If we denote the difference between these values as $$\Delta R$$, then the net flux in the vertical direction can be approximated by $$\Delta R \Delta x \Delta y$$. However,

$\Delta R \Delta x \Delta y = \left(\frac{\Delta R}{\Delta z}\right) \Delta x \Delta y \Delta z \approx \left(\frac{\partial R}{\partial z}\right) \Delta V.\nonumber$

Therefore, the net flux in the vertical direction can be approximated by $$\left(\frac{\partial R}{\partial z}\right)\Delta V$$. Similarly, the net flux in the x-direction can be approximated by$$\left(\frac{\partial P}{\partial x}\right)\Delta V$$and the net flux in the y-direction can be approximated by$$\left(\frac{\partial Q}{\partial y}\right)\Delta V$$. Adding the fluxes in all three directions gives an approximation of the total flux out of the box:

$Total \, flux \approx \left(\frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} + \frac{\partial R}{\partial z} \right) \Delta V = div \, F \Delta V. \nonumber$

This approximation becomes arbitrarily close to the value of the total flux as the volume of the box shrinks to zero.

The sum of $$div F \Delta V$$ over all the small boxes approximating E is approximately $$\iiint_E div \, F dV$$. On the other hand, the sum of $$div \, F \Delta V$$ over all the small boxes approximating E is the sum of the fluxes over all these boxes. Just as in the informal proof of Stokes’ theorem, adding these fluxes over all the boxes results in the cancelation of a lot of the terms. If an approximating box shares a face with another approximating box, then the flux over one face is the negative of the flux over the shared face of the adjacent box. These two integrals cancel out. When adding up all the fluxes, the only flux integrals that survive are the integrals over the faces approximating the boundary of E. As the volumes of the approximating boxes shrink to zero, this approximation becomes arbitrarily close to the flux over S.

$$\Box$$

Example $$\PageIndex{1}$$: Verifying the Divergence Theorem

Verify the divergence theorem for vector field $$F = \langle x - y, \, x + z, \, z - y \rangle$$ and surface S that consists of cone $$x^2 + y^2 = z^2, \, 0 \leq z \leq 1$$, and the circular top of the cone (see the following figure). Assume this surface is positively oriented. Solution

Let E be the solid cone enclosed by S. To verify the theorem for this example, we show that

$\iiint_E div \, F dV = \iint_S F \cdot dS\nonumber$

by calculating each integral separately.

To compute the triple integral, note that $$div \, F = P_x + Q_y + R_z = 2$$, and therefore the triple integral is

\begin{align*} \iiint_E div \, dV = 2 \iiint_E dV \\[4pt] = 2 \, (volume \, of \, E). \end{align*}

The volume of a right circular cone is given by $$\pi r^2 \frac{h}{3}$$. In this case, $$h = r = 1$$. Therefore,

$\iiint_E div \, F dV = 2 \, (volume \, of \, E) = \frac{2\pi}{3}.\nonumber$

To compute the flux integral, first note that S is piecewise smooth; S can be written as a union of smooth surfaces. Therefore, we break the flux integral into two pieces: one flux integral across the circular top of the cone and one flux integral across the remaining portion of the cone. Call the circular top $$S_1$$ and the portion under the top $$S_2$$. We start by calculating the flux across the circular top of the cone. Notice that $$S_1$$ has parameterization

$r(u,v) = \langle u \, \cos v, \, u \, \sin v, \, 1 \rangle, \, 0 \leq u \leq 1, \, 0 \leq v \leq 2\pi.\nonumber$

Then, the tangent vectors are $$t_u = \langle \cos v, \, \sin v, \, 0 \rangle$$ and $$t_v = \langle -u \, \cos v, \, u \, \sin v, 0 \rangle$$. Therefore, the flux across $$S_1$$ is

\begin{align*} \iint_{S_1} F \cdot dS &= \int_0^1 \int_0^{2\pi} F ( r ( u,v)) \cdot (t_u \times t_v) \, dA \\[4pt] &= \int_0^1 \int_0^{2\pi} \langle u \, \cos v - u \, \sin v, \, u \, \cos v + 1, \, 1 - u \, \sin v \rangle \cdot \langle 0,0,v \rangle \, dv\, du \\[4pt] &= \int_0^1 \int_0^{2\pi} u - u^2 \sin v \, dv du \\[4pt] &= \pi. \end{align*}

We now calculate the flux over $$S_2$$. A parameterization of this surface is

$r(u,v) = \langle u \, \cos v, \, u \, \sin v, \, u \rangle, \, 0 \leq u \leq 1, \, 0 \leq v \leq 2\pi.\nonumber$

The tangent vectors are $$t_u = \langle \cos v, \, \sin v, \, 1 \rangle$$ and $$t_v = \langle -u \, \sin v, \, u \, \cos v, 0 \rangle$$, so the cross product is

$t_u \times t_v = \langle - u \, \cos v, \, -u \, \sin v, \, u \rangle.\nonumber$

Notice that the negative signs on the x and y components induce the negative (or inward) orientation of the cone. Since the surface is positively oriented, we use vector $$t_v \times t_u = \langle u \, \cos v, \, u \, \sin v, \, -u \rangle$$ in the flux integral. The flux across $$S_2$$ is then

\begin{align*} \iint_{S_1} F \cdot dS &= \int_0^1 \int_0^{2\pi} F ( r ( u,v)) \cdot (t_u \times t_v) \, dA \\[4pt] &= \int_0^1 \int_0^{2\pi} \langle u \, \cos v - u \, \sin v, \, u \, \cos v + u, \, u \, - \sin v \rangle \cdot \langle u \, \cos v, \, u \, \sin v, \, -u \rangle \\[4pt] &= \int_0^1 \int_0^{2\pi} u^2 \cos^2 v + 2u^2 \sin v - u^2 dv\,du &= \frac{\pi}{3} \end{align*}

The total flux across S is

$\iint_{S_2} F \cdot dS = \iint_{S_1} F \cdot dS + \iint_{S_2} F \cdot dS = \frac{2\pi}{3} = \iiint_E div \, F dV,\nonumber$

and we have verified the divergence theorem for this example.

Exercise $$\PageIndex{1}$$

Verify the divergence theorem for vector field $$F (x,y,z) = \langle x + y + z, \, y, \, 2x - y \rangle$$ and surface S given by the cylinder $$x^2 + y^2 = 1, \, 0 \leq z \leq 3$$ plus the circular top and bottom of the cylinder. Assume that S is positively oriented.

Hint

Calculate both the flux integral and the triple integral with the divergence theorem and verify they are equal.

Both integrals equal $$6\pi$$.

​​Recall that the divergence of continuous field F at point P is a measure of the “outflowing-ness” of the field at P. If F represents the velocity field of a fluid, then the divergence can be thought of as the rate per unit volume of the fluid flowing out less the rate per unit volume flowing in. The divergence theorem confirms this interpretation. To see this, let P be a point and let $$B_{\tau}$$ be a ball of small radius r centered at P (Figure $$\PageIndex{3}$$). Let $$S_{\tau}$$ be the boundary sphere of $$B_{\tau}$$. Since the radius is small and F is continuous, $$div \, F(Q) \approx div \, F(P)$$ for all other points Q in the ball. Therefore, the flux across $$S_{\tau}$$ can be approximated using the divergence theorem:

$\iint_{S_{\tau}} F \cdot dS = \iiint_{B_{\tau}} div \, F dV \approx \iiint_{B_{\tau}} div F(P) dV.\nonumber$

Since $$div F(P)$$ is a constant,

$\iiint_{B_{\tau}} div F(P) dV = div \, F(P) \, V(B_{\tau}).\nonumber$

Therefore, flux $\iint_{S_{\tau}} F \cdot dS$ can be approximated by $$F(P) \, V(B_{\tau})$$. This approximation gets better as the radius shrinks to zero, and therefore

$div \, F(P) = \lim_{\tau \rightarrow 0} \frac{1}{V(B_{\tau})} \iint_{S_{\tau}} F \cdot dS.\nonumber$

This equation says that the divergence at P is the net rate of outward flux of the fluid per unit volume. Figure $$\PageIndex{3}$$: Ball $$B_{\tau}$$ of small radius r centered at P.

## Using the Divergence Theorem

The divergence theorem translates between the flux integral of closed surface S and a triple integral over the solid enclosed by S. Therefore, the theorem allows us to compute flux integrals or triple integrals that would ordinarily be difficult to compute by translating the flux integral into a triple integral and vice versa.

Example $$\PageIndex{2}$$: Applying the Divergence Theorem

Calculate the surface integral

$\iint_S F \cdot dS, \nonumber$

where $$S$$ is cylinder $$x^2 + y^2 = 1, \, 0 \leq z \leq 2$$, including the circular top and bottom, and $$F = \left\langle \frac{x^3}{3} + yz, \, \frac{y^3}{3} - \sin (xz), \, z - x - y \right\rangle$$.

Solution

We could calculate this integral without the divergence theorem, but the calculation is not straightforward because we would have to break the flux integral into three separate integrals: one for the top of the cylinder, one for the bottom, and one for the side. Furthermore, each integral would require parameterizing the corresponding surface, calculating tangent vectors and their cross product..

By contrast, the divergence theorem allows us to calculate the single triple integral

$\iiint_E div \, F \, dV,\nonumber$

where E is the solid enclosed by the cylinder. Using the divergence theorem (Equation \ref{divtheorem}) and converting to cylindrical coordinates, we have

\begin{align*} \iint_S F \cdot dS &= \iiint_E div \, F \, dV, \\[4pt] &= \iiint_E (x^2 + y^2 + 1) \, dV \\[4pt] &= \int_0^{2\pi} \int_0^1 \int_0^2 (r^2 + 1) \, r \, dz \, dr \, d\theta \\[4pt] &= \frac{3}{2} \int_0^{2\pi} d\theta \\[4pt] &= 3\pi. \end{align*}

Exercise $$\PageIndex{2}$$

Use the divergence theorem to calculate flux integral $\iint_S F \cdot dS,\nonumber$ where S is the boundary of the box given by $$0 \leq x \leq 2, \, 1 \leq y \leq 4, \, 0 \leq z \leq 1$$ and $$F = \langle x^2 + yz, \, y - z, \, 2x + 2y + 2z \rangle$$ (see the following figure). Hint

Calculate the corresponding triple integral.

30

Example $$\PageIndex{3}$$: Applying the Divergence Theorem

Let $$v = \left\langle - \frac{y}{z}, \, \frac{x}{z}, \, 0 \right\rangle$$ be the velocity field of a fluid. Let C be the solid cube given by $$1 \leq x \leq 4, \, 2 \leq y \leq 5, \, 1 \leq z \leq 4$$, and let S be the boundary of this cube (see the following figure). Find the flow rate of the fluid across S. Figure $$\PageIndex{4}$$: Vector field $$v = \left\langle - \frac{y}{z}, \, \frac{x}{z}, \, 0 \right\rangle$$.

Solution

The flow rate of the fluid across S is $$\iint_S v \cdot dS$$. Before calculating this flux integral, let’s discuss what the value of the integral should be. Based on Figure $$\PageIndex{4}$$, we see that if we place this cube in the fluid (as long as the cube doesn’t encompass the origin), then the rate of fluid entering the cube is the same as the rate of fluid exiting the cube. The field is rotational in nature and, for a given circle parallel to the xy-plane that has a center on the z-axis, the vectors along that circle are all the same magnitude. That is how we can see that the flow rate is the same entering and exiting the cube. The flow into the cube cancels with the flow out of the cube, and therefore the flow rate of the fluid across the cube should be zero.

To verify this intuition, we need to calculate the flux integral. Calculating the flux integral directly requires breaking the flux integral into six separate flux integrals, one for each face of the cube. We also need to find tangent vectors, compute their cross product. However, using the divergence theorem makes this calculation go much more quickly:

\begin{align*} \iint_S v \cdot dS = \iiint_C div (v) \, dV \\[4pt] = \iiint_C 0 \, dV = 0.\end{align*}

Therefore the flux is zero, as expected.

Exercise $$\PageIndex{3}$$

Let $$v = \left\langle \frac{x}{z}, \, \frac{y}{z}, \, 0 \right\rangle$$ be the velocity field of a fluid. Let C be the solid cube given by $$1 \leq x \leq 4, \, 2 \leq y \leq 5, \, 1 \leq z \leq 4$$, and let S be the boundary of this cube (see the following figure). Find the flow rate of the fluid across S. Hint

Use the divergence theorem and calculate a triple integral

$$9 \, \ln (16)$$

Example illustrates a remarkable consequence of the divergence theorem. Let S be a piecewise, smooth closed surface and let F be a vector field defined on an open region containing the surface enclosed by S. If F has the form $$F = \langle f (y,z), \, g(x,z), \, h(x,y)\rangle$$, then the divergence of F is zero. By the divergence theorem, the flux of F across S is also zero. This makes certain flux integrals incredibly easy to calculate. For example, suppose we wanted to calculate the flux integral $$\iint_S F \cdot dS$$ where S is a cube and

$F = \langle \sin (y) \, e^{yz}, \, x^2z^2, \, \cos (xy) \, e^{\sin x} \rangle.$

Calculating the flux integral directly would be difficult, if not impossible, using techniques we studied previously. At the very least, we would have to break the flux integral into six integrals, one for each face of the cube. But, because the divergence of this field is zero, the divergence theorem immediately shows that the flux integral is zero.

We can now use the divergence theorem to justify the physical interpretation of divergence that we discussed earlier. Recall that if F is a continuous three-dimensional vector field and P is a point in the domain of F, then the divergence of F at P is a measure of the “outflowing-ness” of F at P. If F represents the velocity field of a fluid, then the divergence of F at P is a measure of the net flow rate out of point P (the flow of fluid out of P less the flow of fluid in to P). To see how the divergence theorem justifies this interpretation, let $$B_{\tau}$$ be a ball of very small radius r with center P, and assume that $$B_{\tau}$$ is in the domain of F. Furthermore, assume that $$B_{\tau}$$ has a positive, outward orientation. Since the radius of $$B_{\tau}$$ is small and F is continuous, the divergence of F is approximately constant on $$B_{\tau}$$. That is, ifv $$P'$$ is any point in $$B_{\tau}$$, then $$div \, F(P) \approx div \, F(P')$$. Let $$S_{\tau}$$ denote the boundary sphere of $$B_{\tau}$$. We can approximate the flux across $$S_{\tau}$$ using the divergence theorem as follows:

$\iint_{S_{\tau}} F \cdot dS = \iiint_{B_{\tau}} div \, F \, dV$

$\approx \iiint_{B_{\tau}} div \, F (P) \, dV$

$= div \, F (P) \, V(B_{\tau}).$

As we shrink the radius r to zero via a limit, the quantity $$div \, F (P) \, V(B_{\tau})$$ gets arbitrarily close to the flux. Therefore,

$div \, F(P) = \lim_{\tau \rightarrow 0} \frac{1}{V(B_{\tau})} \iint_{S_{\tau}} F \cdot dS$

and we can consider the divergence at P as measuring the net rate of outward flux per unit volume at P. Since “outflowing-ness” is an informal term for the net rate of outward flux per unit volume, we have justified the physical interpretation of divergence we discussed earlier, and we have used the divergence theorem to give this justification.

## Application to Electrostatic Fields

The divergence theorem has many applications in physics and engineering. It allows us to write many physical laws in both an integral form and a differential form (in much the same way that Stokes’ theorem allowed us to translate between an integral and differential form of Faraday’s law). Areas of study such as fluid dynamics, electromagnetism, and quantum mechanics have equations that describe the conservation of mass, momentum, or energy, and the divergence theorem allows us to give these equations in both integral and differential forms.

One of the most common applications of the divergence theorem is to electrostatic fields. An important result in this subject is Gauss’ law. This law states that if S is a closed surface in electrostatic field E, then the flux of E across S is the total charge enclosed by S (divided by an electric constant). We now use the divergence theorem to justify the special case of this law in which the electrostatic field is generated by a stationary point charge at the origin.

If $$(x,y,z)$$ is a point in space, then the distance from the point to the origin is $$r = \sqrt{x^2 + y^2 + z^2}$$. Let $$F_{\tau}$$ denote radial vector field $$F_{\tau} = \dfrac{1}{\tau^2} \left\langle \dfrac{x}{\tau}, \, \dfrac{y}{\tau}, \, \dfrac{z}{\tau}\right\rangle$$.The vector at a given position in space points in the direction of unit radial vector $$\left\langle \dfrac{x}{\tau}, \, \dfrac{y}{\tau}, \, \dfrac{z}{\tau}\right\rangle$$ and is scaled by the quantity $$1/\tau^2$$. Therefore, the magnitude of a vector at a given point is inversely proportional to the square of the vector’s distance from the origin. Suppose we have a stationary charge of q Coulombs at the origin, existing in a vacuum. The charge generates electrostatic field E given by

$E = \dfrac{q}{4\pi \epsilon_0}F_{\tau},$

where the approximation $$\epsilon_0 = 8.854 \times 10^{-12}$$ farad (F)/m is an electric constant. (The constant $$\epsilon_0$$ is a measure of the resistance encountered when forming an electric field in a vacuum.) Notice that E is a radial vector field similar to the gravitational field described in [link]. The difference is that this field points outward whereas the gravitational field points inward. Because

$E = \dfrac{q}{4\pi \epsilon_0}F_{\tau} = \dfrac{q}{4\pi \epsilon_0}\left(\dfrac{1}{\tau^2} \left\langle \dfrac{x}{\tau}, \, \dfrac{y}{\tau}, \, \dfrac{z}{\tau}\right\rangle\right),$

we say that electrostatic fields obey an inverse-square law. That is, the electrostatic force at a given point is inversely proportional to the square of the distance from the source of the charge (which in this case is at the origin). Given this vector field, we show that the flux across closed surface S is zero if the charge is outside of S, and that the flux is $$q/epsilon_0$$ if the charge is inside of S. In other words, the flux across S is the charge inside the surface divided by constant $$\epsilon_0$$. This is a special case of Gauss’ law, and here we use the divergence theorem to justify this special case.

To show that the flux across S is the charge inside the surface divided by constant $$\epsilon_0$$, we need two intermediate steps. First we show that the divergence of $$F_{\tau}$$ is zero and then we show that the flux of $$F_{\tau}$$ across any smooth surface S is either zero or $$4\pi$$. We can then justify this special case of Gauss’ law.

Example $$\PageIndex{4}$$: The Divergence of $$F_{\tau}$$ is Zero

Verify that the divergence of $$F_{\tau}$$ is zero where $$F_{\tau}$$ is defined (away from the origin).

Solution

Since $$\tau = \sqrt{x^2 + y^2 + z^2}$$, the quotient rule gives us

\begin{align*} \dfrac{\partial}{\partial x} \left( \dfrac{x}{\tau^3} \right) = \dfrac{\partial}{\partial x} \left( \dfrac{x}{(x^2+y^2+z^2)^{3/2}} \right) \\[4pt] = \dfrac{(x^2+y^2+z^2)^{3/2} - x\left[\dfrac{3}{2} (x^2+y^2+z^2)^{1/2}2x\right]}{(x^2+y^2+z^2)^3} \\[4pt] = \dfrac{\tau^3 -3x^2\tau}{\tau^6} = \dfrac{\tau^2 - 3x^2}{\tau^5}. \end{align*}

Similarly,

$\dfrac{\partial}{\partial y} \left( \dfrac{y}{\tau^3} \right) = \dfrac{\tau^2 - 3y^2}{\tau^5} \, and \, \dfrac{\partial}{\partial z} \left( \dfrac{z}{\tau^3} \right) = \dfrac{\tau^2 - 3z^2}{\tau^5}. \nonumber$

Therefore,

\begin{align*} div \, F_{\tau} = \dfrac{\tau^2 - 3x^2}{\tau^5} + \dfrac{\tau^2 - 3y^2}{\tau^5} + \dfrac{\tau^2 - 3z^2}{\tau^5} \\[4pt] = \dfrac{3\tau^2 - 3(x^2+y^2+z^2)}{\tau^5} \\[4pt] = \dfrac{3\tau^2 - 3\tau^2}{\tau^5} = 0. \end{align*}

Notice that since the divergence of $$F_{\tau}$$ is zero and E is $$F_{\tau}$$ scaled by a constant, the divergence of electrostatic field E is also zero (except at the origin).

Flux across a Smooth Surface

Let S be a connected, piecewise smooth closed surface and let $$F_{\tau} = \dfrac{1}{\tau^2} \left\langle \dfrac{x}{\tau}, \, \dfrac{y}{\tau}, \, \dfrac{z}{\tau}\right \rangle$$. Then,

$\iint_SF_{\tau} \cdot dS = \begin{cases}0 \text{if S does not encompass the origin} \\ 4\pi \text{if S encompasses the origin.} \end{cases}$

In other words, this theorem says that the flux of $$F_{\tau}$$ across any piecewise smooth closed surface S depends only on whether the origin is inside of S.

Proof

The logic of this proof follows the logic of [link], only we use the divergence theorem rather than Green’s theorem.

First, suppose that S does not encompass the origin. In this case, the solid enclosed by S is in the domain of $$F_{\tau}$$, and since the divergence of $$F_{\tau}$$ is zero, we can immediately apply the divergence theorem and find that $\iint_S F \cdot dS$ is zero.

Now suppose that S does encompass the origin. We cannot just use the divergence theorem to calculate the flux, because the field is not defined at the origin. Let $$S_a$$ be a sphere of radius a inside of S centered at the origin. The outward normal vector field on the sphere, in spherical coordinates, is

$t_{\phi} \times t_{\theta} = \langle a^2 \cos \theta \, \sin^2 \phi, \, a^2 \sin \theta \, \sin^2 \phi, \, a^2 \sin \phi \, \cos \phi \rangle$

(see [link]). Therefore, on the surface of the sphere, the dot product $$F_{\tau} \cdot N$$ (in spherical coordinates) is

\begin{align*} F_{\tau} \cdot N = \left \langle \dfrac{\sin \phi \, \cos \theta}{a^2}, \, \dfrac{\sin \phi \, \sin \theta}{a^2}, \, \dfrac{\cos \phi}{a^2} \right \rangle \cdot \langle a^2 \cos \theta \, \sin^2 \phi, a^2 \sin \theta \, \sin^2 \phi, \, a^2 \sin \phi \, \cos \phi \rangle \\[4pt] = \sin \phi ( \langle \sin \phi \, \cos \theta, \, \sin \phi \, \sin \theta, \, \cos \phi \rangle \cdot \langle \sin \phi \, \cos \theta, \sin \phi \, \sin \theta, \, \cos \phi \rangle ) \\[4pt] = \sin \phi. \end{align*}

The flux of $$F_{\tau}$$ across $$S_a$$ is

$\iint_{S_a} F_{\tau} \cdot N dS = \int_0^{2\pi} \int_0^{\pi} \sin \phi \, d\phi \, d\theta = 4\pi.$

Now, remember that we are interested in the flux across S, not necessarily the flux across $$S_a$$. To calculate the flux across S, let E be the solid between surfaces $$S_a$$ and S. Then, the boundary of E consists of $$S_a$$ and S. Denote this boundary by $$S - S_a$$ to indicate that S is oriented outward but now $$S_a$$ is oriented inward. We would like to apply the divergence theorem to solid E. Notice that the divergence theorem, as stated, can’t handle a solid such as E because E has a hole. However, the divergence theorem can be extended to handle solids with holes, just as Green’s theorem can be extended to handle regions with holes. This allows us to use the divergence theorem in the following way. By the divergence theorem,

\begin{align*} \iint_{S-S_a} F_{\tau} \cdot dS = \iint_S F_{\tau} \cdot dS - \iint_{S_a} F_{\tau} \cdot dS \\[4pt] = \iiint_E div \, F_{\tau} \, dV \\[4pt] = \iiint_E 0 \, dV = 0. \end{align*}

Therefore,

$\iint_S F_{\tau} \cdot dS = \iint_{S_a} F_{\tau} \cdot dS = 4\pi, \nonumber$

and we have our desired result.

$$\Box$$

Now we return to calculating the flux across a smooth surface in the context of electrostatic field $$E = \dfrac{q}{4\pi \epsilon_0} F_{\tau}$$ of a point charge at the origin. Let S be a piecewise smooth closed surface that encompasses the origin. Then

\begin{align*} \iint_S E \cdot dS = \iint_S \dfrac{q}{4\pi \epsilon_0} F_{\tau} \cdot dS\\[4pt] = \dfrac{q}{4\pi \epsilon_0} \iint_S F_{\tau} \cdot dS \\[4pt] = \dfrac{q}{\epsilon_0}. \end{align*}

If S does not encompass the origin, then

$\iint_S E \cdot dS = \dfrac{q}{4\pi \epsilon_0} \iint_S F_{\tau} \cdot dS = 0. \nonumber$

Therefore, we have justified the claim that we set out to justify: the flux across closed surface S is zero if the charge is outside of S, and the flux is $$q/\epsilon_0$$ if the charge is inside of S.

This analysis works only if there is a single point charge at the origin. In this case, Gauss’ law says that the flux of E across S is the total charge enclosed by S. Gauss’ law can be extended to handle multiple charged solids in space, not just a single point charge at the origin. The logic is similar to the previous analysis, but beyond the scope of this text. In full generality, Gauss’ law states that if S is a piecewise smooth closed surface and Q is the total amount of charge inside of S, then the flux of E across S is $$Q/\epsilon_0$$.

Example $$\PageIndex{5}$$: Using Gauss’ law

Suppose we have four stationary point charges in space, all with a charge of 0.002 Coulombs (C). The charges are located at $$(0,0,1), \, (1,1,4), (-1,0,0)$$, and $$(-2,-2,2)$$. Let E denote the electrostatic field generated by these point charges. If S is the sphere of radius 2 oriented outward and centered at the origin, then find

$\iint_S E \cdot dS. \nonumber$

Solution

According to Gauss’ law, the flux of E across S is the total charge inside of S divided by the electric constant. Since S has radius 2, notice that only two of the charges are inside of S: the charge at $$0,1,1)$$ and the charge at $$(-1,0,0)$$. Therefore, the total charge encompassed by S is 0.004 and, by Gauss’ law,

$\iint_S E \cdot dS = \dfrac{0.004}{8.854 \times 10^{-12}} \approx 4.418 \times 10^9 \, V - m. \nonumber$

Exercise $$\PageIndex{4}$$

Work the previous example for surface S that is a sphere of radius 4 centered at the origin, oriented outward.

Hint

Use Gauss’ law.

$$\approx 6.777 \times 10^9$$

## Key Concepts

• The divergence theorem relates a surface integral across closed surface S to a triple integral over the solid enclosed by S. The divergence theorem is a higher dimensional version of the flux form of Green’s theorem, and is therefore a higher dimensional version of the Fundamental Theorem of Calculus.
• The divergence theorem can be used to transform a difficult flux integral into an easier triple integral and vice versa.
• The divergence theorem can be used to derive Gauss’ law, a fundamental law in electrostatics.

## Key Equations

• Divergence theorem $\iiint_E div \, F \, dV = \iint_S F \cdot dS \nonumber$

## Glossary

divergence theorem
a theorem used to transform a difficult flux integral into an easier triple integral and vice versa
Gauss’ law
if S is a piecewise, smooth closed surface in a vacuum and Q is the total stationary charge inside of S, then the flux of electrostatic field E across S is $$Q/\epsilon_0$$
inverse-square law
the electrostatic force at a given point is inversely proportional to the square of the distance from the source of the charge