Neutron star evolutions with WENO discontinuous Galerkin methods

Discontinuous Galerkin methods combine high accuracy with shock-handling properties. This combination makes them well-suited to solve the PDEs describing relativistic astrophysical scenarios such as neutron star mergers and accretion discs, where fluid shocks and surfaces naturally arise. However, the development of schemes for handling shocks in a discontinuous Galerkin code is ongoing, with recent efforts working toward increasing the robustness and accuracy of these schemes. Using the SpECTRE (https://github.com/sxs-collaboration/spectre) code, we compare different shock-handling schemes: simple slope limiters and more recent, higher-order WENO-based limiters. We evaluate the performance of these different limiters for standard shock problems and neutron star evolutions.