5.4 Surface Integrals, Stokes' Theorem, and the Divergence Theorem

5.4 Surface Integrals, Stokes' Theorem, and the Divergence Theorem

Surface integrals and the classical integral theorems of Stokes and Gauss represent the culmination of multivariable calculus. On the GRE Subject Test in Mathematics, questions frequently evaluate your ability to simplify intractable surface flux calculations into trivial volume integrals or loop circulations, exploiting topological boundary relationships and vector identities.

---

Parametric Surfaces and Surface Area

A two-dimensional surface $S \subset \mathbb{R}^3$ can be parametrized by a vector-valued function of two parameters: $\mathbf{r}(u, v) = \langle x(u, v), \; y(u, v), \; z(u, v) \rangle, \quad (u, v) \in D \subset \mathbb{R}^2.$

Tangent Vectors and the Fundamental Normal

The tangent vectors along the coordinate grid lines are: $\mathbf{r}_u = \frac{\partial \mathbf{r}}{\partial u} = \left\langle \frac{\partial x}{\partial u}, \frac{\partial y}{\partial u}, \frac{\partial z}{\partial u} \right\rangle, \quad \mathbf{r}_v = \frac{\partial \mathbf{r}}{\partial v} = \left\langle \frac{\partial x}{\partial v}, \frac{\partial y}{\partial v}, \frac{\partial z}{\partial v} \right\rangle.$ The surface is regular (smooth) at points where the cross product $\mathbf{r}_u \times \mathbf{r}_v \neq \mathbf{0}$. The vector $\mathbf{N} = \mathbf{r}_u \times \mathbf{r}_v$ is the standard upward normal vector.

Surface Area Element

The differential surface area element is the area of the infinitesimal tangent parallelogram: $dS = \|\mathbf{r}_u \times \mathbf{r}_v\| \, du \, dv.$

Explicit Surfaces $z = g(x, y)$

When $S$ is given as the graph of a function $z = g(x, y)$, choose parameters $u = x$ and $v = y$: $\mathbf{r}(x, y) = \langle x, y, g(x, y) \rangle \implies \mathbf{r}_x \times \mathbf{r}_y = \langle -g_x, -g_y, 1 \rangle.$ The norm yields the standard graph area formula: $dS = \sqrt{1 + \left( \frac{\partial g}{\partial x} \right)^2 + \left( \frac{\partial g}{\partial y} \right)^2} \, dx \, dy.$

Spherical Parameterization

For a sphere of fixed radius $R$ parametrized by colatitude $\varphi$ and azimuth $\theta$: $\mathbf{r}(\varphi, \theta) = \langle R \sin \varphi \cos \theta, \; R \sin \varphi \sin \theta, \; R \cos \varphi \rangle \implies \|\mathbf{r}_\varphi \times \mathbf{r}_\theta\| = R^2 \sin \varphi.$ Hence $dS = R^2 \sin \varphi , d\varphi , d\theta$.

---

Surface Integrals of Scalar and Vector Fields (Flux)

Scalar Surface Integrals

For a continuous scalar field $f(x, y, z)$ on $S$: $\iint_S f(x, y, z) \, dS = \iint_D f(\mathbf{r}(u, v)) \|\mathbf{r}_u \times \mathbf{r}_v\| \, du \, dv.$

Oriented Surfaces and Flux Integrals

A surface $S$ is orientable if there exists a continuous unit normal vector field $\mathbf{n}(x, y, z)$ defined over all of $S$ (excluding one-sided non-orientable surfaces such as the Möbius strip). The flux of a vector field $\mathbf{F}$ through an oriented surface $S$ is: $\iint_S \mathbf{F} \cdot d\mathbf{S} = \iint_S (\mathbf{F} \cdot \mathbf{n}) \, dS.$ Substituting $\mathbf{n} = \frac{\mathbf{r}_u \times \mathbf{r}_v}{|\mathbf{r}_u \times \mathbf{r}_v|}$ and $dS = |\mathbf{r}_u \times \mathbf{r}_v| , du , dv$, the norms cancel identically: $\iint_S \mathbf{F} \cdot d\mathbf{S} = \iint_D \mathbf{F}(\mathbf{r}(u, v)) \cdot (\mathbf{r}_u \times \mathbf{r}_v) \, du \, dv.$ For an explicit surface $z = g(x, y)$ oriented upward (normal with positive $z$-component): $\iint_S \mathbf{F} \cdot d\mathbf{S} = \iint_D (-P g_x - Q g_y + R) \, dx \, dy.$

---

Stokes' Theorem

Stokes' Theorem equates the circulation of a vector field around the boundary of a surface to the flux of its curl across the surface.

Formal Statement

Let $S$ be an oriented, piecewise-smooth surface bounded by a simple, closed, piecewise-smooth boundary curve $C = \partial S$. Let the orientation of $C$ be compatible with that of $S$ via the right-hand rule (if the fingers of your right hand curl in the direction of traversal along $C$, your thumb points in the direction of the normal $\mathbf{n}$). If $\mathbf{F}$ is a $C^1$ vector field: $\oint_{\partial S} \mathbf{F} \cdot d\mathbf{r} = \iint_S (\operatorname{curl} \mathbf{F}) \cdot \mathbf{n} \, dS = \iint_S (\nabla \times \mathbf{F}) \cdot d\mathbf{S}.$

Surface Independence of Curl Flux

A paramount consequence on the GRE: The flux $\iint_S (\nabla \times \mathbf{F}) \cdot d\mathbf{S}$ depends solely on the boundary curve $\partial S$, not on the interior shape of the surface $S$! If $S_1$ and $S_2$ share the same oriented boundary curve $\partial S_1 = \partial S_2$: $\iint_{S_1} (\nabla \times \mathbf{F}) \cdot d\mathbf{S} = \iint_{S_2} (\nabla \times \mathbf{F}) \cdot d\mathbf{S}.$ Whenever faced with an intractable surface (e.g., a wavy paraboloid or high-degree polynomial dome), replace it with the flat planar disk or polygon spanning the same boundary!

Closed Surfaces

If $S$ is a closed surface (such as a full sphere or torus), it has no boundary ($\partial S = \emptyset$). Therefore: $\iint_S (\nabla \times \mathbf{F}) \cdot d\mathbf{S} = 0.$

---

Gauss's Divergence Theorem

The Divergence Theorem (Gauss's Theorem) relates the outward flux of a vector field across a closed surface to the volume integral of its divergence over the enclosed region.

Formal Statement

Let $E \subset \mathbb{R}^3$ be a bounded solid region whose boundary $S = \partial E$ is a closed, piecewise-smooth surface oriented with outward-pointing unit normal vectors $\mathbf{n}$. If $\mathbf{F}$ is a $C^1$ vector field on an open region containing $E$: $\iint_{\partial E} \mathbf{F} \cdot \mathbf{n} \, dS = \iiint_E (\operatorname{div} \mathbf{F}) \, dV = \iiint_E (\nabla \cdot \mathbf{F}) \, dV.$

The Capping Technique for Open Surfaces

When asked for the flux through an open surface $S$ (e.g., an open cylinder, paraboloid bowl, or hemisphere without a base):

1. The Divergence Theorem requires a closed surface. Create one by attaching a simple planar "cap" $S_{\text{cap}}$ (typically a flat disk on a coordinate plane) such that $S \cup S_{\text{cap}} = \partial E$.
2. Apply the Divergence Theorem to the combined closed boundary: $\iint_S \mathbf{F} \cdot \mathbf{n} \, dS + \iint_{S_{\text{cap}}} \mathbf{F} \cdot \mathbf{n}_{\text{cap}} \, dS = \iiint_E (\nabla \cdot \mathbf{F}) \, dV.$
3. Solve for the desired flux: $\iint_S \mathbf{F} \cdot \mathbf{n} \, dS = \iiint_E (\nabla \cdot \mathbf{F}) \, dV - \iint_{S_{\text{cap}}} \mathbf{F} \cdot \mathbf{n}_{\text{cap}} \, dS.$ Because $\mathbf{n}_{\text{cap}}$ is typically a constant coordinate unit vector (such as $-\mathbf{k}$) and a coordinate is constant on the cap (such as $z = 0$), the cap flux integral is almost always trivial.

---

Unified Hierarchy of Integral Theorems

All fundamental theorems of multivariable calculus are specific incarnations of the Generalized Stokes' Theorem for differential forms: $\int_\Omega d\omega = \int_{\partial \Omega} \omega$.

---

Step-by-Step Worked Problems

Problem 1: Stokes' Theorem Surface Replacement

Let $C$ be the circle $x^2 + y^2 = 4$ in the plane $z = 3$, oriented counterclockwise when viewed from above. Evaluate $\oint_C \mathbf{F} \cdot d\mathbf{r}$ for $\mathbf{F} = \langle -y^3, ; x^3, ; z^3 \rangle$.

Solution:

1. By Stokes' Theorem, $\oint_C \mathbf{F} \cdot d\mathbf{r} = \iint_S (\nabla \times \mathbf{F}) \cdot \mathbf{n} , dS$, where $S$ is any oriented surface bounded by $C$.
2. Choose the flat circular disk $S = {(x, y, 3) : x^2 + y^2 \le 4}$ lying in the plane $z = 3$. The upward normal is $\mathbf{n} = \mathbf{k} = \langle 0, 0, 1 \rangle$.
3. Compute $\nabla \times \mathbf{F}$: $\nabla \times \mathbf{F} = \det \begin{pmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \\ \partial_x & \partial_y & \partial_z \\ -y^3 & x^3 & z^3 \end{pmatrix} = \mathbf{i}(0) - \mathbf{j}(0) + \mathbf{k}\left( \frac{\partial}{\partial x}(x^3) - \frac{\partial}{\partial y}(-y^3) \right) = \langle 0, 0, 3x^2 + 3y^2 \rangle.$
4. Dot with normal: $(\nabla \times \mathbf{F}) \cdot \mathbf{n} = 3(x^2 + y^2)$.
5. Convert the double integral over the disk of radius 2 to polar coordinates: $\iint_S 3(x^2 + y^2) \, dA = 3 \int_0^{2\pi} d\theta \int_0^2 r^2 (r \, dr) = 3(2\pi) \left[ \frac{r^4}{4} \right]_0^2 = 6\pi \left( \frac{16}{4} \right) = 24\pi.$

Problem 2: Divergence Theorem with Polynomial Vector Field

Compute the outward flux of $\mathbf{F} = \langle x^3, ; y^3, ; z^3 \rangle$ across the sphere $S: x^2 + y^2 + z^2 = R^2$.

Solution:

1. $S$ is the closed boundary of the solid ball $E = {\mathbf{x} : |\mathbf{x}| \le R}$.
2. Compute the divergence of $\mathbf{F}$: $\nabla \cdot \mathbf{F} = \frac{\partial}{\partial x}(x^3) + \frac{\partial}{\partial y}(y^3) + \frac{\partial}{\partial z}(z^3) = 3(x^2 + y^2 + z^2).$
3. By the Divergence Theorem: $\iint_S \mathbf{F} \cdot \mathbf{n} \, dS = \iiint_E 3(x^2 + y^2 + z^2) \, dV.$
4. Evaluate in spherical coordinates, where $x^2 + y^2 + z^2 = \rho^2$ and $dV = \rho^2 \sin \varphi , d\rho , d\varphi , d\theta$: $\iiint_E 3\rho^2 (\rho^2 \sin \varphi \, d\rho \, d\varphi \, d\theta) = 3 \left( \int_0^{2\pi} d\theta \right) \left( \int_0^\pi \sin \varphi \, d\varphi \right) \left( \int_0^R \rho^4 \, d\rho \right).$
5. Multiply the components: $3 \times (2\pi) \times (2) \times \left( \frac{R^5}{5} \right) = \frac{12\pi R^5}{5}.$

---

GRE Exam Traps & Pitfalls

Trap 1: Forgetting to Check if Surface is Closed in Divergence Theorem The Divergence Theorem applies exclusively to CLOSED surfaces. Applying it directly to an open hemisphere or cylinder without including the base cap is an immediate point loss.

Trap 2: Normal Orientation and Boundary Consistency in Stokes' Theorem Ensure the boundary traversal matches the surface normal via the right-hand rule. If your thumb points along $\mathbf{n}$, your fingers define the positive direction of $\partial S$. If the problem specifies the opposite traversal, multiply the flux by $-1$.

Trap 3: Inward vs. Outward Flux in Gauss's Theorem By convention, Gauss's Divergence Theorem assumes outward-pointing unit normal vectors. If a question requests "inward flux", the result is $-\iiint_E (\nabla \cdot \mathbf{F}) , dV$.