高等等离子体物理(一)线性理论(研究生教材)王晓钢北京大学物理学院2009 年2 月等离子体的流体理论1.等离子体的流体描述1.1 等离子体的双流体模型1.2 Hall 磁流体(Hall-MHD )模型1.3 电子磁流体(E-MHD )模型1.4 理想磁流体力学(MHD)方程组1.5 位力定理1.6 变分原理2.理想磁流体平衡2.1磁场与磁面2.2Z-箍缩与花箍缩2.3一维平衡与螺旋箍缩2.4Grad-Shafra no 方程3.等离子体的理想磁流体稳定性3.1能量原理3.2扭曲模与交换模3.3 一维稳定性,直柱托卡马克4.磁流体力学波4.1线性磁流体(MHD )方程4.2非磁化等离子体中的磁流体波4.3磁化等离子体中的磁流体波5.均匀等离子体中的波(双流体理论)5.1 双流体模型5.2 介电张量与色散关系5.3 静电波简介5.4准静电波与准电磁波5.4电磁波简介1.等离子体的流体描述1.1等离子体的双流体模型等离子体是由大量带电粒子组成的物质状态。
一般意义上的等离子体由带正电的离子和带负电的电子组成。
由于带电粒子之间的Coulomb长程相互作用,等离子体呈整体电中性,即总的正电荷与负电荷相等。
因此,除特殊的非中性(一般是强耦合的)等离子体之外,我们可以用带负电的电子流体和带正电的离子流体组成的“双流体”模型来描述等离子体的宏观行为。
这种近似牵涉到等离子体时空尺度的讨论,我们在后面将进一步详细论述。
基于流体力学的图像及其近似,或者从统计物理的分布函数及其满足的方程(如Vlasov方程或者Fokker-Planek方程等,取决与碰撞项的形式,这里用类Markov过程的碰撞项(f° - f)/ •三丫(f° - f))出发,我们得到“双流体”方程组:连续性方程(统计方程的零阶矩)平' n:=0,(I-01):t动量方程(力平衡方程,统计方程的一阶矩)n:m「u:. U =, u B ,n:.q:. E n.m:m,(1-02)- c状态方程(对统计方程各阶矩的“不封闭链” (Hierarchy )的一种截断)芒u -. 5 . 一u;(I-03)-tCoulomb 定律(Poisson 方程)v E= 4八n:q:., (I-04)Fayraday 定律 12Bc .t 这里:=i,e ;对〉类粒子来说:n :.是粒子数密度, 荷,u :.是流体速度,p :. = n :「.是理想气体近似下的分压强;而’-:是〉类与1类 粒子之间的碰撞频率(当:•二:时为自碰撞)。
E , B , J 则分别是电场强度、磁 感应强度、和等离子体电流密度。
关于状态方程,我们以后会进一步讨论。
这里我们只是指出:参数 的取值决 定等离子体的状态,如等温(isothermal )状态对应卞';不可压缩状态对应 Y T 血;其它的Y 值对应“绝热”状态。
1.2 Hall 磁流体(Hall-MHD )模型一般来说,双流体模型是描述等离子体宏观(大于粒子回旋半径的尺度)运动的有力工具;在高频波段也可以应用,甚至在回旋半径的尺度上也可以得到一 些有用的结果。
但是,由于电子与离子质量之间超过三个数量级的差别,在具体 计算双流体模型的时候,会遇到所谓“刚性”问题:即电子已经完全改变了运动 状态,离子还基本没有动!这使得我们在计算离子时空尺度下的物理问题时,耗费大量的计算机时间。
而且由于code 本身的精度,即使经过长时间运算看到了离 子的运动,其结果或者是看到了很强的数值不稳定性、或者是很难令人相信。
而Ampere 定律J 丄兰 c c jt(1-05)Gaussion 定理 (I-06)(I-07)m 是粒子质量,q.是粒子电的是所谓电子磁流体(Electron MHD )模型。
这个模型也是在电子和离子的运动为了稳定code 引进的数值耗散,则往往带来人为的非物理的效应。
即使对纯理论 的解析推导,不仅过程繁杂,而且得到的物理图像也不清晰。
所以我们经常引进进一步的近似。
因为离子运动的时间尺度远远长于电子的时间尺度(通常在 40倍以上,对于 回旋运动来说则可以大于1840倍),所以我们在主要考虑离子运动时,可以认为 电子响应是“瞬时”的(instantaneously or simultaneously )。
这样,我们可以保持其它方程不变,近似地把(1-02 )中电子的质量趋于零,得到:这里 =4>e /「2e 是所谓的“ Spitzer 电阻”。
利用在“准电中性”近似n 「理下,U e 二* - J /ae ,J = me u i - n e e% :n e e u e 是等离子体电流,这个方程可以写为:必须注意到:这里我们还用到了 u e …U j 的条件。
明显地,u e - u i 要求电子运动与离子运动的分离,即所谓Hall 效应。
所以我们称这个近似模型为 Hall 磁流体模 型;这个方程则称为 Hall 磁流体的广义欧姆定律(Hall MHD Generalized Ohm ' s Law )。
方程中的JXB /n e ec 明显地就是我们在电磁学课程里熟知的 Hall 电场项。
1.3电子磁流体(E-MHD )模型而另一方面,我们在主要考虑电子运动时,可以认为离子响应是“无穷慢”的,或者说离子可以看成是保持总体电中性的 “背景”。
或者说,把离子看成是“稳态”的(::/::t =0,但是可以有ui = 0 )。
将/ ::t =0的近似带入离子的方程,得到 分离的情况下得到的,适用于比Hall MHD 模型更小的空间尺度和更快的时间尺度 的问题。
泌 - n °e u e , n e e(l-02e)E Ur_Bc J B n e ec IPe n °e J 汉B (J 「n i e u i): n ec1.4理想磁流体(MHD )方程组如果不仅整体等离子体呈电中性,而且在非常小的局部也呈电中性,我们可以把这个局部取做流体元,则有n e =n j = n。
上面(I-04)的右边等于零,而(1-02)的不同电荷粒子方程相加可以消去小尺度下(即流体元)的电场。
这样可以在很多情况下使问题得到简化。
这个图像,我们称为磁流体(magnetohydrodynamics, MHD)近似(或称“磁流体力学”近似)。
在这个近似下,宏观的“大尺度”电场满足的方程可以由(I-02 )两式之差(得到的Ohm定律)来计算。
这一节里,我们详细讨论这一近似。
等离子体过程的时空尺度研究物理问题时首要的是讨论时空尺度。
经典的宏观(大空间尺度)、低速(慢时间尺度)下的牛顿力学与相对论(快时间尺度)、量子力学(小空间尺度)的适用范围就是典型例子。
在等离子体中存在着很多的运动模式,我们无法、也没有必要同时考虑所有这些运动模式。
那么哪一种(或者几种)运动模式是主导的、起着决定作用的?要回答这个问题,就要进行时空尺度分析:我们关心的是哪个时空尺度下的物理问题,在这个时空尺度下存在哪几种运动模式?所以,对于等离子体这样的存在大量运动模式的连续介质来说,时空尺度分析尤其重要。
磁流体(MHD )理论的基本假设磁流体理论本质上来说是一种与流体力学相类似的连续介质的理论。
因为考虑宏观的大尺度问题,其特征长度L H(或L o)一般可以看成是所研究等离子体区域的大小,比如柱形等离子体的横截面(的半径)。
而特征时间尺度沖(或者特征频率-.H ~1/ -H)则可以用一个特征信号穿越这一尺度的时间来表征。
这等于特征尺度L H与特征信号在等离子体这一介质中传播的速度之比。
在流体理论里,这1/2显然是声速C s二p/^m ,这里'm是质量密度。
但是在磁化等离子体中,对于大尺度的MHD问题来说,这一特征速度是所谓 Alfv e速度2 B/ 4二亠1/2。
当然,如果所研究的等离子体可以看成是一个驱动(driven)系统,那么其特征时间尺度应该由驱动频率给出。
磁流体(MHD)理论基于下列假设:—1 E* 非相对论假设:■ ■ / k ~ L H / ■ H ~V typical I: c ■0yp C c t*流体假设:1)局域热力学平衡(局域 Maxwellian分布)假设:F ::::::・H ~・H(要求较咼碰撞频率:压强是标量,n—■ 0);2)忽略有限Larmor半径(FLR)效应:• 'H 2」ci :We,O/L H—【/L H X 1;3)单流体(准电中性)假设(即Debye球内有大量粒子,也称等离子体假设):n 「::’D I:L H,‘H I:• 'pe,= '•> 0。
a我们会发现,局域Maxwellian分布的假设对于“无碰撞”理想(ideal)等离子体(其平均自由时间-ee, -ii, -e^ ■ -H )来说的不是一个好的假设。
我们需要进一步讨论:1)粒子间“碰撞” (collision )和关联(correlation )之间的关系,以及长程碰撞的“集体”(collective )效应和短程碰撞之间的关系;2)以及导向中心理论的回旋动理学(gyrokinetic )和漂移动理学(drift-kinetic) 近似。
磁流体(MHD )方程组如果我们利用m二n i二n ,定义小的等离子体元的“单流体”物理量:质量密nm u i n e m e U e m e度:P m三n (m+讥)茫nm,流体速度 U =------------------ = u i +— u e纭u i,等离子nm + n e m e m i体压强p = n (Te+T i ),等离子体电流J= ne( u i -u e),将(1-01 )、( 1-02 )分别对不同电荷分量求和得到连续性方程丄 i n u =0,:t动量方程(并利用Ampere定律)而电子的动量方程(I-02 ) (:•二e时)可以写成这里、e包括了电子自碰撞>ee及电子一离子碰撞Ji ;或者J B ' P enec ne(I-08)nm i 兰u::t5「卩亠7丿 c 4(I-09)E u^_Bc亠込u ene em eueJ i u e这个方程我们一般称为广义Ohm定律,其中='e m e / ne?为经典的等离子体(1-14)Spitzer 电阻率。
一般来说电流主要是电子携带的,可以忽略最后两项得到常用的 广义Ohm 定律(1-10)所谓理想磁流体近似,在“无碰撞”近似下忽略方程右边的电阻项,在大尺度磁 流体近似下忽略方程右边其它各项(Hall 电场项由离子惯性尺度、电子压强梯度 项由“离子声回旋半径”尺度、电子惯性项由电子惯性尺度表征),得到理想磁流 体的Ohm 定律u KB E0。