Skip to content

散度定理

Divergence theorem · Gauss theorem · Gauss–Ostrogradsky theorem

区域内部散度的体积分等于其边界上向外法向通量的积分。

形式陈述

ΩRn 是有界开集,其边界 Ω 至少分片 C1ν 表示几乎处处定义的向外单位法向量。若向量场 F=(F1,,Fn)Ω 的某邻域上属于 C1,则

ΩdivFdx=ΩFνdS,

其中

divF=j=1nFjxj.

左侧按 n 维体积积分,右侧按边界的 (n1) 维面积积分。边界定向固定为向外;若改用向内法向量,右侧必须整体变号。

该结论是广义 Stokes 定理的 Euclidean 特例。令标准体积形式为

ω=dx1dxn,

并以 ιFω 表示向量场 F 与体积形式的内乘,则

d(ιFω)=(divF)ω.

应用 Ωdα=Ωα,边界上的拉回恰为 (Fν)dS

更弱的版本允许 Lipschitz 边界与 Sobolev 向量场,并用迹算子或测度论法向迹解释边界通量;正则性条件随版本改变,不能在不说明积分与边界意义时直接删去。

直觉

divF(x) 测量点 x 附近每单位体积的净产生率。把区域切成许多小控制体后,相邻小块共享面的流出与流入符号相反,内部通量两两抵消,最后只剩整个区域外边界上的净流出。这是定理的几何图像,但把极限做严谨需要边界与向量场正则性。

定理不是说每一点的流量由散度单独决定,而是把局部源汇累积成全局边界账本。无散场可以在区域内部高速循环,只要进入与流出总体抵消,其边界净通量仍为零。

例子与计算

Ω=BR(0)R3F(x)=x=(x1,x2,x3)。此时 divF=3,故

ΩdivFdx=34πR33=4πR3.

在球面上 ν=x/R,所以 Fν=R,右侧为

BRRdS=R4πR2=4πR3.

两侧相等;半径增加时,虽然边界面积只按 R2 增长,边界处法向速度还多出一个 R 因子。

F(x,y)=(y,x)Ω 是平面上的单位圆盘,则 divF=0。圆周上 F 与切向平行、与外法向垂直,故边界通量也为零;这说明零散度允许旋转,并不要求 F=0

边界与失败情形

奇点不能被悄悄放进区域。三维场

F(x)=xx3

x0 处散度为零,但任意以原点为中心的球面通量等于 4π。它不满足在 Ω 邻域内 C1 的条件;分布意义下散度为 4πδ0。若只把经典散度零代入体积分,就会得到错误结论。

有孔区域的边界含多个分支。对环形区域,外边界取向外,内孔边界的“向外”却指向孔内;所有分支都沿区域的外法向定向,不能把每条曲线都按视觉上的逆时针统一处理。

无界区域需要截断后控制无穷远边界项。即使 divF 可积,若 F 的衰减不足,通量极限也未必存在。分片光滑或 Lipschitz 边界可以处理棱角,但高度不规则边界要转向有限周长集合等更一般理论。

有限体积法在网格控制体上写出的离散通量平衡受散度定理启发,却不是连续定理本身。数值通量、重构与时间离散会引入一致性和稳定性问题,不能用一句“由散度定理”替代离散误差分析。

推论与应用

对任意常数 c,取 F=cej 可恢复边界法向分量的整体平衡;取 F=uej 则得到分部积分公式

Ωjudx=ΩuνjdS.

更一般地对标量函数 u 和向量场 F 应用 div(uF)=uF+udivF,可推导 Green 恒等式及偏微分方程的弱形式。

流体质量守恒、电磁学中的 Gauss 定律、连续介质平衡律与有限体积离散都使用“体内源项等于边界净通量”这一接口。每个应用仍需说明向量场、区域边界与通量的具体正则性。

参考资料
  • Michael Spivak, Calculus on Manifolds, Westview Press, 1965,Ch. 4–5。
  • John M. Lee, Introduction to Smooth Manifolds, 2nd ed., Springer, 2013,Ch. 16。
  • Lawrence C. Evans, Partial Differential Equations, 2nd ed., AMS, 2010,Appendix C 与 Ch. 5。