Formulation of discontinuous Galerkin methods for relativistic astrophysics

Abstract

The DG algorithm is a powerful method for solving pdes, especially for evolution equations in conservation form. Since the algorithm involves integration over volume elements, it is not immediately obvious that it will generalize easily to arbitrary time-dependent curved spacetimes. We show how to formulate the algorithm in such spacetimes for applications in relativistic astrophysics. We also show how to formulate the algorithm for equations in non-conservative form, such as Einstein's field equations themselves. We find two computationally distinct formulations in both cases, one of which has seldom been used before for flat space in curvilinear coordinates but which may be more efficient. We also give a new derivation of the ALE algorithm (Arbitrary Lagrangian-Eulerian) using 4-vector methods that is much simpler than the usual derivation and explains why the method preserves the conservation form of the equations. The various formulations are explored with some simple numerical experiments that also explore the effect of the metric identities on the results.

References in corpus (7)

• Solving Einstein's Equations With Dual Coordinate Frames
• 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
• Astrophysical hydrodynamics with a high-order discontinuous Galerkin scheme and adaptive mesh refinement
• Finite Element, Discontinuous Galerkin, and Finite Difference Evolution Schemes in Spacetime
• Discontinuous Galerkin method for the spherically reduced BSSN system with second-order operators
• Bound-Preserving Discontinuous Galerkin Methods for Conservative Phase Space Advection in Curvilinear Coordinates

Cited by in corpus (24)

• SpECTRE: A Task-based Discontinuous Galerkin Code for Relativistic Astrophysics
• Massively Parallel Simulations of Binary Black Hole Intermediate-Mass-Ratio Inspirals
• Solving 3D relativistic hydrodynamical problems with WENO discontinuous Galerkin methods
• General-relativistic neutron star evolutions with the discontinuous Galerkin method
• The new discontinuous Galerkin methods based numerical relativity program Nmesh
• Realizability-Preserving DG-IMEX Method for the Two-Moment Model of Fermion Transport
• An Operator-Based Local Discontinuous Galerkin Method Compatible With the BSSN Formulation of the Einstein Equations
• Provably Stable Flux Reconstruction High-Order Methods on Curvilinear Elements
• A Physical-Constraint-Preserving Finite Volume WENO Method for Special Relativistic Hydrodynamics on Unstructured Meshes
• A DG-IMEX method for two-moment neutrino transport: Nonlinear solvers for neutrino-matter coupling
• Simulating magnetized neutron stars with discontinuous Galerkin methods
• A high-order shock capturing discontinuous Galerkin-finite-difference hybrid method for GRMHD
• thornado-hydro: a discontinuous Galerkin method for supernova hydrodynamics with nuclear equations of state
• A hp-adaptive discontinuous Galerkin solver for elliptic equations in numerical relativity
• Unified discontinuous Galerkin scheme for a large class of elliptic equations
• Simulating binary black hole mergers using discontinuous Galerkin methods
• A new first-order formulation of the Einstein equations exploiting analogies with electrodynamics
• General relativistic force-free electrodynamics with a discontinuous Galerkin-finite difference hybrid method
• High-order discontinuous Galerkin schemes with subcell finite volume limiter and adaptive mesh refinement for a monolithic first-order BSSNOK formulation of the Einstein-Euler equations
• Discontinuous Galerkin scheme for elliptic equations on extremely stretched grids
• Neutron star evolution by combining discontinuous Galerkin and finite volume methods
• A polytopal discrete de Rham scheme for the exterior calculus Einstein's equations
• Parameter control for binary black hole initial data
• A Discontinuous Galerkin Method for General Relativistic Hydrodynamics in thornado