散度定理的推导过程-散度定理推导过程

在物理学中，它直观地描述了流体或电流在闭合区域内的“源点”效应；在数学分析中，它架起了微积分从点集论到拓扑学的桥梁，使得高维空间中的推导变得严谨而优美。

尽管从直观理解到纯数学证明的跨越需要深厚的功底，但通过系统梳理步骤，我们可以清晰地看到这一奇偶函数关系的内在逻辑。本文将从界域职考网xinlishi.cc专家视角出发，结合经典案例，为您详细拆解散度定理的推导过程。

核心散度定理，格林公式，斯托克斯公式，高斯公式

当我们观察一个静止的不可压缩流体在容器内流动时，单位时间内流出的流体总量，应该等于该容器内各点流速的总和。这个总和就是体积分。而单位时间内流入容器外部流体的总量，则是穿过容器表面的积分。两者相等，即为散度定理的雏形。

为了使其更具普遍性，我们将研究对象推广为任意二维向量场，不再局限于流体。假设有一个定义在区域 $D$ 上的二维向量场 $mathbf{F} = (P, Q)$，其中 $P$ 和 $Q$ 是定义在 $D$ 上的光滑函数，$D$ 为平面上的有界区域。现在，我们考虑一个包围区域 $D$ 的简单闭曲线 $C$，方向为逆时针方向（正向）。

假设曲线 $C$ 将平面区域 $D$ 划分为两部分：一部分是内部区域 $D$，另一部分是外部区域 $D^c$。如果我们从 $C$ 的起点沿内部区域 $D$ 走，再沿曲线 $C$ 回到终点，最后从终点沿外部区域 $D^c$ 的边界反走，就形成了一个完全封闭的回路。在这个封闭回路内，向量场 $mathbf{F}$ 的散度 $nabla cdot mathbf{F}$ 处处存在且连续。

根据散度定理（或称高斯 - 斯托克斯定理在二维的对应形式），该封闭回路上的向量场的通量与回路围成的区域 $D$ 上的散度积分相等。用数学符号表示，就是：

$int_{C} mathbf{F} cdot dmathbf{r} = iint_{D} (nabla cdot mathbf{F}) , dA$

这里，左边是向量场沿边界 $C$ 的线积分，右边是散度在区域 $D$ 上的体积分。

这个等式揭示了“场”与“源”之间的数量守恒关系。左边关注的是边界上的流动情况，右边关注的是体内的生成情况。当我们把曲线 $C$ 的外边界 $C^text{out}$ 和内部区域 $D$ 的外边界 $C^text{in}$ 看作同一顺逆方向时，若曲线方向为顺时针，则 $C^text{out}$ 为正方向，$C^text{in}$ 为负方向，此时 $C$ 的积分等于 $C^text{out}$ 的全积分减去 $C^text{in}$ 的全积分，最终结论即为上述等式。

我们在数学化过程中遇到了一大难题：如何计算这个封闭曲线 $C$ 上的线积分？如果我们能找到一条从点 $(x_1, y_1)$ 到点 $(x_2, y_2)$ 的替代路径，使得积分值相同，那么原曲线 $C$ 上的值也就随之确定了。

由此，我们需要建立一个更强大的工具，能够计算任意光滑路径的积分。这个工具就是格林公式（Green's Theorem）。格林公式建立了平面上的曲线积分与区域上的二重积分之间的深刻联系，它将原本难以计算的线积分转化为我们熟悉的重积分形式。
因此，格林公式不仅是连接微积分基本定理的桥梁，更是推导散度定理的关键枢纽。

我们将引入格林公式，通过层层推导，最终完成散度定理的证明。

回顾格林公式的内容：在平面有界区域 $D$ 上，若函数 $mathcal{P}$ 和 $mathcal{Q}$ 具有一阶连续偏导数，且 $C$ 是 $D$ 的周围简单闭曲线，方向为正向，则有：

$oint_{C} (P , dx + Q , dy) = iint_{D} left( frac{partial Q}{partial x} - frac{partial P}{partial y} right) , dA$

注意，格林公式中的积分形式是 $oint (P dx + Q dy)$，这实际上就是向量场 $mathbf{F} = (P, Q)$ 沿曲线 $C$ 的线积分，即 $oint_{C} mathbf{F} cdot dmathbf{r}$。

既然我们已经知道曲线积分等于 $iint_{D} (Q_x - P_y) , dA$，那么直接代入散度定理等式 $int_{C} mathbf{F} cdot dmathbf{r} = iint_{D} (nabla cdot mathbf{F}) , dA$ 中似乎已经完成了证明？不，这只是个等式。我们需要证明两边代表的量在本质上是一致的，或者更准确地说，我们需要利用格林公式来消除线积分中的变量，从而显式地写出散度函数的积分形式。

让我们重新审视等式右边。我们将线积分展开为累次积分形式：

$oint_{C} P , dx = int_{x_1}^{x_2} P , dx + int_{y_1}^{y_2} P , dy$

注意，这里的变量 $x$ 和 $y$ 是在标量域上变化的，即 $int_{x_1}^{x_2} P(x, y) , dx$ 表示对 $x$ 积分。同理，$Q$ 对 $y$ 的积分写为 $int_{y_1}^{y_2} Q(x, y) , dy$。

现在，我们将这两个项代入线积分表达式：

$oint_{C} P , dx + Q , dy = int_{x_1}^{x_2} P , dx + int_{y_1}^{y_2} Q , dy$

这里的关键在于，积分变量 $x$ 和 $y$ 的范围是由曲线 $C$ 决定的。为了消除这个依赖，我们引入一条特殊的辅助路径来构造一个闭合回路。假设曲线 $C$ 的参数方程为 $mathbf{r}(t) = (x(t), y(t))$，其中 $t$ 从 $a$ 变到 $b$，且 $a < b$。

我们可以构造一条从起点 $(x_1, y_1)$ 到终点 $(x_2, y_2)$ 的辅助路径 $C_0$。这条路径 $C_0$ 由两段组成：

1.一段是沿 $C$ 从起点到终点的路径；

2.另一段是水平线段 $L_1$，从终点 $(x_2, y_2)$ 水平向左移动到 $(x_1, y_2)$，此时 $y$ 保持 $y_2$ 不变；

3.最后一段是垂直线段 $L_2$，从 $(x_1, y_2)$ 垂直向下移动到起点 $(x_1, y_1)$，此时 $x$ 保持 $x_1$ 不变。

现在，我们将这三段路径记为 $C$（逆时针方向）、$L_1$（向右）和 $L_2$（向上），这样这三段路径共同构成了一个闭合回路 $C_{closed}$。根据格林公式，对于向量场 $mathbf{F}=(P,Q)$ 在闭合回路上的积分，等于区域 $D_{closed}$ 上散度的积分：

$oint_{C_{closed}} mathbf{F} cdot dmathbf{r} = iint_{D_{closed}} (nabla cdot mathbf{F}) , dA$

而闭合回路 $C_{closed}$ 上的积分可以分解为三部分：

$oint_{C_{closed}} mathbf{F} cdot dmathbf{r} = int_{C} mathbf{F} cdot dmathbf{r} + int_{L_1} mathbf{F} cdot dmathbf{r} + int_{L_2} mathbf{F} cdot dmathbf{r}$

让我们分别计算这三段的积分：

1.对于线段 $L_1$：从 $(x_2, y_2)$ 到 $(x_1, y_2)$，此时 $y = y_2$ 为常数，$dy = 0$。积分变为 $int_{x_2}^{x_1} P(x, y_2) , dx$。

2.对于线段 $L_2$：从 $(x_1, y_2)$ 到 $(x_1, y_1)$，此时 $x = x_1$ 为常数，$dx = 0$。积分变为 $int_{y_2}^{y_1} Q(x_1, y) , dy$。

因此，闭合回路的积分整体变成了：

$oint_{C_{closed}} mathbf{F} cdot dmathbf{r} = int_{C} (P , dx + Q , dy) + int_{x_2}^{x_1} P(x, y_2) , dx + int_{y_2}^{y_1} Q(x_1, y) , dy$

现在，我们利用格林公式对闭合回路 $C_{closed}$ 进行积分计算。区域 $D_{closed}$ 是 $D$ 与 $L_1$、$L_2$ 围成的区域。根据格林公式：

$iint_{D_{closed}} (frac{partial Q}{partial x} - frac{partial P}{partial y}) , dA = iint_{D} (frac{partial Q}{partial x} - frac{partial P}{partial y}) , dA + iint_{D_{L1}} (frac{partial Q}{partial x} - frac{partial P}{partial y}) , dA + iint_{D_{L2}} (frac{partial Q}{partial x} - frac{partial P}{partial y}) , dA$

其中 $D_{L1}$ 是由 $L_1$ 和 $C$ 围成的区域，$D_{L2}$ 是由 $L_2$ 和 $C$ 围成的区域。

我们首先计算 $D_{L1}$ 上的积分。在 $D_{L1}$ 上，$y$ 保持不变（等于 $y_2$），$dy = 0$，故 $dA = dx , dy = dx$。于是积分变为：

$iint_{D_{L1}} frac{partial Q}{partial x} , dx = int_{x_2}^{x_1} left[ int_{y_1}^{y_2} frac{partial Q}{partial x} , dy right] dx$

交换积分次序：$int_{y_1}^{y_2} left[ int_{x_2}^{x_1} frac{partial Q}{partial x} , dx right] dy$

根据分部积分法，$int frac{partial Q}{partial x} , dx = Q$，所以内层积分为 $[Q(x, y_2) - Q(x_2, y_2)]$。

接着计算 $D_{L2}$ 上的积分。在 $D_{L2}$ 上，$x$ 保持不变（等于 $x_1$），$dx = 0$，故 $dA = dy , dx = dy$。于是积分变为：

$iint_{D_{L2}} frac{partial P}{partial y} , dy = int_{y_1}^{y_2} left[ int_{x_1}^{x_2} frac{partial P}{partial y} , dx right] dy$

这里有一个符号需要注意。格林公式中的二阶偏导数项是 $-frac{partial P}{partial y}$，但在计算 $D_{L2}$ 时，我们使用的是 $-frac{partial P}{partial y}$，而 $P$ 对 $x$ 的导数是 $frac{partial P}{partial x}$。由于 $x$ 固定，$frac{partial P}{partial x} = 0$？不对，格林公式是 $iint (nabla cdot mathbf{F}) dA = iint (Q_x - P_y) dx dy$。在 $D_{L2}$ 上，$y$ 变化，$dx=0$，所以 $dA=dy$。项是 $Q_x - P_y$。由于 $P$ 不随 $y$ 变，$P_y=0$，项为 $Q_x$？

让我们重新梳理格林公式的应用。格林公式 $oint_{C} P dx + Q dy = iint_{D} (Q_x - P_y) dx dy$。在 $D_{L2}$ 上，$x=x_1$，$dx=0$，$y in [y_1, y_2]$。$iint_{D_{L2}} (Q_x - P_y) dx dy = int_{y_1}^{y_2} left( int_{x_1}^{x_2} (Q_x - P_y) dx right) dy$。这里的内层积分 $int_{x_1}^{x_2} Q_x dx = Q(x_2, y) - Q(x_1, y) = P(x_2, y) - P(x_1, y)$？这似乎不对，因为我们在计算 $iint (Q_x - P_y)$。

让我们回到微积分基本定理。$int_{x_1}^{x_2} frac{partial Q}{partial x} dx = Q(x_2, y) - Q(x_1, y)$。同理 $int_{x_1}^{x_2} frac{partial P}{partial y} dx = frac{partial P}{partial y}$ 对 $x$ 积分？不对，$P$ 是 $x$ 的函数，$P_y$ 是对 $y$ 的偏导。$int_{x_1}^{x_2} P_y(y) dx$ 是常数 $P_y$ 乘以 $(x_2-x_1)$？

正确推导如下：

1.计算 $Q_x$ 的积分部分：$int_{x_1}^{x_2} frac{partial Q}{partial x} dx$。这里 $Q$ 视为 $y$ 的函数。根据微积分基本定理，$int_{x_1}^{x_2} frac{partial Q}{partial x} dx = Q(x_2, y) - Q(x_1, y)$。但这包含了 $Q_y$ 项吗？不，这是 $Q(x_2, y) - Q(x_1, y)$。

实际上，我们直接对 $x$ 积分：$int_{x_1}^{x_2} frac{partial Q}{partial x} dx = left[ Q right]_{x_1}^{x_2} = Q(x_2, y) - Q(x_1, y)$。但这还没涉及 $P$。

让我们重新看格林公式右边的项：$Q_x - P_y$。

计算 $iint_{D_{L1}} (Q_x - P_y) dx dy = int_{y_1}^{y_2} left[ int_{x_2}^{x_1} Q_x dx - P_y right] dy$。注意积分限 $x_2$ 到 $x_1$ 是负的，或者写成 $int_{x_1}^{x_2} -Q_x dx - P_y dx$。$int_{x_1}^{x_2} Q_x dx = Q(x_2, y) - Q(x_1, y)$。$int_{x_1}^{x_2} P_y dx$？不对，$P_y$ 是 $y$ 的函数，对 $x$ 积分是 $P_y(x_2, y) - P_y(x_1, y)$ 吗？不对，$Q$ 和 $P$ 都是 $x,y$ 的函数。

纠正：$iint_{D_{L1