Chapter 9: Vector Calculus
Section 9.9: Stokes' Theorem
Example 9.9.1
Apply Stokes' theorem to the vector field F=z i−x j−y k; the curve C, the unit circle with center at the origin; and the upper hemisphere as the capping surface S.
Solution
Mathematical Solution
Figure 9.9.1(a) shows arrows of the curl field of F, along with the upper hemisphere S. The black arrow at the "north pole" is a representative normal taken on S.
The curl of F:
∇×F=ijk∂x∂y∂zz−x−y = −11−1
A unit normal on S can be obtained by normalizing Rx×Ry , where R is a position-vector representation of S.
use Student:-VectorCalculus in module() local F,p; F:=VectorField(<z,-x,-y>); p:=Flux(Curl(F),Surface(<x,y,sqrt(-x^2-y^2+1)>,x=-1..1,y=-sqrt(-x^2+1)..sqrt(-x^2+1)),output=plot,fieldoptions=[grid=[3,3,3]],scaling=constrained,caption="",tickmarks=[3,3,3],axes=frame,orientation=[155,85,0]); print(p); end module: end use:
Figure 9.9.1(a) ∇×F and the surface S
Let Z=1−x2−y2 be a Cartesian representation of the upper hemisphere S. It follows that
R=xyZ ⇒Rx×Ry = ijk10−x/Z01−y/Z = x/Zy/Z1 and N=xyZ
The element of surface area can be obtained from
1+Zx2+Zy2=1+x/Z2+y/Z2 = 1/Z
or from Rx×Ry=1/Z. In either event,
∇×F·N dσ=−11−1·xyZ⋅1Z dA=yZ−xZ−1dA
so that ∫∫S∇×F·N dσ can be implemented in Cartesian coordinates as
∫−11∫−1−x21−x2y1−x2−y2−x1−x2−y2−1 dy dx=−π
or in polar coordinates as
∫01∫02 πrr sinθ1−r2−r cosθ1−r2−1 dθ dr= −π
The line integral around C, the unit circle centered at the origin, given by ∳CF·dr, can be evaluated if C is parametrized by the position vector r=cost i+sint j+0 k so that on C
F·dr=0−cos(t)−sin(t)·−sin(t) dtcos(t) dt0= −cos2t dt
Consequently, ∳CF·dr = −∫02 πcos2t dt= −π.
The parametrization chosen for C induces a counterclockwise traverse of the circle, an orientation consistent with the choice of an outward normal on S.
Maple Solution - Interactive
Initialize
Tools≻Load Package: Student Vector Calculus
Loading Student:-VectorCalculus
Tools≻Tasks≻Browse: Calculus - Vector≻ Vector Algebra and Settings≻ Display Format for Vectors
Press the Access Settings button and select "Display as Column Vector"
Display Format for Vectors
Define the vector field F
Enter a free vector whose components are those of F. Context Panel: Evaluate and Display Inline
Context Panel: Student Vector Calculus≻Conversions≻To Vector Field
Context Panel: Assign to a Name≻F
z,−x,−y = →to Vector Field →assign to a nameF
Obtain ∇×F on the surface S
Common Symbols palette: Del and cross-product operators
Context Panel: Evaluate and Display Inline
Context Panel: Assign to a Name≻curlF
∇×F = →assign to a namecurlF
Table 9.9.1(a) contains a task template with which the flux of ∇×F through S is computed. Should the "Clear All and Reset" button in the Task Template be pressed, all the data that has been input to the template will be lost. In that event, the reader should simply re-launch the example to recover the appropriate inputs to the template.
Tools≻Tasks≻Browse: Calculus - Vector≻Integration≻Flux≻3-D≻Through a Surface Defined over a Disk
Flux through a Surface Defined over a Disk
For the Vector Field:
Select Coordinate SystemCartesian [x,y,z]Cartesian - othercylindricalsphericalbipolarcylindricalbisphericalcardioidalcardioidcylindricalcasscylindricalconicalellcylindricalhypercylindricalinvcasscylindricallogcylindricallogcoshcylindricaloblatespheroidalparaboloidalparacylindricalprolatespheroidalrosecylindricalsixspheretangentcylindricaltangentspheretoroidal
Table 9.9.1(a) Flux of ∇×F through S
Table 9.9.1(b) contains the calculation of the line integral ∳CF·dr , where C is the circle capped by S.
Write the name F. Context Panel: Evaluate and Display Inline
Context Panel: Student Vector Calculus≻Line Integral (Complete the Line Integral Domain dialog as per Figure 9.9.1(b).)
Context Panel: Evaluate Integral
Figure 9.9.1(b) Line Integral Domain dialog
F = →line integral∫02⁢π∫01−rⅆrⅆt=−π
Table 9.9.1(b) Line integral: tangential component of F along the circle C
The astute reader will realize that Maple has evaluated a line integral as an iterated double-integral by invoking Stokes' theorem! Consequently, a validation of Stokes' theorem demands that the line integral be evaluated from first principles. This is easily done if C is parametrized by the position vector
r=cost i+sint j+0 k so that on C
and hence, ∫CF·dr = −∫02 πcos2t ⅆt = −π.
Maple Solution - Coded
Install the Student VectorCalculus package.
withStudent:-VectorCalculus:
Set display of vectors with BasisFormat command.
BasisFormatfalse:
Define F with the VectorField command.
F≔VectorFieldz,−x,−y:
Define the upper hemisphere as Z=zx,y.
Z≔1−x2−y2:
Use the Curl and Flux commands to obtain the flux of ∇×F through S
FluxCurlF,Surfacex,y,Z,x,y=Circle0,0,1,r,θ,output=integral
∫01∫02⁢π−r⁢r⁢cos⁡θ−r⁢sin⁡θ+−r2+1−r2+1ⅆθⅆr
FluxCurlF,Surfacex,y,Z,x,y=Circle0,0,1,r,θ = −π
Use the LineInt command to form and evaluate ∳CF·dr
LineIntF,Circle3D0,0,0,1,0,0,1,output=integral
∫02⁢π∫01−rⅆrⅆt
LineIntF,Circle3D0,0,0,1,0,0,1 = −π
The astute reader will realize that Maple has evaluated a line integral as an iterated double-integral by invoking Stokes' theorem! Consequently, a validation of Stokes' theorem demands that the line integral be evaluated from first principles. This has already been done twice, in the previous two sections.
Figure 9.9.1(a) can be obtained with the Flux command, provided the integration is implemented in Cartesian coordinates with the following syntax. The actual options applied in the figure can be seen in the code hidden in the table cell containing the figure.
FluxCurlF,Surfacex,y,Z,x=−1..1,y=−1−x2..1−x2,output=plot:
<< Chapter Overview Section 9.9 Next Example >>
© Maplesoft, a division of Waterloo Maple Inc., 2024. All rights reserved. This product is protected by copyright and distributed under licenses restricting its use, copying, distribution, and decompilation.
For more information on Maplesoft products and services, visit www.maplesoft.com
Download Help Document