第五章 太阳风.pptVIP

2002-3-28 《等离子体天体物理》课程讲义(3) 第五章 太阳风 5.1 引言 5.2 Parker模型 5.3 磁场效应 第六章 磁重联 6.1 引言 6.2 磁湮灭（Sonnerup Priest 1975） 6.3 磁重联的定性效应 6.4 电流片的形成 6.5 线性磁重联 6.6 快速稳态磁重联 6. 7 磁重联在天体物理中的应用 将这个结果叠加到无穷远处的场（Be）中，可以得到： 入流区域由稍微弯曲的磁力线、均匀水平磁场（Bex0）、和一个势场三部分构成。 该势场是在相当远处为零、扩散区为零、而在激波阵面上具有法向分量BN的Laplace方程的解。 作为近似，忽略激波的弯曲，于是问题成为寻找上半平面的解。如图所示。 这个磁场应该在无穷远处为零，在x-轴的Le与L之间为2BN，相应地，在负x-轴的-Le与-L之间为-2BN。采用§3.2节介绍的势场计算方法，可得在原点处的磁场为： 但是，在激波阵面上，根据(4.37)式，慢激波是以法向磁场对应的阿尔芬速度运动, 则: 请更正 其中ve=BN / ?(??)，Me=ve / vAe 。上面这个Bi 就是我们所寻找的结果。 因为Me 1，而 Bi ?Be ，于是尺度律(6.38)和(6.39)成为： 因为log(Rme)变化很慢，上述重联率的典型数值是0.01。可见对于Rme的典型值，该重联率远大于Sweet-Parker模型的结果。 在Petschek重联模型提出之后的20多年里，人们以为快速重联问题已被彻底解决。直到1986年，根据高分辨率的数值试验，一类新的重联模型被发现了，而Petschek模型只是其中一个特例。 上式表明，当磁雷诺数（Rme）或者重联率（Me）增加时，中心区域的尺寸要下降。Petschek指出，当Bi 变得太小时，这个机制会阻止重联的进行。于是，Petschek在(6.42)中设Bi=1/2Be，估计了一个最大重联率（Me*） (c)准均匀场快速磁重联的统一理论（Priest Forbes, 1986） 数值试验表明，常常会出现与Petschek理论不一致的结果。例如，会出现发散而不是收敛的入流；扩散区域大小随重联率变化会产生不同的尺度律。 　 因此，目标是分析入流区域，确定 Bi/Be 的比值，进而确定 Mi 随重联率 Me 变化的规律。 考虑二维稳态理想MHD方程： 选择如下的解形式，其中Me1 是展开因子，使解为均匀场（Be=Bex0）的小扰动： 其中： ，而E为常数。 忽略压力梯度项，则（6.45）的最低阶近似给出： 如右图所示，边界条件为顶部B1x=0，侧边?B1y/?x=0，而在底部x轴上： 或若 在上述条件下，可得到分离变量解为： 其中： 从上述表达式我们可以计算Bi/Be并由(6.37)带入(6.42)，最终得到根据给定的Rme值而画出的Me随Mi的变化曲线，如图(6.10)所示。 由图可见，对于给定的Rme值，确实存在一个最大重联率（ Me*）； 另外， Me*随Rme值的变化规律（b=0曲线）确实靠近Petschek的估计（虚线）。 参数b的含义如后文所述。 Priest教授当时是加进了压力项：即在方程（6.47）和（6.48）的右端分别增加了dp1/dy和-(?/Be)dp1/dy。但对解的影响只是在（6.49）的和式里每项增加了个常数项：anb。 注意培养良好的科研习惯：经常以自己的方式重新推导别人的结果，使其成为自己的东西，并深入理解其前提。所谓 融会贯通。 新参数b的变化导致了不同的重联机制： （i） b=0 为Petschek模型（弱快激波膨胀）； （ii） b=1 为类Sonnerup模型（弱慢激波膨胀）； （iii）b0 为慢波模压缩型； （iv）0b1 为快、慢波杂交模膨胀型； （v） b1 为慢波模的通量积聚型。 如右图所示。 在这里： 压缩是指当等离子体进入扩散区时，压力增加； 膨胀是指当等离子体进入扩散区时，压力减少； 慢波模的作用是使压力和磁场具有不同的变化趋势； 快波模的作用是使压力和磁场具有相同的变化趋势。 (d)数值实验 Biskamp (1986) 完成了一系列引人注目的数值实验来研究快速稳态磁重联。 他发现当Rme(或Me)增加时，扩散区的长度和宽度也增加，与Petschek机制相反，于是得出结论：对于大的Rme值，Petschek机制及快速磁重联不存在。 但是Priest Forbes (1992)用新的模型通过进一步工作表明， Biskamp的特殊尺度律是由于边