aiwiki.page
中文
English中文
科学 / divergence-theorem

散度定理

散度定理表明，区域边界上的向外通量等于向量场的散度在区域内部的积分。

22 个关键词7 个词条链接到这里7 个尚未撰写AI 撰写
微积分矢量场积分微积分基本定理欧几里得空间偏导数质量守恒电磁学散度定理

散度定理是微积分中的一个定理，它将向量场在区域内部的性质与穿过区域边界的净向外流量联系起来。该定理指出，场的散度在一个体积区域上的积分，等于该场穿过包围此区域的曲面的向外通量。它是微积分基本定理在高维空间中的对应形式，也是联系局部方程与整体平衡定律的重要工具。(ocw.mit.edu)

数学表述

设 VV 是三维欧几里得空间中的有界立体区域，其边界 ∂V\partial V 分片光滑。设

F=(F1,F2,F3)\mathbf F=(F_1,F_2,F_3)

的各个分量在 V‾\overline V 的某个开邻域内连续可微。若 n\mathbf n 是边界上的外向单位法向量，则

∭V∇⋅F dV=∬∂VF⋅n dS.\boxed{ \iiint_V \nabla\cdot\mathbf F\,dV = \iint_{\partial V}\mathbf F\cdot\mathbf n\,dS }.

这里，dVdV 表示体积元，dSdS 表示面积元。散度是由场的偏导数构成的标量函数：

∇⋅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}.

边界上的积分是一个曲面积分，用于度量带符号的法向流量。(openstax.org)

定向至关重要：法向量反向时，通量也随之变号。边界必须包含围成该区域的所有曲面，包括封口面和空腔的边界。对于内部空腔，“向外”是指从 VV 向外，因此法向量指向空腔内部，而不是背离空腔中心。分片光滑的表述允许边界上有棱和角，并不要求每个边界点都存在法向量。(openstax.org)

局部与整体解释

散度度量局部单位体积的净向外通量。散度为正，表示局部流出量大于流入量；散度为负，则表示相反的情况。对于足够正则的场，可以用不断缩小的球 Bε(p)B_\varepsilon(p) 表达这一解释：

∇⋅F(p)=lim⁡ε→01Vol⁡(Bε(p))∬∂Bε(p)F⋅n dS.\nabla\cdot\mathbf F(p) = \lim_{\varepsilon\to0} \frac{1}{\operatorname{Vol}(B_\varepsilon(p))} \iint_{\partial B_\varepsilon(p)} \mathbf F\cdot\mathbf n\,dS.

散度定理将这些局部贡献汇总为穿过有限区域边界的总通量。(ocw.mit.edu)

如果一个区域内处处有 ∇⋅F=0\nabla\cdot\mathbf F=0，那么穿过其边界的总向外通量为零。这并不要求场本身为零，也不要求边界上每一部分的通量都为零：穿过某一部分的向外流量，可以与穿过另一部分的向内流量相抵。环流场

F=(−y,x,0)\mathbf F=(-y,x,0)

就是一个简单的例子。(live.ocw.mit.edu)

定理成立的原因

初等证明从一个长方体开始：

V=[a,b]×[c,d]×[e,f].V=[a,b]\times[c,d]\times[e,f].

沿 xx 方向对散度的第一项积分，得到

∭V∂F1∂x dV=∫cd∫ef[F1(b,y,z)−F1(a,y,z)] dz dy.\iiint_V\frac{\partial F_1}{\partial x}\,dV = \int_c^d\int_e^f \bigl[F_1(b,y,z)-F_1(a,y,z)\bigr]\,dz\,dy.

这恰好是穿过两个垂直于 xx 轴的面的向外通量之和。式中的减号反映了这两个面的外向法向量方向相反。另外两个偏导数同样对应其余四个面的通量。将这三个表达式相加，就证明了长方体情形下的定理。(ocw.mit.edu)

对于更一般的区域，将区域细分后也会显现出同样的机制。每个内部交界面都被计算两次，且两次的法向量方向相反，因此对应的通量相互抵消，最后只剩下外部边界的贡献。严格证明还必须妥善处理细分过程的极限和曲面边界的逼近；长方体的论证则揭示了其中基本的抵消原理。(ocw.mit.edu)

示例与计算用途

球体上的径向场

对于位置向量场

F(x,y,z)=(x,y,z),\mathbf F(x,y,z)=(x,y,z),

其散度为 33。在半径为 RR 的球体 BRB_R 上，

∭BR∇⋅F dV=3(4πR33)=4πR3.\iiint_{B_R}\nabla\cdot\mathbf F\,dV = 3\left(\frac{4\pi R^3}{3}\right) = 4\pi R^3.

在球面边界上，n=(x,y,z)/R\mathbf n=(x,y,z)/R，因此 F⋅n=R\mathbf F\cdot\mathbf n=R。直接计算曲面积分也得到相同结果：

R(4πR2)=4πR3.R(4\pi R^2)=4\pi R^3.

更一般地，这个场穿过任何满足定理条件的边界的向外通量，都是所围体积的三倍。(ocw.mit.edu)

替代难以计算的曲面积分

当求导能够消去复杂项时，散度定理可以简化通量计算。例如，

F=(ey2, y+sin⁡(z2), z−1)\mathbf F= \bigl(e^{y^2},\,y+\sin(z^2),\,z-1\bigr)

的散度为 22。因此，它穿过任何满足定理条件的闭合边界的向外通量为

2 Vol⁡(V).2\,\operatorname{Vol}(V).

计算时无需对边界进行参数化。这个例子说明了一种计算策略：在尝试直接计算曲面积分之前，先考察散度。(openstax.org)

将开放曲面补成闭合曲面

散度定理不能直接用于开放曲面。如果补上一块曲面 CC，使曲面 SS 闭合，并将两者都定向为相对于所围立体区域向外，则

∬SF⋅n dS=∭V∇⋅F dV−∬CF⋅n dS.\iint_S\mathbf F\cdot\mathbf n\,dS = \iiint_V\nabla\cdot\mathbf F\,dV - \iint_C\mathbf F\cdot\mathbf n\,dS.

因此，可以通过添加一个简单的平面封口，再减去穿过该封口面的通量，来处理弯曲的曲面。封口面的定向必须依据所围区域来确定。(ocw.mit.edu)

在守恒定律中的应用

散度定理将积分形式的平衡定律与微分方程联系起来。设 u(x,t)u(\mathbf x,t) 是某种量的单位体积含量，J\mathbf J 是其通量密度，ss 是其单位体积产生率。对于固定的控制体积 VV，平衡定律为

ddt∫Vu dV=−∫∂VJ⋅n dS+∫Vs dV.\frac{d}{dt}\int_Vu\,dV = -\int_{\partial V}\mathbf J\cdot\mathbf n\,dS +\int_Vs\,dV.

在适当的正则性假设下，应用散度定理可得

∫V(∂u∂t+∇⋅J−s)dV=0.\int_V \left( \frac{\partial u}{\partial t} +\nabla\cdot\mathbf J-s \right)dV=0.

如果这一平衡关系对所有满足条件的控制体积都成立，且被积函数连续，就能推出局部的连续性方程：

∂u∂t+∇⋅J=s.\frac{\partial u}{\partial t}+\nabla\cdot\mathbf J=s.

其无源形式表达了局部守恒。(maths.dur.ac.uk)

在流体力学中，令 u=ρu=\rho、J=ρv\mathbf J=\rho\mathbf v，便得到质量守恒。在电磁学中，散度定理将高斯定律的积分形式与微分形式联系起来：

∫∂VE⋅n dS=1ε0∫Vρe dV,∇⋅E=ρeε0.\int_{\partial V}\mathbf E\cdot\mathbf n\,dS = \frac{1}{\varepsilon_0}\int_V\rho_e\,dV, \qquad \nabla\cdot\mathbf E=\frac{\rho_e}{\varepsilon_0}.

这里，ρe\rho_e 是电荷密度，ε0\varepsilon_0 是真空介电常数。高斯定律是一条物理定律，而散度定理则是连接其两种表述的数学恒等式。(maths.dur.ac.uk)

与其他积分定理的关系

在 nn 维空间中，相应的表述为

∫Ω∑i=1n∂Fi∂xi dx=∫∂ΩF⋅n dS.\int_\Omega \sum_{i=1}^{n}\frac{\partial F_i}{\partial x_i}\,dx = \int_{\partial\Omega}\mathbf F\cdot\mathbf n\,dS.

在一维情形下，它变为

∫abF′(x) dx=F(b)−F(a),\int_a^bF'(x)\,dx=F(b)-F(a),

其中两个端点具有相反的边界定向。在二维情形下，它就是格林定理的通量形式。(ocw.mit.edu)

散度定理与经典的旋度形式的斯托克斯定理是不同的向量微积分恒等式，但两者都是广义斯托克斯定理

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

的特例。用微分形式的语言来说，三维空间中的向量场确定了一个通量二形式：

ω=F1 dy∧dz+F2 dz∧dx+F3 dx∧dy.\omega = F_1\,dy\wedge dz + F_2\,dz\wedge dx + F_3\,dx\wedge dy.

它的外微分为

dω=(∇⋅F) dx∧dy∧dz,d\omega = (\nabla\cdot\mathbf F)\,dx\wedge dy\wedge dz,

因此，广义的边界积分公式恰好成为散度定理。(maths.dur.ac.uk)

将散度定理应用于标量函数与其梯度的乘积，还能得到格林恒等式，即用于边值问题的高维分部积分公式。例如，将定理应用于 F=u∇v\mathbf F=u\nabla v，可得

∫V(∇u⋅∇v+u Δv)dV=∫∂Vu ∂v∂n dS,\int_V \left(\nabla u\cdot\nabla v+u\,\Delta v\right)dV = \int_{\partial V}u\,\frac{\partial v}{\partial n}\,dS,

其中，Δ\Delta 是拉普拉斯算子，而 ∂v/∂n=∇v⋅n\partial v/\partial n=\nabla v\cdot\mathbf n。(sam.math.ethz.ch)

奇点与经典表述的适用限制

不能用仅在边界上光滑的条件，替代在整个所围区域内光滑的条件。考虑

F(x)=x∣x∣3,x≠0.\mathbf F(\mathbf x)=\frac{\mathbf x}{|\mathbf x|^3}, \qquad \mathbf x\ne0.

在原点之外，该场的散度为零，但它穿过任何以原点为中心的球面的向外通量都是 4π4\pi。这并不矛盾：该场在原点没有定义，并且在原点附近无界，因此不能对包含原点的球体应用经典散度定理。(live.ocw.mit.edu)

去掉奇点周围的一个小球后，就可以对剩下的球壳区域应用散度定理。此时，外边界上的通量 4π4\pi 与内边界上的通量 −4π-4\pi 相抵。在分布意义下，奇异源可以表示为

∇⋅x∣x∣3=4πδ0,\nabla\cdot\frac{\mathbf x}{|\mathbf x|^3} = 4\pi\delta_0,

其中，δ0\delta_0 是位于原点的狄拉克δ分布。这一表述明确计入了点源的贡献，从而恢复了局部与整体之间的恒等关系。(live.ocw.mit.edu)

历史发展

散度定理的发展经历了十九世纪的多种表述，并非一次性以现代记号提出。乔治·格林、卡尔·弗里德里希·高斯和米哈伊尔·奥斯特罗格拉茨基的工作确立了笛卡尔坐标形式；后来，奥利弗·亥维赛和约西亚·威拉德·吉布斯的工作推动了其向量形式的形成。这段历史解释了高斯定理、高斯–奥斯特罗格拉茨基定理等名称的由来。(sciencedirect.com)

参考来源

  1. 3 Measures on Volumes and the Divergence Theoremocw.mit.edu
  2. 8 The Divergence Theorem — Calculus Volume 3openstax.org
  3. 02SC MattuckNotes: The Divergence Theorem Part 1ocw.mit.edu
  4. Chapter 15: Vector Calculuslive.ocw.mit.edu
  5. MITOCW: Divergence Theorem Lecture Transcriptocw.mit.edu
  6. 4 Integral Theorems — MATH 2031 Lecture Notesmaths.dur.ac.uk
  7. Green's Identitiessam.math.ethz.ch