三、矢量场的通量和散度 高斯定理
矢量场A通过一个截面S的通量ΦAΦA∬SA⋅dS∬SAcosθdS1式中θ为A与面元dS的法线en之间夹角,dSendS。如流速场中的流量,电场和磁场中的电场强度通量、磁通量,都属于“通量”的概念。令S为一闭合曲面,它包含的体积为ΔV,设想S面逐渐缩小到空间某点P。用ΦA代表矢量场A在闭合面SΦA∬SA⋅dS2当ΔV→0时,ΦA也趋于0。
1、定义
矢量场
A
\pmb{A}
A 通过一个截面
S
S
S 的通量
Φ
A
{\Phi_A}
ΦA 定义为下列面积分:
Φ
A
=
∬
S
A
⋅
d
S
=
∬
S
A
cos
θ
d
S
(1)
{\Phi_A} = \iint_S \pmb{A}\pmb{\cdot}{d}\pmb{S}=\iint_S {A}\cos\theta{d}S\tag{1}
ΦA=∬SA⋅dS=∬SAcosθdS(1)
式中
θ
\theta
θ 为
A
\pmb{A}
A 与面元
d
S
{d}\pmb{S}
dS 的法线
e
n
\pmb e_n
en 之间夹角,
d
S
=
e
n
d
S
d\pmb{S}=\pmb{e_n}d{S}
dS=endS。如流速场中的流量,电场和磁场中的电场强度通量、磁通量,都属于“通量”的概念。
令
S
S
S 为一闭合曲面,它包含的体积为
Δ
V
\Delta V
ΔV,设想
S
S
S 面逐渐缩小到空间某点
P
P
P。用
Φ
A
\Phi_A
ΦA 代表矢量场
A
\pmb{A}
A 在闭合面
S
S
S 上的通量:
Φ
A
=
∬
S
A
⋅
d
S
(2)
\Phi_A=\iint_S \pmb{A} \pmb{\cdot}d\pmb{S}\tag{2}
ΦA=∬SA⋅dS(2)
当
Δ
V
→
0
\Delta V\rightarrow0
ΔV→0 时,
Φ
A
\Phi_A
ΦA 也趋于0。若两者之比有一极限,则这个极限值为矢量场
A
\pmb A
A 在
P
P
P 点的散度,记作
d
i
v
\rm{div}
div
A
\pmb{A}
A 或
▽
⋅
A
\triangledown\pmb{\cdot A}
▽⋅A:
▽
⋅
A
≡
lim
Δ
V
→
0
Φ
A
Δ
V
=
lim
Δ
V
→
0
∬
S
A
⋅
d
S
Δ
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≡ΔV→0limΔVΦA=ΔV→0limΔV∬SA⋅dS(3)
矢量场的散度是个标量场。
2、散度的坐标表示式
上述散度的定义式(3)是与坐标的选取无关的,下面我们来研究它的直角坐标表示式。如图所示
以
P
P
P 点为中心取一个棱边分别与
x
,
y
,
z
x,y,z
x,y,z 轴平行的平行六面体,设边长分别为
Δ
x
,
Δ
y
,
Δ
z
\Delta x,\Delta y,\Delta z
Δx,Δy,Δz。现在来计算通过这平行六面体表面的通量。
先考虑与
x
x
x 轴垂直的一对表面。它们的面积为
Δ
y
Δ
z
\Delta y\Delta z
ΔyΔz。设
P
P
P 点的坐标为
(
x
,
y
,
z
)
(x,y,z)
(x,y,z),则这一对表面的
x
x
x 坐标分别为
x
−
Δ
/
2
x-\Delta/2
x−Δ/2 和
x
+
Δ
/
2
x+\Delta/2
x+Δ/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)。在计算通量的时候,只有与表面垂直的分量,即
A
x
A_x
Ax 分量起作用,它们在两表面上的数值分别是
A
x
(
x
−
Δ
x
/
2
,
y
,
z
)
\pmb{A_x}(x-\Delta x/2,y,z)
Ax(x−Δx/2,y,z) 和
A
x
(
x
+
Δ
x
/
2
,
y
,
z
)
\pmb{A_x}(x+\Delta x/2,y,z)
Ax(x+Δx/2,y,z),于是穿过这一对表面的通量分别是
A
x
(
x
−
Δ
x
/
2
,
y
,
z
)
Δ
y
Δ
z
\pmb{A_x}(x-\Delta x/2,y,z)\Delta y\Delta z
Ax(x−Δx/2,y,z)ΔyΔz 和
A
x
(
x
+
Δ
x
/
2
,
y
,
z
)
Δ
y
Δ
z
\pmb{A_x}(x+\Delta x/2,y,z)\Delta y\Delta z
Ax(x+Δx/2,y,z)ΔyΔz,二者一进一出,它们的代数和为
Φ
1
=
A
x
(
x
+
Δ
x
/
2
,
y
,
z
)
Δ
y
Δ
z
−
A
x
(
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Δz−Ax(x−Δx/2,y,z)ΔyΔz(4)
围绕
P
P
P 点将
A
x
A_x
Ax 按泰勒级数展开:
A
x
(
x
±
Δ
x
2
,
y
,
z
)
=
A
x
(
x
,
y
,
z
)
±
∂
A
x
∂
x
Δ
x
2
+
高次项
(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)±∂x∂Ax2Δx+高次项(5)
于是
Φ
1
=
[
A
x
(
x
,
y
,
z
)
+
∂
A
x
∂
x
Δ
x
2
]
Δ
y
Δ
z
−
[
A
x
(
x
,
y
,
z
)
−
∂
A
x
∂
x
Δ
x
2
]
Δ
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)+∂x∂Ax2Δx]ΔyΔz−[Ax(x,y,z)−∂x∂Ax2Δx]ΔyΔz+高次项(6)
即
Φ
1
=
∂
A
x
∂
x
Δ
x
Δ
y
Δ
z
+
高次项
(7)
\Phi_1=\frac{\partial A_x}{\partial x}\Delta x\Delta y\Delta z + 高次项\tag{7}
Φ1=∂x∂AxΔxΔyΔz+高次项(7)
同理可以得到穿过与
y
y
y 轴和
z
z
z 轴垂直的两对表面的通量代数和分别为
Φ
2
=
∂
A
y
∂
y
Δ
x
Δ
y
Δ
z
+
高次项
Φ
3
=
∂
A
z
∂
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=∂y∂AyΔxΔyΔz+高次项Φ3=∂z∂AzΔxΔyΔz+高次项(8)
最后我们得到穿过平行六面体六个表面的通量代数总和为
Φ
=
∯
平行六面体表面
A
⋅
d
S
=
Φ
1
+
Φ
2
+
Φ
3
=
(
∂
A
x
∂
x
+
∂
A
y
∂
y
+
∂
A
z
∂
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}
Φ=∬平行六面体表面A⋅dS=Φ1+Φ2+Φ3=(∂x∂Ax+∂y∂Ay+∂z∂Az)ΔxΔyΔz+高次项(9)
因为平行六面体的体积
Δ
=
Δ
x
Δ
y
Δ
z
\Delta = \Delta x\Delta y\Delta z
Δ=ΔxΔyΔz,按照散度的定义式(3)得
▽
⋅
A
=
lim
Δ
V
→
0
∬
S
A
⋅
d
S
Δ
V
=
lim
Δ
x
→
0
,
Δ
y
→
0
,
Δ
z
→
0
(
∂
A
x
∂
x
+
∂
A
y
∂
y
+
∂
A
z
∂
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=ΔV→0limΔV∬SA⋅dS=Δx→0,Δy→0,Δz→0limΔxΔyΔz(∂x∂Ax+∂y∂Ay+∂z∂Az)ΔxΔyΔz+高次项(10)
即
▽
⋅
A
=
∂
A
x
∂
x
+
∂
A
y
∂
y
+
∂
A
z
∂
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=∂x∂Ax+∂y∂Ay+∂z∂Az(11)
这就是散度的直角坐标表示式,下面我们不加推导地写出散度在其他常用坐标中的表达式,以备参考。
柱坐标
▽
⋅
A
=
1
ρ
∂
∂
ρ
(
ρ
A
ρ
)
+
1
ρ
∂
A
φ
∂
φ
+
∂
A
z
∂
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φ+∂z∂Az(12)
球坐标
▽
⋅
A
=
1
γ
2
[
∂
∂
γ
(
γ
2
A
γ
)
]
+
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) 中取任意闭合面
S
S
S,用
V
V
V 代表它所包围的体积。如图(a),用一曲面
D
D
D(下面叫它“隔板”)把体积
V
V
V 及其表面
S
S
S 分为两部分:
V
1
V_1
V1 和
V
2
V_2
V2,以及
S
1
′
S_1^\prime
S1′ 和
S
2
′
S_2^\prime
S2′,这里
V
1
+
V
2
=
V
,
S
1
′
+
S
2
′
=
S
V_1+V_2=V,S_1^\prime+S_2^\prime=S
V1+V2=V,S1′+S2′=S。体积
V
1
V_1
V1 的全部表面为
S
1
′
+
D
=
S
1
S_1^\prime +D=S_1
S1′+D=S1,体积
V
2
V_2
V2 的全部表面为
S
2
′
+
D
=
S
2
S_2^\prime+D=S_2
S2′+D=S2。穿过
S
1
S_1
S1 和
S
2
S_2
S2 的通量分别是
Φ
1
=
∯
S
1
A
⋅
d
S
1
=
∬
S
1
′
A
⋅
d
S
1
+
∬
D
A
⋅
d
S
1
Φ
2
=
∯
S
2
A
⋅
d
S
2
=
∬
S
2
′
A
⋅
d
S
2
+
∬
D
A
⋅
d
S
2
(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=∬S1A⋅dS1=∬S1′A⋅dS1+∬DA⋅dS1Φ2=∬S2A⋅dS2=∬S2′A⋅dS2+∬DA⋅dS2(14)
在上右式中右端的第二项虽然都是矢量场
A
\pmb A
A 穿过“隔板”
D
D
D 的通量,但是对于闭合高斯面
S
1
S_1
S1 和
S
2
S_2
S2 来说,在
D
D
D 上的外法线
n
1
\pmb {n_1}
n1 和
n
2
\pmb n_2
n2 方向相反,所以这两项绝对值相等,但符号正好相反。于是
Φ
1
+
Φ
2
=
∬
S
1
′
A
⋅
d
S
1
+
∬
S
2
′
A
⋅
d
S
2
=
∯
S
A
⋅
d
S
=
Φ
(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=∬S1′A⋅dS1+∬S2′A⋅dS2=∬SA⋅dS=Φ(15)
也就是说,将闭合曲面
S
S
S 所包围的空间用“隔板”隔开后,穿过两部分的通量的代数和不变,它仍等于穿过
S
S
S 的总通量
Φ
\Phi
Φ。
以上结论不难推广到把
V
V
V 分割成更多块的情况。这时我们有
Φ
=
∑
i
=
1
n
Φ
i
(16)
\Phi=\sum_{i=1}^n \Phi_i\tag{16}
Φ=i=1∑nΦi(16)
如果把体积
V
V
V 无限分割下去,使每块体积
Δ
V
i
\Delta V_i
ΔVi 都趋于 0,则按照散度的定义,
Φ
i
=
∯
S
i
A
⋅
d
S
i
=
(
▽
⋅
A
)
i
Δ
V
i
(17)
\Phi_i=\oiint_{S_i}\pmb{A\cdot}d\pmb{S_i}=(\triangledown\pmb{\cdot A})_i\Delta V_i\tag{17}
Φi=∬SiA⋅dSi=(▽⋅A)iΔVi(17)
其中
(
Δ
⋅
A
)
i
(\Delta\pmb{\cdot A})_i
(Δ⋅A)i 是
A
A
A 的散度在体积元
Δ
V
i
\Delta V_i
ΔVi 内的数值。把上式代入(16):
Φ
=
∯
S
A
⋅
d
S
≈
∑
i
=
1
n
(
▽
⋅
A
)
i
Δ
V
i
(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}
Φ=∬SA⋅dS≈i=1∑n(▽⋅A)iΔVi(18)
取极限后右端变为体积分:
∯
S
i
A
⋅
d
S
=
∭
V
▽
⋅
A
d
V
(19)
\oiint_{S_i}\pmb{A\cdot}d\pmb{S}=\iiint_V\triangledown\pmb{\cdot A}dV\tag{19}
∬SiA⋅dS=∭V▽⋅AdV(19)
式(19)表明:矢量场通过任意闭合曲面
S
S
S 的通量等于它所包围的体积
V
V
V 内散度的积分。这就是矢量场论中的高斯定理。
更多推荐
所有评论(0)