论文部分内容阅读
多电流片较为广泛地存在于空间等离子体和实验等离子体中,如太阳风、日冕以及有反剪切磁位型的托卡马克中。它们往往可以发展产生双撕裂模不稳定性。在经历撕裂模各自独立发展的非线性演化阶段后,两个有理面上的磁岛会相互驱动发展从而爆发更加剧烈的磁场重联,磁场能量迅速转化成等离子体的动能,这种不稳定性严重影响着实验装置的稳定。而在高磁雷诺数下,重联区域的电流片由于受到两个有理面磁岛的相互驱动而变得细长。当电流片的横纵比大于一定值时,电流片变得不稳定从而产生二级磁岛。二级磁岛的产生影响着磁场的拓扑结构,同时也会加速磁场能量的释放。模拟研究高磁雷诺数下双撕裂模的演化,高精度的数值算法至关重要。我们从守恒磁流体方程出发,推导得到扰动形式的守恒磁流体方程,采用矢通量分裂(FVS)并结合加权本质无振荡格式(WENO)的有限差分方法,发展了高精度的磁流体数值模拟程序。通过此程序既可以模拟磁流体的线性发展过程,也可以用来研究磁流体的非线性演化。二维磁流体数值模拟程序计算量较大,为了使程序在计算机集群上高效并行运算,我们利用共享存储并行编程(OpenMP)和信息传递接口(MPI)的程序库设计了主从结构的并行程序。程序在节点间采用MPI进行通讯,在计算节点内部采用OpenMP进行并行化计算。同时程序会根据每个计算节点的计算能力对计算负载进行均衡化处理,从而降低各计算节点等待同步的时间,提高并行计算的效率。通过一系列的算例测试验证了程序的高精度以及模拟磁流体不稳定性的准确性,并将程序运行在计算集群多个节点上,检验了其并行加速的效果。剪切流在空间等离子体中和实验室等离子体中都普遍存在,对于磁场重联有着重要影响。我们在模拟中加入了反对称分布的极向剪切流,用数值模拟的方法探究了高磁雷诺数下剪切流对双撕裂模演化及二级磁岛发生过程的影响,研究发现:弱的极向剪切流对双撕裂模的线性发展和非线性发展都具有致稳的作用,且流速越大、流剪切越强,致稳效果越明显。当发生二级磁岛不稳定性时,剪切流会抑制磁场重联区域的入流驱动,使得形成的二级磁岛数目随着剪切流的大小以及剪切率的增加而下降。当剪切流到达一定强度时,二级磁岛不稳定性被完全抑制。因此在具有反剪切磁场位型的托卡马克中加入剪切流可以抑制双撕裂模不稳定性的发展。