Calc3.6

The Divergence

Definition

We've discussed partial derivatives, and we've even discussed a way to take a sort of "complete" derivative of multivariable functions (the gradient). But the gradient only works for scalar functions; how about a "complete" kind of derivative for vector fields?

In two dimensions, there is just one, but in three dimensions, there are, for reasons which won't appear until much later in your mathematical education. The first of these is called the divergence, which we investigate now.

Consider a vector field ${\vec {F}}=f_{1}{\vec {i}}+f_{2}{\vec {j}}+f_{3}{\vec {k}}$, which we interpret as representing the velocity of a fluid occupying space, or a very thin layer of flowing fluid on a plane. Suppose we wished to find out if there are any sources of fluid (like a faucet), or drains - how might we go about this?

Well, consider some specific point, $(x_{0},y_{0},z_{0})={\vec {p}}$. If we surround ${\vec {p}}$ with a sphere, call it S, then obviously

$\iint \limits _{S}\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\;\;\;\subset \!\supset {\vec {F}}dS$

will contain some contribution from whatever "flow" is coming out of, or going into, ${\vec {p}}$. It will, of course, also contains contributions from the various other points in its interior, but, if we let the radius r of the sphere S go to 0, then this integral will ONLY measure the flow at ${\vec {p}}$. Hence we define the divergence of a vector field ${\vec {F}}$ as

$\operatorname {div} \ {\vec {F}}=\lim _{r\to 0}\left({\frac {1}{{\frac {4}{3}}\pi r^{3}}}\iint \limits _{S}\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\;\;\;\subset \!\supset {\vec {F}}\cdot {\vec {n}}_{S}dS\right)$

where ${\vec {n}}_{S}$ is the outward pointing surface normal of unit length to S.

Computation

See the proof page for details.

This is an unwieldy definition, so we frequently use the equivalent statement

$\operatorname {div} {\vec {F}}={\frac {\partial f_{1}}{\partial x}}+{\frac {\partial f_{2}}{\partial y}}+{\frac {\partial f_{3}}{\partial z}}$

or in two dimensions, just ${\frac {\partial f_{1}}{\partial x}}+{\frac {\partial f_{2}}{\partial y}}$. Note that if we were to define a "vector" $\nabla ={\vec {i}}\partial _{x}+{\vec {j}}\partial _{y}+{\vec {k}}\partial _{z}$, we would have $\operatorname {div} \ {\vec {F}}=\nabla \cdot {\vec {F}}$. Note that this isn't technically a vector, since the component are not number but differential operators, and the kind of "multiplication" we're doing isn't really multiplication, since you can't multiply an operator by a function; instead we're differentiating the components. Nevertheless, it's a handy notational convenience, and it actually does point to a deeper truth we won't begin to address in this course, so we'll continue to use $\nabla \cdot {\vec {F}}$ for the divergence.

The Curl

Since we've taken the dot product of our new differential operator $\nabla$ with a vector field, it only makes sense to see what happens when we take the cross product as well.

$\nabla \times {\vec {F}}={\begin{vmatrix}{\vec {i}}&{\vec {j}}&{\vec {k}}\\\partial _{x}&\partial _{y}&\partial _{z}\\f_{1}&f_{2}&f_{3}\end{vmatrix}}=\left({\frac {\partial f_{3}}{\partial y}}-{\frac {\partial f_{2}}{\partial z}}\right){\vec {i}}-\left({\frac {\partial f_{3}}{\partial x}}-{\frac {\partial f_{1}}{\partial z}}\right){\vec {j}}+\left({\frac {\partial f_{2}}{\partial x}}-{\frac {\partial f_{1}}{\partial y}}\right){\vec {k}}$

$=\left({\frac {\partial f_{3}}{\partial y}}-{\frac {\partial f_{2}}{\partial z}}\right){\vec {i}}+\left({\frac {\partial f_{1}}{\partial z}}-{\frac {\partial f_{3}}{\partial x}}\right){\vec {j}}+\left({\frac {\partial f_{2}}{\partial x}}-{\frac {\partial f_{1}}{\partial y}}\right){\vec {k}}$.

Of course, this is useless if we don't know what it measures. Let's investigate:

Here, we see space with the positive x-axis poking out of the screen towards us. Note that the difference in values and direction of the two partial derivatives is cause by the fact that the direction of the lines coming out of their common basepoint is rotating around in the neighborhood of that point.

The first component, $\left(\partial _{y}f_{3}-\partial _{z}f_{2}\right)$, gives us the extent of the vector $\nabla \times {\vec {F}}$ in the x-axis direction. The expressions gives us, basically, how fast the fluid moving "up" along the z-axis "tips" in the direction of the y-axis, minus how much the fluid moving in the direction of the y axis tips in the direction of the z axis.