论文部分内容阅读
本文的主要目的是研究求解磁流体力学方程的高精度高分辨率无振荡的高效数值方法。磁流体力学（MHD）方程的数值解法在天体物理、受控热核反应、雷达系统通信、发电系统以及流动控制等领域都具有广泛的应用。本文的工作主要包括两个方面：一、基于磁流体力学方程与双曲型守恒律的紧密联系，把求解双曲型守恒律的两类高效差分方法推广应用于求解磁流体力学方程；二、改进已有的一类求解磁流体力学方程的交错型中心差分格式。其主要内容包括以下几个方面：1.将求解双曲型守恒律的MmB（Maximum and minimum Bounded）差分格式推广应用于求解磁流体力学方程。基于通量分裂和单元平均的分片线性重构，通过适当选取数值导数，并利用Runge-Kutta TVD时间离散方法，得到了一类求解一维磁流体力学方程的二阶精度高分辨率无振荡的MmB格式，然后利用维数分裂（dimension-by-dimension）方法将格式推广到二维情形。最后，给出了一系列典型数值实验，验证了格式的有效性。2.将Kurganov和Levy提出的求解双曲型守恒律的三阶半离散CWENO（Central weighted essentially non-oscillatory）格式推广应用于求解磁流体力学方程。基于三阶CWENO重构，得到了一类求解一维磁流体力学方程的三阶半离散CWENO格式。然后利用维数分裂方法将格式推广到二维情形。最后，给出了一系列典型数值实验，验证了格式的有效性。3.将Balbas、Tadmor和Wu提出的求解一维和二维磁流体力学方程的交错型无振荡中心差分格式进行改进，得到了一类求解一维和二维磁流体力学方程的二阶和三阶非交错型高分辨率无振荡中心差分格式。最后，给出了一系列典型数值实验，验证了格式的有效性。4.对本文所得到的三类格式进行了比较，指出了它们各自的优缺点，并提出了有待今后进一步开展的工作。