Physical Remedy

Gauss 定理的积分式可以从球面 Coulomb 定理推知,即

EdS=ρdVε0EdS=ρdVε0

对左侧使用微积分的 Gauss 定理,有

(E)dV=ρdVε0(E)dV=ρdVε0

令体积趋于 00 即由连续性得

E=ρε0E=ρε0

即 Gauss 定律微分式。


试探电荷能否在静电场中保持 力学平衡(mechanical equilibrium)?

  • 即,该点处电场为零;给定微扰后,有 回复力(restoring force)令其向位移的反方向运动。

这等效于,所有电场线均指向该点。使用 Gauss 定律分析后,得出唯一的可能性是:在与该点重合的地方,有一个电性相反的点电荷。【但是这个反电性的电荷又要靠谁束缚呢?】

这可以等效于,纯靠电磁力不可能存在力学平衡态。使用管道限制鞍点的鞍向,可以存在力学平衡态;同理,使用交变电场可以产生电阱束缚之;同理,使用光照可以束缚之。


对于电偶极子(electric dipole),位矢 RR 处的电势为 φ=pR4πε0R3φ=pR4πε0R3;电场为 E=p4πε0=3cos2θ+1R3E=p4πε0=3cos2θ+1R3,其中 θθRRpp 的夹角。

  • 也可以写成一坨的 E=14πε0(3(pR)RR5pR3)E=14πε0(3(pR)RR5pR3)。特别地,中垂线上的 RR,前一项为零,式子退化为 E=14πε0pR3E=14πε0pR3

电偶极子的另一种理解方式:

初始两电荷重合抵消,全空间无电场、无电势。

现在,令正电荷沿 zz 方向移动了小位移 ΔzΔz。在某点产生的电势,相当于该场点在正电荷电场中移动小位移 ΔzΔz 的电势变化。

于是,将电偶极子看作正电荷移动 l2l2,负电荷移动 l2l2,二者产生的电势变化并非抵消而是同号叠加。于是

Δφ=Δφ++Δφ=q4πε01R(l2)q4πε01Rl2=14πε0p1R

注意,此处推导依赖于 Rl 的前提,此时 R 附近两电荷的电场均可以直接用 q4πε01R 近似表达。

1R=RR3,于是

φ=14πε0pRR3


计算带电体 VR 处的电势,为

φ=Vρ(r)dV4πε0r

其中 r 是被积点的位矢,r=|Rr|

Rrr,则 Taylor 展开可知 1r=1Rr1R+12!i,jrirj2rirj1R+,于是

φ=14πε0ρ(r)(1Rr1R+12!i,jrirj2rirj1R+)dV

其中,零次项即将其当作点电荷,一次项当作电偶极子,二次项是电四极子……

即,

φ=14πε0QRp1R+16i,jDi,j2rirj1R+

其中

Q=ρ(r)dVp=ρ(r)rdVDi,j=3rirjρ(r)dV

电四极子衡量球面分布,球对称的电荷分布没有四极矩。两个重合的偶极子反接重合,无电场;然后偶极子各自平移,产生电场即为四极子。四极矩可以由偶极子电势微分得到。

H2O 等极性分子,可以被当作偶极子(总电荷为零,但偶极矩非零);CO2 等非极性分子,可以当成两个偶极子重叠构成的四极子(总电荷、偶极矩均为零,只能考虑四极矩)。


导体的场合。导体产生的电场不好研究,因为导体电子在外电场影响下会重分布。

  • 咋办?猜分布、验证是否是等势面!由唯一性定理,只要猜中就是唯一解!

无导体时的唯一性定理:

在体积 V 中,若已知电荷分布 ρ(x)【注意:ρ(x) 是固定电荷分布,不是导体中的流动电荷】,且给定下列二者与边界有关信息之一:

  • 边界电势 φ|S
  • 边界法向电势导数(法向电场强度)φn|S

则内部电场即可唯一确定。

引理:由 Gauss 定理,

(C)dV=CdS

现在,若满足 C=Φφ,其中 Φ,φ 是任意标量场,则

(C)dV=CdS(Φ2φ+Φφ)dV=(Φφ)dS=ΦφndS

上式称作 Green 定理。

因为 E=ρε0,而 E=φ,因此有 2φ=ρε0 的 Poisson 方程。

现在,假设存在两个满足唯一性定理条件的分布 φ,φ,则在内部均满足 Poisson 方程。

U=φφ,则在 V 内部有 2U=0 的条件;边界上条件为:

  • 电势处处为零。
  • 电势法向导数为零

二者之一。

现在,在 Green 定理中,令 φ=Φ=U,则

(U2U+UU)dV=UUndS

因为或者 U=0,或者 Un=0,所以 RHS 为零。

而因为在内部 2U 处处为零,所以 (U)2dV=0。于是可知 U 处处为零,即其在空间内部为常数。

若给定的边界条件是边界电势,则易知 U=0,可知解唯一;若给定的边界条件是边界梯度,则解之间相差常数。


现在,扩充唯一性定理至有导体的场合。

则:

  • 在无导体的空间上,电荷分布给定。
  • 在每个导体的电势/电荷给定。
  • 边界的电势/电势梯度给定。

此时电场唯一分布。


Image charge 法:

在导体内部幻想一些电荷,保持导体的边界条件不变。于是可以用幻想出的电荷来替代导体中的诱导电荷。

例:无限接地面。对一侧镜像到另一侧,则接地面上电场处处为零。

例:接地球。对接地球外一点,其镜像点即为与其互动形成阿氏球恰为该接地球的另一点。

  • 如果是非接地球呢?在球心放一个与幻想电荷电量相反之电荷即可。

Image Charge 法不一定能用。更常见的场合,在分析导体内部的电场时,还是只能处理如下的问题,即有边界条件的 Laplace 方程

2φ=0

三维的场合,常常没有解析解;但是,如果是二维或近似二维,即电场在 z 轴方向没有或几乎没有变化的场合,则 Laplace 方程退化为 Poisson 方程

2xx2+2yy2=0

此时可以使用复变函数方法解决。

对于解析函数,其满足额外的 C-R 方程条件:

ux=vyuy=vx

再求导可知,解析函数的实部与虚部均满足 Laplace 方程,即

uxx+uyy=0vxx+vyy=0

  • 除此之外,还有额外条件 uv=0,即 u,v 的等势线处处正交。

则,如果可以将 φ 扩充为一个解析函数的实部或虚部之一,则其即满足条件。

而二维的场合,取一个为等势线,一个为电场线,则其满足要求!

首先,有

E=V=VxiVy

然后,

W=BAEdl=BA(Vx+iVy)(dx+idy)=BAVxdx+Vydy=V(B)V(A)

其中, 被看作是复数之间的点乘,实部虚部分别相乘;dl 是路径 AB 上的长度微元。

而令

N=BAEds=BA(Vx+iVy)(dy+idx)=BA(VxdyVydx)=:BAUydy+Uxdx=U(B)U(A)

这里的 ds 是法向量与路径垂直的面积微元的法向量;U 是一个人为定义的通量函数(类似电流密度函数的概念);沿电场线走,通量 U 处处不变。

这么搞之后,有

E=VxiVy=VxiUx=(i)ˉf(z)

f(z)=(i)ˉE


假如我们已知 f(z),则其实部即为电通量 U,虚部即为电势 VU 等势即为电场线,V 等势即为电势线;f(z) 被称为 复势

带电长直导线的电场为

E=λ2πε0r2r

用复数表示就是

λ2πε0ˉz

于是知

f(z)=(i)λ2πε0z

积分即知

f(z)=(i)λ2πε0lnz

这即为长直导线复势公式。


在所有 xZ,z=0 处,放置着无限长垂直 y 轴导线。则在 x,z 两维上满足 Laplace 方程

2φx2+2φz2=0

通过 Fourier 展开,猜想每一项是

φn=Fn(z)cos2πnxa

的形式。代入得到应有

d2Fn(z)dz24π2n2a2Fn(z)=0Fn(z)=Anez/z0,z0=a/2πn

因此,φn 关于 z0 为比例快速衰减。同时,有 F0(z)=E0z+C,于是有总式子

φ=E0z+Anez/z0(n)cosxz0+C

zz0(n) 时,有 E=E0。这说明,数根铁棒编织成的网络就足以起到静电屏蔽之作用,密封性不用卡太死。


电场中存储着电势能。

电势能

U=12φdq=12φρdV

而由 Gauss 定理,2φ=ρ/ε0,因此有

U=ε02φ2φdV=ε02((φφ)φφ)dV=ε02E2dVε02(φφ)dS

注意,U 的积分是全空间的,因此后一项的 S 应该膨胀趋于无穷;然后因为电荷分布是局域的而非全空间的,在充分远处有 φ1r,而 φ1r2,但 Sr2,因此 S 时最后的积分趋于 0。因此即有

U=12ε0E2dV

其中 12ε0E2 即为电磁场的局域能量。

  • 这个模型对于点电荷的场合无效:计算可知理想点电荷的能量趋于无穷。因此存在经典的电子半径 rc2.82×1015m,小于该半径则电荷无法被当成点电荷。
  • 因此,当任意形式下牵涉到点电荷时,再用全空间电场积分就不对了!
  • 当同时牵扯到点电荷和连续电荷时,我们唯一能使用的公式就是 12(φ1dq+φ2q);其中,φ1 是空间中该点的电场,φ2 是空间中除 q 外其它物体生成的电场。

忽略边缘效应,则两极板间的电容为

C=QV=Aε0d

其中 A 是极板面积、d 是距离。

现在往中间插一块厚度为 b 的导体板,易知电容变为原本的 ddb 倍。

现在往中间插一块绝缘体板。在靠近正电荷的一侧,感应出少量负电荷(但不如导体板的场合那么多);在靠近负电荷的一侧同理。在绝缘体板与电容板间,电场不变;在绝缘体板内部,电场变小,但没有像导体版的场合一样衰减至零。因此,插入绝缘体板后,电势差变小,但没有像导体的场合一样小,因此电容变大,但不会到 ddb 的限度。


电极化强度 P=Nqδ。因此,P 与电荷面密度 σ 有相同量纲,并且有 σpol=Pen,其中 σpol 是极化电荷面密度。

PdS=qin=qout

并且

EdS=qin+qinε0

因此

(εE+P)dS=qin

DdS=qin

其中 D=εE+P,和面电荷有相同量纲。当 P=χε0E 时,有 D=(1+χ)ε0E=κε0E

极化强度的 Gauss 定理:ρpol=P

越能导的物体,介电常数越大;金属的常数是无穷,因为它可以把电子移到无穷远,因此 P=


现在,提出电介质问题的方法论。

PEE0 互相影响。这不牛!

但是 D 只与自由电荷 (free charge) 有关,因此可以列 D 的 Gauss 定律等算 D,然后由 D=εE 的关系算得 E,然后由 E 即可知 D

例:宽为 d 的无限大电介质板,相对介电系数 εr,内部均匀分布电荷 ρ0

取圆筒,筒的两个底面均在介质板外,则有

2DS=ρ0dSD=E=ρ0d2

取圆筒,筒的两个底面均在介质板内,则有

2DS=ρ02xSD=ρ0xE=ρ0xεrε0P=(εr1)ρ0xεrε0

其中 x 是圆筒距介质板中心的距离。

可以按照 ρ=P 的方式计算极化电荷。在这个问题中,介质板表面分布着均匀的正极化电荷,内部分布着均匀的负极化电荷;实质是正的 ρ0 电荷密度周围吸附了均匀的负极化电荷。

介质中的点电荷会被极化电荷屏蔽,电场除以 εr


两导体极板带 ±Q 的电荷,面积为 S,距离为 d,内部是 ε0 的真空。现在往一半的面积中充入 εrε0 的电介质,各数据如何变化?

首先,因为两极板都是导体,所以充入电介质后其仍然等势,因此有无电介质不影响极板间的电场处处匀强。因此 ±Q 的电荷会以 1:εr 的比率分配给 ε0 周围和 εrε0 的周围,把电介质中的电场强行撑大。

  • 模型其实也可以被看作是两电容并联,两端电压相等,总电量会在二者间分布。有 C=C1+C2。因此 C1=ε0S/2d,C2=εrε0S/2d
  • 原理是,电介质抵消了电场,但为了保持匀强电场所以电子不得不移动过来。

电流密度 j 是单位时间通过单位面积的电荷量,即有 ΔQ=jΔsΔtj 的单位是 Cm2s1,或是 Am2

v 为载流子速度,则有 j=ρv=Nqv,其中 N 为载流子密度。I 是电流密度通量,即

I=Sjds

电流守恒定律:

Sjds=˙Q

令体积趋于零,可知

j=˙ρ

由 Lorentz 公式 F=qv×B 可知,dF=NdV(qv×B)=j×BdV。现在令讨论区域是导线,则 dV=Ad,于是

dF=I×BddFd=I×B

其中,I 是电流矢量,被认为垂直于积分面。此乃 Ampere 力公式。


静磁学:

  • ρ 恒定。
  • 电流是恒定电流。(这意味着 j 恒定)
  • 电场、磁场不随时间改变。

磁场源为电流密度,即 ×B=μ0j

Ampere 定律:

Bdl=(×B)da=μ0jda=μ0I

长螺线管内部磁场比外界强很多。对螺线管附近取一条回路,使用环路定理可知,B0=μ0nI,其中 n=N/LN 是总匝数,L 是总长度。


F=qv×B,j=Nqv。我们的 v 是相对于哪个系的速度?

答曰:取任何惯性系均可。

但是不同惯性系下理应遵从同一物理规律。这暗示着电磁场是统一、可以相互转换的。

导线中有电流 I 在从右往左流。负电荷以 v0 的速度在导线上方从左往右跑。

在导线上方的磁场是垂直纸面向里的磁场。因为是负电荷,电荷会受指向导线的力。

现在,在负电荷静止、导线运动的系中分析。这个系中,电荷静止,因此不能受到磁场力;受到的是电场力,因为运动的导体不再处处电中性,会有电场。

具体推导?不妨令导线中载流子即为电子。在导线静止系中,ρ+ 静止,ρv 的速度从左往右跑。只要满足处处 ρ+=ρ,则导线即为电中性(注意,这里的 ρ 已经经过相对论修正)。有

F=qv0×BF=qv0μ0I2πr=qv0μ0ρvA2πr

B=μ0I2πr 由 Ampere 定理保证)其中 I=ρvA

现在考虑 v0=v 的特殊状况。因为满足 c2=1μ0ε0 的好性质,所以有

F=qρA2πrε0v2c2

S 坐标系下,ρ 是静止的,ρ+v 的速度移动;ρ+ 会生成磁场,但是没有磁场力。

现在,问题在于由相对论效应,运动的长度会变短,因此运动的 ρ+ 会比静止的 ρ+ 更大,有 ρ+=ρ+1v2/c2;同理,ρ 是静止坐标系,因此运动的 ρ 会比静止的 ρ 更大,因此 ρ=ρ1v2/c2

因此,总 ρ=ρv2/c21v2/c2,乃是一带正电导线;计算

F=qρ+A2πε0rv2/c21v2/c2F=F1v2/c2

这个区别在于,力在 Lorentz 变换下也会得到改变,改变关系恰为上述关系。


电荷密度和电流密度共同组成了一个四维矢量

ρ=ρ01v2/c2j=ρv=ρ0v1v2/c2jμ=(j,icρ)xμ=(x,ict)

jμxt=0j+ρt=0


磁感应强度 B 有对应的磁矢势 A,满足 B=×A。因为梯度的旋度永远为零,所以可以有 A=A+φ

Coulomb 规范 (Coulomb Gauge) 要求 A=0。虽然此时 A 仍然不能唯一确定。例如,令 B=B0ez,则 A1=xB0ey,A2=yB0ex 以及它们的 {t,1t} 线性组合都是满足 Coulomb 规范的 A

因为 ×B=μ0j,因此有 ×(×A)=μ0j;因此 (A)2A=μ0j;而由 Coulomb 规范有 A=0,因此有 2A=μ0j

注意:对于矢量场的 Laplace 方程,其等效于每一维分别的 2Ai=μ0ji,其中 i{x,y,z}。于是,即有每一维上分开来类似于静电场的 Poisson 方程,每个 ji 等效于一个 ρ(事实上,它们之间的区别是把一个长度量纲换成了时间量纲)。

在全空间上,静电场 Poisson 方程的解即为

φ(r)=14πε0ρ(r)dVrr

因为 A 的每一维均遵循和 φ 仅有常系数区别的方程,所以可类推得到

A(r)=μ04πj(r)dVrr

因此,由电流分布求 B 的方式,也可以遵照上式,先由 jA,再由 B=×AB

既然 jρ 相似,那么 I 就与 λ 相似;面电流 J 与面电荷密度 σ​ 则相似。


例如,考虑沿 z 轴正方向一股电流 I,有

j=Iπa2ez

jx=jy=0,而 jz 就像一根 ρIπa2 的带电棒。有 φ=λ2πε0lnr,则 Az=μ0I2πlnr。求旋度即得,B=μ0I2πr×er,与通过 Ampere 环路定理推出结果相同。


同理,考虑长螺线管 solenoid。其只在表面有电流密度 j,则 j 在表面绕圈,有 jx=jsinθ,jy=jcosθ,jz=0。表面的 j 就可以类比表面的 σ,即 σsinθ 的圆柱,可以看作两根 ±ρ 的圆柱重叠。推一坨即可。


平面上有一个正方形逆时针电流元。则 jz=0,于是 Az=0;在 x,y 方向上,它各自等效于两段相对电流。这其实类似于电偶极子。

具体而言,有

φ=14πε0peRR2Ax=μ04πIabR2yR

(注意,负号是因为逆时针等效的电偶极子是沿 y 轴负方向的)

总结为,

Ax=μ04πIabR2yR,Ay=μ04πIabR2xRA=μ04πIabR2ez×er

因此,对于面 S,其磁偶极矩 μ=IS,其中 S=SeS。则有

A=μ04πμ×erR2=μ04πμ×RR3Bx=Ayz=μ0μ4π3xzR5By=μ0μ4π3yzR5Bz=μ0μ4π(3z2R51R3)B=μ04π(μR3+3(μR)RR5)

E=14πε0(pR3+3(pR)RR5) 类似。即,环电流元与依右手螺旋定则穿过其之电偶极子类似。

虽然在源周围,电场和磁场遵循不同的原理;但是在远场处,二者的散度和旋度均可以近似为零。这就是为何二者有相似的形式。


对式子

A(r)=μ04πj(r)dVrr

求旋度可以得到 B(r)。有

B=r×A=μ04πr×j(r)rrdV=μ04π(r1rr×j(r)+r×j(r)rr)dV

这里,注明 r 意为关于场点坐标 r 求梯度,而非源点坐标 r。因此,r×j(r)=0

=μ04π(errrr2×j(r))dV=μ04πj×e12r212dV

这个式子与电场中

E=14πε0ρe12dVr212

的式子类似。

并且,因为 jdV=Id,因此即有

B=μ04πId×e12r212

此乃 Biot-Savart 定理。


导线附近的磁场是

B=μ04πIr0(cosθ1cosθ2)

因此,如果是无限长导线,则其为

B=μ02πIr0

与 Ampere 定理结果吻合。


在均匀磁场 B 中的磁偶极子受力。虽然合力为零,但其受力矩 τ=μ×B,与电场中电偶极子受力矩 p×E 相类。由虚功原理可知,U=μB。但是注意,这并非一个电流回路所具有的全体能量:例如,维持一个电流回路需要某些能量。因此,我们只把这部分能量称作 Umech,因为它其实是一种势能(机械能)。即有 Umech=μB。但是,这个能量与 pE 不同,它不是真实的能量(前者是把两个电荷按在一块所需要的能量)。

对于一个矩形回路而言,Umech 同时也上把它从无磁场区域移至磁场中时,非电磁力所做功。同时,在非均匀磁场中,电磁力非零,且有 F=(μB)

Umech 虽然并非总能量,但它是计算受力等场合一个好用的原理,前提是所有电流都是恒定电流。

通过对 Lorentz 力等分析,我们可知:磁场中的导线受外力移动。对导线所做机械功总是恰等于对电流源所做电功。这是因为,F=q(E+v×B);如果没有 E,则 Fv=0。一个变磁场会产生电场,因此此分析只适用于恒定磁场中移动的导线。

虚功原理缺了什么呢?缺了正在产生磁场的电流能量。

现在假设有两个线圈。把两个线圈从无穷远处摆在一块。

固定线圈 A,把 B 移过来。B 移过来的过程中,其电流在运动,磁场会变动,于是产生电场,其对 A 中的电荷会做功,这个能量未被考虑。这一能量记作 UB,则 UB+Umech=0。同理,如果看作是固定 B、移来 A,有 UA+Umech=0

总能量 Ut=UA+UB+Umech=Umech。因此,有磁偶极子的真实能量 Utotal=Umech=(μB)=μB。只有在假定一切电流都是恒定电流时,才可以应用 Umech 求机械力。

这一部分可以类比:首先有电容器能量为 12Q2C。应用虚功,可知 ΔU=Q22ΔCC2

但是,假如我们人为定义一个值为 12CV2 的能量(它与真实能量相反),声称其变化等于机械功,并且忽略维持电压不变由电源做的功,也可正确获得结果。

上述推导类比的结果是,把磁偶极子扔进磁场,忽视磁偶极子对磁场的影响,则能量为 μB;但是,磁偶极子会对磁场源有一些影响,这些影响又需要额外的其它东西来抵消;考虑到这部分影响后,总能量即为 μB​。


对于任意形状的单一电流回路,其总能量

U=μB=IBΔS=I平面BdS=I边界Ad

现在,复杂回路可以看作众多上述单一回路的总和。于是因 Id=jdV 即有

U=12全空间jAdV

与静电学中 U=12ρφdV 类似。


在电动力学中,电磁场的传播需要时间。因此,φA 的公式改变为延迟势(retarded potential)

φ(x,t)=14πε0ρ(x,tr/c)dVrA(x,t)=μ04πj(x,tr/c)dVr

影响会以 c 的速度传播。同时,B,E 的关系也略有修改

B=×AE=φ˙A

其中,˙A=At 的项即来自 Faraday 电磁感应公式。

但是不考这个捏。


一段回路的电动势(electromotive force)是单位电荷绕其一圈运动时,力(不分种类)做功的量。即,

qE=Fd

这里为什么可以是不分种类的力?因为静电场环量为零。

均匀磁场上两根平行导轨电路间滑动杆子产生的电动势 E=Blv

  • ……所以,你的量纲根本不是力的量纲而是电场的量纲,为什么要叫 em force 呵喂!

通量法则:环路上的 E 是环路内部磁通量变化率的相反数。即,

E=˙Φ

在杆子滑动的场合,电动势是 Lorentz 力分力产生的。但是,在只有磁场本身变化时,不会有 Lorentz 力,因此电动势只有可能由电场力提供,即电场本身有电势。

于是有 Faraday 定律:

×E=˙B

验证其结果,即为

Ed=(×E)dS=˙BdS=tBdS=˙Φ

通量定律 E=˙Φ 总是正确的。考虑电磁力公式 F=q(E+v×B),其中环路中的 E 如果总功非零,则必然来自于变化的磁场,v 则是导线切割磁感线的影响。


在匀强磁场环境中放置一个环形轨道。

2πrE=ddt(Bavgπr2)E=r2˙Bavgdp=qr2dBavgΔp=qr2ΔBavg

同时,圆周运动会有向心力,即

qvBorbit=dptdtdptdt=ωp=vprqBorbit=prΔp=qrΔBorbitΔBavg=2ΔBorbit

如果满足这样的关系,则不需要轨道就能约束电子。


在圆盘的圆周上固定一堆带电金属球,圆盘的圆心周围通过一个磁场。现在,撤去磁场,按照 Faraday 电磁感应定律,在圆周上应该感应出电动势,于是金属球受力,圆周应该转动。

一个东西凭空开转了,是否违背角动量守恒?

答曰:并非违背,此乃电磁场角动量。


匀强磁场中有一导线环。其通量 Φ=NBScosθ。若它以 ω 开转,则 E=NBSωsinωt。此乃交流发电机。

在远离发电机的部分,电磁场几乎不再变化,此时回到静电场的场合,可以定义电势 V,且满足 V=E

其功率

dP=vF(nd)P=vFnd=vnFd=EI=NISBωsinωt


两根嵌套螺线管。B=μ0N1I1/,则 E2=N2S˙B=μ0N1N2S˙I1=M21˙I1,其中互感系数 M21=μ0N1N2S。由对称性,两互感系数相同。

事实上,任两个线圈间都有互感。

E1=ddt1BdS=ddt1AdA=μ04π2I2dr12E1=μ04πddt12I2d1d2r12=M12˙I2

事实上,一根线圈也会有自感。自感必然为负。

总的

E1=M12˙I2+M11˙I1E2=M21˙I1+M22˙I2

其中 M11=L1,M22=L2

单个线圈的

E=L˙I

即,自感总是遏制电流的变化。

综上,有或许没啥用的

M=μ04π12d1d2r12

的公式。

对比:

V=L˙IF=m˙vU=12LI2W=12mv2

其中,第二行是由

˙W=EI

积分得到 W=U=12LI2 得到的。这是一个电路的自能。

电路系统总能量

U=12L1I21+12L2I22+MI1I2

考虑建立电路的过程:I1 先被建立,然后建立 I2 时,为了遏制 I1 中的互感,需要额外付出 MI1I2 的能量。


自感系数如何计算?考虑 U=12LI2U=12jAdV 两个公式,得到

L=1I2jAdV

与电流无关。

U=12μ0(×B)AdV=12μ0((B×A)+×AB)dV=12μ0(B×A)dS+12μ0B2dV

其中,有 B1r2,A1r,因此若电流是局域的,让被积分区域遍及全空间,则前一项积分为零。于是有

U=12μ0全空间B2dV

这个公式与静电场的 ε02E2dV 类似。


Ampere 定律是 ×B=μ0j,推论为 j=0,这个定律和电荷守恒定律,即 j+˙ρ=0 不协调。

因此,Ampere 定律必须得到修正,才能与电荷守恒定律适谐。

相反,在考虑“位移电流”后,其修正为

×B=μ0j+μ0ε0˙E

其求散度后则与电荷守恒定律适谐。

如何理解位移电流?考虑一个不断向外辐射电子的电荷元,则其相当于有向各个方向的电流,则其会有磁场。可知:所有电流的磁场彼此抵消,全空间内没有磁场。

˙Q=4πr2j,同时电场也会变化,使用 Gauss 定律分析知 ˙E=jε0

因为全空间没有磁场,所以须有 ×B=μ0j+ε0μ0˙E=0


考虑一个正在靠无穷远处电源充电的电容器。则有 I=˙Q

现在,选取 +Q 盘上方包含 I 的一个回路 γ,以其为边界选取面 Γ,则若导体是理想导体,则电流中不应有电场,直接即有

γBd=μ0Γjds=μ0I

但是,如果换成通过电容区域的 Γ,此时 μ0I 不再通过 Γ,这一部分理应由 ˙E 所弥补。有

ε0μ0˙EdS=ε0μ0ddtQε0=μ0I

给出相同的结果。

即,当空间中的电场 E 改变时,应视作有一个 ε0˙E 的等效 j 流过。

对比:

μ0ΓjdS=μ0I=μ0ε0Γ˙EdS


y-z 平面上有一个平面。现在,突然加之一个面电流 J 沿 y 轴正方向。会有什么影响?

假如 J 是恒定电流,那么会产生 B=μ0J/2,方向在 x 轴正方向侧沿 z 轴负方向,负方向侧沿正方向。

显然,这个磁场不可能瞬间充满全空间。场仅能以有限速率传播。在其波前传播到某一位置时,其左侧有 μ0J/2 的磁场,右侧有 0 的磁场。在波前处,会有 EL=BLv,即 E=Bv 的感应电场。在 x 轴正方向上有磁场的区域,其同时会有方向沿 y 轴负方向的电场。

随着波的传播,不仅是磁场在 build up 并感应出电场,这个电场同样会 build up 并感应出磁场。有 BL=1c2ELv,则 B=1c2Ev。这个磁场应该恰为感应出电场的磁场,由此可得须有 v=c

综上,在打开电流 J 指向 y+ 时,沿 x+ 方向会有:

  • B=μ0J/2 沿 z 方向。
  • E=Bv 沿 y 方向。

的电磁场传播。

同理,倘若电流在 T 时刻被切断,则相当于又叠加了一个 J,会出现 B,E 的电磁场在传播。最终,会有一个长为 vT 的波包向外界传播,这个波包即不再依赖波源。

一切 TEM (Transverse Electric Magnetic) 波都满足:

  • E=Bc
  • ε0E22=B22μ0
  • 波传播方向 k、电场方向 B、磁场方向 E 三者互相垂直。

TE 波只有电场垂直于传播方向,而 TM 波只有磁场。


B=0 恒成立,可知总是可以通过磁矢势 A 来表示 B=×A

因此,由 ×E=˙B,可知 ×(E+˙A)=0,因此可以定义 E+˙A=φ,则有 E=φ˙A

A=A+ϕ 时,这个操作不改变 B 但是改变了 E;当 φ=φ˙ϕ 时,B,E 均不会改变。

另外两个 Maxwell 方程用 φ,A 重写后是

2φ˙A=ρε0c22A+c2(A)+˙φ+¨A=jε0

使用 Lorentz 规范 A=1c2˙φ,其即化简为

2A1c2¨A=μ0j2φ1c2¨φ=ρ/ε0

通过 Lorentz 规范,Aφ 即解耦,并且二者即具有相似模式。通过 ρ,j 解出 φ,A 再分别得到 E,B 往往比直接解方程要容易。

在远场环境下,ρ,j 都为零,此时二者具有相似的波动方程模式

2ψ1c2¨ψ=0

此时可以验证,BE 二者亦满足相同的波动方程模式。即,远场时有

2B1c2¨B=02E1c2¨E=0

因此,远场时不需要依赖 φ,A,可以直接解 B,E


Maxwell 方程有八个方程(有两个是矢量方程)。为方便,我们假定 B,E 均只与 x,t 相关,且仅考虑远场场合,此时

E=0Exx=0B=0Bxx=0(×E)x=0Bxt=0(×B)x=01c2Ext=0

这表明,此时 E,B 均须是匀强的。这是平面波的场合。

现在,假设 E 只有 Ey 分量、Ez=0

(×E)y=ExzEzx=Byt=0(×E)z=EyxExy=Bzt(×B)z=ByxBxy=1c2Ezt=0(×B)y=BxzBzx=1c2EytBzt=Eyx1c2Eyt=Bzx

综上,只有 Ey 分量的 E 激起只有 Bz 方向的 B,只有 Ez 方向的 E 激起只有 By 方向的 B

对上述方程整理得到满足

2Eyx21c2¨Ey=02Bzx21c2¨Bz=0

其解即为 ψ=f(xct) 的通解。也即,其相当于一个波动 fc 的速度沿 x 轴正方向传播。

综上,通解为:

E=(0,Ey,Ez)B=(0,By,Bz)Ey=f(xct)+g(x+ct)Ez=F(xct)+G(x+ct)cBy=F(xct)+G(x+ct)cBz=f(xct)g(x+ct)

传播方向(能量流方向)k 沿着 E×B 的方向(这是所谓 Poynting vector)


一切电磁波都可以被描述为平面波的叠加。但是,一些场合使用球面波会更好!

此时,满足 ψ=ψ(r,t)。处理可得,

2ψ=ψ+2rψ=1r(rψ)

(其中,此处 意为对 r 求导)

于是,波动方程 2ψ1c2¨ψ=0 退化为

(rψ)=1c22t2(rψ)

解为 rψ=f(rct)。即,ψ=f(rct)r,振幅会随着距离反比衰减。因为能量正比于振幅平方,波面正比于距离平方,所以这个关系保证能量守恒。

  • 注:其有另一解 ψ=g(r+ct)r。但是这另一解对应的波是向中心传播的内向波,不符合物理实际,应舍去。

f(rct)rr0 时发散。但是这是合理的:在源点处 Maxwell 方程不再退化为波动方程。

r0 时,方程解退化为 ψ(r)=f(t)r 的形式。同时,1c2¨ψ 一项较之 2ψ 也可以忽略(r0 时,小尺度下,空间变化率相对于时间变化率是极大的)。考虑源后,方程退化为 2ψ=s

这个方程类似于 2φ=ρ/ε0 的静电学式子。类比知

ψ=s(r,tr/c)dV4πr

这个势是 retarded potential(推迟势)。


交流电路!如果 c/f=λl,其中 l 是元件尺度,则此时原件可以被看作是理想的 lumped elements(集总元件)。而短波电路的场合,元件则应被看作 distributed elements(分布元件)

time harmonic(时谐)的量满足 sinusoidal(正弦)关系,例如

V=V0exp(iωt),I=I0exp(iωt),E=E0exp(iωt)

Kirchhoff 定理:环路的 Ed=Vn=0;节点的 In=0

一段不含电源的电路总是可以被当作一个阻抗。任何电源都可以被当成一个理想电源和一个阻抗。

Norton 定理:等效电流源时,电流是两接口短路的电流,阻抗是把所有电流源断路、电压源当成短路的阻抗。

Thevenin 定理:等效电压源时,电压是两接口断路的电流,阻抗是把所有电流源断路、电压源短路的阻抗。

阻抗 Z=R+iX,因此可以被看作一个纯电阻 R 和一个电抗 reactence X 的串联。则有

E=I(R+iX)=I0Rexp(iωt)+iI0Xexp(iωt)ReE=I0RcosωtI0XsinωtReI=I0cosωtˉP=1TT0(Reω)(ReI)dt=12I20R

因此,只有实部 R​ 会耗能,虚部 X​ 不耗能。


梯型网络

-z_1-+-z_1-+-z_1-+-...
     |     |     |
    z_2   z_2   z_2
     |     |     |
-----+-----+-----+-...

把整个网络看成阻抗 z0,则有

z0=z1+z2z0z2+z0

解得 z0=z12+z21/4+z1z2

  • 注意,不能取减号,否则会出现负电阻。

现在,考虑 z1 是感抗 iωLz2 是容抗 iωC 的场合。此时 z0=iωL2+L/CL2ω2/4

发现,如果 ω>2/LC,则 z0 是纯虚数,不耗能;小于的场合,会耗能,这是因为能量会向着远处传播,传播系数为

α=L/Cω2L2/4iωL/2L/Cω2L2/4+iωL/2

实质是,ω<2/LC 传播系数为 1,传播不会衰减;大于时系数小于 1,其会衰减,迅速衰减至 0

这可以被当作一个低频滤波器,只允许低频波传输,高频波会被反射(当成短路)。

同理,如果把 LC 换位,则其变为高频滤波器,只允许高频波传输,低频波会被当成断路。

导线单位长度也会有其电容和电感。其即可被梯形网络建模。分析可得,满足

2xVL0C0¨V=02xIL0C0¨I=0

的波动方程,相速度 v=1/L0C0,波长 λ=2π/ωL0C0,特征阻抗 z0=iωL2+L/Cω2L2/4L0/C0,且滤波频率是正无穷(也即不会滤波)

即,导线中传输的信号,阻抗与相速度均只与导线自身相关,并且可以传输一切频率的波。


交流电路中的电容,其极板间的电场是交变的,其同样会 induce 一个交变磁场,这个交变磁场又会 induce 交变电场……

E1=E0exp(iωt),则 2πrB1=1c2πr2˙E1,于是 B1=iωr2c2exp(iωt)

又有 hE2=r0ω2R2c2exp(iωt)hdR,解得 E2(r)=r2ω24c2E0exp(iωt),其与 E1 方向是相反的。二者叠加后,E 场即不再是匀强电场。

……但是,这个 E 又会感应出另一个磁场……

最终发现,解可以写成一个 Bessel 函数

J0(x)=i=0(1)i(i!)2(x/2)2i

的形式,满足 E=E0exp(iωt)J0(ωr/c)

这是柱坐标系下 Laplace 方程的一个解。

x2.405 时,J0(x)=0。因此,在一个圆柱体满足 ωr/c=2.405 时,其内部的电磁场可以与外界隔绝,此乃 resonator 谐振器……

……但是,在边界上,虽然电场等于零,但是磁场却不等于 0。这意味着谐振器的侧面上会有电流(电容中,电流会从外电路走;但是谐振器是一个导体圆筒,所以电流会从筒壁走),而上下底面上有电荷。

在没有电阻时,谐振器可以一直振。但是如果考虑电阻,能量则会耗散。当频率接近某值时,震荡效果最好。

……并且,Bessel 函数的根不止 2.405 附近的那个,它有多个根。不仅如此,J0 仅仅是一种 Bessel 方程,还有多种 Jn(x) 的 Bessel 函数系。因此,谐振器的模式可以有很多种。


真空中 Maxwell 方程是

E=ρ/ε0B=0×E=˙B×B=μ0j+ε0μ0˙E

现在考虑在介质面,面两侧的电导率分别是 ε1,ε2,磁导率均为 μ0,则电磁场如何变化?

ΔS×Δd 的界面两侧的小体。则 ΔS(E2E1)+SΔd(E2E1)=Qε0

Δd0。可知:en(E2E1)=σ/ε0。而,有 en(B2B1)=0

注意,此处的 σ 是总电荷密度 σtot。假如换成 D 的语言的话,就可以换成自由电荷密度

en(D2D1)=σfree

Δl×Δd 的界面两侧的小环。则 Δl(E2E1)+Δd(E2+E1)=ΔlΔd˙B.

Δd0,可知 E2=E1。也即,en×(E2E1)=0。同理对磁场分析可知,Δl(B2E1)=μ0I,由此有 en×(B2B1)=μ0J,其中 J 是面电流密度。

综上,界面两侧的 Maxwell 方程为

en(E2E1)=σtot/ε0en(B2B1)=0en×(E2E1)=0en×(B2B1)=μ0J

因此,电场的突变只能发生在界面的法向(例如导体法向的电场会突变到零,这是面电荷导致的),切向永远不能发生电场突变;磁场的突变只能发生在界面的切向(这是面电流导致的),法线则不能发生磁场突变。

电场强度 E 有对应的电位移 D=εE,而磁感应强度 B 也有对应的磁场强度 H 满足 B=μH


在远场场合,电磁场满足方程

×E=˙B×B=μ0ε0˙E

×(×E)=t(×B)=μ0ε0¨E

因为 E=0,分析旋度性质可知

2E=1c2¨E

磁场同理有

2B=1c2¨B

若我们认为有 E(x,y,z,t)=E(x,y,z)exp(iωt) 的时谐波(即,电场的时空相对分离)的话,则有

2E(x,y,z)+ω2c2E(x,y,z)=0

因此,若令波数 wave number k=ωc(亦有 k=2πλ),则得到 Helmholtz 方程

2E+k2E=0E=0


在柱坐标系 (ρ,φ,z) 下,Helmholtz 方程同样成立。此时,有

2=2ρ2+1ρρ+1ρ22φ2+2z2

若认为其在 z 方向亦是谐波 Ez=ez(ρ,φ)exp(jβz)=R(ρ)P(φ)exp(jβz),代入 Helmholtz 并推了依托知

ρ2R+ρR+(k2β2)ρ2RR=PP

因为在柱坐标系下,P 需要是 2π 的整周期函数,所以知上式必须等于 n2

kc=k2β2,一通乱推后值得

R(kcρ)=AJn(kcρ)

其中 Jn 是第一类 Bessel 级数。

Bessel 级数(第一类 Jα,第二类 Yα)是方程

x2y+xy+(x2α2)y=0

的解,满足

y=AJα(x)+BYα(x)

其中,Jα(0) 有定义,而 Yα(x0)

因此综上,有

Ez(ρ,φ,z)=AJn(kcρ)(c1sinnφ+c2cosnφ)exp(jβz)

其中,k 的效果即在于影响 kc 进而影响 Bessel 函数的零点位置。零点是重要的,因为在导体筒的边界条件上,Ez 是切线分量,而导体筒内部是无电场的,所以边界条件须是 Ez=0


考虑 a×b×c 导体盒中的电场。其仍然满足 Helmholtz 方程,因此含时的部分已经被消去。

与柱的场合不同,此时 x,y,z 三维是等地位的。对于某一维 i{x,y,z},假设其三维独立,即 Ei=Xi(x)Yi(y)Zi(z),代入可知

XYZ+XYZ+XYZ+k2XYZ=0X/X+Y/Y+Z/Z+k2=0

这需要 X/X=k2x,Y/Y=k2y,Z/Z=k2z。为什么是负号?因为各个的方程 X+k2xX=0,其解为 X=Cxcoskxx+Dxsinkxx,而如果取成负号的时候,它的解就会是指数形式而非三角函数形式,在导体盒中就不合理了。

因此,总结,则解为

{X=Cxcoskxx+DxsinkxxY=Cycoskyy+DysinkyyZ=Czcoskzz+Dzsinkzzk2x+k2y+k2z=k2Ei=XiYiZi

不妨认为我们在考虑 Ex。此时在 y=0/b,z=0/c 两处,Ex 会是无法突变的切分量,因此需要 Y(0)=Y(b)=0,则要求对于某个整数 m,需要 ky=mπb,且 Cy=0kz=nπc,且 Cz=0

同时,又因为 E=0,则有 Ett+Enn=0;但是边界处有 Et=0,所以则有 Enn=0,即在 x=0/a 处,应有 Exx=0。因此,需要 X(0)=X(a)=0,即 kx=tπa,且 Dx=0

综上,x,y,z 彼此分离的解,就是

{Ex(z,y,z)=Acoskxxsinkyysinkzzkx=mπaky=nπbkz=pπck2x+k2y+k2z=k2

这是以固定频率在铁盒子里振的方式。因为 kx,ky,kz 彼此能取的值有限,因此 k 也只能有限,所以含时的频率 ω 也只能取若干离散的值。因此,铁盒里的频率是离散的,而非连续的。

考虑 x,y,z 三维,则有

Ex=AxcoskxxxsinkxyysinkxzzEy=AysinkyxxcoskyyysinkyzzEz=Azsinkzxxsinkzyycoskzzz

但是,由于 E=0 的强性质,这相当于有无穷多个方程,这要求我们必须有 kxx=kyx=kzxy,z 同理。因此综合而言就直接有

Ex=AxcoskxxsinkyysinkzzEy=AysinkxxcoskyysinkzzEz=Azsinkxxsinkyycoskzz

现在,求散度后,知须有 Axkx+Ayky+Azkz=0。而因为 kx=mπa,ky=nπb,kz=pπc,可知 m,n,p 中至多只有一个为零,否则电场直接就不存在了。

因为 k=ω/c,所以即得时间频率

ω2mnp=(m2π2/a2+n2π2/b2+p2π2/c2)c2

不妨认为 a,c>b,则 ω101=π2/a2+π2/c2c 是基模(fundamental mode),其为该铁盒中能谐振的最低频波,低于该频的波都不可能谐振。该波为

Ex=Ez=0Ey=Aysinxaπsinzcπ

因此,基波总是沿最短边的波,其在正中心的振幅最大。

×E=˙B=iωB。因此

B=iω×E={Bx=iωAyπcsinxaπcoszcπBy=0Bz=iωAyπacosxaπsinzcπ

相似的结果也可以通过解 B 的 Helmholtz 方程,只不过满足与 E 不同的边界条件,来解出。

在矩形谐振腔的 TE101 模式中,会有平行于最短边的电场,和旋转的电场;电荷会在最短边的两面上周期性地富集,同时会有环绕的电流。


趋肤效应(skin effect):在交变电流的场合,无电荷但有电流,有 E=0,×B=μ0j+1c2˙E,×E=˙B,B=0。对第三个式子两侧求旋度后,得到

-\nabla^2\b E=-\dfrac{\p}{\p t}(\nabla\times\b B)
\\=-\mu_0\sigma\dot{\b E}-\dfrac1{c^2}\ddot{\b E}&(\b j=\sigma\b E)
\\=-\i\mu_0\omega\b E+\dfrac{\omega^2}{c^2}\b E

因此有

2E+ω2μ0ε0(1iσωε0)E=0

在良导体中,传导电流 σE 远大于位移电流 ω2c2E,因此其退化为

2EiωμσE=0

有特征根 λ=±ωμσ/2(1+i)。则其有震荡项 E0 和 decay 项,decay 项由 δ=2/ωμσ 决定,满足

E=E0exp(z/δ)

则其关于距离指数衰减。


波导:沿着 z 轴截面不变的中空金属筒。

波导分析时,不适宜用磁感应强度 E,而是用磁场强度 H。有 B=μH

将电场分割为沿波导方向(z)方向的电场和垂直波导方向的电场。则有

E(x,y,z)=(ˉe(x,y)+ez(x,y)ez)exp(jβz)exp(jωt)

其中,e(x,y) 是垂直波导 (transverse,横) 分量,ez(x,y) 是沿波导 (longitudinal,纵) 分量。

同理,

H(x,y,z)=(ˉh(x,y)+hz(x,y)ez)exp(jβz)exp(jωt)

Ex(x,y,z)=ex(x,y)exp(jβz)Ey(x,y,z)=ey(x,y)exp(jβz)Hx(x,y,z)=hx(x,y)exp(jβz)Hy(x,y,z)=hy(x,y)exp(jβz)Ez(x,y,z)=ez(x,y)exp(jβz)Hz(x,y,z)=hz(x,y)exp(jβz)

假设波导里无电荷或电流。于是有

×E=μ0˙H=jωμ0ˉH

其中,ˉH 是分离出的不含时 H,即 H=ˉHexp(jωt)

同理,

×H=ε0˙E=jωε0ˉE

易知:二者的空间部分都满足 Helmholtz 方程,即

2ˉE+k2ˉE=0ˉE=0

各自把 Helmholtz 的分量展开来,两个 Laplacian 生成了六个方程。

总之而言,得到

Hx=jk2c(ωεEzyβHzx)

其中 k2c=k2β2,k2=ω2/c2

  1. 考虑 TEM 波的场合。此时 Hz=Ez=0。回到 Helmholtz 的六个方程,得到需要 kc=0 以使得 Hx0。最后推了依托后,得到:ˉe(x,y)=tΦ(x,y),其中 t=(x,y)。并且,有 E=te(x,y)+Ezz=0,由此得到:2tΦ(x,y)=0。那么,净电场和净磁场的方程即与静电场的方程相同,于是只需得知波导壁上的电磁场分布,即可根据边界条件推出内部的电磁场。但是在纯导体环境下,导体壁上不应有电场分布,因此得到结论:波导中没有 TEM 波!!!但是相反,在 transmission line 的场合,因为内外管壁存在电场差,因此这个场合存在 TEM 波,其不含时分量与双筒电容的电场分布相同。【综上,TEM 只在至少有两个导体的场合才可能存在】
    此时,定义传输线的微波阻抗(wave impedence) ZTEM=ExHy=ωMβ=μ0ε0=:η。【VI=EdlHdl
  • TE 波:Ez=0 但是 Hz0
posted @   Troverld  阅读(50)  评论(1编辑  收藏  举报
相关博文:
阅读排行:
· 物流快递公司核心技术能力-地址解析分单基础技术分享
· 单线程的Redis速度为什么快?
· 展开说说关于C#中ORM框架的用法!
· Pantheons:用 TypeScript 打造主流大模型对话的一站式集成库
· SQL Server 2025 AI相关能力初探
点击右上角即可分享
微信分享提示