HOW-MHD: A High-Order WENO-Based Magnetohydrodynamic Code with a High-Order Constrained Transport Algorithm for Astrophysical Applications
arXiv:2304.04360 · doi:10.3847/1538-4357/acdf4b
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.
arXiv comment22 pages, 17 figures, submitted to ApJ
References in corpus (37)
- PLUTO: a Numerical Code for Computational Astrophysics
- Interstellar Turbulence I: Observations and Processes
- Athena: A New Code for Astrophysical MHD
- Density, Velocity, and Magnetic Field Structure in Turbulent Molecular Cloud Models
- Smoothed Particle Hydrodynamics and Magnetohydrodynamics
- The Statistics of Supersonic Isothermal Turbulence
- The Anisotropy of MHD Alfvénic Turbulence
- An Unsplit Godunov Method for Ideal MHD via Constrained Transport
- The Density Probability Distribution in Compressible Isothermal Turbulence: Solenoidal versus Compressive Forcing
- The Star Formation Rate of Supersonic MHD Turbulence
- 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
- Numerical Magnetohydrodynamics in Astrophysics: Algorithm and Tests for One-Dimensional Flow
- A Divergence-Free Upwind Code for Multidimensional Magnetohydrodynamic Flows
- High Order Upwind Schemes for Multidimensional Magnetohydrodynamics
- MHD simulations of jet acceleration from Keplerian accretion disks: the effects of disk resistivity
- From coronal observations to MHD simulations, the building blocks for 3D models of solar flares
- A Solution Accurate, Efficient and Stable Unsplit Staggered Mesh Scheme for Three Dimensional Magnetohydrodynamics
- Global Simulations of Accretion Disks I: Convergence and Comparisons with Local Models
- Numerical Magnetohydrodynamics in Astrophysics: Algorithm and Tests for Multi-Dimensional Flow
- Extended MHD modeling of the steady solar corona and the solar wind
- Cosmological MHD simulation of a cooling flow cluster
- Modeling the Large Scale Structures of Astrophysical Jets in the Magnetically Dominated Limit
- Vorticity, Shocks and Magnetic Fields in Subsonic, ICM-like Turbulence
- Finite Difference Weighted Essentially Non-Oscillatory Schemes with Constrained Transport for Ideal Magnetohydrodynamics
- A Multi-dimensional Code for Isothermal Magnetohydrodynamic Flows in Astrophysics
- Improved Accuracy of High-Order WENO Finite Volume Methods on Cartesian Grids
- A moving mesh unstaggered constrained transport scheme for magnetohydrodynamics
- Shock Waves and Energy Dissipation in Magnetohydrodynamic Turbulence
- Magnetohydrodynamic Simulation Code CANS+: Assessments and Applications
- 3D Simulations of MHD Jet Propagation Through Uniform and Stratified External Environments
- Systematic construction of upwind constrained transport schemes for MHD
- Turbulence Dynamo in Galaxy Clusters
- A high-order weighted finite difference scheme with a multi-state approximate Riemann solver for divergence-free magnetohydrodynamic simulations
- As A Matter of State: The role of thermodynamics in magnetohydrodynamic turbulence
- 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