On the Divergence-free Condition and Conservation Laws in Numerical Simulations for Supersonic Magnetohydrodynamical Flows

An approach to maintain exactly the eight conservation laws and the divergence-free condition of magnetic fields is proposed for numerical simulations of multidimensional magnetohdyrodynamic (MHD) equations. The approach is simple and may be easily applied to both dimensionally split and unsplit Godunov schemes for supersonic MHD flows. The numerical schemes based on the approach are second-order accurate in both space and time if the original Godunov schemes are. As an example of such schemes, a scheme based on the approach and an approximate MHD Riemann solver is presented. The Riemann solver is simple and is used to approximately calculate the time-averaged flux. The correctness, accuracy, and robustness of the scheme are shown through numerical examples. A comparison in numerical solutions between the proposed scheme and a Godunov scheme without the divergence-free constraint implemented is presented.

[1]  Jean-Pierre Croisille,et al.  Numerical simulation of the MHD equations by a kinetic-type method , 1995 .

[2]  Paul R. Woodward,et al.  An approximate Riemann solver for ideal magnetohydrodynamics , 1994 .

[3]  Phillip Colella,et al.  A higher-order Godunov method for the equations of ideal magnetohydrodynamics , 1992 .

[4]  T. Jones,et al.  Cosmic-Ray Acceleration during the Impact of Shocks on Dense Clouds , 1993 .

[5]  J. Michael Picone,et al.  Evolution of the Orszag-Tang vortex system in a compressible medium , 1991 .

[6]  Paul R. Woodward,et al.  Extension of the Piecewise Parabolic Method to Multidimensional Ideal Magnetohydrodynamics , 1994 .

[7]  P. Colella.  Multidimensional upwind methods for hyperbolic conservation laws , 1990 .

[8]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[9]  M. Brio,et al.  An upwind differencing scheme for the equations of ideal magnetohydrodynamics , 1988 .

[10]  J. Hawley,et al.  Simulation of magnetohydrodynamic flows: A Constrained transport method , 1988 .

[11]  K. Zahnle,et al.  Explosion of comet Shoemaker-Levy 9 on entry into the Jovian atmosphere. , 1994, The Astrophysical journal.

[12]  Dongsu Ryu,et al.  Numerical magnetohydrodynamics in astrophysics: Algorithm and tests for multidimensional flow , 1995 .

[13]  J. M. Picone,et al.  Evolution of the Orszag-Tang vortex system in a compressible medium. I: Initial average subsonic flow , 1989 .

[14]  M. Norman,et al.  ZEUS-2D : a radiation magnetohydrodynamics code for astrophysical flows in two space dimensions. II : The magnetohydrodynamic algorithms and tests , 1992 .

[15]  J. Saltzman,et al.  An unsplit 3D upwind method for hyperbolic conservation laws , 1994 .

[16]  M. Norman,et al.  Shock interactions with magnetized interstellar clouds. 1: Steady shocks hitting nonradiative clouds , 1994 .

[17]  Paul R. Woodward,et al.  A simple Riemann solver and high-order Godunov schemes for hyperbolic systems of conservation laws , 1995 .

[18]  J. Brackbill,et al.  The Effect of Nonzero ∇ · B on the numerical solution of the magnetohydrodynamic equations☆ , 1980 .

[19]  S. Orszag,et al.  Small-scale structure of two-dimensional magnetohydrodynamic turbulence , 1979, Journal of Fluid Mechanics.

[20]  Phillip Colella,et al.  A Higher-Order Godunov Method for Multidimensional Ideal Magnetohydrodynamics , 1994, SIAM J. Sci. Comput..

[21]  Department of Materials Science,et al.  The Magnetohydrodynamics of Supersonic Gas Clouds: MHD Cosmic Bullets and Wind-swept Clumps , 1996, astro-ph/9606106.