Divergence theorem

Template:Short description Template:Redirect Template:Redirect {{#invoke:sidebar|collapsible | class = plainlist | titlestyle = padding-bottom:0.25em; | pretitle = Part of a series of articles about | title = Calculus | image = <math>\int_{a}^{b} f'(t) \, dt = f(b) - f(a)</math> | listtitlestyle = text-align:center; | liststyle = border-top:1px solid #aaa;padding-top:0.15em;border-bottom:1px solid #aaa; | expanded = vector | abovestyle = padding:0.15em 0.25em 0.3em;font-weight:normal; | above =

| list2name = differential | list2titlestyle = display:block;margin-top:0.65em; | list2title = Template:Bigger | list2 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | heading1 = Definitions
 | content1 =

 | heading2 = Concepts
 | content2 =

 | heading3 = Rules and identities
 | content3 =

}}

| list3name = integral | list3title = Template:Bigger | list3 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | content1 =

| heading2 = Definitions

 | content2 =

 | heading3 = Integration by
 | content3 =

}}

| list4name = series | list4title = Template:Bigger | list4 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | content1 =

 | heading2 = Convergence tests
 | content2 =

}}

| list5name = vector | list5title = Template:Bigger | list5 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | content1 =

 | heading2 = Theorems
 | content2 =

}}

| list6name = multivariable | list6title = Template:Bigger | list6 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | heading1 = Formalisms
 | content1 =

 | heading2 = Definitions
 | content2 =

}}

| list7name = advanced | list7title = Template:Bigger | list7 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | content1 =

}}

| list8name = specialized | list8title = Template:Bigger | list8 =

| list9name = miscellanea | list9title = Template:Bigger | list9 =

}}

In vector calculus, the divergence theorem, also known as Gauss's theorem or Ostrogradsky's theorem,<ref name="Katz">Template:Cite journal reprinted in Template:Cite book</ref> is a theorem relating the flux of a vector field through a closed surface to the divergence of the field in the volume enclosed.

More precisely, the divergence theorem states that the surface integral of a vector field over a closed surface, which is called the "flux" through the surface, is equal to the volume integral of the divergence over the region enclosed by the surface. Intuitively, it states that "the sum of all sources of the field in a region (with sinks regarded as negative sources) gives the net flux out of the region".

The divergence theorem is an important result for the mathematics of physics and engineering, particularly in electrostatics and fluid dynamics. In these fields, it 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.

Explanation using liquid flowEdit

Template:See also Vector fields are often illustrated using the example of the velocity field of a fluid, such as a gas or liquid. A moving liquid has a velocity—a speed and a direction—at each point, which can be represented by a vector, so that the velocity of the liquid at any moment forms a vector field. Consider an imaginary closed surface S inside a body of liquid, enclosing a volume of liquid. The flux of liquid out of the volume at any time is equal to the volume rate of fluid crossing this surface, i.e., the surface integral of the velocity over the surface.

Since liquids are incompressible, the amount of liquid inside a closed volume is constant; if there are no sources or sinks inside the volume then the flux of liquid out of S is zero. If the liquid is moving, it may flow into the volume at some points on the surface S and out of the volume at other points, but the amounts flowing in and out at any moment are equal, so the net flux of liquid out of the volume is zero.

However, if a source of liquid is inside the closed surface, such as a pipe through which liquid is introduced, the additional liquid will exert pressure on the surrounding liquid, causing an outward flow in all directions. This will cause a net outward flow through the surface S. The flux outward through S equals the volume rate of flow of fluid into S from the pipe. Similarly if there is a sink or drain inside S, such as a pipe which drains the liquid off, the external pressure of the liquid will cause a velocity throughout the liquid directed inward toward the location of the drain. The volume rate of flow of liquid inward through the surface S equals the rate of liquid removed by the sink.

If there are multiple sources and sinks of liquid inside S, the flux through the surface can be calculated by adding up the volume rate of liquid added by the sources and subtracting the rate of liquid drained off by the sinks. The volume rate of flow of liquid through a source or sink (with the flow through a sink given a negative sign) is equal to the divergence of the velocity field at the pipe mouth, so adding up (integrating) the divergence of the liquid throughout the volume enclosed by S equals the volume rate of flux through S. This is the divergence theorem.<ref>Template:Cite book</ref>

The divergence theorem is employed in any conservation law which states that the total volume of all sinks and sources, that is the volume integral of the divergence, is equal to the net flow across the volume's boundary.<ref>Template:Citation</ref>

Mathematical statementEdit

File:Divergence theorem.svg

A region Template:Mvar bounded by the surface Template:Mvar = ∂Template:Mvar with the surface normal Template:Mvar

Suppose Template:Mvar is a subset of <math>\mathbb{R}^n</math> (in the case of Template:Math represents a volume in three-dimensional space) which is compact and has a piecewise smooth boundary Template:Mvar (also indicated with <math>\partial V = S</math>). If Template:Math is a continuously differentiable vector field defined on a neighborhood of Template:Mvar, then:<ref name="Wiley">Template:Cite book</ref><ref name="Kreyszig">Template:Cite book</ref>

Template:Oiint)\,\mathrm{d}S .</math>

}} The left side is a volume integral over the volume Template:Mvar, and the right side is the surface integral over the boundary of the volume Template:Mvar. The closed, measurable set <math>\partial V</math> is oriented by outward-pointing normals, and <math>\mathbf{\hat{n}}</math> is the outward pointing unit normal at almost each point on the boundary <math>\partial V</math>. (<math>\mathrm{d} \mathbf{S}</math> may be used as a shorthand for <math>\mathbf{n} \mathrm{d} S</math>.) 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 Template:Mvar.

Informal derivationEdit

The divergence theorem follows from the fact that if a volume Template:Mvar is partitioned into separate parts, the flux out of the original volume is equal to the algebraic sum of the flux out of each component volume.<ref name="Benford">{{#invoke:citation/CS1|citation |CitationClass=web }}</ref><ref name="Purcell">Template:Cite book</ref> This is true despite the fact that the new subvolumes have surfaces that were not part of the original volume's surface, because these surfaces are just partitions between two of the subvolumes and the flux through them just passes from one volume to the other and so cancels out when the flux out of the subvolumes is summed.

File:Divergence theorem 1 - split volume.png

A volume divided into two subvolumes. At right the two subvolumes are separated to show the flux out of the different surfaces.

See the diagram. A closed, bounded volume Template:Mvar is divided into two volumes Template:Math and Template:Math by a surface Template:Math (green). The flux Template:Math out of each component region Template:Math is equal to the sum of the flux through its two faces, so the sum of the flux out of the two parts is

where Template:Math and Template:Math are the flux out of surfaces Template:Math and Template:Math, Template:Math is the flux through Template:Math out of volume 1, and Template:Math is the flux through Template:Math out of volume 2. The point is that surface Template:Math is part of the surface of both volumes. The "outward" direction of the normal vector <math>\mathbf{\hat n}</math> is opposite for each volume, so the flux out of one through Template:Math is equal to the negative of the flux out of the other so these two fluxes cancel in the sum.

<math>\Phi_\text{31} = \iint_{S_3} \mathbf{F} \cdot \mathbf{\hat n} \; \mathrm{d}S = -\iint_{S_3} \mathbf{F} \cdot (-\mathbf{\hat n}) \; \mathrm{d}S = -\Phi_\text{32}</math>

Therefore:

Since the union of surfaces Template:Math and Template:Math is Template:Mvar

File:Divergence theorem 2 - volume partition.png

The volume can be divided into any number of subvolumes and the flux out of V is equal to the sum of the flux out of each subvolume, because the flux through the green surfaces cancels out in the sum. In (b) the volumes are shown separated slightly, illustrating that each green partition is part of the boundary of two adjacent volumes

This principle applies to a volume divided into any number of parts, as shown in the diagram.<ref name="Purcell" /> Since the integral over each internal partition (green surfaces) appears with opposite signs in the flux of the two adjacent volumes they cancel out, and the only contribution to the flux is the integral over the external surfaces (grey). Since the external surfaces of all the component volumes equal the original surface.

<math>\Phi(V) = \sum_{V_\text{i}\subset V} \Phi(V_\text{i})</math>

File:Divergence theorem 3 - infinitesimals.png

</math> approaches <math>\operatorname{div} \mathbf{F}</math>

The flux Template:Math out of each volume is the surface integral of the vector field Template:Math over the surface

<math>\iint_{S(V)} \mathbf{F} \cdot \mathbf{\hat n} \; \mathrm{d}S = \sum_{V_\text{i}\subset V} \iint_{S(V_\text{i})} \mathbf{F} \cdot \mathbf{\hat n} \; \mathrm{d}S</math>

The goal is to divide the original volume into infinitely many infinitesimal volumes. As the volume is divided into smaller and smaller parts, the surface integral on the right, the flux out of each subvolume, approaches zero because the surface area Template:Math approaches zero. However, from the definition of divergence, the ratio of flux to volume, <math>\frac{\Phi(V_\text{i})}{|V_\text{i}|} = \frac{1}{|V_\text{i}|} \iint_{S(V_\text{i})} \mathbf{F} \cdot \mathbf{\hat n} \; \mathrm{d}S</math>, the part in parentheses below, does not in general vanish but approaches the divergence Template:Math as the volume approaches zero.<ref name="Purcell" />

<math>\iint_{S(V)} \mathbf{F} \cdot \mathbf{\hat n} \; \mathrm{d}S = \sum_{V_\text{i} \subset V} \left(\frac{1}{|V_\text{i}|} \iint_{S(V_\text{i})} \mathbf{F} \cdot \mathbf{\hat n} \; \mathrm{d}S\right) |V_\text{i}|</math>

As long as the vector field Template:Math has continuous derivatives, the sum above holds even in the limit when the volume is divided into infinitely small increments

<math>\iint_{S(V)} \mathbf{F} \cdot \mathbf{\hat n} \; \mathrm{d}S = \lim_{|V_\text{i}|\to 0}\sum_{V_\text{i}\subset V} \left(\frac{1}{|V_\text{i}|}\iint_{S(V_\text{i})} \mathbf{F} \cdot \mathbf{\hat n} \; \mathrm{d}S\right) |V_\text{i}|</math>

As <math>|V_\text{i}|</math> approaches zero volume, it becomes the infinitesimal Template:Math, the part in parentheses becomes the divergence, and the sum becomes a volume integral over Template:Mvar Template:Equation box 1 Since this derivation is coordinate free, it shows that the divergence does not depend on the coordinates used.

ProofsEdit

For bounded open subsets of Euclidean spaceEdit

We are going to prove the following:Template:Citation needed Template:Math theorem

Proof of Theorem. <ref name="Alt 2016 p. ">Template:Cite book</ref> Template:Ordered list \int_{-\infty}^{\infty}u_{x_i}(x)\,dx_i\,dx' = 0</math> by the fundamental theorem of calculus, and <math>\int_{\partial \Omega}u\nu_i\,dS = 0</math> since <math>u</math> vanishes on a neighborhood of <math>\partial \Omega</math>. Thus the theorem holds for <math>u</math> with compact support in <math>\Omega</math>. Thus we have reduced to the case where <math>u</math> has compact support in some <math>U_j</math>.

| 3 = So assume <math>u</math> has compact support in some <math>U_j</math>. The last step now is to show that the theorem is true by direct computation. Change notation to <math>U = U_j</math>, and bring in the notation from (2) used to describe <math>U</math>. Note that this means that we have rotated and translated <math>\Omega</math>. This is a valid reduction since the theorem is invariant under rotations and translations of coordinates. Since <math>u(x) = 0</math> for <math>|x'| \geq r</math> and for <math>|x_n - g(x')| \geq h</math>, we have for each <math>i \in \{1, \dots, n\}</math> that <math display="block"> \begin{align}

   \int_{\Omega}u_{x_i}\,dV &= \int_{|x'| < r}\int_{g(x') - h}^{g(x')}u_{x_i}(x', x_n)\,dx_n\,dx' \\
   &= \int_{\mathbb{R}^{n - 1}}\int_{-\infty}^{g(x')}u_{x_i}(x', x_n)\,dx_n\,dx'.

\end{align} </math> For <math>i = n</math> we have by the fundamental theorem of calculus that <math display="block">\int_{\mathbb{R}^{n - 1}}\int_{-\infty}^{g(x')}u_{x_n}(x', x_n)\,dx_n\,dx' = \int_{\mathbb{R}^{n - 1}}u(x', g(x'))\,dx'.</math> Now fix <math>i \in \{1, \dots, n - 1\}</math>. Note that <math display="block">\int_{\mathbb{R}^{n - 1}}\int_{-\infty}^{g(x')}u_{x_i}(x', x_n)\,dx_n\,dx' = \int_{\mathbb{R}^{n - 1}}\int_{-\infty}^{0}u_{x_i}(x', g(x') + s)\,ds\,dx'</math> Define <math>v : \mathbb{R}^{n} \to \mathbb{R}</math> by <math>v(x', s) = u(x', g(x') + s)</math>. By the chain rule, <math display="block">v_{x_i}(x', s) = u_{x_i}(x', g(x') + s) + u_{x_n}(x', g(x') + s)g_{x_i}(x').</math> But since <math>v</math> has compact support, we can integrate out <math>dx_i</math> first to deduce that <math display="block">\int_{\mathbb{R}^{n - 1}}\int_{-\infty}^{0}v_{x_i}(x', s)\,ds\,dx' = 0.</math> Thus <math display="block"> \begin{align}

   \int_{\mathbb{R}^{n - 1}}\int_{-\infty}^{0}u_{x_i}(x', g(x') + s)\,ds\,dx' &= \int_{\mathbb{R}^{n - 1}}\int_{-\infty}^{0}-u_{x_n}(x', g(x') + s)g_{x_i}(x')\,ds\,dx' \\
   &= \int_{\mathbb{R}^{n - 1}}-u(x', g(x'))g_{x_i}(x')\,dx'.

\end{align} </math> In summary, with <math>\nabla u = (u_{x_1}, \dots, u_{x_n})</math> we have <math display="block">\int_{\Omega}\nabla u\,dV = \int_{\mathbb{R}^{n - 1}}\int_{-\infty}^{g(x')}\nabla u\,dV = \int_{\mathbb{R}^{n - 1}}u(x', g(x'))(-\nabla g(x'), 1)\,dx'.</math> Recall that the outward unit normal to the graph <math>\Gamma</math> of <math>g</math> at a point <math>(x', g(x')) \in \Gamma</math> is <math>\nu(x', g(x')) = \frac{1}{\sqrt{1 + |\nabla g(x')|^2}}(-\nabla g(x'), 1)</math> and that the surface element <math>dS</math> is given by <math display="inline">dS = \sqrt{1 + |\nabla g(x')|^2}\,dx'</math>. Thus <math display="block">\int_{\Omega}\nabla u\,dV = \int_{\partial \Omega}u\nu\,dS.</math> This completes the proof. }}

For compact Riemannian manifolds with boundaryEdit

We are going to prove the following:Template:Citation needed Template:Math theorem

Proof of Theorem. <ref name="Taylor 2011 p. "> Template:Cite book </ref> We use the Einstein summation convention. By using a partition of unity, we may assume that <math>u</math> and <math>X</math> have compact support in a coordinate patch <math>O \subset \overline{\Omega}</math>. First consider the case where the patch is disjoint from <math>\partial \Omega</math>. Then <math>O</math> is identified with an open subset of <math>\mathbb{R}^n</math> and integration by parts produces no boundary terms: <math display="block"> \begin{align}

   (\operatorname{grad} u, X) &= \int_{O}\langle \operatorname{grad} u, X \rangle \sqrt{g}\,dx \\
   &= \int_{O}\partial_j u X^j \sqrt{g}\,dx \\
   &= -\int_{O}u \partial_j(\sqrt{g}X^j)\,dx \\
   &= -\int_{O} u \frac{1}{\sqrt{g}}\partial_j(\sqrt{g}X^j)\sqrt{g}\,dx \\
   &= (u, -\frac{1}{\sqrt{g}}\partial_j(\sqrt{g}X^j)) \\
   &= (u, -\operatorname{div} X).

\end{align} </math> In the last equality we used the Voss-Weyl coordinate formula for the divergence, although the preceding identity could be used to define <math>-\operatorname{div}</math> as the formal adjoint of <math>\operatorname{grad}</math>. Now suppose <math>O</math> intersects <math>\partial \Omega</math>. Then <math>O</math> is identified with an open set in <math>\mathbb{R}_{+}^n = \{x \in \mathbb{R}^n : x_n \geq 0\}</math>. We zero extend <math>u</math> and <math>X</math> to <math>\mathbb{R}_+^n</math> and perform integration by parts to obtain <math display="block"> \begin{align}

   (\operatorname{grad} u, X) &= \int_{O}\langle \operatorname{grad} u, X \rangle \sqrt{g}\,dx \\
   &= \int_{\mathbb{R}_+^n}\partial_j u X^j \sqrt{g}\,dx \\
   &= (u, -\operatorname{div} X) - \int_{\mathbb{R}^{n - 1}}u(x', 0)X^n(x', 0)\sqrt{g(x', 0)}\,dx',

\end{align} </math> where <math>dx' = dx_1 \dots dx_{n - 1}</math>. By a variant of the straightening theorem for vector fields, we may choose <math>O</math> so that <math>\frac{\partial}{\partial x_n}</math> is the inward unit normal <math>-N</math> at <math>\partial \Omega</math>. In this case <math>\sqrt{g(x', 0)}\,dx' = \sqrt{g_{\partial \Omega}(x')}\,dx' = dS</math> is the volume element on <math>\partial \Omega</math> and the above formula reads <math display="block"> (\operatorname{grad} u, X) = (u, -\operatorname{div} X) + \int_{\partial \Omega}u\langle X, N \rangle \,dS. </math> This completes the proof.

CorollariesEdit

By replacing Template:Math in the divergence theorem with specific forms, other useful identities can be derived (cf. vector identities).<ref name=spiegel>Template:Cite book</ref>

With <math>\mathbf{F}\rightarrow \mathbf{F}g</math> for a scalar function Template:Mvar and a vector field Template:Math,

Template:Oiint

A special case of this is <math>\mathbf{F} = \nabla f</math>, in which case the theorem is the basis for Green's identities.

With <math>\mathbf{F}\rightarrow \mathbf{F}\times \mathbf{G}</math> for two vector fields Template:Math and Template:Math, where <math>\times</math> denotes a cross product,

Template:Oiint

With <math>\mathbf{F}\rightarrow \mathbf{F}\cdot \mathbf{G}</math> for two vector fields Template:Math and Template:Math, where <math>\cdot </math> denotes a dot product,

Template:Oiint

With <math>\mathbf{F}\rightarrow f\mathbf{c}</math> for a scalar function Template:Math and vector field c:<ref name=mathworld>MathWorld</ref>

Template:Oiint

The last term on the right vanishes for constant <math>\mathbf{c}</math> 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 <math>\mathbf{c}</math> to be constant:

Template:Oiint

With <math>\mathbf{F}\rightarrow \mathbf{c}\times\mathbf{F}</math> for vector field Template:Math and constant vector c:<ref name=mathworld/>

Template:Oiint

By reordering the triple product on the right hand side and taking out the constant vector of the integral,

Template:Oiint

Hence,

Template:Oiint

ExampleEdit

File:Vector Field on a Sphere.png

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

File:SurfacesWithAndWithoutBoundary.svg

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 we wish to evaluate

Template:Oiint

where Template:Mvar is the unit sphere defined by

<math>S = \left \{ (x,y, z) \in \mathbb{R}^3 \ : \ x^2+y^2+z^2 = 1 \right \},</math>

and Template:Math is the vector field

<math>\mathbf{F} = 2x\mathbf{i}+y^2\mathbf{j}+z^2\mathbf{k}.</math>

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:

<math>\iiint_W (\nabla \cdot \mathbf{F})\,\mathrm{d}V = 2\iiint_W (1 + y + z)\, \mathrm{d}V = 2\iiint_W \mathrm{d}V + 2\iiint_W y\, \mathrm{d}V + 2\iiint_W z\, \mathrm{d}V,</math>

where Template:Mvar is the unit ball:

<math>W = \left \{ (x,y, z) \in \mathbb{R}^3 \ : \ x^2+y^2+z^2\leq 1 \right \}.</math>

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:

<math>\iiint_W y\, \mathrm{d}V = \iiint_W z\, \mathrm{d}V = 0.</math>

Therefore,

Template:Oiint

because the unit ball Template:Mvar has volume Template:Math.

ApplicationsEdit

Differential and integral forms of physical lawsEdit

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 equationsEdit

{{#invoke:Labelled list hatnote|labelledList|Main article|Main articles|Main page|Main pages}}

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).<ref name="C.B. Parker 1994">Template:Cite book</ref>

Inverse-square lawsEdit

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.<ref name="C.B. Parker 1994"/>

HistoryEdit

Joseph-Louis Lagrange introduced the notion of surface integrals in 1760 and again in more general terms in 1811, in the second edition of his Mécanique Analytique. Lagrange employed surface integrals in his work on fluid mechanics.<ref name=":0">Template:Cite book</ref> He discovered the divergence theorem in 1762.<ref>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.</ref>

Carl Friedrich Gauss was also using surface integrals while working on the gravitational attraction of an elliptical spheroid in 1813, when he proved special cases of the divergence theorem.<ref>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.</ref><ref name=":0" /> He proved additional special cases in 1833 and 1839.<ref name=":32">Template:Cite journal</ref> But it was Mikhail Ostrogradsky, who gave the first proof of the general theorem, in 1826, as part of his investigation of heat flow.<ref>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," Template:Webarchive Mathematics Magazine, 52(3): 146–156 {{#invoke:doi|main}}; for Ostragradsky's proof of the divergence theorem, see pages 147–148.</ref> Special cases were proven by George Green in 1828 in An Essay on the Application of Mathematical Analysis to the Theories of Electricity and Magnetism,<ref>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.</ref><ref name=":32" /> Siméon Denis Poisson in 1824 in a paper on elasticity, and Frédéric Sarrus in 1828 in his work on floating bodies.<ref>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.</ref><ref name=":32" />

Worked examplesEdit

Example 1Edit

To verify the planar variant of the divergence theorem for a region <math>R</math>:

<math>R = \left \{ (x, y) \in \mathbb{R}^2 \ : \ x^2 + y^2 \leq 1 \right \},</math>

and the vector field:

<math> \mathbf{F}(x,y)= 2 y\mathbf{i} + 5x \mathbf{j}.</math>

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

such that <math>0 \leq s \leq 2\pi</math> where <math>s</math> units is the length arc from the point <math>s = 0</math> to the point <math>P</math> on <math>C</math>. Then a vector equation of <math>C</math> is

<math>C(s) = \cos(s)\mathbf{i} + \sin(s)\mathbf{j}.</math>

At a point <math>P</math> on <math>C</math>:

<math> P = (\cos(s),\, \sin(s)) \, \Rightarrow \, \mathbf{F} = 2\sin(s)\mathbf{i} + 5\cos(s)\mathbf{j}.</math>

Therefore,

<math>\begin{align}

\oint_C \mathbf{F} \cdot \mathbf{n}\, \mathrm{d}s &= \int_0^{2\pi} (2 \sin(s) \mathbf{i} + 5 \cos(s) \mathbf{j}) \cdot (\cos(s) \mathbf{i} + \sin(s) \mathbf{j})\, \mathrm{d}s\\ &= \int_0^{2\pi} (2 \sin(s) \cos(s) + 5 \sin(s) \cos(s))\, \mathrm{d}s\\ &= 7\int_0^{2\pi} \sin(s) \cos(s)\, \mathrm{d}s\\ &= 0. \end{align}</math>

Because <math>M = \mathfrak{Re}(\mathbf{F}) = 2y</math>, we can evaluate Template:Nowrap and because Template:Nowrap <math>\frac{\partial N}{\partial y} = 0</math>. Thus

<math>\iint_R \, \mathbf{\nabla}\cdot\mathbf{F} \, \mathrm{d}A = \iint_R \left (\frac{\partial M}{\partial x} + \frac{\partial N}{\partial y} \right) \, \mathrm{d}A = 0. </math>

Example 2Edit

Let's say we wanted to evaluate the flux of the following vector field defined by <math> \mathbf{F}=2x^2 \textbf{i} +2y^2 \textbf{j} +2z^2\textbf{k} </math> bounded by the following inequalities:

<math>\left\{0\le x \le 3\right\}, \left\{-2\le y \le 2\right\}, \left\{0\le z \le 2\pi\right\}</math>

By the divergence theorem,

Template:Oiint

We now need to determine the divergence of <math>\textbf{F}</math>. If <math>\mathbf{F}</math> is a three-dimensional vector field, then the divergence of <math>\textbf{F}</math> is given by <math display="inline">\nabla \cdot \textbf{F} = \left( \frac{\partial}{\partial x}\textbf{i} + \frac{\partial}{\partial y}\textbf{j} + \frac{\partial}{\partial z}\textbf{k} \right) \cdot \textbf{F}</math>.

Thus, we can set up the following flux integral Template:Oiint as follows:

<math>

\begin{align} I &=\iiint_V \nabla \cdot \mathbf{F} \, \mathrm{d}V\\[6pt] &=\iiint_V \left( \frac{\partial\mathbf{F_x}}{\partial x}+\frac{\partial\mathbf{F_y}}{\partial y}+\frac{\partial\mathbf{F_z}}{\partial z} \right) \mathrm{d}V\\[6pt] &=\iiint_V (4x+4y+4z) \, \mathrm{d}V\\[6pt] &=\int_0^3 \int_{-2}^2 \int_0^{2\pi} (4x+4y+4z) \, \mathrm{d}V \end{align} </math>

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

<math>\begin{align}

\int_0^3 \int_{-2}^2 \int_0^{2\pi} (4x+4y+4z) \, \mathrm{d}V &=\int_{-2}^2 \int_0^{2\pi} (12y+12z+18) \, \mathrm{d}y \, \mathrm{d}z\\[6pt] &=\int_0^{2\pi} 24 (2z+3)\, \mathrm{d}z\\[6pt] &=48\pi(2\pi+3) \end{align}

</math>

GeneralizationsEdit

Multiple dimensionsEdit

One can use the generalised Stokes' theorem to equate the Template:Mvar-dimensional volume integral of the divergence of a vector field Template:Math over a region Template:Mvar to the Template:Math-dimensional surface integral of Template:Math over the boundary of Template:Mvar:

<math> \underbrace{ \int \cdots \int_U }_n \nabla \cdot \mathbf{F} \, \mathrm{d}V = \underbrace{ \oint_{} \cdots \oint_{\partial U} }_{n-1} \mathbf{F} \cdot \mathbf{n} \, \mathrm{d}S </math>

This equation is also known as the divergence theorem.

When Template:Math, this is equivalent to Green's theorem.

When Template:Math, it reduces to the fundamental theorem of calculus, part 2.

Tensor fieldsEdit

{{#invoke:Labelled list hatnote|labelledList|Main article|Main articles|Main page|Main pages}}

Writing the theorem in Einstein notation:

Template:Oiint

suggestively, replacing the vector field Template:Math with a rank-Template:Mvar tensor field Template:Mvar, this can be generalized to:<ref>Template:Cite book</ref>

Template:Oiint{\partial x_{i_q}} \mathrm{d}V=</math>

| intsubscpt = <math>\scriptstyle S</math> | integrand = <math>T_{i_1i_2\cdots i_q\cdots i_n}n_{i_q}\, \mathrm{d}S .</math> }}

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<ref>see for example:
Template:Cite book, and
Template:Cite book</ref>).

ReferencesEdit

Template:Reflist

External linksEdit

Template:Sister project

Template:Springer
Differential Operators and the Divergence Theorem at MathPages
The Divergence (Gauss) Theorem by Nick Bykov, Wolfram Demonstrations Project.
{{#invoke:Template wrapper|{{#if:|list|wrap}}|_template=cite web

|_exclude=urlname, _debug, id |url = https://mathworld.wolfram.com/{{#if:DivergenceTheorem%7CDivergenceTheorem.html}} |title = Divergence Theorem |author = Weisstein, Eric W. |website = MathWorld |access-date = |ref = Template:SfnRef }} – This article was originally based on the GFDL article from PlanetMath at https://web.archive.org/web/20021029094728/http://planetmath.org/encyclopedia/Divergence.html

Template:Calculus topics

Divergence theorem

Contents

Explanation using liquid flowEdit

Mathematical statementEdit

Informal derivationEdit

ProofsEdit

For bounded open subsets of Euclidean spaceEdit

For compact Riemannian manifolds with boundaryEdit

CorollariesEdit

ExampleEdit

ApplicationsEdit

Differential and integral forms of physical lawsEdit

Continuity equationsEdit

Inverse-square lawsEdit

HistoryEdit

Worked examplesEdit

Example 1Edit

Example 2Edit

GeneralizationsEdit

Multiple dimensionsEdit

Tensor fieldsEdit

See alsoEdit

ReferencesEdit

External linksEdit