# Divergence

For other uses, see Divergence (disambiguation).

In vector calculus, divergence is a vector operator that measures the magnitude of a vector field's source or sink at a given point, in terms of a signed scalar. More technically, the divergence represents the volume density of the outward flux of a vector field from an infinitesimal volume around a given point.

For example: consider air as it is heated or cooled. The relevant vector field for this example is the velocity of the moving air at a point. If air is heated in a region it will expand in all directions such that the velocity field points outward from that region. Therefore the divergence of the velocity field in that region would have a positive value, as the region is a source. If the air cools and contracts, the divergence has a negative value, as the region is a sink.

## Definition of divergence

In physical terms, the divergence of a three-dimensional vector field is the extent to which the vector field flow behaves like a source or a sink at a given point. It is a local measure of its "outgoingness"—the extent to which there is more exiting an infinitesimal region of space than entering it. If the divergence is nonzero at some point then there must be a source or sink at that position.[1] (Note that we are imagining the vector field to be like the velocity vector field of a fluid (in motion) when we use the terms flow, sink and so on.)

More rigorously, the divergence of a vector field F at a point p is defined as the limit of the net flow of F across the smooth boundary of a three-dimensional region V divided by the volume of V as V shrinks to p. Formally,

$\operatorname{div}\,\mathbf{F}(p) = \lim_{V \rightarrow \{p\}} \iint_{S(V)} {\mathbf{F}\cdot\mathbf{n} \over |V| } \; dS$

where |V | is the volume of V, S(V) is the boundary of V, and the integral is a surface integral with n being the outward unit normal to that surface. The result, div F, is a function of p. From this definition it also becomes explicitly visible that div F can be seen as the source density of the flux of F.

In light of the physical interpretation, a vector field with constant zero divergence is called incompressible or solenoidal – in this case, no net flow can occur across any closed surface.

The intuition that the sum of all sources minus the sum of all sinks should give the net flow outwards of a region is made precise by the divergence theorem.

### Application in Cartesian coordinates

Let x, y, z be a system of Cartesian coordinates in 3-dimensional Euclidean space, and let ijk be the corresponding basis of unit vectors.

The divergence of a continuously differentiable vector field F = U i + V j + W k is equal to the scalar-valued function:

$\operatorname{div}\,\mathbf{F} = \nabla\cdot\mathbf{F} =\frac{\partial U}{\partial x} +\frac{\partial V}{\partial y} +\frac{\partial W}{\partial z }.$

Although expressed in terms of coordinates, the result is invariant under orthogonal transformations, as the physical interpretation suggests.

The common notation for the divergence ∇ · F is a convenient mnemonic, where the dot denotes an operation reminiscent of the dot product: take the components of ∇ (see del), apply them to the components of F, and sum the results. Because applying an operator is different from multiplying the components, this is considered an abuse of notation.

The divergence of a continuously differentiable second-order tensor field $\underline{\underline{\epsilon}}$ is a first-order tensor field:[2]

$\overrightarrow{\operatorname{div}}\,(\mathbf{\underline{\underline{\epsilon}}}) = \begin{bmatrix} \frac{\partial \epsilon_{xx}}{\partial x} +\frac{\partial \epsilon_{yx}}{\partial y} +\frac{\partial \epsilon_{zx}}{\partial z} \\[6pt] \frac{\partial \epsilon_{xy}}{\partial x} +\frac{\partial \epsilon_{yy}}{\partial y} +\frac{\partial \epsilon_{zy}}{\partial z} \\[6pt] \frac{\partial \epsilon_{xz}}{\partial x} +\frac{\partial \epsilon_{yz}}{\partial y} +\frac{\partial \epsilon_{zz}}{\partial z} \end{bmatrix}$

### Cylindrical coordinates

For a vector expressed in cylindrical coordinates as

$\mathbf F = \mathbf e_r F_r + \mathbf e_\theta F_{\theta} + \mathbf e_z F_z,$

where ea is the unit vector in direction a, the divergence is[3]

$\operatorname{div}\, \mathbf F = \nabla\cdot\mathbf F = \frac1r \frac{\partial}{\partial r} (rF_r) + \frac1r \frac{\partial F_\theta}{\partial\theta} + \frac{\partial F_z}{\partial z}\, .$

### Spherical coordinates

In spherical coordinates, with $\theta$ the angle with the z axis and $\phi$ the rotation around the z axis, the divergence reads[4]

$\operatorname{div}\, \mathbf F = \nabla\cdot\mathbf F = \frac1{r^2} \frac{\partial}{\partial r}(r^2 F_r) + \frac1{r\sin\theta} \frac{\partial}{\partial \theta} (\sin\theta\, F_\theta) + \frac1{r\sin\theta} \frac{\partial F_\phi}{\partial \phi}.$

## Decomposition theorem

It can be shown that any stationary flux v(r) which is at least two times continuously differentiable in ${\mathbb R}^3$ and vanishes sufficiently fast for |r| → ∞ can be decomposed into an irrotational part E(r) and a source-free part B(r). Moreover, these parts are explicitly determined by the respective source-densities (see above) and circulation densities (see the article Curl):

For the irrotational part one has

$\mathbf E=-\nabla \Phi(\mathbf r)\, ,$

with

$\Phi (\mathbf{r})=\int_{\mathbb R^3}\,{\rm d}^3\mathbf r'\;\frac{\operatorname{div}\,\mathbf{v}(\mathbf{r}')}{4\pi|\mathbf{r}-\mathbf{r}'|}\,.$

The source-free part, B, can be similarly written: one only has to replace the scalar potential Φ(r) by a vector potential A(r) and the terms −Φ by +×A, and the source-density div v by the circulation-density ×v.

This "decomposition theorem" is in fact a by-product of the stationary case of electrodynamics. It is a special case of the more general Helmholtz decomposition which works in dimensions greater than three as well.

## Properties

The following properties can all be derived from the ordinary differentiation rules of calculus. Most importantly, the divergence is a linear operator, i.e.

$\operatorname{div}( a\mathbf{F} + b\mathbf{G} ) = a\;\operatorname{div}( \mathbf{F} ) + b\;\operatorname{div}( \mathbf{G} )$

for all vector fields F and G and all real numbers a and b.

There is a product rule of the following type: if $\varphi$ is a scalar valued function and F is a vector field, then

$\operatorname{div}(\varphi \mathbf{F}) = \operatorname{grad}(\varphi) \cdot \mathbf{F} + \varphi \;\operatorname{div}(\mathbf{F}),$

or in more suggestive notation

$\nabla\cdot(\varphi \mathbf{F}) = (\nabla\varphi) \cdot \mathbf{F} + \varphi \;(\nabla\cdot\mathbf{F}).$

Another product rule for the cross product of two vector fields F and G in three dimensions involves the curl and reads as follows:

$\operatorname{div}(\mathbf{F}\times\mathbf{G}) = \operatorname{curl}(\mathbf{F})\cdot\mathbf{G} \;-\; \mathbf{F} \cdot \operatorname{curl}(\mathbf{G}),$

or

$\nabla\cdot(\mathbf{F}\times\mathbf{G}) = (\nabla\times\mathbf{F})\cdot\mathbf{G} - \mathbf{F}\cdot(\nabla\times\mathbf{G}).$

The Laplacian of a scalar field is the divergence of the field's gradient:

$\operatorname{div}(\nabla\varphi) = \Delta\varphi.$

The divergence of the curl of any vector field (in three dimensions) is equal to zero:

$\nabla\cdot(\nabla\times\mathbf{F})=0$

If a vector field F with zero divergence is defined on a ball in R3, then there exists some vector field G on the ball with F = curl(G). For regions in R3 more complicated than this, the latter statement might be false (see Poincaré lemma). The degree of failure of the truth of the statement, measured by the homology of the chain complex

$\{\mbox{scalar fields on }U\} \;$
$\to\{\mbox{vector fields on }U\} \;$
$\to\{\mbox{vector fields on }U\} \;$
$\to\{\mbox{scalar fields on }U\} \;$

(where the first map is the gradient, the second is the curl, the third is the divergence) serves as a nice quantification of the complicatedness of the underlying region U. These are the beginnings and main motivations of de Rham cohomology.

## Relation with the exterior derivative

One can express the divergence as a particular case of the exterior derivative, which takes a 2-form to a 3-form in R3. Define the current two form

$j =F_1\ dy\wedge dz + F_2\ dz\wedge dx + F_3\ dx\wedge dy$.

It measures the amount of "stuff" flowing through a surface per unit time in a "stuff fluid" of density $\rho = 1 dx\wedge dy\wedge dz$ moving with local velocity F. Its exterior derivative $d j$ is then given by

$d j = \left( \frac{\partial F_1}{\partial x} +\frac{\partial F_2}{\partial y} +\frac{\partial F_3}{\partial z} \right) dx\wedge dy\wedge dz = (\nabla\cdot \mathbf{F}) \rho$

Thus, the divergence of the vector field F can be expressed as:

$\nabla \cdot \mathbf{F} = \star {\mathrm d} {\star {\mathbf{F}^\flat}}$

Here the superscript $\flat$ is one of the two musical isomorphisms, and $\star$ is the Hodge dual. Note however that working with the current two form itself and the exterior derivative is usually easier than working with the vector field and divergence, because unlike the divergence, the exterior derivative commutes with a change of (curvilinear) coordinate system.

## Generalizations

The divergence of a vector field can be defined in any number of dimensions. If

$\mathbf{F}=(F_1, F_2, \dots, F_n),$

in a Euclidean coordinate system where $\mathbf{x}=(x_1, x_2, \dots, x_n)$ and $d\mathbf{x}=(dx_1, dx_2, \dots, dx_n)$, define

$\operatorname{div}\,\mathbf{F} = \nabla\cdot\mathbf{F} =\frac{\partial F_1}{\partial x_1} +\frac{\partial F_2}{\partial x_2}+\cdots +\frac{\partial F_n}{\partial x_n}.$

The appropriate expression is more complicated in curvilinear coordinates.

In the case of one dimension, a "vector field" is simply a regular function, and the divergence is simply the derivative.

For any n, the divergence is a linear operator, and it satisfies the "product rule"

$\nabla\cdot(\varphi \mathbf{F}) = (\nabla\varphi) \cdot \mathbf{F} + \varphi \;(\nabla\cdot\mathbf{F}).$

for any scalar-valued function $\varphi$.

The divergence can be defined on any manifold of dimension n with a volume form (or density) $\mu$ e.g. a Riemannian or Lorentzian manifold. Generalising the construction of a two form for a vector field on $\mathbb{R}^3$, on such a manifold a vector field X defines a n−1 form $j = i_X \mu$ obtained by contracting X with $\mu$. The divergence is then the function defined by

$d j = \operatorname{div}(X) \mu$

Standard formulas for the Lie derivative allow us to reformulate this as

$\mathcal{L}_X \mu = \operatorname{div}(X) \mu$

This means that the divergence measures the rate of expansion of a volume element as we let it flow with the vector field.

On a Riemannian or Lorentzian manifold the divergence with respect to the metric volume form can be computed in terms of the Levi Civita connection $\nabla$

$\operatorname{div}(X) = \nabla\cdot X = X^a_{;a}$

where the second expression is the contraction of the vector field valued 1-form $\nabla X$ with itself and the last expression is the traditional coordinate expression used by physicists.

An equivalent expression without using connection is

$\operatorname{div}(X) = \frac{1}{\sqrt{\operatorname{det} g}} \partial_a (\sqrt{\operatorname{det} g} X^a)$

where $g$ is the metric and $\partial_a$ denotes partial derivative with respect to coordinate $x^a$.

Divergence can also be generalised to tensors. In Einstein notation, the divergence of a contravariant vector $F^\mu$ is given by

$\nabla\cdot\mathbf{F}=\nabla_\mu F^\mu$

where $\nabla_\mu$ is the covariant derivative.

Equivalently, some authors define the divergence of any mixed tensor by using the "musical notation #":

If T is a (p,q)-tensor (p for the contravariant vector and q for the covariant one), then we define the divergence of T to be the (p,q−1)-tensor

$(\operatorname{div} T)(Y_1,...,Y_{q-1}) = \operatorname{trace}(X\mapsto \#( \nabla T)(X,\cdot,Y_1,...,Y_{q-1} ))$

that is we trace the covariant derivative on the first two covariant indices.