10 Parametric Surfaces and Surface Integrals

A progressive guide to parametrizing surfaces, finding tangent planes and normals, computing surface area and scalar surface integrals, and evaluating flux through oriented surfaces.

Parametrizing a Surface

A surface in three-dimensional space is two-dimensional, so its points generally require two parameters. A parametrized surface is written as

r(u,v)=⟨x(u,v),y(u,v),z(u,v)⟩,\mathbf r(u,v)=\langle x(u,v),y(u,v),z(u,v)\rangle,

where (u,v)(u,v) ranges over a region DD in the parameter plane. The region DD is the parameter domain, and each point of DD maps to a point on the surface.

Holding one parameter constant produces coordinate curves. Setting uu equal to a constant lets vv vary, while setting vv equal to a constant lets uu vary. These curves describe how the parameter domain traces the surface.

Common choices include:

  • For a graph z=f(x,y)z=f(x,y), use r(x,y)=⟨x,y,f(x,y)⟩\mathbf r(x,y)=\langle x,y,f(x,y)\rangle.

  • For the lateral surface of a cylinder of radius aa, use r(θ,z)=⟨acos⁡θ,asin⁡θ,z⟩\mathbf r(\theta,z)=\langle a\cos\theta,a\sin\theta,z\rangle.

  • For a sphere of radius aa, use r(ϕ,θ)=⟨asin⁡ϕcos⁡θ,asin⁡ϕsin⁡θ,acos⁡ϕ⟩\mathbf r(\phi,\theta)=\langle a\sin\phi\cos\theta,a\sin\phi\sin\theta,a\cos\phi\rangle, with 0≤ϕ≤π0\le\phi\le\pi and 0≤θ≤2π0\le\theta\le 2\pi.

A is one for which the two parameter directions are not parallel, so they determine a genuine tangent plane.

Takeaway: Begin every surface calculation by identifying a parametrization and its parameter domain.

Tangent Planes and Normals

The tangent directions of r(u,v)\mathbf r(u,v) are the fields

ru=∂r∂u,rv=∂r∂v.\mathbf r_u=\frac{\partial\mathbf r}{\partial u},\qquad \mathbf r_v=\frac{\partial\mathbf r}{\partial v}.

Their cross product gives a :

n=ru×rv.\mathbf n=\mathbf r_u\times\mathbf r_v.

A parametrization is regular at a point when

ru×rv≠0.\mathbf r_u\times\mathbf r_v\ne\mathbf 0.

At P=r(u0,v0)=(x0,y0,z0)P=\mathbf r(u_0,v_0)=(x_0,y_0,z_0), the tangent plane is

(ru×rv)∣(u0,v0)⋅⟨x−x0,y−y0,z−z0⟩=0.\left.(\mathbf r_u\times\mathbf r_v)\right|_{(u_0,v_0)}\cdot\langle x-x_0,y-y_0,z-z_0\rangle=0.

The two unit normal choices are

N=±ru×rv∥ru×rv∥.\mathbf N=\pm\frac{\mathbf r_u\times\mathbf r_v}{\|\mathbf r_u\times\mathbf r_v\|}.

For a graph z=f(x,y)z=f(x,y),

rx×ry=⟨−fx,−fy,1⟩.\mathbf r_x\times\mathbf r_y=\langle-f_x,-f_y,1\rangle.

This vector points upward because its third component is positive. The downward normal is the negative of the upward normal.

Takeaway: Use the cross product to obtain both the tangent-plane normal and the oriented area direction.

Surface Area

A small parameter rectangle with side lengths Δu\Delta u and Δv\Delta v is approximated on the surface by a parallelogram spanned by ruΔu\mathbf r_u\Delta u and rvΔv\mathbf r_v\Delta v. Its approximate area is

∥ru×rv∥ Δu Δv.\|\mathbf r_u\times\mathbf r_v\|\,\Delta u\,\Delta v.

Therefore, the is

dS=∥ru×rv∥ du dv,dS=\|\mathbf r_u\times\mathbf r_v\|\,du\,dv,

and the total area is

Area⁡(S)=∬D∥ru×rv∥ dA.\operatorname{Area}(S)=\iint_D\|\mathbf r_u\times\mathbf r_v\|\,dA.

For a graph z=f(x,y)z=f(x,y), this becomes

Area⁡(S)=∬D1+fx2+fy2 dA.\operatorname{Area}(S)=\iint_D\sqrt{1+f_x^2+f_y^2}\,dA.

For example, consider z=x+2yz=x+2y above 0≤x≤10\le x\le1 and 0≤y≤20\le y\le2. Since fx=1f_x=1 and fy=2f_y=2,

Area⁡(S)=∫01∫026 dy dx=26.\operatorname{Area}(S)=\int_0^1\int_0^2\sqrt{6}\,dy\,dx=2\sqrt{6}.

Takeaway: Surface area uses the magnitude ∥ru×rv∥\|\mathbf r_u\times\mathbf r_v\|, not the oriented cross product itself.

Scalar Surface Integrals

A has the form ∬Sg dS\iint_S g\,dS, where g(x,y,z)g(x,y,z) is a scalar function. Under the parametrization r(u,v)\mathbf r(u,v), compute it by

∬Sg dS=∬Dg(r(u,v))∥ru×rv∥ du dv.\iint_S g\,dS=\iint_D g\bigl(\mathbf r(u,v)\bigr)\|\mathbf r_u\times\mathbf r_v\|\,du\,dv.

For a graph z=f(x,y)z=f(x,y), the formula is

∬Sg dS=∬Dg(x,y,f(x,y))1+fx2+fy2 dA.\iint_S g\,dS=\iint_D g\bigl(x,y,f(x,y)\bigr)\sqrt{1+f_x^2+f_y^2}\,dA.

If g=1g=1, the result is surface area. If gg is a surface density, the result is total mass.

For the plane z=x+yz=x+y over the unit square 0≤x,y≤10\le x,y\le1, with density ρ(x,y,z)=z\rho(x,y,z)=z, one has fx=fy=1f_x=f_y=1 and dS=3 dx dydS=\sqrt{3}\,dx\,dy. Substituting z=x+yz=x+y gives

M=∫01∫01(x+y)3 dy dx=3.M=\int_0^1\int_0^1(x+y)\sqrt{3}\,dy\,dx=\sqrt{3}.

Because only a magnitude appears in the area factor, reversing the does not change a .

Takeaway: Substitute the surface into the scalar function, then multiply by the .

Through an Oriented Surface

Let F(x,y,z)=⟨P(x,y,z),Q(x,y,z),R(x,y,z)⟩\mathbf F(x,y,z)=\langle P(x,y,z),Q(x,y,z),R(x,y,z)\rangle be a vector field on an oriented surface. Its is

∬SF⋅N dS.\iint_S\mathbf F\cdot\mathbf N\,dS.

Using a parametrization, the computational formula is

∬SF⋅N dS=∬DF(r(u,v))⋅(ru×rv) du dv.\iint_S\mathbf F\cdot\mathbf N\,dS=\iint_D\mathbf F\bigl(\mathbf r(u,v)\bigr)\cdot(\mathbf r_u\times\mathbf r_v)\,du\,dv.

The cross product must have the required . If the opposite is needed, use rv×ru\mathbf r_v\times\mathbf r_u, which is the negative of ru×rv\mathbf r_u\times\mathbf r_v.

For a graph z=f(x,y)z=f(x,y) with upward ,

N dS=⟨−fx,−fy,1⟩ dx dy.\mathbf N\,dS=\langle-f_x,-f_y,1\rangle\,dx\,dy.

Thus,

∬SF⋅N dS=∬DF(x,y,f(x,y))⋅⟨−fx,−fy,1⟩ dA.\iint_S\mathbf F\cdot\mathbf N\,dS=\iint_D\mathbf F\bigl(x,y,f(x,y)\bigr)\cdot\langle-f_x,-f_y,1\rangle\,dA.

For example, let F=⟨x,y,z⟩\mathbf F=\langle x,y,z\rangle and z=1−x−yz=1-x-y over D={(x,y):x≥0, y≥0, x+y≤1}D=\{(x,y):x\ge0,\ y\ge0,\ x+y\le1\}. Since fx=fy=−1f_x=f_y=-1, the upward is ⟨1,1,1⟩ dx dy\langle1,1,1\rangle\,dx\,dy. On the surface,

F⋅⟨1,1,1⟩=x+y+(1−x−y)=1.\mathbf F\cdot\langle1,1,1\rangle=x+y+(1-x-y)=1.

Therefore, the upward equals the area of the triangular domain:

∬D1 dA=12.\iint_D1\,dA=\frac{1}{2}.

Takeaway: uses the oriented cross product directly; replacing it with its magnitude loses the direction information.

A Reliable Solution Strategy

A reliable workflow keeps the geometry and the integration factor separate:

  1. Describe the surface and specify its parameter domain.

  2. Choose r(u,v)\mathbf r(u,v), or use a graph representation when the surface is given by z=f(x,y)z=f(x,y).

  3. Compute ru\mathbf r_u and rv\mathbf r_v.

  4. Compute ru×rv\mathbf r_u\times\mathbf r_v.

  5. Check whether its direction matches the required ; reverse it if necessary.

  6. Substitute the parametrization into the integrand.

  7. Use ∥ru×rv∥\|\mathbf r_u\times\mathbf r_v\| for a and ru×rv\mathbf r_u\times\mathbf r_v itself for .

  8. Integrate over the parameter domain.

The central distinction is:

  • Scalar surface integrals measure quantities distributed over area and are independent of .

  • integrals measure passage through a surface and change sign when is reversed.

Final takeaway: Identify the surface, compute its tangent directions, choose the correct area factor, and verify before integrating.