一类TVD型组合差分方法及其在磁流体数值计算中的应用

根据太阳风数值模拟的特点,考虑到算法的质量(收敛速度、稳定性、精度等),结合磁流体数值计算的特性,对三维球坐标下磁流体动力学(MHD)方程组中的流体部分采用一种修正Lax-Friedrichs差分法而对磁场部分采用MacCormack格式,发展了一类快捷的具有TVD特性的组合数值新方法.作为格式的检验,在一维情况下,将其与PPM格式进行了比较,对一维快慢磁流体激波问题得到了与PPM格式精度相同的结果,然后将其应用到定态太阳风的数值模拟上,在不同等离子体β情形下,可得到理想的太阳风定态结构,为今后将此数值模式应用到具有复杂磁场位型或三维真实太阳风暴的数值模拟研究奠定了基础.