The del•B=0 constraint in shock-capturing magnetohydrodynamics codes

被引:751
作者
Tóth, G [1 ]
机构
[1] Dept Atom Phys, H-1117 Budapest, Hungary
关键词
numerical approximation; stability and convergence of difference methods; magnetohydrodynamics;
D O I
10.1006/jcph.2000.6519
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Seven schemes to maintain the del . B=0 constraint numerically are compared. All these algorithms can be combined with shock-capturing Godunov type base schemes. They fall into three categories: the eight-wave formulation maintains the constraint to truncation error, the projection scheme enforces the constraint in some discretization by projecting the magnetic field, while the five different versions of the constrained transport/central difference type schemes conserve del . B to machine accuracy in some discretization for every grid cell. It is shown that the three constrained transport algorithms, which have been introduced recently, can be recast into pure finite volume schemes, and the staggered representation of the magnetic field is unnecessary. Another two new and simple central difference based algorithms are introduced. The properties of the projection scheme are discussed in some detail, and I prove that it has the same order of accuracy as the base scheme even for discontinuous solutions. I describe a flexible and efficient implementation of the projection scheme using conjugate gradient type iterative methods. Generalizations to resistive MHD, to axial symmetry, and to non-Cartesian grids are given for all schemes. The theoretical discussion is followed by numerical tests, where the robustness, accuracy, and efficiency of the seven schemes and the base scheme can be directly compared. All simulations are done with the Versatile Advection Code, in which several shock-capturing base schemes are implemented. Although the eight-wave formulation usually works correctly, one of the numerical tests demonstrates that its non-conservative nature can occasionally produce incorrect jumps across strong discontinuities. Based on a large number of tests, the projection scheme, one of the new central difference based schemes, and one of the constrained transport schemes are found to be the most accurate and reliable among the examined methods. (C) 2000 Academic Press.
引用
收藏
页码:605 / 652
页数:48
相关论文
共 52 条
[31]   Transport of magnetic flux in an arbitrary coordinate ALE code [J].
Peterkin, RE ;
Frese, MH ;
Sovinec, CR .
JOURNAL OF COMPUTATIONAL PHYSICS, 1998, 140 (01) :148-171
[32]   Nonlinear mhd simulations of wave dissipation in flux tubes [J].
Poedts, S ;
Toth, G ;
Belien, AJC ;
Goedbloed, JP .
SOLAR PHYSICS, 1997, 172 (1-2) :45-52
[33]  
POWELL KG, 1994, 9424 ICASE
[34]   Notes on the eigensystem of magnetohydrodynamics [J].
Roe, PL ;
Balsara, DS .
SIAM JOURNAL ON APPLIED MATHEMATICS, 1996, 56 (01) :57-67
[35]  
ROE PL, 1981, J COMPUT PHYS, V43, P357, DOI [10.1016/0021-9991(81)90128-5, 10.1006/jcph.1997.5705]
[36]   A divergence-free upwind code for multidimensional magnetohydrodynamic flows [J].
Ryu, DS ;
Miniati, F ;
Jones, TW ;
Frank, A .
ASTROPHYSICAL JOURNAL, 1998, 509 (01) :244-255
[37]   NUMERICAL MAGNETOHYDRODYNAMICS IN ASTROPHYSICS - ALGORITHM AND TESTS FOR MULTIDIMENSIONAL FLOW [J].
RYU, DS ;
JONES, TW ;
FRANK, A .
ASTROPHYSICAL JOURNAL, 1995, 452 (02) :785-+
[38]  
Schutgens NAJ, 1999, ASTRON ASTROPHYS, V345, P1038
[39]   ZEUS-2D - A RADIATION MAGNETOHYDRODYNAMICS CODE FOR ASTROPHYSICAL FLOWS IN 2 SPACE DIMENSIONS .2. THE MAGNETOHYDRODYNAMIC ALGORITHMS AND TESTS [J].
STONE, JM ;
NORMAN, ML .
ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES, 1992, 80 (02) :791-818
[40]   ON CONSTRUCTION AND COMPARISON OF DIFFERENCE SCHEMES [J].
STRANG, G .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1968, 5 (03) :506-+