# Divergence theorem

{{#invoke:Hatnote|hatnote}}Template:Main other {{#invoke:Hatnote|hatnote}}Template:Main other {{#invoke: Sidebar | collapsible }}

In vector calculus, the divergence theorem, also known as Gauss's theorem or Ostrogradsky's theorem, is a result that relates the flow (that is, flux) of a vector field through a surface to the behavior of the vector field inside the surface.

More precisely, the divergence theorem states that the outward flux of a vector field through a closed surface is equal to the volume integral of the divergence over the region inside the surface. Intuitively, it states that the sum of all sources minus the sum of all sinks gives the net flow out of a region.

The divergence theorem is an important result for the mathematics of engineering, in particular in electrostatics and fluid dynamics.

In physics and engineering, the divergence theorem is usually applied in three dimensions. However, it generalizes to any number of dimensions. In one dimension, it is equivalent to the fundamental theorem of calculus. In two dimensions, it is equivalent to Green's theorem.

The theorem is a special case of the more general Stokes' theorem.

## Intuition

If a fluid is flowing in some area, then the rate at which fluid flows out of a certain region within that area can be calculated by adding up the sources inside the region and subtracting the sinks. The fluid flow is represented by a vector field, and the vector field's divergence at a given point describes the strength of the source or sink there. So, integrating the field's divergence over the interior of the region should equal the integral of the vector field over the region's boundary. The divergence theorem says that this is true.

The divergence theorem is thus a conservation law which states that the volume total of all sinks and sources, that is the volume integral of the divergence, is equal to the net flow across the volume's boundary.

## Mathematical statement The divergence theorem can be used to calculate a flux through a closed surface that fully encloses a volume, like any of the surfaces on the left. It can not directly be used to calculate the flux through surfaces with boundaries, like those on the right. (Surfaces are blue, boundaries are red.)

Suppose Template:Mvar is a subset of $\mathbb {R} ^{n}$ (in the case of n = 3, V represents a volume in 3D space) which is compact and has a piecewise smooth boundary Template:Mvar (also indicated with V = S). If F is a continuously differentiable vector field defined on a neighborhood of Template:Mvar, then we have:

Template:Oiint

The left side is a volume integral over the volume Template:Mvar, the right side is the surface integral over the boundary of the volume Template:Mvar. The closed manifold V is quite generally the boundary of Template:Mvar oriented by outward-pointing normals, and n is the outward pointing unit normal field of the boundary V. (dS may be used as a shorthand for ndS.) The symbol within the two integrals stresses once more that V is a closed surface. In terms of the intuitive description above, the left-hand side of the equation represents the total of the sources in the volume Template:Mvar, and the right-hand side represents the total flow across the boundary V.

### Corollaries

By applying the divergence theorem in various contexts, other useful identities can be derived (cf. vector identities).

• Applying the divergence theorem to the product of a scalar function Template:Mvar and a vector field F, the result is
Template:Oiint
A special case of this is F = ∇ f, in which case the theorem is the basis for Green's identities.
• Applying the divergence theorem to the cross-product of two vector fields F × G, the result is
Template:Oiint
• Applying the divergence theorem to the product of a scalar function, f, and a non-zero constant vector c, the following theorem can be proven:
Template:Oiint
• Applying the divergence theorem to the cross-product of a vector field F and a non-zero constant vector c, the following theorem can be proven:
Template:Oiint

## Example

Suppose we wish to evaluate

Template:Oiint

where Template:Mvar is the unit sphere defined by

$S=\left\{x,y,z\in \mathbb {R} ^{3}\ :\ x^{2}+y^{2}+z^{2}=1\right\}.$ and F is the vector field

${\mathbf {F} }=2x{\mathbf {i} }+y^{2}{\mathbf {j} }+z^{2}{\mathbf {k} }.$ The direct computation of this integral is quite difficult, but we can simplify the derivation of the result using the divergence theorem, because the divergence theorem says that the integral is equal to:

$\iiint _{W}(\nabla \cdot \mathbf {F} )\,dV=2\iiint _{W}(1+y+z)\,dV=2\iiint _{W}dV+2\iiint _{W}y\,dV+2\iiint _{W}z\,dV.$ where Template:Mvar is the unit ball:

$W=\left\{x,y,z\in \mathbb {R} ^{3}\ :\ x^{2}+y^{2}+z^{2}\leq 1\right\}.$ Since the function Template:Mvar is positive in one hemisphere of Template:Mvar and negative in the other, in an equal and opposite way, its total integral over Template:Mvar is zero. The same is true for Template:Mvar:

$\iiint _{W}y\,dV=\iiint _{W}z\,dV=0.$ Therefore,

Template:Oiint

because the unit ball Template:Mvar has volume {{ safesubst:#invoke:Unsubst||B=4π/3}}. ## Applications ### Differential form and integral form of physical laws As a result of the divergence theorem, a host of physical laws can be written in both a differential form (where one quantity is the divergence of another) and an integral form (where the flux of one quantity through a closed surface is equal to another quantity). Three examples are Gauss's law (in electrostatics), Gauss's law for magnetism, and Gauss's law for gravity. #### Continuity equations {{#invoke:main|main}} Continuity equations offer more examples of laws with both differential and integral forms, related to each other by the divergence theorem. In fluid dynamics, electromagnetism, quantum mechanics, relativity theory, and a number of other fields, there are continuity equations that describe the conservation of mass, momentum, energy, probability, or other quantities. Generically, these equations state that the divergence of the flow of the conserved quantity is equal to the distribution of sources or sinks of that quantity. The divergence theorem states that any such continuity equation can be written in a differential form (in terms of a divergence) and an integral form (in terms of a flux). ### Inverse-square laws Any inverse-square law can instead be written in a Gauss' law-type form (with a differential and integral form, as described above). Two examples are Gauss' law (in electrostatics), which follows from the inverse-square Coulomb's law, and Gauss' law for gravity, which follows from the inverse-square Newton's law of universal gravitation. The derivation of the Gauss' law-type equation from the inverse-square formulation (or vice versa) is exactly the same in both cases; see either of those articles for details. ## History The theorem was first discovered by Lagrange in 1762, then later independently rediscovered by Gauss in 1813, by Ostrogradsky, who also gave the first proof of the general theorem, in 1826, by Green in 1828, etc. Subsequently, variations on the divergence theorem are correctly called Ostrogradsky's theorem, but also commonly Gauss's theorem, or Green's theorem. ## Examples To verify the planar variant of the divergence theorem for a region Template:Mvar: $R=\left\{x,y\in \mathbb {R} ^{2}\ :\ x^{2}+y^{2}\leq 1\right\},$ and the vector field: $\mathbf {F} (x,y)=2y\mathbf {i} +5x\mathbf {j} .$ The boundary of Template:Mvar is the unit circle, Template:Mvar, that can be represented parametrically by: $x=\cos(s),\quad y=\sin(s)$ such that 0 ≤ s ≤ 2π where Template:Mvar units is the length arc from the point s = 0 to the point Template:Mvar on Template:Mvar. Then a vector equation of Template:Mvar is $C(s)=\cos(s)\mathbf {i} +\sin(s)\mathbf {j} .$ At a point Template:Mvar on Template:Mvar: $P=(\cos(s),\,\sin(s))\,\Rightarrow \,\mathbf {F} =2\sin(s)\mathbf {i} +5\cos(s)\mathbf {j} .$ Therefore, {\begin{aligned}\oint _{C}\mathbf {F} \cdot \mathbf {n} \,ds&=\int _{0}^{2\pi }(2\sin(s)\mathbf {i} +5\cos(s)\mathbf {j} )\cdot (\cos(s)\mathbf {i} +\sin(s)\mathbf {j} )\,ds\\&=\int _{0}^{2\pi }(2\sin(s)\cos(s)+5\sin(s)\cos(s))\,ds\\&=7\int _{0}^{2\pi }\sin(s)\cos(s)\,ds\\&=0.\end{aligned}} Because M = 2y, {{ safesubst:#invoke:Unsubst||B=M/x}} = 0, and because N = 5x, {{ safesubst:#invoke:Unsubst||\$B=N/y}} = 0. Thus

$\iint _{R}\,\operatorname {div} \mathbf {F} \,dA=\iint _{R}\left({\frac {\partial M}{\partial x}}+{\frac {\partial N}{\partial y}}\right)\,dA=0.$ ## Generalizations

### Multiple dimensions

One can use the general Stokes' Theorem to equate the Template:Mvar-dimensional volume integral of the divergence of a vector field F over a region Template:Mvar to the (n − 1)-dimensional surface integral of F over the boundary of Template:Mvar:

$\int _{U}\nabla \cdot {\mathbf {F} }\,dV_{n}=\oint _{\partial U}{\mathbf {F} }\cdot {\mathbf {n} }\,dS_{n-1}$ This equation is also known as the Divergence theorem.

When n = 2, this is equivalent to Green's theorem.

When n = 1, it reduces to the Fundamental theorem of calculus.

### Tensor fields

{{#invoke:main|main}}

Writing the theorem in Einstein notation:

Template:Oiint

suggestively, replacing the vector field F with a rank-Template:Mvar tensor field Template:Mvar, this can be generalized to:

Template:Oiint

where on each side, tensor contraction occurs for at least one index. This form of the theorem is still in 3d, each index takes values 1, 2, and 3. It can be generalized further still to higher (or lower) dimensions (for example to 4d spacetime in general relativity).