22.3 Surface area calculation

Surface area calculation extends the idea of arc length for curves to two-dimensional surfaces in 3D space. By parameterizing a surface with two variables and using the cross product of partial derivatives, you can set up a double integral that computes the exact area of curved, complex shapes. This technique is foundational for flux integrals, divergence theorems, and many applications in physics and engineering.

Surface Parameterization

Representing Surfaces with Parameters

A parametric surface is defined by a vector-valued function $\mathbf{r}(u, v) = \langle x(u,v),\; y(u,v),\; z(u,v) \rangle$ that maps a region $D$ in the $uv$-plane to a surface in 3D space. The two independent parameters $u$ and $v$ let you describe complex, curved surfaces using a flat coordinate domain.

The partial derivatives $\frac{\partial \mathbf{r}}{\partial u}$ and $\frac{\partial \mathbf{r}}{\partial v}$ are tangent vectors to the surface at each point. Their cross product produces a normal vector whose magnitude measures how much the parameterization locally stretches or compresses area. This magnitude acts as the scaling factor that converts a small rectangle $du \, dv$ in parameter space into the corresponding area element on the surface.

A change of variables converts a surface integral into a double integral over $D$. The cross-product magnitude plays the same role here that the Jacobian determinant plays in substitution for multiple integrals: it accounts for the geometric distortion introduced by the parameterization.

Coordinate Systems for Surface Representation

Different coordinate systems simplify parameterization depending on the surface's symmetry.

• Cartesian coordinates $(x, y, z)$: Each point is located by distances along three perpendicular axes. Best for planes, graphs of functions $z = f(x,y)$, and surfaces without special symmetry.
• Cylindrical coordinates $(r, \theta, z)$: Here $r$ is the distance from the $z$-axis, $\theta$ is the angle in the $xy$-plane, and $z$ is the height. Natural for cylinders, cones, and any surface with circular cross-sections. A typical parameterization sets $u = \theta$ and $v = z$ (or $v = r$).
• Spherical coordinates $(\rho, \theta, \varphi)$: $\rho$ is the distance from the origin, $\theta$ is the azimuthal angle, and $\varphi$ is the polar angle from the positive $z$-axis. Ideal for spheres, ellipsoids, and surfaces with radial symmetry. For a sphere of radius $R$, you parameterize with $u = \theta$ and $v = \varphi$, giving $\mathbf{r}(\theta, \varphi) = \langle R\sin\varphi\cos\theta,\; R\sin\varphi\sin\theta,\; R\cos\varphi \rangle$.

Surface Integrals

Calculating Surface Area

The surface area of $\mathbf{r}(u,v)$ over a parameter domain $D$ is given by:

$A = \iint_D \left\|\frac{\partial \mathbf{r}}{\partial u} \times \frac{\partial \mathbf{r}}{\partial v}\right\| \, du \, dv$

Here's how to set up and evaluate this integral:

1. Write the parameterization $\mathbf{r}(u,v)$ and identify the domain $D$ with its bounds on $u$ and $v$.
2. Compute the partial derivatives $\frac{\partial \mathbf{r}}{\partial u}$ and $\frac{\partial \mathbf{r}}{\partial v}$.
3. Take their cross product $\frac{\partial \mathbf{r}}{\partial u} \times \frac{\partial \mathbf{r}}{\partial v}$. This vector is normal to the surface at each point.
4. Find the magnitude $\left\|\frac{\partial \mathbf{r}}{\partial u} \times \frac{\partial \mathbf{r}}{\partial v}\right\|$. This scalar is the differential area element $dS$.
5. Integrate over $D$ to get the total surface area.

The differential area element (scalar form) is:

$dS = \left\|\frac{\partial \mathbf{r}}{\partial u} \times \frac{\partial \mathbf{r}}{\partial v}\right\| \, du \, dv$

The vector surface element retains the direction of the normal:

$d\mathbf{S} = \left(\frac{\partial \mathbf{r}}{\partial u} \times \frac{\partial \mathbf{r}}{\partial v}\right) du \, dv$

You use $dS$ (scalar) for area calculations and scalar surface integrals, and $d\mathbf{S}$ (vector) for flux integrals.

Evaluating Surface Integrals

Scalar surface integrals integrate a function $f(x,y,z)$ over a surface:

$\iint_S f(x,y,z) \, dS = \iint_D f\!\left(\mathbf{r}(u,v)\right) \left\|\frac{\partial \mathbf{r}}{\partial u} \times \frac{\partial \mathbf{r}}{\partial v}\right\| du \, dv$

You substitute the parameterization into $f$ so that everything is expressed in terms of $u$ and $v$, then evaluate the resulting double integral over $D$.

Vector surface integrals (flux integrals) measure the flow of a vector field $\mathbf{F}$ through a surface:

$\iint_S \mathbf{F} \cdot d\mathbf{S} = \iint_D \mathbf{F}\!\left(\mathbf{r}(u,v)\right) \cdot \left(\frac{\partial \mathbf{r}}{\partial u} \times \frac{\partial \mathbf{r}}{\partial v}\right) du \, dv$

The dot product picks out the component of $\mathbf{F}$ that passes through the surface (the component along the normal). The sign of the result depends on the orientation you choose for the normal, so be consistent with the direction of the cross product.

Steps for evaluating either type:

1. Parameterize the surface and determine $D$.
2. Compute $\frac{\partial \mathbf{r}}{\partial u} \times \frac{\partial \mathbf{r}}{\partial v}$ (and its magnitude, if needed).
3. Substitute $\mathbf{r}(u,v)$ into $f$ or $\mathbf{F}$.
4. Set up the double integral over $D$ and evaluate.

Vector Calculus Fundamentals

Vector Products and Their Applications

The cross product $\mathbf{a} \times \mathbf{b}$ produces a vector perpendicular to both $\mathbf{a}$ and $\mathbf{b}$. Its magnitude equals the area of the parallelogram spanned by $\mathbf{a}$ and $\mathbf{b}$. In surface area calculations, you cross the two tangent vectors $\frac{\partial \mathbf{r}}{\partial u}$ and $\frac{\partial \mathbf{r}}{\partial v}$ to get both the normal direction and the local area scaling.

The dot product $\mathbf{a} \cdot \mathbf{b}$ returns a scalar equal to $\|\mathbf{a}\|\|\mathbf{b}\|\cos\theta$, where $\theta$ is the angle between the vectors. In flux integrals, the dot product of $\mathbf{F}$ with the surface element $d\mathbf{S}$ extracts how much of the field passes through the surface rather than along it.

These two products work together in surface integrals: the cross product builds the surface element, and the dot product projects the vector field onto it.

Gradient, Divergence, and Curl

These differential operators connect to surface integrals through the major theorems of vector calculus.

• Gradient $\nabla f = \left\langle \frac{\partial f}{\partial x},\; \frac{\partial f}{\partial y},\; \frac{\partial f}{\partial z}\right\rangle$: Points in the direction of steepest increase of $f$, with magnitude equal to the rate of increase. For a surface defined implicitly by $F(x,y,z) = c$, the gradient $\nabla F$ is normal to the surface, which provides another route to finding normal vectors.
• Divergence $\nabla \cdot \mathbf{F} = \frac{\partial F_x}{\partial x} + \frac{\partial F_y}{\partial y} + \frac{\partial F_z}{\partial z}$: A scalar that measures the net outward flux per unit volume. The Divergence Theorem relates a flux integral over a closed surface to a triple integral of divergence over the enclosed volume.
• Curl $\nabla \times \mathbf{F} = \left\langle \frac{\partial F_z}{\partial y} - \frac{\partial F_y}{\partial z},\; \frac{\partial F_x}{\partial z} - \frac{\partial F_z}{\partial x},\; \frac{\partial F_y}{\partial x} - \frac{\partial F_x}{\partial y}\right\rangle$: A vector that captures the local rotational tendency of $\mathbf{F}$. Stokes' Theorem relates the flux of the curl through a surface to the circulation of $\mathbf{F}$ around the surface's boundary.

Understanding these operators now will pay off directly when you reach the Divergence Theorem and Stokes' Theorem, where surface area computation and flux integrals are the central tools.