Related papers: Solving 3D relativistic hydrodynamical problems wi…
Discontinuous Galerkin (DG) methods provide a means to obtain high-order accurate solutions in regions of smooth fluid flow while, with the aid of limiters, still resolving strong shocks. These and other properties make DG methods…
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…
We have developed the formalism necessary to employ the discontinuous-Galerkin approach in general-relativistic hydrodynamics. The formalism is firstly presented in a general 4-dimensional setting and then specialized to the case of…
This paper develops $P^K$-based non-central and central Runge-Kutta discontinuous Galerkin (DG) methods with WENO limiter for the one- and two-dimensional special relativistic magnetohydrodynamical (RMHD) equations, $K=1,2,3$. The…
We present a new approach to stabilizing high-order Runge-Kutta discontinuous Galerkin (RKDG) schemes using weighted essentially non-oscillatory (WENO) reconstructions in the context of hyperbolic conservation laws. In contrast to RKDG…
In astrophysics, the two main methods traditionally in use for solving the Euler equations of ideal fluid dynamics are smoothed particle hydrodynamics and finite volume discretization on a stationary mesh. However, the goal to efficiently…
Simulations of relativistic hydrodynamics often need both high accuracy and robust shock-handling properties. The discontinuous Galerkin method combines these features --- a high order of convergence in regions where the solution is smooth…
We present a new class of high-order accurate numerical algorithms for solving the equations of general-relativistic ideal magnetohydrodynamics in curved spacetimes. In this paper we assume the background spacetime to be given and static,…
A discontinuous Galerkin (DG) method suitable for large-scale astrophysical simulations on Cartesian meshes as well as arbitrary static and moving Voronoi meshes is presented. Most major astrophysical fluid dynamics codes use a finite…
We propose a way to maintain strong consistency and facilitate error analysis in the context of dissipation-based WENO stabilization for continuous and discontinuous Galerkin discretizations of conservation laws. Following Kuzmin and Vedral…
Weighted essentially non-oscillatory (WENO) reconstruction schemes are presented that preserve cylindrical symmetry for radial flows on an equal-angle polar mesh. These new WENO schemes are used with a Lagrangian discontinuous Galerkin (DG)…
A new scheme for communication between overset grids using subcells and Weighted Essentially Non Oscillatory (WENO) reconstruction for two-dimensional problems has been proposed. The effectiveness of this procedure is demonstrated using the…
We present a new perspective on the use of weighted essentially nonoscillatory (WENO) reconstructions in high-order methods for scalar hyperbolic conservation laws. The main focus of this work is on nonlinear stabilization of continuous…
This is the first in a series of papers on the construction and validation of a three-dimensional code for general relativistic hydrodynamics, and its application to general relativistic astrophysics. This paper studies the consistency and…
Hydrodynamical numerical methods that converge with high-order hold particular promise for astrophysical studies, as they can in principle reach prescribed accuracy goals with higher computational efficiency than standard second- or…
Solving the Euler equations of ideal hydrodynamics as accurately and efficiently as possible is a key requirement in many astrophysical simulations. It is therefore important to continuously advance the numerical methods implemented in…
Modern astrophysical simulations aim to accurately model an ever-growing array of physical processes, including the interaction of fluids with magnetic fields, under increasingly stringent performance and scalability requirements driven by…
Discontinuous Galerkin methods are popular because they can achieve high order where the solution is smooth, because they can capture shocks while needing only nearest-neighbor communication, and because they are relatively easy to…
This paper is concerned with developing accurate and efficient numerical methods for fully nonlinear second order elliptic and parabolic partial differential equations (PDEs) in multiple spatial dimensions. It presents a general framework…
This paper develops three high-order accurate discontinuous Galerkin (DG) methods for the one-dimensional (1D) and two-dimensional (2D) nonlinear Dirac (NLD) equations with a general scalar self-interaction. They are the Runge-Kutta DG…