Conservative Forces and Scalar Potentials: Needs "Vectoranalysis'" Curl X 2, 2 Y Z, Y 2, Cartesian X, Y, Z
Conservative Forces and Scalar Potentials: Needs "Vectoranalysis'" Curl X 2, 2 Y Z, Y 2, Cartesian X, Y, Z
Conservative Forces and Scalar Potentials: Needs "Vectoranalysis'" Curl X 2, 2 Y Z, Y 2, Cartesian X, Y, Z
SCALAR POTENTIALS
In our study of vector fields, we have encountered several types of conservative forces. If a force is conservative, it has a number of important properties. It is important to note that any one of the properties listed below
implies all the others; in other words, if one of these properties is true for a vector field, then they are all true.
If a force is conservative :
F = 0
the value of the line integral F dl does not depend on the path, but only on the endpoints
b
You are familiar with conservative forces; the gravitational force between point masses and the
electrical force between point charges are both examples of conservative forces. This means that
the work done in moving a particle from one height to another in a gravitational field depends
only the mass of the object and the difference in heights; the actual path taken between the starting and end points is irrelevant.
We can easily determine whether a force is conservative by computing its curl. For instance,
your last homework asked you to determine the line integral of the function :
`
`
`
v = x2 x + 2 y z y + y2 z
Over several paths from (0, 0, 0) to (1, 1, 1). If we knew a priori that the force were conservative, we would only have to calculate one line integral between these two points, knowing that
the line integral for all other paths would be equal. The simplest way to determine if a force is
conservative is to take its curl :
Needs"VectorAnalysis`"
Curlx ^ 2, 2 y z, y ^ 2, Cartesianx, y, z
0, 0, 0
1,1,1
Knowing that the curl is zero, we now know that the value of 0,0,0 v dl will be the same for all
possible paths between these limits.
scalarpotentials.nb
We would like to be able to figure out the scalar potential that generates the vector field of the
force. If we could do this, calculating line integrals becomes almost trivial, requiring only the
simplest integrations. Let's review why this is the case. As we know, if F is conservative, then
it can be derived from a scalar potential such that:
F = -f = so that the line integral becomes :
W = F dl = -
b
= -
f
x
dx +
f
y
f `
f `
f `
xyz
x
y
z
f `
f `
f `
`
`
`
x+
y+
z dx x + dy y + dz z
x
y
z
dy +
f
z
dz = - df = - f b - f a
b
The next to last equality occurs because the total derivative of a multivariate function f (x, y, z)
is simply :
df =
f
x
dx +
f
y
dy +
f
z
dz
Let's now consider how to determine the scalar potential that generates a conservative force. If
we assume that such a potential exists, then F = -f, and we can equate the components of these
two vectors:
f ` f `
f `
`
`
`
x+
y+
z
F = Fx x + Fy y + Fz z and f =
x
y
z
Equating components, we get the three equations :
Fx = -
f
x
; Fy = -
f
y
; Fz = -
f
z
We will use the relationships to determine the scalar potential generating the function
`
`
`
F = x2 x + 2 y z y + y2 z
For this given F, we know :
Fx= x2 = - f x; F y = 2 y z = - f/y ; Fz = y2 = - f z
We are ready to begin calculating the scalar function f that generates F. We will start with the
scalarpotentials.nb
f
x
= x2 df = - x2 dx f =
- x3
3
+ g y, z
(1)
It is important to understand the appearance of the g (y, z) term. In eq. (1) above, we integrate
with respect to x to determine the x depence of f. You learned in Calc I that indefinite integration produces a constant. In this case, g(y,z) plays that role, since g(y,z) is constant with respect
to x. When we take the partial integral with respect to x, we cannot rule out the presence of any
functions that depend only on y or z, so our constant of integration is written most generally as a
function depending on y and z, i.e., g(y,z).
Our next step is to take the potential function as we have it in eq. (1), differentiate that with
respect to y and compare this expression of f/y with the expression for f/y found in the
equation for the force F. Taking this derivative:
f =
g y, z
-x3
3
+ g y, z
f
y
g y, z
y
= -2 y z
This last equality allows us to partially integrate with respect to y and find :
y
= -2 y z g y, z = - 2 y z dy = - y2 z + h z
(2)
As before, our indefinite integral produces a constant of integration; since we integrated with
respect to y, the most general way to write the constant of integration is to treat it as a function
solely of z.
When we now combine eq. (1) and eq. (2), we see the scalar potential becomes :
f=
- x3
3
- y2 z + h z
(3)
Now, we take the potential function as we have it in (3), differentiate it with respect to z, and
compare that result to the expression for f/z in the original force equation. Taking the partial
derivative of eq. (3) with respect to z yields:
f =
-x3
3
- y2 z + h z
f
z
= - y2 +
dh z
dz
= -y2
This last differentiation shows that dh(z)/dz = 0, so that h is at most a numerical constant. We
have now fully determined that the scalar potential for this force is written as:
scalarpotentials.nb
f =
-x3
3
-y z = 2
x3
3
+ y2 z
Let' s check if this actually reproduces our force by taking the negative of the gradient of this
function :
`
`
`
F = -f = -x2 x + 2 y z y + y2 z
and we see that the force is in fact derived from the potential in eq. (4)
Now that we know the scalar potential, we can see its utility in calculating work integrals. As
shown above, the line integral for a conservative can be written as :
W = - df = -f b - f a
b
For the vector we have been using in this example (and the vector you used in homework), we
can calculate the line integral simply by evaluating the scalar potential as written in eq. (5) at the
end points. In other words :
df = -f 1, 1, 1 - f 0, 0, 0 = b
4
- 12 1 - 0 =
3
3
-1
A Worked Example :
Let' s consider the function :
`
`
`
F = 3 x2 y2 x + 2 x3 y + cos z y - y sin z z
First, we will determine whether it is conservative. If it is, we will compute the scalar potential
that generates it, and then calculate the line integral of this force over the path from (0, 0, 0) to
(1, 1, p)
I recommend you do this by hand as an exercise :
F =
`
x
`
y
`
z
3 x2 y2 2 x3 y + cos z -y sin z
`
`
`
x -sin z - -sinz - y 0 - 0 + z 6 x2 y - 6 x2 y = 0, 0, 0
scalarpotentials.nb
Since the curl of F is zero, we know this is a conservative force. Now, we equate the x component of F with the x component of f, and find:
3 x2 y2 = -
f
x
f = 3 x2 y2 dx = -x3 y2 + g y, z
Now, we differentiate the expression for f immediately above with respect to y and equate to the
y component of the force F:
f
y
This implies :
g y, z
y
= - 2 x3 y +
g y, z
y
= -2 x3 y + cos z
f = -2 x3 y - y cos z + h z
We now differentiate this expression for f with respect to z, and equate the result to the z component of the force F:
f
z
= y sin z +
dh z
dz
= y sin z
this last relationship makes it clear h (z) is zero (or some other numerical constant), and we can
write our scalar potential as :
f = -2 x3 y + y cos z
With knowledge of the potential, we can calculate the work done between any two points simply
by evaluating the potential at those points. For us, the work done in exerting this force from (0,
0, 0) to (1, 1, p) is :
W = - f 1, 1, p - f 0, 0, 0
Let' s use this opportunity to become acquainted with another interesting Mathematica operation.
Look carefully at the input below which will evaluate our potential at one of the points :
Clearx, y, z
2x3 y y Cosz;
. x 1, y 1, z
After we cleared variables, we defined our scalar function. Then we made use of the "slash-dot"
operation; this is a rule that says to replace x by 1, and y by 1 and z by p. Similarly, we can find
scalarpotentials.nb
. x 0, y 0, z 0
And you obtain zero as you could easily find by inspection. Clearly, the value of the line integral is simply - 1 - 0. Since we never redefined f, we can continue to use its original defintion in
our succeeding "/." assignments.
Finally, we show how we could have computed the potential at the two limits in one command :
. x 1, y 1, z , x 0, y 0, z 0
1, 0