Math Core

Lesson 5.6 · Vector Calculus

Parametric surfaces and surface integrals

A curve is described by one parameter, r(t)\mathbf{r}(t). A surface needs two, r(u,v)\mathbf{r}(u, v). Once you can parametrize a surface, you can find its area, integrate a function over it (a surface integral), and measure how much of a vector field flows through it (a flux integral). These are the two-dimensional cousins of arc length and line integrals, and they set up Stokes' theorem and the divergence theorem.

Parametric surfaces

A parametric surface is the set of points

r(u,v)=⟨x(u,v), y(u,v), z(u,v)⟩,(u,v) in a region D.\mathbf{r}(u, v) = \langle x(u, v),\ y(u, v),\ z(u, v)\rangle, \qquad (u, v) \text{ in a region } D.

Some standard parametrizations:

surfaceparametrizationparameter domain
graph z=g(x,y)z = g(x, y)⟨x, y, g(x,y)⟩\langle x,\ y,\ g(x, y)\ranglethe region below the surface
cylinder x2+y2=a2x^2 + y^2 = a^2⟨acos⁡θ, asin⁡θ, z⟩\langle a\cos\theta,\ a\sin\theta,\ z\rangle0≤θ≤2π0 \le \theta \le 2\pi, zz in some interval
sphere x2+y2+z2=a2x^2 + y^2 + z^2 = a^2⟨asin⁡ϕcos⁡θ, asin⁡ϕsin⁡θ, acos⁡ϕ⟩\langle a\sin\phi\cos\theta,\ a\sin\phi\sin\theta,\ a\cos\phi\rangle0≤ϕ≤π0 \le \phi \le \pi, 0≤θ≤2π0 \le \theta \le 2\pi
plane through P0P_0 containing a\mathbf{a}, b\mathbf{b}P0+ua+vbP_0 + u\mathbf{a} + v\mathbf{b}all u,vu, v

The sphere uses the spherical angles: ϕ\phi is measured down from the positive zz-axis and θ\theta is the usual polar angle.

Tangent vectors, normals and surface area

Holding vv fixed and varying uu traces a curve on the surface with tangent vector ru=∂r/∂u\mathbf{r}_u = \partial\mathbf{r}/\partial u; similarly rv\mathbf{r}_v. These two tangent vectors span the tangent plane, so their cross product ru×rv\mathbf{r}_u\times\mathbf{r}_v is a normal vector to the surface.

A small rectangle Δu\Delta u by Δv\Delta v in the parameter domain maps to a small patch on the surface that is nearly a parallelogram with sides ru Δu\mathbf{r}_u\,\Delta u and rv Δv\mathbf{r}_v\,\Delta v. Its area is ∣ru×rv∣ Δu Δv\lvert\mathbf{r}_u\times\mathbf{r}_v\rvert\,\Delta u\,\Delta v.

Surface area and the surface element

The area element on the surface is dS=∣ru×rv∣ dAdS = \lvert\mathbf{r}_u\times\mathbf{r}_v\rvert\,dA, and

Area(S)=∬D∣ru×rv∣ dA.\text{Area}(S) = \iint_D \lvert\mathbf{r}_u\times\mathbf{r}_v\rvert\,dA.

For a graph z=g(x,y)z = g(x, y): rx×ry=⟨−gx, −gy, 1⟩\mathbf{r}_x\times\mathbf{r}_y = \langle -g_x,\ -g_y,\ 1\rangle, so dS=1+gx2+gy2 dAdS = \sqrt{1 + g_x^2 + g_y^2}\,dA.

For a sphere of radius aa: dS=a2sin⁡ϕ dϕ dθdS = a^2\sin\phi\,d\phi\,d\theta.

Worked example: A cylinder

Find the area of the cylinder x2+y2=4x^2 + y^2 = 4 for 0≤z≤30 \le z \le 3.

With r(θ,z)=⟨2cos⁡θ, 2sin⁡θ, z⟩\mathbf{r}(\theta, z) = \langle 2\cos\theta,\ 2\sin\theta,\ z\rangle: rθ=⟨−2sin⁡θ, 2cos⁡θ, 0⟩\mathbf{r}_\theta = \langle -2\sin\theta,\ 2\cos\theta,\ 0\rangle and rz=⟨0,0,1⟩\mathbf{r}_z = \langle 0, 0, 1\rangle. Then rθ×rz=⟨2cos⁡θ, 2sin⁡θ, 0⟩\mathbf{r}_\theta\times\mathbf{r}_z = \langle 2\cos\theta,\ 2\sin\theta,\ 0\rangle, which has length 22. So

Area=∫03 ⁣ ⁣∫02π2 dθ dz=12π,\text{Area} = \int_0^3\!\!\int_0^{2\pi} 2\,d\theta\,dz = 12\pi,

which matches circumference times height, 4π⋅34\pi\cdot 3.

Worked example: A paraboloid

Find the area of the part of z=x2+y2z = x^2 + y^2 that lies below the plane z=2z = 2.

Here gx=2xg_x = 2x, gy=2yg_y = 2y, so dS=1+4x2+4y2 dAdS = \sqrt{1 + 4x^2 + 4y^2}\,dA. The region below is the disk x2+y2≤2x^2 + y^2 \le 2. In polar coordinates, with u=1+4r2u = 1 + 4r^2:

Area=∫02π ⁣ ⁣∫021+4r2 r dr dθ=2π⋅112[(1+4r2)3/2]02=π6(27−1)=13π3.\text{Area} = \int_0^{2\pi}\!\!\int_0^{\sqrt{2}} \sqrt{1 + 4r^2}\,r\,dr\,d\theta = 2\pi\cdot\frac{1}{12}\Big[(1 + 4r^2)^{3/2}\Big]_0^{\sqrt{2}} = \frac{\pi}{6}(27 - 1) = \frac{13\pi}{3}.

Surface integrals of functions

Just as ∫Cf ds\int_C f\,ds weights a curve by ff, a surface integral weights a surface by ff:

∬Sf(x,y,z) dS=∬Df(r(u,v)) ∣ru×rv∣ dA.\iint_S f(x, y, z)\,dS = \iint_D f\big(\mathbf{r}(u, v)\big)\,\lvert\mathbf{r}_u\times\mathbf{r}_v\rvert\,dA.

If ff is a density (mass per unit area) of a thin sheet, this is the sheet's mass. With f=1f = 1 it is the surface area.

Worked example: Integrating over a hemisphere

Evaluate ∬Sz dS\displaystyle\iint_S z\,dS, where SS is the upper half of the unit sphere.

On the unit sphere z=cos⁡ϕz = \cos\phi and dS=sin⁡ϕ dϕ dθdS = \sin\phi\,d\phi\,d\theta, with 0≤ϕ≤π20 \le \phi \le \tfrac{\pi}{2}:

∬Sz dS=∫02π ⁣ ⁣∫0π/2cos⁡ϕsin⁡ϕ dϕ dθ=2π⋅12=π.\iint_S z\,dS = \int_0^{2\pi}\!\!\int_0^{\pi/2}\cos\phi\sin\phi\,d\phi\,d\theta = 2\pi\cdot\frac{1}{2} = \pi.

Oriented surfaces and flux

To measure flow through a surface you need to say which side is which. A surface is oriented when you choose a continuous unit normal n\mathbf{n} at every point. For a closed surface such as a sphere, the standard choice is the outward normal. For a graph z=g(x,y)z = g(x, y), "upward" means the normal ⟨−gx,−gy,1⟩\langle -g_x, -g_y, 1\rangle, whose zz-component is positive. (A Möbius strip has only one side and can't be oriented.)

Definition

Flux integral

The flux of F\mathbf{F} across the oriented surface SS is

∬SF⋅dS=∬SF⋅n dS=∬DF(r(u,v))⋅(ru×rv) dA,\iint_S \mathbf{F}\cdot d\mathbf{S} = \iint_S \mathbf{F}\cdot\mathbf{n}\,dS = \iint_D \mathbf{F}\big(\mathbf{r}(u, v)\big)\cdot(\mathbf{r}_u\times\mathbf{r}_v)\,dA,

where the parametrization is chosen so that ru×rv\mathbf{r}_u\times\mathbf{r}_v points in the direction of n\mathbf{n}.

For a graph z=g(x,y)z = g(x, y) with upward orientation and F=⟨P,Q,R⟩\mathbf{F} = \langle P, Q, R\rangle:

∬SF⋅dS=∬D(−P gx−Q gy+R)dA.\iint_S \mathbf{F}\cdot d\mathbf{S} = \iint_D \left(-P\,g_x - Q\,g_y + R\right) dA.

If F\mathbf{F} is the velocity of a fluid, the flux is the volume of fluid crossing SS per unit time, in the direction of n\mathbf{n}. Notice that the ∣ru×rv∣\lvert\mathbf{r}_u\times\mathbf{r}_v\rvert factors cancel: n dS=(ru×rv) dA\mathbf{n}\,dS = (\mathbf{r}_u\times\mathbf{r}_v)\,dA.

Worked example: Flux through a paraboloid

Find the flux of F=⟨y, x, z⟩\mathbf{F} = \langle y,\ x,\ z\rangle upward through the part of z=4−x2−y2z = 4 - x^2 - y^2 above the xyxy-plane.

Here gx=−2xg_x = -2x and gy=−2yg_y = -2y, so

−P gx−Q gy+R=2xy+2xy+(4−x2−y2)=4xy+4−r2.-P\,g_x - Q\,g_y + R = 2xy + 2xy + (4 - x^2 - y^2) = 4xy + 4 - r^2.

Over the disk r≤2r \le 2, the term 4xy4xy integrates to 00 by symmetry. The rest is

∫02π ⁣ ⁣∫02(4−r2) r dr dθ=2π(8−4)=8π.\int_0^{2\pi}\!\!\int_0^2 (4 - r^2)\,r\,dr\,d\theta = 2\pi\left(8 - 4\right) = 8\pi.

Common mistake

Scalar surface integrals use ∣ru×rv∣\lvert\mathbf{r}_u\times\mathbf{r}_v\rvert (a length); flux integrals use ru×rv\mathbf{r}_u\times\mathbf{r}_v itself (a vector, dotted with F\mathbf{F}). And for flux, check that your normal points the way the problem asks. If it points the wrong way, the answer changes sign.

Tip

On a sphere of radius aa centered at the origin, the outward unit normal is simply n=⟨x,y,z⟩a\mathbf{n} = \dfrac{\langle x, y, z\rangle}{a}. So for F=⟨x,y,z⟩\mathbf{F} = \langle x, y, z\rangle, F⋅n=a\mathbf{F}\cdot\mathbf{n} = a everywhere and the flux is a⋅4πa2=4πa3a\cdot 4\pi a^2 = 4\pi a^3.

Practice

Practice 1

For the plane r(u,v)=⟨u, v, 3−u−2v⟩\mathbf{r}(u, v) = \langle u,\ v,\ 3 - u - 2v\rangle, compute the normal vector ru×rv\mathbf{r}_u\times\mathbf{r}_v.

Enter a point like (2, -3)

Practice 2

Find the area of the part of the plane z=1+2x+2yz = 1 + 2x + 2y that lies above the rectangle 0≤x≤10 \le x \le 1, 0≤y≤20 \le y \le 2.

Enter a number. Fractions like 3/4 and sqrt(2) are OK.

Practice 3

Use dS=a2sin⁡ϕ dϕ dθdS = a^2\sin\phi\,d\phi\,d\theta to find the surface area of a sphere of radius 33.

Enter a number. Fractions like 3/4 and sqrt(2) are OK.

Practice 4

Evaluate ∬Sz dS\displaystyle\iint_S z\,dS, where SS is the cylinder x2+y2=1x^2 + y^2 = 1, 0≤z≤20 \le z \le 2 (the side only, no top or bottom).

Enter a number. Fractions like 3/4 and sqrt(2) are OK.

Practice 5

Find the flux of the constant field F=⟨0,0,3⟩\mathbf{F} = \langle 0, 0, 3\rangle upward through the disk x2+y2≤4x^2 + y^2 \le 4 in the plane z=1z = 1.

Enter a number. Fractions like 3/4 and sqrt(2) are OK.

Practice 6

Find the outward flux of F=⟨x,y,z⟩\mathbf{F} = \langle x, y, z\rangle across the sphere x2+y2+z2=4x^2 + y^2 + z^2 = 4.

Enter a number. Fractions like 3/4 and sqrt(2) are OK.

Practice 7

Find the flux of F=⟨x,y,0⟩\mathbf{F} = \langle x, y, 0\rangle upward through the part of the paraboloid z=1−x2−y2z = 1 - x^2 - y^2 above the xyxy-plane.

Enter a number. Fractions like 3/4 and sqrt(2) are OK.

Practice 8

Which quantity does ∬SF⋅dS\displaystyle\iint_S \mathbf{F}\cdot d\mathbf{S} change sign under, while ∬Sf dS\displaystyle\iint_S f\,dS does not?