Div, Grad, and Curl: 7.1 The Operator
Div, Grad, and Curl: 7.1 The Operator
Div, Grad, and Curl: 7.1 The Operator
= 0 on D = constant. (7.3)
Moreover, for any scalar c, we have:
1
7.2 Divergence
Let f : D Rn , D Rn , be a dierentiable vector eld. (Note that both
spaces are n-dimensional.) Let f1 , f2 , . . . , fn be the component (scalar) elds
of f . The divergence of f is dened to be
n
fj
div(f ) = f = . (7.5)
j=1
x j
n
tr(A) = aii .
i=1
f = tr(Df ). (7.7)
Let be a twice dierentiable scalar eld. Then its Laplacian is dened
to be
2 = (). (7.8)
It follows from (7.1),(7.5),(7.6) that
2 2 2
=
2
+ 2 + + 2 . (7.9)
x21 x2 xn
One says that is harmonic i 2 = 0. Note that we can formally
consider the dot product
2 n
=( ,..., )( ,..., )= . (7.10)
x1 xn x1 xn j=1
x2j
Then we have
2
2 = ( ). (7.11)
These last two examples are special cases of the fact, which we mention
without proof, that for any function f : D C which is dierentiable in the
complex sense, the real and imaginary part, (f ) and (f ), are harmonic
functions. Here f is dierentiable in the complex sense if its total derivative
Df at a point z D, a priori a R-linear map from C to itself, is in fact given
by multiplication with a complex number, which we then call f (z). More
concretely,
( a b) this means that the matrix of Df in the basis 1, i is of the form
b a for some real numbers a, b. We then have f (z) = a + bi. There is
a large supply of such functions since any f given (locally) by a convergent
power series in z is complex dierentiable.
In (i) we can take f (z) = ez = ex+iy = ex cos(y) + iex sin(y) and in
(ii) we can take f (z) = log(z) = log(rei ) = log(r) + i but we must be
careful about the domain. To have a well dened argument for all z D
we must make a cut in the plane and can only dene f on, for example,
D = {z = x + iy| y = 0 x > 0} or D = {z = x + iy| y = 0 x < 0}. But
the union of D and D is C {0} as in (ii) .
(iii) D = Rn {0}; (x1 , x2 , ..., xn ) = (x21 + x22 + + x2n )/2 = r for some
xed R.
Then x i
= r1 xri = r2 xi , and
2
= ( 2)r4 xi xi + r2 1.
x2i
Hence 2 = ni=1 (( 2)r4 x2i + r2 ) = ( 2 + n)r2 .
So is harmonic for = 0 or = 2 n ( = 1 for n = 3).
3
7.3 Cross product in R3
The three-dimensional space is very special in that it admits a vector prod-
uct, often called the cross product. Let i,j,k denote the standard basis of
R3 . Then, for all pairs of vectors v = xi + yj + zk and v = x i + y j + z k,
the cross product is dened by
( i j k
)
v v = det x y z = (yz y z)i (xz x z)j + (xy x y)k. (7.12)
x y z
Corollary: v v = 0.
4
More generally, the same argument shows that the (signed) volume of the
parallelepiped spanned by any three vectors u, v, v is u (v v ).
(a) (h) = 0.
(b) ( f ) = 0.
( ) ( ) ( )
2h 2h 2h 2h 2h 2h
= i+ j+ k.
yz zy zx xz xy yx
5
( ) ( ) ( )
2R 2Q 2R 2P 2Q 2P
= + + + .
xy xz yx yz zx zy
2R 2R
Again, since f is C 2 , xy
= yx
, etc., and we get the assertion.
Done.
Warning: There exist twice dierentiable scalar (resp. vector) elds
h (resp. f ), which are not C 2 , for which (a) (resp. (b)) does not hold.
When the vector eld f represents uid ow, it is often called irrota-
tional when its curl is 0. If this ow describes the movement of water in a
stream, for example, to be irrotational means that a small boat being pulled
by the ow will not rotate about its axis. We will see later in this chapter
the condition f = 0 occurs naturally in a purely mathematical setting
as well.
( ) ( ) ( ( ) ( ))
x y x y
= i+ j+ k
z x + y2
2 z x2 + y 2 x x2 + y 2 y x2 + y 2
[ ]
(x2 + y 2 ) + 2x2 (x2 + y 2 ) 2y 2
= k = 0.
(x2 + y 2 )2 (x2 + y 2 )2
(ii) Let m be any integer = 3, D = R3 {0}, and
f (x, y, z) = r1m (xi + yj + zk), where r = x2 + y 2 + z 2 . Show that f is not
the curl of another vector eld. Indeed, suppose f = g. Then, since f
is C 1 , g will be C 2 , and by the Proposition proved above,
f = ( g) would be zero. But,
( ) (
x y z )
f = , , m, m, m
x y z r r r
6
1 ( ) 1
= 2m
3r m
m(x 2
+ y 2
+ z 2 m2
)r = m (3 m).
r r
This is non-zero as m = 3. So f is not a curl.
{(x, y, z) R(
3
|(x, y) )D} by setting f (x, y, z) = g(x, y) = P i + Qj. Then
i j k ( )
f = det x
y z = Q
y
P
x
k, because P
z
= Q
z
= 0. Thus we
P Q 0
get:
Q P
( f ) k = . (7.15)
y x
And Greens theorem becomes:
H
Theorem 1 ( f ) k dx dy = C P dx + Q dy
7
Proposition 1 Let g : D R2 , D R2 open and simply connected,
g = (P, Q), be a C 1 vector field. Set f (x, y, z) = g(x, y), for all (x, y, z) R3
with (x, y) D. Suppose f = 0. Then g is conservative on D.
H
Proof: Since f = 0, Theorem 1 implies that C P dx + Q dy = 0 for all
Jordan
H curves C contained in D. In fact, f = 0 also implies that
C
P dx + Q dy = 0 for all closed curves but we wont prove this. Hence f is
conservative. Done.
Again, dene f (x, y, z) to be g(x, y) for all (x, y, z) in R3 such that (x, y)
D. Since g is evidently C 1 , f will be C 1 as well. By the Proposition above, it
will suce to check if f is irrotational, i.e., f = 0, on D R. This was
already shown in Example (i) of section 4 of this chapter. So g is conservative.