无碰撞激波上游低频波的激发及其对激波耗散机制的影响[]陆全明郭.doc

无碰撞激波上游低频波的激发及其对激波耗散机制的影响* 陆全明 郭帆 王水 （中国科学技术大学地球和空间科学学院 合肥 230026） 摘 要 由于很大一部分来自激波上游的粒子被激波面所反射，因此在准平行无碰撞激波的上游存在着等离子体束流。本文通过一维混合模拟方法计算了束流密度较小（）和较大两种不同情形下的等离子体束流不稳定性。结果表明在束流密度较小时，束流激发的主要是平行于背景磁场方向传播的右旋波，此波动只能对束流粒子产生影响，包括减速和加热。在束流密度较大时，束流可同时激发平行和反平行于背景磁场方向传播的右旋波，它们除能对束流粒子产生影响外，还可通过非共振作用加速和加热背景粒子。它们对准平行无碰撞激波耗散机制的影响也进行了讨论。 关键词：束流不稳定性；波动；混合模拟；无碰撞激波 1．引言 等离子体束流不稳定性是两束不同流速的等离子体相互作用后激发等离子体波动的过程，在此过程中等离子体的集体作用起着主导作用。束流在等离子体中，尤其是空间等离子体中是一种普遍存在的物理现象[1-4]。本文主要讨论离子束流激发的低频波动，一个比较典型的例子是在无碰撞准平行激波的上游被激波反射的粒子与太阳风粒子形成的束流[5]。无碰撞激波的结构与激波上游的参数有着密切的关系，其中最重要的参数是激波上游等离子体沿激波法向的Mach数及磁场方向与激波法向的夹角。按照磁场方向与无碰撞激波法向的夹角，可将无碰撞激波分为准平行和准垂直两类：的激波为准垂直激波，而的激波为准平行激波[6]。在准平行无碰撞激波中，当激波上游的数足够大时，部分来自激波上游的粒子可被激波反射，由于其速度方向与磁场方向的夹角很小，粒子可以运动到激波面前几个甚至十个回旋半径的地方，从而激发起等离子体束流不稳定性[7-10]。对等离子体束流不稳定性的研究有助于我们更加深入地了解在激波上游观测到的大振幅低频波和激波的耗散机制。 在太阳风中,一般认为离子束流可激发低频电磁波，反过来低频电磁波可通过共振散射并加热束流离子。Akimoto利用一维混合模拟方法研究了离子束流激发低频电磁波的非线性过程[11]。Onsager[12]等发现束流激发的电磁波有两种：共振和非共振的右旋波，在束流密度很小时，主要激发共振模，当束流密度增大到一定程度时，非共振模也可以被激发起来，近期的研究表明它是分别是平行和反平行于背景磁场方向传播的磁声波[13]。当背景离子的垂直温度大于平行温度时，平行和反平行于背景磁场方向传播的Alfvén波也可被激发[14]。同时，Wu等最近发现离子束流激发的大振幅低频电磁波可加热背景离子[15]，而Chen等的研究也表明大振幅低频电磁波可通过非共振作用加热离子[16]。Lu和Wang的粒子模拟研究也表明离子束流激发的大振幅低频电磁波同样也可加热电子[17]。 在本文中，我们用一维混合模拟的方法研究离子束流激发的低频电磁波以及它们对背景离子的影响，在此基础上讨论了高马赫数情况下平行无碰撞激波的耗散机制。 2 数值方法 在本文中所采用的计算模型是一维混合模拟方法，它着重讨论离子的动力学过程，而将电子视为无质量的流体处理[18,19]。一维混合模拟方法中离子在四维相空间（）内运动，受Lorentz力作用，而电子满足无质量电子流体的动量方程。同时在计算中满足准中性条件，即在空间每一点电子和离子的电荷密度相等。联立求解每个质子的运动方程﹑宏观的电子动量方程和能量方程以及Maxwell方程组即可用混合模拟方法求解物理问题。粒子数密度以平均密度标定，磁场用初始磁场标定，时间尺度用质子回旋周期标定，，其中为质子质量，e为单位电荷。速度用初始时刻的Alfven速度标定，长度用标定。引入磁Lundquist数，定义为，其中c为光速，为波―粒子相互作用引起的反常电阻，在混合模拟中将其视为常数。计算区域的总长度，空间网格距取为，每个网格点中束流和背景的粒子数都是200个。背景和束流粒子不同数密度可通过对粒子取相应的权重而获得。取，。时间步长，边界条件采用周期边界条件。 3.计算结果 我们计算了在不同条件下的等离子体束流不稳定性，包括不同的束流与背景粒子的数密度之比及He2+粒子的影响，但保持束流相对于背景粒子的速度为。 计算中我们选取所选参考系为质心参考系，在此参考系中粒子的总电流电流为0，即： ﹑﹑分别为各成分粒子数密度﹑粒子流体速度和电荷数，i代表不同成分的粒子(包括束流粒子及背景粒子)。 一般来说，束流不稳定性可以激发两种类型的波动，分别是平行及反平行于背景磁场传播的右旋波。在束流密度很小时，主要激发平行背景磁场传播的右旋波；当束流密度增大到一定程度时，反平行于背景磁场传播的右旋波也可以被激发起来[11]。而通过Terasawa的方法可将平行及反平行于背景磁场