高斯公式(散度定理)

高斯公式的意义

高斯公式,又称散度定理,建立了空间闭区域上的三重积分沿区域边界闭曲面的第二类曲面积分之间的联系。它是格林公式在三维空间的推广,将”区域内部的散度累积”与”边界上的通量”联系起来,是向量分析的核心定理之一。

高斯公式

定理1

设空间闭区域 Ω\Omega 由分片光滑的闭曲面 Σ\Sigma 围成,函数 P(x,y,z),Q(x,y,z),R(x,y,z)P(x,y,z), Q(x,y,z), R(x,y,z)Ω\Omega 上具有一阶连续偏导数,则有

ΣPdydz+Qdzdx+Rdxdy=Ω(Px+Qy+Rz)dv\oiint_\Sigma P \, dy\,dz + Q \, dz\,dx + R \, dx\,dy = \iiint_\Omega \left( \frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} + \frac{\partial R}{\partial z} \right) dv

其中 Σ\SigmaΩ\Omega 的整个边界曲面的外侧

几何解释
推论
证明
符号说明

用散度表示

引入散度的概念:F=Px+Qy+Rz\nabla \cdot \vec{F} = \frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} + \frac{\partial R}{\partial z},则高斯公式可以简洁地表示为:

ΣFdS=Ω(F)dv\oiint_\Sigma \vec{F} \cdot d\vec{S} = \iiint_\Omega (\nabla \cdot \vec{F}) \, dv

即:向量场通过闭曲面外侧的通量 = 散度在闭区域上的三重积分

高斯公式的物理意义

高斯公式的应用

应用一:简化闭曲面积分

当直接计算闭曲面上的第二类曲面积分困难时,可以利用高斯公式将其转化为三重积分。特别是当散度 F\nabla \cdot \vec{F} 很简单(如常数、零)时,三重积分往往更容易计算。

应用二:非闭曲面的补面法

如果曲面 Σ\Sigma 不是闭曲面,可以补上一块曲面 Σ1\Sigma_1,使 Σ+Σ1\Sigma + \Sigma_1 成为闭曲面,然后利用高斯公式:

ΣFdS=Σ+Σ1FdSΣ1FdS=Ω(F)dvΣ1FdS\iint_\Sigma \vec{F} \cdot d\vec{S} = \oiint_{\Sigma+\Sigma_1} \vec{F} \cdot d\vec{S} - \iint_{\Sigma_1} \vec{F} \cdot d\vec{S} = \iiint_\Omega (\nabla \cdot \vec{F}) dv - \iint_{\Sigma_1} \vec{F} \cdot d\vec{S}

其中 Σ1\Sigma_1 的侧需要选择使得 Σ+Σ1\Sigma + \Sigma_1 取外侧(或内侧,需一致)。

应用三:散度为零的场

如果向量场 F\vec{F} 的散度 F=0\nabla \cdot \vec{F} = 0(无源场/管状场),则通过任意闭曲面的通量为零:

ΣFdS=0\oiint_\Sigma \vec{F} \cdot d\vec{S} = 0

这意味着对于无源场,通过具有相同边界曲线的任意两块曲面的通量相等(与曲面形状无关,只依赖边界)。

典型例题

例题1:闭曲面上的通量

计算 Σxdydz+ydzdx+zdxdy\oiint_\Sigma x \, dy\,dz + y \, dz\,dx + z \, dx\,dy,其中 Σ\Sigma 是球面 x2+y2+z2=a2x^2+y^2+z^2 = a^2 的外侧。

参考答案(5 个标签)
高斯公式球面外侧通量散度

解题思路: 闭曲面外侧,用高斯公式转化为三重积分。散度为常数3。

详细步骤

  1. P=x,Q=y,R=zP = x, \, Q = y, \, R = z
  2. Px=1,Qy=1,Rz=1\frac{\partial P}{\partial x} = 1, \, \frac{\partial Q}{\partial y} = 1, \, \frac{\partial R}{\partial z} = 1
  3. 散度 F=3\nabla \cdot \vec{F} = 3
  4. 由高斯公式:ΣFdS=Ω3dv=3×球体积=3×43πa3=4πa3\oiint_\Sigma \vec{F} \cdot d\vec{S} = \iiint_\Omega 3 \, dv = 3 \times \text{球体积} = 3 \times \frac{4}{3}\pi a^3 = 4\pi a^3

答案4πa34\pi a^3

例题2:补面法

计算 Σ(x2cosα+y2cosβ+z2cosγ)dS\iint_\Sigma (x^2 \cos\alpha + y^2 \cos\beta + z^2 \cos\gamma) \, dS,其中 Σ\Sigma 是锥面 x2+y2=z2x^2+y^2 = z^2 介于 z=0z=0z=hz=hh>0h>0)之间部分的外侧(取下侧)。

参考答案(5 个标签)
高斯公式补面法锥面第一类曲面积分外侧

解题思路: 这是第一类曲面积分形式,但被积函数是 Fn\vec{F} \cdot \vec{n},等价于第二类曲面积分。锥面不是闭曲面,补面 z=hz=h 的上侧。

详细步骤

  1. 补面 Σ1:z=h,x2+y2h2\Sigma_1: z = h, \, x^2+y^2 \leq h^2,取上侧。则 Σ+Σ1\Sigma + \Sigma_1 构成闭曲面,取外侧。
  2. F=(x2,y2,z2)\vec{F} = (x^2, y^2, z^2),散度 F=2x+2y+2z\nabla \cdot \vec{F} = 2x+2y+2z
  3. 由高斯公式:Σ+Σ1FdS=Ω2(x+y+z)dv\oiint_{\Sigma+\Sigma_1} \vec{F} \cdot d\vec{S} = \iiint_\Omega 2(x+y+z) dv
  4. 由对称性,Ωxdv=0\iiint_\Omega x \, dv = 0Ωydv=0\iiint_\Omega y \, dv = 0
  5. 柱坐标计算 Ωzdv\iiint_\Omega z \, dvΩ:0θ2π,0ρh,ρzh\Omega: 0 \leq \theta \leq 2\pi, \, 0 \leq \rho \leq h, \, \rho \leq z \leq h Ωzdv=02πdθ0hρdρρhzdz=2π0hρ12(h2ρ2)dρ=π0h(h2ρρ3)dρ=π(h42h44)=πh44\iiint_\Omega z \, dv = \int_0^{2\pi} d\theta \int_0^h \rho \, d\rho \int_\rho^h z \, dz = 2\pi \int_0^h \rho \cdot \frac{1}{2}(h^2-\rho^2) d\rho = \pi \int_0^h (h^2\rho - \rho^3) d\rho = \pi(\frac{h^4}{2} - \frac{h^4}{4}) = \frac{\pi h^4}{4}
  6. Σ+Σ1=2×πh44=πh42\oiint_{\Sigma+\Sigma_1} = 2 \times \frac{\pi h^4}{4} = \frac{\pi h^4}{2}
  7. 计算补面 Σ1\Sigma_1 上的积分:z=hz=h,上侧,FdS=z2dxdy=h2dxdy\vec{F} \cdot d\vec{S} = z^2 dx\,dy = h^2 dx\,dy Σ1FdS=x2+y2h2h2dxdy=h2πh2=πh4\iint_{\Sigma_1} \vec{F} \cdot d\vec{S} = \iint_{x^2+y^2 \leq h^2} h^2 dx\,dy = h^2 \cdot \pi h^2 = \pi h^4
  8. 因此 Σ=Σ+Σ1Σ1=πh42πh4=πh42\iint_\Sigma = \oiint_{\Sigma+\Sigma_1} - \iint_{\Sigma_1} = \frac{\pi h^4}{2} - \pi h^4 = -\frac{\pi h^4}{2}

答案πh42-\frac{\pi h^4}{2}

例题3:散度为零的场

证明向量场 F=rr3\vec{F} = \frac{\vec{r}}{r^3}(其中 r=(x,y,z),r=x2+y2+z2\vec{r} = (x,y,z), \, r = \sqrt{x^2+y^2+z^2})在除原点外的区域内散度为零,并计算通过包围原点的任意闭曲面外侧的通量。

参考答案(5 个标签)
高斯公式散度为零点源电场奇点

解题思路: 计算散度,注意原点是奇点。对包围原点的闭曲面,用小球面挖去原点后应用高斯公式。

详细步骤

  1. F=(xr3,yr3,zr3)\vec{F} = (\frac{x}{r^3}, \frac{y}{r^3}, \frac{z}{r^3})
  2. x(xr3)=1r33x2r5\frac{\partial}{\partial x}(\frac{x}{r^3}) = \frac{1}{r^3} - \frac{3x^2}{r^5}
  3. 同理 y(yr3)=1r33y2r5\frac{\partial}{\partial y}(\frac{y}{r^3}) = \frac{1}{r^3} - \frac{3y^2}{r^5}z(zr3)=1r33z2r5\frac{\partial}{\partial z}(\frac{z}{r^3}) = \frac{1}{r^3} - \frac{3z^2}{r^5}
  4. 散度 =3r33(x2+y2+z2)r5=3r33r2r5=0= \frac{3}{r^3} - \frac{3(x^2+y^2+z^2)}{r^5} = \frac{3}{r^3} - \frac{3r^2}{r^5} = 0r0r \neq 0
  5. 对包围原点的闭曲面 Σ\Sigma,作小球面 Σε:x2+y2+z2=ε2\Sigma_\varepsilon: x^2+y^2+z^2 = \varepsilon^2(取内侧,即指向原点),在 Σ\SigmaΣε\Sigma_\varepsilon 之间的区域应用高斯公式: Σ+ΣεFdS=Ω0dv=0\oiint_{\Sigma + \Sigma_\varepsilon^-} \vec{F} \cdot d\vec{S} = \iiint_\Omega 0 \, dv = 0
  6. 因此 ΣFdS=ΣεFdS=Σε+FdS\oiint_\Sigma \vec{F} \cdot d\vec{S} = -\oiint_{\Sigma_\varepsilon^-} \vec{F} \cdot d\vec{S} = \oiint_{\Sigma_\varepsilon^+} \vec{F} \cdot d\vec{S}Σε+\Sigma_\varepsilon^+ 为外侧)
  7. 在小球面上 r=εr = \varepsilonF=rε3\vec{F} = \frac{\vec{r}}{\varepsilon^3},法向量 n=rε\vec{n} = \frac{\vec{r}}{\varepsilon} Fn=rrε4=ε2ε4=1ε2\vec{F} \cdot \vec{n} = \frac{\vec{r} \cdot \vec{r}}{\varepsilon^4} = \frac{\varepsilon^2}{\varepsilon^4} = \frac{1}{\varepsilon^2}
  8. 通量 =Σε1ε2dS=1ε24πε2=4π= \iint_{\Sigma_\varepsilon} \frac{1}{\varepsilon^2} dS = \frac{1}{\varepsilon^2} \cdot 4\pi\varepsilon^2 = 4\pi

答案:通过包围原点的任意闭曲面外侧的通量为 4π4\pi(与曲面形状无关,这就是静电场中的高斯定律)。


练习题

练习1

计算 Σ(xy)dydz+(yz)dzdx+(zx)dxdy\oiint_\Sigma (x-y) \, dy\,dz + (y-z) \, dz\,dx + (z-x) \, dx\,dy,其中 Σ\Sigma 是立方体 0x1,0y1,0z10 \leq x \leq 1, \, 0 \leq y \leq 1, \, 0 \leq z \leq 1 表面的外侧。

参考答案(4 个标签)
高斯公式立方体外侧散度

解题思路: 闭曲面外侧,用高斯公式。散度为常数3。

详细步骤

  1. P=xy,Q=yz,R=zxP = x-y, \, Q = y-z, \, R = z-x
  2. Px=1,Qy=1,Rz=1\frac{\partial P}{\partial x} = 1, \, \frac{\partial Q}{\partial y} = 1, \, \frac{\partial R}{\partial z} = 1,散度 =3= 3
  3. 由高斯公式:=Ω3dv=3×1=3= \iiint_\Omega 3 \, dv = 3 \times 1 = 3

答案33

练习2

计算 Σxz2dydz+(x2yz3)dzdx+(2xy+y2z)dxdy\iint_\Sigma xz^2 \, dy\,dz + (x^2y - z^3) \, dz\,dx + (2xy + y^2z) \, dx\,dy,其中 Σ\Sigma 是上半球面 z=a2x2y2z = \sqrt{a^2-x^2-y^2} 的上侧。

参考答案(4 个标签)
高斯公式补面法上半球上侧

解题思路: 上半球面不是闭曲面,补面 z=0z=0 的下侧,用高斯公式。

详细步骤

  1. 补面 Σ1:z=0,x2+y2a2\Sigma_1: z=0, \, x^2+y^2 \leq a^2,取下侧。Σ+Σ1\Sigma+\Sigma_1 为闭曲面外侧。
  2. P=xz2,Q=x2yz3,R=2xy+y2zP = xz^2, \, Q = x^2y-z^3, \, R = 2xy+y^2z
  3. 散度 =z2+x2+y2= z^2 + x^2 + y^2
  4. 由高斯公式:Σ+Σ1=Ω(x2+y2+z2)dv=02πdθ0π/2sinϕdϕ0ar2r2dr=2π1a55=2πa55\oiint_{\Sigma+\Sigma_1} = \iiint_\Omega (x^2+y^2+z^2) dv = \int_0^{2\pi} d\theta \int_0^{\pi/2} \sin\phi \, d\phi \int_0^a r^2 \cdot r^2 dr = 2\pi \cdot 1 \cdot \frac{a^5}{5} = \frac{2\pi a^5}{5}
  5. 补面上 z=0z=0,下侧,FdS=Rdxdy=2xydxdy\vec{F} \cdot d\vec{S} = R \, dx\,dy = 2xy \, dx\,dy(下侧取负) Σ1=Dxy2xydxdy=0\iint_{\Sigma_1} = -\iint_{D_{xy}} 2xy \, dx\,dy = 0(由对称性)
  6. 因此 Σ=2πa550=2πa55\iint_\Sigma = \frac{2\pi a^5}{5} - 0 = \frac{2\pi a^5}{5}

答案2πa55\frac{2\pi a^5}{5}


总结

本文出现的符号

符号类型读音/说明在本文中的含义
Σ\oiint_\Sigma闭曲面积分closed surface integral沿闭曲面的积分
F\nabla \cdot \vec{F}散度divergence向量场的散度,Px+Qy+Rz\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z}
Ω\Omega空间区域spatial domain闭曲面围成的空间区域
Σ\Sigma边界曲面boundary surface区域 Ω 的边界闭曲面
r\vec{r}位置向量position vector从原点到点(x,y,z)的向量
rr距离distance点到原点的距离,r=rr=\|\vec{r}\|

中英对照

中文术语英文术语音标说明
高斯公式Gauss’s theorem/ɡaʊsɪz ˈθɪərəm/散度定理,三重积分与边界曲面积分的关系
散度定理divergence theorem/daɪˈvɜːrdʒəns ˈθɪərəm/即高斯公式
散度divergence/daɪˈvɜːrdʒəns/向量场的发散程度
通量flux/flʌks/向量场通过曲面的量
source/sɔːrs/散度为正的点,向量场向外发散
sink/sɪŋk/散度为负的点,向量场向内汇聚
补面法surface completion method/ˈsɜːrfɪs kəmˈpliːʃən ˈmɛθəd/补面使曲面闭合以应用高斯公式
无源场source-free field/sɔːrs friː fiːld/散度为零的向量场,又称管状场
管状场solenoidal field/səˈlɛnɔɪdəl fiːld/即无源场,散度为零
奇点singularity/ˌsɪŋɡjəˈlærəti/函数无定义或不光滑的点(如原点)