4 research outputs found
Initial data transients in binary black hole evolutions
We describe a method for initializing characteristic evolutions of the
Einstein equations using a linearized solution corresponding to purely outgoing
radiation. This allows for a more consistent application of the characteristic
(null cone) techniques for invariantly determining the gravitational radiation
content of numerical simulations. In addition, we are able to identify the {\em
ingoing} radiation contained in the characteristic initial data, as well as in
the initial data of the 3+1 simulation. We find that each component leads to a
small but long lasting (several hundred mass scales) transient in the measured
outgoing gravitational waves.Comment: 18 pages, 4 figure
Simulation of Binary Black Hole Spacetimes with a Harmonic Evolution Scheme
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 problem. This includes demonstrations of the effectiveness of
constraint damping, and how the time slicing can be controlled through the use
of a source function evolution equation. In addition, some results from an
ongoing study of binary black hole coalescence, where the black holes are
formed via scalar field collapse, are shown. Scalar fields offer a convenient
route to exploring certain aspects of black hole interactions, and one
interesting, though tentative suggestion from this early study is that behavior
reminiscent of "zoom-whirl" orbits in particle trajectories is also present in
the merger of equal mass, non-spinning binaries, with appropriately fine-tuned
initial conditions.Comment: 16 pages, 14 figures; replaced with published versio
The discrete energy method in numerical relativity: Towards long-term stability
The energy method can be used to identify well-posed initial boundary value
problems for quasi-linear, symmetric hyperbolic partial differential equations
with maximally dissipative boundary conditions. A similar analysis of the
discrete system can be used to construct stable finite difference equations for
these problems at the linear level. In this paper we apply these techniques to
some test problems commonly used in numerical relativity and observe that while
we obtain convergent schemes, fast growing modes, or ``artificial
instabilities,'' contaminate the solution. We find that these growing modes can
partially arise from the lack of a Leibnitz rule for discrete derivatives and
discuss ways to limit this spurious growth.Comment: 18 pages, 22 figure
From Geometry to Numerics: interdisciplinary aspects in mathematical and numerical relativity
This article reviews some aspects in the current relationship between
mathematical and numerical General Relativity. Focus is placed on the
description of isolated systems, with a particular emphasis on recent
developments in the study of black holes. Ideas concerning asymptotic flatness,
the initial value problem, the constraint equations, evolution formalisms,
geometric inequalities and quasi-local black hole horizons are discussed on the
light of the interaction between numerical and mathematical relativists.Comment: Topical review commissioned by Classical and Quantum Gravity.
Discussion inspired by the workshop "From Geometry to Numerics" (Paris, 20-24
November, 2006), part of the "General Relativity Trimester" at the Institut
Henri Poincare (Fall 2006). Comments and references added. Typos corrected.
Submitted to Classical and Quantum Gravit