Foundations to FrontierPart VI: Calculus IIIBezhig

Module 46: The Divergence theorem and a unified view of calculus

8 exercises, each with a complete solution. 1 figure.

On this page
  1. A. Prerequisites and observable goals
  2. B. The central idea: interior outflow accumulates on the boundary
  3. C. The theorem and its domain requirements
  4. D. Why a singular source changes the argument
  5. E. A method-selection procedure
  6. F. Worked examples
  7. Computing bridge: Conservation as a system-level check
  8. G. Independent practice
  9. H. Complete solutions: Module 46
  10. I. Retention and readiness

A. Prerequisites and observable goals#

Before

Triple integrals, divergence, outward flux, and boundary orientation.

After

Convert closed-surface flux to a volume integral, close an open surface with explicit corrections, handle cavities and singularities, and select the appropriate boundary theorem.

B. The central idea: interior outflow accumulates on the boundary#

The Divergence theorem connects the net outward flux of a vector field through a closed surface to the integral of its divergence throughout the enclosed solid. Internal transfers cancel: flow leaving one small cell enters its neighbor. Only the exterior boundary contributes to net escape. This is the three-dimensional outflow counterpart of Green’s planar flux theorem.

C. The theorem and its domain requirements#

Let EE be a suitable bounded solid with piecewise smooth closed boundary ∂E\partial E, oriented outward. If F\vv F has continuous first partial derivatives on an open neighborhood of the entire solid, then

∬∂EF⋅n dS=∭Ediv⁡F dV.\iint_{\partial E}\vv F\cdot\vv n\dd S=\iiint_E\diver\vv F\dd V.

The field must be regular inside the solid, not merely on its boundary. For a solid with a cavity, outward on the inner boundary points into the cavity, away from the material of the solid. Open surfaces are not closed boundaries; add the missing patches and subtract their fluxes afterward. [O3Div]

The fundamental boundary patterns are:

∫abF′ dx=F(b)−F(a),∫C∇ϕ⋅ dr=ϕ(B)−ϕ(A),\int_a^b F'\dd x=F(b)-F(a),\qquad \int_C\grad\phi\cdot\dd\vv r=\phi(B)-\phi(A),
∬Scurl⁡F⋅n dS=∮∂SF⋅ dr,∭Ediv⁡F dV=∬∂EF⋅n dS.\iint_S\curl\vv F\cdot\vv n\dd S=\oint_{\partial S}\vv F\cdot\dd\vv r, \qquad \iiint_E\diver\vv F\dd V=\iint_{\partial E}\vv F\cdot\vv n\dd S.

These share an interior-derivative/boundary-accumulation structure, but their objects, dimensions, and hypotheses differ. The analogy helps choose a theorem; it does not license exchanging curl, divergence, or gradient.

D. Why a singular source changes the argument#

For F=r/∥r∥3\vv F=\vv r/\|\vv r\|^3 away from the origin, direct differentiation gives divergence zero there. On a sphere of radius RR, the outward normal is r/R\vv r/R, so F⋅n=1/R2\vv F\cdot\vv n=1/R^2, giving total flux 4π4\pi. There is no contradiction: the enclosed origin is outside the field’s domain. The theorem cannot integrate zero divergence throughout a solid where the field is not defined everywhere.

On a spherical shell that excludes the origin, the field is smooth. The outer flux is 4π4\pi and the inner outward-from-shell flux is −4π-4\pi, so the total is zero, agreeing with the volume integral. The missing source point and the cavity’s orientation together explain the difference between the two situations.

E. A method-selection procedure#

If asked for flux, determine whether the surface is closed. For a closed surface, compute divergence and compare the volume integral’s complexity with direct surface integration. For an open patch, determine whether adding a simple cap makes the theorem useful; track its flux with the outward orientation of the enclosed solid. Search for singularities anywhere in that solid. For circulation, consider a potential, Green, or Stokes instead; divergence is not a shortcut for every vector integral.

F. Worked examples#

An ellipsoid and a simple divergence. For F=(x,y,z)\vv F=(x,y,z), divergence is 33. Over an ellipsoid with semiaxes a,b,ca,b,c, total outward flux is 3(4πabc/3)=4πabc3(4\pi abc/3)=4\pi abc, provided the entire ellipsoid boundary is included. This avoids parameterizing its surface normals.

Closing an open cap. Suppose a cap and a horizontal base disk bound a solid. Write Φcap+Φbase=∭Ediv⁡F dV\Phi_{\text{cap}}+\Phi_{\text{base}}=\iiint_E\diver\vv F\dd V. Only then solve for the requested cap flux. A downward base normal can make Φbase\Phi_{\text{base}} negative, so subtracting it increases the cap result. This is exactly the same sign discipline as subtracting an entire f(x)f(x) in a difference quotient, now applied to oriented contributions.

Planar analogue: the exact sampled field (x,y)(x,y) points outward through a circle. Its divergence is positive. In a partition, internal edge flux cancels, leaving the boundary flux; the three-dimensional theorem uses the same cancellation across internal faces.

Computing bridge: Conservation as a system-level check#

The cancellation of internal fluxes is a mathematical prototype for a conservation test. A numerical transport model may track change in stored mass plus boundary outflow minus sources; an accounting service may track conserved totals across a transfer. These are different domains, not interchangeable theorems, but both benefit from explicitly naming what can enter, leave, or be created. J03 checks transaction effects; T05 checks the conservation laws appropriate to its mathematical model.

G. Independent practice#

F foundation C standard course work T transfer S optional stretch L runnable lab or project. Attempt the complete set before opening its solutions.

M46.01FSolution

Find outward flux of F=(x,y,z)\vv F=(x,y,z) through the full boundary of [0,1]×[0,2]×[0,3][0,1]\times[0,2]\times[0,3].

M46.02CSolution

Find outward flux of F=(x3,y3,z3)\vv F=(x^3,y^3,z^3) through the unit sphere.

M46.03CSolution

Find outward flux of F=(x,y,z)\vv F=(x,y,z) through x2+y2/4+z2/9=1x^2+y^2/4+z^2/9=1.

M46.04TSolution

Find upward flux of F=(x,y,z+1)\vv F=(x,y,z+1) through the paraboloid cap z=4−x2−y2z=4-x^2-y^2, z≥0z\ge0, by adding its base disk.

M46.05TSolution

Show that the total outward flux of any constant vector field through any closed surface bounding a suitable solid is zero. Explain why individual patches may still have nonzero flux.

M46.06TSolution

For F=r/∥r∥3\vv F=\vv r/\|\vv r\|^3, compute flux through a radius-RR sphere centered at the origin and explain why divergence zero off the origin does not make this flux zero.

M46.07CSolution

For F=(x,y,z)\vv F=(x,y,z), find total outward flux from the shell 1≤x2+y2+z2≤41\le x^2+y^2+z^2\le4. Separately identify outer- and inner-boundary contributions.

M46.08SSolution

For a radial field F=r h(r)\vv F=\vv r\,h(r) with r=x2+y2+z2>0r=\sqrt{x^2+y^2+z^2}>0 and differentiable hh, derive div⁡F=3h(r)+rh′(r)\diver\vv F=3h(r)+rh'(r). Find all such divergence-free radial fields and their fluxes through origin-centered spheres. Which extend continuously to the origin?

Optional small hints

For M46.02, integrate 3ρ23\rho^2 over the unit ball. For M46.04, the base normal is downward and its field component is −1-1. For M46.08, use ∂r/∂x=x/r\partial r/\partial x=x/r and solve the resulting first-order equation for hh.

H. Complete solutions: Module 46#

Compare the first step where your work differs. A corrected answer is not yet an independent success.

M46.01Show solution

Closed-surface flux from constant divergence

The field is smooth throughout the box and its neighborhood. Divergence is 33, and volume is 1⋅2⋅3=61\cdot2\cdot3=6. Thus total outward flux is 18\boxed{18}. The six faces are all included; a single face would not have this flux.

Back to M46.01
M46.02Show solution

A radial divergence makes spherical coordinates efficient

Divergence is 3x2+3y2+3z2=3ρ23x^2+3y^2+3z^2=3\rho^2. Therefore

Φ=∫02π∫0π∫013ρ4sin⁡ϕ dρ dϕ dθ=(2π)(2)(3/5)=12π/5.\Phi=\int_0^{2\pi}\int_0^\pi\int_0^1 3\rho^4\sin\phi\dd\rho\dd\phi\dd\theta =(2\pi)(2)(3/5)=\boxed{12\pi/5}.

The field is polynomial everywhere, so no interior singularity obstructs the theorem.

Back to M46.02
M46.03Show solution

Use volume rather than a complicated normal field

The ellipsoid semiaxes are 1,2,31,2,3, so its volume is (4π/3)(6)=8π(4\pi/3)(6)=8\pi. Since divergence is 33, the outward flux is 24π\boxed{24\pi}. The axis lengths are read from square-root denominators, not from the denominator values 4,94,9 themselves.

Back to M46.03
M46.04Show solution

Subtract the missing disk with its signed contribution

The cap and base disk bound E={0≤z≤4−r2,0≤r≤2}E=\{0\le z\le4-r^2,0\le r\le2\}. Its volume is 8π8\pi, and divergence of the field is 33, so closed flux is 24π24\pi. On the base z=0z=0, the outward normal is downward, giving normal field component −1-1 and base flux −4π-4\pi. Hence Φcap=24π−(−4π)=28π\Phi_{\text{cap}}=24\pi-(-4\pi)=\boxed{28\pi}. A direct check uses upward graph vector area (2x,2y,1) dA(2x,2y,1)\dd A; the dot product becomes r2+5r^2+5, whose disk integral is also 28π28\pi.

Back to M46.04
M46.05Show solution

Zero divergence means no net enclosed source for a regular field

Every component derivative of a constant field is zero. The Divergence theorem therefore gives total outward flux ∭E0 dV=0\iiint_E0\dd V=\boxed0. A patch facing the field can have positive flux and an opposing patch negative flux; the theorem asserts cancellation of the closed total, not zero flow at every boundary point.

Back to M46.05
M46.06Show solution

Compute directly when the enclosed singularity blocks the theorem

On the sphere, F=n/R2\vv F=\vv n/R^2, so the normal component is 1/R21/R^2. Multiplying by surface area 4πR24\pi R^2 gives 4π\boxed{4\pi}. The field’s divergence is zero only for r>0r>0. It is undefined at the enclosed origin, so the smoothness assumption on a neighborhood of the entire solid fails. A volume integral of its classical divergence over that solid is not a valid application of the theorem as stated.

Back to M46.06
M46.07Show solution

Outward from a shell points inward on its cavity

The shell volume is (4π/3)(23−13)=28π/3(4\pi/3)(2^3-1^3)=28\pi/3, so total flux is 33 times this, 28π\boxed{28\pi}. The radius-two outer sphere contributes 32π32\pi. On the radius-one inner sphere, the outward-from-shell normal is the negative radial unit vector, giving flux −4π-4\pi. Their sum is 28π28\pi, checking both the theorem and cavity orientation.

Back to M46.07
M46.08Show solution

Derive and solve the radial divergence equation

The first component is xh(r)xh(r), whose xx derivative is h+xh′(r)(x/r)h+xh'(r)(x/r). Summing the analogous derivatives gives

div⁡F=3h+x2+y2+z2rh′=3h+rh′.\diver\vv F=3h+\frac{x^2+y^2+z^2}{r}h'=3h+rh'.

Setting it to zero on r>0r>0 gives (r3h)′=0(r^3h)'=0, so h(r)=C/r3h(r)=C/r^3 and F=Cr/r3\boxed{\vv F=C\vv r/r^3}. On a radius-RR sphere, the normal component is C/R2C/R^2, giving flux 4πC\boxed{4\pi C}. If C≠0C\ne0, the magnitude is ∣C∣/r2|C|/r^2, unbounded near zero, so no continuous extension exists. Only C=0C=0 extends continuously, as the zero field. The nonzero flux and the singularity are linked by the derived radial structure, not a violation of the Divergence theorem.

Back to M46.08

I. Retention and readiness#

Revisit M15.08, M43.08, M44.08, and M45.06. Explain which boundary theorem fits an unfamiliar integral and list its domain and orientation assumptions before calculating. Complete the Calculus III assessment, then use its rubric to distinguish a setup gap from an arithmetic slip. The goal is independent method selection with justified conditions, not merely recognition of a familiar final formula.