# Surface integrals ​

See also the graded exercises and additional exercises associated to this part of the course. If you want more, Sections 6.4 to 6.8 of OpenStax Calculus Volume 3 are a good option.

In this section we consider surfaces and how to define integral of vector fields over these surfaces. This is similar in many ways to line integrals but a higher dimensional version. Curves (for line integrals) are 1D subsets of higher dimensional space whereas surfaces are 2D subsets of higher dimensional space. Identically to line integrals, the first step is to understand a practical way to represent the surfaces, just like with curves we used paths as the parametric representation of the curve. Once we have clarified the parametric representation of surface we can define the surface integral (of a vector field) and show that it satisfies various properties which we would expect, including that the integral is independent of the choice of parametrization. Similar to how we were able to use a line integral (of a scalar) to calculate the length of a curve we can use a surface integral (of a scalar) to calculate the area of a surface.

We then introduce two important operators that act on vector fields, namely curl and divergence. Using these operators and the surface integral we introduce two theorems, Gauss' Theorem and Stokes' Theorem. These theorems connect line integrals with surface integrals and with volume integrals.

## Representation of a surface ​

Before developing parametric representations of surfaces let's recall an example of parametric representation of a curve (path). For example, the half circle $C=\left\{\left(x,y\right):{x}^{2}+{y}^{2}=1,y\ge 1\right\}$ can be parametrized in many ways, including the following two paths.

$\begin{array}{rl}\mathbit{\alpha }\left(x\right)& =\left(x,\sqrt{1-{x}^{2}}\right),\phantom{\rule{1em}{0ex}}x\in \left[-1,1\right],\\ \mathbit{\alpha }\left(t\right)& =\left(\mathrm{cos}t,\mathrm{sin}t\right),\phantom{\rule{1em}{0ex}}t\in \left[0,\pi \right].\end{array}$

In a similar way, now in 2D we can have a parametric representation of a hemisphere.

The hemisphere $S=\left\{\left(x,y,z\right):{x}^{2}+{y}^{2}+{z}^{2}=1,z\ge 0\right\}$ can be represented parametrically in many ways, including

$\begin{array}{rl}\mathbf{r}\left(x,y\right)& =\left(x,y,\sqrt{1-{x}^{2}-{y}^{2}}\right),\phantom{\rule{1em}{0ex}}\left(x,y\right)\in \left\{{x}^{2}+{y}^{2}\le 1\right\},\\ \mathbf{r}\left(u,v\right)& =\left(\mathrm{cos}u\mathrm{cos}v,\mathrm{sin}u\mathrm{cos}v,\mathrm{sin}v\right),\phantom{\rule{1em}{0ex}}\left(u,v\right)\in \left[0,2\pi \right]×\left[0,\pi /2\right].\end{array}$

Observe that the second form above can be deduced from spherical coordinates (fixed distance from the origin).

The cone $S=\left\{\left(x,y,z\right):{z}^{2}={x}^{2}+{y}^{2},z\in \left[0,1\right]\right\}$ can be represented parametrically in many ways, including

$\begin{array}{rl}\mathbf{r}\left(x,y\right)& =\left(x,y,\sqrt{{x}^{2}+{y}^{2}}\right),\phantom{\rule{1em}{0ex}}\left(x,y\right)\in \left\{{x}^{2}+{y}^{2}\le 1\right\},\\ \mathbf{r}\left(u,v\right)& =\left(v\mathrm{cos}u,v\mathrm{sin}u,v\right),\phantom{\rule{1em}{0ex}}\left(u,v\right)\in \left[0,2\pi \right]×\left[0,1\right].\end{array}$

Observe that the second form can be deduced from spherical coordinates (fixed angle from $z$-axis) or from cylindrical coordinates (fixed ratio between vertical coordinate and distance from the the $z$-axis).

### Fundamental vector product ​

A key notion for parametric surfaces and natural geometric object is the fundamental vector product. Consider the parametric surface, denoted $\mathbf{r}\left(T\right)$, and suppose it has the form

$\mathbf{r}\left(u,v\right)=\left(X\left(u,v\right),Y\left(u,v\right),Z\left(u,v\right)\right),\phantom{\rule{1em}{0ex}}\left(u,v\right)\in T.$

The vector-valued function defined as

$\frac{\partial \mathbf{r}}{\partial u}×\frac{\partial \mathbf{r}}{\partial v}=\left(\begin{array}{c}{\partial }_{u}X\\ {\partial }_{u}Y\\ {\partial }_{u}Z\end{array}\right)×\left(\begin{array}{c}{\partial }_{v}X\\ {\partial }_{v}Y\\ {\partial }_{v}Z\end{array}\right)$

is called the fundamental vector product of the representation $\mathbf{r}$.

By definition, the vector-valued functions $\frac{\partial \mathbf{r}}{\partial u}$ and $\frac{\partial \mathbf{r}}{\partial v}$ are tangent to the surface. As such, assuming that they are linearly independent, the fundamental vector product $\frac{\partial \mathbf{r}}{\partial u}×\frac{\partial \mathbf{r}}{\partial v}$ is normal to the surface (orthogonal to every curve which passes through the surface). Moreover the norm of the vector represents the local scaling of area (small parallelograms).

As always we need to take some care about smoothness of the objects we work with. If $\left(u,v\right)$ is a point in $T$ at which $\frac{\partial \mathbf{r}}{\partial u}$ and $\frac{\partial \mathbf{r}}{\partial v}$ are continuous and the fundamental vector product is non-zero then $\mathbf{r}\left(u,v\right)$ is said to be a regular point for that representation.

A surface $\mathbf{r}\left(T\right)$ is said to be smooth if all its points are regular points.

Just like we saw with paths to represent curves, there are many different ways we can find the parametric representation of a given surface. If the surface $S$ has the form $z=f\left(x,y\right)$ (the surface in written in explicit form) then we can use $x,y$ as the parameters and have the representation

$\mathbf{r}\left(x,y\right)=\left(x,y,f\left(x,y\right)\right),\phantom{\rule{1em}{0ex}}\left(x,y\right)\in T.$

The region $T$ is the projection of $S$ onto the $xy$-plane. For such a surface we compute

$\frac{\partial \mathbf{r}}{\partial x}=\left(\begin{array}{c}1\\ 0\\ {\partial }_{x}f\end{array}\right),\phantom{\rule{1em}{0ex}}\frac{\partial \mathbf{r}}{\partial y}=\left(\begin{array}{c}0\\ 1\\ {\partial }_{y}f\end{array}\right),$

and consequently

$\frac{\partial \mathbf{r}}{\partial x}×\frac{\partial \mathbf{r}}{\partial y}=\left(\begin{array}{c}1\\ 0\\ {\partial }_{x}f\end{array}\right)×\left(\begin{array}{c}0\\ 1\\ {\partial }_{y}f\end{array}\right)=\left(\begin{array}{c}-{\partial }_{x}f\\ -{\partial }_{y}f\\ 1\end{array}\right).$

An example of such a representation is as follows for the hemisphere. Let $T=\left\{{x}^{2}+{y}^{2}\le 1\right\}$, and let

$\mathbf{r}\left(x,y\right)=\left(x,y,\sqrt{1-{x}^{2}-{y}^{2}}\right).$

The surface $\mathbf{r}\left(T\right)$ is the unit hemisphere $\left\{\left(x,y,z\right):{x}^{2}+{y}^{2}+{z}^{2}=1,z\ge 0\right\}$. The fundamental vector product of this representation is

In this case, all points are regular except the equator $\left\{\left(x,y,0\right):{x}^{2}+{y}^{2}=1\right\}$.

Let $T=\left[0,2\pi \right]×\left[0,\pi /2\right]$ and let

$\mathbf{r}\left(u,v\right)=\left(\mathrm{cos}u\mathrm{cos}v,\mathrm{sin}u\mathrm{cos}v,\mathrm{sin}v\right).$

The surface $\mathbf{r}\left(T\right)$ is the unit hemisphere $\left\{\left(x,y,z\right):{x}^{2}+{y}^{2}+{z}^{2}=1\right\}$. This representation is connected to spherical coordinates. We calculate that

$\frac{\partial \mathbf{r}}{\partial u}\left(u,v\right)=\left(\begin{array}{c}-\mathrm{sin}u\mathrm{cos}v\\ \mathrm{cos}u\mathrm{cos}v\\ 0\end{array}\right),\phantom{\rule{1em}{0ex}}\frac{\partial \mathbf{r}}{\partial v}\left(u,v\right)=\left(\begin{array}{c}-\mathrm{cos}u\mathrm{sin}v\\ -\mathrm{sin}u\mathrm{sin}v\\ \mathrm{cos}v\end{array}\right),$

and so the fundamental vector product of this representation is

In this case many points map to the north pole $\left(0,0,1\right)$ and so the north pole is not a regular point. Additionally there are two points which map to each point on the line between equator and north pole $\left\{\left(x,y,z\right)\in \mathbf{r}\left(T\right):x\ge 0,y=0\right\}$.

## Surface integral of scalar field ​

Mirroring the process for line integrals we will define surface integrals both for scalar fields and for vector fields. The surface integral of a scalar field is closely related to the area of a parametric surface, just like the length of a curve is closely related to the line integral of a scalar field.

Definition

The area of the parametric surface $S=\mathbf{r}\left(T\right)$ is defined as the double integral

Observe that the definition is in terms of a multiple integral over the region $T$, and the quantity being integrated is the norm of the fundamental vector product.

Later we will show that $\mathrm{Area}\left(S\right)$ is independent of the choice of representation as we require for such a definition, it would be unreasonable if the area of a surface depended on the choice of representation.

We will check that this definition corresponds to a fact that we already know by computing the surface area of a hemisphere. Let, as before, $T=\left[0,2\pi \right]×\left[0,\pi /2\right]$ and let $\mathbf{r}\left(u,v\right)=\left(\mathrm{cos}u\mathrm{cos}v,\mathrm{sin}u\mathrm{cos}v,\mathrm{sin}v\right)$. The norm of the fundamental vector product (which we computed earlier) is

Taking the definition of area and evaluating the multiple integral, this means that

The surface integral of a scalar field is defined in a way similar to the area of a surface.

Definition (scalar surface integral)

Let $S=\mathbf{r}\left(T\right)$ be a parametric surface and let $f$ be a scalar field defined on $S$. The surface integral of $f$ over $S$ is defined as

whenever the double integral on the right exists.

Observe that, if we choose $f\equiv 1$, that is we choose the scalar field identically equal to $1$, then we obtain the formula for the area of the surface). This is just the same as the line integral of a scalar and the length of the corresponding curve.

From the point of view of applications, we could take $f$ as the density of thin material which has the shape of the surface $S$ and then is the total mass of this piece of material. Extending this idea we could also calculate the centre of mass of this piece of material.

## Change of surface parametrization ​

In order to validate the definition of a surface integral and consequently that of the area of a surface, we will now show that the the value of the evaluated integral doesn't depend on the choice of representation for any given surface.

Suppose that $\mathbf{q}\left(A\right)$ and $\mathbf{r}\left(B\right)$ are both representations of the same surface, and that $\mathbf{r}=\mathbf{q}\circ G$ for some differentiable $G:B\to A$. Then

Since $\mathbf{r}\left(u,v\right)=\mathbf{q}\left(S\left(u,v\right),T\left(u,v\right)\right)$ we calculate (chain rule and vector product) that

$\left[\frac{\partial \mathbf{r}}{\partial u}×\frac{\partial \mathbf{r}}{\partial v}\right]\left(u,v\right)=\left[\left(\frac{\partial \mathbf{q}}{\partial s}×\frac{\partial \mathbf{q}}{\partial t}\right)\left(\frac{\partial S}{\partial u}\frac{\partial T}{\partial v}-\frac{\partial S}{\partial v}\frac{\partial T}{\partial u}\right)\right]\left(S\left(u,v\right),T\left(u,v\right)\right).$

Observe that $\frac{\partial S}{\partial u}\frac{\partial T}{\partial v}-\frac{\partial S}{\partial v}\frac{\partial T}{\partial u}$ is the Jacobian determinant associated to change of variables $\left(u,v\right)↦\left(S\left(u,v\right),T\left(u,v\right)\right)$. Consequently, by the change of variables theorem,

so the definition does make sense.

## Surface integral of a vector field ​

In preparation for defining the surface integral of a vector field we need the notion of the normal vector of a surface. This is a natural geometric notion, for each point in the surface it is the unit vector field which is orthogonal to the surface.

Definition (unit normal)

Let $S=\mathbf{r}\left(T\right)$ be a parametric surface. At each regular point the two unit normals are

${\mathbf{n}}_{1}=\frac{\frac{\partial \mathbf{r}}{\partial u}×\frac{\partial \mathbf{r}}{\partial v}}{‖\frac{\partial \mathbf{r}}{\partial u}×\frac{\partial \mathbf{r}}{\partial v}‖}\phantom{\rule{1em}{0ex}}\text{and}\phantom{\rule{1em}{0ex}}{\mathbf{n}}_{2}=-{\mathbf{n}}_{1}.$

This definition makes $‖{\mathbf{n}}_{1}‖=‖{\mathbf{n}}_{2}‖=1$. That there are two normal vectors is expected because there are two sides to the surface at each point, one is just the opposite direction to the other. When we have two parameterizations of the same surface, they always have the same pair of normals at any regular point, but which one is ${\mathbf{n}}_{1}$ and which one is ${\mathbf{n}}_{2}$ can be different.

If $\mathbf{f}$ is a vector field then $\mathbf{f}\cdot \mathbf{n}$ is the component of the flow in direction of $\mathbf{n}$.

Definition (vector surface integral)

Let $S=\mathbf{r}\left(T\right)$ be a parametric surface and $\mathbf{f}$ a vector field. The integral

is said to be the surface integral of $\mathbf{f}$ with respect to the normal $\mathbf{n}$.

For convenience let $\mathbf{N}=\frac{\partial \mathbf{r}}{\partial u}×\frac{\partial \mathbf{r}}{\partial v}$ and $\mathbf{n}=\mathbf{N}/‖\mathbf{N}‖$. Observe that

and so for evaluating the surface integral of a vector field there is typically no need to evaluate the norm of the fundamental vector product. Also note that because ${\mathbf{n}}_{1}=-{\mathbf{n}}_{2}$. This means that chosing one normal or the other simply corresponds to a minus sign in the evaluated integral. This is the notion that there is a choice of orientation inherent with a surface. As a tangible example imagine that the surface has a flow passing it and this flow is determined by a vector field. Then the surface integral would represent the total flow passing the given surface in a given direction.

## Curl and divergence ​

Suppose that $\mathbf{f}=\left(\begin{array}{c}{f}_{x}\\ {f}_{y}\\ {f}_{z}\end{array}\right)$ is a differentiable vector field.

Definition (curl)

The curl of $\mathbf{f}$ is defined as

$\mathrm{\nabla }×\mathbf{f}=\left(\begin{array}{c}\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}\end{array}\right).$

Definition (divergence)

The divergence of $\mathbf{f}$ is defined as

$\mathrm{\nabla }\cdot \mathbf{f}=\frac{\partial {f}_{x}}{\partial x}+\frac{\partial {f}_{y}}{\partial y}+\frac{\partial {f}_{z}}{\partial z}.$

Often the notation $\mathrm{curl}\mathbf{f}=\mathrm{\nabla }×\mathbf{f}$ and $\mathrm{div}\mathbf{f}=\mathrm{\nabla }\cdot \mathbf{f}$ is used instead. Note that the symbols"$×$" and "$\cdot$" used in the notation for curl and divergence are not truly representing the vector and scalar product but are more a convenient way to remember the definitions. These quantities satisfy the following basic properties which can all be proved by basic calculations.

• If $\mathbf{f}=\mathrm{\nabla }\phi$ then $\mathrm{\nabla }×\mathbf{f}=\mathbf{0}$,
• $\mathrm{\nabla }\cdot \left(\mathrm{\nabla }×\mathbf{f}\right)=0$,
• $\mathrm{\nabla }×\left(\mathrm{\nabla }×\mathbf{f}\right)=\mathrm{\nabla }\left(\mathrm{\nabla }\cdot \mathbf{f}\right)-{\mathrm{\nabla }}^{2}\mathbf{f}$.

The quantity defined as ${\mathrm{\nabla }}^{2}\phi =\mathrm{\nabla }\cdot \left(\mathrm{\nabla }\phi \right)=\frac{{\partial }^{2}\phi }{\partial {x}^{2}}+\frac{{\partial }^{2}\phi }{\partial {y}^{2}}+\frac{{\partial }^{2}\phi }{\partial {z}^{2}}$ is called the Laplacian and occurs in many applications of physics and mathematics.

Some examples:

If