四、矢量场的环量和旋度、斯托克斯定理

1、定义

矢量场 $A$ 沿闭合回路 $L$ 的线积分称为环量，用 ${\Gamma }_A$ 表示环量，则有
$${\Gamma }_A={\oint }_{L}A\cdot dl\text{\hspace{1em}(1)}$$
令 $\Delta S$ 为闭合曲线 $L$ 包围的面积，$e_n$ 为 $\Delta S$ 的右旋单位法线矢量。设想回路 $L$ 逐渐缩小，最后缩到空间某点 $P$。当 $\Delta S\to 0$ 时，${\Gamma }_A$ 也趋于0。若两者之比有一极限，则这极限值为矢量场 $A$ 的旋度（它是个矢量）在 $e_n$ 上的投影。$A$ 的旋度记作 $\mathrm{curl}$$A$ 或 $\mathrm{rot}$$A$，或 $▽\times A$。上述定义可写作
$$(▽\times A{)}_n=\underset{\Delta \to 0}{\lim }\frac{{\Gamma }_A}{\Delta S}=\underset{\Delta \to 0}{\lim }\frac{{\oint }_{L}A\cdot dl}{\Delta S}\text{\hspace{1em}(2)}$$
矢量场的旋度也是个矢量场。

2、旋度的坐标表达式

旋度直角坐标表示式的推导

下面我们来研究旋度的直角坐标表达式。先看旋度的 $x$ 分量。如上图（a），取一个与 $x$ 轴垂直的矩形回路 $L_x$，它的边分别与 $y,z$ 轴平行，边长为 $\Delta $ 和 $\Delta z$，取回路 $L_x$ 的环绕方向，使它的右旋法线 $e_n$ 指向 $+x$ 方向。设回路的中心 $P$ 点的坐标为 $(x,y,z)$，则在 $1,2,3,4$ 四边上的矢量场 $A$ 沿回路元的平行分量是 $A_z(x,y+\Delta y/2,z),-A_y(x,y,z+\Delta z/2),-A_z(x,y-\Delta y/2,z),A_y(x,y,z-\Delta z/2)$。
所以
$$\begin{array}{rl}{\oint }_{L_x}A\cdot dl& =A_z(x,y+\Delta y/2,z)\Delta z\\ & -A_z(x,y,-\Delta y/2)\Delta z\\ & -A_y(x,y,z+\Delta z/2)\Delta y\\ & +A_y(x,y,z-\Delta z/2)\Delta y\end{array}\text{\hspace{1em}(3)}$$
围绕 $P$ 点将 $A_y,A_z$ 按泰勒级数展开：
$$\begin{gathered} A_y(x,y,z\pm \Delta z/2)=A_y(x,y,z)\pm \frac{\partial A_y}{\partial z}\frac{\Delta z}{2}+高次项 \\ {A}_z(x,y\pm \Delta y/2,z)=A_z(x,y,z)\pm \frac{\partial A_z}{\partial y}\frac{\Delta y}{2}+高次项\text{\hspace{1em}(4)} \end{gathered}$$
带入前式，得
$${\oint }_{L_x}A\cdot dl=(\frac{\partial A_z}{\partial y}-\frac{\partial A_y}{\partial z})\Delta y\Delta z+高次项\text{\hspace{1em}(5)}$$
因为回路 $L_x$ 包围的矩形面积为 $\Delta S=\Delta y\Delta z$，按照旋度的定义式（2），得
$$\begin{array}{rl}(▽\times A{)}_x& \equiv \underset{\Delta S\to 0}{\lim }\frac{{\oint }_{L_x}A\cdot dl}{\Delta S}\\ & =\underset{\Delta y\to 0,{\Delta }_z\to 0}{\lim }\frac{(\frac{\partial A_z}{\partial y}-\frac{\partial A_y}{\partial z})\Delta y\Delta z+高次项}{\Delta y\Delta z}\\ & =\frac{\partial A_z}{\partial y}-\frac{\partial A_y}{\partial z}\end{array}\text{\hspace{1em}(6)}$$
同理可以得到旋度的 $y,z$ 两个分量。现将全部分量的直角坐标表示罗列如下：
$$\{\begin{array}{l}(▽\times A{)}_x=\frac{\partial A_z}{\partial y}-\frac{\partial A_y}{\partial z},\\ \\ (▽\times A{)}_y=\frac{\partial A_x}{\partial z}-\frac{\partial A_z}{\partial x},\\ \\ (▽\times A{)}_z=\frac{\partial A_y}{\partial x}-\frac{\partial A_x}{\partial y},\end{array}\text{\hspace{1em}(7)}$$
旋度矢量的直角坐标表达式为
$$\begin{array}{rl}▽\times A& =(\frac{\partial A_z}{\partial y}-\frac{\partial A_y}{\partial z})i+(\frac{\partial A_x}{\partial z}-\frac{\partial A_z}{\partial x})j+(\frac{\partial A_y}{\partial x}-\frac{\partial A_x}{\partial y})k\\ & =\left|\begin{array}{ccc}i& j& k\\ & & \\ \frac{\partial }{\partial x}& \frac{\partial }{\partial y}& \frac{\partial }{\partial z}\\ & & \\ A_x& A_y& A_z\end{array}\right|\end{array}\text{\hspace{1em}(8)}$$
下面不加推导地给出旋度在其他常用坐标中的表达式，必备参考。
柱坐标
$$▽\times A=(\frac{1}{\rho }\frac{\partial A_z}{\partial \phi }-\frac{\partial A_{\phi }}{\partial z})e_{\rho }+(\frac{\partial A_{\rho }}{\partial z}-\frac{\partial A_z}{\partial \rho })e_{\phi }+\frac{1}{\rho }(\frac{\partial (\rho A_{\phi })}{\partial \rho }-\frac{\partial A_{\rho }}{\partial \phi })e_z\text{\hspace{1em}(9)}$$
球坐标
$$▽\times A=\left[\frac{1}{r\sin \theta }\right(\frac{\partial }{\partial \theta }(\sin \theta A_{\phi })-\frac{\partial A_{\theta }}{\partial \phi }\left)\right]e_r+[\frac{1}{r\sin \theta }\frac{\partial A_r}{\partial \phi }-\frac{1}{r}(\frac{\partial }{\partial r}(r{A}_{\phi })\left)\right]e_{\theta }+[\frac{1}{r}\frac{\partial (r{A}_{\theta })}{\partial r}-\frac{1}{r}\frac{\partial A_r}{\partial \theta }]e_{\phi }\text{\hspace{1em}(10)}$$

3、斯托克斯定理

斯托克斯定理的证明

在矢量场 $A(x,y,z)$ 中取任意闭合回路 $L$，现用一条曲线搭在回路 $L$ 上的 $M,N$ 两点之间。$M$ 和 $N$ 把 $L$ 分为 $L_1^{\prime }$ 和 $L_2^{\prime }$ 两部分，如上图（a）。$L_1^{\prime }$ 和 $MN$ 组成新的闭合回路 $L_1$，$L_2^{\prime }$ 和 $MN$ 组成新的闭合回路 $L_2$，$L_1$ 和 $L_2$ 的环绕方向一致，沿 $L_1$ 和 $L_2$ 的环量分别是
$$\begin{gathered} {\Gamma }_1\equiv {\oint }_{L_1}A\cdot dl={\int }_{L_1^{\prime }}A\cdot dl+{\int }_M^{N}A\cdot dl, \\ {\Gamma }_2\equiv {\oint }_{L_2}A\cdot dl={\int }_{L_2^{\prime }}A\cdot dl+{\int }_N^{M}A\cdot dl.\text{\hspace{1em}(11)} \end{gathered}$$
故
$${\Gamma }_1+{\Gamma }_2={\int }_{L_1^{\prime }}A\cdot dl+{\int }_{L_2^{\prime }}A\cdot dl={\oint }_{L}A\cdot dl\equiv \Gamma \text{\hspace{1em}(12)}$$
即矢量场在闭合回路 $L$ 上的环量等于分割出来的两个闭合回路 $L_1$ 和 $L_2$ 上的环量之和。如图(b)，用许多曲线，像织成的网子一样绷在回路 $L$ 的“框架”上，则每个网眼是一个闭合回路 $L_i$，令它们的环绕方向都一致，用 ${\Gamma }_i$ 代表 $L_i$ 上的环量，则有
$$\Gamma =\sum _{i=1}^n{\Gamma }_i\text{\hspace{1em}(13)}$$
这就是说，$L$ 上的环量是由各局部的环量累积起来的。
如果把上述分割过程无限继续下去，使每个小回路的面积 $\Delta S_i$ 都区域0，则按照旋度的定义，
$${\Gamma }_i\equiv {\oint }_{L}A\cdot dl=(▽\times A{)}_{ni}\Delta S_i=(▽\times A)\cdot \Delta S_i\text{\hspace{1em}(14)}$$
这里$(▽\times A{)}_{ni}$代表旋度$(▽\times A)$在 $\Delta S_i$ 的右旋单位法线矢量 $e_{ni}$ 上的投影，$\Delta S_i=e_{ni}\Delta S_i$ 是矢量面元，带入式(13)，得
$$\Gamma \equiv {\oint }_{L}A\cdot dl=\sum _{i=1}^n(▽\times A)\cdot \Delta S_i\text{\hspace{1em}(15)}$$
取极限后，右端变为面积分：
$${\oint }_{L}A\cdot dl={\iint }_S(▽\times A)\cdot dS\text{\hspace{1em}(16)}$$
式(16)表明：矢量场在任何闭合回路 $L$ 上的环量等于以它为边界的曲面 $S$ 上的旋度的积分。这就是斯托克斯定理。

4、备注

我们已经多次使用了符号 “$▽$”，单一直是将它和一个场函数 $\Phi$ 或 $A$ 连起来写，而未说明它单独代表什么。其实 $▽$ 是一个矢量性质的算符，叫做劈型算符或那勃勒算符，也称矢量微分算子，它的直角坐标表达式为
$$▽=i\frac{\partial }{\partial x}+j\frac{\partial }{\partial y}+k\frac{\partial }{\partial z}$$
我们可以把 $▽$ 形式地“乘”在一个标量场 $\Phi$ 上，称为它的梯度 $▽\Phi$，也可以把它形式地“点乘”或“叉乘”在一个矢量场 $A$ 上，成为它的散度 $▽\cdot A$ 或旋度 $▽\times A$，这样做的结果，我们得到的正是前面的式子。