## Curl and Stokes’s theorem

Required math: vectors, calculus

Required physics: electrostatics

Reference: Griffiths, David J. (2007) Introduction to Electrodynamics, 3rd Edition; Prentice Hall – Sec. 1.3.5

One of the main vector derivatives is known as the curl. To see the motivation for it, we will consider a vector field ${\mathbf{A}}$ which, to make things concrete, can represent the flow of a fluid. We can imagine putting a little paddle wheel into the fluid at various points to measure the circulation of the fluid. If there is a vortex in the fluid, if we put a paddle wheel centred around the vortex, the wheel will turn. In general, any difference in tangential speed of the fluid from one side of the wheel to the other would cause the wheel to turn. For example, in a river, frequently the water flows more rapidly in the middle of the stream than it does at the banks, so if we put a paddle wheel (with a vertical axis, not a horizontal axis such as would be found in a water mill) anywhere in the river between the middle and the bank, the difference in speeds of the water across the wheel would cause it to turn.

As always, in calculus, we’re interested in how these ideas work on an infinitesimal scale. To keep things simple, we’ll consider only one component of the circulation: that in the ${xy}$ plane. Suppose we have a small rectangle in this plane with one corner at location ${(x_{0},y_{0})}$ and the diagonally opposite corner at ${(x_{0}+dx,y_{0}+dy)}$. To calculate the circulation around this rectangle, we would need to integrate the component of ${\mathbf{A}}$ that is parallel to each edge, with the integral extending around all four sides of the rectangle. That is, we would need to integrate ${\mathbf{A}\cdot d\mathbf{l}}$ around the rectangle, where ${d\mathbf{l}}$ is the incremental length vector that points along an edge. The direction in which we go around the rectangle will affect the sign of the answer, and by convention the integral is taken counterclockwise. That is, along the bottom edge ${d\mathbf{l}=dx\hat{\mathbf{i}}}$ and along the top edge, ${d\mathbf{l}=-dx\hat{\mathbf{i}}}$, and along the left edge ${d\mathbf{l}=-dy\hat{\mathbf{j}}}$, and right edge ${d\mathbf{l}=dy\hat{\mathbf{j}}}$.

Taking the dot product of ${\mathbf{A}}$ with ${d\mathbf{l}}$ selects out the ${x}$ component of ${\mathbf{A}}$ along the horizontal paths, and the ${y}$ component of ${\mathbf{A}}$ along the vertical paths. Since the components of ${\mathbf{A}}$ are functions of ${x}$ and ${y}$, we can use a first order expansion to get their values on all four sides of the rectangle.

At location ${(x_{0},y_{0})}$, we have ${A_{x}=A_{x}(x_{0},y_{0})}$, while at ${(x_{0},y_{0}+dy)}$ we get ${A_{x}=A_{x}(x_{0},y_{0}+dy)=A_{x}(x_{0},y_{0})+(\partial A_{x}/\partial y)(dy)}$ to first order. So to first order, the contribution to the circulation integral from the top and bottom edges is

$\displaystyle A_{x}(x_{0},y_{0})dx+\left[A_{x}(x_{0},y_{0})+\frac{\partial A_{x}}{\partial y}dy\right](-dx)=-\frac{\partial A_{x}}{\partial y}dxdy \ \ \ \ \ (1)$

By a similar argument we can work out the contributions from the left and right edges and get

$\displaystyle A_{y}(x_{0},y_{0})(-dy)+\left[A_{y}(x_{0},y_{0})+\frac{\partial A_{y}}{\partial x}dx\right]dy=\frac{\partial A_{y}}{\partial x}dxdy \ \ \ \ \ (2)$

Combining the two, the total circulation is

$\displaystyle \left(\frac{\partial A_{y}}{\partial x}-\frac{\partial A_{x}}{\partial y}\right)dxdy \ \ \ \ \ (3)$

Dividing out ${dxdy}$ we get the circulation per unit area as

$\displaystyle \frac{\partial A_{y}}{\partial x}-\frac{\partial A_{x}}{\partial y} \ \ \ \ \ (4)$

We can do a similar analysis in the other two coordinate planes, and we find that the three terms are the components of the curl of ${\mathbf{A}}$ which is defined as

$\displaystyle \nabla\times\mathbf{A}\equiv\left|\begin{array}{ccc} \hat{\mathbf{i}} & \hat{\mathbf{j}} & \hat{\mathbf{k}}\\ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z}\\ A_{x} & A_{y} & A_{z} \end{array}\right| \ \ \ \ \ (5)$

That is, the curl is defined analogously to a vector cross product, except the first ‘vector’ is the differential operator ${\nabla}$.

Stokes’s theorem can now be derived by considering some surface that is bounded by a closed path. Note that the surface itself is not closed (that is, it is not a surface like a sphere that has no edges), but apart from that, the exact form of the surface doesn’t matter. All that matters is the boundary curve.

We can divide the surface up into a number of little patches and perform the analysis above on each patch, then add up the contribution of each patch to the circulaton to get the total circulation around the bounding edge. The key point is that each edge of one patch that is shared with a neighbouring patch will be traversed in opposite directions as we do the integrals over the two patches, so the contributions of all these internal edges cancel out, leaving only the unmatched edges along the boundary. The net result is that the integral of the curl of the vector field over the surface is equal to the line integral of the vector field around the boundary. That is

$\displaystyle \int_{S}(\nabla\times\mathbf{A})\cdot d\mathbf{a}=\int_{L}\mathbf{A}\cdot d\mathbf{l} \ \ \ \ \ (6)$

where it is important to note that the integral on the left is taken over the surface, while that on the right is taken over the bounding curve. This is Stokes’s theorem.

There are several theorems involving the curl, but probably the most important one is that the curl of a gradient of a scalar field ${V}$ is always zero. That is

$\displaystyle \nabla\times(\nabla V)=0 \ \ \ \ \ (7)$

This can be seen most easily by substituting the gradient into the definition of the curl above:

$\displaystyle \nabla\times(\nabla V)=\left|\begin{array}{ccc} \hat{\mathbf{i}} & \hat{\mathbf{j}} & \hat{\mathbf{k}}\\ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z}\\ \frac{\partial V}{\partial x} & \frac{\partial V}{\partial y} & \frac{\partial V}{\partial z} \end{array}\right| \ \ \ \ \ (8)$

Since the order in which partial derivatives are calculated doesn’t matter (for any function of interest to physicists, anyway), this determinant is always zero.

A direct consequence of this in electrostatics is that, since the electric field can be expressed as the gradient of a potential function

$\displaystyle \mathbf{E}=-\nabla V \ \ \ \ \ (9)$

the curl of the electric field (in electrostatics) is always zero:

$\displaystyle \nabla\times\mathbf{E}=0 \ \ \ \ \ (10)$