A Comparison Study of Three CESE Schemes in MHD Simulation
JI Zhen1,2, ZHOU Yu-Fen1
1State Key Laboratory of Space Weather, Center for Space Science and Applied Research, Chinese Academy of Sciences, Beijing 100190 2Graduate School, Chinese Academy of Sciences, Beijing 100049
A Comparison Study of Three CESE Schemes in MHD Simulation
JI Zhen1,2, ZHOU Yu-Fen1
1State Key Laboratory of Space Weather, Center for Space Science and Applied Research, Chinese Academy of Sciences, Beijing 100190 2Graduate School, Chinese Academy of Sciences, Beijing 100049
The space-time conservation element and solution element (CESE) scheme is a new second order numerical scheme based on the concept of space-time conservation integration. In order to further overcome excessive numerical damping due to small Courant-Friedrichs-Lewy (CFL) number and to obtain a high quality solution, a Courant number insensitive (CNIS) scheme and a high-order scheme have been proposed by Chang et al. for fluid mechanics problems recently. In this study, to explore the potential capability of applications of the CNIS CESE scheme and the high-order CESE scheme to magnetohydrodynamics (MHD) equations, several benchmark MHD problems are calculated in one and two dimensions: (i) Brio and Wu's shock tube, (ii) Dai and Woodward's case, (iii) the Orszag-Tang vortex problem, (iv) the Riemann problem. The numerical results just prove that the CNIS scheme is more accurate and can keep the divergence free condition of the magnetic field, even if the CFL number is «1. Meanwhile, the tests show that the high order CESE scheme possesses the ability to solve MHD problems but is sensitive to the Courant number.
The space-time conservation element and solution element (CESE) scheme is a new second order numerical scheme based on the concept of space-time conservation integration. In order to further overcome excessive numerical damping due to small Courant-Friedrichs-Lewy (CFL) number and to obtain a high quality solution, a Courant number insensitive (CNIS) scheme and a high-order scheme have been proposed by Chang et al. for fluid mechanics problems recently. In this study, to explore the potential capability of applications of the CNIS CESE scheme and the high-order CESE scheme to magnetohydrodynamics (MHD) equations, several benchmark MHD problems are calculated in one and two dimensions: (i) Brio and Wu's shock tube, (ii) Dai and Woodward's case, (iii) the Orszag-Tang vortex problem, (iv) the Riemann problem. The numerical results just prove that the CNIS scheme is more accurate and can keep the divergence free condition of the magnetic field, even if the CFL number is «1. Meanwhile, the tests show that the high order CESE scheme possesses the ability to solve MHD problems but is sensitive to the Courant number.
JI Zhen;ZHOU Yu-Fen. A Comparison Study of Three CESE Schemes in MHD Simulation[J]. 中国物理快报, 2010, 27(8): 85201-085201.
JI Zhen, ZHOU Yu-Fen. A Comparison Study of Three CESE Schemes in MHD Simulation. Chin. Phys. Lett., 2010, 27(8): 85201-085201.
[1] Chang S C 1995 J. Comput. Phys. 119 295 [2] Chang S C, Wang X Y and Chow C Y 1998 NASA TM-1998-208843 [3] Chang S C, Wang X Y and Chow C Y 1999 J. Comput. Phys. 156 89 [4] Zhang Z C and Yu S T 1999 AIAA Paper 1999-0904 [5] Zhang Z C, Yu S T and Chang S C 2002 J. Comput. Phys. 175 168 [6] Feng X S, Hu Y Q and Wei F S 2006 Sol. Phys. 239 235 [7] Feng X S, Zhou Y F and Wu S T 2007 Astrophys. J. 655 1110 [8] Zhou Y F, Feng X S and Wu S T 2008 Chin. Phys. Lett. 25 790 [9] Chang S C 2002 AIAA Paper 2002-3890 [10] Yen Joseph C and Wagner Donald A 2005 AIAA Paper 2005-2820 [11] Chang S C 2007 AIAA Paper 2007-5820 [12] Chang S C 2008 NASA TM-2008-215138 [13] Liu K X and Wang J T 2004 Chin. Phys. Lett. 21 2085 [14] Zhang D L, Wang J T and Wang G 2009 Chin. J. Comput. Phys. 26 211 [15] Zhang M J, Henry Lin S C and John Yu S T 2002 AIAA Paper 2002-3888 [16] Jiang G S, Wu C C 1999 J. Comput. Phys. 150 561 [17] Zhang M J, John Yu S T and Henry Lin S C el at 2006 J. Comput. Phys. 214 599 [18] Xiong M, Peng Z, Hu Y Q, Zheng H N 2009 Chin. Phys. Lett. 26 015202 [19] Zheng H N, Zhang Y Y and Wang S et al 2006 Chin. Phys. Lett. 23 399