在这篇文章内,矢量 与标量 分别用粗体 与斜体 显示。例如,位置矢量通常用
r
{\displaystyle \mathbf{r}\,\!}
表示;而其大小则用
r
{\displaystyle r\,\!}
来表示。
两个同性电荷的电场线和等势线。
在静电学 里,电势 ,或作电位 (英语:electric potential )[1] ,是标量 ,定义为处于电场 中某个位置的单位电荷 所具有的电势能 。[2] 其数值不具有绝对意义,只具有相对意义,因此为了便于分析问题,必须设定一个参考位置,并把它设为零,称为零势能点。通常,会把无穷远处的电势设定为零。那么,电势可以定义如下:假设检验电荷从无穷远位置,经过任意路径,克服电场力,以缓慢、没有产生加速度的方式移动到某位置,则在这位置的电势,等于因移动检验电荷所做的功 与检验电荷的电荷量 的比值。在国际单位制 里,电势的单位为伏特 (
V
=
J
/
C
{\displaystyle \scriptstyle{\text{V} = \text{J}/\text{C}}}
)(Volt),它是为了纪念意大利物理学家亚历山德罗·伏打 (Alessandro Volta)而命名。
电势必需满足泊松方程 ,同时符合相关边界条件;假设在某区域内的电荷密度 为零,则泊松方程约化为拉普拉斯方程 ,电势必需满足拉普拉斯方程。
在电动力学 里,当含时电磁场 存在的时候,电势可以延伸为“广义电势”。特别注意,广义电势不能被视为电势能每单位电荷。
简介
处于外电场 的带电粒子 会受到外电场施加的作用力 ,称为电场力 ,促使带电粒子加速运动。对于带正电粒子,电场力与电场 同方向;对于带负电粒子,电场力与电场反方向。电场力的数值大小与电荷量、电场数值大小成正比。
作用力 与势能 之间有非常直接的关系。随着物体朝著作用力的方向的加速运动,物体的动能 变大,势能变小。例如,一个石头在山顶的重力势能 大于在山脚的重力势能。随着物体的滚落,重力势能变小,动能变大。
对于某种特别作用力,科学家可以定义其矢量场 和其位势 ,使得物体因为这矢量场而具有的势能,只与物体位置、参考位置之间的距离有关。称这种作用力为保守力 ,这种矢量场为保守场 。
例如,重力 、静电场的电场力,都是保守力。静电场的标势 称为电势 ,或称为静电势 。
电势和磁矢势 共同形成一个四维矢量 ,称为四维势 。从某一个惯性参考系 观察到的四维势,应用洛伦兹变换 ,可以计算出另外一个惯性参考系所观察到的四维势。
静电学里的电势
在静电学里,电场
E
{\displaystyle \mathbf{E}}
内某位置
r
{\displaystyle \mathbf{r}}
的电势
ϕ
{\displaystyle \phi}
,以方程定义为[2]
ϕ
(
r
)
=
d
e
f
U
E
(
r
)
/
q
{\displaystyle \phi(\mathbf{r})\ \stackrel{def}{=}\ U_\mathrm{E}(\mathbf{r})/q}
;
其中,
U
E
{\displaystyle U_\mathrm{E}}
是在位置
r
{\displaystyle \mathbf{r}}
的检验电荷
q
{\displaystyle q}
所具有的电势能。
电势能的数值是人为设定的,没有绝对意义,只有相对于某参考位置的已设定参考值时才有物理意义。假若要设定电势能在空间任意位置的数值,必须先设定其在某参考位置
r
0
{\displaystyle \mathbf{r}_0}
的数值。为了方便运算,假设其参考数值为0。然后,就可以将在位置
r
{\displaystyle \mathbf{r}}
的电势能
U
E
(
r
)
{\displaystyle U_\mathrm{E}(\mathbf{r})}
定义为从参考位置
r
0
{\displaystyle \mathbf{r}_0}
缓慢地将检验电荷
q
{\displaystyle q}
移动至
r
{\displaystyle \mathbf{r}}
所需做的机械功
W
{\displaystyle W}
:
U
E
(
r
)
=
d
e
f
W
{\displaystyle U_\mathrm{E}(\mathbf{r})\ \stackrel{def}{=}\ W}
。
移动检验电荷时所施加的外力
F
{\displaystyle \mathbf{F}}
,必须恰巧抵消处于电场
E
{\displaystyle \mathbf{E}}
的检验电荷
q
{\displaystyle q}
所感受到的电场力
q
E
{\displaystyle q\mathbf{E}}
,即
F
=
−
q
E
{\displaystyle \mathbf{F}=-q\mathbf{E}}
。其所做机械功等于外力
F
{\displaystyle \mathbf{F}}
的路径积分:
W
=
∫
L
F
⋅
d
ℓ
=
−
q
∫
L
E
⋅
d
ℓ
{\displaystyle W=\int_\mathbb{L} \mathbf{F} \cdot \mathrm{d} \boldsymbol{\ell}=-q\int_\mathbb{L} \mathbf{E} \cdot \mathrm{d} \boldsymbol{\ell}}
;
其中,
L
{\displaystyle \mathbb{L} }
是从参考位置
r
0
{\displaystyle \mathbf{r}_0}
到位置
r
{\displaystyle \mathbf{r}}
的一条任意路径,
d
ℓ
{\displaystyle \mathrm{d} \boldsymbol{\ell}}
是微小线元素。
在静电学里,
∇
×
E
=
0
{\displaystyle \mathbf{\nabla} \times \mathbf{E} = 0}
,电场是保守场,所以,在积分时,可以选择任意路径
L
{\displaystyle \mathbb{L} }
,计算出来的结果都一样。欲知更详尽细节,请参阅条目保守力 。由于这方程右边的路径积分 跟路径
L
{\displaystyle \mathbb{L} }
无关,只跟路径的初始位置
r
0
{\displaystyle \mathbf{r}_0}
、终止位置
r
{\displaystyle \mathbf{r}}
有关,因此若能够假设无穷远位置
∞
{\displaystyle \infty}
的电势能为0,则可以设定参考位置
r
0
{\displaystyle \mathbf{r}_0}
在无穷远位置
∞
{\displaystyle \infty}
:
U
E
(
r
)
=
−
q
∫
∞
r
E
⋅
d
ℓ
{\displaystyle U_\mathrm{E}(\mathbf{r})= - q\int_\infty^\mathbf{r} \mathbf{E} \cdot \mathrm{d} \boldsymbol{\ell}}
。
所以,电势就是从无穷远位置到检验位置对于电场做路径积分所得结果的负值:
ϕ
(
r
)
=
−
∫
∞
r
E
⋅
d
ℓ
{\displaystyle \phi(\mathbf{r})= - \int_\infty^\mathbf{r} \mathbf{E} \cdot \mathrm{d} \boldsymbol{\ell}}
。
在任意两个位置
r
1
{\displaystyle \mathbf{r}_1}
、
r
2
{\displaystyle \mathbf{r}_2}
之间的“电势差”
Δ
ϕ
{\displaystyle \Delta\phi}
为
Δ
ϕ
=
ϕ
(
r
2
)
−
ϕ
(
r
1
)
=
−
∫
r
1
r
2
E
⋅
d
ℓ
{\displaystyle \Delta\phi=\phi(\mathbf{r}_2)-\phi(\mathbf{r}_1)= - \int_{\mathbf{r}_1}^{\mathbf{r}_2} \mathbf{E} \cdot \mathrm{d} \boldsymbol{\ell}}
。
由于电场
E
{\displaystyle \mathbf{E}}
是保守场,电势差也与积分路径无关,只跟积分路径的初始位置与终止位置有关。
点电荷
由点电荷 Q 所产生的电势,在距离 r 时,可表示为
V
=
1
4
π
ε
0
Q
r
{\displaystyle V = \frac{1}{4 \pi \varepsilon_0} \frac{Q}{r} }
其中,ε 0 是真空电容率 。
在无限远处,电势为零。由多个点电荷产生的电势,相等于各点电荷所产生的电势之和。此外,电势场是标量场 ,电场则是矢量场 。
叠加原理
电场遵守叠加原理 :假设在三维空间里,由两组完全不相交的电荷分布所产生的电场分别为
E
1
{\displaystyle \mathbf{E}_1}
、
E
2
{\displaystyle \mathbf{E}_2}
,则总电场为
E
t
=
E
1
+
E
2
{\displaystyle \mathbf{E}_t=\mathbf{E}_1+\mathbf{E}_2}
。
总电势为每单位电荷 克服电场力所做的机械功之和:
ϕ
t
(
r
)
=
−
∫
∞
r
E
t
⋅
d
ℓ
=
−
∫
∞
r
(
E
1
+
E
2
)
⋅
d
ℓ
=
ϕ
1
(
r
)
+
ϕ
2
(
r
)
{\displaystyle \phi_t(\mathbf{r})= - \int_\infty^\mathbf{r} \mathbf{E}_t \cdot \mathrm{d} \boldsymbol{\ell}= - \int_\infty^\mathbf{r} (\mathbf{E}_1+\mathbf{E}_2) \cdot \mathrm{d} \boldsymbol{\ell}=\phi_1(\mathbf{r})+\phi_2(\mathbf{r})}
。
所以,电势也遵守叠加原理。当计算一组电荷分布所产生的电势时,只需要知道在电荷分布的每个源位置的单独电荷所产生在检验位置的电势,就可以应用积分运算,得到整个电荷分布所产生在检验位置的电势。
电势的微分方程
应用积分符号内取微分 方法,电势的梯度 为
∇
ϕ
(
r
)
=
−
∇
∫
∞
r
E
(
r
′
)
⋅
d
ℓ
′
=
−
E
(
r
)
{\displaystyle \mathbf{\nabla} \phi(\mathbf{r})= - \mathbf{\nabla} \int_\infty^\mathbf{r} \mathbf{E}(\mathbf{r}') \cdot \mathrm{d} \boldsymbol{\ell}^{\,\prime}=- \mathbf{E}(\mathbf{r})}
。
所以,电场与电势之间的关系为
E
(
r
)
=
−
∇
ϕ
(
r
)
{\displaystyle \mathbf{E}(\mathbf{r}) = - \mathbf{\nabla} \phi(\mathbf{r})}
。
根据高斯定律 的方程,
∇
⋅
E
=
ρ
/
ϵ
0
{\displaystyle \mathbf{\nabla} \cdot \mathbf{E}=\rho / \epsilon_0}
;
其中,
ρ
{\displaystyle \rho}
是电荷密度 ,
ϵ
0
{\displaystyle \epsilon_0}
是电常数 。
所以,电势满足泊松方程 :
∇
2
ϕ
=
−
ρ
/
ϵ
0
{\displaystyle \nabla^2 \phi= - \rho/ \epsilon_0}
。
假设电荷密度为零,则这方程变为拉普拉斯方程 :
∇
2
ϕ
=
0
{\displaystyle \nabla^2 \phi=0}
。
请注意,假若
∇
×
E
≠
0
{\displaystyle \mathbf{\nabla}\times\mathbf{E} \ne 0}
,也就是说,电场不具保守性 (由于随时间变化的磁场 造成的效应;参阅麦克斯韦方程组 ),则不能使用这些方程。
由于电势乃是标量,而电场是具有三个分量的矢量,所以,很多时候,使用电势来解析问题会省去很多运算工作,带来很大的便利。
拉普拉斯方程的解答
在某空间区域内,假设电荷密度 为零,则电势必须满足拉普拉斯方程 ,并且符合所有相关边界条件 。
边界条件
在静电学里,有三种边界条件:
狄利克雷边界条件 :在所有边界,电势都已良态给定。具有这种边界条件的问题称为狄利克雷问题 。
纽曼边界条件 :在所有边界,电势的法向导数都已良态给定。具有这种边界条件的问题称为纽曼问题 。
混合边界条件:一部分边界的电势都已良态给定,其它边界的电势的法向导数也已良态给定。
根据拉普拉斯方程的唯一性定理 ,对于这些种类的边界条件,拉普拉斯方程的解答都具有唯一性。所以,只要找到一个符合边界条件的解答,则这解答必定为正确解答。
分离变数法
应用分离变数法 来解析拉普拉斯方程,可以将问题的偏微分方程改变为一组较容易解析的常微分方程 。对于一般问题,通常会采用直角坐标系、圆柱坐标系或球坐标系来分离拉普拉斯方程。但是,对于其它比较特别的问题,另外还有八种坐标系可以用来分离拉普拉斯方程。[3] 分离之后,找到每一个常微分方程的通解(通常为一组本征方程的叠加),电势可以表达为这些通解的乘积。将这表达式与边界条件相匹配,就可以设定一般解的系数,从而找到问题的特解。根据拉普拉斯方程的唯一性定理,这特解也是唯一的正确解答。
两个半平面导体案例
被位于
y
=
0
{\displaystyle y=0}
的绝缘线条分隔为处于y+ 、y- -半平面的两个导体的电势分别设定为
+
V
{\displaystyle +V}
、
−
V
{\displaystyle -V}
。
假设在xy-平面的无限平面导体 被一条位于
y
=
0
{\displaystyle y=0}
的绝缘线条分为两半,两个处于y+ 、y- -半平面的导体的电势分别设定为
+
V
{\displaystyle +V}
、
−
V
{\displaystyle -V}
,则计算z+ -半空间任意位置的电势这问题,由于边界条件的几何形状适合用直角坐标来描述,可以以直角坐标
(
x
,
y
,
z
)
{\displaystyle (x,y,z)}
将拉普拉斯方程表示为:
∇
2
ϕ
=
∂
2
ϕ
∂
x
2
+
∂
2
ϕ
∂
y
2
+
∂
2
ϕ
∂
z
2
=
0
{\displaystyle \nabla^2 \phi=\frac{\partial^2 \phi}{\partial x^2 }+\frac{\partial ^2\phi}{\partial y^2 }+\frac{\partial^2\phi}{\partial z^2 }=0}
。
因为这案例与x-坐标无关,方程可以简化为
∇
2
ϕ
(
y
,
z
)
=
∂
2
ϕ
∂
y
2
+
∂
2
ϕ
∂
z
2
=
0
{\displaystyle \nabla^2 \phi(y,z)=\frac{\partial ^2\phi}{\partial y^2 }+\frac{\partial^2\phi}{\partial z^2 }=0}
。
应用分离变数法,猜想解答的形式为
ϕ
(
y
,
z
)
=
Y
(
y
)
Z
(
z
)
{\displaystyle \phi(y,z)=Y(y)Z(z)}
。
将这公式代入拉普拉斯方程,则可得到
1
Y
(
y
)
d
2
Y
(
y
)
d
y
2
+
1
Z
(
z
)
d
2
Z
(
z
)
d
z
2
=
0
{\displaystyle \frac{1}{Y(y)}\ \frac{\mathrm{d}^2 Y(y)}{\mathrm{d}y^2}+\frac{1}{Z(z)}\ \frac{\mathrm{d}^2 Z(z)}{\mathrm{d}z^2}=0}
。
注意到这方程的每一个项目都只含有一个变量,并且跟其它变量无关。所以,每一个项目都等于常数:
1
Y
(
y
)
d
2
Y
(
y
)
d
y
2
=
C
{\displaystyle \frac{1}{Y(y)}\ \frac{\mathrm{d}^2 Y(y)}{\mathrm{d}y^2}=C}
、
1
Z
(
z
)
d
2
Z
(
z
)
d
z
2
=
−
C
{\displaystyle \frac{1}{Z(z)}\ \frac{\mathrm{d}^2 Z(z)}{\mathrm{d}z^2}=-C}
。
这样,一个二次偏微分方程被改变为两个简单的二次常微分方程。解答分别为
Y
(
y
)
=
A
1
e
i
k
y
+
A
2
e
−
i
k
y
{\displaystyle Y(y)=A_1 e^{iky}+A_2 e^{-iky}}
、
Z
(
z
)
=
B
1
e
k
z
+
B
2
e
−
k
z
{\displaystyle Z(z)=B_1 e^{kz}+B_2 e^{-kz}}
;
其中,
A
1
(
k
)
{\displaystyle A_1(k)}
、
A
2
(
k
)
{\displaystyle A_2(k)}
、
B
1
(
k
)
{\displaystyle B_1(k)}
、
B
2
(
k
)
{\displaystyle B_2(k)}
都是系数函数。
当
z
{\displaystyle z}
趋向于无穷大时,
Z
(
z
)
{\displaystyle Z(z)}
趋向于零,所以,
B
1
=
0
{\displaystyle B_1=0}
。综合起来,电势为
ϕ
(
y
,
z
)
=
∫
0
∞
(
A
1
e
i
k
y
+
A
2
e
−
i
k
y
)
e
−
k
z
d
k
{\displaystyle \phi(y,z)=\int_0^{\infty}(A_1 e^{iky}+A_2 e^{-iky})e^{-kz}\mathrm{d}k}
。
由于在
z
=
0
{\displaystyle z=0}
,y+ 、y- -半平面的电势分别为
+
V
{\displaystyle +V}
、
−
V
{\displaystyle -V}
,所以,
当
y
>
0
{\displaystyle y>0}
时,
∫
0
∞
(
A
1
e
i
k
y
+
A
2
e
−
i
k
y
)
d
k
=
+
V
{\displaystyle \int_0^{\infty}(A_1 e^{iky}+A_2 e^{-iky})\mathrm{d}k=+V}
、
当
y
<
0
{\displaystyle y<0}
时,
∫
0
∞
(
A
1
e
i
k
y
+
A
2
e
−
i
k
y
)
d
k
=
−
V
{\displaystyle \int_0^{\infty}(A_1 e^{iky}+A_2 e^{-iky})\mathrm{d}k=-V}
。
应用傅里叶变换,可以得到
A
1
(
k
)
=
V
2
π
(
∫
0
∞
e
−
i
k
y
′
d
y
′
−
∫
−
∞
0
e
−
i
k
y
′
d
y
′
)
{\displaystyle A_1(k)=\frac{V}{2\pi}\left( \int_0^{\infty} e^{-iky'}\mathrm{d}y'-\int_{-\infty}^0 e^{-iky'}\mathrm{d}y' \right)}
、
A
2
(
k
)
=
V
2
π
(
∫
0
∞
e
i
k
y
′
d
y
′
−
∫
−
∞
0
e
i
k
y
′
d
y
′
)
{\displaystyle A_2(k)=\frac{V}{2\pi}\left( \int_0^{\infty} e^{iky'}\mathrm{d}y'-\int_{-\infty}^0 e^{iky'}\mathrm{d}y' \right)}
。
所以,由
A
1
(
k
)
{\displaystyle A_1(k)}
项目贡献出的电势为
ϕ
1
=
V
2
π
∫
0
∞
d
k
{
∫
0
∞
e
i
k
(
y
−
y
′
)
−
k
z
d
y
′
−
∫
−
∞
0
e
i
k
(
y
−
y
′
)
−
k
z
d
y
′
}
=
−
V
2
π
∫
0
∞
d
y
′
i
(
y
−
y
′
)
−
z
+
V
2
π
∫
−
∞
0
d
y
′
i
(
y
−
y
′
)
−
z
{\displaystyle \begin{align}\phi_1 & =\frac{V}{2\pi}\int_0^{\infty}\mathrm{d}k\left\{\int_0^{\infty}e^{ik(y-y')-kz}\mathrm{d}y' - \int_{-\infty}^0e^{ik(y-y')-kz}\mathrm{d}y' \right\} \\
& =-\ \frac{V}{2\pi}\int_0^{\infty}\frac{\mathrm{d}y'}{i(y-y')-z}+\ \frac{V}{2\pi}\int_{-\infty}^0\frac{\mathrm{d}y'}{i(y-y')-z} \\
\end{align}}
。
类似地,由
A
2
(
k
)
{\displaystyle A_2(k)}
项目贡献出的电势为
ϕ
2
=
V
2
π
∫
0
∞
d
k
{
∫
0
∞
e
−
i
k
(
y
−
y
′
)
−
k
z
d
y
′
−
∫
−
∞
0
e
−
i
k
(
y
−
y
′
)
−
k
z
d
y
′
}
=
−
V
2
π
∫
0
∞
d
y
′
−
i
(
y
−
y
′
)
−
z
+
V
2
π
∫
−
∞
0
d
y
′
−
i
(
y
−
y
′
)
−
z
{\displaystyle \begin{align}\phi_2 & =\frac{V}{2\pi}\int_0^{\infty}\mathrm{d}k\left\{\int_0^{\infty}e^{-ik(y-y')-kz}\mathrm{d}y' - \int_{-\infty}^0e^{-ik(y-y')-kz}\mathrm{d}y' \right\} \\
& =-\ \frac{V}{2\pi}\int_0^{\infty}\frac{\mathrm{d}y'}{-i(y-y')-z}+\ \frac{V}{2\pi}\int_{-\infty}^0\frac{\mathrm{d}y'}{-i(y-y')-z} \\
\end{align}}
。
总电势为[4]
ϕ
=
V
z
π
∫
0
∞
d
y
′
(
y
−
y
′
)
2
+
z
2
−
V
z
π
∫
−
∞
0
d
y
′
(
y
−
y
′
)
2
+
z
2
=
2
V
π
arctan
(
y
z
)
{\displaystyle \begin{align}\phi & =\frac{Vz}{\pi}\int_0^{\infty}\frac{\mathrm{d}y'}{(y-y')^2+z^2}-\ \frac{Vz}{\pi}\int_{-\infty}^0\frac{\mathrm{d}y'}{(y-y')^2+z^2} \\
& =\frac{2V}{\pi}\ \arctan{\left(\frac{y}{z}\right)} \\
\end{align}}
。
泊松方程的解答
电荷分布所产生的电势
根据库仑定律 ,一个源位置为
r
′
{\displaystyle \mathbf{r}'}
的点电荷
q
{\displaystyle q}
,所产生在任意位置
r
{\displaystyle \mathbf{r}}
的电场为
E
(
r
)
=
q
4
π
ϵ
0
(
r
−
r
′
)
|
r
−
r
′
|
3
{\displaystyle \mathbf{E}(\mathbf{r}) = \frac{q}{4 \pi \epsilon_0}\ \frac{(\mathbf{r} - \mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|^3}}
。
对于一群点电荷,应用叠加原理 ,总电场等于每一个点电荷所产生的电场的叠加。体积区域
V
′
{\displaystyle \mathbb{V}'}
内部电荷密度为
ρ
(
r
′
)
{\displaystyle \rho(\mathbf{r}')}
的电荷分布,在检验位置
r
{\displaystyle \mathbf{r}}
所产生的电场为
E
(
r
)
=
1
4
π
ϵ
0
∫
V
′
ρ
(
r
′
)
(
r
−
r
′
)
|
r
−
r
′
|
3
d
3
r
′
{\displaystyle \mathbf{E}(\mathbf{r}) =\frac{1}{4 \pi \epsilon_0}\int_{\mathbb{V}'} \rho(\mathbf{r}')\frac{(\mathbf{r} - \mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|^3}\ \mathrm{d}^3 r'}
;
其中,
d
3
r
′
{\displaystyle \mathrm{d}^3 r'}
是微小体积元素。
应用一条矢量恒等式 ,
∇
1
|
r
−
r
′
|
=
−
(
r
−
r
′
)
|
r
−
r
′
|
3
{\displaystyle \nabla\frac{1}{|\mathbf{r} - \mathbf{r}'|}=-\ \frac{(\mathbf{r} - \mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|^3}}
,
可以得到
E
(
r
)
=
−
1
4
π
ϵ
0
∇
∫
V
′
ρ
(
r
′
)
|
r
−
r
′
|
d
3
r
′
{\displaystyle \mathbf{E}(\mathbf{r}) =-\ \frac{1}{4 \pi \epsilon_0}\nabla\int_{\mathbb{V}'} \frac{\rho(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|}\ \mathrm{d}^3 r'}
。
设定在无穷远的电势为参考值0,则在任意位置的电势为
ϕ
(
r
)
=
1
4
π
ϵ
0
∫
V
′
ρ
(
r
′
)
|
r
−
r
′
|
d
3
r
′
{\displaystyle \phi(\mathbf{r}) = \frac{1}{4 \pi \epsilon_0}\int_{\mathbb{V}'} \frac{\rho(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|}\ \mathrm{d}^3 r'}
;(1)
应用一则关于狄拉克δ函数 的矢量恒等式
∇
2
(
1
|
r
−
r
′
|
)
=
−
4
π
δ
(
r
−
r
′
)
{\displaystyle \nabla^2 \left(\frac{1}{|\mathbf{r} - \mathbf{r}'|}\right)
= - 4\pi\delta(\mathbf{r} - \mathbf{r}')}
,
假设检验位置
r
{\displaystyle \mathbf{r}}
在积分体积
V
′
{\displaystyle \mathbb{V}'}
内,则可得到泊松方程:
∇
2
ϕ
(
r
)
=
1
4
π
ϵ
0
∫
V
′
∇
2
(
ρ
(
r
′
)
|
r
−
r
′
|
)
d
3
r
′
=
−
1
ϵ
0
∫
V
′
ρ
(
r
′
)
δ
(
r
−
r
′
)
d
3
r
′
=
−
ρ
(
r
)
ϵ
0
{\displaystyle \nabla^2\phi(\mathbf{r}) = \frac{1}{4 \pi \epsilon_0}\int_{\mathbb{V}'} \nabla^2 \left(\frac{\rho(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|}\right)
\ \mathrm{d}^3 r'=-\ \frac{1}{\epsilon_0}\int_{\mathbb{V}'} \rho(\mathbf{r}')\delta(\mathbf{r} - \mathbf{r}')\ \mathrm{d}^3 r'
=-\ \frac{\rho(\mathbf{r})}{\epsilon_0}}
。
所以,电势的方程(1)为泊松方程的解答。
边界条件
电势的方程(1)只考虑到一群电荷分布所产生的电势。假若遭遇边界条件为电势的静电学问题,就不能使用方程(1),必需使用更具功能的方法。
根据格林第二恒等式 ,对于任意良态函数
ϕ
(
r
)
{\displaystyle \phi(\mathbf{r})}
与
ψ
(
r
)
{\displaystyle \psi(\mathbf{r})}
,[5]
∫
V
(
ϕ
∇
2
ψ
−
ψ
∇
2
ϕ
)
d
3
r
=
∮
S
(
ϕ
∂
ψ
∂
n
−
ψ
∂
ϕ
∂
n
)
d
2
r
{\displaystyle \int_{\mathbb{V}} \left( \phi \nabla^2 \psi - \psi \nabla^2 \phi\right)\ \mathrm{d}^3 r= \oint_{\mathbb{S}} \left( \phi {\partial \psi \over \partial n} - \psi {\partial \phi \over \partial n}\right)\ \mathrm{d}^2 r}
;
其中,
V
{\displaystyle \mathbb{V}}
是积分体积,
S
{\displaystyle \mathbb{S}}
是包住
V
{\displaystyle \mathbb{V}}
的闭表面,
d
2
r
{\displaystyle \mathrm{d}^2 r}
是微小面元素,
∂
ϕ
∂
n
{\displaystyle \partial\phi \over \partial n}
或
∂
ϕ
∂
n
{\displaystyle \partial\phi \over \partial n}
都是取垂直于闭表面
S
{\displaystyle \mathbb{S}}
的法向导数 ,都是从积分体积
V
{\displaystyle \mathbb{V}}
朝外指出。
设定
ϕ
(
r
′
)
{\displaystyle \phi(\mathbf{r}') }
为在
r
′
{\displaystyle \mathbf{r}' }
的电势,
ψ
=
1
|
r
−
r
′
|
{\displaystyle \psi=\frac{1}{|\mathbf{r} - \mathbf{r}'|}}
为
r
′
{\displaystyle \mathbf{r}' }
与
r
{\displaystyle \mathbf{r} }
之间的距离。应用泊松方程
∇
2
ϕ
(
r
)
=
−
ρ
/
ϵ
0
{\displaystyle \nabla^2\phi(\mathbf{r})=-\rho/\epsilon_0 }
,则可得到
∫
V
′
[
ϕ
(
r
′
)
∇
2
(
1
|
r
−
r
′
|
)
+
ρ
(
r
′
)
ϵ
0
|
r
−
r
′
|
]
d
3
r
′
=
∮
S
′
[
ϕ
∂
∂
n
′
(
1
|
r
−
r
′
|
)
−
(
1
|
r
−
r
′
|
)
∂
ϕ
∂
n
′
]
d
2
r
′
{\displaystyle \int_{\mathbb{V}'} \left[ \phi(\mathbf{r}') \nabla^2 \left(\frac{1}{|\mathbf{r} - \mathbf{r}'|}\right) + \frac{\rho(\mathbf{r}')}{\epsilon_0 |\mathbf{r} - \mathbf{r}'|} \right]\mathrm{d}^3 r' = \oint_{\mathbb{S}'} \left[ \phi\ {\partial \over \partial n'}\left(\frac{1}{|\mathbf{r} - \mathbf{r}'|}\right) - \left(\frac{1}{|\mathbf{r} - \mathbf{r}'|}\right) {\partial \phi \over \partial n'}\right]\mathrm{d}^2 r'}
。
再应用矢量恒等式
∇
2
(
1
|
r
−
r
′
|
)
=
−
4
π
δ
(
r
−
r
′
)
{\displaystyle \nabla^2 \left(\frac{1}{|\mathbf{r} - \mathbf{r}'|}\right)=-4\pi\delta(\mathbf{r} - \mathbf{r}')}
。
假设检验位置
r
{\displaystyle \mathbf{r}}
在积分体积
V
′
{\displaystyle \mathbb{V}'}
内,则可得到
ϕ
(
r
)
=
1
4
π
ϵ
0
∫
V
′
ρ
(
r
′
)
|
r
−
r
′
|
d
3
r
′
+
1
4
π
∮
S
′
[
(
1
|
r
−
r
′
|
)
∂
ϕ
∂
n
′
−
ϕ
∂
∂
n
′
(
1
|
r
−
r
′
|
)
]
d
2
r
′
{\displaystyle \phi(\mathbf{r}) =\frac{1}{4 \pi \epsilon_0} \int_{\mathbb{V}'} \frac{\rho(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|}\ \mathrm{d}^3 r' +\frac{1}{4 \pi}\oint_{\mathbb{S}'} \left[\left(\frac{1}{|\mathbf{r} - \mathbf{r}'|}\right) {\partial \phi \over \partial n'}-\phi\ {\partial \over \partial n'}\left(\frac{1}{|\mathbf{r} - \mathbf{r}'|}\right)\right]\mathrm{d}^2 r'}
。
这方程右手边的体积分就是电势的方程(1),而面积分就是因为边界条件 而添加的项目。这是
V
′
{\displaystyle \mathbb{V}'}
体内与体外之间的边界曲面。面积分的第一个项目要求给定在边界曲面的法向电场,即
E
n
′
=
−
∂
ϕ
∂
n
′
{\displaystyle E_{n'}= -{\partial \phi \over \partial n'}}
,也就是面感应电荷密度
σ
=
ϵ
0
E
n
′
{\displaystyle \sigma=\epsilon_0 E_{n'}}
。面积分的第二个项目要求给定在边界曲面的电势
ϕ
{\displaystyle \phi}
。假若能够知道积分体积内的电荷密度、在闭曲面的面电荷密度与电势,就可以计算出在积分体积内任意位置的电势。
根据柯西边界条件 ,有时候,给定在边界曲面的法向电场与电势,可能会因为给定过多边界条件,而造成无法计算出一致的电势的状况。实际而言,只要给定法向电场或电势,两者之一,就可以计算出电势。[5]
假若积分体积为无穷大空间,当
r
′
{\displaystyle r'}
趋向于无穷大时,则面积分的被积分项目会以
1
/
r
′
3
{\displaystyle 1/r'^3}
速率递减,而积分面积会以
r
′
2
{\displaystyle r'^2}
速率递增,所以,面积分项目会趋向于零,这方程约化为先前的电势方程(1)。
格林函数
包括函数
1
/
|
r
−
r
′
|
{\displaystyle 1/|\mathbf{r} - \mathbf{r}'|}
在内,有一类函数
G
(
r
,
r
′
)
{\displaystyle G(\mathbf{r},\mathbf{r}')}
,称为格林函数 ,能够满足方程
∇
2
G
(
r
,
r
′
)
=
−
4
π
δ
(
r
−
r
′
)
{\displaystyle \nabla^2 G(\mathbf{r},\mathbf{r}')= - 4\pi\delta(\mathbf{r} - \mathbf{r}')}
。
另外,假设函数
H
(
r
,
r
′
)
{\displaystyle H(\mathbf{r},\mathbf{r}')}
满足拉普拉斯方程
∇
2
H
(
r
,
r
′
)
=
0
{\displaystyle \nabla^2 H(\mathbf{r},\mathbf{r}')=0}
,
则函数
G
′
(
r
,
r
′
)
=
G
(
r
,
r
′
)
+
H
(
r
,
r
′
)
{\displaystyle G'(\mathbf{r},\mathbf{r}')=G(\mathbf{r},\mathbf{r}')+H(\mathbf{r},\mathbf{r}')}
也是格林函数。
应用这灵活性质,可以更严格地规定格林函数:[5]
对于狄利克雷问题 ,当源位置
r
′
{\displaystyle \mathbf{r}'}
在边界表面
S
′
{\displaystyle {\mathbb{S}'}}
时,规定格林函数
G
D
(
r
,
r
′
)
=
0
{\displaystyle G_D(\mathbf{r},\mathbf{r}')=0}
。这样,从格林第二恒等式,设定
ϕ
(
r
′
)
{\displaystyle \phi(\mathbf{r}') }
为在
r
′
{\displaystyle \mathbf{r}' }
的电势,
ψ
(
r
,
r
′
)
=
G
D
(
r
,
r
′
)
{\displaystyle \psi(\mathbf{r},\mathbf{r}')=G_D(\mathbf{r},\mathbf{r}')}
,则可得到
ϕ
(
r
)
=
1
4
π
ϵ
0
∫
V
′
ρ
(
r
′
)
G
D
(
r
,
r
′
)
d
3
r
′
−
1
4
π
∮
S
′
ϕ
(
r
′
)
∂
G
D
(
r
,
r
′
)
∂
n
′
d
2
r
′
{\displaystyle \phi(\mathbf{r}) =\frac{1}{4 \pi \epsilon_0} \int_{\mathbb{V}'} \rho(\mathbf{r}')G_D(\mathbf{r},\mathbf{r}')\ \mathrm{d}^3 r'
-\ \frac{1}{4 \pi}\oint_{\mathbb{S}'} \phi(\mathbf{r}')\ {\partial G_D(\mathbf{r},\mathbf{r}')\over \partial n'}\mathrm{d}^2 r'}
。(2)
对于满足纽曼问题 ,当源位置
r
′
{\displaystyle \mathbf{r}'}
在边界表面
S
′
{\displaystyle {\mathbb{S}'}}
时,规定格林函数
∮
S
′
∂
G
D
(
r
,
r
′
)
∂
n
′
d
2
r
′
=
−
4
π
S
{\displaystyle \oint_{\mathbb{S}'}\frac{\partial G_D(\mathbf{r},\mathbf{r}')}{\partial n'}\mathrm{d}^2 r'=-\frac{4\pi}{S}}
。
这两种规定都能够唯一地设定格林函数。注意到格林函数是一个几何函数,与整个系统的电荷分布无关。对于任何系统,只要计算出适合其几何形状的格林函数,则不论系统的电荷分布为何,都可以使用同样的格林函数。
无限平面导体案例
位于xy-平面的是一个接地 的无限平面导体。其上方的点电荷
q
{\displaystyle q}
的直角坐标是
(
0
,
0
,
a
)
{\displaystyle (0,\,0,\,a)}
。
假设xy-平面是接地 的无限平面导体 ,则对于z+ 半空间、满足狄利克雷边界条件的格林函数为
G
D
(
r
,
r
′
)
=
1
(
x
−
x
′
)
2
+
(
y
−
y
′
)
2
+
(
z
−
z
′
)
2
−
1
(
x
−
x
′
)
2
+
(
y
−
y
′
)
2
+
(
z
+
z
′
)
2
{\displaystyle \begin{matrix}G_D(\mathbf{r},\mathbf{r}')=\cfrac{1}{\sqrt{(x-x')^2+(y-y')^2+ (z - z')^2}} \\
\qquad\qquad\qquad -\ \cfrac{1}{\sqrt{(x-x')^2+(y-y')^2+(z+z')^2}} \\
\end{matrix}}
;
其中,
(
x
,
y
,
z
)
{\displaystyle (x,y,z)}
、
(
x
′
,
y
′
,
z
′
)
{\displaystyle (x',y',z')}
分别是检验位置
r
{\displaystyle \mathbf{r}}
、源位置
r
′
{\displaystyle \mathbf{r}'}
的直角坐标 。
由于接地导体的电势为零,方程(2)的面积分项目等于零,方程(2)变为
ϕ
(
r
)
=
1
4
π
ϵ
0
∫
V
′
ρ
(
r
′
)
G
D
(
r
,
r
′
)
d
3
r
′
{\displaystyle \phi(\mathbf{r}) =\frac{1}{4 \pi \epsilon_0} \int_{\mathbb{V}'} \rho(\mathbf{r}')G_D(\mathbf{r},\mathbf{r}')\ \mathrm{d}^3 r'
}
。
假设在位置
(
0
,
0
,
a
)
{\displaystyle (0,0,a)}
有点电荷
q
{\displaystyle q}
,则在z+ 半空间任意位置的电势为
ϕ
(
r
)
=
1
4
π
ϵ
0
∫
V
′
ρ
(
r
′
)
(
1
x
2
+
y
2
+
(
z
−
a
)
2
−
1
x
2
+
y
2
+
(
z
+
a
)
2
)
d
3
r
′
=
1
4
π
ϵ
0
(
q
x
2
+
y
2
+
(
z
−
a
)
2
−
q
x
2
+
y
2
+
(
z
+
a
)
2
)
{\displaystyle \begin{align}
\phi(\mathbf{r}) & =\frac{1}{4 \pi \epsilon_0} \int_{\mathbb{V}'} \rho(\mathbf{r}')\left( \frac{1}{\sqrt{x^2+y^2+ (z - a)^2}} - \frac{1}{\sqrt{x^2+y^2+ (z +a)^2}} \right)\ \mathrm{d}^3 r' \\
& =\frac{1}{4 \pi \epsilon_0}\left( \frac{q}{\sqrt{x^2+y^2+ (z - a)^2}} - \frac{q}{\sqrt{x^2+y^2+ (z+a)^2}} \right) \\
\end{align}}
。
仔细检察这方程,右手边第一个项目,是在没有平面导体的状况时,点电荷
q
{\displaystyle q}
所产生的电势;右手边第二个项目,是使用镜像法 时,镜像电荷
−
q
{\displaystyle -q}
所产生的电势。请参阅镜像法 条目的点电荷与无限平面导体 段落。
导引
已知函数
1
/
|
r
−
r
′
|
{\displaystyle 1/|\mathbf{r} - \mathbf{r}'|}
为格林函数
G
(
r
,
r
′
)
{\displaystyle G(\mathbf{r},\mathbf{r}')}
,满足方程
∇
2
G
(
r
,
r
′
)
=
−
4
π
δ
(
r
−
r
′
)
{\displaystyle \nabla^2 G(\mathbf{r},\mathbf{r}')= - 4\pi\delta(\mathbf{r} - \mathbf{r}')}
。
在三维无限空间里,
1
/
|
r
−
r
′
|
{\displaystyle 1/|\mathbf{r} - \mathbf{r}'|}
的傅里叶级数 为[6]
1
|
r
−
r
′
|
≡
1
2
π
2
∫
−
∞
∞
∫
−
∞
∞
∫
−
∞
∞
d
3
k
e
i
k
⋅
(
r
−
r
′
)
k
2
=
1
2
π
2
∫
−
∞
∞
∫
−
∞
∞
d
k
x
d
k
y
e
i
k
x
(
x
−
x
′
)
+
i
k
y
(
y
−
y
′
)
∫
−
∞
∞
d
k
z
e
i
k
z
(
z
−
z
′
)
k
x
2
+
k
y
2
+
k
z
2
{\displaystyle \begin{align}\frac{1}{|\mathbf{r} - \mathbf{r}'|}
& \equiv \frac{1}{2\pi^2}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\mathrm{d}^3 k \frac{e^{i\mathbf{k}\cdot(\mathbf{r}-\mathbf{r}')}}{k^2} \\
& = \frac{1}{2\pi^2}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\mathrm{d}k_x\ \mathrm{d}k_y e^{ik_x(x-x')+ik_y(y-y')}\int_{-\infty}^{\infty}\ \mathrm{d}k_z \frac{e^{ik_z(z-z')}}{k_x^2+k_y^2+k_z^2} \\
\end{align}}
。
现在,必需找到格林函数
G
D
(
r
,
r
′
)
=
G
(
r
,
r
′
)
+
H
(
r
,
r
′
)
{\displaystyle G_D(\mathbf{r},\mathbf{r}')=G(\mathbf{r},\mathbf{r}')+H(\mathbf{r},\mathbf{r}')}
,满足狄利克雷边界条件
G
D
(
(
x
,
y
,
0
)
,
r
′
)
=
0
{\displaystyle G_D((x,y,0),\mathbf{r}')=0}
,同时,函数
H
(
r
,
r
′
)
{\displaystyle H(\mathbf{r},\mathbf{r}')}
满足拉普拉斯方程
∇
2
H
(
r
,
r
′
)
=
0
{\displaystyle \nabla^2 H(\mathbf{r},\mathbf{r}')=0}
。
对于z+ 半空间,
H
(
r
,
r
′
)
{\displaystyle H(\mathbf{r},\mathbf{r}')}
以傅里叶级数 扩张为
H
(
r
,
r
′
)
=
1
2
π
2
∫
−
∞
∞
∫
−
∞
∞
d
k
x
d
k
y
e
i
k
x
(
x
−
x
′
)
+
i
k
y
(
y
−
y
′
)
∫
−
∞
∞
d
k
z
[
B
(
k
,
z
′
)
e
i
k
z
z
+
C
(
k
,
z
′
)
e
−
i
k
z
z
]
{\displaystyle \begin{align}H(\mathbf{r},\mathbf{r}')
& = \frac{1}{2\pi^2}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\ \mathrm{d}k_x\ \mathrm{d}k_y e^{ik_x(x-x')+ik_y(y-y')}\int_{-\infty}^{\infty}\ \mathrm{d}k_z \left[B(\mathbf{k},z')e^{ik_z z} +C(\mathbf{k},z')e^{-ik_z z}\right]\\
\end{align}}
。
对于x-坐标与对于y-坐标的傅里叶级数 扩张,
H
{\displaystyle H}
函数与
G
{\displaystyle G}
函数的形式相同。这是因为对于无限空间案例与无限平面导体案例,两种案例的x-边界条件与y-边界条件都相同,只有z-边界条件稍有改变。将
H
{\displaystyle H}
函数的方程代如,
G
D
(
r
,
r
′
)
{\displaystyle G_D(\mathbf{r},\mathbf{r}')}
变为
G
D
(
r
,
r
′
)
=
1
2
π
2
∫
−
∞
∞
∫
−
∞
∞
d
k
x
d
k
y
e
i
k
x
(
x
−
x
′
)
+
i
k
y
(
y
−
y
′
)
∫
−
∞
∞
d
k
z
[
e
i
k
z
(
z
−
z
′
)
k
x
2
+
k
y
2
+
k
z
2
+
B
(
k
,
z
′
)
e
i
k
z
z
+
C
(
k
,
z
′
)
e
−
i
k
z
z
]
{\displaystyle G_D(\mathbf{r},\mathbf{r}')
= \frac{1}{2\pi^2}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\ \mathrm{d}k_x\ \mathrm{d}k_y e^{ik_x(x-x')+ik_y(y-y')}\int_{-\infty}^{\infty}\ \mathrm{d}k_z \left[\frac{e^{ik_z(z-z')}}{k_x^2+k_y^2+k_z^2} +B(\mathbf{k},z')e^{ik_z z} +C(\mathbf{k},z')e^{-ik_z z}\right]}
;
其中,
B
(
k
,
z
′
)
{\displaystyle B(\mathbf{k},z')}
与
C
(
k
,
z
′
)
{\displaystyle C(\mathbf{k},z')}
都是系数函数。
由于
G
D
(
(
x
,
y
,
0
)
,
r
′
)
=
0
{\displaystyle G_D((x,y,0),\mathbf{r}')=0}
,对于任意
k
{\displaystyle \mathbf{k}}
与
z
′
{\displaystyle z'}
,
B
(
k
,
z
′
)
{\displaystyle B(\mathbf{k},z')}
与
C
(
k
,
z
′
)
{\displaystyle C(\mathbf{k},z')}
之间的关系为
e
−
i
k
z
z
′
k
x
2
+
k
y
2
+
k
z
2
+
B
(
k
,
z
′
)
+
C
k
,
z
′
)
=
0
{\displaystyle \frac{e^{-ik_z z'}}{k_x^2+k_y^2+k_z^2} +B(\mathbf{k},z')+C\mathbf{k},z')=0}
、
B
(
k
,
z
′
)
=
B
0
e
−
i
k
z
z
′
k
x
2
+
k
y
2
+
k
z
2
{\displaystyle B(\mathbf{k},z')=\frac{B_0 e^{-ik_z z'}}{k_x^2+k_y^2+k_z^2} }
、
C
(
k
,
z
′
)
=
C
0
e
−
i
k
z
z
′
k
x
2
+
k
y
2
+
k
z
2
{\displaystyle C(\mathbf{k},z')=\frac{C_0 e^{-ik_z z'}}{k_x^2+k_y^2+k_z^2} }
;
其中,
B
0
{\displaystyle B_0}
与
C
0
{\displaystyle C_0}
都是系数常数,而且,
B
0
+
C
0
=
−
1
{\displaystyle B_0+C_0=-1}
将这些公式代入
G
D
{\displaystyle G_D}
,可以得到
G
D
(
r
,
r
′
)
=
1
2
π
2
∫
−
∞
∞
∫
−
∞
∞
d
k
x
d
k
y
e
i
k
x
(
x
−
x
′
)
+
i
k
y
(
y
−
y
′
)
∫
−
∞
∞
d
k
z
{
(
1
+
B
0
)
k
2
[
e
i
k
z
(
z
−
z
′
)
−
e
i
k
z
(
z
+
z
′
)
]
}
{\displaystyle G_D(\mathbf{r},\mathbf{r}')
= \frac{1}{2\pi^2}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\ \mathrm{d}k_x\ \mathrm{d}k_y e^{ik_x(x-x')+ik_y(y-y')}\int_{-\infty}^{\infty}\ \mathrm{d}k_z \left\{\frac{(1+B_0)}{k^2} \left[ e^{ik_z(z-z')}-e^{ik_z(z+z')}\right]\right\}}
。
为了满足方程
∇
2
G
D
(
r
,
r
′
)
=
−
4
π
δ
(
r
−
r
′
)
{\displaystyle \nabla^2 G_D(\mathbf{r},\mathbf{r}')= - 4\pi\delta(\mathbf{r} - \mathbf{r}')}
,必需设定
B
0
=
0
{\displaystyle B_0=0}
。所以,
G
D
(
r
,
r
′
)
=
1
2
π
2
∫
−
∞
∞
∫
−
∞
∞
d
k
x
d
k
y
e
i
k
x
(
x
−
x
′
)
+
i
k
y
(
y
−
y
′
)
∫
−
∞
∞
d
k
z
{
1
k
2
[
e
i
k
z
(
z
−
z
′
)
−
e
i
k
z
(
z
+
z
′
)
]
}
=
1
|
r
−
r
′
|
−
1
|
r
−
r
″
|
=
1
(
x
−
x
′
)
2
+
(
y
−
y
′
)
2
+
(
z
−
z
′
)
2
−
1
(
x
−
x
′
)
2
+
(
y
−
y
′
)
2
+
(
z
+
z
′
)
2
{\displaystyle \begin{align}G_D(\mathbf{r},\mathbf{r}') & = \frac{1}{2\pi^2}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\ \mathrm{d}k_x\ \mathrm{d}k_y e^{ik_x(x-x')+ik_y(y-y')}\int_{-\infty}^{\infty}\ \mathrm{d}k_z \left\{\frac{1}{k^2} \left[ e^{ik_z(z-z')}-e^{ik_z(z+z')}\right]\right\} \\
& =\frac{1}{|\mathbf{r} - \mathbf{r}'|}-\frac{1}{|\mathbf{r} - \mathbf{r}''|} \\
& =\cfrac{1}{\sqrt{(x-x')^2+(y-y')^2+ (z - z')^2}} -\ \cfrac{1}{\sqrt{(x-x')^2+(y-y')^2+(z+z')^2}} \\
\end{align}}
;
其中,
r
″
=
(
x
′
,
y
′
,
−
z
′
)
{\displaystyle \mathbf{r}''=(x',y',-z')}
是镜像电荷的位置。
两个半平面导体案例
假设在xy-平面的无限平面导体 被一条位于
y
=
0
{\displaystyle y=0}
的绝缘线条分为两半,两个处于y+ 、y- -半平面的导体的电势分别设定为
+
V
{\displaystyle +V}
与
−
V
{\displaystyle -V}
,则由于
ρ
(
r
′
)
=
0
{\displaystyle \rho(\mathbf{r}')=0}
,方程(2)变为
ϕ
(
r
)
=
−
1
4
π
∮
S
′
ϕ
(
r
′
)
∂
G
D
(
r
,
r
′
)
∂
n
′
d
2
r
′
{\displaystyle \phi(\mathbf{r}) =-\ \frac{1}{4 \pi}\oint_{\mathbb{S}'} \phi(\mathbf{r}')\ {\partial G_D(\mathbf{r},\mathbf{r}')\over \partial n'}\mathrm{d}^2 r'}
。(3)
注意到
V
′
{\displaystyle \mathbb{V}'}
是z+ -半空间,xy-平面是其边界闭曲面的一部分,格林函数在xy-平面的法向导数的方向是朝着负z方向:
∂
G
D
∂
n
′
=
−
∂
G
D
∂
z
′
=
−
z
−
z
′
[
(
x
−
x
′
)
2
+
(
y
−
y
′
)
2
+
(
z
−
z
′
)
2
]
3
/
2
−
z
+
z
′
[
(
x
−
x
′
)
2
+
(
y
−
y
′
)
2
+
(
z
+
z
′
)
2
]
3
/
2
=
−
2
z
[
(
x
−
x
′
)
2
+
(
y
−
y
′
)
2
+
z
2
]
3
/
2
{\displaystyle \begin{align}{\partial G_D\over \partial n'} & =-\ {\partial G_D\over \partial z'} \\
& =-\ \cfrac{z - z'}{[(x-x')^2+(y-y')^2+ (z - z')^2]^{3/2}}\ -\ \cfrac{z + z'}{[(x-x')^2+(y-y')^2+(z+z')^2]^{3/2}} \\
& =-\ \cfrac{2z}{[(x-x')^2+(y-y')^2+ z^2]^{3/2}} \\
\end{align}}
。
V
′
{\displaystyle \mathbb{V}'}
的边界闭曲面在无穷远位置的电势为0,所以,只需要计算xy-平面给出的贡献,就可以得到在
V
′
{\displaystyle \mathbb{V}'}
内部任意位置的电势。将上述方程代入方程(3):[4]
ϕ
(
r
)
=
2
z
4
π
{
∫
0
+
∞
∫
−
∞
∞
V
d
x
′
d
y
′
[
(
x
−
x
′
)
2
+
(
y
−
y
′
)
2
+
z
2
]
3
/
2
+
∫
−
∞
0
−
∫
−
∞
∞
−
V
d
x
′
d
y
′
[
(
x
−
x
′
)
2
+
(
y
−
y
′
)
2
+
z
2
]
3
/
2
}
=
z
V
π
{
∫
0
+
∞
d
y
′
(
y
−
y
′
)
2
+
z
2
−
∫
−
∞
0
−
d
y
′
(
y
−
y
′
)
2
+
z
2
}
=
2
V
π
arctan
(
y
z
)
{\displaystyle \begin{align}\phi(\mathbf{r}) & =\frac{2z}{4 \pi}\left\{\int_{0+}^{\infty}\int_{-\infty}^{\infty}\cfrac{V \mathrm{d}x'\mathrm{d}y'
}{[(x-x')^2+(y-y')^2+ z^2]^{3/2}}+\int_{-\infty}^{0-}\int_{-\infty}^{\infty}\cfrac{-V \mathrm{d}x'\mathrm{d}y'
}{[(x-x')^2+(y-y')^2+ z^2]^{3/2}}\right\} \\
& =\ \frac{zV}{\pi}\left\{\int_{0+}^{\infty}\frac{\mathrm{d}y'}{(y-y')^2+ z^2}
-\int_{-\infty}^{0-} \frac{\mathrm{d}y'}{(y-y')^2+ z^2}\right\} \\
& =\frac{2V}{\pi}\ \arctan{\left(\frac{y}{z}\right)} \\
\end{align}}
。
推广至电动力学
假设磁场含时间(每当电场含时间,则此假设成立。逆过来亦成立),则不能简单地以标势
ϕ
{\displaystyle \phi}
描述电场。因为根据法拉第电磁感应定律 ,
∇
×
E
=
−
∂
B
∂
t
≠
0
{\displaystyle \mathbf{\nabla} \times \mathbf{E}=-\ \frac{\partial \mathbf{B}}{\partial t}\neq 0}
,电场不再具有保守性,
∫
E
⋅
d
ℓ
{\displaystyle \int \mathbf{E}\cdot \mathrm{d} \boldsymbol{\ell}}
跟路径有关。
替代地,在定义标势时,必须引入磁矢势
A
{\displaystyle \mathbf{A}}
,定义为
B
=
d
e
f
∇
×
A
{\displaystyle \mathbf{B}\ \stackrel{def}{=}\ \mathbf{\nabla} \times \mathbf{A}}
;
其中,
B
{\displaystyle \mathbf{B}}
是磁场 。
根据亥姆霍兹定理 [7] (Helmholtz theorem ),假设一个矢量函数
F
(
r
)
{\displaystyle \mathbf{F}(\mathbf{r})}
满足以下两条件:
∇
⋅
F
(
r
)
=
D
(
r
)
{\displaystyle \nabla\cdot\mathbf{F}(\mathbf{r})=D(\mathbf{r})}
、
∇
×
F
(
r
)
=
C
(
r
)
{\displaystyle \nabla\times\mathbf{F}(\mathbf{r})=\mathbf{C}(\mathbf{r})}
;
其中,
D
(
r
)
{\displaystyle D(\mathbf{r})}
是个标量函数,
C
(
r
)
{\displaystyle \mathbf{C}(\mathbf{r})}
是个矢量函数。
再假设
D
(
r
)
{\displaystyle D(\mathbf{r})}
和
C
(
r
)
{\displaystyle \mathbf{C}(\mathbf{r})}
,在无穷远处都足够快速地趋向0,则
F
(
r
)
{\displaystyle \mathbf{F}(\mathbf{r})}
可以用方程表达为
F
(
r
)
=
−
∇
(
1
4
π
∫
V
′
D
(
r
′
)
|
r
−
r
′
|
d
3
r
′
)
+
∇
×
(
1
4
π
∫
V
′
C
(
r
′
)
|
r
−
r
′
|
d
3
r
′
)
{\displaystyle \mathbf{F}(\mathbf{r})= - \nabla\left(\frac{1}{4\pi}\int_{\mathbb{V}'}\frac{D(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|} \mathrm{d}^3 r'\right)+\nabla\times \left(\frac{1}{4\pi}\int_{\mathbb{V}'}\frac{\mathbf{C}(\mathbf{r}')}{ |\mathbf{r} - \mathbf{r}'| }\mathrm{d}^3 r'\right)}
;
在这里,
∇
{\displaystyle \nabla}
只作用于
r
{\displaystyle \mathbf{r}}
,体积分的体积为
V
′
{\displaystyle \mathbb{V}'}
。
采用库仑规范 (Coulomb gauge ),则磁矢势
A
{\displaystyle \mathbf{A}}
遵守
∇
⋅
A
=
0
{\displaystyle \mathbf{\nabla} \cdot \mathbf{A}=0}
。
所以,
A
(
r
)
=
∇
×
(
1
4
π
∫
V
′
B
(
r
′
)
|
r
−
r
′
|
d
3
r
′
)
=
1
4
π
∫
V
′
B
(
r
′
)
×
r
−
r
′
|
r
−
r
′
|
3
d
3
r
′
{\displaystyle \mathbf{A}(\mathbf{r})=\nabla\times \left(\frac{1}{4\pi}\int_{\mathbb{V}'}\frac{\mathbf{B}(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|}\mathrm{d}^3 r'\right) =\frac{1}{4\pi}\int_{\mathbb{V}'}\mathbf{B}(\mathbf{r}')\times\frac{\mathbf{r} - \mathbf{r}'}{|\mathbf{r} - \mathbf{r}'|^3} \mathrm{d}^3 r'}
。
注意到,以上这些推导,并没有涉及时间参数。加入时间参数
t
{\displaystyle t}
,结果也成立。所以,永远可以找到磁矢势
A
{\displaystyle \mathbf{A}}
:
A
(
r
,
t
)
=
1
4
π
∫
V
′
B
(
r
′
,
t
)
×
r
−
r
′
|
r
−
r
′
|
3
d
3
r
′
{\displaystyle \mathbf{A}(\mathbf{r},\,t)=\frac{1}{4\pi}\int_{\mathbb{V}'} \mathbf{B}(\mathbf{r}',\,t)\times\frac{\mathbf{r} - \mathbf{r}'}{|\mathbf{r} - \mathbf{r}'|^3} \mathrm{d}^3 r' }
。
根据法拉第电磁感应定律 ,矢量场
G
=
E
+
∂
A
/
∂
t
{\displaystyle \mathbf{G}=\mathbf{E}+\partial\mathbf{A}/\partial t}
是一个保守场:
∇
×
G
=
∇
×
E
+
∇
×
∂
A
/
∂
t
=
0
{\displaystyle \nabla\times\mathbf{G}=\nabla\times\mathbf{E}+\nabla\times\partial\mathbf{A}/\partial t=0}
。
所以,必定可以找到标势
ϕ
{\displaystyle \phi}
,满足
G
=
−
∇
ϕ
{\displaystyle \mathbf{G} = - \nabla \phi}
。因此,下述方程成立:
E
=
−
∇
ϕ
−
∂
A
∂
t
{\displaystyle \mathbf{E} = - \mathbf{\nabla}\phi - \frac{\partial\mathbf{A}}{\partial t}}
。
静电势只是这含时定义的一个特别案例,在这案例里,磁矢势
A
{\displaystyle \mathbf{A}}
不含时间。从另一方面来说,对于含时矢量场,电场的路径积分与静电学的结果大不相同:
∫
a
b
E
⋅
d
ℓ
≠
ϕ
(
b
)
−
ϕ
(
a
)
{\displaystyle \int_a^b \mathbf{E} \cdot \mathrm{d}\boldsymbol{\ell} \neq \phi(b) - \phi(a)}
。
参阅
参考文献
↑ 電勢 . 中华语文知识库. [2016-03-03 ] (简体中文(中国大陆)) .
↑ 2.0 2.1 Halliday, David; Robert Resnick, Jearl Walker, Fundamental of Physics 7th, USA: John Wiley and Sons, Inc.: pp. 630ff, 2005, ISBN 0-471-23231-9
↑ Jackson 1999 ,第70-72页
↑ 4.0 4.1 Beyer, William, CRC Standard Mathematical Table 28th, CRC Press, 1987, ISBN 0-8493-0628-0 pp. 241, formula #43,
pp. 252, formula#165
↑ 5.0 5.1 5.2 Jackson 1999 ,第35-40页
↑ Jackson 1999 ,第127-129页
↑ Griffiths, David J. Introduction to Electrodynamics 3rd. Prentice Hall. 1998: pp. 555–557. ISBN 0-13-805326-X .
Jackson, John David, Classical Electrodynamic 3rd., USA: John Wiley & Sons, Inc., 1999, ISBN 978-0-471-30932-1
Wangsness, Roald K. Electromagnetic Fields 2nd., Revised, illustrated. Wiley. 1986. ISBN 9780471811862 .