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

I review recent developments in numerical relativity, focussing on progress made in 3D black hole evolution. Progress in development of black hole initial data, apparent horizon boundary conditions, adaptive mesh refinement, and characteristic evolution is highlighted, as well as full 3D simulations of colliding and distorted black holes. For true 3D distorted holes, with Cauchy evolution techniques, it is now possible to extract highly accurate, nonaxisymmetric waveforms fro...

Techniques are developed for projecting the solutions of symmetric hyperbolic evolution systems onto the constraint submanifold (the constraint-satisfying subset of the dynamical field space). These optimal projections map a field configuration to the ``nearest'' configuration in the constraint submanifold, where distances between configurations are measured with the natural metric on the space of dynamical fields. The construction and use of these projections is illustrated ...

We study the stability of three-dimensional numerical evolutions of the Einstein equations, comparing the standard ADM formulation to variations on a family of formulations that separate out the conformal and traceless parts of the system. We develop an implementation of the conformal-traceless (CT) approach that has improved stability properties in evolving weak and strong gravitational fields, and for both vacuum and spacetimes with active coupling to matter sources. Cases ...

The demand to obtain answers to highly complex problems within strong-field gravity has been met with significant progress in the numerical solution of Einstein's equations - along with some spectacular results - in various setups. We review techniques for solving Einstein's equations in generic spacetimes, focusing on fully nonlinear evolutions but also on how to benchmark those results with perturbative approaches. The results address problems in high-energy physics, hologr...

We are entering an era where the numerical construction of generic spacetimes is becoming a reality. The use of computer simulations, in principle, allows us to solve Einstein equations in their full generality and unravel important messages so far hidden in the theory. Despite the problems still found in present applications, significant progress is being achieved in a number of areas where simulations are on the verge of providing useful physical information. This article r...

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 gauge conditions we use is able to drive highly distorted, rotating black holes to an almost static state at late times, with well behaved metric functions, without the need for any special initial conditions or analytically prescribed gauge ...

We propose a re-formulation of the Einstein evolution equations that cleanly separates the conformal degrees of freedom and the non-conformal degrees of freedom with the latter satisfying a first order strongly hyperbolic system. The conformal degrees of freedom are taken to be determined by the choice of slicing and the initial data, and are regarded as given functions (along with the lapse and the shift) in the hyperbolic part of the evolution. We find that there is a two...

The Good-Bad-Ugly-F model is a system of semi-linear wave equations that mimics the asymptotic form of the Einstein field equations in generalized harmonic gauge with specific constraint damping and suitable gauge source functions. These constraint additions and gauge source functions eliminate logarithmic divergences appearing at the leading order in the asymptotic expansion of the metric components. In this work, as a step towards using compactified hyperboloidal slices in ...

We introduce a class of nonlinear hyperbolic conservation laws on a Schwarzschild black hole background and derive several properties satisfied by (possibly weak) solutions. Next, we formulate a numerical approximation scheme which is based on the finite volume methodology and takes the curved geometry into account. An interesting feature of our model is that no boundary conditions is required at the black hole horizon boundary. We establish that this scheme converges to an e...

Motivated by the need to control the exponential growth of constraint violations in numerical solutions of the Einstein evolution equations, two methods are studied here for controlling this growth in general hyperbolic evolution systems. The first method adjusts the evolution equations dynamically, by adding multiples of the constraints, in a way designed to minimize this growth. The second method imposes special constraint preserving boundary conditions on the incoming comp...