Central Runge-Kutta discontinuous Galerkin methods for the special relativistic hydrodynamics

TL;DR

This paper introduces central Runge-Kutta discontinuous Galerkin methods with WENO limiter for solving 1D and 2D special relativistic hydrodynamics equations, achieving stability, accuracy, and efficiency.

Contribution

The paper develops a novel central DG method with WENO limiter for relativistic hydrodynamics, avoiding numerical flux calculations and improving computational efficiency.

Findings

Methods are stable and accurate for complex RHD problems.

WENO limiter effectively reduces computational cost by limiting troubled cells.

Numerical tests confirm robustness and high resolution of the proposed methods.

Abstract

This paper developes Runge-Kutta $P^K$-based central discontinuous Galerkin (CDG) methods with WENO limiter to the one- and two-dimensional special relativistic hydrodynamical (RHD) equations, $K=1,2,3$. Different from the non-central DG methods, the \CDG{} have to find two approximate solutions defined on mutually dual meshes. For each mesh, the CDG approximate solutions on its dual mesh are used to calculate the flux values in the cell and on the cell boundary so that the approximate solutions on mutually dual meshes are coupled with each other, and the use of numerical flux may be avoided. The WENO limiter is adaptively implemented via two steps: the "troubled" cells are first identified by using a modified TVB minmod function, and then the WENO technique is used to locally reconstruct new polynomials of degree $(2K+1)$ replacing the CDG solutions inside the "troubled' cells by the cell average values of the CDG solutions in the neighboring cells as well as the original cell averages of the "troubled" cells. Because the WENO limiter is only employed for finite "troubled" cells, the computational cost can be as little as possible. The accuracy of the CDG without the numerical dissipation is analyzed and calculation of the flux integrals over the cells is also addressed. Several test problems in one and two dimensions are solved by using our \CDG{} with WENO limiter. The computations demonstrate that our methods are stable, accurate, and robust in solving complex RHD problems.