磁流体数值模拟是空间物理研究的重要手段.采用具有TVD(Total Variation Diminishing)特性的LaxFriderichs差分格式求解了GLM-MHD(Generalized Lagrange Multiplier-Magnetohydrodynamics)方程.为降低格式的数值耗散,引入耗散修正系数对算法的通量计算过程进行改进.二维Rotor算例和磁云-电流片相互作用算例的模拟结果表明,GLM-MHD方法可以有效控制磁场散度误差,相对于泊松校正法可以节省一半以上的计算时间.在不破坏格式稳定性基础上,耗散修正系数有效降低了算法的数值耗散.
Magnetohydrodynamics(MHD) numerical simulation is an important tool for space physics research.In this paper,Lax-Eriderchs scheme with TVD property is employed to solve GLMMHD equations.The diffusion turning coefficient is introduced for scheme optimization.Simulation result of 2D rotor test and magnetic cloud current sheet interaction test demonstrates GLM-MHD method's divergence control capability.The simulation consumes less than half of the computational time comparing with simulation utilizing Poisson correction method.While numerical stability is not damaged,numerical diffusion is reduced by the diffusion tuning coefficient.