CosmosDG: An hp-Adaptive Discontinuous Galerkin Code for Hyper-Resolved Relativistic MHD

Abstract

We have extended Cosmos++, a multi-dimensional unstructured adaptive mesh code for solving the covariant Newtonian and general relativistic radiation magnetohydrodynamic (MHD) equations, to accommodate both discrete finite volume and arbitrarily high order finite element structures. The new finite element implementation, called CosmosDG, is based on a discontinuous Galerkin (DG) formulation, using both entropy-based artificial viscosity and slope limiting procedures for regularization of shocks. High order multi-stage forward Euler and strong stability preserving Runge-Kutta time integration options complement high order spatial discretization. We have also added flexibility in the code infrastructure allowing for both adaptive mesh and adaptive basis order refinement to be performed separately or simultaneously in a local (cell-by-cell) manner. We discuss in this report the DG formulation and present tests demonstrating the robustness, accuracy, and convergence of our numerical methods applied to special and general relativistic MHD, though we note an equivalent capability currently also exists in CosmosDG for Newtonian systems.

21 pages, 14 figures, submitted to The Astrophysical Journal

References in corpus (17)

• Numerical simulations of super-critical black hole accretion flows in general relativity
• Primitive Variable Solvers for Conservative General Relativistic Magnetohydrodynamics
• Numerical 3+1 general relativistic magnetohydrodynamics: a local characteristic approach
• A Numerical Method for General Relativistic Magnetohydrodynamics
• Cosmos++: Relativistic Magnetohydrodynamics on Unstructured Grids with Local Adaptive Refinement
• A five-wave HLL Riemann solver for relativistic MHD
• SpECTRE: A Task-based Discontinuous Galerkin Code for Relativistic Astrophysics
• A Second Order Godunov Method for Multidimensional Relativistic Magnetohydrodynamics
• Solving the relativistic magnetohydrodynamics equations with ADER discontinuous Galerkin methods, a posteriori subcell limiting and adaptive mesh refinement
• Discontinuous Galerkin methods for general-relativistic hydrodynamics: formulation and application to spherically symmetric spacetimes
• Numerical Simulations of Optically Thick Accretion onto a Black Hole - I. Spherical Case
• Numerical Simulations of Optically Thick Accretion onto a Black Hole - II. Rotating Flow
• Astrophysical hydrodynamics with a high-order discontinuous Galerkin scheme and adaptive mesh refinement
• Non-Oscillatory Central Difference and Artificial Viscosity Schemes for Relativistic Hydrodynamics
• A discontinuous Galerkin method for solving the fluid and MHD equations in astrophysical simulations
• Multi-Dimensional Astrophysical Structural and Dynamical Analysis I. Development of a Nonlinear Finite Element Approach
• COSMOS: A Radiation-Chemo-Hydrodynamics Code for Astrophysical Problems

Cited by in corpus (14)

• First Sagittarius A* Event Horizon Telescope Results. V. Testing Astrophysical Models of the Galactic Center Black Hole
• The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project
• Constrained transport and adaptive mesh refinement in the Black Hole Accretion Code
• High-order Magnetohydrodynamics for Astrophysics with an Adaptive Mesh Refinement Discontinuous Galerkin Scheme
• Relativistic Tidal Disruption and Nuclear Ignition of White Dwarf Stars by Intermediate Mass Black Holes
• Multi-Frequency General Relativistic Radiation-Hydrodynamics with ${\mathbf{M}}_1$ Closure
• General-relativistic neutron star evolutions with the discontinuous Galerkin method
• Non-linear evolutions of magnetised thick discs around black holes: dependence on the initial data
• General Relativistic Implicit Monte Carlo Radiation-Hydrodynamics
• Nuclear Ignition of White Dwarf Stars by Relativistic Encounters with Rotating Intermediate Mass Black Holes
• Spectral Difference method with a posteriori limiting: Application to the Euler equations in one and two space dimensions
• Divergence-Free Magnetohydrodynamics on Conformally Moving, Adaptive Meshes Using a Vector Potential Method
• Future Simulations of Tidal Disruption Events
• Resolution Study of Thermonuclear Initiation in White Dwarf Tidal Disruption Events