Stokes' theorem

For the equation governing viscous drag in fluids, see Stokes' law.

In vector calculus, and more generally differential geometry, Stokes' theorem (also called the generalized Stokes' theorem) is a statement about the integration of differential forms on manifolds, which both simplifies and generalizes several theorems from vector calculus. Stokes' theorem says that the integral of a differential form ω over the boundary of some orientable manifold Ω is equal to the integral of its exterior derivative dω over the whole of Ω, i.e.,

\int_{\partial \Omega}\omega=\int_\Omega \mathrm {d}\omega.

This modern form of Stokes' theorem is a vast generalization of a classical result. Lord Kelvin communicated it to George Stokes in a letter dated July 2, 1850.[1][2][3] Stokes set the theorem as a question on the 1854 Smith's Prize exam, which led to the result bearing his name, even though it was actually first published by Hermann Hankel in 1861.[3][4] This classical Kelvin–Stokes theorem relates the surface integral of the curl of a vector field F over a surface Σ in Euclidean three-space to the line integral of the vector field over its boundary ∂Σ:

 \iint_{\Sigma} \nabla \times \mathbf{F} \cdot \mathrm{d}\mathbf{\Sigma} = \oint_{\partial\Sigma} \mathbf{F} \cdot \mathrm{d} \mathbf{r}.

This classical statement, along with the classical divergence theorem, fundamental theorem of calculus, and Green's theorem are simply special cases of the general formulation stated above.

Introduction

The fundamental theorem of calculus states that the integral of a function f over the interval [a, b] can be calculated by finding an antiderivative F of f :

\int_a^b f(x)\,\mathrm dx = F(b) - F(a).

Stokes' theorem is a vast generalization of this theorem in the following sense.

In even simpler terms, one can consider that points can be thought of as the boundaries of curves, that is as 0-dimensional boundaries of 1-dimensional manifolds. So, just as one can find the value of an integral (f dx = dF) over a 1-dimensional manifolds ([a, b]) by considering the anti-derivative (F) at the 0-dimensional boundaries ([a, b]), one can generalize the fundamental theorem of calculus, with a few additional caveats, to deal with the value of integrals (dω) over n-dimensional manifolds (Ω) by considering the anti-derivative (ω) at the (n  1)-dimensional boundaries (dΩ) of the manifold.

So the fundamental theorem reads:

\int_{[a, b]} f(x)\,\mathrm{d}x = \int_{[a, b]} \mathrm{d}F = \int_{\{a\}^- \cup \{b\}^+} F = F(b) - F(a).

General formulation

Let Ω be an oriented smooth manifold of dimension n and let α be an n-differential form that is compactly supported on Ω. First, suppose that α is compactly supported in the domain of a single, oriented coordinate chart {U, φ}. In this case, we define the integral of α over Ω as

\int_\Omega \alpha = \int_{\phi(U)} \left(\phi^{-1}\right)^* \alpha,

i.e., via the pullback of α to Rn.

More generally, the integral of α over Ω is defined as follows: Let {ψi} be a partition of unity associated with a locally finite cover {Ui, φi} of (consistently oriented) coordinate charts, then define the integral

\int_\Omega \alpha \equiv \sum_i \int_{U_i} \psi_i \alpha,

where each term in the sum is evaluated by pulling back to Rn as described above. This quantity is well-defined; that is, it does not depend on the choice of the coordinate charts, nor the partition of unity.

Stokes' theorem reads: If ω is an (n  1)-form with compact support on Ω and ∂Ω denotes the boundary of Ω with its induced orientation, then

A "normal" integration manifold (here called D instead of Ω) for the special case n = 2
\int_\Omega \mathrm{d}\omega = \int_{\partial \Omega} \omega \ \ \left( = \oint_{\partial \Omega} \omega\right).

Here d is the exterior derivative, which is defined using the manifold structure only. On the r.h.s., a circle is sometimes used within the integral sign to stress the fact that the (n  1)-manifold ∂Ω is closed.[5] The r.h.s. of the equation is often used to formulate integral laws; the l.h.s. then leads to equivalent differential formulations (see below).

The theorem is often used in situations where Ω is an embedded oriented submanifold of some bigger manifold on which the form ω is defined.

A proof becomes particularly simple if the submanifold Ω is a so-called "normal manifold", as in the figure on the r.h.s., which can be segmented into vertical stripes (e.g. parallel to the xn direction), such that after a partial integration concerning this variable, nontrivial contributions come only from the upper and lower boundary surfaces (coloured in yellow and red, respectively), where the complementary mutual orientations are visible through the arrows.

Topological preliminaries; integration over chains

Let M be a smooth manifold. A smooth singular k-simplex of M is a smooth map from the standard simplex in Rk to M. The free abelian group, Sk, generated by singular k-simplices is said to consist of singular k-chains of M. These groups, together with the boundary map, ∂, define a chain complex. The corresponding homology (resp. cohomology) is called the smooth singular homology (resp. cohomology) of M.

On the other hand, the differential forms, with exterior derivative, d, as the connecting map, form a cochain complex, which defines de Rham cohomology.

Differential k-forms can be integrated over a k-simplex in a natural way, by pulling back to Rk. Extending by linearity allows one to integrate over chains. This gives a linear map from the space of k-forms to the k-th group in the singular cochain, Sk*, the linear functionals on Sk. In other words, a k-form ω defines a functional

I(\omega)(c) = \oint_c \omega

on the k-chains. Stokes' theorem says that this is a chain map from de Rham cohomology to singular cohomology; the exterior derivative, d, behaves like the dual of ∂ on forms. This gives a homomorphism from de Rham cohomology to singular cohomology. On the level of forms, this means:

  1. closed forms, i.e., dω = 0, have zero integral over boundaries, i.e. over manifolds that can be written as \partial \,\sum\nolimits_c M_c, and
  2. exact forms, i.e., ω = dσ, have zero integral over cycles, i.e. if the boundaries sum up to the empty set: \sum\nolimits_c\partial M_c=\varnothing.

De Rham's theorem shows that this homomorphism is in fact an isomorphism. So the converse to 1 and 2 above hold true. In other words, if {ci} are cycles generating the k-th homology group, then for any corresponding real numbers, {ai}, there exist a closed form, ω, such that

\oint_{c_i} \omega = a_i,

and this form is unique up to exact forms.

Underlying principle

To simplify these topological arguments, it is worthwhile to examine the underlying principle by considering an example for d = 2 dimensions. The essential idea can be understood by the diagram on the left, which shows that, in an oriented tiling of a manifold, the interior paths are traversed in opposite directions; their contributions to the path integral thus cancel each other pairwise. As a consequence, only the contribution from the boundary remains. It thus suffices to prove Stokes' theorem for sufficiently fine tilings (or, equivalently, simplices), which usually is not difficult.

Special cases

The general form of the Stokes theorem using differential forms is more powerful and easier to use than the special cases. The traditional versions can be formulated using Cartesian coordinates without the machinery of differential geometry, and thus are more accessible. Further, they are older and their names are more familiar as a result. The traditional forms are often considered more convenient by practicing scientists and engineers but the non-naturalness of the traditional formulation becomes apparent when using other coordinate systems, even familiar ones like spherical or cylindrical coordinates. There is potential for confusion in the way names are applied, and the use of dual formulations.

Kelvin–Stokes theorem

An illustration of the Kelvin–Stokes theorem, with surface Σ, its boundary ∂Σ and the "normal" vector n.

This is a (dualized) 1+1 dimensional case, for a 1-form (dualized because it is a statement about vector fields). This special case is often just referred to as the Stokes' theorem in many introductory university vector calculus courses and as used in physics and engineering. It is also sometimes known as the curl theorem.

The classical Kelvin–Stokes theorem:

 \int_\Sigma \nabla \times \mathbf{F} \cdot \mathrm{d}\mathbf{\Sigma} = \oint_{\partial\Sigma} \mathbf{F} \cdot \mathrm{d} \mathbf{r},

which relates the surface integral of the curl of a vector field over a surface Σ in Euclidean three-space to the line integral of the vector field over its boundary, is a special case of the general Stokes theorem (with n = 2) once we identify a vector field with a 1 form using the metric on Euclidean three-space. The curve of the line integral, ∂Σ, must have positive orientation, meaning that dr points counterclockwise when the surface normal, dΣ, points toward the viewer, following the right-hand rule.

One consequence of the Kelvin–Stokes theorem is that the field lines of a vector field with zero curl cannot be closed contours. The formula can be rewritten as:

\iint_{\Sigma}\left\{\left(\frac{\partial R}{\partial y}-\frac{\partial Q}{\partial z}\right)\mathrm{d}y\mathrm{d}z +\left(\frac{\partial P}{\partial z}-\frac{\partial R}{\partial x}\right)\mathrm{d}z\mathrm{d}x  +\left (\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)\mathrm{d}x\mathrm{d}y\right\}
=\oint_{\partial\Sigma}\left \{P\mathrm{d}x+Q\mathrm{d}y+R\mathrm{d}z\right\},

where P, Q and R are the components of F.

These variants are rarely used:

 \int_{\Sigma} \left [ g \left(\nabla \times \mathbf{F} \right) + \left( \nabla g \right) \times \mathbf{F} \right ] \cdot \mathrm{d}\mathbf{\Sigma} = \oint_{\partial\Sigma} g \mathbf{F} \cdot \mathrm{d} \mathbf{r},
 \int_{\Sigma} \left[ \mathbf{F} \left(\nabla \cdot \mathbf{G} \right) - \mathbf{G}\left(\nabla \cdot \mathbf{F} \right) + \left( \mathbf{G} \cdot \nabla \right) \mathbf{F} - \left(\mathbf{F} \cdot \nabla \right) \mathbf{G} \right ] \cdot \mathrm{d}\mathbf{\Sigma} = \oint_{\partial\Sigma} \left( \mathbf{F} \times \mathbf{G}\right) \cdot \mathrm{d} \mathbf{r}.
 \int_{\Sigma} \nabla (\mathbf{F}\cdot \mathrm{d}\mathbf{\Sigma}) -(\nabla\cdot \mathbf{F}) \mathrm{d}\mathbf{\Sigma} = \oint_{\partial\Sigma}  \mathrm{d} \mathbf{r}\times\mathbf{F}.

Green's theorem

Green's theorem is immediately recognizable as the third integrand of both sides in the integral in terms of P, Q, and R cited above.

In electromagnetism

Two of the four Maxwell equations involve curls of 3-D vector fields and their differential and integral forms are related by the Kelvin–Stokes theorem. Caution must be taken to avoid cases with moving boundaries: the partial time derivatives are intended to exclude such cases. If moving boundaries are included, interchange of integration and differentiation introduces terms related to boundary motion not included in the results below (see Differentiation under the integral sign):

Name Differential form Integral form (using Kelvin–Stokes theorem plus relativistic invariance, \scriptstyle\int \frac{\partial}{\partial t} ...\to \frac{\mathrm d}{\mathrm dt} \int ...)
Maxwell–Faraday equation
Faraday's law of induction:
\ \nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}} {\partial t}\ \ \oint_C \mathbf{E} \cdot \mathrm{d}\mathbf{l} = \iint_S  \nabla \times \mathbf{E} \cdot \mathrm{d}\mathbf{A} = -\,\iint_S \frac{\partial \mathbf{B}}{\partial t} \cdot \mathrm{d}\mathbf{A}\ (with C and S not necessarily stationary)
Ampère's law
(with Maxwell's extension):
\ \ \nabla \times \mathbf{H} = \mathbf{J} + \frac{\partial \mathbf{D}} {\partial t}\ \ \,\ \begin{align}\oint_C \mathbf{H} \cdot \mathrm{d}\mathbf{l} &= \iint_S \nabla \times \mathbf{H} \cdot \mathrm{d} \mathbf{A}\\ &= \iint_S \mathbf{J} \cdot \mathrm{d} \mathbf{A} + \iint_S \frac{\partial \mathbf{D}}{\partial t} \cdot \mathrm{d} \mathbf{A}\end{align}\,\ (with C and S not necessarily stationary)

The above listed subset of Maxwell's equations are valid for electromagnetic fields expressed in SI units. In other systems of units, such as CGS or Gaussian units, the scaling factors for the terms differ. For example, in Gaussian units, Faraday's law of induction and Ampère's law take the forms[6][7]

\nabla \times \mathbf{E} = -\frac{1}{c} \frac{\partial \mathbf{B}} {\partial t},
\nabla \times \mathbf{H} = \frac{1}{c} \frac{\partial \mathbf{D}} {\partial t} + \frac{4\pi}{c} \mathbf{J},

respectively, where c is the speed of light in vacuum.

Divergence theorem

Likewise, the Ostrogradsky–Gauss theorem (also known as the divergence theorem or Gauss's theorem)

\int_{\mathrm{Vol}} \nabla \cdot \mathbf{F} \ \mathrm{d}_\mathrm{Vol} = \oint_{\partial \mathrm{Vol}} \mathbf{F} \cdot \mathrm{d} \mathbf{\Sigma}

is a special case if we identify a vector field with the n−1 form obtained by contracting the vector field with the Euclidean volume form. An application of this is the case \mathbf{F}=f\mathbf{c} where \mathbf c is an arbitrary constant vector. Working out the divergence of the product gives

\mathbf c \cdot \int_{\mathrm{Vol}} \nabla f \ \mathrm{d}_\mathrm{Vol} = \mathbf c \cdot \oint_{\partial \mathrm{Vol}} f \mathrm{d} \mathbf{\Sigma}

Since this holds for all \vec c, we find

\int_{\mathrm{Vol}} \nabla f \ \mathrm{d}_\mathrm{Vol} = \oint_{\partial \mathrm{Vol}} f \mathrm{d} \mathbf{\Sigma}

Notes

  1. See:
    • Victor J. Katz (May 1979), "The history of Stokes' theorem," Mathematics Magazine, 52 (3): 146–156.
    • The letter from Thomson to Stokes appears in: William Thomson and George Gabriel Stokes with David B. Wilson, ed., The Correspondence between Sir George Gabriel Stokes and Sir William Thomson, Baron Kelvin of Largs, Volume 1: 1846–1869 (Cambridge, England: Cambridge University Press, 1990), pages 96–97.
    • Neither Thomson nor Stokes published a proof of the theorem. The first published proof appeared in 1861 in: Hermann Hankel, Zur allgemeinen Theorie der Bewegung der Flüssigkeiten [On the general theory of the movement of fluids] (Göttingen, (Germany): Dieterische University Buchdruckerei, 1861); see pages 34–37. Hankel doesn't mention the author of the theorem.
    • In a footnote, Larmor mentions earlier researchers who had integrated, over a surface, the curl of a vector field. See: George G. Stokes with Sir Joseph Larmor and John Wm. Strutt (Baron Rayleigh), ed.s, Mathematical and Physical Papers by the late Sir George Gabriel Stokes, ... (Cambridge, England: University of Cambridge Press, 1905), vol. 5, pages 320–321.
  2. Olivier Darrigol,Electrodynamics from Ampere to Einstein, p. 146,ISBN 0198505930 Oxford (2000)
  3. 1 2 Spivak (1965), p. vii, Preface.
  4. See:
    • The 1854 Smith's Prize Examination is available on-line at: Clerk Maxwell Foundation. Maxwell took this examination and tied for first place with Edward John Routh in the Smith's Prize examination of 1854. See footnote 2 on page 237 of: James Clerk Maxwell with P. M. Harman, ed., The Scientific Letters and Papers of James Clerk Maxwell, Volume I: 1846–1862 (Cambridge, England: Cambridge University Press, 1990), page 237; see also Wikipedia's article "Smith's prize" or the Clerk Maxwell Foundation.
    • James Clerk Maxwell, A Treatise on Electricity and Magnetism (Oxford, England: Clarendon Press,1873), volume 1, pages 25–27. In a footnote on page 27, Maxwell mentions that Stokes used the theorem as question 8 in the Smith's Prize Examination of 1854. This footnote appears to have been the cause of the theorem's being known as "Stokes' theorem".
  5. For mathematicians this fact is known, therefore the circle is redundant and often omitted. However, one should keep in mind here that in thermodynamics, where frequently expressions as \oint_W \,\{d_{\,total\,}U\} appear (wherein the total derivative, see below, should not be confused with the exterior one), the integration path W is a one-dimensional closed line on a much higher-dimensional manifold. That is, in a thermodynamic application, where U is a function of the temperature \alpha_1:=T, the volume \alpha_2:=V, and the electrical polarization \alpha_3:=P of the sample, one has \{d_{\,total\,}U\}=\sum_{i=1}^3\,\frac{\partial U}{\partial\alpha_i}\mathrm d\alpha_i, and the circle is really necessary, e.g. if one considers the differential consequences of the integral postulate \oint_W\,\{d_{\,total\,}U\}\, \stackrel{!}{=}\,0.
  6. J.D. Jackson, Classical Electrodynamics, 2nd Ed (Wiley, New York, 1975).
  7. M. Born and E. Wolf, Principles of Optics, 6th Ed. (Cambridge University Press, Cambridge, 1980).

Further reading

External links

This article is issued from Wikipedia - version of the Thursday, January 07, 2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.