Related papers: Excision methods for high resolution shock capturi…
We solve Einstein's field equations coupled to relativistic hydrodynamics in full 3+1 general relativity to evolve astrophysical systems characterized by strong gravitational fields. We model rotating, collapsing and binary stars by…
Higher order finite volume schemes for magnetohydrodynamics (MHD) and relativistic magnetohydrodynamics (RMHD) are very valuable because they allow us to carry out astrophysical simulations with very high accuracy. However, astrophysical…
We present a finite volume scheme, first on the Burgers equations, then on the Euler equations, based on a conservative reconstruction of shocks inside each cells of the mesh. Its main features are the following points. First, the scheme is…
We present a method for general relativistic smoothed particle hydrodynamics (GRSPH), based on an entropy-conservative form of the general relativistic hydrodynamic equations for a perfect fluid. We aim to replace approximate treatments of…
High-order reconstruction schemes for the solution of hyperbolic conservation laws in orthogonal curvilinear coordinates are revised in the finite volume approach. The formulation employs a piecewise polynomial approximation to the…
We propose novel less diffusive schemes for conservative one- and two-dimensional hyperbolic systems of nonlinear partial differential equations (PDEs). The main challenges in the development of accurate and robust numerical methods for the…
Central-upwind (CU) schemes are Riemann-problem-solver-free finite-volume methods widely applied to a variety of hyperbolic systems of PDEs. Exact solutions of these systems typically satisfy certain bounds, and it is highly desirable or…
This paper develops the genuinely multidimensional HLL Riemann solver and finite volume scheme for the two-dimensional special relativistic hydrodynamic equations on Cartesian meshes and studies its physical-constraint-preserving (PCP)…
We develop a numerical hydrodynamics code using a pseudo-Newtonian formulation that uses the weak field approximation for the geometry, and a generalized source term for the Poisson equation that takes into account relativistic effects. The…
An implicit algorithm for solving the equations of general relativistic hydrodynamics in conservative form in three-dimensional axi-symmetry is presented. This algorithm is a direct extension of the pseudo-Newtonian implicit radiative…
Pressure Poisson equation (PPE) reformulations of the incompressible Navier-Stokes equations (NSE) replace the incompressibility constraint by a Poisson equation for the pressure and a suitable choice of boundary conditions. This yields a…
When using black hole excision to numerically evolve a fully generic black hole spacetime, most 3-D 3+1 codes use an $xyz$-topology (spatial) grid. In such a grid, an $r = \constant$ excision surface must be approximated by an irregular and…
We present ECHO-QGP, a numerical code for $(3+1)$-dimensional relativistic viscous hydrodynamics designed for the modeling of the space-time evolution of the matter created in high energy nuclear collisions. The code has been built on top…
In this paper we propose the first better than second order accurate method in space and time for the numerical solution of the resistive relativistic magnetohydrodynamics (RRMHD) equations on unstructured meshes in multiple space…
We assess the suitability of a recent high-resolution central scheme developed by Kurganov & Tadmor (2000) for the solution of the relativistic hydrodynamics equations. The novelty of this approach relies on the absence of Riemann solvers…
The second paper of this series presents two robust entropy stable shock-capturing methods for discontinuous Galerkin spectral element (DGSEM) discretizations of the compressible magneto-hydrodynamics (MHD) equations. Specifically, we use…
This paper develops the genuinely multidimensional HLL Riemann solver for the two-dimensional special relativistic hydrodynamic equations on Cartesian meshes and studies its physical-constraint preserving (PCP) property. Based on the…
In the present study, a characteristic-based boundary condition scheme is developed for the compressible magnetohydrodynamic (MHD) equations in the general curvilinear coordinate system, which is an extension of the characteristic boundary…
We investigate how the accuracy and stability of numerical relativity simulations of 1D colliding plane waves depends on choices of equation formulations, gauge conditions, boundary conditions, and numerical methods, all in the context of a…
We compare high-order methods including spectral difference (SD), flux reconstruction (FR), the entropy-stable discontinuous Galerkin spectral element method (ES-DGSEM), modal discontinuous Galerkin methods, and WENO to select the best…