Chapter 9: Vector Calculus
Section 9.9: Stokes' Theorem
Example 9.9.3
Apply Stokes' theorem to the vector field F=y z i+x2z j+x y k; the curve C, the unit circle with center at the origin; and a closed cylinder with lid in the plane z=a>0 as the capping surface S.
Solution
Mathematical Solution
Begin by calculating
∇×F=ijk∂x∂y∂zy zx2zx y = x−x20z (2 x−1)
The normal component of this vector must be integrated over two surfaces: S1, the cylinder; and S2, the unit disk in the plane z=a. Using the parametrization x=cosθ,y=sinθ for the cylinder, an outward unit normal is N1=cosθ i+sinθ j+0 k. Hence, on this surface
∇×F·N1=cosθ−cos2θ cosθ
Represent the cylinder by the position vector R=cosθ i+sinθ j+z k so dσ=Rθ×Rz dθ dz. Now
Rθ×Rz = ijk−sinθcosθ0001 = cos(θ)sin(θ)0
so dσ=Rθ×Rz dθ dz = 1⋅dθ dz. Hence,
∫∫S1∇×F·N1 dσ = ∫02 π∫0acosθcosθ−cos2θdz dθ = π a
On the surface S2, N2=k, so ∇×F·N2=z 2 x−1. Using cylindrical coordinates,
∫01∫02 πa 2 r cosθ−1 r ⅆθ ⅆr= −π a
The net flux of ∇×F through S=S1⋃S2 is π a−π a=0.
The line integral can be evaluated if C is parametrized by the position vector
r=cost i+sint j+0 k so that on C
F·dr=00sin(t)cos(t)·−sin(t) dtcos(t) dt0= 0 dt
Clearly, then, ∫CF·dr = ∫02 π0 ⅆt = 0.
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
y z,x2z,x y = →to Vector Field →assign to a nameF
Obtain ∇×F
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
To evaluate ∫∫S1∇×F·N dσ, where S1 is the vertical wall of the cylinder, use the task template in Table 9.9.3(a) where ∇×F is given in Cartesian coordinates, but the cylinder S1 is given in cylindrical coordinates.
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 Parametric Surface
Flux through a Parametrically Defined Surface
For the Vector Field:
Select Coordinate SystemCartesian [x,y,z]Cartesian - othercylindricalsphericalbipolarcylindricalbisphericalcardioidalcardioidcylindricalcasscylindricalconicalellcylindricalhypercylindricalinvcasscylindricallogcylindricallogcoshcylindricaloblatespheroidalparaboloidalparacylindricalprolatespheroidalrosecylindricalsixspheretangentcylindricaltangentspheretoroidal
Select Coordinate SystemCartesiancylindricalsphericalbipolarcylindricalbisphericalcardioidalcardioidcylindricalcasscylindricalconicalellcylindricalhypercylindricalinvcasscylindricallogcylindricallogcoshcylindricaloblatespheroidalparaboloidalparacylindricalprolatespheroidalrosecylindricalsixspheretangentcylindricaltangentspheretoroidal
Table 9.9.3(a) ∫∫S1∇×F·N dσ, where S1 is the cylinder
To evaluate ∫∫S2∇×F·N dσ, where S2 is the "lid" for the cylinder, use the task template in Table 9.9.3(b) where ∇×F is given in Cartesian coordinates, but the cylinder S1 is given in cylindrical coordinates.
Tools≻Tasks≻Browse: Calculus - Vector≻Integration≻Flux≻3-D≻Through a Surface Defined over a Disk
Flux through a Surface Defined over a Disk
Table 9.9.3(b) ∫∫S2∇×F·N dσ, where S2 is the "lid" for the cylinder
Table 9.9.3(c) accesses the LineInt command through the Context Panel.
Form and evaluate the line integral of F around C
Write the name F. Context Panel: Evaluate and Display Inline
Context Panel: Student Vector Calculus≻Line Integral (Complete the dialog as per Figure 9.9.3(a).)
Context Panel: Evaluate Integral
Figure 9.9.3(a) Line Integral Domain dialog
F = →line integral∫02⁢π∫010ⅆrⅆt=0
Table 9.9.3(c) Evaluation of the line integral of F around 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
Maple Solution - Coded
Let S be the union of two surfaces S1 (the cylinder) and S2 (the disk on top).
Install the Student VectorCalculus package.
withStudent:-VectorCalculus:
Set display of vectors with BasisFormat command.
BasisFormatfalse:
Define F with the VectorField command.
F≔VectorFieldy z,x2z,x y:
Use the Curl and Flux commands to obtain the flux of ∇×F through S1
FluxCurlF,Surface1,θ,z,θ=0..2 π,z=0..a,coords=cylindricalr,θ,z,output=integral
∫02⁢π∫0a−cos⁡θ2+cos⁡θ⁢cos⁡θⅆzⅆθ
FluxCurlF,Surface1,θ,z,θ=0..2 π,z=0..a,coords=cylindricalr,θ,z = π⁢a
Use the Curl and Flux commands to obtain the flux of ∇×F through S2
FluxCurlF,Surfacer,θ,a,r=0..1,θ=0..2 π,coords=cylindricalr,θ,z,output=integral
∫01∫02⁢π2⁢r⁢cos⁡θ⁢a−a⁢cos⁡θ2⁢r+sin⁡θ2⁢rⅆθⅆr
FluxCurlF,Surfacer,θ,a,r=0..1,θ=0..2 π,coords=cylindricalr,θ,z = −π⁢a
Table 9.9.3(d) uses the LineInt command to obtain the value of ∳CF·dr, where C is the unit circle centered at the origin, and lying in the plane z=0.
LineIntF,Circle3D0,0,0,1,0,0,1,output=integral
∫02⁢π∫010ⅆrⅆt
LineIntF,Circle3D0,0,0,1,0,0,1 = 0
Table 9.9.3(d) Line integral of the tangential component of F around 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 has already been done twice, in the previous two sections.
<< Previous Example 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