An Extension of the Athena++ Framework for General Equations of State

Abstract

We present modifications to the Athena++ framework to enable use of general equations of state (EOS). Part of our motivation for doing so is to model transient astrophysics phenomena, as these types of events are often not well approximated by an ideal gas. This necessitated changes to the Riemann solvers implemented in Athena++. We discuss the adjustments made to the HLLC, and HLLD solvers and EOS calls required for arbitrary EOS. We demonstrate the reliability of our code in a number of tests which utilize a relatively simple, but non-trivial EOS based on hydrogen ionization, appropriate for the transition from atomic to ionized hydrogen. Additionally, we perform tests using an electron-positron Helmholtz EOS, appropriate for regimes where nuclear statistical equilibrium is a good approximation. These new complex EOS tests overall show that our modifications to Athena++ accurately solve the Riemann problem with linear convergence and linear-wave tests with quadratic convergence. We provide our test solutions as a means to check the accuracy of other hydrodynamic codes. Our tests and additions to Athena++ will enable further research into (magneto)hydrodynamic problems where realistic treatments of the EOS are required.

arXiv comment32 pages, 12 figures, 18 tables, accepted by ApJS

References in corpus (16)

• The NumPy array: a structure for efficient numerical computation
• The Man Behind the Curtain: X-rays Drive the UV through NIR Variability in the 2013 AGN Outburst in NGC 2617
• First Results from the Catalina Real-time Transient Survey
• The Zwicky Transient Facility: Data Processing, Products, and Archive
• The All-Sky Automated Survey for Supernovae (ASAS-SN) Light Curve Server v1.0
• The disc instability model of dwarf-novae and low-mass X-ray binary transients
• CASTRO: A New Compressible Astrophysical Solver. I. Hydrodynamics and Self-Gravity
• Numerical Magnetohydrodynamics in Astrophysics: Algorithm and Tests for One-Dimensional Flow
• An Extension of the Athena++ Code Framework for GRMHD Based on Advanced Riemann Solvers and Staggered-Mesh Constrained Transport
• General-Relativistic Large-Eddy Simulations of Binary Neutron Star Mergers
• Convection Causes Enhanced Magnetic Turbulence in Accretion Disks in Outburst
• Dwarf Nova Outbursts with Magnetorotational Turbulence
• Impact of convection and resistivity on angular momentum transport in dwarf novae
• Convection Enhances Magnetic Turbulence in AM CVn Accretion Disks
• On the Piecewise Parabolic Method for Compressible Flow with Stellar Equations of State
• Solving the Riemann Problem for Realistic Astrophysical Fluids

Cited by in corpus (6)

• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• The Athena++ Adaptive Mesh Refinement Framework: Multigrid Solvers for Self-Gravity
• The Early Evolution of Magnetar Rotation I: Slowly Rotating "Normal" Magnetars
• Disc Novae: Thermodynamics of Gas Assisted Binary Black Hole Formation in AGN Discs
• Hydrodynamic simulations of black hole evolution in AGN discs I: orbital alignment of highly inclined satellites
• Two-Dimensional Radiation-Hydrodynamic Simulations of Luminous Red Novae