多维理想磁流体流动的数值模拟

摘要

采用修正的四步Runge—Kutta方法求解三维一般曲线坐标系下的理想磁流体方程组，为克服数值振荡，加特征型NND格式进行后处理。特征型NND格式推广到求解三维磁流体(MHD)问题需要知道雅可比矩阵的左右特征矩阵，在具体计算时需要克服矩阵的奇性。本文用三维程序采用推广的特征NND格式计算了多个典型的二维流动，计算结果表明，特征NND格式保持了TVD格式高精度的优点，又具有计算简单的特点，在包括强弱间断等复杂波系的非定常MHD流场数值模拟中是成功的。