散度、高斯散度定理

贡献者： addis

预备知识　全微分，曲面积分，流密度

　　我们在矢量场中取一个闭合曲面 $S$ ，其内部空间记为 $V$ 。以向外为正方向，矢量场 $F (r)$ 在闭合曲面的通量 $Φ$ 可以用以下面积分表示，积分范围默认为 $S$

\begin{matrix} (1) & Φ = \oint F \cdot d s . \end{matrix}

　　现在我们把该曲面以其内部一点 $r$ 为中心按比例不断缩小，若通量与体积 $V$ 的比值存在极限，就把该极限叫做该点的散度（divergence），用 $\nabla \cdot$ 算符¹记为（下文将介绍）

\begin{matrix} (2) & \nabla \cdot F (r) \equiv lim_{V \to 0} \frac{Φ}{V} . \end{matrix}

若场的分布连续且光滑，则该极限处处存在且与曲面的形状无关（证明略），我们就得到了矢量场的散度场（注意是标量场）。

例 1　匀速水流

　　假设密度不变的水以匀速流动，质量的流密度场 $j (r)$ 为恒定场。对于任何一个闭合曲面，流入的流量（负值）和流出的流量（正值）相等，总通量为零。所以该场的散度处处为零。

例 2　变速水流

　　在例 1 中，流密度场随 $x$ 坐标线性增加， $j = (j_{0} + α x) \hat{x}$ （ $α$ 为常数），那么取一个边长为 $h$ 的立方体表面作为闭合曲面，从左侧的流量为 $Φ_{L} = - (j_{0} + α x_{0}) h^{2}$ ，右侧的流量为 $Φ_{R} = [j_{0} + α (x_{0} + h)] h^{2}$ ，其余四个面与 $j$ 平行，没有流量。闭合曲面的总流量为 $Φ = α h^{3} = α V$ 。根据定义，水流的散度处处为 $α$ 。分析可发现该水流中单位体积单位时间必然会凭空产生质量为 $α$ 的水（虽然实际中不可能）。所以散度也叫源密度（source density）。

1. 直角坐标系中的散度

　　若在直角坐标系中给出矢量场

\begin{matrix} (3) & F (x, y, z) = F_{x} (x, y, z) \hat{x} + F_{y} (x, y, z) \hat{y} + F_{z} (x, y, z) \hat{z}, \end{matrix}

令闭合曲面为立方体

[x, y, z]

[x + h, y + h, z + h]

的表面。先来考虑

x

方向两个正方形的通量

Φ_{x}

，在点

r (x, y, z)

附近对

F_{x} (r)

使用全微分近似得（为简便书写，以下的函数值和偏导都默认在

r

处取值）

\begin{matrix} (4) & F_{x} (x + x^{'}, y + y^{'}, z + z^{'}) \approx F_{x} + \frac{\partial F_{x}}{\partial x} x^{'} + \frac{\partial F_{x}}{\partial y} y^{'} + \frac{\partial F_{x}}{\partial z} z^{'} . \end{matrix}

由于只有

x

方向的场分量对

Φ_{x}

有贡献，

\begin{matrix} (5) & \begin{aligned} Φ_{x} & \approx \int_{0}^{h} \int_{0}^{h} d y^{'} d z^{'} \times \\ [(F_{x} + \frac{\partial F_{x}}{\partial x} h + \frac{\partial F_{x}}{\partial y} y^{'} + \frac{\partial F_{x}}{\partial z} z^{'}) - (F_{x} + \frac{\partial F_{x}}{\partial x} 0 + \frac{\partial F_{x}}{\partial y} y^{'} + \frac{\partial F_{x}}{\partial z} z^{'})] \\ = \frac{\partial F_{x}}{\partial x} h \int_{0}^{h} \int_{0}^{h} d y^{'} d z^{'} = \frac{\partial F_{x}}{\partial x} h^{3} = \frac{\partial F_{x}}{\partial x} V . \end{aligned} \end{matrix}

同理可以得到另外四个正方形的通量。六个正方形的总通量为

\begin{matrix} (6) & Φ \approx (\frac{\partial F_{x}}{\partial x} + \frac{\partial F_{y}}{\partial y} + \frac{\partial F_{z}}{\partial z}) V . \end{matrix}

根据定义式 2 ，可得直角坐标中的散度公式

\begin{matrix} (7) & \nabla \cdot F = lim_{V \to 0} \frac{Φ}{V} = \frac{\partial F_{x}}{\partial x} + \frac{\partial F_{y}}{\partial y} + \frac{\partial F_{z}}{\partial z} . \end{matrix}

从形式上，我们可以引入一个

\nabla

算符，在直角坐标系中的形式为

\begin{matrix} (8) & \nabla = \hat{x} \frac{\partial}{\partial x} + \hat{y} \frac{\partial}{\partial y} + \hat{z} \frac{\partial}{\partial z}, \end{matrix}

那么

\nabla \cdot F (r)

从形式上可以看做矢量算符

\nabla

与某点场矢量

F (r)

的 “内积”。根据式 7 ，显然散度是一个线性算符，即多个矢量场的线性组合的散度等于它们分别求散度再线性组合

\begin{matrix} (9) & \nabla \cdot [C_{1} F_{1} (r) + C_{2} F_{2} (r) + . . .] = C_{1} \nabla \cdot F_{1} (r) + C_{2} \nabla \cdot F_{2} (r) + \dots \end{matrix}

例 3　质点引力场的散度

　　令质点 $m$ 在坐标原点，则它的引力场（式 2 ）为

\begin{matrix} (10) & g (r) = - \frac{G m}{r^{3}} r . \end{matrix}

我们在直角坐标系中计算该场的散度。直角坐标系中，有

r = x \hat{x} + y \hat{y} + z \hat{z}

和

r = \sqrt{x^{2} + y^{2} + z^{2}}

，代入上式再求散度，得

\begin{matrix} (11) & \begin{aligned} \nabla \cdot g & = \frac{G m}{(x^{2} + y^{2} + z^{2})^{5 / 2}} [(2 x^{2} - y^{2} - z^{2}) + (2 y^{2} - z^{2} - x^{2}) + (2 z^{2} - x^{2} - y^{2})] \\ = 0, \end{aligned} \end{matrix}

可见引力场的散度为 0。然而需要注意的是，在坐标原点处式 10 的各个方向的偏导都不存在，所以不能用该公式计算散度。以上的结论只适用于原点之外的点。

2. 散度定理

　　我们来考虑一个有限大的闭合曲面并计算通量 $Φ$ 。我们先把曲面内的空间划分成许多体积足够小的微元，第 $i$ 个的体积微元为 $V_{i}$ ，通量为 $Φ_{i} \approx \nabla \cdot F (r_{i}) V_{i}$ 。现在来证明所有小曲面的通量之和等于大曲面的通量。

图 1：证明所有小曲面的通量之和等于大曲面的通量

　　图 1 中所有微元的曲面可划分为两部分，一是相邻两个小曲面的边界（红色），二是小曲面与大曲面重合的部分（黑色）。前者产生的通量之和为零，因为这些边界都是由正方向相反的两块小曲面重合而成，它们产生的通量等大反向，互相抵消。后者产生的通量等于大曲面的通量，这是因为每块黑色边界都是由正方向相同的小曲面和大曲面重合而成，产生的通量等大同向。所以总通量等于

\begin{matrix} (12) & Φ = \sum_{i} Φ_{i} \approx \sum_{i} \nabla \cdot F (r_{i}) V_{i} . \end{matrix}

令微元趋近无穷小，上面的求和变为定积分（积分范围默认为

V

）

\begin{matrix} (13) & \oint F \cdot d s = \int \nabla \cdot F (r) d V . \end{matrix}

这就是高斯散度定理，也简称为高斯定理（Gauss's theorem）或者散度定理（divergence theorem）就是，矢量场在任意闭合曲面的通量等于矢量场的散度在曲面所围空间的体积分。

未完成：给出一些场和曲面，证明散度定理成立

1. ^ 符号 $\nabla$ 的名字为 nabla，作为算符时读作 del，一些教材也会在上方加矢量箭头，原因见下文。

散度、高斯散度定理

例 1 匀速水流

例 2 变速水流

1. 直角坐标系中的散度

例 3 质点引力场的散度

2. 散度定理

例 1　匀速水流

例 2　变速水流

例 3　质点引力场的散度