Introduction
In this part we’ll cover surface integrals over vector fields as well as two very important theorems.
Surface integrals over vector fields
Definition 1 (Flux) Let S S S be a parametric surface, parameterized by, r ⃗ : D → R 3 \vec{r}: D \rarr \mathbb{R}^3 r : D → R 3 and F : R 3 → R 3 F: \mathbb{R}^3 \rarr \mathbb{R}^3 F : R 3 → R 3 , be a vector field. Then:
∬ S F ⃗ ⋅ d S = ∬ D F ⃗ ( r ⃗ ( u , v ) ) ⋅ ( r u ⃗ × r v ⃗ ) d u d v . \iint_S \vec{F} \cdot dS = \iint_D \vec{F}(\vec{r}(u, v)) \cdot (\vec{r_u} \times \vec{r_v}) \, du \, dv. ∬ S F ⋅ d S = ∬ D F ( r ( u , v )) ⋅ ( r u × r v ) d u d v . We call this the flux of F ⃗ \vec{F} F through S S S .
Example 1 Given, F ⃗ = ⟨ z , y , x ⟩ \vec{F} = \langle z, y, x \rangle F = ⟨ z , y , x ⟩ , find the flux of this vector field across the unit sphere.
So, we have the surface:
S = { ( x , y , z ) ∣ x 2 + y 2 + z 2 = 1 } S = \{(x, y, z) \mid x^2 + y^2 + z^2 = 1 \} S = {( x , y , z ) ∣ x 2 + y 2 + z 2 = 1 } As always, good idea to change into polar coordinates, so a good parameterization is:
r ⃗ ( ϕ , θ ) = ⟨ sin ( ϕ ) cos ( θ ) , sin ( ϕ ) sin ( θ ) , cos ( ϕ ) ⟩ , 0 ≤ ϕ ≤ π , 0 ≤ θ ≤ 2 π . \vec{r}(\phi, \theta) = \langle \sin(\phi) \cos(\theta), \sin(\phi) \sin(\theta), \cos(\phi) \rangle, \quad 0 \leq \phi \leq \pi, 0 \leq \theta \leq 2\pi. r ( ϕ , θ ) = ⟨ sin ( ϕ ) cos ( θ ) , sin ( ϕ ) sin ( θ ) , cos ( ϕ )⟩ , 0 ≤ ϕ ≤ π , 0 ≤ θ ≤ 2 π . NB: This computation is just tedious (and we’ll see how we can reduce this shortly)
r ϕ ⃗ × r θ ⃗ = ⟨ sin 2 ( ϕ ) cos ( θ ) , sin 2 ( ϕ ) sin ( θ ) , sin ( ϕ ) cos ( ϕ ) ⟩ \vec{r_{\phi}} \times \vec{r_{\theta}} = \langle \sin^2(\phi) \cos(\theta), \sin^2(\phi) \sin(\theta), \sin(\phi) \cos(\phi) \rangle r ϕ × r θ = ⟨ sin 2 ( ϕ ) cos ( θ ) , sin 2 ( ϕ ) sin ( θ ) , sin ( ϕ ) cos ( ϕ )⟩ Now, one can solve this integral by taking the dot product and applying our definition:
∫ 0 2 π ∫ 0 π ⟨ cos ( ϕ ) , sin ( ϕ ) sin ( θ ) , sin ( ϕ ) cos ( θ ) ⟩ ⋅ ⟨ sin 2 ( ϕ ) cos ( θ ) , sin 2 ( ϕ ) sin ( θ ) , sin ( ϕ ) cos ( ϕ ) ⟩ d ϕ d θ \int_0^{2\pi} \int_0^{\pi} \langle \cos(\phi), \sin(\phi) \sin(\theta), \sin(\phi) \cos(\theta) \rangle \cdot \langle \sin^2(\phi) \cos(\theta), \sin^2(\phi) \sin(\theta), \sin(\phi) \cos(\phi) \rangle \, d\phi \, d\theta ∫ 0 2 π ∫ 0 π ⟨ cos ( ϕ ) , sin ( ϕ ) sin ( θ ) , sin ( ϕ ) cos ( θ )⟩ ⋅ ⟨ sin 2 ( ϕ ) cos ( θ ) , sin 2 ( ϕ ) sin ( θ ) , sin ( ϕ ) cos ( ϕ )⟩ d ϕ d θ Let’s see now how we can do the same example, with fewer computations, first some theory.
Orientation
For surfaces, we say that positive orientation is “pointing outwards”. For hollow surfaces this is always the case. If our parameterization is given by a single function, then we always have the correct choice.
However, imagine half the unit sphere with a filled base. The base with parameterization: D = { ( r , θ ) ∣ 0 ≤ r ≤ 1 , 0 ≤ θ ≤ 2 π } D = \{(r, \theta) \mid 0 \leq r \leq 1, 0 \leq \theta \leq 2\pi \} D = {( r , θ ) ∣ 0 ≤ r ≤ 1 , 0 ≤ θ ≤ 2 π } . Is going to be a negative integral.
Definition 2 (Divergence) For a vector field, F ⃗ : R 3 → R 3 \vec{F}: \mathbb{R}^3 \rarr \mathbb{R}^3 F : R 3 → R 3 , the divergence is:
div F ⃗ = ∂ P ∂ x + ∂ Q ∂ y + ∂ R ∂ z . \operatorname{div} \vec{F} = \dfrac{\partial P}{\partial x} + \dfrac{\partial Q}{\partial y} + \dfrac{\partial R}{\partial z}. div F = ∂ x ∂ P + ∂ y ∂ Q + ∂ z ∂ R .
Theorem 1 (The Divergence Theorem) Let S S S be a hollow surface and let E E E be the 3-dimensional region it surrounds (∂ E = S \partial E = S ∂ E = S ).
Then:
∬ S F ⃗ ⋅ d S = ∭ E div F ⃗ d V . \iint_S \vec{F} \cdot dS = \iiint_E \operatorname{div} \vec{F} \, dV. ∬ S F ⋅ d S = ∭ E div F d V . Equivalently:
∬ ∂ E F ⃗ ⋅ d S = ∭ E div F ⃗ d V . \iint_{\partial E} \vec{F} \cdot dS = \iiint_E \operatorname{div} \vec{F} \, dV. ∬ ∂ E F ⋅ d S = ∭ E div F d V .
Now let’s take the same example.
Example 2 Find the flux of F ⃗ = ⟨ z , y , x ⟩ \vec{F} = \langle z, y, x \rangle F = ⟨ z , y , x ⟩ on the unit sphere.
div F ⃗ = ∂ P ∂ x + ∂ Q ∂ y + ∂ R ∂ z = 0 + 1 + 0 = 1 \operatorname{div} \vec{F} = \dfrac{\partial P}{\partial x} + \dfrac{\partial Q}{\partial y} + \dfrac{\partial R}{\partial z} = 0 + 1 + 0 = 1 div F = ∂ x ∂ P + ∂ y ∂ Q + ∂ z ∂ R = 0 + 1 + 0 = 1 ∬ S F ⃗ ⋅ d S = ∭ E 1 d V \iint_S \vec{F} \cdot dS = \iiint_E 1 \, dV ∬ S F ⋅ d S = ∭ E 1 d V We could now just transform to polar coordinates, but, this is the volume of the unit sphere, which is just:
4 π 3 \boxed{\dfrac{4\pi}{3}} 3 4 π
Stokes’ Theorem
So far we have seen Green’s Theorem as well as The Divergence theorem. But what is the reason that Green’s Theorem only works in R 2 \mathbb{R}^2 R 2 ?
Well, Stokes’ Theorem is essentially Green’s Theorem but in R 3 \mathbb{R}^3 R 3 .
Definition 3 (Curl) For F ⃗ : R 3 → R 3 \vec{F}: \mathbb{R}^3 \rarr \mathbb{R}^3 F : R 3 → R 3 , we define its curl as:
NB: it is also called rot and sometimes denoted as ∇ × F ⃗ \nabla \times \vec{F} ∇ × F .
curl F ⃗ = ⟨ ∂ R ∂ y − ∂ Q ∂ z , ∂ P ∂ z − ∂ R ∂ x , ∂ Q ∂ x − ∂ P ∂ y ⟩ . \operatorname{curl} \vec{F} = \left\langle \dfrac{\partial R}{\partial y} - \dfrac{\partial Q}{\partial z}, \dfrac{\partial P}{\partial z} - \dfrac{\partial R}{\partial x}, \dfrac{\partial Q}{\partial x} - \dfrac{\partial P}{\partial y} \right\rangle. curl F = ⟨ ∂ y ∂ R − ∂ z ∂ Q , ∂ z ∂ P − ∂ x ∂ R , ∂ x ∂ Q − ∂ y ∂ P ⟩ .
Theorem 2 (Stokes’ Theorem) If S ⊂ R 3 S \subset \mathbb{R}^3 S ⊂ R 3 is a surface, surrounded by a simple, closed curve in R 3 \mathbb{R}^3 R 3 . Then:
∫ C F ⃗ ⋅ d r = ∬ S curl F ⃗ ⋅ d S . \int_C \vec{F} \cdot dr = \iint_S \operatorname{curl} \vec{F} \cdot dS. ∫ C F ⋅ d r = ∬ S curl F ⋅ d S .
Example 3 Compute ∫ C F ⋅ d r \int_C F \cdot dr ∫ C F ⋅ d r for F ⃗ = ⟨ − y 2 , x , z 2 ⟩ \vec{F} = \langle -y^2, x, z^2 \rangle F = ⟨ − y 2 , x , z 2 ⟩ and C C C being the intersection of y + z = 2 y + z = 2 y + z = 2 with x 2 + y 2 = 1 x^2 + y^2 = 1 x 2 + y 2 = 1 .
Let’s find the curl:
curl F ⃗ = ⟨ 0 , 0 , 1 + 2 y ⟩ \operatorname{curl} \vec{F} = \langle 0, 0, 1 + 2y \rangle curl F = ⟨ 0 , 0 , 1 + 2 y ⟩ We need to find the surface, the intersection of y + z = 2 y + z = 2 y + z = 2 and x 2 + y 2 = 1 x^2 + y^2 = 1 x 2 + y 2 = 1 will just be the unit disk but at a slight tilt.
∬ S curl F ⃗ ⋅ d S \iint_S \operatorname{curl} \vec{F} \cdot dS ∬ S curl F ⋅ d S Is just the flux through this region, let’s express this in polar coordinates
D = { ( r , θ ) ∣ 0 ≤ r ≤ 1 , 0 ≤ θ ≤ 2 π } D = \{(r, \theta) \mid 0 \leq r \leq 1, 0 \leq \theta \leq 2\pi \} D = {( r , θ ) ∣ 0 ≤ r ≤ 1 , 0 ≤ θ ≤ 2 π } Meaning:
∬ S curl F ⃗ ⋅ d S = ∬ D curl F ⃗ ⋅ ( r u ⃗ × r v ⃗ ) d A \iint_S \operatorname{curl} \vec{F} \cdot dS = \iint_D \operatorname{curl} \vec{F} \cdot (\vec{r_u} \times \vec{r_v}) \, dA ∬ S curl F ⋅ d S = ∬ D curl F ⋅ ( r u × r v ) d A Let’s find a parameterization of the surface:
r ⃗ ( u , v ) = ⟨ u , v , 2 − v ⟩ , r u ⃗ ( u , v ) = ⟨ 1 , 0 , 0 ⟩ , r v ⃗ ( u , v ) = ⟨ 0 , 1 , − 1 ⟩ , r u ⃗ × r v ⃗ = ⟨ 0 , 1 , 1 ⟩ . \begin{align*}
\vec{r}(u, v) & = \langle u, v, 2 - v \rangle, \newline
\vec{r_u}(u, v) & = \langle 1, 0, 0 \rangle, \newline
\vec{r_v}(u, v) & = \langle 0, 1, -1 \rangle, \newline
\vec{r_u} \times \vec{r_v} & = \langle 0, 1, 1 \rangle.
\end{align*} r ( u , v ) r u ( u , v ) r v ( u , v ) r u × r v = ⟨ u , v , 2 − v ⟩ , = ⟨ 1 , 0 , 0 ⟩ , = ⟨ 0 , 1 , − 1 ⟩ , = ⟨ 0 , 1 , 1 ⟩ . ∬ D curl F ⃗ ⋅ ( r u ⃗ × r v ⃗ ) d A = ∬ D ⟨ 0 , 0 , 1 + 2 y ⟩ ⋅ ⟨ 0 , 1 , 1 ⟩ d A = ∬ D 1 + 2 y d A = ∫ 0 2 π ∫ 0 1 ( 1 + 2 r sin ( θ ) ) r d r d θ = ∫ 0 2 π ∫ 0 1 r + 2 r 2 sin ( θ ) d r d θ = ∫ 0 1 ∫ 0 2 π r + 2 r 2 sin ( θ ) d θ d r \begin{align*}
\iint_D \operatorname{curl} \vec{F} \cdot (\vec{r_u} \times \vec{r_v}) \, dA & = \iint_D \langle 0, 0, 1 + 2y \rangle \cdot \langle 0, 1, 1 \rangle \, dA \newline
& = \iint_D 1 + 2y \, dA \newline
& = \int_0^{2\pi} \int_0^1 (1 + 2r \sin(\theta))r \, dr \, d\theta \newline
& = \int_0^{2\pi} \int_0^1 r + 2r^2 \sin(\theta) \, dr \, d\theta \newline
& = \int_0^1 \int_0^{2\pi} r + 2r^2 \sin(\theta) \, d\theta \, dr
\end{align*} ∬ D curl F ⋅ ( r u × r v ) d A = ∬ D ⟨ 0 , 0 , 1 + 2 y ⟩ ⋅ ⟨ 0 , 1 , 1 ⟩ d A = ∬ D 1 + 2 y d A = ∫ 0 2 π ∫ 0 1 ( 1 + 2 r sin ( θ )) r d r d θ = ∫ 0 2 π ∫ 0 1 r + 2 r 2 sin ( θ ) d r d θ = ∫ 0 1 ∫ 0 2 π r + 2 r 2 sin ( θ ) d θ d r sin ( θ ) \sin(\theta) sin ( θ ) has symmetry around 2 π 2\pi 2 π , therefore the integral of that is 0 0 0 .
∫ 0 1 2 π r d r = π \int_0^1 2\pi r \, dr = \boxed{\pi} ∫ 0 1 2 π r d r = π
Example 4 Compute ∬ S curl F ⃗ ⋅ d S \iint_S \operatorname{curl} \vec{F} \cdot dS ∬ S curl F ⋅ d S where F ⃗ = ⟨ x z , y z , x y ⟩ \vec{F} = \langle xz, yz, xy \rangle F = ⟨ x z , y z , x y ⟩ and S S S is the part of the sphere of radius 2 around ( 0 , 0 , 0 ) (0, 0, 0) ( 0 , 0 , 0 ) that lies inside the cylinder x 2 + y 2 ≤ 1 x^2 + y^2 \leq 1 x 2 + y 2 ≤ 1 and above the x y xy x y -plane.
Let’s actually use Stokes’ theorem, “backwards”. C C C will be a circle with:
{ x 2 + y 2 = 1 x 2 + y 2 + z 2 = 4 \begin{cases}
x^2 + y^2 = 1 \newline
x^2 + y^2 + z^2 = 4
\end{cases} { x 2 + y 2 = 1 x 2 + y 2 + z 2 = 4 Which means that z = 3 z = \sqrt{3} z = 3
We can parameterize this curve as:
r ⃗ ( t ) = ⟨ cos ( t ) , sin ( t ) , 3 ⟩ , 0 ≤ t ≤ 2 π , r ′ ⃗ ( t ) = ⟨ − sin ( t ) , cos ( t ) , 0 ⟩ . \begin{align*}
\vec{r}(t) & = \langle \cos(t), \sin(t), \sqrt{3} \rangle, \quad 0 \leq t \leq 2\pi, \newline
\vec{r^\prime}(t) & = \langle -\sin(t), \cos(t), 0 \rangle.
\end{align*} r ( t ) r ′ ( t ) = ⟨ cos ( t ) , sin ( t ) , 3 ⟩ , 0 ≤ t ≤ 2 π , = ⟨ − sin ( t ) , cos ( t ) , 0 ⟩ . ∬ S curl F ⃗ ⋅ d S = ∫ C F ⃗ ⋅ d r = ∫ 0 2 π ⟨ 3 cos ( t ) , 3 sin ( t ) , cos ( t ) sin ( t ) ⟩ ⋅ ⟨ − sin ( t ) , cos ( t ) , 0 ⟩ d t = ∫ 0 2 π − 3 cos ( t ) sin ( t ) + 3 sin ( t ) cos ( t ) d t = ∫ 0 2 π 0 d t = 0 \begin{align*}
\iint_S \operatorname{curl} \vec{F} \cdot dS & = \int_C \vec{F} \cdot dr \newline
& = \int_0^{2\pi} \langle \sqrt{3} \cos(t), \sqrt{3} \sin(t), \cos(t) \sin(t) \rangle \cdot \langle -\sin(t), \cos(t), 0 \rangle \, dt \newline
& = \int_0^{2\pi} -\sqrt{3} \cos(t) \sin(t) + \sqrt{3} \sin(t) \cos(t) \, dt \newline
& = \int_0^{2\pi} 0 \, dt = \boxed{0}
\end{align*} ∬ S curl F ⋅ d S = ∫ C F ⋅ d r = ∫ 0 2 π ⟨ 3 cos ( t ) , 3 sin ( t ) , cos ( t ) sin ( t )⟩ ⋅ ⟨ − sin ( t ) , cos ( t ) , 0 ⟩ d t = ∫ 0 2 π − 3 cos ( t ) sin ( t ) + 3 sin ( t ) cos ( t ) d t = ∫ 0 2 π 0 d t = 0