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

散度定理

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 的单位外法向量。左端把散度当成体积密度加总,右端把 F⋅n\mathbf{F}\cdot\mathbf{n} 当成穿过边界的通量。等式说的是:内部源的代数和,等于净流出。要让两边有意义,曲面必须取外法向,场必须在包含 $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 的证明类似。一般区域可以分解为有限个简单区域。把体积分里对 $x$ 的偏导交给一元微积分基本定理,内层只剩下前、后两张面上 F1F_1 的差。

前表面外法向的 $x$ 分量为正,后表面为负,差值于是变成闭曲面上的 F1n1F_1 n_1。对 $y$、$z$ 各说一遍同一句话,三项相加就是散度。一般区域先切成有限块简单区域,切割面上的通量成对抵消,只留下真正的外边界。

完整演算不必在这里展开。存在性依赖的是区域可剖分,以及 F\mathbf{F} 的连续可微;边界有有限条棱与有限个角时,分段光滑已经够用。场若有奇点,必须先挖掉奇点再取极限,否则前提不成立。库仑场穿过包围原点的球面通量是 4π4\pi、原点外散度却为零,缺的不是证明技巧,缺的是定理假设。

历史背景

拉格朗日早在 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
  • 高斯定律:散度定理在电磁学中的直接应用
  • 连续性方程:物理学中质量、电荷守恒的数学表述

经典斯托克斯定理把曲面上的旋度通量换成边界环量,散度定理把体积内的散度换成闭曲面通量。二者不是互相抄写的一对公式,而是广义斯托克斯公式在不同维数、不同形式上的特化。格林定理对应平面有界区域:把散度定理压到二维,或把曲线上的斯托克斯再降一维,都会落到格林。

称为高斯公式,强调的是高斯处理引力势时用过的特殊情形。奥斯特罗格拉茨基给出一般向量场的证明,所以文献里也写高斯–奥斯特罗格拉茨基定理。名字记录谱系,陈述的是同一条恒等式。外微分语言里它们塌成 ∫Mdω=∫∂Mω\int_M d\omega=\int_{\partial M}\omega,并不是说三维向量分析可以被一笔勾销:计算通量时,仍要核对定向、奇点与光滑性。

具体示例

设 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=(x,y,z)\mathbf{F}=(x,y,z) 处处光滑,散度为常数,通量退成表面积。换一个在原点爆破的场,就不能再把奇点包在里面假装定理仍可用。数值上两边相等,只能验证假设满足的情形;它不能代替对一般区域的证明,也不能用来给奇点场“开特例后门”。

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

散度 ∇⋅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。嘉当的外微分形式语言将格林定理、斯托克斯定理和散度定理统一为同一个定理的不同特例——这是现代微分几何的基础。

维数升高之后,改变的是“边界”的含义:$n$ 维体积的边界是 $(n-1)$ 维超曲面,通量仍是 F⋅n\mathbf{F}\cdot\mathbf{n} 在边界上的积分。外微分把这项写成 $(n-1)$-形式的积分,于是格林、斯托克斯、散度看起来像三条定理,写出来却是同一条。具体计算仍要选坐标系、核对定向;统一语言消除的是记忆负担,不是分析条件。

跨域连接

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

参考文献

  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).
  6. Tom M. Apostol, Calculus, Volume 2 (2nd ed., Wiley, 1969).

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