Journal
ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES
Volume 149, Issue 2, Pages 447-455Publisher
UNIV CHICAGO PRESS
DOI: 10.1086/378771
Keywords
methods : numerical; MHD
Categories
Ask authors/readers for more resources
We describe a numerical method to solve the magnetohydrodynamic (MHD) equations. The fluid variables are updated along each direction using the flux-conservative, second-order, total variation diminishing (TVD), upwind scheme of Jin & Xin. The magnetic field is updated separately in two-dimensional advection-constraint steps. The electromotive force (EMF) is computed in the advection step using the TVD scheme, and this same EMF is used immediately in the constraint step in order to preserve del . B = 0 without the need to store intermediate fluxes. Operator splitting is used to extend the code to three dimensions, and Runge-Kutta is used to get second-order accuracy in time. The advantages of this code are high-resolution per grid cell, second-order accuracy in space and time, enforcement of the del . B = 0 constraint to machine precision, no memory overhead, speed, and simplicity. A three-dimensional FORTRAN implementation less than 400 lines long is made freely available. We also implemented a fully scalable message-passing parallel MPI version. We present tests of the code on MHD waves and shocks.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available