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

被引:736
作者
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 条
[1]   A staggered mesh algorithm using high order Godunov fluxes to ensure solenoidal magnetic fields in magnetohydrodynamic simulations [J].
Balsara, DS ;
Spicer, DS .
JOURNAL OF COMPUTATIONAL PHYSICS, 1999, 149 (02) :270-292
[2]   Total variation diminishing scheme for adiabatic and isothermal magnetohydrodynamics [J].
Balsara, DS .
ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES, 1998, 116 (01) :133-153
[3]   Linearized formulation of the Riemann problem for adiabatic and isothermal magnetohydrodynamics [J].
Balsara, DS .
ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES, 1998, 116 (01) :119-131
[4]   Shock-capturing approach and nonevolutionary solutions in magnetohydrodynamics [J].
Barmin, AA ;
Kulikovskiy, AG ;
Pogorelov, NV .
JOURNAL OF COMPUTATIONAL PHYSICS, 1996, 126 (01) :77-90
[5]   A 2ND-ORDER PROJECTION METHOD FOR THE INCOMPRESSIBLE NAVIER STOKES EQUATIONS [J].
BELL, JB ;
COLELLA, P ;
GLAZ, HM .
JOURNAL OF COMPUTATIONAL PHYSICS, 1989, 85 (02) :257-283
[6]   FLUX-CORRECTED TRANSPORT .1. SHASTA, A FLUID TRANSPORT ALGORITHM THAT WORKS [J].
BORIS, JP ;
BOOK, DL .
JOURNAL OF COMPUTATIONAL PHYSICS, 1973, 11 (01) :38-69
[7]   THE EFFECT OF NONZERO-DEL.B ON THE NUMERICAL-SOLUTION OF THE MAGNETO-HYDRODYNAMIC EQUATIONS [J].
BRACKBILL, JU ;
BARNES, DC .
JOURNAL OF COMPUTATIONAL PHYSICS, 1980, 35 (03) :426-430
[8]   AN UPWIND DIFFERENCING SCHEME FOR THE EQUATIONS OF IDEAL MAGNETOHYDRODYNAMICS [J].
BRIO, M ;
WU, CC .
JOURNAL OF COMPUTATIONAL PHYSICS, 1988, 75 (02) :400-422
[9]   NUMERICAL SOLUTION OF NAVIER-STOKES EQUATIONS [J].
CHORIN, AJ .
MATHEMATICS OF COMPUTATION, 1968, 22 (104) :745-&
[10]   A simple finite difference scheme for multidimensional magnetohydrodynamical equations [J].
Dai, WL ;
Woodward, PR .
JOURNAL OF COMPUTATIONAL PHYSICS, 1998, 142 (02) :331-369