Stokes公式

Date：2021/11/24

Author：自倚修行

一、基础知识

1.1 Hamilton算子

定义：

\[\nabla = (\frac{\part}{\part x}, \frac{\part}{\part y}, \frac{\part}{\part z}) \]

形式上，可以将Hamilton算子看作向量。同时，Hamilton算子还有二阶形式 \(\nabla = (\frac{\part}{\part x},\frac{\part}{\part y})\)。

利用Hamilton算子，可以定义如下概念：

梯度：

\[\mathbf{grad}\ F = \nabla F = (\frac{\part F}{\part x}, \frac{\part F}{\part y}, \frac{\part F}{\part z}) \]

散度：

\[\mathbf{div}\ F = \nabla\cdot F = \frac{\part F}{\part x} + \frac{\part F}{\part y} + \frac{\part F}{\part z} \]

旋度：

\[\mathbf{rot}\ F = \nabla\times F = \left|\begin{matrix} i & j & k \\ \frac{\part}{\part x}& \frac{\part}{\part y}& \frac{\part}{\part z} \\ P & Q & R \end{matrix}\right|,\ F = P\ \mathbf{i} + Q\ \mathbf{j} + R\ \mathbf{k} \]

1.2 外微分

外微分是对 \(k\) - 形式微分的特殊称呼。

\(k\) - 形式通常如下定义：

\(k\) - 形式：

\[\omega = \sum\limits_{1\le i_1

举例如下：

\(0\) - 形式：

\[\omega = f(x,y,z) \]

\(1\) - 形式：

\[\omega = f(x,y,z)\mathbf{d}x + g(x,y,z)\mathbf{d}y + r(x,y,z)\mathbf{d}z \]

\(2\) - 形式：

\[\omega = f(x,y,z)\mathbf{d}y\wedge \mathbf{d}z + g(x,y,z)\mathbf{d}z\wedge \mathbf{d}x + r(x,y,z)\mathbf{d}x\wedge \mathbf{d}y \]

一般如下定义外微分：

外微分定义：

\[\begin{aligned} \mathbf{d} \omega &= \sum\limits_{1\le i_1<\cdots

不加证明地给出如下性质（\(\omega\) 是 \(k\) - 形式，\(\eta\) 是 \(l\) - 形式）：

\[\mathbf{d} x_i \wedge \mathbf{d}x_i = 0 \\ \mathbf{d}^2 \omega = \mathbf{d}(\mathbf{d} \omega) = 0 \\ \mathbf{d}(\omega \wedge \eta) = \mathbf{d}\omega \wedge \eta + (-1)^k \omega \wedge \mathbf{d}\eta \\ \]

二、曲线积分

2.1 第一类曲线积分

第一类曲线积分是对曲线质量的积分，可以用于计算曲线长度。

首先考虑第一类曲线积分中的弧长微元，由勾股定理：

\[\begin{aligned} \Delta s &= \sqrt{[\Delta x]^2+[\Delta y]^2+[\Delta y]^2} \\ \mathrm{d}s &= \sqrt{[x'(t)]^2+[y'(t)]^2+[z'(t)]^2} \mathrm{d}t \\ \end{aligned} \]

自然的，计算第一类曲线积分：

\[\int_L f(x,y,z) \mathrm{d}s = \int_\alpha^\beta f(x(t),y(t),z(t))\sqrt{[x'(t)]^2+[y'(t)]^2+[z'(t)]^2} \mathrm{d}t \]

2.2 第二类曲线积分

第二类曲线积分是求力沿路径所做的功，写法灵活。

第二类曲线积分中的路径微元是如下向量：

\[\mathrm{d}\mathbf{s} = (\mathrm{d}x, \mathrm{d}y, \mathrm{d}z) \]

第二类曲线积分有如下写法：

\[\begin{aligned} \int_L F\cdot \mathrm{d}\mathbf{s} &= \int_L P\mathrm{d}x+Q\mathrm{d}y+R\mathrm{d}z \\ &= \int_L [P\cos\alpha+Q\cos\beta+R\cos\gamma]\mathrm{d}s \\ &= \int_L F \cdot \tau\ \mathrm{d}s \end{aligned} \]

其中，\(\tau\) 是曲线在各点处的单位切向量，\(\alpha\ \beta \ \gamma\) 为 \(\tau=(\cos\alpha,\cos\beta,\cos\gamma)\) 中定义。

计算第二类曲线积分主要是利用 \(1\) - 形式：

\[\begin{aligned} \int_L F\cdot \mathrm{d}\mathbf{s} &= \int_L P\mathrm{d}x+Q\mathrm{d}y+R\mathrm{d}z \\ &= \int_\alpha^\beta [Px'(t)+Qy'(t)+Rz'(t)]\mathrm{d}t \end{aligned} \]

三、曲面积分

3.1 第一类曲面积分

第一类曲面积分是对曲面质量的积分，可以用于计算曲面面积。

首先考虑第一类曲面积分中的面积微元，利用面积在各坐标平面的投影：

\[\mathrm{d}S = \cos\alpha \mathrm{d}y\mathrm{d}z + \cos\beta \mathrm{d}z\mathrm{d}x+\cos\gamma \mathrm{d}x\mathrm{d}y \]

计算第一类曲面积分需要利用向量叉积将 \(x,y,z\) 空间的面积微元投射到 \(u,v\) 空间。

设 \(x = x(u,v),y = y(u,v), z = z(u,v), r=(x,y,z)\)，则：

\[\mathrm{d}S = \| r_u \times r_v \| \mathrm{d}u\mathrm{d}v \]

其中的变换参数：

\[\begin{aligned} \| r_u \times r_v \| &= \left|\begin{matrix} i & j & k \\ x_u & y_u & z_u \\ x_v & y_v & z_v \end{matrix}\right| \\ &= \sqrt{\left[\frac{\partial(y, z)}{\partial(u, v)}\right]^{2}+\left[\frac{\partial(z, x)}{\partial(u, v)}\right]^{2}+\left[\frac{\partial(x, y)}{\partial(u, v)}\right]^{2}} \\ &= \sqrt{EG-F^2} \end{aligned} \]

此时有 \(E = r_u\cdot r_u,\ F = r_u \cdot r_v,\ G = r_v\cdot r_v\) 称为 \(\text{Gauss}\) 系数。

自然的，计算第一类曲面积分：

\[\int_\Sigma f(x,y,z) \mathrm{d}S = \int_{\Sigma'} f(x(u,v), y(u,v), z(u,v))\sqrt{EG-F^2}\ \mathrm{d}u\mathrm{d}v \]

其中，\(\Sigma'\) 是曲面 \(\Sigma\) 在 \(u,v\) 空间上的投影。

3.2 第二类曲面积分

第二类曲面积分是对过曲面的流量进行积分，写法灵活。

第二类曲面积分中的面积微元是如下向量：

\[\mathrm{d}\mathbf{S} = (\mathrm{d}y\mathrm{d}z, \mathrm{d}z\mathrm{d}x, \mathrm{d}x\mathrm{d}y) \]

第二类曲面积分有如下写法：

\[\begin{aligned} \int_\Sigma F\cdot \mathrm{d}\mathbf{S} &= \int_\Sigma P \mathrm{d}y\mathrm{d}z+Q\mathrm{d}z\mathrm{d}x+R\mathrm{d}x\mathrm{d}y \\ &= \int_\Sigma [P\cos\alpha+Q\cos\beta+R\cos\gamma] \mathrm{d}S \\ &= \int_\Sigma F\cdot \mathbf{n}\ \mathrm{d}S \end{aligned} \]

其中，\(\mathbf{n}\) 是曲面微元的单位外法向量，\(\alpha\ \beta \ \gamma\) 为 \(\mathbf{n}=(\cos\alpha,\cos\beta,\cos\gamma)\) 中定义。

计算第二类曲面积分，通常会使用被积曲面的各类对称性。一般地，求解时还是利用 \(x,y,z\) 空间的曲面微元投射到 \(u,v\) 空间，然后对 \(2\) - 形式进行积分。

\[\begin{aligned} \mathrm{d}\mathbf{S} &= (\mathrm{d}y\mathrm{d}z, \mathrm{d}z\mathrm{d}x, \mathrm{d}x\mathrm{d}y) \\ &= \left( \frac{\part(y,z)}{\part(u,v)},\frac{\part(z,x)}{\part(u,v)},\frac{\part(x,y)}{\part(u,v)} \right)\mathrm{d}u\mathrm{d}v \\ &= \pm\sqrt{EG-F^2}\ \mathbf{n}\ \mathrm{d}u\mathrm{d}v \\ &= \pm \mathbf{n}\ \mathrm{d}S \\ \end{aligned} \]

所以有：

\[\begin{aligned} \int_\Sigma F\cdot \mathrm{d}\mathbf{S} &= \int_\Sigma P \mathrm{d}y\mathrm{d}z+Q\mathrm{d}z\mathrm{d}x+R\mathrm{d}x\mathrm{d}y \\ &= \pm\int_{\Sigma'} \left[P\frac{\part(y,z)}{\part(u,v)}+Q\frac{\part(z,x)}{\part(u,v)}+R\frac{\part(x,y)}{\part(u,v)}\right] \mathrm{d}u\mathrm{d}v \\ \end{aligned} \]

式子中的 \(\pm\) 可以利用外法方向向量 \(\mathbf{n}\) 和有向被积曲面 \(\Sigma\) 方向是否一致来判断（一致则取正）。上式中 \(\Sigma'\) 是曲面 \(\Sigma\) 在 \(u,v\) 空间上的投影。

四、常见公式

一维 \(1\) - 形式 \(\Leftrightarrow\) \(0\) - 形式：

\(\text{Newton-Leibniz}\) 公式

\[\int_a^b f'(x) \mathrm{d}x = f(x)|_a^b \]

二维 \(2\) - 形式 \(\Leftrightarrow\) \(1\) - 形式：

\(\text{Green}\) 公式

\[\iint_{D}\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right) \mathrm{d} x \mathrm{d} y=\oint_{\part D} P \mathrm{d} x+Q \mathrm{d} y \]

其中左端为标准 \(2\) - 形式，可利用外微分对右端 \(1\) - 形式求导获得。

三维 \(3\) - 形式 \(\Leftrightarrow\) \(2\) - 形式：

\(\text{Gauss}\) 公式

\[\iiint_{V} \nabla \cdot \boldsymbol{F}\ \mathrm{d} V=\iint_{\partial V} \boldsymbol{F} \cdot \mathrm{d} \boldsymbol{S} \]

三维 \(2\) - 形式 \(\Leftrightarrow\) \(1\) - 形式：

\(\text{Stokes}\) 公式

\[\iint_{\Sigma}\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 = \oint_{\part\Sigma} \boldsymbol{F} \cdot \mathrm{d}\mathbf{s} \]

其中左端为标准 \(2\) - 形式，可利用外微分对右端 \(1\) - 形式求导获得。

五、广义Stokes公式

广义 \(\text{Stokes}\) 公式：

\[\int_{\partial \Omega} \omega=\int_{\Omega} \mathrm{d} \omega \]

其中的微分算符为外微分。

利用广义 \(\text{Stokes}\) 公式，可以导出第四节中全部公式。