A Novel High-Order, Entropy Stable, 3D AMR MHD Solver with Guaranteed Positive Pressure

Abstract

We describe a high-order numerical magnetohydrodynamics (MHD) solver built upon a novel non-linear entropy stable numerical flux function that supports eight travelling wave solutions. By construction the solver conserves mass, momentum, and energy and is entropy stable. The method is designed to treat the divergence-free constraint on the magnetic field in a similar fashion to a hyperbolic divergence cleaning technique. The solver described herein is especially well-suited for flows involving strong discontinuities. Furthermore, we present a new formulation to guarantee positivity of the pressure. We present the underlying theory and implementation of the new solver into the multi-physics, multi-scale adaptive mesh refinement (AMR) simulation code $\text{FLASH}$ (http://flash.uchicago.edu). The accuracy, robustness and computational efficiency is demonstrated with a number of tests, including comparisons to available MHD implementations in $\text{FLASH}$.

arXiv comment42 pages

References in corpus (17)

• A Multi-Code Analysis Toolkit for Astrophysical Simulation Data
• The SILCC (SImulating the LifeCycle of molecular Clouds) project: I. Chemical evolution of the supernova-driven ISM
• Extensible Component Based Architecture for FLASH, A Massively Parallel, Multiphysics Simulation Code
• Kinetic energy preserving and entropy stable finite volume schemes for compressible Euler and Navier-Stokes equations
• Numerical Magnetohydrodynamics in Astrophysics: Algorithm and Tests for One-Dimensional Flow
• High Order Upwind Schemes for Multidimensional Magnetohydrodynamics
• A robust numerical scheme for highly compressible magnetohydrodynamics: Nonlinear stability, implementation and tests
• Modelling CO formation in the turbulent interstellar medium
• Modelling the supernova-driven ISM in different environments
• A Solution Accurate, Efficient and Stable Unsplit Staggered Mesh Scheme for Three Dimensional Magnetohydrodynamics
• A Direct Multigrid Poisson Solver for Oct-Tree Adaptive Meshes
• An Intercomparison Between Divergence-Cleaning and Staggered Mesh Formulations for Numerical Magnetohydrodynamics
• CosmoMHD: A Cosmological Magnetohydrodynamics Code
• Molecular cooling in the diffuse interstellar medium
• Relations between WENO3 and Third-order Limiting in Finite Volume Methods
• The Consistent Multi-fluid Advection method
• High-Order Discontinuous Galerkin Finite Element Methods with Globally Divergence-Free Constrained Transport for Ideal MHD

Cited by in corpus (19)

• Ideal GLM-MHD: About the entropy consistent nine-wave magnetic field divergence diminishing ideal magnetohydrodynamics equations
• From parallel to perpendicular -- On the orientation of magnetic fields in molecular clouds
• Comparison of some Entropy Conservative Numerical Fluxes for the Euler Equations
• A purely hyperbolic discontinuous Galerkin approach for self-gravitating gas dynamics
• High-order Magnetohydrodynamics for Astrophysics with an Adaptive Mesh Refinement Discontinuous Galerkin Scheme
• A Novel Averaging Technique for Discrete Entropy-Stable Dissipation Operators for Ideal MHD
• On the resolution requirements for modelling molecular gas formation in solar neighbourhood conditions
• Formulation of Entropy-Stable schemes for the multicomponent compressible Euler equations
• Entropy Symmetrization and High-Order Accurate Entropy Stable Numerical Schemes for Relativistic MHD Equations
• I. Jet Formation and Evolution due to 3D Magnetic Reconnection
• Unravelling the structure of magnetised molecular clouds with SILCC-Zoom: sheets, filaments and fragmentation
• Magnetic fields in star-forming systems (II): examining dust polarization, the Zeeman effect, and the Faraday rotation measure as magnetic field tracers
• Magnetic helicity dissipation and production in an ideal MHD code
• A Discontinuous Galerkin Solver in the FLASH Multi-Physics Framework
• Characteristic boundary conditions for magnetohydrodynamic equations
• Numerical Methods for Simulating Star Formation
• Compact Runge-Kutta flux reconstruction methods for non-conservative hyperbolic equations
• A Divergence-Free and $H(div)$-Conforming Embedded-Hybridized DG Method for the Incompressible Resistive MHD equations
• Linking simulations and observations of star-forming regions I. Masses and temperatures of ALMAGAL-like cores