Related papers: Excision methods for high resolution shock capturi…
Excision techniques are used in order to deal with black holes in numerical simulations of Einstein equations and consist in removing a topological sphere containing the physical singularity from the numerical domain, applying instead…
A new algorithm for solving the general relativistic MHD equations is described in this paper. We design our scheme to incorporate black hole excision with smooth boundaries, and to simplify solving the combined Einstein and MHD equations…
Multidimensional shock-capturing numerical schemes for special relativistic hydrodynamics (RHD) are computationally more expensive than their correspondent Euler versions, due to the nonlinear relations between conservative and primitive…
We present numerical hydrodynamical evolutions of rapidly rotating relativistic stars, using an axisymmetric, nonlinear relativistic hydrodynamics code. We use four different high-resolution shock-capturing (HRSC) finite-difference schemes…
Fallback supernovae and the collapsar scenario for long-gamma ray burst and hypernovae have received considerable interest as pathways to black-hole formation and extreme transient events. Consistent simulations of these scenarios require a…
We analyze the excision strategy for simulating black holes. The problem is modeled by the propagation of quasi-linear waves in a 1-dimensional spatial region with timelike outer boundary, spacelike inner boundary and a horizon in between.…
We present a numerical method to solve the equations of general relativistic hydrodynamics in a given external gravitational field. The method is based on a generalization of Roe's approximate Riemann solver for the non relativistic Euler…
We present a new numerical code, ECHO, based on an Eulerian Conservative High Order scheme for time dependent three-dimensional general relativistic magnetohydrodynamics (GRMHD) and magnetodynamics (GRMD). ECHO is aimed at providing a…
We generalize the closest point method (CPM) to solve surface partial differential equations with general boundary conditions. The proposed extrapolation method provides a unified framework for treating a broad class of inhomogeneous…
We present a new three-dimensional general-relativistic hydrodynamic evolution scheme coupled to dynamical spacetime evolutions which is capable of efficiently simulating stellar collapse, isolated neutron stars, black hole formation, and…
We introduce a new scheme adaption strategy for one- and two-dimensional hyperbolic systems of conservation laws. The proposed approach builds upon the adaptive framework introduced in [S. Chu, A. Kurganov, and I. Menshov, Appl. Numer.…
A recent paper by Lucas-Serrano et al. indicates that a high-resolution central (HRC) scheme is robust enough to yield accurate hydrodynamical simulations of special relativistic flows in the presence of ultrarelativistic speeds and strong…
We report a new implementation for axisymmetric simulation in full general relativity. In this implementation, the Einstein equations are solved using the Nakamura-Shibata formulation with the so-called cartoon method to impose an…
We present some recent developments on shock capturing methods for nonlinear hyperbolic systems of balance laws, whose prototype is the Euler system of compressible fluid flows, and especially discuss {structure-preserving} techniques. The…
The use of limiting methods for high-order numerical approximations of hyperbolic conservation laws generally requires defining an admissible region/bounds for the solution. In this work, we present a novel approach for computing solution…
We use rigorous techniques from numerical analysis of hyperbolic equations in bounded domains to construct stable finite-difference schemes for Numerical Relativity, in particular for their use in black hole excision. As an application, we…
In this chapter, we aim at presenting the basic techniques necessary to go beyond the widely accepted paradigm of second-order numerics. We specifically focus on finite-volume schemes for hyperbolic conservation laws occuring in fluid…
The paper develops high-order accurate physical-constraints-preserving finite difference WENO schemes for special relativistic hydrodynamical (RHD) equations, built on the local Lax-Friedrich splitting, the WENO reconstruction, the…
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…
This paper studies the physical-constraints-preserving (PCP) Lagrangian finite volume schemes for one- and two-dimensional special relativistic hydrodynamic (RHD) equations. First, the PCP property (i.e. preserving the positivity of the…