Related papers: Black hole evolution by spectral methods
We have recently constructed a numerical code that evolves a spherically symmetric spacetime using a hyperbolic formulation of Einstein's equations. For the case of a Schwarzschild black hole, this code works well at early times, but…
We describe a numerical code that solves Einstein's equations for a Schwarzschild black hole in spherical symmetry, using a hyperbolic formulation introduced by Choquet-Bruhat and York. This is the first time this formulation has been used…
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 consider unconstrained evolution schemes for the hyperboloidal initial value problem in numerical relativity as a promising candidate for the optimally efficient numerical treatment of radiating compact objects. Here, spherical symmetry…
We present a new many-parameter family of hyperbolic representations of Einstein's equations, which we obtain by a straightforward generalization of previously known systems. We solve the resulting evolution equations numerically for a…
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.…
Numerical relativity has traditionally been pursued via finite differencing. Here we explore pseudospectral collocation (PSC) as an alternative to finite differencing, focusing particularly on the solution of the Hamiltonian constraint (an…
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…
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…
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 present a code for solving the coupled Einstein-hydrodynamics equations to evolve relativistic, self-gravitating fluids. The Einstein field equations are solved in generalized harmonic coordinates on one grid using pseudospectral…
A numerical solution scheme for the Einstein field equations based on generalized harmonic coordinates is described, focusing on details not provided before in the literature and that are of particular relevance to the binary black hole…
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…
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…
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 report on the successful numerical evolution of the compactified hyperboloidal initial value problem in general relativity using generalized harmonic gauge. We work in spherical symmetry, using a massless scalar field to drive dynamics.…
Partial differential equations exhibiting an anisotropic scaling between space and time -- such as those of Horava-Lifshitz gravity -- have a dispersive nature. They contain higher-order spatial derivatives, but remain second order in time.…
In this paper we consider the single patch pseudo-spectral scheme for tensorial and spinorial evolution problems on the 2-sphere presented in [3,4] which is based on the spin-weighted spherical harmonics transform. We apply and extend this…
This is the first of a series of papers describing a numerical implementation of the conformally rescaled Einstein equation, an implementation designed to calculate asymptotically flat spacetimes, especially spacetimes containing black…
Three dimensional (3D) numerical evolutions of static black holes with excision are presented. These evolutions extend to about 8000M, where M is the mass of the black hole. This degree of stability is achieved by using growth-rate…