# 16.8: Stokes' Theorem - Mathematics We are searching data for your request:

Forums and discussions:
Manuals and reference books:
Data from registers:
Wait the end of the search in all databases.
Upon completion, a link will appear to access the found materials.

So far the only types of line integrals which we have discussed are those along curves in (mathbb{R}^ 2). But the definitions and properties which were covered in Sections 4.1 and 4.2 can easily be extended to include functions of three variables, so that we can now discuss line integrals along curves in (mathbb{R}^ 3).

Definition (PageIndex{1}): Line Integrals

For a real-valued function (f (x, y, z)) and a curve (C) in (mathbb{R}^ 3), parametrized by (x = x(t), y = y(t), z = z(t), a ≤ t ≤ b), the line integral of (f (x, y, z) extbf{ along }C) with respect to arc length (s) is

[int_C f (x, y, z),ds = int_a^b f (x(t), y(t), z(t))sqrt{x ′ (t)^2 + y ′ (t)^2 + z ′ (t)^2},dt.label{Eq4.34}]

The line integral of (f (x, y, z)) along (C) with respect to (x) is

[int_C f (x, y, z),dx =int_a^b f (x(t), y(t), z(t)) x ′ (t),dt.label{Eq4.35}]

The line integral of (f (x, y, z)) along (C) with respect to (y) is

[int_C f (x, y, z),d y =int_a^b f (x(t), y(t), z(t)) y ′ (t),dt .label{Eq4.36}]

The line integral of (f (x, y, z)) along (C) with respect to (z) is

[int_C f (x, y, z),dz =int_a^b f (x(t), y(t), z(t)) z ′ (t),dt .label{Eq4.37}]

Similar to the two-variable case, if (f (x, y, z) ≥ 0) then the line integral (int_C f (x, y, z),ds) can be thought of as the total area of the “picket fence” of height (f (x, y, z)) at each point along the curve (C) in (mathbb{R}^ 3).

Vector fields in (mathbb{R}^ 3) are defined in a similar fashion to those in (mathbb{R}^ 2), which allows us to define the line integral of a vector field along a curve in (mathbb{R}^ 3).

Definition (PageIndex{2})

For a vector field ( extbf{f}(x, y, z) = P(x, y, z) extbf{i}+Q(x, y, z) extbf{j}+R(x, y, z) extbf{k}) and a curve (C) in (mathbb{R}^ 3) with a smooth parametrization (x = x(t), y = y(t), z = z(t), a ≤ t ≤ b), the line integral of (f) along (C) is

[egin{align} int_C extbf{f} cdot d extbf{r} &=int_C P(x, y, z),dx +int_C Q(x, y, z),d y + int_C R(x, y, z),dz label{Eq4.38} [4pt] &=int_a^b extbf{f}(x(t), y(t), z(t))cdot extbf{r} ′ (t)dt ,label{Eq4.39} [4pt] end{align}]

where ( extbf{r}(t) = x(t) extbf{i}+ y(t) extbf{j}+ z(t) extbf{k}) is the position vector for points on (C).

Similar to the two-variable case, if ( extbf{f}(x, y, z)) represents the force applied to an object at a point ((x, y, z)) then the line integral (int_C extbf{f}cdot d extbf{r}) represents the work done by that force in moving the object along the curve (C) in (mathbb{R}^ 3).

Some of the most important results we will need for line integrals in (mathbb{R}^ 3) are stated below without proof (the proofs are similar to their two-variable equivalents).

Theorem (PageIndex{1})

For a vector field ( extbf{f}(x, y, z) = P(x, y, z) extbf{i} +Q(x, y, z) extbf{j} + R(x, y, z) extbf{k}) and a curve (C) with a smooth parametrization (x = x(t), y = y(t), z = z(t), a ≤ t ≤ b) and position vector ( extbf{r}(t) = x(t) extbf{i}+ y(t) extbf{j}+ z(t) extbf{k}),

[int_C extbf{f}cdot d extbf{r} = int_C extbf{f}cdot extbf{T},ds,label{Eq4.40}]

where ( extbf{T}(t) = dfrac{ extbf{r} ′ (t)}{ lVert extbf{r} ′ (t) Vert }) is the unit tangent vector to (C) at ((x(t), y(t), z(t))).

Theorem (PageIndex{2}): Chain Rule

If (w = f (x, y, z)) is a continuously differentiable function of (x, y, ext{ and }z, ext{ and }x = x(t), y = y(t) ext{ and }z = z(t)) are differentiable functions of (t, ext{ then }w) is a differentiable function of (t), and

[dfrac{dw}{ dt} = dfrac{∂w}{ ∂x}dfrac{ dx}{ dt} + dfrac{∂w}{ ∂y}dfrac{ d y}{ dt} + dfrac{∂w }{∂z}dfrac{ dz}{ dt}.label{Eq4.41}]

Also, if (x = x(t_1 ,t_2), y = y(t_1 ,t_2) ext{ and }z = z(t_1 ,t_2)) are continuously differentiable function of ((t_1 ,t_2)), then

[dfrac{∂w}{ ∂t_1} = dfrac{∂w}{ ∂x}dfrac{ ∂x}{ ∂t_1} + dfrac{∂w}{ ∂y}dfrac{ ∂y}{ ∂t_1} +dfrac{∂w}{ ∂z}dfrac{ ∂z}{ ∂t_1} label{Eq4.42}]

and

[dfrac{∂w}{ ∂t_2} = dfrac{∂w}{ ∂x}dfrac{ ∂x}{ ∂t_2} + dfrac{∂w}{ ∂y}dfrac{ ∂y}{ ∂t_2} +dfrac{∂w}{ ∂z}dfrac{ ∂z}{ ∂t_2} label{Eq4.43}]

Theorem (PageIndex{3}): Potential

Let ( extbf{f}(x, y, z) = P(x, y, z) extbf{i} + Q(x, y, z) extbf{j} + R(x, y, z) extbf{k}) be a vector field in some solid (S), with (P, Q ext{ and }R) continuously differentiable functions on (S). Let (C) be a smooth curve in (S) parametrized by (x = x(t), y = y(t), z = z(t), a ≤ t ≤ b). Suppose that there is a real-valued function (F(x, y, z)) such that (∇F = extbf{f} ext{ on }S). Then

[int_C extbf{f}cdot d extbf{r} = F(B) − F(A) ,label{Eq4.44}]

where (A = (x(a), y(a), z(a)) ext{ and }B = (x(b), y(b), z(b))) are the endpoints of (C).

Corollary

If a vector field ( extbf{f}) has a potential in a solid (S), then (oint_C extbf{f}cdot d extbf{r} = 0) for any closed curve (C) in (S) (i.e. (oint_C ∇Fcdot d extbf{r} = 0) for any real-valued function (F(x, y, z))).

Example (PageIndex{1})

Let (f (x, y, z) = z) and let (C) be the curve in (mathbb{R}^ 3) parametrized by

[ onumber x = tsin t ,quad y = t cos t ,quad z = t ,quad 0 ≤ t ≤ 8π .]

Evaluate (int_C f (x, y, z),ds). (Note: (C) is called a conical helix. See Figure 4.5.1).

Solution

Since (x ′ (t) = sin t+ t cos t, y ′ (t) = cos t− tsin t, ext{ and }z ′ (t) = 1), we have

[ onumber egin{align}x ′ (t)^ 2 + y ′ (t)^ 2 + z ′ (t)^ 2 &= (sin^2 t+2tsin t cos t+ t^ 2 cos^2 t)+(cos^2 t−2tsin t cos t+ t^ 2 sin^2 t)+1 [4pt] onumber &=t^ 2 (sin^2 t+cos^2 t)+sin^2 t+cos^2 t+1 [4pt] onumber &=t^2 +2, [4pt] end{align}]

so since (f (x(t), y(t), z(t)) = z(t) = t) along the curve (C), then

[ onumber egin{align} int_C f (x, y, z),ds &= int_0^{8pi} f (x(t), y(t), z(t))sqrt{x ′ (t)^ 2 + y ′ (t)^ 2 + z ′ (t)^ 2},dt [4pt] onumber &=int_0^{8pi} tsqrt{t^2+2},dt [4pt] onumber &= left ( dfrac{1}{3} (t^ 2 +2)^{3/2} ight ) Big |_0^{8pi} = dfrac{1}{3} left ( (64π^ 2 +2)^{3/2} −2 sqrt{ 2} ight ). [4pt] end{align}] Example (PageIndex{2})

Let ( extbf{f}(x, y, z) = x extbf{i} + y extbf{j}+ 2z extbf{k}) be a vector field in (mathbb{R}^ 3). Using the same curve (C) from Example 4.12, evaluate (int_C extbf{f}cdot d extbf{r}).

Solution:

It is easy to see that (F(x, y, z) = dfrac{x^ 2}{ 2} + dfrac{y^ 2}{ 2} + z^ 2) is a potential for ( extbf{f}(x, y, z)) (i.e. (∇F = extbf{f})).

So by Theorem 4.12 we know that

[ onumber egin{align} int_C extbf{f}cdot d extbf{r} &=F(B) − F(A) , ext{ where }A = (x(0), y(0), z(0)) ext{ and }B = (x(8π), y(8π), z(8π)), ext{ so} [4pt] onumber &=F(8πsin 8π,8πcos 8π,8π) − F(0sin 0,0cos 0,0) [4pt] onumber &=F(0,8π,8π) − F(0,0,0) [4pt] onumber &= 0+ dfrac{(8π)^ 2}{ 2} +(8π)^ 2 −(0+0+0) = 96π^ 2 . [4pt] end{align}]

We will now discuss a generalization of Green’s Theorem in (mathbb{R}^ 2) to orientable surfaces in (mathbb{R}^ 3), called Stokes’ Theorem. A surface (Σ) in (mathbb{R}^ 3) is orientable if there is a continuous vector field N in (mathbb{R}^ 3) such that N is nonzero and normal to (Σ) (i.e. perpendicular to the tangent plane) at each point of (Σ). We say that such an N is a normal vector field.

For example, the unit sphere (x^ 2+y^ 2+z^ 2 = 1) is orientable, since the continuous vector field ( extbf{N}(x, y, z) = x extbf{i}+ y extbf{j}+z extbf{k}) is nonzero and normal to the sphere at each point. In fact, (− extbf{N}(x, y, z)) is another normal vector field (see Figure 4.5.2). We see in this case that ( extbf{N}(x, y, z)) is what we have called an outward normal vector, and (− extbf{N}(x, y, z)) is an inward normal vector. These “outward” and “inward” normal vector fields on the sphere correspond to an “outer” and “inner” side, respectively, of the sphere. That is, we say that the sphere is a two-sided surface. Roughly, “two-sided” means “orientable”. Other examples of two-sided, and hence orientable, surfaces are cylinders, paraboloids, ellipsoids, and planes. You may be wondering what kind of surface would not have two sides. An example is the Möbius strip, which is constructed by taking a thin rectangle and connecting its ends at the opposite corners, resulting in a “twisted” strip (see Figure 4.5.3). If you imagine walking along a line down the center of the Möbius strip, as in Figure 4.5.3(b), then you arrive back at the same place from which you started but upside down! That is, your orientation changed even though your motion was continuous along that center line. Informally, thinking of your vertical direction as a normal vector field along the strip, there is a discontinuity at your starting point (and, in fact, at every point) since your vertical direction takes two different values there. The Möbius strip has only one side, and hence is nonorientable.

For an orientable surface (Σ) which has a boundary curve (C), pick a unit normal vector n such that if you walked along (C) with your head pointing in the direction of n, then the surface would be on your left. We say in this situation that n is a positive unit normal vector and that (C) is traversed n-positively. We can now state Stokes’ Theorem:

Theorem (PageIndex{4}): Stoke's Theorem

Let (Σ) be an orientable surface in (mathbb{R}^ 3) whose boundary is a simple closed curve (C), and let ( extbf{f}(x, y, z) = P(x, y, z) extbf{i} +Q(x, y, z) extbf{j} + R(x, y, z) extbf{k}) be a smooth vector field defined on some subset of (mathbb{R}^ 3) that contains (Σ). Then

[oint_C extbf{f}cdot d extbf{r} = iintlimits_Σ ( ext{curl } extbf{f})cdot extbf{n},dσ , label{Eq4.45}]

where

[ ext{curl } extbf{f} = left ( dfrac{∂R}{ ∂y}-dfrac{∂Q}{ ∂z} ight ) extbf{i} + left ( dfrac{∂P}{ ∂z} - dfrac{∂R}{ ∂x} ight ) extbf{j} + left ( dfrac{∂Q}{ ∂x} -dfrac{∂P}{ ∂y} ight ) extbf{k},label{Eq4.46}]

n is a positive unit normal vector over (Σ), and (C) is traversed n-positively.

Proof: As the general case is beyond the scope of this text, we will prove the theorem only for the special case where (Σ) is the graph of (z = z(x, y)) for some smooth real-valued function (z(x, y), ext{ with }(x, y)) varying over a region (D) in (mathbb{R}^ 2).

Projecting (Σ) onto the (x y)-plane, we see that the closed curve (C) (the boundary curve of (Σ)) projects onto a closed curve (C_D) which is the boundary curve of (D) (see Figure 4.5.4). Assuming that (C) has a smooth parametrization, its projection (C_D) in the (x y)-plane also has a smooth parametrization, say [ onumber C_D, :, x = x(t) , y = y(t) , a ≤ t ≤ b ,]

and so (C) can be parametrized (in (mathbb{R}^ 3) ) as

[ onumber C ,:, x = x(t) ,, y = y(t) ,, z = z(x(t), y(t)) ,, a ≤ t ≤ b ,]

since the curve (C) is part of the surface (z = z(x, y)). Now, by the Chain Rule (Theorem 4.4 in Section 4.2), for (z = z(x(t), y(t)) ext{ as a function of }t), we know that

[ onumber z ′ (t) = dfrac{∂z}{ ∂x} x ′ (t) + dfrac{∂z}{ ∂y} y ′ (t) ,]

and so

[ onumber egin{align} oint_C extbf{f}cdot d extbf{r} &= int_C P(x, y, z),dx+Q(x, y, z),d y+ R(x, y, z),dz [4pt] onumber &=int_a^b left ( P x′ (t)+Q y′ (t)+ R left ( dfrac{∂z}{ ∂x} x ′ (t)+ dfrac{∂z}{ ∂y} y ′ (t) ight ) ight ) ,dt [4pt] onumber &=int_a^b left ( left ( P + R dfrac{∂z}{ ∂x} ight ) x ′ (t)+ left ( Q + R dfrac{∂z}{ ∂y} ight ) y ′ (t) ight ) ,dt [4pt] onumber &=int_{C_D} ilde P (x, y),dx+ ilde Q (x, y),d y , [4pt] end{align}]

where

[ onumber egin{align} & ilde P (x, y) = P(x, y, z(x, y)) + R(x, y, z(x, y)) dfrac{∂z}{ ∂x} (x, y) , ext{ and} [4pt] onumber & ilde Q (x, y) = Q(x, y, z(x, y)) + R(x, y, z(x, y)) dfrac{∂z}{ ∂y} (x, y) [4pt] end{align}]

for ((x, y) ext{ in }D). Thus, by Green’s Theorem applied to the region (D), we have

[oint_C extbf{f}cdot d extbf{r} = iintlimits_D left ( dfrac{∂ ilde Q}{ ∂x} − dfrac{∂ ilde P}{ ∂y} ight ) ,dA.label{Eq4.47}]

Thus,

[ onumber egin{align} dfrac{∂ ilde Q}{ ∂x} &= dfrac{∂}{ ∂x} left ( Q(x, y, z(x, y))+ R(x, y, z(x, y)) dfrac{∂z}{ ∂y} (x, y) ight ), ext{ so by the Product Rule we get} [4pt] onumber &= dfrac{∂}{ ∂x} (Q(x, y, z(x, y)))+left (dfrac{∂}{ ∂x} R(x, y, z(x, y)) ight ) dfrac{∂z}{ ∂y} (x, y)+R(x, y, z(x, y)) dfrac{∂}{ ∂x} left ( dfrac{∂z}{ ∂y} (x, y) ight ) [4pt] end{align}]

Now, by Equation ef{Eq4.42} in Theorem 4.11, we have

[ onumber egin{align} dfrac{∂}{ ∂x} (Q(x, y, z(x, y)))&= dfrac{∂Q}{ ∂x}dfrac{ ∂x}{ ∂x}+dfrac{∂Q}{ ∂y}dfrac{ ∂y}{ ∂x}+dfrac{∂Q}{ ∂z}dfrac{ ∂z}{ ∂x} [4pt] onumber &= dfrac{∂Q}{ ∂x} cdot 1 + dfrac{∂Q}{ ∂y} cdot 0 + dfrac{∂Q}{ ∂z}dfrac{ ∂z}{ ∂x} [4pt] onumber &= dfrac{∂Q}{ ∂x} + dfrac{∂Q}{∂z}dfrac{ ∂z}{ ∂x} . [4pt] end{align}]

Similarly,

[ onumber dfrac{∂}{ ∂x} (R(x, y, z(x, y))) = dfrac{∂R}{ ∂x} + dfrac{∂R}{ ∂z}dfrac{ ∂z}{ ∂x} .]

Thus,

[ onumber egin{align}dfrac{∂ ilde Q}{ ∂x} &=dfrac{∂Q}{ ∂x} + dfrac{∂Q}{ ∂z}dfrac{ ∂z}{ ∂x} + left ( dfrac{∂R}{ ∂x} + dfrac{∂R}{ ∂z}dfrac{ ∂z}{ ∂x} ight ) dfrac{∂z }{∂y}+R(x, y, z(x, y)) dfrac{∂^2 z}{ ∂x∂y} [4pt] onumber &=dfrac{∂Q}{ ∂x} + dfrac{∂Q}{ ∂z}dfrac{ ∂z}{ ∂x} + dfrac{∂R}{ ∂x}dfrac{ ∂z}{ ∂y}+dfrac{∂R}{ ∂z}dfrac{ ∂z}{ ∂x}dfrac{ ∂z}{ ∂y} + R dfrac{∂^2 z}{ ∂x∂y} . [4pt] end{align}]

In a similar fashion, we can calculate

[ onumber dfrac{∂ ilde P}{ ∂y} = dfrac{∂P}{ ∂y} + dfrac{∂P}{ ∂z}dfrac{ ∂z}{ ∂y} + dfrac{∂R}{ ∂y}dfrac{ ∂z}{ ∂x} + dfrac{∂R}{ ∂z}dfrac{ ∂z}{ ∂y}dfrac{ ∂z}{ ∂x} +Rdfrac{∂^2 z}{ ∂y∂x} .]

So subtracting gives

[dfrac{∂ ilde Q}{ ∂x} - dfrac{∂ ilde P}{ ∂y} = left ( dfrac{∂Q}{ ∂z} − dfrac{∂R}{ ∂y} ight ) dfrac{∂z}{ ∂x} + left ( dfrac{∂R}{ ∂x} − dfrac{∂P}{ ∂z} ight ) dfrac{∂z}{ ∂y} + left ( dfrac{∂Q}{ ∂x} − dfrac{∂P}{ ∂y} ight ) label{Eq4.48}]

since (dfrac{∂^2 z}{ ∂x∂y} = dfrac{∂^2 z}{ ∂y∂x}) by the smoothness of (z = z(x, y)). Hence, by Equation ef{Eq4.47},

[oint_C extbf{f}cdot d extbf{r} = iint_D left ( - left (dfrac{∂R}{ ∂y} − dfrac{∂Q}{ ∂z} ight ) dfrac{∂z}{ ∂x} - left ( dfrac{∂P}{ ∂z} − dfrac{∂R}{ ∂x} ight ) dfrac{∂z}{ ∂y} + left ( dfrac{∂Q}{ ∂x} − dfrac{∂P}{ ∂y} ight ) ight ) ,dA label{Eq4.49}]

after factoring out a −1 from the terms in the first two products in Equation ef{Eq4.48}.

Now, recall from Section 2.3 (see p.76) that the vector ( extbf{N} = − dfrac{∂z}{ ∂x} extbf{ i}− dfrac{∂z}{ ∂y} extbf{j}+ extbf{k}) is normal to the tangent plane to the surface (z = z(x, y)) at each point of (Σ). Thus,

[ onumber extbf{n} = dfrac{ extbf{N}}{lVert extbf{N} Vert } = dfrac{− dfrac{∂z}{ ∂x} extbf{ i}− dfrac{∂z}{ ∂y} extbf{j}+ extbf{k}}{sqrt{1+left ( dfrac{∂z}{ ∂x} ight ) ^2+left ( dfrac{∂z}{ ∂y} ight )^2}}]

is in fact a positive unit normal vector to (Σ) (see Figure 4.5.4). Hence, using the parametrization ( extbf{r}(x, y) = x extbf{i} + y extbf{j} + z(x, y) extbf{k}, ext{ for }(x, y) ext{ in }D), of the surface (Σ), we have (dfrac{∂ extbf{r}}{ ∂x} = extbf{i} + dfrac{∂z}{ ∂x} extbf{k}) and (dfrac{∂ extbf{r}}{ ∂y} = extbf{j}+ dfrac{∂z}{ ∂y} extbf{k}), and so (lVert dfrac{ ∂ extbf{r}}{ ∂x} imes dfrac{∂ extbf{r}}{ ∂y} Vert = sqrt{ 1+ left ( dfrac{∂z}{ ∂x} ight )^2 + left ( dfrac{∂z}{ ∂y} ight )^2}). So we see that using Equation ef{Eq4.46} for curl f, we have

[ onumber egin{align} iintlimits_Σ ( ext{curl } extbf{f})cdot extbf{n} ,dσ &=iintlimits_D ( ext{curl } extbf{f}) Big lVert dfrac{∂ extbf{r}}{ ∂x} imes dfrac{∂ extbf{r}}{ ∂y} Big Vert , dA [4pt] onumber &=iintlimits_D left ( left ( dfrac{∂R}{ ∂y} − dfrac{∂Q}{ ∂z} ight ) extbf{i} + left ( dfrac{∂P}{ ∂z} − dfrac{∂R}{ ∂x} ight ) extbf{j} + left ( dfrac{∂Q}{ ∂x} − dfrac{∂P}{ ∂y} ight ) extbf{k} ight ) cdot left ( -dfrac{∂z}{ ∂x} extbf{i}− dfrac{∂z}{ ∂y} extbf{j}+ extbf{k} ight ) ,dA [4pt] onumber &=iintlimits_D left ( - left ( dfrac{∂R}{ ∂y} − dfrac{∂Q}{ ∂z} ight ), dfrac{∂z}{ ∂x} - left ( dfrac{∂P}{ ∂z} − dfrac{∂R}{ ∂x} ight ) , dfrac{∂z}{ ∂y} + left ( dfrac{∂Q}{ ∂x} − dfrac{∂P}{ ∂y} ight ) ight ) , dA, [4pt] end{align}]

which, upon comparing to Equation ef{Eq4.49}, proves the Theorem.

( ag{( extbf{QED})})

Note: The condition in Stokes’ Theorem that the surface (Σ) have a (continuously varying) positive unit normal vector n and a boundary curve (C) traversed n-positively can be expressed more precisely as follows: if ( extbf{r}(t)) is the position vector for (C) and ( extbf{T}(t) = extbf{r} ′ (t)/ Vert extbf{r} ′ (t) Vert) is the unit tangent vector to (C), then the vectors T, n, T ( imes) n form a right-handed system.

Also, it should be noted that Stokes’ Theorem holds even when the boundary curve (C) is piecewise smooth.

Example (PageIndex{3})

Verify Stokes’ Theorem for ( extbf{f}(x, y, z) = z extbf{i} + x extbf{j} + y extbf{k}) when (Σ) is the paraboloid (z = x^ 2 + y^ 2) such that (z ≤ 1) (see Figure 4.5.5). Solution:

The positive unit normal vector to the surface (z = z(x, y) = x^ 2 + y^ 2) is

[ onumber extbf{n} = dfrac{−dfrac{∂z}{ ∂x} extbf{i}− dfrac{∂z}{ ∂y} extbf{j}+ extbf{k}}{sqrt{1+ left ( dfrac{∂z}{ ∂x} ight )^2+left ( dfrac{∂z}{ ∂y} ight )^2}}= dfrac{−2x extbf{i}−2y extbf{j}+ extbf{k}}{sqrt{1+4x^2 +4y^2}},]

and curl f = (1−0)i+(1−0)j+(1−0)k = i+j+k, so

[ onumber ( ext{curl } extbf{f})cdot extbf{n} = (−2x−2y+1)/sqrt{1+4x^ 2 +4y^ 2} .]

Since (Σ) can be parametrized as ( extbf{r}(x, y) = x extbf{i} + y extbf{j} + (x^ 2 + y^ 2 ) extbf{k} ext{ for }(x, y) ext{ in the region }D = {(x, y) : x^ 2 + y^ 2 ≤ 1}), then

[ onumber egin{align} iintlimits_Σ ( ext{curl } extbf{f})cdot extbf{n},dσ &=iintlimits_D ( ext{curl } extbf{f}) Big lVert dfrac{∂ extbf{r}}{ ∂x} imes dfrac{∂ extbf{r}}{ ∂y} Big Vert , dA [4pt] onumber &= iintlimits_D dfrac{−2x−2y+1}{sqrt{1+4x^ 2 +4y^ 2}}sqrt{1+4x^ 2 +4y^ 2}, d A [4pt] onumber &= iintlimits_D (−2x−2y+1),d A, ext{ so switching to polar coordinates gives} [4pt] onumber &=int_0^{2pi} int_0^1 (−2r cos θ −2r sin θ +1),r ,dr, dθ [4pt] onumber &=int_0^{2pi} int_0^1 (−2r^ 2 cos θ −2r^ 2 sin θ + r),dr, dθ [4pt] onumber &=int_0^{2pi} left (− dfrac{2r^ 3}{ 3} cos θ − dfrac{2r^ 3}{ 3} sin θ + dfrac{r^ 2}{ 2} Big |_{r=0}^{r=1} ight ) ,dθ [4pt] onumber &= int_0^{2pi} left ( − dfrac{2}{ 3} cos θ − dfrac{2}{ 3} sin θ + dfrac{1}{ 2} ight ) , dθ [4pt] onumber &= − dfrac{2}{ 3} sin θ +dfrac{2}{ 3} cos θ + dfrac{1}{ 2} θ Big |_0^{2pi} = pi . [4pt] end{align}]

The boundary curve (C) is the unit circle (x^ 2 + y^ 2 = 1) laying in the plane (z = 1) (see Figure 4.5.5), which can be parametrized as (x = cos t, y = sin t, z = 1 ext{ for }0 ≤ t ≤ 2π). So

[ onumber egin{align} oint_C extbf{f}cdot d extbf{r} &= int_0^{2pi} ((1)(−sin t)+(cos t)(cos t)+(sin t)(0)),dt [4pt] onumber &= int_0^{2pi} left ( −sin t+ dfrac{1+cos 2t}{ 2} ight ) ,dt quad left ( ext{here we used } cos^2 t = dfrac{1+cos 2t}{ 2} ight ) [4pt] onumber &= cos t+ dfrac{t}{ 2} + dfrac{sin 2t}{ 4}Big |_0^{2pi} = π . [4pt] end{align}]

So we see that (oint_C extbf{f}cdot d extbf{r} =iintlimits_Σ ( ext{curl } extbf{f})cdot extbf{n}dσ), as predicted by Stokes’ Theorem.

The line integral in the preceding example was far simpler to calculate than the surface integral, but this will not always be the case.

Example (PageIndex{4})

Let (Σ) be the elliptic paraboloid (z =dfrac{x^ 2}{ 4} + dfrac{y^ 2}{ 9} ext{ for }z ≤ 1), and let (C) be its boundary curve. Calculate ( oint_C extbf{f}cdot d extbf{r} ext{ for } extbf{f}(x, y, z) = (9xz + 2y) extbf{i} + (2x + y^ 2 ) extbf{j} + (−2y^ 2 + 2z) extbf{k}), where (C) is traversed counterclockwise

Solution

The surface is similar to the one in Example (PageIndex{3}), except now the boundary curve (C) is the ellipse (dfrac{x^ 2}{ 4} + dfrac{y^ 2}{ 9} = 1) laying in the plane (z = 1). In this case, using Stokes’ Theorem is easier than computing the line integral directly. As in Example 4.14, at each point ((x, y, z(x, y))) on the surface (z = z(x, y) = dfrac{x^ 2}{ 4} + dfrac{y^ 2}{ 9}) the vector

[ onumber extbf{n} = dfrac{-dfrac{∂z}{ ∂x} extbf{i}− dfrac{∂z}{ ∂y} extbf{j}+ extbf{k}}{sqrt{1 + left ( dfrac{∂z}{ ∂x} ight )^2 + left ( dfrac{∂z}{ ∂y} ight )^2}}=dfrac{− dfrac{x}{ 2} extbf{i}− dfrac{2y}{ 9} extbf{j}+ extbf{k}}{sqrt{1+ dfrac{x^ 2}{ 4} + dfrac{4y^ 2}{ 9}}},]

is a positive unit normal vector to (Σ). And calculating the curl of f gives

[ onumber ext{curl } extbf{f} = (−4y−0) extbf{i} + (9x−0) extbf{j} + (2−2) extbf{k} = −4y extbf{i} + 9x extbf{j} + 0 extbf{k} ,]

so

[ onumber ( ext{curl } extbf{f})cdot extbf{n} = dfrac{(−4y)(− dfrac{x}{ 2} )+(9x)(− dfrac{2y}{ 9} )+(0)(1)}{sqrt{1+ dfrac{x^ 2}{ 4} + dfrac{4y^ 2}{ 9}}} = dfrac{2x y−2x y+0}{sqrt{1+ dfrac{x^ 2}{ 4} + dfrac{4y^ 2}{ 9}}} = 0,]

and so by Stokes’ Theorem

[ onumber oint_C extbf{f}cdot d extbf{r} = iintlimits_Σ ( ext{curl } extbf{f})cdot extbf{n},dσ = iintlimits_Σ 0, dσ = 0.]

In physical applications, for a simple closed curve (C) the line integral (oint_C extbf{f}cdot d extbf{r}) is often called the circulation of f around (C). For example, if E represents the electrostatic field due to a point charge, then it turns out that curl ( extbf{E}= extbf{0}), which means that the circulation (oint_C extbf{E}cdot d extbf{r} = 0) by Stokes’ Theorem. Vector fields which have zero curl are often called irrotational fields.

In fact, the term curl was created by the 19th century Scottish physicist James Clerk Maxwell in his study of electromagnetism, where it is used extensively. In physics, the curl is interpreted as a measure of circulation density. This is best seen by using another definition of curl f which is equivalent to the definition given by Equation ef{Eq4.46}. Namely, for a point ((x, y, z) ext{ in }mathbb{R}^ 3),

[ extbf{n}cdot ( ext{curl } extbf{f}(x, y, z) = limlimits_{S o 0} dfrac{1}{S}oint_C extbf{f}cdot d extbf{r},label{Eq4.50}]

where (S) is the surface area of a surface (Σ) containing the point ((x, y, z)) and with a simple closed boundary curve (C) and positive unit normal vector n at ((x, y, z)). In the limit, think of the curve (C) shrinking to the point ((x, y, z)), which causes (Σ), the surface it bounds, to have smaller and smaller surface area. That ratio of circulation to surface area in the limit is what makes the curl a rough measure of circulation density (i.e. circulation per unit area). An idea of how the curl of a vector field is related to rotation is shown in Figure 4.5.6. Suppose we have a vector field ( extbf{f}(x, y, z)) which is always parallel to the (x y)-plane at each point ((x, y, z)) and that the vectors grow larger the further the point ((x, y, z)) is from the (y)-axis. For example, ( extbf{f}(x, y, z) = (1+ x^ 2 ) extbf{j}). Think of the vector field as representing the flow of water, and imagine dropping two wheels with paddles into that water flow, as in Figure 4.5.6. Since the flow is stronger (i.e. the magnitude of f is larger) as you move away from the (y)-axis, then such a wheel would rotate counterclockwise if it were dropped to the right of the (y)-axis, and it would rotate clockwise if it were dropped to the left of the (y)-axis. In both cases the curl would be nonzero (curl ( extbf{f}(x, y, z) = 2x extbf{k}) in our example) and would obey the right-hand rule, that is, curl ( extbf{f}(x, y, z)) points in the direction of your thumb as you cup your right hand in the direction of the rotation of the wheel. So the curl points outward (in the positive (z)-direction) if (x > 0) and points inward (in the negative (z)-direction) if (x < 0). Notice that if all the vectors had the same direction and the same magnitude, then the wheels would not rotate and hence there would be no curl (which is why such fields are called irrotational, meaning no rotation).

Finally, by Stokes’ Theorem, we know that if (C) is a simple closed curve in some solid region (S) in (mathbb{R}^ 3) and if ( extbf{f}(x, y, z)) is a smooth vector field such that curl ( extbf{f} = 0 ext{ in }S), then

[ onumber oint_C extbf{f}cdot d extbf{r} = iintlimits_Σ ( ext{curl } extbf{f}cdot extbf{n},dσ = iintlimits_Σ extbf{0}cdot extbf{n}, dσ = iintlimits_Σ 0,dσ = 0,]

where (Σ) is any orientable surface inside (S) whose boundary is (C) (such a surface is sometimes called a capping surface for (C)). So similar to the two-variable case, we have a threedimensional version of a result from Section 4.3, for solid regions in (mathbb{R}^ 3) which are simply connected (i.e. regions having no holes):

The following statements are equivalent for a simply connected solid region (S) in (mathbb{R}^ 3) :

1. ( extbf{f}(x, y, z) = P(x, y, z) extbf{i}+Q(x, y, z) extbf{j}+ R(x, y, z) extbf{k}) has a smooth potential (F(x, y, z) ext{ in }S)
2. (int_C extbf{f}cdot d extbf{r}) is independent of the path for any curve (C) in (S)
3. (oint_C extbf{f}cdot d extbf{r} = 0) for every simple closed curve (C) in (S)
4. (dfrac{∂R}{ ∂y} = dfrac{∂Q}{ ∂z} , dfrac{∂P}{ ∂z} = dfrac{∂R}{ ∂x} , ext{ and }dfrac{∂Q}{ ∂x} = dfrac{∂P}{ ∂y}) in (S) (i.e. curl ( extbf{f} = extbf{0} ext{ in }S))

Part (d) is also a way of saying that the differential form (P ,dx+Q ,d y+ R, dz) is exact.

Example (PageIndex{5})

Determine if the vector field ( extbf{f}(x, y, z) = x yz extbf{i}+xz extbf{j}+x y extbf{k}) has a potential in (mathbb{R}^ 3).

Solution

Since (mathbb{R}^ 3) is simply connected, we just need to check whether curl f = 0 throughout (mathbb{R}^ 3), that is,

[ onumber dfrac{∂R}{ ∂y} = dfrac{∂Q}{ ∂z} ,quad dfrac{∂P}{ ∂z} = dfrac{∂R}{ ∂x} ,quad ext{and }dfrac{∂Q}{ ∂x} = dfrac{∂P}{ ∂y} ]

throughout (mathbb{R}^ 3), where (P(x, y, z) = x yz, Q(x, y, z) = xz, ext{ and }R(x, y, z) = x y). But we see that

[ onumber dfrac{∂P}{ ∂z} = x y ,, dfrac{∂R}{ ∂x} = y quad Rightarrow quad dfrac{∂P}{ ∂z} eq dfrac{∂R}{ ∂x} ext{ for some }(x, y, z) ext{ in }mathbb{R}^ 3 .]

Thus, (f(x, y, z)) does not have a potential in (mathbb{R}^ 3).