等离子体物理讲义06磁流体力学及静平衡12汇总.docx
《等离子体物理讲义06磁流体力学及静平衡12汇总.docx》由会员分享,可在线阅读,更多相关《等离子体物理讲义06磁流体力学及静平衡12汇总.docx(39页珍藏版)》请在冰点文库上搜索。
等离子体物理讲义06磁流体力学及静平衡12汇总
等离子体物理学讲义
No.6
马石庄
2012.03.07.北京
第6讲MHD方程与静力平衡
教学目的:
建立等离子体的磁流体模型,在拟稳态近似下,建立磁流体动力学方程。
依据磁Reynolds数,掌握理想MHD的磁冻结定理和拓扑不变量;无力平衡和有力平衡。
主要内容:
§1MHD方程(3
1.1导心理论引出(3
1.2MHD近似(9
1.3磁应力张量(12
§2电磁感应方程(15
2.1磁冻结定理(16
2.2拓扑不变量(21
2.3磁场扩散(26
§3MHD静平衡(28
3.1维里定理(30
3.2无力平衡(34
3.3有力平衡(36
习题6(44
在研究等离子体的宏观运动时,通常可以近似地把它当作导电流体来处理。
这种模型适合于缓慢变化的等离子体现象。
所谓缓慢变化是指等离子体的特征长度和特征时间远大子等离子体粒子的平均自由程和平均碰撞时间。
在这种情况下,等离子体可以近似地看作处于局部热平衡状态,因而可以像通常的流体力学中那样定义流体的速度,压强,密度,温度等流体力学及热力学参量并用这些宏观参量来描述等离子体的宏观运动。
§1MHD方程
当导电流体在电磁场中运动时,流体内感生出电场从而产生电流。
这个电流一方面与磁场相互作用,产生机械力,对流体运动产生重大影响;另一方面感应出改变原有电磁场的磁场。
于是就形成了电磁现象和流体动力学现象相互作用的复杂图像。
这些现象必须要用电磁场方程和流体动力学方程的联立方程组来进行研究。
1.1导心理论引出
等离子体中的带电粒子在电磁场中的运动可以看作是围绕磁力
线回转的粒子引导中心的漂移叠加,下面探讨微观单个粒子的行为与宏观流体行为之间的关系,给出一种物理直观图象。
如图1所示,基本思路是计算导心运动导致的流过等离子体中任意开曲面的垂直电
流密度,考察这个电流与等离子体压强梯度和惯性力之间的联系。
取曲面的法向与磁场正交,仔细考虑回转半径扩张的影响。
首先考虑粒子运动的主要贡献是来自圆周回转运动,每个粒子进出曲面的方向相反,对电流没有贡献,如图1(b。
换言之,在一个回转周期中,没有净电荷流动。
垂直电流由两种不同的机制产生。
一个是导心垂直漂移产生的穿过曲面的电荷流,如图1(c;还有一种曲面边界附近的回转运动,如图1(d,所谓磁化电流。
粒子的导心漂移速度由漂移,B漂移,曲率漂移和极化漂移构成
EB
2
d
d
EB
2
其中/
·
是磁力线的曲率半径。
在磁流体力学尺度内,通常漂移比电子和离子的任何其它漂移都要大量级,粒子的垂直速度和平行速度与热运动速度相当,即~~。
所以,在流体模型中,无论电子还是离子,导心的主要运动都是漂移运动。
由于电子和离子以相同的速度漂移,因此可以引入垂直方向上的宏观速度~~,这里
改写为
相当于垂直方向上的Ohm定律。
现在考察导心漂移运动产生的垂直电流。
尽管漂移对粒子漂移的贡献最大,由于电中性条件以及电子离子的漂移方向相同,漂移产生的净电流也为零。
再考察其它漂移的贡献,将电子的漂移与离子的漂移相减,然后对所有的粒子求和
~d
代入漂移速度,漂移自动消去
d
d
d
假设粒子服从稳定Maxwell分布,在速度空间的局域直角坐标系中
容易计算得到
d
d
其中,。
磁化电流源自边界附近只穿过曲面一次的带电粒子,它们对垂直电流的贡献要通过计算每个粒子携带电流与所有只穿过一次的粒子数的乘积来估计。
在一个回转周期内,载荷的粒子的产生的平均电流为
2
负号表明电流沿逆磁方向,趋于抵消原来的磁场。
这类带电粒子可能具有四种不同的情形,如图2所示。
这些轨道中心的轨迹组成一个以回旋半径L为半径的圆,任何带电粒子只要导心位于所选曲面边界上的点为中心,以L为半径的圆内,就属于仅穿越曲面一次的粒子。
包含这些粒子的体积元为d·d,其中
L,d是曲面边界元。
一般说来回转轨道的法向沿方向,与d并不平行,如图3所示,其投影减少了单次穿越所选曲面的粒子数
目。
体积元内速度为的粒子数为ddd,相应的电流为ddd,通过对速度空间和沿曲面边界的所有体积元积分,
即得总磁化电流
M2d·d对于局域Maxwell平衡分布,积分并引入M,得到
·dMM·d·d·d比较得到
M
合并导心漂移电流和磁化电流M
dd
化简各项,得到
dd
其中用到恒等式·0和··/2,合并消去,从导心理论得到动量方程
d
这个形式与理想磁流体力学得到的结果一致,说明流体力学方法和导心理论方法两种处理是等价的,只是对宏观磁流体力学行为的自洽描述不同。
1.2MHD近似
电磁现象的一般规律满足Maxwell方程
∂∂··0
本构方程是Ohm定律
这里假定介质是静止的(对参考系而言,,等就是在这个参考系内定义的。
特别是电流与电场之间的关系,一般说来只适用于静止的导体。
为了求出运动导电流体内电流和电场的关系式,从参考系变换到另一个以速度相对于运动的参考系,其中导电流体在所考虑的时刻是静止的。
在这个参考系内,有,其中是内的电场强度。
根据相对论关于场变换的公式,准确到/的量级,用系内的场表示为
于是得到
为了得到参考系中的表达式,考虑电流密度的一般定义和,其中和分别是电荷在座标系和内的速度。
在和的相对速度远小于光速的非相对论情形,速度按Galilio公式变换
由此得到
即
下面将表明,磁流体力学范围内运流电流与传导电流比较可以略去。
因此,从一个参考系变换到另一个参考系时,电流实际上保持不变。
物质内的交变电磁场的特征,主要决定于物质的种类和场频率的量级。
在磁流体力学范围内,通常研究的是在外加交变磁场内的大导电流体中产生的现象。
这时可以假定,场的变化速度不很大,满足这样两个条件。
首先,假定相应场频率的波长~/,大于流体运动的特征长度,即
或
1
其中~1/为场变化特征时间。
其次,假定电导率和场频率之比满足:
1或
1
即场变化的特征时间远大于粒子碰撞时间,由于等离子体是良导体,这一条件实际上总是满足的,这样的电磁场和电流是准静态的,称为MHD近似。
当MHD近似条件满足时,位移电流∂⁄,运流电流和电
场力可以忽略。
利用,·及∂/,1/,并且假定,,可以作如下估计
|∂⁄||||⁄|||1||||·1
||||||
1|||||·|||1
||||||||
1于是运动导电流体内的MHD近似电磁场方程为
∂Faraday定律Ampere定律··0
Ohm定律
在MHD近似下,进一步从Ampere定律和Ohm定律中消去电流,得到
两边取旋度
运用恒等式
··
由于·0和Faraday定律,得到描写磁场演化的电磁感应方程
∂其中磁扩散系数
1
电磁感应方程说明磁场随时间的演化受两项控制。
第一项与流体速度有关,称为对流项;第二项与电导率有关,称为扩散项。
1.3磁应力张量
对于电磁场中运动导电流体,仍然要从质量守恒、动量守恒和能量守恒定律得到动力学方程。
现在运动方程中的外加质量力是电磁力和重力。
通常情况下电磁力远大于重力,可略去不计。
考虑了电磁力的流体动力学方程如下
dd·0d·d····电磁场通过Lorentz力公式
对流体产生作用。
从而。
磁流体动力学方程组为
dd·0d·∂d····如果流体是无粘性,不传热的和理想导电,称为理想导电流体,方程组为
dd
·0d·dd0∂
称为理想MHD方程。
进一步分析Lorentz力的作用性质。
将代入得到
1
121·事实上,由···,体积力可以写成
1·12
·这里
112
为Maxwell电磁应力张量的磁场部分。
现在解释给此力以物理意义。
考虑闭合曲面任意体积,在此体积上积分得到合力,体积分可以变换成面积积分,于是合力可以解释成一组等效的表面力
d·d
·
d
d
其中
1
·
1
2
是作用在单位法线矢量为的面元上的等效表面力,称为磁应力。
设平行于磁场的单位矢量为,则体积力和磁应力的关系进一步写成
d·
d
2
d
2
cosd
其中为和
之间的夹角。
由此可见,磁力等效于大小为2⁄的各向同性磁压力和沿磁力线方向,大小为⁄的张力之和,或者可以说,磁力等效于与磁力线相垂直方向上的压力2⁄与沿磁力线方向的张力⁄的之和。
实际上,利用矢量公式
···
和,还可以得到
11·11
·2·2
·其中,最右边第一项是各向同性的磁压力;第二项在平行方向抵消磁压力梯度力,因为本来就没有平行磁场的分量;第三项是磁张力引起,与磁张力和曲率成正比,是弯曲磁力线的恢复力,指向曲率中心。
由此可见,用形象语言说,磁力线很像拉紧的橡皮筋,沿着磁力线方向是张力,磁场增强也就意味着张力增大。
如果磁力线是弯曲的,这个张力可产生指向磁力线曲率中心的恢复力。
磁力线被扰动就像张紧的橡皮筋被拨动,扰动可以传播,就是所谓的磁流体波。
§2电磁感应方程
取是流体空间变化的特征尺度,是磁场的特征尺度,是速度的特征尺度,把用估计,将各物理量标度化
,
其中各无量纲物理量的量级||~~||~1,类比于流体力学中的Reynolds数,引入无量纲磁Reynolds数
||||
得到无量纲电磁感应方程
∂
1
对流项与扩散项量级之比的估计
||
||
·
||
||
~
如果1,对流项为主;如果1,扩散项为主。
2.1磁冻结定理
大磁Reynolds数意味着当1时,电导率,流体速度和特征尺度足够大,使得1,扩散项比起对流项足够小,可以忽略不计,电磁感应方程简化为
∂
这个方程可以用来证明称为“冻结磁场定理”。
这个定理有两个等价的表述,分别有不同的证明。
表述1:
穿过任一随流体一起运动的闭合曲线的磁通量为常数。
在时刻,对于有闭合曲线围合的曲面,磁通量定义为
Φ·d
设曲线随流体一起运动,时间后,曲线上的各点移动到新位置d,则通过曲线的磁通量的变化率为
dΦ
∂
·d
·d
第一项是由于磁场对时间的显式依赖引起的磁通量变化,第二项是由于曲线的移动引起的磁通量变化。
第二项的矢量三重积可以交换位置
dΦd
∂
·d
d·
根据Stokes定理及理想电磁感应方程,得到
dΦ
∂
·d
·d
因此
dΦ
表述1有个推论,考虑两条闭合曲线和在时刻为磁力线所连接,这些磁力线形成一条磁通量为Φ的磁通管。
随着流体运动,既然在组成磁通管的曲面上dΦ/d0,在以后任何时刻,这条磁通管的磁通量保持不变。
由于磁通管的两端可以收缩为无穷小,磁通管就
变成一条磁力线,因此Alfven可以说“磁力线冻结在流体之中”。
必
须指出,磁冻结定理是是说流体的运动不能穿越磁力线,沿磁力线方
向的流体运动不违反磁冻结定理。
表述2:
如果随流体运动的一曲线初始时为磁力线,那么随后就
一直保持。
为了证明这个命题,需要引入通量坐标(fluxcoordinates。
由
于磁场是无散场,只有两个自由度,因此总是可以用两个标量场可以
表示,使用Clebsch变量,,也称为磁Euler势
其中,是沿磁力线是保持不变的常数
··0
对于空间一点,磁力线可以表示为
用表示从沿的曲线长度。
三维曲线正交坐标系,称为,,,在磁力线附近,0,自动满足
····0
其实,是沿磁力线保持不变的常数,因为
··0
证明表述2,基本的思路是计算沿给定的流体质点轨迹,计算证明·的时间变化率为零。
为此
d
·∂
···
其中第一项
∂
·∂
··
∂
根据电磁感应方程,有
∂
···设恒定的曲面随流体一起运动,即要求
d
∂
·0,
∂
·
分别代回,得到
dd
···········
方程等号右端的第一,三和四项可以简化
······
··
因此
d···表明,如果在0时·0,那么对所有的时间·0。
如果在0以通量坐标,的常数曲面标志一条磁力线,那么此后一直标识这条磁力线。
2.2拓扑不变量
在理想MHD
中,任何有限体积中包含无数多条磁通管。
关于磁
通量管,具有拓扑不变量,在MHD湍流和快发电机理论有重要应用。
对于有闭合曲线围合的曲面,磁通量定义为
Φ·d
·d
d
考虑积分
·d
其中是中的第条磁通管的体积。
磁通管以速度与流体一道运动,则的时间变化为
d
∂
·d·
∂
d·
d
d
其中最后一项
d
dd
ddddddddd·d
使用Faraday定律,有
d
·d
·d
··d
其中是标量势。
由于
···第二个积分可以写为
·d·d
·d
·d
·d
类似地,第一个积分可以写为
·d·d·d
因为·0和在上·0,是磁通管。
因此
dd2·d·d··d
理想MHD有,就有
dd····d
0因为在在上·0,·0。
因此,对于系统中每一根磁通管const.。
称为Woltjer不变量1,它取决于理想MHD中和在上·0,·0。
这两个条件中,第一个是磁通管的性质,第二个也是理想MHD的结果,磁通管随流体运动。
可以给Wöltjer不变量以
物理解释2。
考虑如图所示的联
通磁通管。
通量管包含磁通
Φ,通量管包含磁通Φ,通
量管的Wöltjer不变量为
·d
对这条通量管,有
1L.Wöltjer,Proc.Nat.Acad.Sciences44,4891958.
2H.K.Moffatt,MagneticFieldGenerationinElectricallyConductingFluids,CambridgeUniversityPress,Cambridge,UK1978.
dddd
dddddd
ddd·dd
从而
·d·d·d
第一个积分刚好就是通量管内包含的磁通Φ;第二积分是曲线围合的磁通,如果通量管有“右手”连通的,磁通是Φ;如果通量管有“左手”连通的,磁通是Φ。
因此
ΦΦ
类似地
ΦΦ
如果通量管缠绕了圈,那么有
ΦΦ
同样的结果也可以对单结通量管得到,如下图。
综上所述,Wöltjer不变量是通量
管的连通性,即拓扑(topology的度
量,成为磁螺旋度(magnetichelicity。
因为理想MHD中const.,
这意味着,通量管的拓扑不变并且永远保持。
这个
性质因此是一个拓扑不变量,确实是磁场与流体一道运动的另外一种表述,是理想MHD的Ohm定律0的结果,对流体的可
能运动给以很强的约束。
交叉螺旋度(crosshelicity
·d
时间变化率为
d
·∂
·
∂
d
··1
1
·d
积分号下的第三项涉及Lorentz力,对积分的贡献为零。
如果假设绝热过程~,由于·0,第二项为
1
··
运用矢量恒等式
·
2
和
···第一和第四项结合得到散度形式
···
运用散度定理,的时间变化率为
d
·d
在边界上,当··0时,上述积分为零。
在这些条件下,
交互螺旋度是理想MHD的不变量。
2.3磁场扩散
当1时,意味着等离子体的电阻效应占主导地位。
从磁通量Φ变化的角度,显然有
dΦd·d
1
·d
磁冻结定理不再成立,流体质点可以与磁力线分道扬镳,磁力线可以穿越流体质点移动。
由此,等离子体中的磁力线才会发生MHD重联(reconnetion和MHD弛豫(relaxation现象.
若取是流体空间变化的特征尺度,是磁场的特征尺度,是速度的特征尺度,特征时间为,从无量纲电磁感应方程
∂
1
当∞时,方程
∂
1
有解指数衰减的解
exp
exp
exp
其中磁场扩散特征时间为
~
例如,对于地核,衰减时间为10sec~10年,相对于地球几十亿年的存在,如果有原生磁场,早已衰减殆尽。
维持磁场,需要发电机作
用。
对于磁约束,等离子体的有限电导
率也是一个需要克服的难点。
如图,设
等离子体被磁场所约束,处于准平衡状
态(流体的惯性项可忽略不计,等离
子体主要处在“内部”,磁场主要在“外
侧”。
等离子体的连续性方程为
∂
··
其中,离子质量密度为,动量方程为
由定态Ohm定律
其中为标量势,可以把电场是外加的。
定态Ohm定律与动量方程合,得到
1
代入连续性方程
∂·
·
在等温过程中,,const.,等离子体密度的演化依据
∂
··
其中扩散系数
且漂移速度
第一项表示穿越约束磁力线的扩散,第二项表示向内的运流作用。
扩散的特征时间为
因此,如果没有外加电场,等离子体在这个时间尺度上将损失。
换言之,等离子体的电阻对磁约束有负面效应。
由于有限电阻存在,磁能不再守恒,而螺旋度可能是不变量。
如下图操作,总有
2Φ
Figure:
Asequenceofhelicityconservingoperationsontwosinglylinkeduxtubes(a.Wecarefullymakeacutin(bandreconnectanddeformin(ctoshowtwo360_twistsandfurtherdeformin(dtoshowtwocrosses.
§3MHD静平衡
对于流体而言,动量方程可以写为Newton第二定律的形式
d
d
静力平衡(stationaryequilibirium要求
dd
非静力平衡(non–stationaryequilibirium要求
0·
选择坐标系统,使得流体是静止的,在MHD方程中,静力平衡(要求
0,
因此静力平衡
Ampere定律
和
·0
成立。
显然有推论
·0,·0
电流和磁场必定处于等压强曲面上。
前者说,沿着磁场方向压强是不变的;后者说,要保持一个与磁场不平行的电流需要存在压强梯度,反之亦然。
在理想MHD中,压强梯度力可以被Lorentz力平衡,这就是磁约束的基础。
3.1维里定理
在流体力学中,静力平衡是相对简单的。
考虑重力场中的流体,流体静力平衡的条件是
在一维情形,有
即得
d
如果密度剖面给定或者存在关系,就可以得到压强的分布。
在后面的情形中,
dddd
dd其中是声速。
在等温条件下,const.,得到解
expg对于非静力平衡,有
·
因此,有限大小的流体流动可以导致非均匀压强。
在一维情形中,有
d
1
即
1
const.
就是Bernoulli定理。
对于更一般的情形3,在MHD中,动量演化依据
·
平衡的条件表示为
·0
即
其中,由于0,
2
1
是总应力张量。
把重写为
2
比较方便。
定义等效的正交压强和平行压强
3ThediscussionoftheVirialTheoremfollowsthatofV.D.Shafranov,“PlasmaEquilibriuminaMagneticField”,inReviewsofPlasmaPhysics,M.A.Leontovich(ed.,ConsultantsBureau,NewYork(1965.
因此
总应力张量可以表示为
即
trace
引入置于磁场中的孤立流体,在平衡态,则有trace0成立。
在体积上积分并应用Gauss定理,得到
d
d
d
其中是包围体积的闭合曲面,未必是流体的体积及其表面。
求积分,得到
3
d
·
·d
等号左端是正定的,现在估计右端的符号和大小。
考虑所有的电流密度和压强都完全在流体内部,闭合曲面完全在流体之外,磁场只是有磁流体内部的电流产生,如图所示。
现在
取∞,使得上式右端只与磁场有关。
磁流体内部的电流引起的矢量势为
||
d
如果远离,那么可以做展开
1||
1
||
·
||
代回矢量势的表示式,既然·0,以及
d
发现
~,~
且
·d~
因此,积分的右端也是,当
0时趋于零。
与左边的保持正的相矛盾,意味着条件不可能为平衡所满足。
这就是维里定理(TheVirialTheorem。
表明,磁化流体不可能在内部电流产生的力作用下处于MHD平衡。
换言之,MHD平衡一定是为外部电流维持的。
在实验室中,维里定理一定满足,因为都是用外部线圈产生磁场的。
在天体物理条件下维里定理未必满足,而宇宙恰恰是充满动机的地方。
维里定理
RudolfClausius1822‐1888在1870年创立,又称均
功定理:
若粒子受力的作用,在某时间的位移有限,
则在外力长时间∞的作用下,所作功平均值·等
于粒子平均动能负值的两倍
·
·称为维力。
若系统由个粒子组成,则有
·
对于保守力场中的经典粒子,若作用位势是粒子位矢次幂的齐次函数,则力,由维里定理,可得
2
对重力场和静电场来说,因1,则
1
2
对于具有磁场和自转的稳定体系,维里定理表示为:
220
其中是系统转动的动能,是系统的总磁能,是系统的总位能。
3.2无力平衡
在稀薄等离子体中,压强梯度力可以忽略不计,称为无力平衡(force‐freeequilibria。
有
因此存在标量函数,使得,即
这样的是矢量场,称为Beltrami场。
显然
0····
由于·0,因此,
·0
标量函数沿磁力线方向保持恒定。
引入参数
再对运算,
·0
理论上,两个方程可以联立求解未知量和,实际上可能遇到微妙的数学问题4。
一个常见的问题是,const.的柱位形等离子体,则·0自动满足;0,则关于和分量分别满足方程
d
d
d
d
10
和
d
d
d
d
都具有Bessel方程的形式,解分别为
其中,和分别是零阶和一阶Bessel函数。
注意到当,2.4048,符号改变