Related papers: A 3+1 Computational Scheme for Dynamic Spherically…
We report on a new 3D numerical code designed to solve the Einstein equations for general vacuum spacetimes. This code is based on the standard 3+1 approach using cartesian coordinates. We discuss the numerical techniques used in developing…
When using the black hole exclusion (horizon boundary condition) technique, $K$ is usually nonzero and spatially variable, so none of the special cases of York's conformal-decomposition algorithm apply, and the full 4-vector nonlinear York…
We present a new technique for time-domain numerical evolution of the scalar field generated by a pointlike scalar charge orbiting a black hole. Time-domain evolution offers an efficient way for calculating black hole perturbations,…
We consider the numerical evolution of dynamic black hole initial data sets with a full 3D, nonlinear evolution code. These data sets consist of single black holes distorted by strong gravitational waves, and mimic the late stages of…
Spacetime singularities in numerical relativity can be avoided by excising a region of the computational domain from inside the apparent horizon. We report on results of such a scheme that is based on using ({\it i}) a horizon locking…
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 evolve a scalar field in a fixed Kerr-Schild background geometry to test simple $(3+1)$-dimensional algorithms for singularity excision. We compare both centered and upwind schemes for handling the shift (advection) terms, as well as…
It was recently shown that spacetime singularities in numerical relativity could be avoided by excising a region inside the apparent horizon in numerical evolutions. In this paper we report on the details of the implementation of this…
We address the hyperboloidal initial value problem in the context of Numerical Relativity, motivated by its evolution on hyperboloidal slices: smooth spacelike slices that reach future null infinity, the "location" in spacetime where…
We report new results which establish that the accurate 3-dimensional numerical simulation of generic single-black-hole spacetimes has been achieved by characteristic evolution with unlimited long term stability. Our results cover a…
We consider the numerical evolution of black hole initial data sets, consisting of single black holes distorted by strong gravitational waves, with a full 3D, nonlinear evolution code. These data sets mimic the late stages of coalescing…
The scheduled launch of the LISA Mission in the next decade has called attention to the gravitational self-force problem. Despite an extensive body of theoretical work, long-time numerical computations of gravitational waves from…
Motivated by the goal for high accuracy modeling of gravitational radiation emitted by isolated systems, recently, there has been renewed interest in the numerical solution of the hyperboloidal initial value problem for Einstein's field…
Numerical relativity has faced the problem that standard 3+1 simulations of black hole spacetimes without singularity excision and with singularity avoiding lapse and vanishing shift fail after an evolution time of around 30-40M due to the…
Current methods of evolving a spacetime containing one or more black holes are plagued by instabilities that prohibit long-term evolution. Some of these instabilities may be due to the numerical method used, traditionally finite…
When using black hole excision to numerically evolve a black hole spacetime with no continuous symmetries, most 3+1 finite differencing codes use a Cartesian grid. It's difficult to do excision on such a grid, because the natural $r =…
Solving the 4-d Einstein equations as evolution in time requires solving equations of two types: the four elliptic initial data (constraint) equations, followed by the six second order evolution equations. Analytically the constraint…
We present a new pseudo-spectral code for the simulation of evolution systems that are second order in space. We test this code by evolving a non-linear scalar wave equation. These non-linear waves can be stably evolved using very simple…
We present long term evolutions of a single black hole of mass $M$ with the BSSN system using pseudospectral methods. For our simulations we use the SGRID code where the BSSN system is implemented in its standard second order in space form.…
Spherically symmetric (1D) black-hole spacetimes are considered as a test for numerical relativity. A finite difference code, based in the hyperbolic structure of Einstein's equations with the harmonic slicing condition is presented.…