1.无碰撞激波的MHD模拟方法.doc

无碰撞激波的理想MHD模拟* 陆启明，杨维纮 （中国科学技术大学近代物理系，合肥230027） ［摘 要］ 本文提出了通过求解理想MHD方程模拟空间无碰撞激波的方法，并且使用该方法模拟了垂直无碰撞激波与行星际反向磁场结构和高密度等离子体团的相互作用过程，同时与粒子模拟的结果进行比对。两者的结果非常类似，而且在MHD模拟中得到了一些粒子模拟中没有观察到的现象。模拟结果表明了理想MHD模拟是准确且可行的，同时相对于粒子模拟又有很好的计算效率，便于扩展至高维的情形。 ［关键词］ 无碰撞激波　MHD方程 数值模拟 0. 引言 行星际空间和星际空间中充满着完全电离的稀薄等离子体，粒子平均自由程非常大，经典库仑碰撞效应往往可以略去。这些无碰撞的等离子体通常以超声速运动，形成太阳风和星风。当太阳风和星风遇到存在磁场的行星和恒星的阻挡时，在界面处将形成各种间断面，如地球的磁顶层和弓激波、太阳系的日球顶层等[1]。以弓激波为例，观察资料表明，弓激波是无碰撞的激波，上游是未扰动的超声速太阳风，而下游的等离子体以亚声速绕过地球的磁顶层[2]。无碰撞激波是等离子体物理、空间物理和天体物理学中的重要基础性课题，对它的深入研究有助于了解激波本身的产生、演化、耗散机制以及各种行星际结构与激波的相互作用问题。 同时，随着计算机技术的不断发展，使得通过计算机来模拟无碰撞激波成为可能。对无碰撞激波的数值研究通常可分