HOW-MHD: A High-Order WENO-Based Magnetohydrodynamic Code with a High-Order Constrained Transport Algorithm for Astrophysical Applications

Abstract

Due to the prevalence of magnetic fields in astrophysical environments, magnetohydrodynamic (MHD) simulation has become a basic tool for studying astrophysical fluid dynamics. To further advance the precision of MHD simulations, we have developed a new simulation code that solves ideal adiabatic or isothermal MHD equations with high-order accuracy. The code is based on the finite-difference weighted essentially non-oscillatory (WENO) scheme and the strong stability-preserving Runge-Kutta (SSPRK) method. Most of all, the code implements a newly developed, high-order constrained transport (CT) algorithm for the divergence-free constraint of magnetic fields, completing its high-order competence. In this paper, we present the version in Cartesian coordinates, which includes a fifth-order WENO and a fourth-order five-stage SSPRK, along with extensive tests. With the new CT algorithm, fifth-order accuracy is achieved in convergence tests involving the damping of MHD waves in three-dimensional space. And substantially improved results are obtained in magnetic loop advection and magnetic reconnection tests, indicating a reduction in numerical diffusivity. In addition, the reliability and robustness of the code, along with its high accuracy, are demonstrated through several tests involving shocks and complex flows. Furthermore, tests of turbulent flows reveal the advantages of high-order accuracy, and show the adiabatic and isothermal codes have similar accuracy. With its high-order accuracy, our new code would provide a valuable tool for studying a wide range of astrophysical phenomena that involve MHD processes.

22 pages, 17 figures, submitted to ApJ

References in corpus (12)

• PLUTO: a Numerical Code for Computational Astrophysics
• Athena: A New Code for Astrophysical MHD
• Smoothed Particle Hydrodynamics and Magnetohydrodynamics
• The Statistics of Supersonic Isothermal Turbulence
• The Density Probability Distribution in Compressible Isothermal Turbulence: Solenoidal versus Compressive Forcing
• An Unsplit Godunov Method for Ideal MHD via Constrained Transport in Three Dimensions
• ECHO: an Eulerian Conservative High Order scheme for general relativistic magnetohydrodynamics and magnetodynamics
• MHD simulations of jet acceleration from Keplerian accretion disks: the effects of disk resistivity
• Cosmological MHD simulation of a cooling flow cluster
• Modeling the Large Scale Structures of Astrophysical Jets in the Magnetically Dominated Limit
• WENO-Wombat: Scalable Fifth-Order Constrained-Transport Magnetohydrodynamics for Astrophysical Applications
• Effects of Forcing on Shocks and Energy Dissipation in Interstellar and Intracluster Turbulences

Cited by in corpus (3)

• Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows
• Physical Constraint Preserving Higher Order Finite Volume Schemes for Divergence-Free Astrophysical MHD and RMHD
• An Alternative Finite Difference WENO-like Scheme with Physical Constraint Preservation for Divergence-Preserving Hyperbolic Systems