1、定义

矢量场 A\pmb{A}A 通过一个截面 SSS 的通量 ΦA{\Phi_A}ΦA 定义为下列面积分:
ΦA=∬SA⋅dS=∬SAcos⁡θdS(1) {\Phi_A} = \iint_S \pmb{A}\pmb{\cdot}{d}\pmb{S}=\iint_S {A}\cos\theta{d}S\tag{1} ΦA=SAdS=SAcosθdS(1)
式中 θ\thetaθA\pmb{A}A 与面元 dS{d}\pmb{S}dS 的法线 en\pmb e_nen 之间夹角,dS=endSd\pmb{S}=\pmb{e_n}d{S}dS=endS。如流速场中的流量,电场和磁场中的电场强度通量、磁通量,都属于“通量”的概念。
SSS 为一闭合曲面,它包含的体积为 ΔV\Delta VΔV,设想 SSS 面逐渐缩小到空间某点 PPP。用 ΦA\Phi_AΦA 代表矢量场 A\pmb{A}A 在闭合面 SSS 上的通量:
ΦA=∬SA⋅dS(2) \Phi_A=\iint_S \pmb{A} \pmb{\cdot}d\pmb{S}\tag{2} ΦA=SAdS(2)
ΔV→0\Delta V\rightarrow0ΔV0 时,ΦA\Phi_AΦA 也趋于0。若两者之比有一极限,则这个极限值为矢量场 A\pmb AAPPP 点的散度,记作 div\rm{div}divA\pmb{A}A▽⋅A\triangledown\pmb{\cdot A}A
▽⋅A≡lim⁡ΔV→0ΦAΔV=lim⁡ΔV→0∬SA⋅dSΔV(3) \triangledown\pmb{\cdot A}\equiv\lim_{\Delta V\rightarrow 0}\frac{\Phi_A}{\Delta V}=\lim_{\Delta V\rightarrow 0}\frac{\iint_{S} \pmb{A\cdot}d\pmb{S}}{\Delta V}\tag{3} AΔV0limΔVΦA=ΔV0limΔVSAdS(3)
矢量场的散度是个标量场。

2、散度的坐标表示式

上述散度的定义式(3)是与坐标的选取无关的,下面我们来研究它的直角坐标表示式。如图所示
在这里插入图片描述

散度直角坐标表示式的推导

PPP 点为中心取一个棱边分别与 x,y,zx,y,zx,y,z 轴平行的平行六面体,设边长分别为 Δx,Δy,Δz\Delta x,\Delta y,\Delta zΔx,Δy,Δz。现在来计算通过这平行六面体表面的通量。
先考虑与 xxx 轴垂直的一对表面。它们的面积为 ΔyΔz\Delta y\Delta zΔyΔz。设 PPP 点的坐标为 (x,y,z)(x,y,z)(x,y,z),则这一对表面的 xxx 坐标分别为 x−Δ/2x-\Delta/2xΔ/2x+Δ/2x+\Delta/2x+Δ/2,从而在这一对表面上的矢量场分别为 A(x−Δx/2,y,z)\pmb{A}(x-\Delta x/2,y,z)A(xΔx/2,y,z)A(x+Δx/2,y,z)\pmb{A}(x+\Delta x/2,y,z)A(x+Δx/2,y,z)。在计算通量的时候,只有与表面垂直的分量,即 AxA_xAx 分量起作用,它们在两表面上的数值分别是 Ax(x−Δx/2,y,z)\pmb{A_x}(x-\Delta x/2,y,z)Ax(xΔx/2,y,z)Ax(x+Δx/2,y,z)\pmb{A_x}(x+\Delta x/2,y,z)Ax(x+Δx/2,y,z),于是穿过这一对表面的通量分别是 Ax(x−Δx/2,y,z)ΔyΔz\pmb{A_x}(x-\Delta x/2,y,z)\Delta y\Delta zAx(xΔx/2,y,z)ΔyΔzAx(x+Δx/2,y,z)ΔyΔz\pmb{A_x}(x+\Delta x/2,y,z)\Delta y\Delta zAx(x+Δx/2,y,z)ΔyΔz,二者一进一出,它们的代数和为
Φ1=Ax(x+Δx/2,y,z)ΔyΔz−Ax(x−Δx/2,y,z)ΔyΔz(4) \Phi_1=\pmb{A_x}(x+\Delta x/2,y,z)\Delta y\Delta z - \pmb{A_x}(x-\Delta x/2,y,z)\Delta y\Delta z\tag{4} Φ1=Ax(x+Δx/2,y,z)ΔyΔzAx(xΔx/2,y,z)ΔyΔz(4)
围绕 PPP 点将 AxA_xAx 按泰勒级数展开:
Ax(x±Δx2,y,z)=Ax(x,y,z)±∂Ax∂xΔx2+高次项(5) \pmb{A_x}(x\pm\frac{\Delta x}{2},y,z)=\pmb{A_x}(x,y,z)\pm\frac{\partial A_x}{\partial x}\frac{\Delta x}{2}+高次项\tag{5} Ax(x±2Δx,y,z)=Ax(x,y,z)±xAx2Δx+高次项(5)
于是
Φ1=[Ax(x,y,z)+∂Ax∂xΔx2]ΔyΔz−[Ax(x,y,z)−∂Ax∂xΔx2]ΔyΔz+高次项(6) \Phi_1=\big[A_x(x,y,z)+\frac{\partial A_x}{\partial x}\frac{\Delta x}{2}\big]\Delta y\Delta z-\big[A_x(x,y,z)-\frac{\partial A_x}{\partial x}\frac{\Delta x}{2}\big]\Delta y\Delta z+高次项\tag{6} Φ1=[Ax(x,y,z)+xAx2Δx]ΔyΔz[Ax(x,y,z)xAx2Δx]ΔyΔz+高次项(6)

Φ1=∂Ax∂xΔxΔyΔz+高次项(7) \Phi_1=\frac{\partial A_x}{\partial x}\Delta x\Delta y\Delta z + 高次项\tag{7} Φ1=xAxΔxΔyΔz+高次项(7)
同理可以得到穿过与 yyy 轴和 zzz 轴垂直的两对表面的通量代数和分别为
Φ2=∂Ay∂yΔxΔyΔz+高次项Φ3=∂Az∂zΔxΔyΔz+高次项(8) \Phi_2=\frac{\partial A_y}{\partial y}\Delta x\Delta y\Delta z + 高次项\\ \quad\\ \Phi_3=\frac{\partial A_z}{\partial z}\Delta x\Delta y\Delta z + 高次项\tag{8} Φ2=yAyΔxΔyΔz+高次项Φ3=zAzΔxΔyΔz+高次项(8)
最后我们得到穿过平行六面体六个表面的通量代数总和为
Φ=∯平行六面体表面A⋅dS=Φ1+Φ2+Φ3=(∂Ax∂x+∂Ay∂y+∂Az∂z)ΔxΔyΔz+高次项(9) \begin{aligned} \Phi &=\oiint_{平行六面体表面}\pmb{A\cdot}d\pmb{S}=\Phi_1+\Phi_2+\Phi_3\\ &=\big(\frac{\partial A_x}{\partial x}+\frac{\partial A_y}{\partial y}+\frac{\partial A_z}{\partial z}\big)\Delta x\Delta y\Delta z+高次项\tag{9} \end{aligned} Φ= 平行六面体表面AdS=Φ1+Φ2+Φ3=(xAx+yAy+zAz)ΔxΔyΔz+高次项(9)
因为平行六面体的体积 Δ=ΔxΔyΔz\Delta = \Delta x\Delta y\Delta zΔ=ΔxΔyΔz,按照散度的定义式(3)得
▽⋅A=lim⁡ΔV→0∬SA⋅dSΔV=lim⁡Δx→0,Δy→0,Δz→0(∂Ax∂x+∂Ay∂y+∂Az∂z)ΔxΔyΔz+高次项ΔxΔyΔz(10) \begin{aligned} \triangledown\pmb{\cdot A}&=\lim_{\Delta V\rightarrow 0}\frac{\iint_{S} \pmb{A\cdot}d\pmb{S}}{\Delta V}\\ &=\lim_{\Delta x\rightarrow 0,\Delta y\rightarrow 0,\Delta z\rightarrow 0}\frac{\big(\frac{\partial A_x}{\partial x}+\frac{\partial A_y}{\partial y}+\frac{\partial A_z}{\partial z}\big)\Delta x\Delta y\Delta z+高次项}{\Delta x\Delta y\Delta z} \end{aligned}\tag{10} A=ΔV0limΔVSAdS=Δx0,Δy0,Δz0limΔxΔyΔz(xAx+yAy+zAz)ΔxΔyΔz+高次项(10)

▽⋅A=∂Ax∂x+∂Ay∂y+∂Az∂z(11) \triangledown\pmb{\cdot A}=\frac{\partial A_x}{\partial x}+\frac{\partial A_y}{\partial y}+\frac{\partial A_z}{\partial z}\tag{11} A=xAx+yAy+zAz(11)
这就是散度的直角坐标表示式,下面我们不加推导地写出散度在其他常用坐标中的表达式,以备参考。
柱坐标
▽⋅A=1ρ∂∂ρ(ρAρ)+1ρ∂Aφ∂φ+∂Az∂z(12) \triangledown\pmb{\cdot A}=\frac{1}{\rho}\frac{\partial}{\partial \rho}(\rho A_\rho)+\frac{1}{\rho}\frac{\partial A_\varphi}{\partial \varphi}+\frac{\partial A_z}{\partial z}\tag{12} A=ρ1ρ(ρAρ)+ρ1φAφ+zAz(12)
球坐标
▽⋅A=1γ2[∂∂γ(γ2Aγ)]+1γsin⁡θ[∂∂θ(sin⁡θAθ)]+1γsin⁡θ∂Aφ∂φ(13) \triangledown\pmb{\cdot A}=\frac{1}{\gamma^2}\big[\frac{\partial}{\partial\gamma}(\gamma^2A_\gamma)\big]+\frac{1}{\gamma\sin\theta}\big[\frac{\partial}{\partial\theta}(\sin\theta A_\theta)\big]+\frac{1}{\gamma\sin\theta}\frac{\partial A_\varphi}{\partial\varphi}\tag{13} A=γ21[γ(γ2Aγ)]+γsinθ1[θ(sinθAθ)]+γsinθ1φAφ(13)

3、高斯定理

在这里插入图片描述

高斯定理的证明

在矢量场 A(x,y,z)\pmb{A}(x,y,z)A(x,y,z) 中取任意闭合面 SSS,用 VVV 代表它所包围的体积。如图(a),用一曲面 DDD(下面叫它“隔板”)把体积 VVV 及其表面 SSS 分为两部分:V1V_1V1V2V_2V2,以及 S1′S_1^\primeS1S2′S_2^\primeS2,这里 V1+V2=V,S1′+S2′=SV_1+V_2=V,S_1^\prime+S_2^\prime=SV1+V2=V,S1+S2=S。体积 V1V_1V1 的全部表面为 S1′+D=S1S_1^\prime +D=S_1S1+D=S1,体积 V2V_2V2 的全部表面为 S2′+D=S2S_2^\prime+D=S_2S2+D=S2。穿过 S1S_1S1S2S_2S2 的通量分别是
Φ1=∯S1A⋅dS1=∬S1′A⋅dS1+∬DA⋅dS1Φ2=∯S2A⋅dS2=∬S2′A⋅dS2+∬DA⋅dS2(14) \Phi_1=\oiint_{S_1}\pmb{A\cdot}d\pmb{S_1}=\iint_{S_1^\prime}\pmb{A\cdot}d\pmb{S_1}+\iint_{D}\pmb{A\cdot}d\pmb{S_1}\\ \Phi_2=\oiint_{S_2}\pmb{A\cdot}d\pmb{S_2}=\iint_{S_2^\prime}\pmb{A\cdot}d\pmb{S_2}+\iint_{D}\pmb{A\cdot}d\pmb{S_2}\tag{14} Φ1= S1AdS1=S1AdS1+DAdS1Φ2= S2AdS2=S2AdS2+DAdS2(14)
在上右式中右端的第二项虽然都是矢量场 A\pmb AA 穿过“隔板” DDD 的通量,但是对于闭合高斯面 S1S_1S1S2S_2S2 来说,在 DDD 上的外法线 n1\pmb {n_1}n1n2\pmb n_2n2 方向相反,所以这两项绝对值相等,但符号正好相反。于是
Φ1+Φ2=∬S1′A⋅dS1+∬S2′A⋅dS2=∯SA⋅dS=Φ(15) \begin{aligned} \Phi_1+\Phi_2&=\iint_{S_1^\prime}\pmb{A\cdot}d\pmb{S_1}+\iint_{S_2^\prime}\pmb{A\cdot}d\pmb{S_2}\\ &=\oiint_S\pmb{A\cdot}d\pmb{S}=\Phi \end{aligned}\tag{15} Φ1+Φ2=S1AdS1+S2AdS2= SAdS=Φ(15)
也就是说,将闭合曲面 SSS 所包围的空间用“隔板”隔开后,穿过两部分的通量的代数和不变,它仍等于穿过 SSS 的总通量 Φ\PhiΦ
以上结论不难推广到把 VVV 分割成更多块的情况。这时我们有
Φ=∑i=1nΦi(16) \Phi=\sum_{i=1}^n \Phi_i\tag{16} Φ=i=1nΦi(16)
如果把体积 VVV 无限分割下去,使每块体积 ΔVi\Delta V_iΔVi 都趋于 0,则按照散度的定义,
Φi=∯SiA⋅dSi=(▽⋅A)iΔVi(17) \Phi_i=\oiint_{S_i}\pmb{A\cdot}d\pmb{S_i}=(\triangledown\pmb{\cdot A})_i\Delta V_i\tag{17} Φi= SiAdSi=(A)iΔVi(17)
其中 (Δ⋅A)i(\Delta\pmb{\cdot A})_i(ΔA)iAAA 的散度在体积元 ΔVi\Delta V_iΔVi 内的数值。把上式代入(16):
Φ=∯SA⋅dS≈∑i=1n(▽⋅A)iΔVi(18) \Phi=\oiint_{S}\pmb{A\cdot}d\pmb{S}\approx\sum_{i=1}^n(\triangledown\pmb{\cdot A})_i\Delta V_i\tag{18} Φ= SAdSi=1n(A)iΔVi(18)
取极限后右端变为体积分:
∯SiA⋅dS=∭V▽⋅AdV(19) \oiint_{S_i}\pmb{A\cdot}d\pmb{S}=\iiint_V\triangledown\pmb{\cdot A}dV\tag{19} SiAdS=VAdV(19)
式(19)表明:矢量场通过任意闭合曲面 SSS 的通量等于它所包围的体积 VVV 内散度的积分。这就是矢量场论中的高斯定理。

Logo

腾讯云面向开发者汇聚海量精品云计算使用和开发经验,营造开放的云计算技术生态圈。

更多推荐