DISPATCH methods: an approximate, entropy-based Riemann solver for ideal magnetohydrodynamics

Abstract

With the advance of supercomputers we can now afford simulations with very large ranges of scales. In astrophysical applications, e.g. simulating Solar, stellar and planetary atmospheres, interstellar medium, etc; physical quantities, like gas pressure, density, temperature, plasma $\beta$, Mach, Reynolds numbers can vary by orders of magnitude. This requires a robust solver, which can deal with a very wide range of conditions and be able to maintain hydrostatic equilibrium where it is applicable. We reformulate a Godunov-type HLLD Riemann solver that it would be suitable to maintain hydrostatic equilibrium in atmospheric applications in a range of Mach numbers, regimes where kinetic and magnetic energies dominate over thermal energy without any ad-hoc corrections. We change the solver to use entropy instead of total energy as the primary thermodynamic variable in the system of MHD equations. The entropy is not conserved, it increases when kinetic and magnetic energy is converted to heat, as it should. We propose using an approximate entropy - based Riemann solver as an alternative to already widely used Riemann solver formulations where it might be beneficial. We conduct a series of standard tests with varying conditions and show that the new formulation for the Godunov type Riemann solver works and is promising.

13 pages, 17 figures, 3 page appendix with numerical method; accepted in A&A

References in corpus (18)

• GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
• Athena: A New Code for Astrophysical MHD
• An Unsplit Godunov Method for Ideal MHD via Constrained Transport
• A High Order Godunov Scheme with Constrained Transport and Adaptive Mesh Refinement for Astrophysical MHD
• High Order Upwind Schemes for Multidimensional Magnetohydrodynamics
• The Magnetic Rayleigh-Taylor Instability in Three Dimensions
• A Well-Posed Kelvin-Helmholtz Instability Test and Comparison
• rpSPH: a novel Smoothed Particle Hydrodynamics Algorithm
• Generation solar-like differential rotation
• Pebble dynamics and accretion onto rocky planets. I. Adiabatic and convective models
• High-order Magnetohydrodynamics for Astrophysics with an Adaptive Mesh Refinement Discontinuous Galerkin Scheme
• DISPATCH: A Numerical Simulation Framework for the Exa-scale Era. I. Fundamentals
• A Low Mach Number Solver: Enhancing Applicability
• Pebble dynamics and accretion onto rocky planets. II. Radiative models
• Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers
• AZEuS: An Adaptive Zone Eulerian Scheme for Computational MHD
• Curvilinear Grids for WENO Methods in Astrophysical Simulations
• Toward fully compressible numerical simulations of stellar magneto-convection with the RAMSES code