P04 · 第 2 章 · 第一编 静电学

电势、Poisson 方程与边值问题

由静电场的路径无关性定义电势和电势能,推导 Poisson 与 Laplace 方程,并以导体边界、唯一性定理、分离变量和镜像法求解典型静电边值问题。

报告页面错误
预备知识Coulomb 定律、电场与 Gauss 定律梯度偏微分方程

本章目标

  1. 从静电场的线积分定义电势差,保持路径方向、电荷符号与功的符号一致。
  2. 由点电荷势和连续电荷积分求电势,再用负梯度恢复带方向的电场。
  3. 从 Gauss 定律推导 Poisson 与 Laplace 方程,核对电势、场和电荷密度的 SI 单位。
  4. 陈述导体静电边界、Dirichlet 与 Neumann 条件,并用能量积分证明相应唯一性结论。
  5. 使用球对称解、镜像法和一维 Poisson 模型求解可复核的边值问题,说明方法的适用区域。
页面阅读位置0% · 仅保存在此浏览器
章节未开始
本册完成进度0/6 章 · 0%
本页目录

从力做功转向标量势

直接叠加 Coulomb 场时,每个源都贡献一个矢量。静电场还有更简洁的结构:固定电荷产生的场是保守场,闭合路径上的环流为零,

CEd=0.\oint_C\boldsymbol E\cdot\mathrm d\boldsymbol\ell=0.

因此两点之间的线积分与路径无关,可以用一个标量函数记录。以下约定路径从初点 AA 指向末点 BBd\mathrm d\boldsymbol\ell 沿积分方向;空间坐标采用右手系,曲面法向在每个边值问题中另行声明。

电势差与电势

静电势差定义为

V(B)V(A)=ABEd.V(B)-V(A)=-\int_A^B\boldsymbol E\cdot\mathrm d\boldsymbol\ell.

电势单位为伏特,1V=1JC11\,\mathrm V=1\,\mathrm{J\,C^{-1}}。选择参考点 r0\boldsymbol r_0 和参考值 V(r0)V(\boldsymbol r_0) 后,其他点的电势由上式确定。整体加同一常数不会改变电场。

负号表达正电荷在电场方向运动时电势降低。局部取一段位移 d\mathrm d\boldsymbol\ell,有 dV=Ed\mathrm dV=-\boldsymbol E\cdot\mathrm d\boldsymbol\ell,所以

E=V.\boldsymbol E=-\nabla V.

电势的等值面与电场正交,电场指向电势下降最快的方向。V\nabla V 的单位是 Vm1\mathrm{V\,m^{-1}},与 NC1\mathrm{N\,C^{-1}} 相同。

电势能、外力功与符号

把电荷 qq 放入给定外场,其电势能定义为 U=qVU=qV。从 AA 缓慢移到 BB 时,电场做功与势能变化满足

WE,AB=ΔU=q[V(A)V(B)].W_{E,A\to B}=-\Delta U=q[V(A)-V(B)].

若移动足够缓慢、动能变化可忽略,外力做功为 Wext=ΔUW_{\mathrm{ext}}=\Delta U。对负电荷,U=qVU=qV 会反转高低关系:电子自发加速的方向仍由 F=qE\boldsymbol F=q\boldsymbol E 决定,不能只看“电势高低”而遗漏电荷符号。

点电荷 qq 取无穷远为零势时,

V(r)=14πε0qrr.V(\boldsymbol r)=\frac{1}{4\pi\varepsilon_0} \frac{q}{\lVert\boldsymbol r-\boldsymbol r'\rVert}.

多个点电荷的电势直接作代数和。连续体电荷则给出

V(r)=14πε0ρ(r)rrdV.V(\boldsymbol r)=\frac{1}{4\pi\varepsilon_0} \iiint\frac{\rho(\boldsymbol r')}{\lVert\boldsymbol r-\boldsymbol r'\rVert} \,\mathrm dV'.

标量积分通常比逐分量求场容易,之后再以 V-\nabla V 恢复方向。对无限线电荷或无限平面,令无穷远电势为零并不成立,因为相应积分发散;此时只定义两个有限位置之间的势差。

例 1:电势为正而电场指向负 x

q1=+4.00nCq_1=+4.00\,\mathrm{nC} 位于 x=0.0300mx=-0.0300\,\mathrm mq2=1.00nCq_2=-1.00\,\mathrm{nC} 位于 x=0.0600mx=0.0600\,\mathrm m。求 x=0.120mx=0.120\,\mathrm m 处的电势和 xx 分量。两距离分别为 0.150m0.150\,\mathrm m0.0600m0.0600\,\mathrm m,因此

V=ke(4.00×1090.1501.00×1090.0600)89.9V.V=k_e\left(\frac{4.00\times10^{-9}}{0.150} -\frac{1.00\times10^{-9}}{0.0600}\right) \approx89.9\,\mathrm V.

电势为正不能决定场方向。该点位于两个电荷右侧,正电荷场向 +x+x,负电荷场向 x-x

Ex=ke(4.00×1090.15021.00×1090.06002)8.99×102NC1.E_x=k_e\left(\frac{4.00\times10^{-9}}{0.150^2} -\frac{1.00\times10^{-9}}{0.0600^2}\right) \approx-8.99\times10^2\,\mathrm{N\,C^{-1}}.

负电荷更近,主导局部斜率,所以电场向左。若放入 q0=+2.00nCq_0=+2.00\,\mathrm{nC},其电势能为 U=q0V1.80×107JU=q_0V\approx1.80\times10^{-7}\,\mathrm J;这仍不等于该点所受力,力要用 q0Eq_0\boldsymbol E 计算。

Poisson 方程与 Laplace 方程

E=V\boldsymbol E=-\nabla V 代入真空中的局部 Gauss 定律 E=ρ/ε0\nabla\cdot\boldsymbol E=\rho/\varepsilon_0,得到

2V=ρε0.\boxed{\nabla^2V=-\frac{\rho}{\varepsilon_0}}.

这就是 Poisson 方程。无体电荷区域满足

2V=0,\boxed{\nabla^2V=0},

称为 Laplace 方程。右侧 ρ/ε0\rho/\varepsilon_0 的单位为 Vm2\mathrm{V\,m^{-2}},与电势的二阶空间导数一致。无源不等于零场:平行板之间 ρ=0\rho=0,电势可以线性变化,电场仍为非零常量。

Laplace 方程的解称为调和函数。它在内部不能出现严格的孤立极大值或极小值,极值由边界控制。这个性质说明纯静电场不能在三维自由空间中构成稳定的正电荷势能最低点;若加入时间变化场、磁场或机械约束,结论的前提才会改变。

例 2:一维均匀电荷层中的 Poisson 解

0xL0\le x\le L 的平板区域内,体电荷密度为常量 ρ>0\rho>0,两端边界固定为 V(0)=V(L)=0VV(0)=V(L)=0\,\mathrm V。若场只依赖 xx,方程为

d2Vdx2=ρε0.\frac{\mathrm d^2V}{\mathrm dx^2}=-\frac{\rho}{\varepsilon_0}.

积分两次并代入边界条件,

V(x)=ρ2ε0x(Lx),Ex=dVdx=ρε0(xL2).V(x)=\frac{\rho}{2\varepsilon_0}x(L-x), \qquad E_x=-\frac{\mathrm dV}{\mathrm dx} =\frac{\rho}{\varepsilon_0}\left(x-\frac L2\right).

中心处电场为零而电势最大;左半区场向 x-x,右半区场向 +x+x,符合正电荷向两侧产生场的对称性。取 ρ=1.00μCm3\rho=1.00\,\mathrm{\mu C\,m^{-3}}L=0.0200mL=0.0200\,\mathrm m,中心电势为

V(L/2)=ρL28ε05.65V,V(L/2)=\frac{\rho L^2}{8\varepsilon_0}\approx5.65\,\mathrm V,

两端内侧场强大小均为 ρL/(2ε0)1.13×103Vm1\rho L/(2\varepsilon_0)\approx1.13\times10^3\,\mathrm{V\,m^{-1}}。边界外部由何种导体或外加电源维持零势,不由这段区域方程单独说明。

边界数据决定哪一个解

微分方程只限制区域内部,边界条件从可能解族中选出实际解。静电问题常见三类数据:

  • Dirichlet 条件给定边界上的 VV,例如接地导体满足 V=0V=0
  • Neumann 条件给定外法向导数 V/n=Vn^\partial V/\partial n=\nabla V\cdot\widehat{\boldsymbol n},等价于给定 En=V/nE_n=-\partial V/\partial n
  • Robin 条件给定 aV+bV/naV+b\,\partial V/\partial n,可描述某些有效界面或数值边界。

这里的 n^\widehat{\boldsymbol n} 默认从所求区域指向区域外部。若研究“导体外部真空区域”,其内边界上的区域外法向可能指向导体内部;而表面电荷公式常取从导体指向真空的法向。两种法向相反时,公式符号也相反,必须在图或文字中声明。

纯 Neumann 问题还需要相容条件。对 Poisson 方程积分并用散度定理,

ΩVndA=1ε0ΩρdV.\oint_{\partial\Omega}\frac{\partial V}{\partial n}\,\mathrm dA =-\frac{1}{\varepsilon_0}\iiint_\Omega\rho\,\mathrm dV.

边界法向导数若不满足这个总通量约束,就不存在解。即使存在,VV 仍可整体加常数;电场则唯一。

导体边界与表面电荷

静电平衡导体材料内部 E=0\boldsymbol E=0,所以导体为等势体。表面切向场连续且为零;真空侧的法向场由表面电荷决定。若法向 n^c\widehat{\boldsymbol n}_c 从导体指向真空,

Eoutn^c=σε0,Voutnc=σε0.\boldsymbol E_{\mathrm{out}}\cdot\widehat{\boldsymbol n}_c =\frac{\sigma}{\varepsilon_0}, \qquad \frac{\partial V_{\mathrm{out}}}{\partial n_c} =-\frac{\sigma}{\varepsilon_0}.

因此导体表面给定电势后,解出外部 VV 并求法向导数,就能反推出非均匀的 σ\sigma。曲率大的尖端附近常出现较大的法向梯度和表面电荷密度,但理想尖点可能导致数学发散,真实材料的有限曲率和场致发射会限制该模型。

在两种普通线性介质交界面,电势通常连续,切向 E\boldsymbol E 连续;法向电位移 D\boldsymbol D 的跳变由自由面电荷决定。不能把真空中的 σ/ε0\sigma/\varepsilon_0 公式不加区分地用于任意介质总电荷。

唯一性定理:同一边界为何不会给出两种场

V1V_1V2V_2 在同一区域满足同一 Poisson 方程,并在整个边界具有相同 Dirichlet 数据。令 u=V1V2u=V_1-V_2,则 2u=0\nabla^2u=0 且边界上 u=0u=0。由乘积求散度,

(uu)=u2+u2u.\nabla\cdot(u\nabla u)=|\nabla u|^2+u\nabla^2u.

在区域上积分并用散度定理:

Ωu2dV=ΩuundA=0.\iiint_\Omega|\nabla u|^2\,\mathrm dV =\oint_{\partial\Omega}u\frac{\partial u}{\partial n}\,\mathrm dA=0.

体积分的被积函数非负,故 u=0\nabla u=0uu 是常数,再由边界 u=0u=0u0u\equiv0。因此 Dirichlet 问题的电势唯一。纯 Neumann 条件下同样可证明 u=0\nabla u=0,但常数不能由法向导数确定,所以只有电场唯一。

这一证明也是检验构造法的依据:镜像法并未真实地在导体另一侧安放电荷;只要构造出的势在所求区域满足 Poisson 方程和正确边界数据,唯一性保证它就是该区域的解。

分离变量如何把边界形状变成模态

规则矩形区域常可用分离变量求 Laplace 方程。以半无限条带 0<x<L0<x<Ly>0y>0 为例,令两侧 x=0,Lx=0,L 接地,远离底边时 V0V\to0,底边给定 V(x,0)=V0sin(πx/L)V(x,0)=V_0\sin(\pi x/L)。设 V=X(x)Y(y)V=X(x)Y(y),代入 Vxx+Vyy=0V_{xx}+V_{yy}=0,可把同一分离常数写成

X+k2X=0,Yk2Y=0.X''+k^2X=0, \qquad Y''-k^2Y=0.

两侧零势选择 X=sin(nπx/L)X=\sin(n\pi x/L),远处衰减排除增长指数。当前边界只含 n=1n=1 模态,所以

V(x,y)=V0sin(πxL)eπy/L.V(x,y)=V_0\sin\left(\frac{\pi x}{L}\right) e^{-\pi y/L}.

由负梯度得到

Ex=πV0Lcos(πxL)eπy/L,Ey=πV0Lsin(πxL)eπy/L.E_x=-\frac{\pi V_0}{L}\cos\left(\frac{\pi x}{L}\right)e^{-\pi y/L}, \qquad E_y=\frac{\pi V_0}{L}\sin\left(\frac{\pi x}{L}\right)e^{-\pi y/L}.

直接求二阶导数可核对 2V=0\nabla^2V=0,代入三处边界也逐项成立。一般底边函数可展开为正弦级数,每个横向波数 nπ/Ln\pi/L 都以 enπy/Le^{-n\pi y/L} 衰减;边界上的细小尺度对应较大的 nn,离开边界后消失得更快。这一结论来自给定矩形几何,换成圆柱或球面时应改用相应坐标的本征函数。

球形电容器:用边界选出 Laplace 解

例 3:同心球导体之间的电势与电容

内导体半径 a=0.0200ma=0.0200\,\mathrm m,外导体内半径 b=0.0500mb=0.0500\,\mathrm m。真空区 a<r<ba<r<b 内无体电荷,边界为 V(a)=100VV(a)=100\,\mathrm VV(b)=0VV(b)=0\,\mathrm V。球对称 Laplace 方程为

1r2ddr(r2dVdr)=0,\frac{1}{r^2}\frac{\mathrm d}{\mathrm dr} \left(r^2\frac{\mathrm dV}{\mathrm dr}\right)=0,

积分得 V=A+B/rV=A+B/r。边界给出

B=V(a)V(b)1/a1/b=3.333Vm,A=Bb=66.67V.B=\frac{V(a)-V(b)}{1/a-1/b} =3.333\,\mathrm{V\,m}, \qquad A=-\frac Bb=-66.67\,\mathrm V.

所以 E=(B/r2)e^r\boldsymbol E=(B/r^2)\widehat{\boldsymbol e}_r,从高势内球指向接地外球。在内球外侧 E(a)=8.33×103Vm1E(a)=8.33\times10^3\,\mathrm{V\,m^{-1}},在外球内侧 E(b)=1.33×103Vm1E(b)=1.33\times10^3\,\mathrm{V\,m^{-1}}。由 E=Q/(4πε0r2)E=Q/(4\pi\varepsilon_0r^2)

Q=4πε0B3.71×1010C,C=Q100V3.71pF.Q=4\pi\varepsilon_0B\approx3.71\times10^{-10}\,\mathrm C, \qquad C=\frac{Q}{100\,\mathrm V}\approx3.71\,\mathrm{pF}.

同一几何的一般电容为 C=4πε0ab/(ba)C=4\pi\varepsilon_0ab/(b-a)。电容只依赖几何和介质,不随这里选择的 100V100\,\mathrm V 改变,前提是材料响应保持线性且没有击穿。

镜像法:满足边界的等效构造

例 4:点电荷与接地无限导体平面

真空区域为 z>0z>0,接地导体占 z0z\le0。点电荷 q=+2.00nCq=+2.00\,\mathrm{nC} 位于 (0,0,h)(0,0,h)h=0.0400mh=0.0400\,\mathrm m。在所求区域内,用位于 (0,0,h)(0,0,-h) 的虚拟电荷 q-q 构造

V(ρ,z)=q4πε0[1ρ2+(zh)21ρ2+(z+h)2].V(\rho,z)=\frac{q}{4\pi\varepsilon_0} \left[ \frac{1}{\sqrt{\rho^2+(z-h)^2}} -\frac{1}{\sqrt{\rho^2+(z+h)^2}} \right].

z=0z=0 两项相消,满足接地边界;z>0z>0 除真实电荷位置外满足 Laplace 方程。唯一性定理保证它给出导体外部真实电势。真实电荷所受力可以按它与镜像电荷之间距离 2h2h 计算:

F=e^zkeq2(2h)25.62×106e^zN.\boldsymbol F =-\widehat{\boldsymbol e}_z \frac{k_eq^2}{(2h)^2} \approx-5.62\times10^{-6}\widehat{\boldsymbol e}_z\,\mathrm N.

负号表示吸向导体。导体表面取从导体指向真空的法向 +e^z+\widehat{\boldsymbol e}_z,由外侧法向场可得

σ(ρ)=qh2π(ρ2+h2)3/2.\sigma(\rho)=-\frac{qh}{2\pi(\rho^2+h^2)^{3/2}}.

σ\sigma 处处为负,积分总量为 q-q。镜像电荷只用于 z>0z>0 的数学构造;不能拿该势描述导体内部,也不能把镜像的自能直接计入真实空间。

电场能量与边值解的物理读法

电容系统从无电荷缓慢充到 QQ,若电压随电荷线性变化为 V=q/CV=q/C,外力做功储存在场中:

U=0QqCdq=Q22C=12CV2=12QV.U=\int_0^Q\frac qC\,\mathrm dq =\frac{Q^2}{2C}=\frac12CV^2=\frac12QV.

真空静电场的能量密度为

uE=12ε0E2,u_E=\frac12\varepsilon_0E^2,

单位为 Jm3\mathrm{J\,m^{-3}}。积分 uEdV\iiint u_E\,\mathrm dV 可与电容公式相互核对。点电荷的自能积分在 r0r\to0 时发散,表明经典点粒子模型不能给出有限自能;对有限尺寸带电体,需从实际半径开始积分。

常见失误与适用范围

常见误区

V=0V=0 的地方电场一定为零。”电势值依赖参考零点,电场取决于空间梯度。接地导体表面 V=0V=0,紧邻表面的法向场可以很强。

常见误区

“无电荷区满足 Laplace 方程,所以电势处处为常数。”线性函数也满足二阶导数为零;边界电势不同就会维持非零场。

常见误区

“镜像电荷是真实感应电荷集中成的点。”真实导体上的感应电荷是连续面分布。镜像只在指定求解区域复现相同边界电势。

本章采用静电近似、理想导体和真空介电常数。时变场满足 ×E=B/t\nabla\times\boldsymbol E=-\partial\boldsymbol B/\partial t,一般不能用单一全局标量势表示全部电场;非线性介质、空间变化介电常数和有限电导还会改变方程与界面条件。

练习:电势、方程与边界

练习

匀强电场 E=250e^xVm1\boldsymbol E=250\widehat{\boldsymbol e}_x\,\mathrm{V\,m^{-1}}。把 q=3.00nCq=-3.00\,\mathrm{nC}x=0.10mx=0.10\,\mathrm m 缓慢移到 x=0.34mx=0.34\,\mathrm m,求势差、电势能变化和电场做功。

查看提示
沿匀强场方向移动时 EdE\cdot dℓ 为正;先求 ΔV\Delta V,再分别使用 ΔU=qΔV\Delta U=q\Delta VWE=ΔUW_E=-\Delta U
查看解答
场为 E=250Vm1E=250\,\mathrm{V}\cdot m^{-1} 沿 +x,从 x=0.10mx=0.10\,\mathrm{m} 到 0.34 m 有 ΔV=EΔx=60.0V\Delta V=-E\Delta x=-60.0\,\mathrm{V}。对 q=3.00nCq=-3.00\,\mathrm{nC}ΔU=qΔV=+1.80×107J\Delta U=q\Delta V=+1.80\times 10^{-7} J,电场做功 WE=1.80×107JW_E=-1.80\times 10^{-7} J。负电荷被外力沿场方向缓慢移动时势能增加。
练习

给定 V(x,y,z)=α(x2+y2)E0zV(x,y,z)=\alpha(x^2+y^2)-E_0z。求 E\boldsymbol Eρ\rho,并说明横向等势线和场方向。写出 α\alphaE0E_0 应具有的 SI 单位。

查看提示
分别计算三个偏导,再整体加负号;等势面法向与梯度同向。
查看解答
V=(2αx,2αy,E0)\nabla V=(2\alpha x,2\alpha y,-E_0),故 E=(2αx,2αy,E0)E=(-2\alpha x,-2\alpha y,E_0)。在圆柱面 x2+y2=x^{2}+y^{2}=常量上,横向场沿径向;z 项给均匀 +z 场。代入 Laplace 算子得 2V=4α\nabla^{2}V=4\alpha,所以区域含 ρ=4ϵ0α\rho=-4\epsilon_0\alpha。量纲为 [α]=Vm2[\alpha]=V\cdot m^{-2}[E0]=Vm1[E_0]=V\cdot m^{-1}
练习

0<x<L0<x<L 的无电荷区域,边界为 V(0)=V0V(0)=V_0V(L)=VLV(L)=V_L。求 V(x)V(x)E\boldsymbol E,并根据 V0>VLV_0>V_L 判断方向。

查看提示
无源一维区域满足 V=0V''=0,先写一般线性解,再代入两端电势。
查看解答
V(x)=V0+(VLV0)x/LV(x)=V_0+(V_L-V_0)x/LEx=(VLV0)/L=(V0VL)/LE_x=-(V_L-V_0)/L=(V_0-V_L)/L,为常量。若 V0>VLV_0>V_L,场沿 +x。任意整体加常数会同时改变两端给定值,因此 Dirichlet 数据固定后不能再自由平移。
练习

同心球电容器 a=1.00cma=1.00\,\mathrm{cm}b=3.00cmb=3.00\,\mathrm{cm},两球电势差为 500V500\,\mathrm V。求电容、电荷量和场能,并说明内球高势时的场方向。

查看提示
先用一般式 C=4πϵ0ab/(ba)C=4\pi \epsilon_0ab/(b-a),再以 U=CV2/2U=CV^{2}/2;半径必须换成米。
查看解答
a=0.010ma=0.010\,\mathrm{m}b=0.030mb=0.030\,\mathrm{m}C=4πϵ0(0.010×0.030)/(0.020)=1.67pFC=4\pi \epsilon_0(0.010\times 0.030)/(0.020)=1.67 pF。电压 500 V 时 Q=CV=8.34×1010CQ=CV=8.34\times 10^{-10} CU=CV2/2=2.09×107JU=CV^{2}/2=2.09\times 10^{-7} J。内球电势较高时内球带正电,场沿径向向外。
练习

正点电荷 qq 距接地无限平面 hh。求其受力,并写出从无穷远缓慢移到该处时系统能量。解释为何不能把真实电荷与镜像电荷的配对能原样当作系统能量。

查看提示
镜像电荷到真实电荷的距离是 2h;势能不能直接使用两个真实电荷体系的完整配对能。
查看解答
接地平面镜像为 -q,真实电荷受力大小 keq2/(4h2)k_eq^{2}/(4h^{2}),方向垂直指向平面。把电荷从无穷远准静态移到 h,系统能量为 keq2/(4h)-k_eq^{2}/(4h),其对 h 的负导数给出同一吸引力;它是 q 乘镜像势的一半。
练习

证明纯 Neumann Poisson 问题必须满足总通量相容条件。再判断:无体电荷的有限区域若在整个边界指定同一个正常数 V/n=g>0\partial V/\partial n=g>0,是否存在解。

查看提示
对 Poisson 方程在区域上积分,并把 Laplace 算子改写为边界外法向导数积分。
查看解答
必要条件是 ∮Ω(V/n)dA=(1/ϵ0)\partial \Omega(\partial V/\partial n)dA=-(1/\epsilon_0)ΩρdV\Omega \rho dV。若 ρ=0\rho=0 而边界给出的外法向导数处处为正常数 g>0,则左侧为 gAboundary>0gA_{boundary}>0,不可能满足零右侧,因此无解。满足条件时 V 仍只确定到加法常数。

知识连接与继续学习

课程 · 2007

Physics II: Electricity and Magnetism

用于核对 P04 的场与势定义、积分定律方向、边界条件、电路暂态、Maxwell 方程和波动例题。

打开官方来源

MIT OpenCourseWare 8.02 的电势与导体单元覆盖势差、电势能、等势面、导体、电容和典型边界问题,可用于核对本章的物理符号与算例路线。继续学习时应先掌握稳恒电流和磁场,再把静电势的边界方法与磁矢势、Faraday 感应及完整 Maxwell 方程并列比较。