Related papers: Black Hole Excision with Multiple Grid Patches
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…
It is expected that the realization of a convergent and long-term stable numerical code for the simulation of a black hole inspiral collision will depend greatly upon the construction of stable algorithms capable of handling smooth and,…
This is the second in a series of papers describing a 3+1 computational scheme for the numerical simulation of dynamic black hole spacetimes. We discuss the numerical time-evolution of a given black-hole-containing initial data slice in…
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…
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.…
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 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…
We describe an explicit in time, finite-difference code designed to simulate black holes by using the excision method. The code is based upon the harmonic formulation of the Einstein equations and incorporates several features regarding the…
We perform both distorted black hole evolutions and binary black hole head on collisions and compare the results of using a full grid to results obtained by excising the black hole interiors. In both cases the evolutions are found to run…
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…
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 describe a generic infrastructure for time evolution simulations in numerical relativity using multiple grid patches. After a motivation of this approach, we discuss the relative advantages of global and patch-local tensor bases. We…
We present a singularity excision algorithm appropriate for numerical simulations of black holes moving throughout the computational domain. The method is an extension of the excision procedure previously used to obtain stable simulations…
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…
Simulations of binary black hole systems using the Spectral Einstein Code (SpEC) are done on a computational domain that excises the regions inside the black holes. It is imperative that the excision boundaries are outflow boundaries with…
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…
We extend previous work on 3D black hole excision to the case of distorted black holes, with a variety of dynamic gauge conditions that are able to respond naturally to the spacetime dynamics. We show that the combination of excision and…
Recent demonstrations of unexcised black holes traversing across computational grids represent a significant advance in numerical relativity. Stable and accurate simulations of multiple orbits, and their radiated waves, result. This…
We present numerical simulations of orbiting black holes for around twelve cycles, using a high-order multipatch approach. Unlike some other approaches, the computational speed scales almost perfectly for thousands of processors. Multipatch…
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…