A Free, Fast, Simple and Efficient TVD MHD Code

Abstract

We describe a numerical method to solve the magnetohydrodynamic (MHD) equations. The fluid variables are updated along each direction using the flux conservative, 2nd order, total variation diminishing (TVD), upwind scheme of Jin and 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 \grad.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 \grad.B=0 constraint to machine precision, no memory overhead, speed, and simplicity. A 3-D 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.

17 pages, 6 figures, submitted to the astrophysical journal

References in corpus (2)

• Local Three-Dimensional Simulations of Magneto-Rotational Instability in Radiation-Dominated Accretion Disks
• A Primer on Eulerian Computational Fluid Dynamics for Astrophysics

Cited by in corpus (33)

• Athena: A New Code for Astrophysical MHD
• An Unsplit Godunov Method for Ideal MHD via Constrained Transport
• An Unsplit Godunov Method for Ideal MHD via Constrained Transport in Three Dimensions
• Axisymmetric simulations of magneto--rotational core collapse: dynamics and gravitational wave signal
• High Performance P3M N-body code: CUBEP3M
• Global galactic dynamo driven by cosmic-rays and exploding magnetized stars
• GAMER: a GPU-Accelerated Adaptive Mesh Refinement Code for Astrophysics
• Cosmic rays can drive strong outflows from gas-rich high-redshift disk galaxies
• The influence of model parameters on the prediction of gravitational wave signals from stellar core collapse
• Simulations of magnetic fields in filaments
• Modeling Bell's Non-resonant Cosmic Ray Instability
• Gravitational waves from 3D MHD core collapse simulations
• The Fate of Non-Radiative Magnetized Accretion Flows: Magnetically-Frustrated Convection
• Numerical Parameter Survey of Nonradiative Black Hole Accretion -- Flow Structure and Variability of the Rotation Measure
• Three-dimensional hydrodynamical models of wind and outburst-related accretion in symbiotic systems
• Implementation of CR Energy SPectrum (CRESP) algorithm in PIERNIK MHD code. I. Spectrally resolved propagation of CR electrons on Eulerian grids
• FISH: A 3D parallel MHD code for astrophysical applications
• Numerical simulations of buoyancy instabilities in galaxy cluster plasmas with cosmic rays and anisotropic thermal conduction
• PIERNIK mhd code - a multi-fluid, non-ideal extension of the relaxing-TVD scheme (I)
• PIERNIK mhd code - a multi-fluid, non-ideal extension of the relaxing-TVD scheme (II)
• Stellar evolution of massive stars at very low metallicities
• The angular momentum structure of cosmic ray driven galactic outflows triggered by stream accretion
• Magnetohydrodynamics on Heterogeneous architectures: a performance comparison
• The phase structure of cosmic ray driven outflows in stream fed disc galaxies
• PIERNIK mhd code - a multi-fluid, non-ideal extension of the relaxing-TVD scheme (III)
• Detecting and reconstructing gravitational waves from the next Galactic core-collapse supernova in the Advanced Detector Era
• PIERNIK mhd code - a multi-fluid, non-ideal extension of the relaxing-TVD scheme (IV)
• Cosmic-ray driven dynamo in galactic disks
• Approaching the dynamics of hot nucleons in supernovae
• Fast magnetic reconnection in three dimensional MHD simulations
• MHD Generation Code
• Global simulations of galactic dynamo driven by cosmic-rays and exploding magnetized stars
• Cosmic-ray driven dynamo in galaxies