Chapter 9: Vector Calculus
Section 9.8: Divergence Theorem
Example 9.8.3
Apply the Divergence theorem to the vector field F=x i+y j+z k and R, the rectangular box defined by x∈1,2,y∈3,4,z∈5,6.
Solution
Mathematical Solution
Figure 9.8.3(a) shows the box with outward normals on each face, and arrows of the field F.
The divergence of F:
∇·F=∂xx+∂yy+∂zz=1+1+1=3
The integral of ∇·F over the interior of the box:
∫12∫34∫563 dz dy dx = 3
use Student:-VectorCalculus in module() local F,p; F:=VectorField(<x,y,z>); p:=Flux(F,Box(1..2,3..4,5..6),output=plot,caption="",axes=frame,labels=[x,y,z],tickmarks=[2,2,2],orientation=[-35,75,0],fieldoptions=[grid=[5,5,5]]); print(p); end module: end use:
Figure 9.8.3(a) Box and field F
The flux must be computed through each of the six faces of the box. To begin, note that
Face
dσ
N
F·N
x=1
dy dz
−i
−1
x=2
i
2
y=3
dx dz
−j
−3
y=4
j
4
z=5
dx dy
−k
−5
z=6
k
6
so that, by pairing opposing faces, the following three integrals represent the flux through the faces of the box.
∫34∫562−1 dz dy+∫12∫564−3 dz dx+∫12∫346−5 dy dx = 3
Maple Solution - Interactive
The Student VectorCalculus package is needed for calculating the divergence, but it then conflicts with any multidimensional integral set from the Calculus palette. Hence, the Student MultivariateCalculus package is installed to gain Context Panel access to the MultiInt command.
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 the components of F in a free vector. Context Panel: Evaluate and Display Inline
Context Panel: Student Vector Calculus≻Conversions≻To Vector Field
Context Panel: Assign to a Name≻F
x,y,z = →to Vector Field →assign to a nameF
Obtain ∇·F, the divergence of F
Common Symbols palette: Del and dot-product operators
Context Panel: Evaluate and Display Inline
∇·F = 3
Integrate ∇·F over the interior of the box
Tools≻Load Package: Student Multivariate Calculus
Loading Student:-MultivariateCalculus
Write the numeral 3
Context Panel: Student Multivariate Calculus≻Integrate≻Iterated Complete the dialogs as per the figures below.
Context Panel: Evaluate Integral
3→MultiInt∫12∫34∫563ⅆzⅆyⅆx=3
Use a task template to compute the flux of F through the boundaries of R
Tools≻Tasks≻Browse: Calculus - Vector≻Integration≻Flux≻3-D≻Through a Box
Flux through a Box
For the Vector Field:
Select Coordinate SystemCartesian [x,y,z]Cartesian - othercylindricalsphericalbipolarcylindricalbisphericalcardioidalcardioidcylindricalcasscylindricalconicalellcylindricalhypercylindricalinvcasscylindricallogcylindricallogcoshcylindricaloblatespheroidalparaboloidalparacylindricalprolatespheroidalrosecylindricalsixspheretangentcylindricaltangentspheretoroidal
In the Divergence theorem, the volume integral on the left and the flux on the right both have the value 3.
Maple Solution - Coded
Install the Student VectorCalculus package.
withStudent:-VectorCalculus:
Set display of vectors via BasisFormat.
BasisFormatfalse:
Define F via the VectorField command.
F≔VectorFieldx,y,z:
Invoke the Divergence command.
DivergenceF = 3
Use the int command to integrate the divergence of F over R
int3,x,y,z=Parallelepiped1..2,3..4,5..6,output=integral
∫56∫34∫123ⅆxⅆyⅆz
int3,x,y,z=Parallelepiped1..2,3..4,5..6 = 3
Use the Flux command to obtain the flux of F through the boundaries of R
FluxF,Box1..2,3..4,5..6,output=integral
∫56∫341ⅆsⅆt+∫56∫121ⅆsⅆt+∫34∫121ⅆsⅆt
FluxF,Box1..2,3..4,5..6 = 3
Figure 9.8.3(a) is drawn with he Flux command by including the option "output = plot". The actual code for the figure is hidden behind the table cell in which the figure appears.
<< Previous Example Section 9.8 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
