thornado-hydro: a discontinuous Galerkin method for supernova hydrodynamics with nuclear equations of state

Abstract

This paper describes algorithms for non-relativistic hydrodynamics in the toolkit for high-order neutrino radiation hydrodynamics (thornado), which is being developed for multiphysics simulations of core-collapse supernovae (CCSNe) and related problems with Runge-Kutta discontinuous Galerkin (RKDG) methods. More specifically, thornado employs a spectral type nodal collocation approximation, and we have extended limiters - a slope limiter to prevent non-physical oscillations and a bound-enforcing limiter to prevent non-physical states - from the standard RKDG framework to be able to accommodate a tabulated nuclear equation of state (EoS). To demonstrate the efficacy of the algorithms with a nuclear EoS, we first present numerical results from basic test problems in idealized settings in one and two spatial dimensions, employing Cartesian, spherical-polar, and cylindrical coordinates. Then, we apply the RKDG method to the problem of adiabatic collapse, shock formation, and shock propagation in spherical symmetry, initiated with a 15 solar mass progenitor. We find that the extended limiters improve the fidelity and robustness of the RKDG method in idealized settings. The bound-enforcing limiter improves robustness of the RKDG method in the adiabatic collapse application, while we find that slope limiting in characteristic fields is vulnerable to structures in the EoS - more specifically, in the phase transition from nuclei and nucleons to bulk nuclear matter. The success of these applications marks an important step toward applying RKDG methods to more realistic CCSN simulations with thornado in the future.

References in corpus (15)

• Array Programming with NumPy
• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• GW190814: Gravitational Waves from the Coalescence of a 23 M${}_{\odot }$ Black Hole with a 2.6 M${}_{\odot }$ Compact Object
• Origin of the heavy elements in binary neutron-star mergers from a gravitational wave event
• Nucleosynthesis and Remnants in Massive Stars of Solar Metallicity
• Neutron-Rich Nuclei in Heaven and Earth
• The Role of Turbulence in Neutrino-Driven Core-Collapse Supernova Explosions
• Hydrodynamics of core-collapse supernovae and their progenitors
• A Second Relativistic Mean Field and Virial Equation of State for Astrophysical Simulations
• Three-dimensional Boltzmann-Hydro code for core-collapse in massive stars I. special relativistic treatments
• Discontinuous Galerkin methods for general-relativistic hydrodynamics: formulation and application to spherically symmetric spacetimes
• A 3D Simulation of a Neutrino-Driven Supernova Explosion Aided By Convection and Magnetic Fields
• Bound-Preserving Discontinuous Galerkin Methods for Conservative Phase Space Advection in Curvilinear Coordinates
• On the Piecewise Parabolic Method for Compressible Flow with Stellar Equations of State
• On the character of turbulence in self-consistent models of core-collapse supernovae