Splitting Based Scheme for Three-dimensional MHD with Dual Time Stepping  

作　　者：FU Huazheng FENG Xueshang 

机构地区：[1]State Key Laboratory for Space Weather,Center for Space Science and Applied Research,Chinese Academy of Sciences [2]University of Chinese Academy of Sciences

出　　处：《空间科学学报》2015年第1期9-17,共9页Chinese Journal of Space Science

基　　金：Supported by the National Basic Research Program of China(2012CB825601);the National Natural Science Foundationof China(41031066,41231068,41274192,41074121,41204127);the Knowledge Innovation Program of the ChineseAcademy of Sciences(KZZD-EW-01-4);the Specialized Research Fund for State Key Laboratories

摘　　要：A new hybrid numerical scheme of combining an E-CUSP(Energy-Convective Upwind and Split Pressure) method for the fluid part and the Constrained Transport(CT) for the magnetic induction part is proposed.In order to avoid the occurrence of negative pressure in the reconstructed profiles and its updated value,a positivity preserving method is provided.Furthermore,the MHD equations are solved at each physical time step by advancing in pseudo time.The use of dual time stepping is beneficial in the computation since the use of dual time stepping allows the physical time step not to be limited by the corresponding values in the smallest cell and to be selected based on the numerical accuracy criterion.This newly established hybrid scheme combined with positivity preserving method and dual time technique has demonstrated the accurateness and robustness through numerical experiments of benchmark problems such as the 2D Orszag-Tang vortex problem and the3 D shock-cloud interaction problem.A new hybrid numerical scheme of combining an E-CUSP(Energy-Convective Upwind and Split Pressure) method for the fluid part and the Constrained Transport(CT) for the magnetic induction part is proposed.In order to avoid the occurrence of negative pressure in the reconstructed profiles and its updated value,a positivity preserving method is provided.Furthermore,the MHD equations are solved at each physical time step by advancing in pseudo time.The use of dual time stepping is beneficial in the computation since the use of dual time stepping allows the physical time step not to be limited by the corresponding values in the smallest cell and to be selected based on the numerical accuracy criterion.This newly established hybrid scheme combined with positivity preserving method and dual time technique has demonstrated the accurateness and robustness through numerical experiments of benchmark problems such as the 2D Orszag-Tang vortex problem and the3 D shock-cloud interaction problem.

关 键 词：SPLITTING BASED SCHEME THREE-DIMENSIONAL MHD Dual time STEPPING 

分 类 号：P353[天文地球—空间物理学]