Simulating Accretion Disks With Discontinous Galerkin Methods

Discontinuous Galerkin methods have recently seen a lot of interest as a new
method for relativistic astrophysics simulations. Many test problems for the
GRMHD equations have been studied, including 3D TOV stars, blast waves, and
single black hole evolutions. A natural next problem to study is accretion disks
on a stationary background. We present results of accretion disk simulations
using discontiuous Galerkin methods using the SpECTRE
(https://github.com/sxs-collaboration/spectre) code and compare the results to
the Spectral Einstein Code.