# Divergence theorem

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

More precisely, the divergence theorem states that the outward flux of a tensor 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 (with sinks regarded as negative sources) gives the net flux out of a region.

The divergence theorem is an important result for the mathematics of physics and 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.

## 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 first order tensor (or 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.[2]

The divergence theorem is employed in any 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.[3]

## Mathematical statement

A region V bounded by the surface S = ∂V with the surface normal n
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 V is a subset of ${\displaystyle \mathbb ^}$ (in the case of n = 3, V represents a volume in three-dimensional space) which is compact and has a piecewise smooth boundary S (also indicated with V = S). If F is a continuously differentiable vector field defined on a neighborhood of V, then we have:[4]

${\displaystyle \iiint _\left(\mathbf {\nabla } \cdot \mathbf \right)\,dV=}$ ${\displaystyle \scriptstyle S}$ ${\displaystyle (\mathbf \cdot \mathbf )\,dS.}$

The left side is a volume integral over the volume V, the right side is the surface integral over the boundary of the volume V. The closed manifold V is quite generally the boundary of V 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 V, and the right-hand side represents the total flow across the boundary S.

### Corollaries

By replacing ${\displaystyle \mathbf }$ in the divergence theorem with specific forms, other useful identities can be derived (cf. vector identities).[4]

• With ${\displaystyle \mathbf \rightarrow \mathbf g}$ for a scalar function g and a vector field F,
${\displaystyle \iiint _\left[\mathbf \cdot \left(\nabla g\right)+g\left(\nabla \cdot \mathbf \right)\right]dV=}$ ${\displaystyle \scriptstyle S}$ ${\displaystyle g\mathbf \cdot \mathbf dS.}$
A special case of this is F = ∇ f, in which case the theorem is the basis for Green's identities.
• With ${\displaystyle \mathbf \rightarrow \mathbf \times \mathbf }$ for two vector fields F and G,
${\displaystyle \iiint _\left[\mathbf \cdot \left(\nabla \times \mathbf \right)-\mathbf \cdot \left(\nabla \times \mathbf \right)\right]\,dV=}$ ${\displaystyle \scriptstyle S}$ ${\displaystyle (\mathbf \times \mathbf )\cdot \mathbf d\mathbf .}$
• With ${\displaystyle \mathbf \rightarrow f\mathbf }$ for a scalar function f and vector field c:[5]
${\displaystyle \iiint _\mathbf \cdot \nabla f\,dV=}$ ${\displaystyle \scriptstyle S}$ ${\displaystyle (\mathbf f)\cdot \mathbf dS-\iiint _f(\nabla \cdot \mathbf )\,dV.}$
The last term on the right vanishes for constant ${\displaystyle \mathbf }$ or any divergence free (solenoidal) vector field, e.g. Incompressible flows without sources or sinks such as phase change or chemical reactions etc. In particular, taking ${\displaystyle \mathbf }$ to be constant:
${\displaystyle \iiint _\nabla f\,dV=}$ ${\displaystyle \scriptstyle S}$ ${\displaystyle fd\mathbf .}$
• With ${\displaystyle \mathbf \rightarrow \mathbf \times \mathbf }$ for vector field F and constant vector c:[5]
${\displaystyle \iiint _\mathbf \cdot (\nabla \times \mathbf )\,dV=}$ ${\displaystyle \scriptstyle S}$ ${\displaystyle (\mathbf \times \mathbf )\cdot \mathbf d\mathbf .}$
By reordering the triple product on the right hand side and taking out the constant vector of the integral,
${\displaystyle \iiint _(\nabla \times \mathbf )\,dV\cdot \mathbf =}$ ${\displaystyle \scriptstyle S}$ ${\displaystyle (d\mathbf \times \mathbf )\cdot \mathbf .}$
Hence,
${\displaystyle \iiint _(\nabla \times \mathbf )\,dV=}$ ${\displaystyle \scriptstyle S}$ ${\displaystyle \mathbf \times \mathbf dS.}$

## Example

The vector field corresponding to the example shown. Vectors may point into or out of the sphere.

Suppose we wish to evaluate

${\displaystyle \scriptstyle S}$ ${\displaystyle \mathbf \cdot \mathbf \,dS,}$

where S is the unit sphere defined by

${\displaystyle S=\left\{(x,y,z)\in \mathbb ^\ :\ x^+y^+z^=1\right\},}$

and F is the vector field

${\displaystyle \mathbf =2x\mathbf +y^\mathbf +z^\mathbf .}$

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:

${\displaystyle \iiint _(\nabla \cdot \mathbf )\,dV=2\iiint _(1+y+z)\,dV=2\iiint _dV+2\iiint _y\,dV+2\iiint _z\,dV,}$

where W is the unit ball:

${\displaystyle W=\left\{(x,y,z)\in \mathbb ^\ :\ x^+y^+z^\leq 1\right\}.}$

Since the function y is positive in one hemisphere of W and negative in the other, in an equal and opposite way, its total integral over W is zero. The same is true for z:

${\displaystyle \iiint _y\,dV=\iiint _z\,dV=0.}$

Therefore,

${\displaystyle \scriptstyle S}$ ${\displaystyle \mathbf \cdot \mathbf \,S=2\iiint _\,dV={\frac },}$

because the unit ball W has volume 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

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).[6]

### Inverse-square laws

Any inverse-square law can instead be written in a Gauss's law-type form (with a differential and integral form, as described above). Two examples are Gauss's law (in electrostatics), which follows from the inverse-square Coulomb's law, and Gauss's law for gravity, which follows from the inverse-square Newton's law of universal gravitation. The derivation of the Gauss's 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.[6]

## History

The theorem was first discovered by Lagrange in 1762,[7] then later independently rediscovered by Gauss in 1813,[8] by Ostrogradsky, who also gave the first proof of the general theorem, in 1826,[9] by Green in 1828,[10] Simeon-Denis Poisson in 1824 and Frédéric Sarrus in 1828.[11]

## Examples

To verify the planar variant of the divergence theorem for a region R:

${\displaystyle R=\left\{(x,y)\in \mathbb ^\ :\ x^+y^\leq 1\right\},}$

and the vector field:

${\displaystyle \mathbf (x,y)=2y\mathbf +5x\mathbf .}$

The boundary of R is the unit circle, C, that can be represented parametrically by:

${\displaystyle x=\cos(s),\quad y=\sin(s)}$

such that 0 ≤ s ≤ 2π where s units is the length arc from the point s = 0 to the point P on C. Then a vector equation of C is

${\displaystyle C(s)=\cos(s)\mathbf +\sin(s)\mathbf .}$

At a point P on C:

${\displaystyle P=(\cos(s),\,\sin(s))\,\Rightarrow \,\mathbf =2\sin(s)\mathbf +5\cos(s)\mathbf .}$

Therefore,

${\displaystyle {\begin\oint _\mathbf \cdot \mathbf \,ds&=\int _^(2\sin(s)\mathbf +5\cos(s)\mathbf )\cdot (\cos(s)\mathbf +\sin(s)\mathbf )\,ds\\&=\int _^(2\sin(s)\cos(s)+5\sin(s)\cos(s))\,ds\\&=7\int _^\sin(s)\cos(s)\,ds\\&=0.\end}}$

Because M = 2y, M/x = 0, and because N = 5x, N/y = 0. Thus

${\displaystyle \iint _\,\mathbf {\nabla } \cdot \mathbf \,dA=\iint _\left({\frac {\partial M}{\partial x}}+{\frac {\partial N}{\partial y}}\right)\,dA=0.}$

## Applied Example

Let's say we wanted to evaluate the flux of the following vector field defined by ${\displaystyle \mathbf =2x^+2y^+2z^}$ bounded by the following inequalities:

${\displaystyle \left\\left\{-2\leq y\leq 2\right\}\left\}$

We know from the Divergence Theorem that :${\displaystyle \iiint _\left(\mathbf {\nabla } \cdot \mathbf \right)\,dV=}$ ${\displaystyle \scriptstyle S}$ ${\displaystyle (\mathbf \cdot \mathbf )\,dS.}$

We need to determine ${\displaystyle \nabla \cdot \mathbf }$

The divergence of a three dimensional vector field, ${\displaystyle \mathbf }$, is defined as ${\displaystyle {\frac {\partial \mathbf } }{\partial }}+{\frac {\partial \mathbf } }{\partial }}+{\frac {\partial \mathbf } }{\partial }}}$

Thus, we can set up the following integrals:

${\displaystyle {\scriptstyle S}}$ ${\displaystyle \mathbf \cdot \mathbf \,dS}$

${\displaystyle {\begin&=\iiint _\nabla \cdot \mathbf \,dV\\&=\iiint _{\frac {\partial \mathbf } }{\partial }}+{\frac {\partial \mathbf } }{\partial }}+{\frac {\partial \mathbf } }{\partial }}\,dV\\&=\iiint _4x+4y+4z\,dV\\&=\int _^\int _{-2}^\int _^4x+4y+4z\,dV\end}}$

Now that we have set up the integral, we can evaluate it.

${\displaystyle {\begin&\int _^\int _{-2}^\int _^4x+4y+4z\,dV\\&=\int _{-2}^\int _^12y+12z+18\,dy\,dz\\&=\int _^24(2z+3)\,dz\\&=48\pi \cdot (2\pi +3)\\&\approx \end}}$

## Generalizations

### Multiple dimensions

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

${\displaystyle \underbrace {\int \cdots \int _} _\nabla \cdot \mathbf \,dV=\underbrace {\oint \cdots \oint _{\partial U}} _\mathbf \cdot \mathbf \,dS}$

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

Writing the theorem in Einstein notation:

${\displaystyle \iiint _{\dfrac {\partial \mathbf _}{\partial x_}}dV=}$ ${\displaystyle \scriptstyle S}$ ${\displaystyle \mathbf _n_\,dS}$

suggestively, replacing the vector field F with a rank-n tensor field T, this can be generalized to:[12]

${\displaystyle \iiint _{\dfrac {\partial T_i_\cdots i_\cdots i_}}{\partial x_}}}dV=}$ ${\displaystyle \scriptstyle S}$ ${\displaystyle T_i_\cdots i_\cdots i_}n_}\,dS.}$

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[13]).

## Notes

1. ^ Katz, Victor J. (1979). "The history of Stokes's theorem". Mathematics Magazine. 52 (3): 146–156. doi:10.2307/2690275. JSTOR 2690275. reprinted in Anderson, Marlow (2009). Who Gave You the Epsilon?: And Other Tales of Mathematical History. Mathematical Association of America. pp. 78–79. ISBN 978-0883855690.
2. ^ R. G. Lerner; G. L. Trigg (1994). Encyclopaedia of Physics (2nd ed.). VHC. ISBN 978-3-527-26954-9.
3. ^ Byron, Frederick; Fuller, Robert (1992), Mathematics of Classical and Quantum Physics, Dover Publications, p. 22, ISBN 978-0-486-67164-2
4. ^ a b M. R. Spiegel; S. Lipschutz; D. Spellman (2009). Vector Analysis. Schaum’s Outlines (2nd ed.). USA: McGraw Hill. ISBN 978-0-07-161545-7.
5. ^ a b MathWorld
6. ^ a b C.B. Parker (1994). McGraw Hill Encyclopaedia of Physics (2nd ed.). McGraw Hill. ISBN 978-0-07-051400-3.
7. ^ In his 1762 paper on sound, Lagrange treats a special case of the divergence theorem: Lagrange (1762) "Nouvelles recherches sur la nature et la propagation du son" (New researches on the nature and propagation of sound), Miscellanea Taurinensia (also known as: Mélanges de Turin ), 2: 11 – 172. This article is reprinted as: "Nouvelles recherches sur la nature et la propagation du son" in: J.A. Serret, ed., Oeuvres de Lagrange, (Paris, France: Gauthier-Villars, 1867), vol. 1, pages 151–316; on pages 263–265, Lagrange transforms triple integrals into double integrals using integration by parts.
8. ^ C. F. Gauss (1813) "Theoria attractionis corporum sphaeroidicorum ellipticorum homogeneorum methodo nova tractata," Commentationes societatis regiae scientiarium Gottingensis recentiores, 2: 355–378; Gauss considered a special case of the theorem; see the 4th, 5th, and 6th pages of his article.
9. ^ Mikhail Ostragradsky presented his proof of the divergence theorem to the Paris Academy in 1826; however, his work was not published by the Academy. He returned to St. Petersburg, Russia, where in 1828–1829 he read the work that he'd done in France, to the St. Petersburg Academy, which published his work in abbreviated form in 1831.
• His proof of the divergence theorem – "Démonstration d'un théorème du calcul intégral" (Proof of a theorem in integral calculus) – which he had read to the Paris Academy on February 13, 1826, was translated, in 1965, into Russian together with another article by him. See: Юшкевич А.П. (Yushkevich A.P.) and Антропова В.И. (Antropov V.I.) (1965) "Неопубликованные работы М.В. Остроградского" (Unpublished works of MV Ostrogradskii), Историко-математические исследования (Istoriko-Matematicheskie Issledovaniya / Historical-Mathematical Studies), 16: 49–96; see the section titled: "Остроградский М.В. Доказательство одной теоремы интегрального исчисления" (Ostrogradskii M. V. Dokazatelstvo odnoy teoremy integralnogo ischislenia / Ostragradsky M.V. Proof of a theorem in integral calculus).
• M. Ostrogradsky (presented: November 5, 1828 ; published: 1831) "Première note sur la théorie de la chaleur" (First note on the theory of heat) Mémoires de l'Académie impériale des sciences de St. Pétersbourg, series 6, 1: 129–133; for an abbreviated version of his proof of the divergence theorem, see pages 130–131.
• Victor J. Katz (May1979) "The history of Stokes' theorem," Archived April 2, 2015, at the Wayback Machine Mathematics Magazine, 52(3): 146–156; for Ostragradsky's proof of the divergence theorem, see pages 147–148.
10. ^ George Green, An Essay on the Application of Mathematical Analysis to the Theories of Electricity and Magnetism (Nottingham, England: T. Wheelhouse, 1838). A form of the "divergence theorem" appears on pages 10–12.
11. ^ Other early investigators who used some form of the divergence theorem include:
• Poisson (presented: February 2, 1824 ; published: 1826) "Mémoire sur la théorie du magnétisme" (Memoir on the theory of magnetism), Mémoires de l'Académie des sciences de l'Institut de France, 5: 247–338; on pages 294–296, Poisson transforms a volume integral (which is used to evaluate a quantity Q) into a surface integral. To make this transformation, Poisson follows the same procedure that is used to prove the divergence theorem.
• Frédéric Sarrus (1828) "Mémoire sur les oscillations des corps flottans" (Memoir on the oscillations of floating bodies), Annales de mathématiques pures et appliquées (Nismes), 19: 185–211.
12. ^ K.F. Riley; M.P. Hobson; S.J. Bence (2010). Mathematical methods for physics and engineering. Cambridge University Press. ISBN 978-0-521-86153-3.
13. ^ see for example:
J.A. Wheeler; C. Misner; K.S. Thorne (1973). Gravitation. W.H. Freeman & Co. pp. 85–86, §3.5. ISBN 978-0-7167-0344-0., and
R. Penrose (2007). The Road to Reality. Vintage books. ISBN 978-0-679-77631-4.