广义欧姆定律的推导

——郑春开 《等离子体物理》4.7节 P86

 

 

假设电子和离子密度相同ne=ni=n,并保持电中性,应用双流体方程组:

(1){nαt+(nαuα)=0nαmαduαdt=pα+nαFα+Rα32nαdTαdt=pαuα+Qα

其中α=i,e;在等离子体物理中的摩擦系数Rα通常可以写为布拉金斯基给出的形式:

(2)Rei=Rie=meneνei(ueui)

式子中νei为电子-离子碰撞频率。注意,同类粒子碰撞时因为其动量守恒,对摩擦阻力没有贡献,所以Ri=Rie, Re=Rei而且Rei=Rie,于是(1)式变为:

(3){nmi(uit+uiui)+pi=en(E+ui×B)+meneνei(ueui)nmi(uet+ueue)+pe=en(E+ue×B)+meneνei(ueui)

e/mi乘以(3)中的1式,e/me乘以(3)中的2式,然后两式相减:

(4)ne(uituet)+ne(uiuiueue)+e(1mipi1mepe)=e2n(1me+1mi)E+e2n(1miui+1meue)×Bmeneνei(1mi+1me)(uiue)

电流密度

(5)j=ne(uiue)

使用以下近似化简(4):

(6)memi, 1me±1mi1me,u=nimiui+nemeuenimi+neme=miui+meuemi+memiui+meuemi

因为pi=niTi, pe=neTe, ni=ne, TiTe

(4)中的第一项:

(7)net(uiue)=t[ne(uiue)](uiue)ent=jt(uiue)entjt

上式忽略了(uiue)ent这个二级小量项目,是因为平均速度都很像且速度与密度随时间变化环面。同理(4)的左边第二项ne(uiuiueue),右方第二项的因子可以变为:

(8)e2n(1miui+1meue)=e2nme(meui+miuemi)=e2nme[miui+meuemi(miuimeuimimiuemeuemi)]=e2nmeeme[ne(1memi)(uiue)]=e2nmeuemej

于是(4)整理为:

(9)mene2jt=(E+u×B)1nej×B+1nepe1σcj

其中σc=e2nmeνei是等离子体电导率。

整理成空间物理学经常用的形式:

(10)E+u×B=1σcj+1nej×B1nepe+mene2jt

其中1σcj表示为电阻项;1nej×B是霍尔项;1nepe是电子压强梯度项;mene2jt代表电子惯性项。