A CLASS OF TWO-STEP TVD MACCORMACK TYPE NUMERICAL SCHEME FOR MHD EQUATIONS  

作　　者：FENG Xueshang WEI Fengsi ZHONG Dingkun(Center for Space Science and Applied Research, The Chinese Academy of Sciences, Beijing 100080) 

出　　处：《空间科学学报》2003年第6期401-412,共12页Chinese Journal of Space Science

基　　金：the National Natural Science Foundation of China(No.49925412,49990450),the National Basic Research Science Foundation(No.G2000078405)

摘　　要：In this paper, a new numerical scheme of Total Variation Diminishing (TVD) MacCormack type for MagnetoHydroDynamic (MHD) equations is proposed by taking into account of the characteristics such as convergence, stability, resolution. This new scheme is established by solving the MHD equations with a TVD modified MacCormack scheme for the purpose of developing a scheme of quick convergence as well as of TVD property. To show the validation, simplicity and practicability of the scheme for modelling MHD problems, a self-similar Cauchy problem with the discontinuous initial data consisting of constant states, and the collision of two fast MHD shocks, and two-dimensional Orszag and Tang's MHD vortex problem are discussed with the numerical results conforming to the existing results obtained by the Roe type TVD, the high-order Godunov scheme,and Weighted Essentially Non-Oscillatory (WENO) scheme. The numerical tests show that this two-step TVD MacCormack numerical scheme for MHD system is of robust operation in the presence of very strong waves, thin shock fronts, thin contact and slip surface discontinuities.In this paper, a new numerical scheme of Total Variation Diminishing (TVD) Mac-Cormack type for MagnetoHydroDynamic (MHD) equations is proposed by taking into account of the characteristics such as convergence, stability, resolution. This new scheme is established by solving the MHD equations with a TVD modified MacCormack scheme for the purpose of developing a scheme of quick convergence as well as of TVD property. To show the validation, simplicity and practicability of the scheme for modelling MHD problems, a self-similar Cauchy problem with the discontinuous initial data consisting of constant states, and the collision of two fast MHD shocks, and two-dimensional Orszag and Tang's MHD vortex problem are discussed with the numerical results conforming to the existing results obtained by the Roe type TVD, the high-order Godunov scheme, and Weighted Essentially Non-Oscillatory (WENO) scheme. The numerical tests show that this two-step TVD MacCormack numerical scheme for MHD system is of robust operation in the presence of very strong waves, thin shock fronts, thin contact and slip surface discontinuities.

关 键 词：磁液体力学方程 两步TVD-MacCormack图解法 收敛性 稳定性 数值模拟 

分 类 号：O175.2[理学—数学]