Maxwell 方程边界条件
Maxwell equation:
∂tEi−ϵjki∂jBk∂tBi−ϵjki∂jEk=0=0(1)(2)
双曲对称的
记 uα={Ex,Ey,Ez,Bx,By,Bz}(symbol)
将方程(1)(2)写成矩阵形式可得:
∂tUα+Akαβ∂kUβ=0(3)
引进对称化子:
ds2=Sαβduαduβ=δαβduαduβ(4)
它可以将矩阵 A 对称化:
设
nk
是边界的外法向量,
nkAkab
是特征矩阵,
ea^b
是特征矩阵的左特征向量
ea^bnkAkba=va^ea^a(symbol 6)
方程对角化后的新变量(姑且叫作特征场吧)
ua^=ea^bub(7)
分解特征场:(至于为什么这么分解?忘记了)
u0^±=nkEk±nkBk(Q8)
u1^±i=PijEj∓ϵijknjBk(Q9)
- 特征矩阵的特征根
⎛⎝⎜⎜⎜⎜⎜⎜⎜⎜0011−1−1⎞⎠⎟⎟⎟⎟⎟⎟⎟⎟
⎛⎝⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜000n1n3n2n31n1n3n2n31000n3n3n2n1n22+n21n1−n2n110−n2n23+n21n1−n2n3n1−n3n101n2−n23+n21n1n2n3n1−n3n101−n3−n3n2n1n22+n21n1−n2n110⎞⎠⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟
改装以后的特征向量:
⎛⎝⎜⎜⎜⎜⎜⎜⎜⎜⎜000n1n2n3n1n2n3000n3n1n3n2n22+n21−n2n10−n2n1n23+n21−n2n3−n30n1n2n1−n23−n21n2n3−n30−n1−n3n1−n3n2n22+n21−n2n10⎞⎠⎟⎟⎟⎟⎟⎟⎟⎟⎟
注记:
-
A×B:=ϵijkAjBk
-
n⋅(n×B)=niεijkniBk=0
ϵijk=⎧⎩⎨10−1(i,j,k)=(1,2,3;2,3,1;3,1,2)(other)(i,j,k)=(3,2,1;2,1,3;1,3,2)
- ϵijkϵilm=δjlδkm−δjmδkl
- 0=niPij=ni(δij−ninj)=nj−nj=0
- nj=gijni
Definition projector
PijPijPij=gij−ninj=δij−ninj=gij−ninjExample Pijti=tj , Pijnj=gijnj−ninjnj=0 就是说, Pij 可以在切空间上给出一个度量。
Exercise proof that
nkSk=ϵijkniEjBk=Piju1+^iu1^+j−Piju1^−iu1^−jSniSi=E×B=ϵijkniEjBk=point out that niu1^∓i=0
niu1^±i=ni(PijEj∓ϵijknjBk)=niPijEj∓niϵijknjBk=0(Remark)平面波的例子
Consider very short wavelength plane-wave perturbations that propagate normal to a particular boundary
δuα=Vαe−iwt+iknaXa(Q)
k 在这里充当的是波数?k 足够大的时候
∂tδuα+Akαβ∂kδuβ≈i(wδαβ−knkAkαβ)Vβ≈0(QQ)