Skip to content

Curves and Surfaces

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=∫ab∥r"(t)∥ dt=∫ab(dxdt)2+(dydt)2+(dzdt)2 dtL = \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)=∫at∥r′(τ)∥ 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)=(etcos⁡t, etsin⁡t, et)\mathbf{r}(t) = (e^t \cos t,\, e^t \sin t,\, e^t) for 0≤t≤ln⁡20 \leq t \leq \ln 2.

Solution

r′(t)=(etcos⁡t−etsin⁡t, etsin⁡t+etcos⁡t, 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(cos⁡t−sin⁡t)2+e2t(sin⁡t+cos⁡t)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[(cos⁡2t−2sin⁡tcos⁡t+sin⁡2t)+(sin⁡2t+2sin⁡tcos⁡t+cos⁡2t)+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)∥=3 et\lVert \mathbf{r}'(t) \rVert = \sqrt{3}\, e^t

L=∫0ln⁡23 et dt=3 [et]0ln⁡2=3(2−1)=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)=(cos⁡t, sin⁡t, t)\mathbf{r}(t) = (\cos t,\, \sin t,\, t) for 0≤t≤4π0 \leq t \leq 4\pi.

Solution

r′(t)=(−sin⁡t, cos⁡t, 1)\mathbf{r}'(t) = (-\sin t,\, \cos t,\, 1)So ∥r′(t)∥=sin⁡2t+cos⁡2t+1=2\lVert \mathbf{r}'(t) \rVert = \sqrt{\sin^2 t + \cos^2 t + 1} = \sqrt{2}.

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

■\blacksquare

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, T⋅T=1\mathbf{T} \cdot \mathbf{T} = 1So T′⋅T=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 B′⋅B=0\mathbf{B}' \cdot \mathbf{B} = 0. Also B⋅T=0\mathbf{B} \cdot \mathbf{T} = 0So B′⋅T+B⋅T′=0\mathbf{B}' \cdot \mathbf{T} + \mathbf{B} \cdot \mathbf{T}' = 0 Giving B′⋅T=−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′=(N′⋅T) T+(N′⋅N) N+(N′⋅B) 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 N⋅T=0\mathbf{N} \cdot \mathbf{T} = 0: N′⋅T=−N⋅T′=−κ\mathbf{N}' \cdot \mathbf{T} = -\mathbf{N} \cdot \mathbf{T}' = -\kappa. From N⋅N=1\mathbf{N} \cdot \mathbf{N} = 1: N′⋅N=0\mathbf{N}' \cdot \mathbf{N} = 0. From N⋅B=0\mathbf{N} \cdot \mathbf{B} = 0: N′⋅B=−N⋅B′=τ\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)=(acos⁡t, asin⁡t, bt)\mathbf{r}(t) = (a\cos t,\, a\sin t,\, bt) where a,b>0a, b \gt 0.

Solution

r′(t)=(−asin⁡t, acos⁡t, b)\mathbf{r}'(t) = (-a\sin t,\, a\cos t,\, b)

r′′(t)=(−acos⁡t, −asin⁡t, 0)\mathbf{r}''(t) = (-a\cos t,\, -a\sin t,\, 0)

r′′′(t)=(asin⁡t, −acos⁡t, 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′′=∣ijk−asin⁡tacos⁡tb−acos⁡t−asin⁡t0∣=(absin⁡t, −abcos⁡t, 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′′′=absin⁡t⋅asin⁡t+(−abcos⁡t)(−acos⁡t)+a2⋅0=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.

A parametric surface is a C1C^1 map r:D⊆R2→R3\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=∂r∂u,rv=∂r∂v\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×rv∥ru×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}}.

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

A(S)=∬D∥ru×rv∥ du dvA(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+y2≤4x^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+y2≤44x2+4y2+1 dx dyA = \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, 0≤r≤20 \leq r \leq 2, 0≤θ≤2π0 \leq \theta \leq 2\pi.

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

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

=2π⋅18∫117u du=π4[2u3/23]117=π6(173/2−1)= 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

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

∬Sf dS=∬Df(r(u,v)) ∥ru×rv∥ du dv\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

∬SF⋅dS=∬DF(r(u,v))⋅(ru×rv) du dv\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 ∬SF⋅dS\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, z≥0z \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 SF⋅dS=∭E∇⋅F dV=∭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+2∭Ez dV= 2V + 2\iiint_E z\, dV

Where V=12⋅43π(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

∭Ez dV=zˉ⋅V=34⋅16π3=4π\iiint_E z\, dV = \bar{z} \cdot V = \frac{3}{4} \cdot \frac{16\pi}{3} = 4\pi

=2⋅16π3+2⋅4π=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+y2≤4x^2 + y^2 \leq 4 (with downward normal −k-\mathbf{k}): ∬diskF⋅(−k) dS=∬disk0 dS=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 ∬Sz dS\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=4−2x−2yz = 4 - 2x - 2y where x≥0x \geq 0, y≥0y \geq 0, z≥0z \geq 0 I.e., 2x+2y≤42x + 2y \leq 4 or x+y≤2x + y \leq 2.

r(x,y)=(x, y, 4−2x−2y)\mathbf{r}(x,y) = (x,\, y,\, 4 - 2x - 2y), D=(x,y):x≥0, y≥0, x+y≤2D = \\{(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

∬Sz dS=∬D(4−2x−2y)⋅3 dx dy=3∫02∫02−x(4−2x−2y) dy dx\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

=3∫02[(4−2x)y−y2]02−x dx=3∫02(2−x)(2−x) 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

=3∫02(2−x)2 dx=3[−(2−x)33]02=3⋅83=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

:::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).

:::