磁流体力学不稳定性

《磁流体力学不稳定性》由会员分享，可在线阅读，更多相关《磁流体力学不稳定性（13页珍藏版）》请在金锄头文库上搜索。

1、第6章磁流体力学不稳定性6.1 概论等离子体能够被磁场约束并处于力学平衡状态。一个处于力学平衡状态的等离子体位形, 当它受到某种扰动，偏离平衡态时，等离子体将如何反应？是越来越偏离平衡态，最后导致 平衡态被破坏呢，还是很快将扰动抑制住回到平衡态.前者是不稳定平衡，后者是稳定平衡.但当磁流体处在非热力学平衡态，其内部存在着可以转换成扰动能量的自由能时，在合 适的条件下有些扰动就可能发展成为在大范围、长时间、能量超过热噪声水平的大幅度集体 运动.这种集体运动就称为不稳定的模式，相应现象就称为磁流体的不稳定性.研究等离子体的各种不稳定性，阐明其物理机制，探索抑制不稳定性的方法，一直是受 控核聚变研究

2、的重要课题.磁约束等离子体可以处于力学平衡状态，但它不是完全的热力学 平衡态.等离子体处于非热力学平衡状态意味着等离子体具有较高的自由能，因而必然会产 生从较高能量状态过渡到较低能量状态的宏观或微观运动.等离子体偏离热力学平衡态大体 有两类方式.一类是等离子体宏观参数如密度、温度、压强或其它热力学量的空间局域性和 不均匀性；另一类是等离子体的速度空间分布函数偏离麦克斯韦分布.由于前一种原因产生 不稳定性时，等离子体通常以整体形式在空间改变其形状，因而称为宏观不稳定性。由后一 种原因产生的不稳定性称为微观不稳定性.宏观不稳定性通常用磁流体力学方程进行分析， 因而也称为磁流体力学不稳定性，而微观不

3、稳定性则用动力论方程进行分析，因而也叫动力 学不稳定性.由于磁流体力学不稳定性在磁约束核聚变等离子体中具有更重要的地位，处理 方法也相对地比较容易，因此本节仅讨论磁流体力学不稳定性.下面我们将首先从分析流体 的瑞利泰勒不稳定性(Rayleigh Taylor instability )入手，这样做物理图像清晰，易于理 解.然后讨论在分析磁流体力学不稳定性中得到广泛应用的能量原理.在这基础上分析几种 主要的宏观不稳定性，最后讨论等离子体电阻对不稳定性的影响.下面是几种典型的磁流体不稳定模式.例1 .瑞利-泰勒（Rayleigh-Taylor ）不稳定性（图4. 1）；例2 .开尔文亥姆霍兹（Ke

4、lvin Helmholtz ）不稳定性（图 4. 2）；例3 .腊肠型不稳定性（图4.3）;例4.弯曲型不稳定性（图4.4）;例5.磁岛（图4.5）;例6.磁重联（图4. 6）.图磁岛116第g学 班温体力学不鹤定性阳4, 6股重联每种不稳定的扰动在其演化过程中都会依次经历下面三个阶段：线性阶段、非线性阶段 及饱和阶段.在线性阶段，扰动的幅度较小，不同类型的扰动彼此之间并不相互作用，扰动 对它所处的平衡态也无影响，这时扰动的幅度是随时间指数增长的.在非线性阶段，扰动幅 度增大到会反过来使原有的平衡量作一定调整 （因此改变了自己得以不稳定增长的初始条件, 使馈入的自由能量减少），并达到开始和其

5、他扰动模式相互作用（从而彼此间交换能量）的程 度，从而使增长率木断下降.这时扰动幅度是依次随时间的不同幕次（一般是从高幕到低幕 次）而增长的.当时间的幕次最后降低到零时，就达到了演化的终点一一扰动的幅度不再随 时间增加，而一直保持极大值，这就是饱和.本章只讨论磁流体的线性不稳定性.线性不稳定性的基本描述方法（1）简正模法先将描述所研究对象的状态量写成平衡量（零级量）和扰动量（一级小量）之和，然后 把它们代入所用的磁流体方程组，从中减去平衡方程并略去二级小量就得到了线性化的方程 组.对这些方程作（时间）拉氏变换和（空间）傅氏变换 A（r ,t）=Ak exp（i k r-it）后可能出 现下列几

6、种情况：全部空间坐标都能进行傅氏变换.这样线性微分方程组就变成了线性的齐次代数方程组，它的有非平凡解的条件（系数行列式为零）就给出了关于 o=6（k）的色散关系.例如上一章中平板几何位形下的阿尔文波的色放关系正是由这种方式得到的.（ii）只有部分空间坐标能进行傅氏变换，剩余的坐标构成了约化的微分方程组.这时要设法先得到它的通解，然后利用边条件或连接条件也可以得到（o=Mk）的色散关系.例如上一章中，柱坐标下阿尔文波的色散关系就是这样求得的.（iii）所得出的约化微分方程如果是奇异的，如上一章中连续谱阿尔文波所满足的方程（2）能量原理（仅对理想磁流体适用）6.2 瑞利泰勒不稳定性这是一种经典的流

7、体不稳定性.因为这种不稳定性是由重力驱动的，故又称重力不稳定 性.让我们来研究图3.25所示的一个容器.该容器内盛有两种不同质量密度的液体，上面的 液体质量密度大，下面的质量密度小.两种流体之间有明显的分界线.显然，质量密度梯度VP 由下向上，受到的重力由上向下，用 -VG来表示.液体的平衡方程是图3.25流体的瑞利-泰勒不稳定性三（：u） = 0二 t:业W Gdt式中u是流体元的速度.流体达到平衡 u = 0.现在假定在交界面上出现了一个微扰动，其形式为Pi = Pi（x）e,3,ui = ui （x）e,s这样，密度和流体速度便可写成:二：01, u = u0ui = ui从这里开始，参

8、数下标为0表示平衡量，参数下标为i表示扰动量.将（4）式代入平(5)衡方程（3）,我们得到质量守恒方程（：词）=ui 7Q =0 5=0.将（3）式代人（5）式便得到已表.t在整理上式时，已考虑到流体是不可压缩的， 达式:同样可以得到扰动后的动量方程和 U1的表达式:：0包!= 一个G dtI-q G0将(6)式和(8)相结合使得到如下的方程:(6)(8)(9 )(9)式说明，当流体的密度梯度方向跟受到的重力方向相反时就会产生不稳定性，此时。20,这就是说重流体在上面轻流体在下面的这种平衡是不稳定的.只要有微扰(轻轻晃动)，就会破坏原来的平衡状态，直到达到另一种新的平衡态为止.这时重流体在下，

9、轻流体 在上，正好跟原来交换了位置，所以这种不稳定性也叫做 交换不稳定性.现在我们采用类比的方法来研究约束在磁场中的等离子体.假定磁场与等离子体之间达到了平衡，中间有明显的分界面.就是说在等离子体中没有磁场，在磁场中没有等离子体.这时，等离子体除了受到重力之外，还受到磁场的作用力，包括磁场梯度引起的力NVB和磁场的弯曲引起的力my2(b Sb.当然这是指单个粒子受到的力，我们把它们当作等效重力(跟 流体情况作类比)，记作VGeff , Geff= WB + m%b7) b(10)(b)b -B2mv _ W2B B以及粒子能量w=w,w代入上式并对整个麦克斯韦速度分布函数积分，我们可以得到作为

10、流体元的等效重力：B %B)PoGeff T P.+(11)I B B )对干各向同性等离子体，VB&VB, &也,因此叩、Geff - 2 P-B=. = 2Rc所以, Geff2PRe(12)Of将（12）式代入（9）式便得到描述瑞利一泰勒不稳定性的方程2P5A:。(13)因为在低P情况下上式说明，当磁场曲率Re与等离子体密度梯度VP。方向相反，即Re PP0 0.如图3.26 （a）所示.从图中可以看出，这图3.26 ,坏曲率才与“好曲率时的磁力线是凹向等离子体的.这种曲率被称为“ 坏曲率：图3. 26 （b）画出了稳定的磁场位形.此时，磁场曲率 Rc与等离子体压强梯度 （或密度梯度VP

11、。）同向.磁力线凸向等离子体，这种磁场位形的曲率被称为“好曲率在实际的磁场位形中，曲率矢量如往往不断改变方向.也就是说，在某个地方是“好曲率；在另一个地方则变成“坏曲率”.如在简单磁镜场中，在中心部位是“坏曲率”，而在“咽喉”部位则 是“好曲率”因此，有必要引入“平均曲率”的概念.定义： 磁力线管的比容U ,它是磁力线管的几何体积 每V与管内的磁通量。力的比值:6V = JSdl, g=B6S= const,V =、Sdl = B平均曲率的定义为匚建二-ndl 7B因此，平均曲率半径为前面得到的稳定条件（好曲率）是曲率与 VP同向，即VP Rc0,在聚变等离子体中，一般都是中心密度大，即PL

12、cP/cr 0;因此稳定条件要求Rc 0 .这就相当于要求 其中y表示比热比。设每一个变量均为平衡量和扰动量的叠加，即 f = fo + fi+.。为简化 起见，不考虑平衡流，即u0=0。(如果u0#0可以讨论)则将方程(1) (6)线性化之后 可得关于一阶扰动量的微分方程组d prdl 剂LL b剂FU ；：2V:二 0其中VW）为磁面包围的体积.因此，即 V仪）有极大值，其中必有磁场极小值，这相当于平均磁阱.这说明位于磁阱的等离子体是稳定的.与之相反，位于磁山“磁山”的等离子体是不稳定的,6.2等离子体的能量原理不考虑离子和电子的效应，可将等离子体作为单流体来处理。采用理想磁流体力学方程组

13、作为出发点dP v : u =0-:t:du = 一 p J B gd r p _000dt(D(2)(3)(4)丑E-:t(5)E u B =0(6)2 +口(PUi )=0(7)什：0如-R J1 B 0 J 0 B1(8)dt= - ;/p0, U1 -(U1 17 )p0(9)r1 = Du (u 黑 B0 )(10)ft.1,一一、J1 = - 1- B1(11)J 0令相对于流体元平衡位置r0的扰动位移=r-r。为一阶小量，则有她=一盯0)(；：t将上式分别代入方程(7)、(9)和(10),对时间积分，可将扰动密度、扰动压强和扰动磁场 均用扰动位移来表示(13)Pl =-?P0E-

14、( 3)P0(14)BlW (暧 B0)(15)将这些表达式代入方程(8),并利用方程(11),则可得到关于扰动位移 己的二阶微分方程.2 P。=F( 0(16)二 t1 .1 .F(9二P1 +正(寸x B1 /B。+7(X B0 y B1J 0J 0- II； P0V r c P0: B B B0J 0+(w B 0)mW (洛 B 0 )(17)显而易见，F (目相当于由扰动位移所引起的作用在单位流体体积上的力。在适当的边界条件 下解此方程，可以确定位形的稳定性。根据能量守恒原理，扰动位移引起的系统总能量的变化为零，即动能和势能的变化之和 为零-P0怛dr + BW = 02将上式对时间微商可得d:山图rW =0dt利用扰动方程（16）和函数F（E）的自伴性，即可将方程（19）中的第一项写成：。也= F dr:口2计dr1F己dr22 dt则由方程（19）可得扰动势能的变化1_、W = - 士 F 士 dr(18)(19)(20)(21)(22)从直观上来说，线性系统在力F（E）的方向上做位移士时，所做的功为3F（叮/2。根据能量守恒原理，这个功之可能是以消耗势能为代价，因此可得方程（22）。1-2;。+2JlflBoO B 於