Overview

Part 20 - Divergence and Stokes' Theorem

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 SS be a parametric surface, parameterized by, r⃗:D→R3\vec{r}: D \rarr \mathbb{R}^3 and F:R3→R3F: \mathbb{R}^3 \rarr \mathbb{R}^3, be a vector field. Then:

∬SF⃗⋅dS=∬DF⃗(r⃗(u,v))⋅(ru⃗×rv⃗) du dv.\iint_S \vec{F} \cdot dS = \iint_D \vec{F}(\vec{r}(u, v)) \cdot (\vec{r_u} \times \vec{r_v}) \, du \, dv.

We call this the flux of F⃗\vec{F} through SS.

Example 1

Given, F⃗=⟨z,y,x⟩\vec{F} = \langle z, y, x \rangle, find the flux of this vector field across the unit sphere.

So, we have the surface:

S={(x,y,z)∣x2+y2+z2=1}S = \{(x, y, z) \mid 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.

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

Now, one can solve this integral by taking the dot product and applying our definition:

∫02π∫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

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 \}. Is going to be a negative integral.

Definition 2 (Divergence)

For a vector field, F⃗:R3→R3\vec{F}: \mathbb{R}^3 \rarr \mathbb{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}.
Theorem 1 (The Divergence Theorem)

Let SS be a hollow surface and let EE be the 3-dimensional region it surrounds (∂E=S\partial E = S).

Then:

∬SF⃗⋅dS=∭Ediv⁡F⃗ dV.\iint_S \vec{F} \cdot dS = \iiint_E \operatorname{div} \vec{F} \, dV.

Equivalently:

∬∂EF⃗⋅dS=∭Ediv⁡F⃗ dV.\iint_{\partial E} \vec{F} \cdot dS = \iiint_E \operatorname{div} \vec{F} \, dV.

Now let’s take the same example.

Example 2

Find the flux of F⃗=⟨z,y,x⟩\vec{F} = \langle z, y, x \rangle 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∬SF⃗⋅dS=∭E1 dV\iint_S \vec{F} \cdot dS = \iiint_E 1 \, dV

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}}

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 R2\mathbb{R}^2?

Well, Stokes’ Theorem is essentially Green’s Theorem but in R3\mathbb{R}^3.

Definition 3 (Curl)

For F⃗:R3→R3\vec{F}: \mathbb{R}^3 \rarr \mathbb{R}^3, we define its curl as:

NB: it is also called rot and sometimes denoted as ∇×F⃗\nabla \times \vec{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.
Theorem 2 (Stokes’ Theorem)

If S⊂R3S \subset \mathbb{R}^3 is a surface, surrounded by a simple, closed curve in R3\mathbb{R}^3. Then:

∫CF⃗⋅dr=∬Scurl⁡F⃗⋅dS.\int_C \vec{F} \cdot dr = \iint_S \operatorname{curl} \vec{F} \cdot dS.
Example 3

Compute ∫CF⋅dr\int_C F \cdot dr for F⃗=⟨−y2,x,z2⟩\vec{F} = \langle -y^2, x, z^2 \rangle and CC being the intersection of y+z=2y + z = 2 with x2+y2=1x^2 + y^2 = 1.

Let’s find the curl:

curl⁡F⃗=⟨0,0,1+2y⟩\operatorname{curl} \vec{F} = \langle 0, 0, 1 + 2y \rangle

We need to find the surface, the intersection of y+z=2y + z = 2 and x2+y2=1x^2 + y^2 = 1 will just be the unit disk but at a slight tilt.

∬Scurl⁡F⃗⋅dS\iint_S \operatorname{curl} \vec{F} \cdot dS

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 \}

Meaning:

∬Scurl⁡F⃗⋅dS=∬Dcurl⁡F⃗⋅(ru⃗×rv⃗) dA\iint_S \operatorname{curl} \vec{F} \cdot dS = \iint_D \operatorname{curl} \vec{F} \cdot (\vec{r_u} \times \vec{r_v}) \, dA

Let’s find a parameterization of the surface:

r⃗(u,v)=⟨u,v,2−v⟩,ru⃗(u,v)=⟨1,0,0⟩,rv⃗(u,v)=⟨0,1,−1⟩,ru⃗×rv⃗=⟨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*}∬Dcurl⁡F⃗⋅(ru⃗×rv⃗) dA=∬D⟨0,0,1+2y⟩⋅⟨0,1,1⟩ dA=∬D1+2y dA=∫02π∫01(1+2rsin⁡(θ))r dr dθ=∫02π∫01r+2r2sin⁡(θ) dr dθ=∫01∫02πr+2r2sin⁡(θ) dθ dr\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*}

sin⁡(θ)\sin(\theta) has symmetry around 2π2\pi, therefore the integral of that is 00.

∫012πr dr=π\int_0^1 2\pi r \, dr = \boxed{\pi}
Example 4

Compute ∬Scurl⁡F⃗⋅dS\iint_S \operatorname{curl} \vec{F} \cdot dS where F⃗=⟨xz,yz,xy⟩\vec{F} = \langle xz, yz, xy \rangle and SS is the part of the sphere of radius 2 around (0,0,0)(0, 0, 0) that lies inside the cylinder x2+y2≤1x^2 + y^2 \leq 1 and above the xyxy-plane.

Let’s actually use Stokes’ theorem, “backwards”. CC will be a circle with:

{x2+y2=1x2+y2+z2=4\begin{cases} x^2 + y^2 = 1 \newline x^2 + y^2 + z^2 = 4 \end{cases}

Which means that z=3z = \sqrt{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*}∬Scurl⁡F⃗⋅dS=∫CF⃗⋅dr=∫02π⟨3cos⁡(t),3sin⁡(t),cos⁡(t)sin⁡(t)⟩⋅⟨−sin⁡(t),cos⁡(t),0⟩ dt=∫02π−3cos⁡(t)sin⁡(t)+3sin⁡(t)cos⁡(t) dt=∫02π0 dt=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*}

Search

Essays, lecture notes, and research

Type to search the site.