Surface Integral of a Scalar-Valued Function

Learning Objectives

Use a surface integral to calculate the area of a given surface.

Now that we can parameterize surfaces and we can calculate their surface areas, we are able to define surface integrals. First, let’s look at the surface integral of a scalar-valued function. Informally, the surface integral of a scalar-valued function is an analog of a scalar line integral in one higher dimension. The domain of integration of a scalar line integral is a parameterized curve (a one-dimensional object); the domain of integration of a scalar surface integral is a parameterized surface (a two-dimensional object). Therefore, the definition of a surface integral follows the definition of a line integral quite closely. For scalar line integrals, we chopped the domain curve into tiny pieces, chose a point in each piece, computed the function at that point, and took a limit of the corresponding Riemann sum. For scalar surface integrals, we chop the domain region (no longer a curve) into tiny pieces and proceed in the same fashion.

Let [latex]S[/latex] be a piecewise smooth surface with parameterization [latex]{\bf{r}}(u,v)=\langle{x}(u,v),y(u,v)z(u,v)\rangle[/latex] with parameter domain [latex]D[/latex] and let [latex]f(x, y, z)[/latex] be a function with a domain that contains [latex]S[/latex]. For now, assume the parameter domain [latex]D[/latex] is a rectangle, but we can extend the basic logic of how we proceed to any parameter domain (the choice of a rectangle is simply to make the notation more manageable). Divide rectangle [latex]D[/latex] into subrectangles [latex]D_{ij}[/latex] with horizontal width [latex]\Delta{u}[/latex] and vertical length [latex]\Delta{v}[/latex]. Suppose that [latex]i[/latex] ranges from [latex]1[/latex] to [latex]m[/latex] and [latex]j[/latex] ranges from [latex]1[/latex] to [latex]n[/latex] so that [latex]D[/latex] is subdivided into [latex]mn[/latex] rectangles. This division of [latex]D[/latex] into subrectangles gives a corresponding division of [latex]S[/latex] into pieces [latex]S_{ij}[/latex]. Choose point [latex]P_{ij}[/latex] in each piece [latex]S_{ij}[/latex], evaluate [latex]P_{ij}[/latex] at [latex]f[/latex], and multiply by area [latex]S_{ij}[/latex] to form the Riemann sum

[latex]\large{\displaystyle\sum_{i=1}^m\displaystyle\sum_{j=1}^nf(P_{ij})\Delta{S_{ij}}}[/latex].

To define a surface integral of a scalar-valued function, we let the areas of the pieces of [latex]S[/latex] shrink to zero by taking a limit.

definition

The surface integral of a scalar-valued function of [latex]f[/latex] over a piecewise smooth surface [latex]S[/latex] is

[latex]\large{\displaystyle\iint_Sf(x,y,z)dS=\displaystyle\lim_{m,n\to\infty}\displaystyle\sum_{i=1}^m\displaystyle\sum_{j=1}^nf(P_{ij})\Delta{S_{ij}}}[/latex].

Again, notice the similarities between this definition and the definition of a scalar line integral. In the definition of a line integral we chop a curve into pieces, evaluate a function at a point in each piece, and let the length of the pieces shrink to zero by taking the limit of the corresponding Riemann sum. In the definition of a surface integral, we chop a surface into pieces, evaluate a function at a point in each piece, and let the area of the pieces shrink to zero by taking the limit of the corresponding Riemann sum. Thus, a surface integral is similar to a line integral but in one higher dimension.

The definition of a scalar line integral can be extended to parameter domains that are not rectangles by using the same logic used earlier. The basic idea is to chop the parameter domain into small pieces, choose a sample point in each piece, and so on. The exact shape of each piece in the sample domain becomes irrelevant as the areas of the pieces shrink to zero.

Scalar surface integrals are difficult to compute from the definition, just as scalar line integrals are. To develop a method that makes surface integrals easier to compute, we approximate surface areas [latex]\Delta{S_{ij}}[/latex] with small pieces of a tangent plane, just as we did in the previous subsection. Recall the definition of vectors [latex]{\bf{t}}_u[/latex] and [latex]{\bf{t}}_v[/latex]:

[latex]\large{{\bf{t}}_u=\left\langle\frac{\partial{x}}{\partial{u}},\frac{\partial{y}}{\partial{u}},\frac{\partial{z}}{\partial{u}}\right\rangle\text{ and }{\bf{t}}_v=\left\langle\frac{\partial{x}}{\partial{v}},\frac{\partial{y}}{\partial{v}},\frac{\partial{z}}{\partial{v}}\right\rangle}[/latex].

From the material we have already studied, we know that

[latex]\large{\Delta{S}_{ij}\approx||{\bf{t}}_u(P_{ij})\times{\bf{t}}_v(P_{ij})||\Delta{u}\Delta{v}}[/latex].

Therefore,

[latex]\large{\displaystyle\iint_Sf(x,y,z)dS\approx\displaystyle\lim_{m,n\to\infty}\displaystyle\sum_{i=1}^m\displaystyle\sum_{j=1}^nf(P_{ij})||{\bf{t}}_u(P_{ij})\times{\bf{t}}_v(P_{ij})||\Delta{u}\Delta{v}}[/latex].

This approximation becomes arbitrarily close to [latex]\displaystyle\lim_{m,n\to\infty}\displaystyle\sum_{i=1}^m\displaystyle\sum_{j=1}^nf(P_{ij})\Delta{S}_{ij}[/latex] as we increase the number of pieces [latex]S_{ij}[/latex] by letting [latex]m[/latex] and [latex]n[/latex] go to infinity. Therefore, we have the following equation to calculate scalar surface integrals:

[latex]\large{\displaystyle\iint_Sf(x,y,z)dS=\displaystyle\iint_Df({\bf{r}}(u,v)||{\bf{t}}_u\times{\bf{t}}_v||}dA[/latex].

The surface integral of a scalar-valued function allows us to calculate a surface integral by transforming it into a double integral. This equation for surface integrals is analogous to the surface integral equation for line integrals:

[latex]\large{\displaystyle\iint_Cf(x,y,z)ds=\displaystyle\int_a^bf({\bf{r}}(t))||{\bf{r}}^\prime(t)||dt}[/latex].

In this case, vector [latex]{\bf{t}}_u\times{\bf{t}}_v[/latex] is perpendicular to the surface, whereas vector [latex]{\bf{r}}^\prime[/latex] is tangent to the curve.

Example: Calculating a surface integral

Calculate surface integral [latex]\displaystyle\iint_S5dS[/latex], where [latex]S[/latex] is the surface with parameterization [latex]{\bf{r}}(u,v)=\langle{v},u^2,v\rangle[/latex] for [latex]0\leq{u}\leq2[/latex] and [latex]0\leq{v}\leq{u}[/latex].

Show Solution

Example: calculating the surface integral of a cylinder

Calculate surface integral [latex]\displaystyle\iint_S(x+y^2)dS[/latex], where [latex]S[/latex] is cylinder [latex]x^2+y^2=4[/latex], [latex]0\leq{z}\leq3[/latex] (Figure 1).

<img src="/apps/archive/20220422.171947/resources/2de01bbf14da890ca1a4c3a6c3b1a93dca62f7da" data-media-type="image/jpeg" alt="A graph in three dimensions of a cylinder. The base of the cylinder is on the (x,z) plane, with center on the y axis. It stretches along the y axis." id="39">

Figure 1. Integrating function [latex]f(x,y,z)=x+y^2[/latex] over a cylinder.

Show Solution

try it

Calculate [latex]\displaystyle\iint_S(x^2-z)dS[/latex], where [latex]S[/latex] is the surface with parameterization [latex]{\bf{r}}(u,v)=\langle{u},u^2+v^2,1\rangle[/latex], [latex]0\leq{u}\leq2, \ 0\leq{v}\leq3[/latex].

Show Solution

Example: calculating the surface integral of a piece of a sphere

Calculate surface integral [latex]\displaystyle\iint_Sf(x,y,z)dS[/latex], where [latex]f(x, y, z)=z^{2}[/latex] and [latex]S[/latex] is the surface that consists of the piece of sphere [latex]x^{2}+y^{2}+z^{2}=4[/latex] that lies on or above plane [latex]z=1[/latex] and the disk that is enclosed by intersection plane [latex]z=1[/latex] and the given sphere (Figure 2).

<img src="/apps/archive/20220422.171947/resources/8060a589ddea4ad92316c95589f1b486dffea464" data-media-type="image/jpeg" alt="A diagram in three dimensions of the upper half of a sphere. The center is at the origin, and the radius is 2. The top part above the plane z=1 is cut off and shaded; the rest is simply an outline of the hemisphere. The top section has center at (0,0,1) and radius of radical three." id="42">

Figure 2. Calculating a surface integral over surface [latex]S[/latex].

Show Solution

Solution

Notice that [latex]S[/latex] is not smooth but is piecewise smooth; [latex]S[/latex] can be written as the union of its base [latex]S_1[/latex] and its spherical top [latex]S_2[/latex], and both [latex]S_1[/latex] and [latex]S_2[/latex] are smooth. Therefore, to calculate [latex]\displaystyle\iint_Sz^2dS[/latex], we write this integral as [latex]\displaystyle\iint_Sz^2d{S_1}+\displaystyle\iint_{S_2}z^2dS[/latex] and we calculate integrals [latex]\displaystyle\iint_{S_1}z^2dS[/latex] and [latex]\displaystyle\iint_{S_2}z^2dS[/latex].

First, we calculate [latex]\displaystyle\iint_{S_1}z^2dS[/latex]. To calculate this integral we need a parameterization of [latex]S_1[/latex]. This surface is a disk in plane [latex]z=1[/latex] centered at [latex](0, 0, 1)[/latex]. To parameterize this disk, we need to know its radius. Since the disk is formed where plane [latex]z=1[/latex] intersects sphere [latex]x^{2}+y^{2}+z^{2}=4[/latex], we can substitute [latex]z=1[/latex] into equation [latex]x^{2}+y^{2}+z^{2}=4[/latex]:

[latex]x^2+y^2+1=4\Rightarrow{x}^2+y^2=3[/latex].

Therefore, the radius of the disk is [latex]\sqrt3[/latex] and a parameterization of [latex]S_1[/latex] is [latex]{\bf{r}}(u,v)=\langle{u}\cos{v},u\sin{v},1\rangle[/latex], [latex]0\leq{u}\leq\sqrt3[/latex],[latex]0\leq{v}\leq2\pi[/latex]. The tangent vectors are [latex]{\bf{t}}_u=\langle\cos{v},\sin{v},0\rangle[/latex] and [latex]{\bf{t}}_v=\langle-u\sin{v},u\cos{v},0\rangle[/latex], and thus

[latex]{\bf{t}}_u\times{\bf{t}}_v=\begin{vmatrix}{\bf{i}}&{\bf{j}}&{\bf{k}} \\ \cos{v}&\sin{v}&0 \\-u\sin{v}&u\cos{v}&0\end{vmatrix}=\langle0,0,u\cos^2v+u\sin^2v\rangle=\langle0,0,u\rangle[/latex].

The magnitude of this vector is [latex]u[/latex]. Therefore,

[latex]\begin{aligned} \displaystyle\iint_{S_1}z^2dS&=\displaystyle\int_0^{\sqrt3}\displaystyle\int_0^{2\pi}f({\bf{r}}(u,v))||{\bf{t}}_u\times{\bf{t}}_v|| \ dvdu \\ &=\displaystyle\int_0^{\sqrt3}\displaystyle\int_0^{2\pi}u \ dvdu \\ &=2\pi\displaystyle\int_0^{\sqrt3} \ udu \\ &=3\pi \end{aligned}[/latex].

Now we calculate [latex]\displaystyle\iint_{S_2} \ dS[/latex]. To calculate this integral, we need a parameterization of [latex]S_2[/latex]. The parameterization of full sphere [latex]x^{2}+y^{2}+z^{2}=4[/latex] is

[latex]{\bf{r}}(\phi,\theta)=\langle2\cos\theta\sin\phi,2\sin\theta\sin\phi,2\cos\phi\rangle, \ 0\leq\theta\leq2\pi, \ 0\leq\phi\leq\pi[/latex].

Since we are only taking the piece of the sphere on or above plane [latex]z=1[/latex], we have to restrict the domain of [latex]\phi[/latex]. To see how far this angle sweeps, notice that the angle can be located in a right triangle, as shown in Figure 3 (the [latex]\sqrt3[/latex] comes from the fact that the base of [latex]S[/latex] is a disk with radius [latex]\sqrt3[/latex]). Therefore, the tangent of [latex]\phi[/latex] is [latex]\sqrt3[/latex], which implies that [latex]\phi[/latex] is [latex]\pi/6[/latex]. We now have a parameterization of [latex]S_2[/latex]:

[latex]{\bf{r}}(\phi,\theta)=\langle2\cos\theta\sin\phi,2\sin\theta\sin\phi,2\cos\phi\rangle, \ 0\leq\theta\leq2\pi, \ 0\leq\phi\leq\pi/3[/latex].

<img src="/apps/archive/20220422.171947/resources/476497f01653538dbc08bef585c08834b0c1278f" data-media-type="image/jpeg" alt="A diagram of a plane within the three-dimensional coordinate system. Two points are marked on the z axis: (0,0,2) and (0,0,1). The distance from the origin to (0,0,1) is marked as 1, the horizontal distance between the point (0,0,1) and a point of the sphere is labeled radical three, and the angle between the origin and the point on the sphere is theta. There is a line drawn from the origin to the point on the sphere, and this forms a triangle." id="43">

Figure 3. The maximum value of [latex]\phi[/latex] has a tangent value of [latex]{\sqrt3}[/latex]

The tangent vectors are

[latex]{\bf{t}}_\phi=\langle2\cos\theta\cos\phi,2\sin\theta\cos\phi,-2\sin\phi\rangle\text{ and }{\bf{t}}_\theta=\langle-2\sin\theta\sin\phi,u\cos\theta\sin\phi,0\rangle[/latex],

and thus

[latex]\begin{aligned} {\bf{t}}_\phi\times{\bf{t}}_\theta&=\begin{vmatrix} {\bf{i}}&{\bf{j}}&{\bf{k}} \\ 2\cos\theta\cos\phi&2\sin\theta\cos\phi&-2\sin\phi \\ -2\sin\theta\sin\phi&2\cos\theta\sin\phi&0 \end{vmatrix} \\ &=\langle4\cos\theta\sin^2\phi,4\sin\theta\sin^2\phi,4\cos^2\theta\sin\phi+4\sin^2\theta\cos\phi\sin\phi\rangle \\ &=\langle4\cos\theta\sin^2\phi,4\sin\theta\sin^2\phi,4\cos\phi\sin\phi\rangle \end{aligned}[/latex].

The magnitude of this vector is

[latex]\begin{aligned} ||{\bf{t}}_\phi\times{\bf{t}}_\theta||&=\sqrt{16\cos^2\theta\sin^4\phi+16\sin^2\theta\sin^4\phi+16\cos^2\phi\sin^2\phi} \\ &=4\sqrt{\sin^4\phi+\cos^2\phi\sin^2\phi} \end{aligned}[/latex].

Therefore,

[latex]\begin{aligned} \displaystyle\iint_{S_2}z^2dS&=\displaystyle\int_0^{\pi/3}\displaystyle\int_0^{2\pi}f({\bf{r}}(\phi,\theta))||{\bf{t}}_\phi\times{\bf{t}}_\theta|| \ d\theta{d}\phi \\ &=\displaystyle\int_0^{\pi/3}\displaystyle\int_0^{2\pi}16\cos^2\phi\sqrt{\sin^4\phi+\cos^2\phi\sin^2\phi}d\theta{d}\phi \\ &=32\pi\displaystyle\int_0^{\pi/3}\cos^2\phi\sqrt{\sin^4\phi+\cos^2\phi\sin^2\phi}d\phi \\ &=32\pi\displaystyle\int_0^{\pi/3}\cos^2\phi\sin\phi\sqrt{\sin^2\phi+\cos^2\phi}d\phi \\ &=32\pi\displaystyle\int_0^{\pi/3}\cos^2\phi\sin\phi{d}\phi \\ &=32\pi\left[-\frac{\cos^3\phi}3\right]_0^{\pi/6}=32\pi\left[\frac13-\frac{\sqrt3}8\right]=\frac{28\pi}3 \end{aligned}[/latex].

Since

[latex]\displaystyle\iint_Sz^2 \ dS=\displaystyle\iint_{S_1}z^2 \ dS+\displaystyle\iint_{S_2}z^2 \ dS=3\pi\frac{28\pi}3=\frac{37\pi}3[/latex].

Analysis

In this example we broke a surface integral over a piecewise surface into the addition of surface integrals over smooth subsurfaces. There were only two smooth subsurfaces in this example, but this technique extends to finitely many smooth subsurfaces.

try it

Calculate line integral [latex]\displaystyle\iint_S(x-y)dS[/latex], where [latex]S[/latex] is cylinder [latex]x^2+y^2=1[/latex],[latex]0\leq{z}\leq2[/latex], including the circular top and bottom.

Show Solution

Watch the following video to see the worked solution to the above Try It

You can view the transcript for “CP 6.55” here (opens in new window).Scalar surface integrals have several real-world applications. Recall that scalar line integrals can be used to compute the mass of a wire given its density function. In a similar fashion, we can use scalar surface integrals to compute the mass of a sheet given its density function. If a thin sheet of metal has the shape of surface [latex]S[/latex] and the density of the sheet at point [latex](x, y, z)[/latex] is [latex]\rho(x,y,z)[/latex], then mass [latex]m[/latex] of the sheet is [latex]m=\displaystyle\iint_S\rho(x,y,z)dS[/latex].

Example: calculating the mass of a sheet

A flat sheet of metal has the shape of surface [latex]z=1+x+2y[/latex] that lies above rectangle [latex]0\leq{x}\leq4[/latex] and [latex]0\leq{y}\leq2[/latex]. If the density of the sheet is given by [latex]\rho(x,y,z)=x^2yz[/latex], what is the mass of the sheet?

Show Solution

try it

A piece of metal has a shape that is modeled by paraboloid [latex]z=x^2+y^2[/latex], [latex]0\leq{z}\leq4[/latex], and the density of the metal is given by [latex]\rho(x,y,z)=z+1[/latex]. Find the mass of the piece of metal.

Show Solution

Module 6: Vector Calculus