#### Definition 16.6.2. Surface Integral.

Let \(G(x,y,z)\) be a continuous function defined on a surface \(\surfaceS\text{.}\) The surface integral of \(G\) on \(\surfaceS\) is

\begin{equation*}
\iint_\surfaceS G(x,y,z)\, dS\text{.}
\end{equation*}

Consider a smooth surface \(\surfaceS\) that represents a thin sheet of metal. How could we find the mass of this metallic object?

If the density of this object is constant, then we can find mass via “mass\(=\) density × surface area,” and we could compute the surface area using the techniques of the previous section.

What if the density were not constant, but variable, described by a function \(\delta(x,y,z)\text{?}\) We can describe the mass using our general integration techniques as

\begin{equation*}
\text{ mass } = \iint_\surfaceS \, dm\text{,}
\end{equation*}

where \(dm\) represents “a little bit of mass.” That is, to find the total mass of the object, sum up lots of little masses over the surface.

How do we find the “little bit of mass” \(dm\text{?}\) On a small portion of the surface with surface area \(\Delta S\text{,}\) the density is approximately constant, hence \(dm \approx \delta(x,y,z)\Delta S\text{.}\) As we use limits to shrink the size of \(\Delta S\) to 0, we get \(dm = \delta(x,y,z)dS\text{;}\) that is, a little bit of mass is equal to a density times a small amount of surface area. Thus the total mass of the thin sheet is

\begin{equation}
\text{ mass } =\iint_\surfaceS \delta(x,y,z)\, dS\text{.}\tag{16.6.1}
\end{equation}

To evaluate the above integral, we would seek \(\vec r(u,v)\text{,}\) a smooth parametrization of \(\surfaceS\) over a region \(R\) of the \(u\)-\(v\) plane. The density would become a function of \(u\) and \(v\text{,}\) and we would integrate \(\iint_R \delta(u,v)\snorm{\vec r_u\times \vec r_v}\, dA\text{.}\)

The integral in Equation (16.6.1) is a specific example of a more general construction defined below.

Let \(G(x,y,z)\) be a continuous function defined on a surface \(\surfaceS\text{.}\) The surface integral of \(G\) on \(\surfaceS\) is

\begin{equation*}
\iint_\surfaceS G(x,y,z)\, dS\text{.}
\end{equation*}

Surface integrals can be used to measure a variety of quantities beyond mass. If \(G(x,y,z)\) measures the static charge density at a point, then the surface integral will compute the total static charge of the sheet. If \(G\) measures the amount of fluid passing through a screen (represented by \(\surfaceS\)) at a point, then the surface integral gives the total amount of fluid going through the screen.

Find the mass of a thin sheet modeled by the plane \(2x+y+z=3\) over the triangular region of the \(xy\)-plane bounded by the coordinate axes and the line \(y=2-2x\text{,}\) as shown in Figure 16.6.4, with density function \(\delta(x,y,z) = x^2+5y+z\text{,}\) where all distances are measured in cm and the density is given as gm/cm\(^2\text{.}\)

We begin by parametrizing the planar surface \(\surfaceS\text{.}\) Using the techniques of the previous section, we can let \(x=u\) and \(y=v(2-2u)\text{,}\) where \(0\leq u\leq 1\) and \(0\leq v\leq 1\text{.}\) Solving for \(z\) in the equation of the plane, we have \(z=3-2x-y\text{,}\) hence \(z = 3-2u-v(2-2u)\text{,}\) giving the parametrization \(\vec r(u,v) = \langle u, v(2-2u), 3-2u-v(2-2u)\rangle\text{.}\)

We need \(dS=\snorm{\vec r_u\times \vec r_v}dA\text{,}\) so we need to compute \(\vec r_u\text{,}\) \(\vec r_v\) and the norm of their cross product. We leave it to the reader to confirm the following:

\begin{equation*}
\vec r_u = \langle 1,-2v,2v-2\rangle, \vec r_v = \langle 0,2-2u, 2u-2\rangle\text{,}
\end{equation*}

\begin{equation*}
\vec r_u\times \vec r_v = \langle 4-4u,2-2u,2-2u\rangle \text{ and } \snorm{\vec r_u\times \vec r_v} = 2\sqrt{6}\sqrt{(u-1)^2}\text{.}
\end{equation*}

We need to be careful to not “simplify” \(\snorm{\vec r_u\times \vec r_v} = 2\sqrt{6}\sqrt{(u-1)^2}\) as \(2\sqrt{6}(u-1)\text{;}\) rather, it is \(2\sqrt{6}|u-1|\text{.}\) In this example, \(u\) is bounded by \(0\leq u\leq 1\text{,}\) and on this interval \(|u-1| = 1-u\text{.}\) Thus \(dS = 2\sqrt{6}(1-u)dA\text{.}\)

The density is given as a function of \(x\text{,}\) \(y\) and \(z\text{,}\) for which we’ll substitute the corresponding components of \(\vec r\) (with the slight abuse of notation that we used in previous sections):

\begin{align*}
\delta(x,y,z) \amp = \delta\big(\vec r(u,v)\big)\\
\amp = u^2 + 5v(2-2u)+3-2u-v(2-2u)\\
\amp = u^2-8uv-2u+8v+3\text{.}
\end{align*}

Thus the mass of the sheet is:

\begin{align*}
M \amp = \iint_\surfaceS\, dm\\
\amp = \iint_R \delta\big(\vec r(u,v)\big)\snorm{\vec r_u\times \vec r_v}\, dA\\
\amp = \int_0^1\int_0^1 \big(u^2-8uv-2u+8v+3\big)\big(2\sqrt{6}(1-u)\big)\, du\, dv\\
\amp = \frac{31}{\sqrt{6}} \approx 12.66 \text{ gm. }
\end{align*}

Let a surface \(\surfaceS\) lie within a vector field \(\vec F\text{.}\) One is often interested in measuring the *flux* of \(\vec F\) across \(\surfaceS\text{;}\) that is, measuring “how much of the vector field passes across \(\surfaceS\text{.}\)” For instance, if \(\vec F\) represents the velocity field of moving air and \(\surfaceS\) represents the shape of an air filter, the flux will measure how much air is passing through the filter per unit time.

As flux measures the amount of \(\vec F\) passing across \(\surfaceS\text{,}\) we need to find the “amount of \(\vec F\) orthogonal to \(\surfaceS\text{.}\)” Similar to our measure of flux in the plane, this is equal to \(\vec F\cdot \vec n\text{,}\) where \(\vec n\) is a unit vector normal to \(\surfaceS\) at a point. We now consider how to find \(\vec n\text{.}\)

Given a smooth parametrization \(\vec r(u,v)\) of \(\surfaceS\text{,}\) the work in the previous section showing the development of our method of computing surface area also shows that \(\vec r_u(u,v)\) and \(\vec r_v(u,v)\) are tangent to \(\surfaceS\) at \(\vec r(u,v)\text{.}\) Thus \(\vec r_u\times \vec r_v\) is orthogonal to \(\surfaceS\text{,}\) and we let

\begin{equation*}
\vec n = \frac{\vec r_u\times \vec r_v}{\snorm{\vec r_u\times \vec r_v}}\text{,}
\end{equation*}

which is a unit vector normal to \(\surfaceS\) at \(\vec r(u,v)\text{.}\)

The measurement of flux across a surface is a surface integral; that is, to measure total flux we sum the product of \(\vec F\cdot\vec n\) times a small amount of surface area: \(\vec F\cdot \vec n\, dS\text{.}\)

A nice thing happens with the actual computation of flux: the \(\snorm{\vec r_u\times \vec r_v}\) terms go away. Consider:

\begin{align*}
\text{ Flux } \amp = \iint_\surfaceS \vec F\cdot \vec n\, dS\\
\amp = \iint_R \vec F\cdot \frac{\vec r_u\times \vec r_v}{\snorm{\vec r_u\times \vec r_v}}\snorm{\vec r_u\times \vec r_v}\, dA\\
\amp = \iint_R \vec F\cdot (\vec r_u\times \vec r_v)\, dA\text{.}
\end{align*}

The above only makes sense if \(\surfaceS\) is orientable; the normal vectors \(\vec n\) must vary continuously across \(\surfaceS\text{.}\) We assume that \(\vec n\) does vary continuously. (If the parametrization \(\vec r\) of \(\surfaceS\) is smooth, then our above definition of \(\vec n\) will vary continuously.)

Let \(\vec F\) be a vector field with continuous components defined on an orientable surface \(\surfaceS\) with normal vector \(\vec n\text{.}\) The flux of \(\vec F\) across \(\surfaceS\) is

\begin{equation*}
\text{ Flux } = \iint_\surfaceS \vec F\cdot \vec n\, dS\text{.}
\end{equation*}

If \(\surfaceS\) is parametrized by \(\vec r(u,v)\text{,}\) which is smooth on its domain \(R\text{,}\) then

\begin{equation*}
\text{ Flux } = \iint_R \vec F\big(\vec r(u,v)\big)\cdot (\vec r_u\times \vec r_v)\, dA\text{.}
\end{equation*}

Since \(\surfaceS\) is orientable, we adopt the convention of saying one passes from the “back” side of \(\surfaceS\) to the “front” side when moving across the surface parallel to the direction of \(\vec n\text{.}\) Also, when \(\surfaceS\) is closed, it is natural to speak of the regions of space “inside” and “outside” \(\surfaceS\text{.}\) We also adopt the convention that when \(\surfaceS\) is a closed surface, \(\vec n\) should point to the outside of \(\surfaceS\text{.}\) If \(\vec n = \vec r_u\times\vec r_v\) points inside \(\surfaceS\text{,}\) use \(\vec n = \vec r_v\times \vec r_u\) instead.

When the computation of flux is positive, it means that the field is moving from the back side of \(\surfaceS\) to the front side; when flux is negative, it means the field is moving opposite the direction of \(\vec n\text{,}\) and is moving from the front of \(\surfaceS\) to the back. When \(\surfaceS\) is not closed, there is not a “right” and “wrong” direction in which \(\vec n\) should point, but one should be mindful of its direction to make full sense of the flux computation.

We demonstrate the computation of flux, and its interpretation, in the following examples.

Let \(\surfaceS\) be the surface given in Example 16.6.3, where \(\surfaceS\) is parametrized by \(\vec r(u,v) = \langle u, v(2-2u),3-2u-v(2-2u)\rangle\) on \(0\leq u\leq 1\text{,}\) \(0\leq v\leq 1\text{,}\) and let \(\vec F = \langle 1, x,-y\rangle\text{,}\) as shown in Figure 16.6.8. Find the flux of \(\vec F\) across \(\surfaceS\text{.}\)

Using our work from the previous example, we have \(\vec n = \vec r_u\times\vec r_v = \langle 4-4u,2-2u,2-2u\rangle\text{.}\) We also need \(\vec F\big(\vec r(u,v)\big) = \langle 1, u, -v(2-2u)\rangle\text{.}\)

Thus the flux of \(\vec F\) across \(\surfaceS\) is:

\begin{align*}
\text{ Flux } \amp = \iint_\surfaceS \vec F\cdot \vec n\, dS\\
\amp = \iint_R \langle 1,u,-v(2-2u)\rangle\cdot\langle 4-4u,2-2u,2-2u\rangle\, dA\\
\amp = \int_0^1\int_0^1 \big(-4u^2v-2u^2+8uv-2u-4v+4\big)\, du\, dv\\
\amp = 5/3\text{.}
\end{align*}

To make full use of this numeric answer, we need to know the direction in which the field is passing across \(\surfaceS\text{.}\) The graph in Figure 16.6.8 helps, but we need a method that is not dependent on a graph.

Pick a point \((u,v)\) in the interior of \(R\) and consider \(\vec n(u,v)\text{.}\) For instance, choose \((1/2,1/2)\) and look at \(\vec n(1/2,1/2) = \langle 2,1,1\rangle/\sqrt{6}\text{.}\) This vector has positive \(x\text{,}\) \(y\) and \(z\) components. Generally speaking, one has *some* idea of what the surface \(\surfaceS\) looks like, as that surface is for some reason important. In our case, we know \(\surfaceS\) is a plane with \(z\)-intercept of \(z=3\text{.}\) Knowing \(\vec n\) and the flux measurement of positive \(5/3\text{,}\) we know that the field must be passing from “behind” \(\surfaceS\text{,}\) i.e., the side the origin is on, to the “front” of \(\surfaceS\text{.}\)

Let \(\surfaceS_1\) be the unit disk in the \(xy\)-plane, and let \(\surfaceS_2\) be the paraboloid \(z=1-x^2-y^2\text{,}\) for \(z\geq 0\text{,}\) as graphed in Figure 16.6.11. Note how these two surfaces each have the unit circle as a boundary.

Let \(\vec F_1 = \langle 0,0,1\rangle\) and \(\vec F_2 = \langle 0,0,z\rangle\text{.}\) Using normal vectors for each surface that point “upward,” i.e., with a positive \(z\)-component, find the flux of each field across each surface.

We begin by parametrizing each surface.

The boundary of the unit disk in the \(xy\)-plane is the unit circle, which can be described with \(\langle \cos u,\sin u,0\rangle\text{,}\) \(0\leq u\leq 2\pi\text{.}\) To obtain the interior of the circle as well, we can scale by \(v\text{,}\) giving

\begin{equation*}
\vec r_1(u,v) = \langle v\cos u,v\sin u, 0\rangle, 0\leq u\leq 2\pi 0\leq v\leq 1\text{.}
\end{equation*}

As the boundary of \(\surfaceS_2\) is also the unit circle, the \(x\) and \(y\) components of \(\vec r_2\) will be the same as those of \(\vec r_1\text{;}\) we just need a different \(z\) component. With \(z = 1-x^2-y^2\text{,}\) we have

\begin{equation*}
\vec r_2(u,v) = \langle v\cos u,v\sin u, 1-v^2\cos^2u-v^2\sin^2u\rangle = \langle v\cos u,v\sin u, 1-v^2\rangle\text{,}
\end{equation*}

where \(0\leq u\leq 2\pi\) and \(0\leq v\leq 1\text{.}\)

We now compute the normal vectors \(\vec n_1\) and \(\vec n_2\text{.}\)

For \(\vec n_1\text{:}\) \(\vec r_{1u}= \langle -v\sin u, v\cos u,0\rangle\text{,}\) \(\vec r_{1v} = \langle \cos u,\sin u,0\rangle\text{,}\) so

\begin{equation*}
\vec n_1 = \vec r_{1u}\times \vec r_{1v} = \langle 0,0,-v\rangle\text{.}
\end{equation*}

As this vector has a negative \(z\)-component, we instead use

\begin{equation*}
\vec n_1 = \vec r_{1v}\times \vec r_{1u} = \langle 0,0,v\rangle\text{.}
\end{equation*}

Similarly, \(\vec n_2\text{:}\) \(\vec r_{2u}= \langle -v\sin u, v\cos u,0\rangle\text{,}\) \(\vec r_{2v} = \langle \cos u,\sin u,-2v\rangle\text{,}\) so

\begin{equation*}
\vec n_2 = \vec r_{2u}\times \vec r_{2v} = \langle -2v^2\cos u,-2v^2\sin u,-v\rangle\text{.}
\end{equation*}

Again, this normal vector has a negative \(z\)-component so we use

\begin{equation*}
\vec n_2 = \vec r_{2v}\times \vec r_{2u} = \langle 2v^2\cos u,2v^2\sin u,v\rangle\text{.}
\end{equation*}

We are now set to compute flux. Over field \(\vec F_1=\langle 0,0,1\rangle\text{:}\)

\begin{align*}
\text{ Flux across } \surfaceS_1 \amp = \iint_{\surfaceS_1} \vec F_1\cdot \vec n_1\, dS\\
\amp = \iint_R\langle 0,0,1\rangle\cdot\langle 0,0,v\rangle\, dA\\
\amp = \int_0^1\int_0^{2\pi} (v)\, du\, dv\\
\amp = \pi\text{.}
\end{align*}

\begin{align*}
\text{ Flux across } \surfaceS_2 \amp = \iint_{\surfaceS_2} \vec F_1\cdot \vec n_2\, dS\\
\amp = \iint_R\langle 0,0,1\rangle\cdot\langle 2v^2\cos u,2v^2\sin u,v\rangle\, dA\\
\amp = \int_0^1\int_0^{2\pi} (v)\, du\, dv\\
\amp = \pi\text{.}
\end{align*}

These two results are equal and positive. Each are positive because both normal vectors are pointing in the positive \(z\)-directions, as does \(\vec F_1\text{.}\) As the field passes through each surface in the direction of their normal vectors, the flux is measured as positive.

We can also intuitively understand why the results are equal. Consider \(\vec F_1\) to represent the flow of air, and let each surface represent a filter. Since \(\vec F_1\) is constant, and moving “straight up,” it makes sense that all air passing through \(\surfaceS_1\) also passes through \(\surfaceS_2\text{,}\) and vice-versa.

If we treated the surfaces as creating one piecewise-smooth surface \(\surfaceS\text{,}\) we would find the total flux across \(\surfaceS\) by finding the flux across each piece, being sure that each normal vector pointed to the outside of the closed surface. Above, \(\vec n_1\) does not point outside the surface, though \(\vec n_2\) does. We would instead want to use \(-\vec n_1\) in our computation. We would then find that the flux across \(\surfaceS_1\) is \(-\pi\text{,}\) and hence the total flux across \(\surfaceS\) is \(-\pi + \pi = 0\text{.}\) (As \(0\) is a special number, we should wonder if this answer has special significance. It does, which is briefly discussed following this example and will be more fully developed in the next section.)

We now compute the flux across each surface with \(\vec F_2=\langle 0,0,z\rangle\text{:}\)

\begin{align*}
\text{ Flux across } \surfaceS_1 \amp = \iint_{\surfaceS_1} \vec F_2\cdot \vec n_1\, dS .\\
\end{align*}

Over \(\surfaceS_1\text{,}\) \(\vec F_2 = \vec F_2\big(\vec r_2(u,v)\big) = \langle 0,0,0\rangle\text{.}\) Therefore,

\begin{align*} \amp = \iint_R\langle 0,0,0\rangle\cdot\langle 0,0,v\rangle\, dA\\ \amp = \int_0^1\int_0^{2\pi} (0)\, du\, dv\\ \amp = 0\text{.} \end{align*}
\begin{align*}
\text{ Flux across } \surfaceS_2 \amp = \iint_{\surfaceS_2} \vec F_2\cdot \vec n_2\, dS .\\
\end{align*}

Over \(\surfaceS_2\text{,}\) \(\vec F_2 = \vec F_2\big(\vec r_2(u,v)\big) = \langle 0,0,1-v^2\rangle\text{.}\) Therefore,

\begin{align*} \amp = \iint_R\langle 0,0,1-v^2\rangle\cdot\langle 2v^2\cos u,2v^2\sin u,v\rangle\, dA\\ \amp = \int_0^1\int_0^{2\pi} (v^3-v)\, du\, dv\\ \amp = \pi/2\text{.} \end{align*}This time the measurements of flux differ. Over \(\surfaceS_1\text{,}\) the field \(\vec F_2\) is just \(\vec 0\text{,}\) hence there is no flux. Over \(\surfaceS_2\text{,}\) the flux is again positive as \(\vec F_2\) points in the positive \(z\) direction over \(\surfaceS_2\text{,}\) as does \(\vec n_2\text{.}\)

In the previous example, the surfaces \(\surfaceS_1\) and \(\surfaceS_2\) form a closed surface that is piecewise smooth. That the measurement of flux across each surface was the same for some fields (and not for others) is reminiscent of a result from Section 16.4, where we measured flux across curves. The quick answer to why the flux was the same when considering \(\vec F_1\) is that \(\divv \vec F_1 = 0\text{.}\) In the next section, we’ll see the second part of the Divergence Theorem, which will more fully explain this occurrence. We will also explore Stokes’ Theorem, the spatial analogue to Green’s Theorem.

The videos in Figure 16.6.12–16.6.13 present some additional examples involving surface integrals of vector fields. Note that computing flux a cross a cube requires us to consider all six faces. This is a case where the Divergence Theorem can greatly simplify matters. For integrals over spheres, there are often simplifications possible, especially for “radial” vector fields (those parallel to the position vector \(\vec{r}(x,y,z)=\la x,y,z\ra\)).

In the plane, flux is a measurement of how much of the vector field passes across a ; in space, flux is a measurement of how much of the vector field passes across a .

When computing flux, what does it mean when the result is a negative number?

When \(\surfaceS\) is a closed surface, we choose the normal vector so that it points to the of the surface.

If \(\surfaceS\) is a plane, and \(\vec F\) is always parallel to \(\surfaceS\text{,}\) then the flux of \(\vec F\) across \(\surfaceS\) will be .

In the following exercises, a surface \(\surfaceS\) that represents a thin sheet of material with density \(\delta\) is given. Find the mass of each thin sheet.

\(\surfaceS\) is the plane \(z=x+y\) on \(-2\leq x\leq 2\text{,}\) \(-3\leq y\leq 3\text{,}\) with \(\delta(x,y,z) = z+10\text{.}\)

\(\surfaceS\) is the unit sphere, with \(\delta(x,y,z) = x+y+z+10\text{.}\)

In the following exercises, a surface \(\surfaceS\) and a vector field \(\vec F\) are given. Compute the flux of \(\vec F\) across \(\surfaceS\text{.}\) (If \(\surfaceS\) is not a closed surface, choose \(\vec n\) so that it has a positive \(z\)-component, unless otherwise indicated.)

\(\surfaceS\) is the plane \(z = 3x+y\) on \(0\leq x\leq 1\text{,}\) \(1\leq y\leq 4\text{;}\) \(\vec F = \langle x^2,-z,2y\rangle\text{.}\)

\(\surfaceS\) is the plane \(z = 8-x-y\) over the triangle with vertices at \((0,0)\text{,}\) \((1,0)\) and \((1,5)\text{;}\) \(\vec F = \langle 3,1,2\rangle\text{.}\)

\(\surfaceS\) is the paraboloid \(z = x^2+y^2\) over the unit disk; \(\vec F = \langle 1,0,0\rangle\text{.}\)

\(\surfaceS\) is the unit sphere; \(\vec F = \langle y-z,z-x,x-y\rangle\text{.}\)

\(\surfaceS\) is the square in space with corners at \((0,0,0)\text{,}\) \((1,0,0)\text{,}\) \((1,0,1)\) and \((0,0,1)\) (choose \(\vec n\) such that it has a positive \(y\)-component); \(\vec F = \langle 0,-z,y\rangle\text{.}\)

\(\surfaceS\) is the disk in the \(yz\)-plane with radius 1, centered at \((0,1,1)\) (choose \(\vec n\) such that it has a positive \(x\)-component); \(\vec F = \langle y,z,x\rangle\text{.}\)

\(\surfaceS\) is the closed surface composed of \(\surfaceS_1\text{,}\) whose boundary is the ellipse in the \(xy\)-plane described by \(\frac{x^2}{25}+\frac{y^2}9 = 1\) and \(\surfaceS_2\text{,}\) part of the elliptical paraboloid \(f(x,y) = 1-\frac{x^2}{25}-\frac{y^2}9\) (see graph); \(\vec F = \langle 5,2,3\rangle\text{.}\)

\(\surfaceS\) is the closed surface composed of \(\surfaceS_1\text{,}\) part of the unit sphere and \(\surfaceS_2\text{,}\) part of the plane \(z=1/2\) (see graph); \(\vec F = \langle x,-y,z\rangle\text{.}\)