Skip to content

Curves and Surfaces

5.1 Parametric Curves

A parametric curve in R3\mathbb{R}^3 is a C1C^1 function r:[a,b]R3\mathbf{r} : [a, b] \to \mathbb{R}^3 Written r(t)=(x(t),y(t),z(t))\mathbf{r}(t) = (x(t),\, y(t),\, z(t)).

Definition. The arc length of r\mathbf{r} over [a,b][a, b] is

L=abr"(t)dt=ab(dxdt)2+(dydt)2+(dzdt)2dtL = \int_a^b \lVert \mathbf{r}"(t) \rVert\, dt = \int_a^b \sqrt{\left(\frac{dx}{dt}\right)^2 + \left(\frac{dy}{dt}\right)^2 + \left(\frac{dz}{dt}\right)^2}\, dt

Proposition 5.1. The arc length function s(t)=atr(τ)dτs(t) = \int_a^t \lVert \mathbf{r}'(\tau) \rVert\, d\tau Satisfies dsdt=r(t)\frac{ds}{dt} = \lVert \mathbf{r}'(t) \rVertAnd reparametrising by arc length gives a Unit-speed curve: drds=1\lVert \frac{d\mathbf{r}}{ds} \rVert = 1.

Proof. By the Fundamental Theorem of Calculus, dsdt=r(t)\frac{ds}{dt} = \lVert \mathbf{r}'(t) \rVert. If we reparametrise by ssI.e., write r(s)=r(t(s))\mathbf{r}(s) = \mathbf{r}(t(s))Then by the chain rule drds=r(t)dtds\frac{d\mathbf{r}}{ds} = \mathbf{r}'(t) \cdot \frac{dt}{ds}So drds=r(t)dtds=1\lVert \frac{d\mathbf{r}}{ds} \rVert = \lVert \mathbf{r}'(t) \rVert \cdot \left\lvert \frac{dt}{ds} \right\rvert = 1. \blacksquare

Problem. Find the arc length of the curve r(t)=(etcost,etsint,et)\mathbf{r}(t) = (e^t \cos t,\, e^t \sin t,\, e^t) for 0tln20 \leq t \leq \ln 2.

Solution

r(t)=(etcostetsint,etsint+etcost,et)\mathbf{r}'(t) = (e^t \cos t - e^t \sin t,\, e^t \sin t + e^t \cos t,\, e^t)

r(t)2=e2t(costsint)2+e2t(sint+cost)2+e2t\lVert \mathbf{r}'(t) \rVert^2 = e^{2t}(\cos t - \sin t)^2 + e^{2t}(\sin t + \cos t)^2 + e^{2t}

=e2t[(cos2t2sintcost+sin2t)+(sin2t+2sintcost+cos2t)+1]= e^{2t}[(\cos^2 t - 2\sin t \cos t + \sin^2 t) + (\sin^2 t + 2\sin t \cos t + \cos^2 t) + 1]

=e2t[1+1+1]=3e2t= e^{2t}[1 + 1 + 1] = 3e^{2t}

r(t)=3et\lVert \mathbf{r}'(t) \rVert = \sqrt{3}\, e^t

L=0ln23etdt=3[et]0ln2=3(21)=3L = \int_0^{\ln 2} \sqrt{3}\, e^t\, dt = \sqrt{3}\, [e^t]_0^{\ln 2} = \sqrt{3}(2 - 1) = \sqrt{3}

\blacksquare

Problem. Find the arc length of the helix r(t)=(cost,sint,t)\mathbf{r}(t) = (\cos t,\, \sin t,\, t) for 0t4π0 \leq t \leq 4\pi.

Solution

r(t)=(sint,cost,1)\mathbf{r}'(t) = (-\sin t,\, \cos t,\, 1)So r(t)=sin2t+cos2t+1=2\lVert \mathbf{r}'(t) \rVert = \sqrt{\sin^2 t + \cos^2 t + 1} = \sqrt{2}.

L=04π2dt=42πL = \int_0^{4\pi} \sqrt{2}\, dt = 4\sqrt{2}\,\pi

\blacksquare

5.2 Curvature and Torsion

Definition. Let r(s)\mathbf{r}(s) be a unit-speed curve (r(s)=1\lVert \mathbf{r}'(s) \rVert = 1). Define:

  • Unit tangent vector: T(s)=r(s)\mathbf{T}(s) = \mathbf{r}'(s)
  • Curvature: κ(s)=T(s)=r(s)\kappa(s) = \lVert \mathbf{T}'(s) \rVert = \lVert \mathbf{r}''(s) \rVert
  • Principal normal: N(s)=T(s)T(s)\mathbf{N}(s) = \frac{\mathbf{T}'(s)}{\lVert \mathbf{T}'(s) \rVert} (when κ0\kappa \neq 0)
  • Binormal: B(s)=T(s)×N(s)\mathbf{B}(s) = \mathbf{T}(s) \times \mathbf{N}(s)
  • Torsion: τ(s)=B(s)N(s)\tau(s) = -\mathbf{B}'(s) \cdot \mathbf{N}(s)

The vectors T\mathbf{T}, N\mathbf{N}, B\mathbf{B} form the Frenet—Serret frame, an orthonormal Basis that moves with the curve.

Theorem 5.2 (Frenet—Serret Formulas).

T=κN,N=κT+τB,B=τN\mathbf{T}' = \kappa\, \mathbf{N}, \quad \mathbf{N}' = -\kappa\, \mathbf{T} + \tau\, \mathbf{B}, \quad \mathbf{B}' = -\tau\, \mathbf{N}

Proof. Since T\mathbf{T} is a unit vector, TT=1\mathbf{T} \cdot \mathbf{T} = 1So TT=0\mathbf{T}' \cdot \mathbf{T} = 0. Therefore T\mathbf{T}' is orthogonal to T\mathbf{T}So T\mathbf{T}' is parallel to N\mathbf{N} (when κ0\kappa \neq 0). This gives T=κN\mathbf{T}' = \kappa\,\mathbf{N}.

Similarly, B=T×N\mathbf{B} = \mathbf{T} \times \mathbf{N} is a unit vector, so BB=0\mathbf{B}' \cdot \mathbf{B} = 0. Also BT=0\mathbf{B} \cdot \mathbf{T} = 0So BT+BT=0\mathbf{B}' \cdot \mathbf{T} + \mathbf{B} \cdot \mathbf{T}' = 0 Giving BT=BκN=0\mathbf{B}' \cdot \mathbf{T} = -\mathbf{B} \cdot \kappa\,\mathbf{N} = 0. So B\mathbf{B}' is Parallel to N\mathbf{N}Giving B=τN\mathbf{B}' = -\tau\,\mathbf{N}.

For N\mathbf{N}': since {T,N,B}\{\mathbf{T}, \mathbf{N}, \mathbf{B}\} is an orthonormal basis, N=(NT)T+(NN)N+(NB)B\mathbf{N}' = (\mathbf{N}' \cdot \mathbf{T})\,\mathbf{T} + (\mathbf{N}' \cdot \mathbf{N})\,\mathbf{N} + (\mathbf{N}' \cdot \mathbf{B})\,\mathbf{B}. From NT=0\mathbf{N} \cdot \mathbf{T} = 0: NT=NT=κ\mathbf{N}' \cdot \mathbf{T} = -\mathbf{N} \cdot \mathbf{T}' = -\kappa. From NN=1\mathbf{N} \cdot \mathbf{N} = 1: NN=0\mathbf{N}' \cdot \mathbf{N} = 0. From NB=0\mathbf{N} \cdot \mathbf{B} = 0: NB=NB=τ\mathbf{N}' \cdot \mathbf{B} = -\mathbf{N} \cdot \mathbf{B}' = \tau. This gives N=κT+τB\mathbf{N}' = -\kappa\,\mathbf{T} + \tau\,\mathbf{B}. \blacksquare

Intuition. The curvature κ\kappa measures how sharply the curve bends (deviation from a straight line). The torsion τ\tau measures how sharply the curve twists out of the osculating plane (deviation from a Plane curve). A curve lies in a plane if and only if τ=0\tau = 0 everywhere.

For a curve parameterised by an arbitrary parameter tt (not necessarily unit-speed):

κ=r(t)×r(t)r(t)3\kappa = \frac{\lVert \mathbf{r}'(t) \times \mathbf{r}''(t) \rVert}{\lVert \mathbf{r}'(t) \rVert^3}

τ=[r(t)×r(t)]r(t)r(t)×r(t)2\tau = \frac{[\mathbf{r}'(t) \times \mathbf{r}''(t)] \cdot \mathbf{r}^{\prime\prime\prime}(t)}{\lVert \mathbf{r}'(t) \times \mathbf{r}''(t) \rVert^2}

Problem. Find the curvature and torsion of the helix r(t)=(acost,asint,bt)\mathbf{r}(t) = (a\cos t,\, a\sin t,\, bt) where a,b>0a, b \gt 0.

Solution

r(t)=(asint,acost,b)\mathbf{r}'(t) = (-a\sin t,\, a\cos t,\, b)

r(t)=(acost,asint,0)\mathbf{r}''(t) = (-a\cos t,\, -a\sin t,\, 0)

r(t)=(asint,acost,0)\mathbf{r}^{\prime\prime\prime}(t) = (a\sin t,\, -a\cos t,\, 0)

r=a2+b2\lVert \mathbf{r}' \rVert = \sqrt{a^2 + b^2}

r×r=ijkasintacostbacostasint0=(absint,abcost,a2)\mathbf{r}' \times \mathbf{r}'' = \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \\ -a\sin t & a\cos t & b \\ -a\cos t & -a\sin t & 0 \end{vmatrix} = (ab\sin t,\, -ab\cos t,\, a^2)

r×r=a2b2+a4=aa2+b2\lVert \mathbf{r}' \times \mathbf{r}'' \rVert = \sqrt{a^2 b^2 + a^4} = a\sqrt{a^2 + b^2}

κ=aa2+b2(a2+b2)3/2=aa2+b2\kappa = \frac{a\sqrt{a^2 + b^2}}{(a^2 + b^2)^{3/2}} = \frac{a}{a^2 + b^2}

For the torsion:

(r×r)r=absintasint+(abcost)(acost)+a20=a2b(\mathbf{r}' \times \mathbf{r}'') \cdot \mathbf{r}^{\prime\prime\prime} = ab\sin t \cdot a\sin t + (-ab\cos t)(-a\cos t) + a^2 \cdot 0 = a^2 b

τ=a2ba2(a2+b2)=ba2+b2\tau = \frac{a^2 b}{a^2(a^2 + b^2)} = \frac{b}{a^2 + b^2}

\blacksquare

Remark. The helix has constant curvature and constant torsion, reflecting its uniform geometry.

5.3 Parametric Surfaces

A parametric surface is a C1C^1 map r:DR2R3\mathbf{r} : D \subseteq \mathbb{R}^2 \to \mathbb{R}^3 r(u,v)=(x(u,v),y(u,v),z(u,v))\mathbf{r}(u, v) = (x(u,v),\, y(u,v),\, z(u,v)).

The tangent plane at r(u0,v0)\mathbf{r}(u_0, v_0) is spanned by the tangent vectors

ru=ru,rv=rv\mathbf{r}_u = \frac{\partial \mathbf{r}}{\partial u}, \quad \mathbf{r}_v = \frac{\partial \mathbf{r}}{\partial v}

The unit normal to the surface is

n=ru×rvru×rv\mathbf{n} = \frac{\mathbf{r}_u \times \mathbf{r}_v}{\lVert \mathbf{r}_u \times \mathbf{r}_v \rVert}

Examples of parametric surfaces:

  • Sphere (spherical coordinates): r(θ,ϕ)=(ρsinϕcosθ,ρsinϕsinθ,ρcosϕ)\mathbf{r}(\theta, \phi) = (\rho\sin\phi\cos\theta,\, \rho\sin\phi\sin\theta,\, \rho\cos\phi)
  • Cylinder: r(θ,z)=(rcosθ,rsinθ,z)\mathbf{r}(\theta, z) = (r\cos\theta,\, r\sin\theta,\, z)
  • Graph of z=f(x,y)z = f(x,y): r(x,y)=(x,y,f(x,y))\mathbf{r}(x, y) = (x,\, y,\, f(x,y))

For the graph z=f(x,y)z = f(x,y)The normal is n=(fx,fy,1)1+fx2+fy2\mathbf{n} = \frac{(-f_x,\, -f_y,\, 1)}{\sqrt{1 + f_x^2 + f_y^2}}.

5.4 Surface Area

Definition. The area of a parametric surface r:DR3\mathbf{r} : D \to \mathbb{R}^3 is

A(S)=Dru×rvdudvA(S) = \iint_D \lVert \mathbf{r}_u \times \mathbf{r}_v \rVert\, du\, dv

Derivation. Partition DD into small rectangles DijD_{ij} of area ΔuΔv\Delta u\, \Delta v. The image r(Dij)\mathbf{r}(D_{ij}) is approximately a parallelogram spanned by ruΔu\mathbf{r}_u\, \Delta u and rvΔv\mathbf{r}_v\, \Delta vWith area ru×rvΔuΔv\lVert \mathbf{r}_u \times \mathbf{r}_v \rVert\, \Delta u\, \Delta v. Summing and taking the limit gives the formula. \blacksquare

Problem. Find the surface area of the part of the paraboloid z=x2+y2z = x^2 + y^2 that lies below The plane z=4z = 4.

Solution

Parametrise by r(x,y)=(x,y,x2+y2)\mathbf{r}(x, y) = (x,\, y,\, x^2 + y^2) where x2+y24x^2 + y^2 \leq 4.

rx=(1,0,2x)\mathbf{r}_x = (1,\, 0,\, 2x), ry=(0,1,2y)\mathbf{r}_y = (0,\, 1,\, 2y).

rx×ry=ijk102x012y=(2x,2y,1)\mathbf{r}_x \times \mathbf{r}_y = \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \\ 1 & 0 & 2x \\ 0 & 1 & 2y \end{vmatrix} = (-2x,\, -2y,\, 1)

rx×ry=4x2+4y2+1\lVert \mathbf{r}_x \times \mathbf{r}_y \rVert = \sqrt{4x^2 + 4y^2 + 1}

A=x2+y244x2+4y2+1dxdyA = \iint_{x^2+y^2 \leq 4} \sqrt{4x^2 + 4y^2 + 1}\, dx\, dy

Use polar coordinates: x=rcosθx = r\cos\theta, y=rsinθy = r\sin\theta, 0r20 \leq r \leq 2, 0θ2π0 \leq \theta \leq 2\pi.

A=02π024r2+1rdrdθA = \int_0^{2\pi} \int_0^2 \sqrt{4r^2 + 1}\, r\, dr\, d\theta

Let u=4r2+1u = 4r^2 + 1, du=8rdrdu = 8r\, dr:

=2π18117udu=π4[2u3/23]117=π6(173/21)= 2\pi \cdot \frac{1}{8} \int_1^{17} \sqrt{u}\, du = \frac{\pi}{4}\left[\frac{2u^{3/2}}{3}\right]_1^{17} = \frac{\pi}{6}(17^{3/2} - 1)

\blacksquare

5.5 Surface Integrals

Definition (Scalar surface integral). The integral of a scalar function ff over a parametric Surface SS is

SfdS=Df(r(u,v))ru×rvdudv\iint_S f\, dS = \iint_D f(\mathbf{r}(u,v))\, \lVert \mathbf{r}_u \times \mathbf{r}_v \rVert\, du\, dv

Definition (Vector surface integral / flux). The flux of a vector field F\mathbf{F} through an Oriented surface SS is

SFdS=DF(r(u,v))(ru×rv)dudv\iint_S \mathbf{F} \cdot d\mathbf{S} = \iint_D \mathbf{F}(\mathbf{r}(u,v)) \cdot (\mathbf{r}_u \times \mathbf{r}_v)\, du\, dv

Where the orientation is determined by the choice of normal ru×rv\mathbf{r}_u \times \mathbf{r}_v vs. rv×ru\mathbf{r}_v \times \mathbf{r}_u.

Problem. Evaluate SFdS\iint_S \mathbf{F} \cdot d\mathbf{S} where F=(x,y,z2)\mathbf{F} = (x,\, y,\, z^2) and SS is the hemisphere x2+y2+z2=4x^2 + y^2 + z^2 = 4, z0z \geq 0With Upward orientation.

Solution

Use the divergence theorem on the closed hemisphere plus the disk at z=0z = 0. Let EE be the solid hemisphere. Then:

closed SFdS=EFdV=E(1+1+2z)dV\iint_{\mathrm{closed\ S} \mathbf{F} \cdot d\mathbf{S} = \iiint_E \nabla \cdot \mathbf{F}\, dV = \iiint_E (1 + 1 + 2z)\, dV}

=2V+2EzdV= 2V + 2\iiint_E z\, dV

Where V=1243π(23)=16π3V = \frac{1}{2} \cdot \frac{4}{3}\pi(2^3) = \frac{16\pi}{3}.

By symmetry, the centroid of a hemisphere of radius R=2R = 2 is at z=3R/8=3/4z = 3R/8 = 3/4So

EzdV=zˉV=3416π3=4π\iiint_E z\, dV = \bar{z} \cdot V = \frac{3}{4} \cdot \frac{16\pi}{3} = 4\pi

=216π3+24π=32π3+8π=56π3= 2 \cdot \frac{16\pi}{3} + 2 \cdot 4\pi = \frac{32\pi}{3} + 8\pi = \frac{56\pi}{3}

The flux through the disk z=0z = 0, x2+y24x^2 + y^2 \leq 4 (with downward normal k-\mathbf{k}): diskF(k)dS=disk0dS=0\iint_{\mathrm{disk} \mathbf{F} \cdot (-\mathbf{k})\, dS = \iint_{\mathrm{disk} 0\, dS = 0}}.

So the flux through the hemisphere alone is 56π3\frac{56\pi}{3}. \blacksquare

Problem. Evaluate SzdS\iint_S z\, dS where SS is the part of the plane 2x+2y+z=42x + 2y + z = 4 in the First octant.

Solution

Parametrise the surface. Solve for z=42x2yz = 4 - 2x - 2y where x0x \geq 0, y0y \geq 0, z0z \geq 0 I.e., 2x+2y42x + 2y \leq 4 or x+y2x + y \leq 2.

r(x,y)=(x,y,42x2y)\mathbf{r}(x,y) = (x,\, y,\, 4 - 2x - 2y), D=(x,y):x0,y0,x+y2D = \\{(x,y) : x \geq 0,\, y \geq 0,\, x + y \leq 2\\}.

rx=(1,0,2)\mathbf{r}_x = (1,\, 0,\, -2), ry=(0,1,2)\mathbf{r}_y = (0,\, 1,\, -2).

rx×ry=(2,2,1)\mathbf{r}_x \times \mathbf{r}_y = (2,\, 2,\, 1)

rx×ry=4+4+1=3\lVert \mathbf{r}_x \times \mathbf{r}_y \rVert = \sqrt{4 + 4 + 1} = 3

SzdS=D(42x2y)3dxdy=30202x(42x2y)dydx\iint_S z\, dS = \iint_D (4 - 2x - 2y) \cdot 3\, dx\, dy = 3 \int_0^2 \int_0^{2-x} (4 - 2x - 2y)\, dy\, dx

=302[(42x)yy2]02xdx=302(2x)(2x)dx= 3 \int_0^2 \left[(4-2x)y - y^2\right]_0^{2-x}\, dx = 3 \int_0^2 (2-x)(2-x)\, dx

=302(2x)2dx=3[(2x)33]02=383=8= 3 \int_0^2 (2-x)^2\, dx = 3\left[-\frac{(2-x)^3}{3}\right]_0^2 = 3 \cdot \frac{8}{3} = 8

\blacksquare

5.6 Common Pitfalls

:::caution Common Pitfalls

  • Parameterisation domain. Always verify that the parameterisation covers the entire surface and that the map is one-to-one (except possibly on the boundary).
  • Normal orientation. The cross product ru×rv\mathbf{r}_u \times \mathbf{r}_v determines the orientation. Swapping the order changes the sign of the flux integral.
  • Surface area vs. Flux. Surface area uses ru×rv\lVert \mathbf{r}_u \times \mathbf{r}_v \rVert (scalar), while flux uses ru×rv\mathbf{r}_u \times \mathbf{r}_v (vector, oriented).

:::