Curl & Stokes' Theorem

⚡ Key Takeaways

Stokes' Theorem equates the circulation of a vector field around a closed curve to the flux of its curl through any surface bounded by that curve: ∮∂SF⋅dr=∬S(∇×F)⋅dS\oint_{\partial S} \mathbf{F}\cdot d\mathbf{r} = \iint_S (\nabla\times\mathbf{F})\cdot d\mathbf{S}. In words: the net "swirl" of a field over a surface equals the total circulation around its rim. Applies to: electromagnetism, fluid dynamics, differential geometry, any oriented surface with a piecewise-smooth boundary. Who/when: stated by William Thomson (Lord Kelvin) in 1850, popularized as an exam question by George Stokes at Cambridge in 1854.

💡 The Core Insight

The curl ∇×F\nabla\times\mathbf{F} measures the infinitesimal circulation per unit area at a point — imagine a tiny paddlewheel placed in the field. Stokes' Theorem says these microscopic circulations, when summed (integrated) over a surface, do not cancel only along the outer boundary: every interior edge is shared by two adjacent cells with opposite orientations, so their contributions annihilate. What survives is the circulation around the rim. This is the same telescoping cancellation that powers all generalized Stokes-type theorems.

Formal Statement

DefinitionCurl of a Vector Field

For a vector field F=(P,Q,R)\mathbf{F} = (P, Q, R) with continuously differentiable components, the curl is the vector field

1
∇×F=(∂R∂y−∂Q∂z)i+(∂P∂z−∂R∂x)j+(∂Q∂x−∂P∂y)k\quad \nabla\times\mathbf{F} = \left(\frac{\partial R}{\partial y}-\frac{\partial Q}{\partial z}\right)\mathbf{i} + \left(\frac{\partial P}{\partial z}-\frac{\partial R}{\partial x}\right)\mathbf{j} + \left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)\mathbf{k}

Derivation hint: expand the formal determinant ∇×F=det⁡∣ijk∂x∂y∂zPQR∣\nabla\times\mathbf{F} = \det\begin{vmatrix}\mathbf{i}&\mathbf{j}&\mathbf{k}\\ \partial_x&\partial_y&\partial_z\\ P&Q&R\end{vmatrix} along the top row. Here ∇=(∂x,∂y,∂z)\nabla = (\partial_x, \partial_y, \partial_z) is the del operator. The curl is itself a vector; its direction is the axis of maximal rotation (by the right-hand rule) and its magnitude is twice the local angular velocity of an infinitesimal fluid element.

DefinitionCirculation and Flux

The circulation of F\mathbf{F} around a closed curve CC is the line integral

2
Γ=∮CF⋅dr=∮C(P dx+Q dy+R dz)\quad \Gamma = \oint_C \mathbf{F}\cdot d\mathbf{r} = \oint_C (P\,dx + Q\,dy + R\,dz)

Derivation hint: parametrize CC by r(t)\mathbf{r}(t), a≤t≤ba\le t\le b, so dr=r′(t) dtd\mathbf{r} = \mathbf{r}'(t)\,dt and the integrand becomes F(r(t))⋅r′(t) dt\mathbf{F}(\mathbf{r}(t))\cdot\mathbf{r}'(t)\,dt. The flux of the curl through an oriented surface SS with unit normal n^\hat{\mathbf{n}} is

3
Φ=∬S(∇×F)⋅dS=∬S(∇×F)⋅n^ dS\quad \Phi = \iint_S (\nabla\times\mathbf{F})\cdot d\mathbf{S} = \iint_S (\nabla\times\mathbf{F})\cdot\hat{\mathbf{n}}\,dS

Derivation hint: dS=n^ dSd\mathbf{S} = \hat{\mathbf{n}}\,dS; for a parametrized surface r(u,v)\mathbf{r}(u,v), dS=(ru×rv) du dvd\mathbf{S} = (\mathbf{r}_u\times\mathbf{r}_v)\,du\,dv.

TheoremStokes' Theorem

Let SS be a piecewise-smooth oriented surface in R3\mathbb{R}^3 bounded by a piecewise-smooth, simple, closed curve ∂S\partial S with the induced (right-hand-rule) orientation. Let F\mathbf{F} be a vector field whose components have continuous partial derivatives on an open region containing SS. Then

4
∮∂SF⋅dr=∬S(∇×F)⋅dS\quad \oint_{\partial S} \mathbf{F}\cdot d\mathbf{r} = \iint_S (\nabla\times\mathbf{F})\cdot d\mathbf{S}

Derivation hint: reduces to <a class="wiki-link" href="/wiki/ndli-iit-kgp-ma10002-mathematics-ii-greens-theorem">Green&#39;s Theorem</a> when SS is flat in the xyxy-plane; the general case follows by pulling SS back to a parameter domain.

Orientation convention: if you walk along ∂S\partial S with your head pointing in the direction of n^\hat{\mathbf{n}}, the surface SS lies on your left. Equivalently, curl the fingers of your right hand along ∂S\partial S; the thumb gives n^\hat{\mathbf{n}}.

Figure
Stokes' Theorem: circulation around the rim equals flux of curl through the cap
n̂ ∂S S ∇×F ∮∂S F·dr = ∬S (∇×F)·dS

Why It Holds / Derivation

ProofFrom Green's Theorem to the General Surface

We prove the theorem by reducing to Green's Theorem and then assembling the general case.

Step 1 — Flat case. Suppose SS lies in the xyxy-plane with upward normal n^=k\hat{\mathbf{n}}=\mathbf{k}, and F=(P,Q,0)\mathbf{F}=(P,Q,0). Then (∇×F)⋅k=Qx−Py(\nabla\times\mathbf{F})\cdot\mathbf{k} = Q_x - P_y, and

5
∬S(Qx−Py) dA=∮∂S(P dx+Q dy)\quad \iint_S (Q_x - P_y)\,dA = \oint_{\partial S}(P\,dx + Q\,dy)

Derivation hint: this is exactly the circulation form of Green's Theorem.

Step 2 — Parametrized surface. Let r(u,v)\mathbf{r}(u,v), (u,v)∈D(u,v)\in D, parametrize SS with ∂D\partial D mapping to ∂S\partial S. The line integral pulls back:

6
∮∂SF⋅dr=∮∂D(F⋅ru) du+(F⋅rv) dv\quad \oint_{\partial S}\mathbf{F}\cdot d\mathbf{r} = \oint_{\partial D}\Big(\mathbf{F}\cdot\mathbf{r}_u\Big)\,du + \Big(\mathbf{F}\cdot\mathbf{r}_v\Big)\,dv

Derivation hint: along ∂D\partial D, dr=ru du+rv dvd\mathbf{r} = \mathbf{r}_u\,du + \mathbf{r}_v\,dv; substitute into F⋅dr\mathbf{F}\cdot d\mathbf{r}.

Step 3 — Apply Green in the parameter plane. Set P~=F⋅ru\tilde P = \mathbf{F}\cdot\mathbf{r}_u, Q~=F⋅rv\tilde Q = \mathbf{F}\cdot\mathbf{r}_v. Green's Theorem in DD gives ∮∂DP~ du+Q~ dv=∬D(Q~u−P~v) du dv\oint_{\partial D}\tilde P\,du + \tilde Q\,dv = \iint_D (\tilde Q_u - \tilde P_v)\,du\,dv. Computing Q~u−P~v\tilde Q_u - \tilde P_v by the product and chain rules, the mixed second-derivative terms involving ruv\mathbf{r}_{uv} cancel by equality of mixed partials, leaving

7
Q~u−P~v=(∇×F)⋅(ru×rv)\quad \tilde Q_u - \tilde P_v = (\nabla\times\mathbf{F})\cdot(\mathbf{r}_u\times\mathbf{r}_v)

Derivation hint: expand ∂u(F⋅rv)−∂v(F⋅ru)\partial_u(\mathbf{F}\cdot\mathbf{r}_v) - \partial_v(\mathbf{F}\cdot\mathbf{r}_u); the surviving terms are first derivatives of F\mathbf{F} contracted against ru×rv\mathbf{r}_u\times\mathbf{r}_v.

Step 4 — Recognize the flux. Since dS=(ru×rv) du dvd\mathbf{S} = (\mathbf{r}_u\times\mathbf{r}_v)\,du\,dv,

8
∬D(∇×F)⋅(ru×rv) du dv=∬S(∇×F)⋅dS\quad \iint_D (\nabla\times\mathbf{F})\cdot(\mathbf{r}_u\times\mathbf{r}_v)\,du\,dv = \iint_S(\nabla\times\mathbf{F})\cdot d\mathbf{S}

Derivation hint: substitute the surface-element formula directly. Combining Steps 2–4 yields ∮∂SF⋅dr=∬S(∇×F)⋅dS\oint_{\partial S}\mathbf{F}\cdot d\mathbf{r} = \iint_S(\nabla\times\mathbf{F})\cdot d\mathbf{S}. General surfaces are decomposed into such patches; internal edges cancel pairwise. ■\blacksquare

∎
GG
George Gabriel Stokes
1819–1903

Stokes was Lucasian Professor of Mathematics at Cambridge for over fifty years. He did not discover "his" theorem — Lord Kelvin communicated it to him in a letter of July 1850. Stokes' lasting contribution was pedagogical: he set the identity as a problem on the Smith's Prize Examination in February 1854, ensuring generations of physicists (including James Clerk Maxwell, who sat the exam) encountered it. Stokes himself made monumental contributions to fluid dynamics (the Navier–Stokes equations), optics (fluorescence, Stokes shift), and the theory of light. The naming is one of mathematics' famous historical accidents.

Building Intuition

🧠 Intuition

Picture a surface SS tiled into tiny squares. In each square, the curl gives the local circulation per unit area. Now sum the circulations around the boundaries of all these little squares. Every interior edge is shared by exactly two squares, traversed in opposite directions — so those contributions cancel exactly. The only edges not cancelled are those on the outer rim ∂S\partial S. Therefore the total interior swirl (the flux of curl) must equal the leftover boundary circulation. This is the discrete heart of the theorem; refining the tiles to zero size makes it exact. A second intuition: the curl points along the axis a paddlewheel would spin about, and its length is the spin rate. Flux of curl counts how much "spin-axis" pokes through SS. The theorem says that total is fixed by what the flow does along the rim — the interior detail is irrelevant.

Worked Examples

ExampleHemisphere with a Rotational Field

Given: F=(−y,x,z)\mathbf{F} = (-y, x, z), and SS is the upper unit hemisphere x2+y2+z2=1x^2+y^2+z^2=1, z≥0z\ge 0, oriented with outward (upward) normal. Compute ∬S(∇×F)⋅dS\iint_S(\nabla\times\mathbf{F})\cdot d\mathbf{S}.

Step 1 — Use Stokes to convert to a line integral. The boundary ∂S\partial S is the unit circle x2+y2=1x^2+y^2=1 in the plane z=0z=0, oriented counterclockwise (viewed from above).

Step 2 — Parametrize the boundary. r(t)=(cos⁡t,sin⁡t,0)\mathbf{r}(t)=(\cos t,\sin t, 0), 0≤t≤2π0\le t\le 2\pi, so dr=(−sin⁡t,cos⁡t,0) dtd\mathbf{r}=(-\sin t,\cos t,0)\,dt.

Step 3 — Evaluate F\mathbf{F} on ∂S\partial S. F=(−sin⁡t,cos⁡t,0)\mathbf{F}=(-\sin t,\cos t,0).

Step 4 — Dot and integrate.

F⋅dr=(sin⁡2t+cos⁡2t) dt=dt.\mathbf{F}\cdot d\mathbf{r} = (\sin^2 t + \cos^2 t)\,dt = dt.
∮∂SF⋅dr=∫02π1 dt=2π.\oint_{\partial S}\mathbf{F}\cdot d\mathbf{r} = \int_0^{2\pi} 1\,dt = 2\pi.

Check: ∇×F=(0,0,2)\nabla\times\mathbf{F}=(0,0,2); the flat disk of radius 1 has flux ∬2 dA=2⋅π=2π\iint 2\,dA = 2\cdot\pi = 2\pi. The hemisphere gives the same value — surface independence confirmed.

∬S(∇×F)⋅dS=2π\boxed{\iint_S(\nabla\times\mathbf{F})\cdot d\mathbf{S} = 2\pi}
ExampleTriangle Boundary, Direct Line Integral Avoided

Given: F=(y2,x,z2)\mathbf{F}=(y^2, x, z^2). CC is the triangle with vertices (0,0,0)(0,0,0), (1,0,0)(1,0,0), (0,1,0)(0,1,0) oriented counterclockwise in the xyxy-plane. Compute ∮CF⋅dr\oint_C\mathbf{F}\cdot d\mathbf{r}.

Step 1 — Choose surface SS. The flat triangular region TT in z=0z=0, normal n^=k\hat{\mathbf{n}}=\mathbf{k}.

Step 2 — Curl.

∇×F=(∂yz2−∂zx, ∂zy2−∂xz2, ∂xx−∂yy2)=(0, 0, 1−2y).\nabla\times\mathbf{F}=\big(\partial_y z^2 - \partial_z x,\ \partial_z y^2 - \partial_x z^2,\ \partial_x x - \partial_y y^2\big) = (0,\ 0,\ 1 - 2y).

Step 3 — Flux through TT. (∇×F)⋅k=1−2y(\nabla\times\mathbf{F})\cdot\mathbf{k} = 1-2y. The region: 0≤x≤1−y0\le x \le 1-y, 0≤y≤10\le y\le 1.

∬T(1−2y) dA=∫01 ⁣ ⁣∫01−y(1−2y) dx dy=∫01(1−2y)(1−y) dy.\iint_T (1-2y)\,dA = \int_0^1\!\!\int_0^{1-y}(1-2y)\,dx\,dy = \int_0^1 (1-2y)(1-y)\,dy.

Step 4 — Integrate.

∫01(1−3y+2y2) dy=[y−32y2+23y3]01=1−32+23=16.\int_0^1 (1 - 3y + 2y^2)\,dy = \Big[y - \tfrac{3}{2}y^2 + \tfrac{2}{3}y^3\Big]_0^1 = 1 - \tfrac{3}{2} + \tfrac{2}{3} = \tfrac{1}{6}.
∮CF⋅dr=16\boxed{\oint_C\mathbf{F}\cdot d\mathbf{r} = \tfrac{1}{6}}

Common Misconceptions

⚠️ Orientation Must Match

The most common error is mismatched orientation between ∂S\partial S and n^\hat{\mathbf{n}}. If you reverse the surface normal but keep the same boundary orientation, the answer flips sign. Always verify with the right-hand rule: fingers curl along ∂S\partial S, thumb gives n^\hat{\mathbf{n}}. A second pitfall: Stokes' Theorem requires F\mathbf{F} to be C1C^1 on all of SS — a singularity anywhere on the chosen surface invalidates the equality, even if the boundary is fine.

CounterexampleCurl-Free but Nonzero Circulation

Consider the "vortex" field F=(−yx2+y2,xx2+y2,0)\mathbf{F} = \left(\dfrac{-y}{x^2+y^2}, \dfrac{x}{x^2+y^2}, 0\right). A direct computation gives ∇×F=0\nabla\times\mathbf{F}=\mathbf{0} everywhere it is defined. Yet around the unit circle,

∮CF⋅dr=2π≠0.\oint_C\mathbf{F}\cdot d\mathbf{r} = 2\pi \ne 0.

Stokes' Theorem is not violated: every surface spanning the unit circle must pass through the zz-axis, where F\mathbf{F} is undefined (the curl has a delta-function singularity there). The hypothesis "F\mathbf{F} is C1C^1 on SS" fails. This example shows circulation can be "trapped" by a hole — the seed of De Rham Cohomology.

Connections

connection

Stokes' Theorem sits at the center of a web of results. Green's Theorem is its flat, two-dimensional special case (and supplies the proof). The Divergence Theorem is its sibling, relating flux through a closed surface to a volume integral of divergence. All three are unified by the Generalized Stokes Theorem for differential forms: ∫Mdω=∫∂Mω\int_M d\omega = \int_{\partial M}\omega. The vanishing identity ∇⋅(∇×F)=0\nabla\cdot(\nabla\times\mathbf{F})=0 explains surface independence, and connects to De Rham Cohomology — exact circulation depends only on the boundary's homology class. In physics, Stokes' Theorem converts Maxwell's integral and differential forms into one another via Maxwell's Equations.

Real-World Applications

🔬 Where this appears

Historical Development

⏳ History
- 1760s–1780s — Lagrange and others develop line and surface integrals in the study of mechanics and potential theory.
- 1828 — George Green publishes his Essay on the Application of Mathematical Analysis to the Theories of Electricity and Magnetism, introducing Green's Theorem.
- 1831 — Faraday discovers electromagnetic induction, motivating circulation–flux relations physically.
- July 1850 — William Thomson (Lord Kelvin) states the curl theorem in a letter to Stokes.
- February 1854 — Stokes sets the theorem as Question 8 on the Smith's Prize Examination at Cambridge; James Clerk Maxwell is among the candidates.
- 1861–1873 — Maxwell employs the theorem extensively in his treatment of electromagnetism, embedding it in A Treatise on Electricity and Magnetism.
- 1880s–1890s — Heaviside and Gibbs cast the result in modern vector-calculus notation with ∇×\nabla\times.
- 1899 — Élie Cartan develops the exterior calculus of differential forms.
- 1930s–1940s — Cartan and Georges de Rham complete the Generalized Stokes Theorem ∫Mdω=∫∂Mω\int_M d\omega = \int_{\partial M}\omega, unifying all classical integral theorems.

Subtleties & Limits

open-problem

Stokes' Theorem is fully proven for smooth manifolds with boundary, but its boundaries (literal and figurative) remain rich. Regularity: the classical statement requires C1C^1 fields and piecewise-smooth surfaces; extending to rough/fractal boundaries and to currents (in the sense of geometric measure theory) is the domain of Federer–Fleming theory and is technically delicate. Topology: when surfaces have holes or are non-orientable (e.g., the Möbius band), the naïve theorem fails — orientability is essential, and the failure is measured by De Rham Cohomology and homology. Singular fields: as the vortex counterexample shows, isolated singularities act like distributional sources; making this rigorous leads to currents and the theory of distributions. Whether the cleanest formulation lives in the smooth category, the PL category, or geometric measure theory depends on the application, and the interplay remains an active interface between analysis and topology.

See Also

Green's Theorem · Divergence Theorem · Generalized Stokes Theorem · Differential Forms · De Rham Cohomology · Maxwell's Equations · Vector Fields and the Del Operator · Vorticity and Circulation