等离子体物理学二ppt课件

1、等离子体物理学二李毅2021.10 等离子体中，电场、磁场、速度、密度、压力、温度等任何一个物理量 普通会随空间和时间变化。 扰动量原那么上它可以分解为各个平面波的叠加，即： 其中 为波的幅度，是物理量的Fourior分解： 对于其中恣意一支平面波来说，k为波矢，w为频率。 这里我们用复数来表示波是方便的，取其实部就是实践的值。等离子体中的线性波()( , )( ,)itkted d k xxkk()1( ,)( , )2itkt ed dt k xkxx( ,)kk 波的速度可以用相速度和群速度来描画。相速度是波在坚持相位不变的情况下的运动速度。相位为： 相位不变的条件下： 得到相速度：波的

2、相速度pdxvdtk()0ddkxtdtdttk x 波的群速度描画波包整体运动的速度，而波包是由满足一定色散关系的各种频率的波组成。假设该波包的色散关系为 只需频率满足此关系的波才存在，可以表示为： 因此由式积分，在波沿x方向直线传播情况下得：波包( ) )( , )( )i k xk tkx tk edk( ,)( )( )kkk kk( )k 假设波包的主要波数为k0，对应的频率 近似有： 其中群速度定义为： 代入可得： 可见波包的包络以群速度vg的速度前进。波的相速度可以超越光速。但群速度一定不能超越光速，由于群速度可以传送信息和能量，否那么会违背爱因斯坦的狭义相对论原理。 波的群速度

3、0000()()()( , )( )(, )gggik x v ti k vtki k vtgx tk edk exv t t e0gdvdk00()k00( )()gkvkk 一支波沿x方向传播，在y、z两个垂直方向上，电场矢量的分量Ey和Ez普通可以表示成： 其中， Ey0和Ez0 ，a，b均为常数。 在yz平面上的电场分量满足：波的旋转与偏振00cos(),cos()yyzzEEkxtEEkxt02222000cos()cos() sin()sin()cos()sin ()(1)zzyyzzyyEkxtkxtEEEEEEE 这阐明，电场矢量端点在yz平面内的轨迹是椭圆二次曲线中只需椭圆离

4、原点间隔有限，因此是椭圆偏振。 特殊情况下，可以是线偏振a=b或|a-b|=p，偏振方向与y轴夹角为 也可以是圆偏振Ey0=Ez0且|a-b|=p /2 。 当 a-b=p /2 时，例如a=0 而 b =-p /2 ，此时 随着波沿着x方向前进，相位添加，E矢量做右手旋转。所以是波是右旋的。波的旋转与偏振00cos(),sin()yyzzEEkxtEEkxt00arctan(/)yzEE 当 a-b=-p /2 时，例如，当 a=0 而 b =p /2 时： 随着波沿着x方向前进，E矢量按左手旋转。所以这时波是左旋的。 普通情况下，无妨取|a-b|p， 当a-b 0时，是右旋；而a-b0 时

5、，是左旋；a-b =0或p时，是线偏振。波的旋转与偏振00cos(),sin()yyzzEEkxtEEkxtxyz 将等离子体中的扰动作Fourior分解，也即化为多个平面波的线性叠加。假设方程组是线性的，对于一切满足方程组的平面波来说，其线性叠加也满足方程组。因此，从研讨最简单的平面波入手，我们就可以研讨扰动在等离子体中的传播和开展。方程组中的非线性项应该被忽略，这是由方程的线性特性所决议的。另外，非线性项都是二阶或二阶以上的小量，在解线性动摇问题时，可以忽略。波的线性化和平面波分解()( , )( ,)itkted d k xxkk 普通来说，对于等离子体中的动摇来说，其频率和波长有一定的

6、对应关系。或者说，对于一个给定的频率，只需对应波长的动摇才干存在。这种对应关系即为波的色散关系： 波的群速度的计算需求用到波的色散关系： 更重要的是有了色散关系，就知道了初始的扰动 在随后的开展变化：线性波的色散关系 ( ,)0Dk( )gdvkdk( ) )01( )( ),( , )( )2iitkkedted k xk xkkxxxkk0( )x 等离子体中，电子的运动会引起电荷分别，使得等离子体偏离电中性，从而产生静电场。在这个静电场的作用下，电子会改动运动形状，力图使等离子体恢复电中性，但是在等离子体恢复了电中性之后，电子依然具有一定的动能，其运动又会使等离子体产生非电中性。我们称电

7、子的这种振荡为电子静电波，也叫Langmuir波。 这种波维持了等离子体的准电中性。电子静电波 在冷等离子体中，这种动摇可以用一维方程组描画： 将方程组进展线性化和平面波分解，得到方程组：冷等离子体中的电子静电波00()0,()() ,eeeeeenn vtxe nndvEn mne Edtx1011100/eeei nikn vi m veEikEen 经过化简成为： 这阐明，假设要 ，即动摇存在，必需有 在电子热压力不可忽略的情况下，方程改写为 这里g为多方指数，而对于电子做1维运动的电子静电波情况，取 g=3。而对于普通电子做3维运动的情况，取我们熟知的 g=5/3。电子静电波的频率20

8、1201(1)0en eikEm0E 200peen em2(),0,eeeeeepdvdDn mne Ep ndtxdtD 方程组经过线性化和平面波分解，成为： 得到色散关系，也即这组方程存在非0解的条件为： 这里vse是电子的声波速度。热等离子体中的电子静电波10101010100/0eeeeeeei nikn vi m veEikpnikEenppnn 2222,epseseeTk vvm 在 的冷等离子体近似的条件下，回到冷等离子体时的电子静电波色散表达式，此时对应的Langmuir波的群速度为0，因此是不传播的局域震荡。而在热等离子体中，Langmuir波的群速度与电子热运动速度可达

9、同样的量级，类似于电子压力引起的纵波。 电子静电波的频率必需不小于电子等离子体频率 wpe，通常这是较高的频率。在这个频率下，离子由于其质量远大于电子质量，它来不及呼应这么高的频率变化。其运动可以忽略。 热电子静电波的讨论0eT 对于长波情况，色散关系可近似为 其群速度远小于电子的热速度vthe ： 对于短波情况，当 时， 群速度为 与电子热速度相当，这时会产生剧烈的波与电子的相互作用，需求用动力学才干加以研讨。 热电子静电波的讨论223(1)2pDek3()gDethethedvkvvdk1Dek2pe1.5gthevv 离子的运动可以产生频率较低的动摇。在研讨较低频率的等离子体动摇时，需求

10、同时思索电子和离子的运动其中，a代表等离子体中的一切粒子，即电子和各种离子 ：思索离子成分时的静电波10101001100()011()0nkvnn vntxikpdvpim vq En mn q EndtxEq nikEq nxdpnp ndtpn 将以上做过线性化和平面波分解之后的方程组再进展消元化简，得到色散方程 ： 其中wpa和vsa分别是a类粒子对应的等离子体振荡频率和声速。 由于离子质量远大于电子质量，那么 因此在高频时wwpe，色散关系公式中的求和的各项中，离子项远小于电子项，因此可以忽略。只保管电子项，此色散关系回到电子静电波的色散关系式。 思索离子时的静电波色散关系22220

11、222010,ppssnqTvk vmmpipe 思索低频情况为简化分析起见，无妨假设只需一种氢离子成份。 离子声波： 对于低频长波，klDe1，色散关系公式中的电子项和离子项均远大于1由于它们的分母均很接近于0，因此可以忽略第一项常数1，得到离子声波色散关系：离子声波,eeiissiTTkvvm 这很像在普通气体中传播的声波。由于波长很长，在这种长尺度条件下等离子体可以很好地坚持电中性，因此引起的扰动类似于中性气体中产生的紧缩波。但由于离子和电子必需坚持电中性，当离子运动时，电子必需跟随，两者牢牢地结合在一同。这时电子的压力影响也经过这种结合传送给离子，即使离子温度为0，由于有电子压力的存在

12、，也可以产生离子声波。现实上，在以后的动力论中我们知道，假设离子热运动速度与离子声波的速度相当的时候，会产生阻尼景象，离子声波不能存在，因此离子声波大多在Tivs情况vAvA情况下，可忽略位移电流，与MHD结果一致。1222(1)AAvkvc22 2(1)piik c平行传播的右旋圆偏振波 右旋圆偏振波。色散关系为： 偏振关系为右旋圆偏振： 频率较高时，右旋圆偏振波的色散关系又能写为 截止频率为也是垂直传播的X模的截止频率221(|)peekc 22222121()11(|)()ppepieikcn 0,0 xyzEiEE2224eeRpe平行传播的右旋圆偏振波讨论 接近电子盘旋频率、波长较短

13、时，是电子盘旋波： 在电子盘旋频率上共振。共振时，电子可继续从右旋圆偏振波中获得或失去能量。 频率极低时，成为右旋圆偏振Alfven波： 在cvA情况下，可忽略位移电流，与MHD结果一致。1222(1)AAvkvc22 2|(1)peek c 哨声波 右旋偏振波在频率远低于电子盘旋频率但又远高于离盘旋频率时，成为哨声波： 它的群速度为会随着频率升高： 当扰动发生时，高频成分的波群速度较快，会先被观测到，而低频成分的波随后才干被观测到。在地球外表，雷电引起的电磁脉冲扰动在电离层激发低频的哨声波，可以收到沿地球磁场传播到另一端，可以听到由高到低的类似哨声的信号。2 22|epek c222|2|e

14、gepepekcdcvdkAlfven波 线性的Alfven波是左旋和右旋圆偏振波的叠加： 反之，圆偏振波也能看成是两个相互垂直、相位差为90度、振幅一样的线偏振波的叠加。 左旋或右旋圆偏振Alfven波，线偏振Alfven波，它们都具有一样的色散关系。 Alfven波是沿磁场传播的左旋和右旋偏振波的低频长波极限情况。在频率趋于0时，波数k也趋于0，但与普通的截止情况不同，这不是通带与阻带的分界点。()()()000()()22i kzti kzti kztyxyxyEEEeieieEeeeee平行磁场和垂直磁场传播的波的色散关系图kcikcHHLHRLvARLvAkB0k | B0哨声波离子

15、盘旋波电子盘旋波电磁波X模O模AlfvenAlfven电磁波X模高混杂低混杂法拉第旋转 高频电磁波在等离子体中沿磁场传播时，左旋的波和右旋的波遵守不同的色散关系，他们的相速度也不同。假设初始时一个线偏振的电磁波，可以分解为左旋波和右旋波的叠加。设电磁波在z=0处为 ，那么在等离子体中传播之后，波场为 这两支波在等离子体中传播一样一段间隔之后，由于旋转的角度不同，从而两者叠加之后的线偏振波的偏振方向就会有所改动。0(0, )i txztE eEe()()00( , )()()22LRi k zti k ztxyxyEEz teieiEeeee法拉第旋转 其中，kL和kR分别是左旋波和右旋波的波矢

16、，记 ： 那么 可知角度j即为偏振方向转动的角度。而 ： 在磁场知的情况下，可以用丈量电磁波偏振方向的旋转角度来得到等离子体的密度。 ()0( , )(cossin)i k ztxyz tE e Eee(),22LRLRkkkk zk22222222222()1122(|)()(|)()(|)|2 ()()2 ()pepeLReieipeeipeeeiekkzzczzcc 法拉第旋转 当等离子体密度随空间变化时，可以丈量等离子体电子密度的线积分： 电离层电子密度的线积分是地球空间物理的一项重要的探测内容。220202()|zeeeemcn dze动实际动摇实际 动力论方程是描画等离子体分布函数

17、的变化的。对于处置热效应、波与粒子相互作用、多种速度成分的带电粒子等景象时，用磁流膂力学的描画显然是不全面的，这时应该选用动力论来处置。 对于动摇的问题，磁流体描画适宜于冷等离子体，且波与粒子相互作用较弱的情况，除此之外，用动力论方程来研讨动摇问题更加准确全面，且能得到一些流体的动摇实际中没有的结果。 对于空间等离子体，无碰撞的Vlasov动实际方程为：()0tqfffmxvvEvB一维静电动摇 思索一维静电波扰动，线性化之后为： 留意此时v仅仅是坐标变量，它不是1阶小量。经Fourier变换解得 ： 色散关系为： 1100txvqfvfEfm 101xEqf dv01()viqEffm kv

18、200( , )10()pvfDkdvnk kv 一维静电波色散方程 思索长波，k很小，相速度很大，展开： 取平衡时的分布函数为Maxwellian分布： 进而色散方程为：2023011()().0pvfkvkvkvdvnk20000exp(),22vmmvmvfnffTTT 223001()1()().0pm vkvkvkvfdvn kT积分展开项 因对称性，只需偶次项积分不为0： 运用定积分公式： 及 ，递推得 因此22340301.0pmkvk vfdvn kT222122111222xnnxnxnnnnIex dxxdexedxI0I2(21)!2mmmI22220(21)!10mmp

19、mmmmk Tm主值积分的结果 取头两项，并思索主要是电子的奉献，那么 ： 进一步近似可得 ： 这个结果和流体实际得到的色散关系完全一致。而朗道以为，以上运算过程中，积分存在奇点问题，即在速度等于波的相速度时，积分的分母为0。以上的处置方法只是主值积分，正确的计算需求沿着奇点下方的途径进展。22221(13)0peeek Tm2222(1 3)peDek朗道积分围道 假设按照朗道指出的途径积分，结果为： 这时，w不再为实数，而是含有虚部的复数。普通对于形如 其中虚部是小量，那么 这里下标r、i对应为实部和虚部。22202220/1(13)0pepeveeeevkfk Timn k( , )(

20、, )0riDkiDkIm(v)Re(v)积分围道(, )(, )irirrDkDk朗道阻尼率 运用到此处，可以得到电子静电波的阻尼率 ： 至于为什么要运用朗道围道进展积分，还需求从问题本来的物理过程看。2240023300/423233222/2exp()22213exp()822rpevepeerieereevkeeerreeepeDeDefmfn kn kTmmmkTTT kkk 静电动摇的物理过程 假设最初有扰动，可以对Vlasov方程进展时间t的Laplace变换求得以后的扰动电场，而不是做Fourier变换空间上依然做Fourier变换： 积分可得 代入电场方程1100()0ptt

21、xvqefvfEfdtm 01110( )( )( )0vtqE pfpfpfikvfpm0110( )1( )vtqE pffpfpikvm01110( )( )( )0vtqE pfpfpfikvfpm朗道围道的数学解释 可解出电场 ： 经Laplace反变换 这里s是充分大的实数，使得一切积分奇点都在积分线路的左边。而这些奇点使分母为0就是经过对应 之后的色散方程。而s是充分大的条件对应于色散方程对速度 v 的积分中，奇点 具有充分大的虚部，因此积分是从奇点下方经过的。 2110000( )() ()ptvvvfdvf dvqE pikpikvnpikv1( )( )2iptiE te

22、E p dpi pi /kip k朗道阻尼的物了解释 有趣的是，朗道阐明了要处理这个物理问题，积分围道需从下方绕过，才干满足数学上的要求。而一些数学家那么以为这种做法只是纯数学的东西，没有物理意义。直到后来实验和模拟都证明了朗道阻尼的存在，朗道的处置方法才被普遍的接受。 从物理上看，朗道阻尼其实是波与电子的共振相互作用。当电子的运动速度与波的相速度相差不大时，电子就被波的势阱捕获，从而与波一同运动。开场时速度小于波速的粒子得到加速，而开场时速度大于波速的粒子被减速，最后被捕获的粒子平均速度都与波的速度一样。朗道围道的物了解释 对于Maxwellian分布，运动速度在波速附近的粒子中，速度慢的比

23、速度快的粒子更多。从而获得加速的粒子多于减速的粒子。总体看来，波失去能量而粒子获得能量。波的幅度就会逐渐减小，构成阻尼。f(v)v通行粒子和波的势阱中捕获粒子被波加速的粒子绿色和被波减速的粒子蓝色磁化等离子体中动摇的动实际 设磁场沿z方向，而波矢为 ： 对于动力论方程 线性化之后 在速度空间的柱坐标系 中，有zzxkkkee0ffftm Fvxv011110()()()0ffqqiifmm k vEv Bv Bvv(, ,)zvv01111()()0ffqifm k vEv Bv磁化等离子体波中的分布函数 对于Maxwellian的初始分布 ： 代入得 解此微分方程： 其中，232020002

24、mvTfmmfnefTT vv11100()0fimffTB E vk v1010sin ,igigzzmfed e fBTk vk vgabab v E磁化等离子体波中的分布函数 利用Fourior展开式 ： 以及此式对 q 和 b 的偏导： 可积分得 其中，c是对q积分产生的常数，虽然它可以是 的函数，但由于对q 必需有 2p 的周期，应该取为0。sin( )ibinnneJb esinsincos( ),sin( )ibinibinnnnnneJ b eeiJ b eb()0110( )( )( )()i a nigz znxnynnf menfevJ bvJ bviJ bcBTi a

25、nbeeeE,zv v动实际伯恩斯坦波 利用Fourior展开式 ： 动实际伯恩斯坦波是磁化等离子体中的电子静电波。对于静电波，电场满足 ： 为求色散方程，做k空间的平面波分解，得 在对速度空间的q做一个周期的积分之后， 的项积分为0，对求和只需 一项不为0。()0110( ) ()( )( )()i n nnzzxnynnnf menfJ bvvJ bviJ bBTi a nb eeeE1110,qf dEEv2()0201( )( )()0i nnnnzznnq f denJb J b k vk vkTanb vnnnn动实际伯恩斯坦波 因此化简为： 利用定积分公式： 这里 In( ) 是

26、n阶修正Bessel函数。进一步得 其中 220202( )1()0znzznq f v dv dvJbk vnkTan 222()/401()()()22sxpqsnnnpqeJpx J qx xdxeIss22221( )02zmvpezzTznnmk vnmdv eIkTTan 22kTm垂直磁场传播的静电波 思索垂直磁场传播的静电波 ，此时， 色散方程为： 利用 的性质化为 在极限情况下， 用此调查色散关系的极限情况。0,zkkk22( )10pennmenIk Tan( )( )nnI xIx2222212( )10pennen In( )( /2) / !0( )/ 2nnnInI

27、e垂直磁场传播的静电波 因此，在 的极限下，有 在 的极限下，有 这是动实际伯恩斯坦波。22(),2 ,3 ,.HHpe ,2 ,3 ,. 0/khh/12345斜向传播的静电波 假设 ，色散方程中，积分会出现奇点。好像朗道阻尼的处置，需求用朗道围道做积分。即在满足共振条件的情况下，色散方程中会出现虚部。 化简 其中， 是等离子体色散函数。,( )2zmnZTk 221( ) ( ) 102pennzmmeIZk TkT22221( )02zmvpeTnznmmeIdvekTTan 0zk 斜向传播的静电波 等离子体色散函数的定义如下： 其中，s = 0, 1, 2 分别对应着积分围道取从奇点

28、之上对应于主值积分，穿过奇点，以及从奇点之下穿过朗道围道。 伯恩斯坦波没有阻尼，而对于斜向传播的静电波，可产生阻尼。221( ),0,1,2xeZdxisesx等离子体的平衡 平衡时，等离子体不运动，满足 这阐明，磁力线和电流线都是在等压面内。 是沿着等压面的法线方向，j和B都与它垂直，因此他们都平行于等压面 对于平直的磁力线，在垂直方向，有 对于柱等离子体，由于对称性，等压面就是柱的同心圆面。磁场既有轴向也有径向分量：pjB20costant2Bpp( )( )zzB rB rBeejB柱形等离子体的平衡 柱对称平衡时，磁力线具有一定的曲率： 其中 因此，径向的平衡方程为：200()02Bp

29、BB2200()02BdBpdrr211( )( )( )( )rB rB rBrr BBee无力场的平衡 在低b等离子体中，磁场力占主导位置，热压力梯度力可以忽略不计。平衡时，必需电流与磁场平行，满足： 其中a为常数时，是线性无力场，系统到达整体势能最小的平衡形状。 做旋度，得Helmholtz方程： 可分别解其中的三个分量。太阳低日冕中常用无力场模型。 当a不为常数时，是非线性无力场，求解更困难些。 BB22()0 B平衡系统的MHD不稳定性 MHD不稳定性能量原理： 假设系统中的等离子体各处都有微小位移，假设系统整体的能量有所上升，那么系统处于稳定平衡形状，反之，那么不平衡。如下图： 设

30、位移 速度可以表示为dddtttdttvv( , )( )i tte r r位移之后势能添加，系统稳定位移之后势能减小，系统不稳定扰动对系统的影响 延续性方程得到： 作为随体运动，从绝热方程得到 从磁场冻结方程得到 从牛顿方程得到 由于F与x成正比，系统能量的变化为10() BB10()100ppp 2010011p FjBjB220010011112212Wd dVd dVdVdVpdV F F jBjB系统能量的变化 其中： 留意到 ，即磁场与边境面平行，否那么等离子领会沿着磁力线运动，边境无法静止。因此101001100022011100110011()()()11()()()dVdVd

31、VdB dVddB dV jBBB BBB B BBS B BSB B S00nBe201101110201110110011()()22111()222nWdB dVppdVpdSB dVpdV BB SjB BBjB 扰动满足的边境条件 在边境两边总压力是一样的，即 思索到扰动时，其随体的变化为： 在等离子体和真空的边境，在跟随等离子体运动的坐标系中，电场切向分量延续，且都为0，即知边境附近的真空中 其中运用了磁场与边境面平行的边境条件。上式可简化为 并留意到真空的外表和等离子体的外表的矢量方向相反，且真空中没有电流，即 22010010100000()()22eeen nn nBBppB

32、 BBBee1100()()0neennett AeEvBeeB10nne eAB10e B220022eBBp金属壁包围等离子体的系统 运用 化简， 上式中，在完全导电的金属壁上，电场切向分量为零，因此在真空中的金属壁边境附近 因此可以不思索真空与金属壁边境面上的积分。11()00nnt AeeA11e BA0111111121111111neeneeevaceeeevacvacvacdSdSdddVdVB dV BBeABBASBASBABAAB系统能量的变化的各项 最后化简能量公式得 分析式中各项的作用，第一项总是正的，起稳定作用。等离子体形变时，引起真空区间的磁场变化，而真空中的磁场本

33、来处于能量最低形状，任何改动真空磁场形状的变化都要破费能量，因此这项抑制等离子体形变。222201001000001122222100100000100()222212()22221122neenenenevacBBBWpdSdVpdVBBBBdVpddVdVppdV BBejBSjB系统能量的变化各项的作用 第二项，假设从等离子体到真空时，总压力磁压加动压增大，那么这项为正，起稳定作用，反之那么是起不稳定作用。因此要想稳定，需求真空中的磁压力要越来越大。第三项，是等离子体内部的磁场起的稳定作用。第四项，经常引起不稳定，即电流驱动的不稳定性。第五项中的前一部分，经常引起压力不均匀导致的不稳定性

34、，后一部分总是正的，是等离子体的可紧缩性提供的稳定作用对于不可紧缩等离子体，这项为0。交换不稳定性 普通流体中的Rayleigh-Taylor不稳定性，在等离子体中表现为交换不稳定性。以磁场来平衡等离子体所遭到的非电磁力例如重力，就能够会出现这种不稳定性。 想象磁通一样的两束相邻的流管1与2，交换过后体积互变，体积改动的同时，压力变化满足绝热条件，以流管1为例，体积由V1变为V2，压力从P1变为P1，设 内能的表达式为2121,VVVPPP11000000111VVVPV VPV VPVUPdVPV V dV交换不稳定性 交换之后，内能变化为111 22 11 12 21 122 211 12

35、 211 11111121 1211221111111(1)(1)(1)1 (1)(1)1(1)()1 (1)12(1)(1)(1)2UPVPVPVPVPV VPV VPVPVPVVPVPVVPVPVPVVVVVPVVPVV 11112221 112211111112)(1)()(1)122PVP VPVP VPVVVP VP VVPP VVVP VP VV 交换不稳定性 由于上式中的第一项总是正的，因此不稳定性仅在于当第二项小于0时才有能够出现。思索等离子体边境处，等离子体通量管和它相邻的真空通量管交换时，总有dP0，因此稳定条件为dV0 ，即同样磁通量的真空通量管的体积比相邻的等离子体通量

36、管小，因此 式中，F是磁通量，B是磁场强度。00dldlVSdldlBBB Rayleigh-Taylor 不稳定性 思索有重力的情况，在垂直方向y方向，流体密度不均匀，且上面的流体密度大于下面的密度，可以引发Rayleigh-Taylor不稳定性。 设初始扰动： 且流体是不可紧缩的，即 思索系统扰动之后的受力变化为：()( , )(0,( ),( )i kztyztyye r0yzdikdy 1001000211000()( )pppppp F ggRayleigh-Taylor 不稳定性 将方程做旋度运算以便消去压力项，并留意到平衡时的0阶密度只与y有关，那么 利用不可紧缩条件消去得到 在

37、每一种介质里面，密度是常数，方程变为：2002000()()()yzyyikikg g222000()yyykkg 02000( )0kyyyykyeykyeyxykRayleigh-Taylor 不稳定性 在边境y=0附近区间对原方程积分得到 假设上面的密度大于下面的密度，那么得到虚的频率，引发不稳定性，其增长率为 对于有磁场的情况，留意到磁场散度为0，以及边境处磁场法线方向分量为0，因此B0y=0，磁场只需x和z分量。 22002()()ululluulkkgkg ()/()ululkgRayleigh-Taylor 不稳定性 旋度作用后的受力方程中添加磁场力一项为 而 而 代入下式： 2

38、220001()() yyyxkkgik jB1010000,yzyzzzyzzyyzBikBBikBBBB2011111222201100100() (/2 () ) () ) () )() )() ()xxxzyzzyzzyyzzyBikikB BB Bk B BikBk B j BBBBBBBBB100000()()()zyikBB BBBBRayleigh-Taylor 不稳定性 因此， 或 同样在分界面附近积分， 得到频率的表达式： 磁场项永远不小于0，在这里起稳定作用。直观地看，磁力线有张力而趋于平直，它与等离子体冻结，起固定等离子体的作用，抑制不稳定性。2222222000000

39、()( ()/yyyyzz ykkg kBk B 22220000()2/)()0ulzulk Bkkg 22002()()luululkg k B222222220000000(/) )(/)0zyzyyk Bkk Bkg 柱形等离子体不稳定性 假设等离子体分开平衡位置的位移可表示为： 且等离子体是不可紧缩的，即 思索柱形等离子体内部磁场和压力是均匀且电流存在于等离子体柱外表的情况， 内部磁场 外部磁场 受力方程为( , )( )( )( )ikz imi trrzztrrre reee0 00izBBe00( ) ,( )1/ezzBB rB rrBee200()()2dBpdtvBB柱形

40、等离子体不稳定性 记总压强 ： 保管一阶扰动小量时 而从磁感应方程及不可紧缩条件式可知 进而可得 取散度，由不可紧缩条件可得202Bpp2010120BptzB100()ikBB B2220010()k Bp 210p柱形等离子体不稳定性 化为Bessell方程 ： 可解得满足r=0处有限的解： 这里Im(x)是修正Bessell函数。对压力梯度式取方向分量，可得2221221()0ddmkpdrr drr11()( )()ikz imi tmmIkrpp aeIka2220010()() ( )( )()mrmk BkIkaap aIka 柱形等离子体不稳定性 在等离子体外部，由于没有电流，故磁场可以用势场表示： 由于磁场的散度总是为0，即有 同样满足与类似的Bessell方程，其解为在无穷远处为0，在这个边境条件下， 这里Km(x)是修正Bessell函数。因此， 1eB20()( )()ikz imi tmmKkraeKka1()( )( )()merzmkKkamaaiikKkaaBeee柱形等离子体不稳定性 等离子体内部和外部的扰动场如今都曾经知道，还需求用边境条件