跳转到内容
← 返回定理
向量分析中等约 17 分钟阅读

散度定理

Divergence Theorem

gauss·1813
向量分析散度高斯电磁学

定理陈述

散度定理（Divergence Theorem），也称为高斯定理（Gauss's Theorem），建立了体积分与面积分之间的关系。

设 $V$ 是三维空间中的有界闭区域，边界 ∂V\partial V 分段光滑且取外法向。设 F\mathbf{F} 是在包含 $V$ 的开集上连续可微的向量场，则

∭V(∇⋅F) dV=∯∂VF⋅n dS\iiint_V (\nabla \cdot \mathbf{F})\,dV = \oiint_{\partial V} \mathbf{F} \cdot \mathbf{n}\,dS

其中 ∇⋅F=∂F1∂x+∂F2∂y+∂F3∂z\nabla \cdot \mathbf{F} = \frac{\partial F_1}{\partial x} + \frac{\partial F_2}{\partial y} + \frac{\partial F_3}{\partial z} 是 F\mathbf{F} 的散度，n\mathbf{n} 是 ∂V\partial V 的单位外法向量。

破除误解：它不是"三维的微积分基本定理"那么简单，但也不比它更神秘

散度定理常被介绍成"微积分基本定理的高维版本"。这个类比是对的，而且值得说透——因为说透之后，斯托克斯定理、格林定理、散度定理会塌缩成同一句话。微积分基本定理说：

∫abf′(x) dx=f(b)−f(a)\int_a^b f'(x)\,dx = f(b) - f(a)

左边是导数在区间内部的积分，右边是原函数在边界（两个端点）上的值。散度定理说的是同一件事，只是"区间"换成三维区域、"两个端点"换成一张闭曲面、"导数"换成散度。用外微分形式写出来，四条定理只有一条：

∫Mdω=∫∂Mω\int_M d\omega = \int_{\partial M} \omega

取 $M$ 为区间、ω\omega 为函数，得微积分基本定理；取二维区域，得格林定理；取三维区域与 $(n-1)$-形式，得散度定理；取曲面与 $1$-形式，得斯托克斯定理。它们不是四个定理，是一个定理的四个特例——这个统一由嘉当（É. Cartan）的外微分语言完成。

但要警惕两处过度简化：

第一，"通量 = 内部源强度"需要条件。 定理要求 F\mathbf{F} 在包含 $V$ 的开集上连续可微。若场在 $V$ 内部有奇点，等式就不再成立——最典型的例子是点电荷的库仑场 F=r^/r2\mathbf{F} = \hat{r}/r^2，它在原点之外散度恒为零，但穿过任何包围原点的闭曲面的通量都是 4π4\pi。矛盾吗？不矛盾：原点处 F\mathbf{F} 不可微，定理不适用。要把这种情形纳入，需要引入狄拉克 δ 函数（分布意义下 ∇⋅(r^/r2)=4πδ3(r)\nabla \cdot (\hat{r}/r^2) = 4\pi\delta^3(\mathbf{r})）。高斯定律里那个 QencQ_{\text{enc}} 之所以只数闭曲面内部的电荷，根源就在这里。

第二，"散度是通量密度"是结论而非定义。 散度的定义是偏导数之和，看不出任何几何含义；它等于"单位体积的净流出量"这一事实，恰恰是散度定理在无穷小体积上的推论。次序反过来讲会让人误以为定理是同义反复。

直觉理解

散度定理说的是：从一个区域流出的总量等于内部源的总强度。想象一个充满水的容器。散度 ∇⋅F\nabla \cdot \mathbf{F} 度量每个点"产生"或"吸收"水的速率。从容器表面流出的总水量等于内部所有水源的总出水量减去总吸水量。在电学中：从闭合曲面流出的电通量等于曲面包围的总电荷量（高斯定律）。这就是散度定理在物理学中最著名的应用。

散度定理也可以理解为微积分基本定理在三维的推广：FTC 说"内部的导数 = 边界上的值"，散度定理说的是同样的话，只是维度更高。

证明思路

证明概要

步骤 1：将向量场 F=(F1,F2,F3)\mathbf{F} = (F_1, F_2, F_3) 分解为三个分量，分别证明：

∭V∂F1∂x dV=∯∂VF1n1 dS\iiint_V \frac{\partial F_1}{\partial x}\,dV = \oiint_{\partial V} F_1 n_1\,dS

步骤 2：假设 $V$ 是 $x$-型区域（V={(x,y,z):a≤x≤b, (y,z)∈D(x)}V = \{(x,y,z) : a \leq x \leq b,\, (y,z) \in D(x)\}）。

步骤 3：将体积分化为累次积分：

∭V∂F1∂x dV=∬D[∫x1(y,z)x2(y,z)∂F1∂x dx]dA\iiint_V \frac{\partial F_1}{\partial x}\,dV = \iint_{D} \left[\int_{x_1(y,z)}^{x_2(y,z)} \frac{\partial F_1}{\partial x}\,dx\right]dA

=∬D[F1(x2(y,z),y,z)−F1(x1(y,z),y,z)] dA= \iint_D [F_1(x_2(y,z), y, z) - F_1(x_1(y,z), y, z)]\,dA

步骤 4：右端恰好等于 ∯∂VF1n1 dS\oiint_{\partial V} F_1 n_1\,dS——因为前表面的法向量 n1>0n_1 > 0，后表面的 n1<0n_1 < 0。

对 F2F_2 和 F3F_3 的证明类似。一般区域可以分解为有限个简单区域。

历史背景

拉格朗日早在 1762 年研究流体力学时就已写下这一结果，但没有给出证明。高斯在 1813 年研究椭球体的引力时证明了散度定理的若干特殊情形。真正给出一般定理第一个证明的是俄国数学家奥斯特罗格拉茨基（Mikhail Ostrogradsky）——他在 1826 年研究热传导时完成证明，1828 年宣读、1831 年正式发表。因此在俄罗斯和法国的数学传统中，这一定理常被称为高斯—奥斯特罗格拉茨基定理。此后格林（1828）、汤姆森（开尔文勋爵）等人也在各自的物理研究中独立或重新推导了相关公式。

在嘉当的外微分形式语言中，散度定理是斯托克斯定理 ∫Mdω=∫∂Mω\int_M d\omega = \int_{\partial M} \omega 取 ω\omega 为 $(n-1)$-形式的特殊情形。

应用

  1. 电磁学：高斯定律 ∯E⋅dS=Qenc/ε0\oiint \mathbf{E} \cdot d\mathbf{S} = Q_{\text{enc}}/\varepsilon_0——电场通量等于包围电荷
  2. 流体力学：质量守恒的连续性方程 ∂ρ∂t+∇⋅(ρv)=0\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0
  3. 热力学：热传导方程的推导——流入热量等于内部热源
  4. 引力场：高斯引力定律——引力通量与质量的关系
  5. 有限体积法：数值求解偏微分方程的基本方法

与其他定理的关系

  • 微积分基本定理：散度定理是 FTC 在三维的推广
  • 斯托克斯定理：散度定理是斯托克斯定理的特例
  • 格林定理：散度定理在二维的对应——∮F⋅n ds=∬∇⋅F dA\oint \mathbf{F} \cdot \mathbf{n}\,ds = \iint \nabla \cdot \mathbf{F}\,dA
  • 高斯定律：散度定理在电磁学中的直接应用
  • 连续性方程：物理学中质量、电荷守恒的数学表述

具体示例

设 F=(x,y,z)\mathbf{F} = (x, y, z)，$V$ 是单位球 x2+y2+z2≤1x^2 + y^2 + z^2 \leq 1。散度 ∇⋅F=1+1+1=3\nabla \cdot \mathbf{F} = 1 + 1 + 1 = 3。体积分：∭V3 dV=3⋅43π=4π\iiint_V 3\,dV = 3 \cdot \frac{4}{3}\pi = 4\pi。面积分：∯∂VF⋅n dS=∯∂V(x,y,z)⋅(x,y,z) dS=∯∂V1 dS=4π\oiint_{\partial V} \mathbf{F} \cdot \mathbf{n}\,dS = \oiint_{\partial V} (x, y, z) \cdot (x, y, z)\,dS = \oiint_{\partial V} 1\,dS = 4\pi。两边相等，验证了散度定理。再看一个非平凡的例子：设 F=(x2,y2,z2)\mathbf{F} = (x^2, y^2, z^2)，则 ∇⋅F=2(x+y+z)\nabla \cdot \mathbf{F} = 2(x + y + z)。由对称性，∭V2x dV=0\iiint_V 2x\,dV = 0（球心对称），故体积分为 0。面积分同理为 0。

从局部到整体：守恒律的数学形式

散度 ∇⋅F\nabla \cdot \mathbf{F} 度量向量场在某点是"源"还是"汇"：

  • ∇⋅F>0\nabla \cdot \mathbf{F} > 0：源，场从该点向外发散（正电荷处的电场）
  • ∇⋅F<0\nabla \cdot \mathbf{F} < 0：汇，场向该点汇聚（负电荷处的电场）
  • ∇⋅F=0\nabla \cdot \mathbf{F} = 0：无源无汇（不可压缩流体的速度场）

散度定理把这些逐点的源强度加总，等同于穿过整体边界的通量。这一"局部 ↔ 整体"的转换，正是物理学写守恒律的标准手法。以任意守恒量为例：设 ρ\rho 是它的密度（质量、电荷、能量皆可），J\mathbf{J} 是流量密度。微分形式的守恒律说

∂ρ∂t+∇⋅J=0\frac{\partial \rho}{\partial t} + \nabla \cdot \mathbf{J} = 0

在区域 $V$ 上积分并用散度定理，立刻得到积分形式：

ddt∫Vρ dV=−∮∂VJ⋅n dS\frac{d}{dt}\int_V \rho\,dV = -\oint_{\partial V} \mathbf{J} \cdot \mathbf{n}\,dS

读作：区域内总量的变化率，等于净流入量。左边是"里面有多少在变"，右边是"边界上进出多少"——两句话说的是同一件事，散度定理就是这个"同一"的内容。

麦克斯韦方程组的积分形式全部由此而来（旋度那两条用斯托克斯定理）：

微分形式积分形式用到的定理
∇⋅E=ρ/ε0\nabla \cdot \mathbf{E} = \rho/\varepsilon_0∯E⋅dS=Q/ε0\oiint \mathbf{E}\cdot d\mathbf{S} = Q/\varepsilon_0散度定理
∇⋅B=0\nabla \cdot \mathbf{B} = 0∯B⋅dS=0\oiint \mathbf{B}\cdot d\mathbf{S} = 0散度定理
∇×E=−∂B/∂t\nabla \times \mathbf{E} = -\partial\mathbf{B}/\partial t∮E⋅dl=−dΦB/dt\oint \mathbf{E}\cdot d\mathbf{l} = -d\Phi_B/dt斯托克斯定理
∇×B=μ0J+μ0ε0∂E/∂t\nabla \times \mathbf{B} = \mu_0\mathbf{J} + \mu_0\varepsilon_0\partial\mathbf{E}/\partial t∮B⋅dl=μ0I+μ0ε0 dΦE/dt\oint \mathbf{B}\cdot d\mathbf{l} = \mu_0 I + \mu_0\varepsilon_0\,d\Phi_E/dt斯托克斯定理

第二行值得单独看一眼：∇⋅B=0\nabla \cdot \mathbf{B} = 0 说磁场没有源，于是穿过任何闭曲面的磁通量恒为零——这正是"磁单极子不存在"的数学表述。散度定理让一个逐点的微分陈述，变成一个可以实际测量的整体陈述。

散度定理的高维推广

散度定理可以推广到任意维度。在 $n$ 维空间中，设 V⊂RnV \subset \mathbb{R}^n 是有界闭区域，F\mathbf{F} 是向量场，则：

∫V∇⋅F dnx=∮∂VF⋅n dn−1S\int_V \nabla \cdot \mathbf{F}\,d^n x = \oint_{\partial V} \mathbf{F} \cdot \mathbf{n}\,d^{n-1}S

在微分形式语言中，散度定理是斯托克斯定理 ∫Mdω=∫∂Mω\int_M d\omega = \int_{\partial M} \omega 的特例。取 ω=F1 dx2∧⋯∧dxn−F2 dx1∧dx3∧⋯∧dxn+⋯\omega = F_1\,dx_2 \wedge \cdots \wedge dx_n - F_2\,dx_1 \wedge dx_3 \wedge \cdots \wedge dx_n + \cdots（F\mathbf{F} 对应的 $(n-1)$-形式），则 dω=(∇⋅F) dx1∧⋯∧dxnd\omega = (\nabla \cdot \mathbf{F})\,dx_1 \wedge \cdots \wedge dx_n。嘉当的外微分形式语言将格林定理、斯托克斯定理和散度定理统一为同一个定理的不同特例——这是现代微分几何的基础。

跨域连接

  • 偏微分方程：把"区域内总量的变化率等于净流入"写在任意小的体积上，就得到守恒律型的偏微分方程。推论是：这类方程的解可以在有限时间内出现间断（激波），那里微分形式已失去意义，必须退回积分形式来定义弱解。
  • 麦克斯韦方程组：磁场散度为零，于是穿过任何闭曲面的磁通恒为零——这正是"没有磁单极子"的数学写法。推论是：这是一条可证伪陈述，一旦某闭曲面上测到非零磁通，方程本身就得改写，而不是补一个修正项了事。同样的写法也说明：电场通量只数闭曲面内部的电荷，外面的电荷再多，净贡献仍为零。
  • 电化学：电极表面的电流密度就是电荷通量，把它对时间与面积积起来，就等于参与反应的物质的量。推论是：没有副反应时总电量与产物的量必须严格成比例，测得的偏离直接指示漏电流或副反应，不必另做分析。
  • 肾脏生理：滤过率靠物质守恒来测——注入一种既不被重吸收也不被分泌的标志物，其排出速率除以血浆浓度即为滤过率。推论是：一旦该物质也被分泌，清除率就系统性高估滤过率，这正是不同标志物给出不同数值的原因。
  • 碳循环：碳收支相当于对整个地表画一个闭曲面，源与汇之差必须等于大气存量的变化率。推论是：三项中若只有两项可直接测，第三项只能由闭合残差反推，因此它的不确定度是前两项误差之和，而不是独立估计。

参考文献

  1. Carl Friedrich Gauss, Theoria attractionis corporum sphaeroidicorum (1813).
  2. H.M. Schey, Div, Grad, Curl, and All That (4th ed., 2005).
  3. Jerrold E. Marsden & Anthony J. Tromba, Vector Calculus (6th ed., 2012).
  4. 梁昆淼, 《数学物理方法》, 高等教育出版社.
  5. M. Spivak, Calculus on Manifolds (1965).

散度定理 ∭V(∇⋅F) dV=∯∂VF⋅dS\iiint_V(\nabla\cdot\mathbf{F})\,dV=\oiint_{\partial V}\mathbf{F}\cdot d\mathbf{S} 把向量场散度的体积分等同于穿过边界曲面的通量。它是高斯电场定律、流体连续性方程与热传导方程的数学基础。