A Discontinuous Galerkin Method for General Relativistic Hydrodynamics in thornado

Abstract

Discontinuous Galerkin (DG) methods provide a means to obtain high-order accurate solutions in regions of smooth fluid flow while, with the aid of limiters, still resolving strong shocks. These and other properties make DG methods attractive for solving problems involving hydrodynamics; e.g., the core-collapse supernova problem. With that in mind we are developing a DG solver for the general relativistic, ideal hydrodynamics equations under a 3+1 decomposition of spacetime, assuming a conformally-flat approximation to general relativity. With the aid of limiters we verify the accuracy and robustness of our code with several difficult test-problems: a special relativistic Kelvin--Helmholtz instability problem, a two-dimensional special relativistic Riemann problem, and a one- and two-dimensional general relativistic standing accretion shock (SAS) problem. We find good agreement with published results, where available. We also establish sufficient resolution for the 1D SAS problem and find encouraging results regarding the standing accretion shock instability (SASI) in 2D.

arXiv comment14 pages, 6 figures

References in corpus (25)

• A Multi-Code Analysis Toolkit for Astrophysical Simulation Data
• Stability of Standing Accretion Shocks, With an Eye Toward Core Collapse Supernovae
• Radiation hydrodynamics with neutrinos: Variable Eddington factor method for core-collapse supernova simulations
• A Finite Difference Representation of Neutrino Radiation Hydrodynamics in Spherically Symmetric General Relativistic Space-Time
• Truncated Moment Formalism for Radiation Hydrodynamics in Numerical Relativity
• An HLLC Riemann Solver for Relativistic Flows: I. Hydrodynamics
• An Open-Source Neutrino Radiation Hydrodynamics Code for Core-Collapse Supernovae
• The Status of Multi-Dimensional Core-Collapse Supernova Models
• THC: a new high-order finite-difference high-resolution shock-capturing code for special-relativistic hydrodynamics
• The Spherical Accretion Shock Instability in the Linear Regime
• An efficient shock-capturing central-type scheme for multidimensional relativistic flows. I. Hydrodynamics
• A new multi-dimensional general relativistic neutrino hydrodynamics code for core-collapse supernovae. I. Method and code tests in spherical symmetry
• Relativistic simulations of rotational core collapse. I. Methods, initial models, and code tests
• Fornax: a Flexible Code for Multiphysics Astrophysical Simulations
• A new multidimensional, energy-dependent two-moment transport code for neutrino-hydrodynamics
• Non-linear axisymmetric pulsations of rotating relativistic stars in the conformal flatness approximation
• 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 New Spherical Harmonics Scheme for Multi-Dimensional Radiation Transport I: Static Matter Configurations
• Conservative 3+1 General Relativistic Variable Eddington Tensor Radiation Transport Equations
• Impact of Neutrino Opacities on Core-Collapse Supernova Simulations
• Formulation of discontinuous Galerkin methods for relativistic astrophysics
• Realizability-Preserving DG-IMEX Method for the Two-Moment Model of Fermion Transport
• Discontinuous Galerkin finite element methods for radiative transfer in spherical symmetry
• Computational Models of Stellar Collapse and Core-Collapse Supernovae