Extending the lifetime of 3D black hole computations with a new hyperbolic system of evolution equations

We review recent efforts to re-formulate the Einstein equations for fully relativistic numerical simulations. The so-called numerical relativity (computational simulations in general relativity) is a promising research field matching with ongoing astrophysical observations such as gravitational wave astronomy. Many trials for longterm stable and accurate simulations of binary compact objects have revealed that mathematically equivalent sets of evolution equations show differe...

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 differencing. In this paper, we explore the use of a pseudospectral collocation (PSC) method for the evolution of a spherically symmetric black hole spacetime in one dimension using a hyperbolic formulation of Einstein's equations. We demonstrate that...

We discuss a successful three-dimensional cartesian implementation of the Bona-Mass\'o hyperbolic formulation of the 3+1 Einstein evolution equations in numerical relativity. The numerical code, which we call ``Cactus,'' provides a general framework for 3D numerical relativity, and can include various formulations of the evolution equations, initial data sets, and analysis modules. We show important code tests, including dynamically sliced flat space, wave spacetimes, and bla...

Einstein's equations are known to lead to the formation of black holes and spacetime singularities. This appears to be a manifestation of the mathematical phenomenon of finite-time blowup: a formation of singularities from regular initial data. We present a simple hyperbolic system of two semi-linear equations inspired by the Einstein equations. We explore a class of solutions to this system which are analogous to static black-hole models. These solutions exhibit a black-hole...

The first attempts at solving a binary black hole spacetime date back to the 1960s, with the pioneering works of Hahn and Lindquist. In spite of all the computational advances and enormous efforts by several groups, the first stable, long-term evolution of the orbit and merger of two black holes was only accomplished over 40 years later, in 2005. Since then, the field of Numerical Relativity has matured, and been extensively used to explore and uncover a plethora of physical ...

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 present 3D simulations of a scalar field propagating in a Schwarzschild black hole background.

A new constraint suppressing formulation of the Einstein evolution equations is presented, generalizing the five-parameter first-order system due to Kidder, Scheel and Teukolsky (KST). The auxiliary fields, introduced to make the KST system first-order, are given modified evolution equations designed to drive constraint violations toward zero. The algebraic structure of the new system is investigated, showing that the modifications preserve the hyperbolicity of the fundamenta...

Numerical codes based on a direct implementation of the standard ADM formulation of Einstein's equations have generally failed to provide long-term stable and convergent evolutions of black hole spacetimes when excision is used to remove the singularities. We show that, for the case of a single black hole in spherical symmetry, it is possible to circumvent these problems by adding to the evolution equations terms involving the constraints, thus adjusting the standard ADM syst...

Numerical relativity is finally coming of age with the development of massively parallel computers. 3D problems, which were completely intractable several years ago due to limited computer power, can now be performed with grid sizes of about $200^3$. We report on several new codes developed for solving the full 3D Einstein equations, and present results of using them to evolve black holes and gravitational waves.

We present new many-parameter families of strongly and symmetric hyperbolic formulations of Einstein's equations that include quite general algebraic and live gauge conditions for the lapse. The first system that we present has 30 variables and incorporates an algebraic relationship between the lapse and the determinant of the three metric that generalizes the densitized lapse prescription. The second system has 34 variables and uses a family of live gauges that generalizes t...