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. 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. 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.

• By the choice of F, $\tfrac{\mathrm{d}F}{\mathrm{d}x}=f(x)$. In the parlance of differential forms, this is saying that f(x) dx is the exterior derivative of the 0-form, i.e. function, F: in other words, that dF = f dx. The general Stokes theorem applies to higher differential forms ω instead of just 0-forms such as F.
• A closed interval [a, b] is a simple example of a one-dimensional manifold with boundary. Its boundary is the set consisting of the two points a and b. Integrating f over the interval may be generalized to integrating forms on a higher-dimensional manifold. Two technical conditions are needed: the manifold has to be orientable, and the form has to be compactly supported in order to give a well-defined integral.
• The two points a and b form the boundary of the open interval. More generally, Stokes' theorem applies to oriented manifolds M with boundary. The boundary ∂M of M is itself a manifold and inherits a natural orientation from that of the manifold. For example, the natural orientation of the interval gives an orientation of the two boundary points. Intuitively, a inherits the opposite orientation as b, as they are at opposite ends of the interval. So, "integrating" F over two boundary points a, b is taking the difference F(b) − F(a).

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 Rⁿ.

More generally, the integral of α over Ω is defined as follows: Let {ψ_i} be a partition of unity associated with a locally finite cover {U_i, φ_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 Rⁿ 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. 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 x_n 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 R^k to M. The free abelian group, S_k, 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 R^k. 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, S_k*, the linear functionals on S_k. 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 {c_i} are cycles generating the k-th homology group, then for any corresponding real numbers, {a_i}, 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

File:Stokes' Theorem.svg

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

$\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}$

Further reading

• Joos, Georg. Theoretische Physik. 13th ed. Akademische Verlagsgesellschaft Wiesbaden 1980. ISBN 3-400-00013-2
• Marsden, Jerrold E., Anthony Tromba. Vector Calculus. 5th edition W. H. Freeman: 2003.
• Lee, John. Introduction to Smooth Manifolds. Springer-Verlag 2003. ISBN 978-0-387-95448-6
• Stewart, James. Calculus: Concepts and Contexts. 2nd ed. Pacific Grove, CA: Brooks/Cole, 2001.
• Stewart, James. Calculus: Early Transcendental Functions. 5th ed. Brooks/Cole, 2003.

External links

• Proof of the Divergence Theorem and Stokes' Theorem
• Calculus 3 – Stokes Theorem from lamar.edu – an expository explanation