18 research outputs found
Stable radiation-controlling boundary conditions for the generalized harmonic Einstein equations
This paper is concerned with the initial-boundary value problem for the
Einstein equations in a first-order generalized harmonic formulation. We impose
boundary conditions that preserve the constraints and control the incoming
gravitational radiation by prescribing data for the incoming fields of the Weyl
tensor. High-frequency perturbations about any given spacetime (including a
shift vector with subluminal normal component) are analyzed using the
Fourier-Laplace technique. We show that the system is boundary-stable. In
addition, we develop a criterion that can be used to detect weak instabilities
with polynomial time dependence, and we show that our system does not suffer
from such instabilities. A numerical robust stability test supports our claim
that the initial-boundary value problem is most likely to be well-posed even if
nonzero initial and source data are included.Comment: 27 pages, 4 figures; more numerical results and references added,
several minor amendments; version accepted for publication in Class. Quantum
Gra
Implementation of higher-order absorbing boundary conditions for the Einstein equations
We present an implementation of absorbing boundary conditions for the
Einstein equations based on the recent work of Buchman and Sarbach. In this
paper, we assume that spacetime may be linearized about Minkowski space close
to the outer boundary, which is taken to be a coordinate sphere. We reformulate
the boundary conditions as conditions on the gauge-invariant
Regge-Wheeler-Zerilli scalars. Higher-order radial derivatives are eliminated
by rewriting the boundary conditions as a system of ODEs for a set of auxiliary
variables intrinsic to the boundary. From these we construct boundary data for
a set of well-posed constraint-preserving boundary conditions for the Einstein
equations in a first-order generalized harmonic formulation. This construction
has direct applications to outer boundary conditions in simulations of isolated
systems (e.g., binary black holes) as well as to the problem of
Cauchy-perturbative matching. As a test problem for our numerical
implementation, we consider linearized multipolar gravitational waves in TT
gauge, with angular momentum numbers l=2 (Teukolsky waves), 3 and 4. We
demonstrate that the perfectly absorbing boundary condition B_L of order L=l
yields no spurious reflections to linear order in perturbation theory. This is
in contrast to the lower-order absorbing boundary conditions B_L with L<l,
which include the widely used freezing-Psi_0 boundary condition that imposes
the vanishing of the Newman-Penrose scalar Psi_0.Comment: 25 pages, 9 figures. Minor clarifications. Final version to appear in
Class. Quantum Grav
Testing outer boundary treatments for the Einstein equations
Various methods of treating outer boundaries in numerical relativity are
compared using a simple test problem: a Schwarzschild black hole with an
outgoing gravitational wave perturbation. Numerical solutions computed using
different boundary treatments are compared to a `reference' numerical solution
obtained by placing the outer boundary at a very large radius. For each
boundary treatment, the full solutions including constraint violations and
extracted gravitational waves are compared to those of the reference solution,
thereby assessing the reflections caused by the artificial boundary. These
tests use a first-order generalized harmonic formulation of the Einstein
equations. Constraint-preserving boundary conditions for this system are
reviewed, and an improved boundary condition on the gauge degrees of freedom is
presented. Alternate boundary conditions evaluated here include freezing the
incoming characteristic fields, Sommerfeld boundary conditions, and the
constraint-preserving boundary conditions of Kreiss and Winicour. Rather
different approaches to boundary treatments, such as sponge layers and spatial
compactification, are also tested. Overall the best treatment found here
combines boundary conditions that preserve the constraints, freeze the
Newman-Penrose scalar Psi_0, and control gauge reflections.Comment: Modified to agree with version accepted for publication in Class.
Quantum Gra
Reducing orbital eccentricity in binary black hole simulations
Binary black hole simulations starting from quasi-circular (i.e., zero radial
velocity) initial data have orbits with small but non-zero orbital
eccentricities. In this paper the quasi-equilibrium initial-data method is
extended to allow non-zero radial velocities to be specified in binary black
hole initial data. New low-eccentricity initial data are obtained by adjusting
the orbital frequency and radial velocities to minimize the orbital
eccentricity, and the resulting ( orbit) evolutions are compared with
those of quasi-circular initial data. Evolutions of the quasi-circular data
clearly show eccentric orbits, with eccentricity that decays over time. The
precise decay rate depends on the definition of eccentricity; if defined in
terms of variations in the orbital frequency, the decay rate agrees well with
the prediction of Peters (1964). The gravitational waveforms, which contain
cycles in the dominant l=m=2 mode, are largely unaffected by the
eccentricity of the quasi-circular initial data. The overlap between the
dominant mode in the quasi-circular evolution and the same mode in the
low-eccentricity evolution is about 0.99.Comment: 27 pages, 9 figures; various minor clarifications; accepted to the
"New Frontiers" special issue of CQ
Towards absorbing outer boundaries in General Relativity
We construct exact solutions to the Bianchi equations on a flat spacetime
background. When the constraints are satisfied, these solutions represent in-
and outgoing linearized gravitational radiation. We then consider the Bianchi
equations on a subset of flat spacetime of the form [0,T] x B_R, where B_R is a
ball of radius R, and analyze different kinds of boundary conditions on
\partial B_R. Our main results are: i) We give an explicit analytic example
showing that boundary conditions obtained from freezing the incoming
characteristic fields to their initial values are not compatible with the
constraints. ii) With the help of the exact solutions constructed, we determine
the amount of artificial reflection of gravitational radiation from
constraint-preserving boundary conditions which freeze the Weyl scalar Psi_0 to
its initial value. For monochromatic radiation with wave number k and arbitrary
angular momentum number l >= 2, the amount of reflection decays as 1/(kR)^4 for
large kR. iii) For each L >= 2, we construct new local constraint-preserving
boundary conditions which perfectly absorb linearized radiation with l <= L.
(iv) We generalize our analysis to a weakly curved background of mass M, and
compute first order corrections in M/R to the reflection coefficients for
quadrupolar odd-parity radiation. For our new boundary condition with L=2, the
reflection coefficient is smaller than the one for the freezing Psi_0 boundary
condition by a factor of M/R for kR > 1.04. Implications of these results for
numerical simulations of binary black holes on finite domains are discussed.Comment: minor revisions, 30 pages, 6 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
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
A minimization problem for the lapse and the initial-boundary value problem for Einstein's field equations
Testing gravitational-wave searches with numerical relativity waveforms: Results from the first Numerical INJection Analysis (NINJA) project
The Numerical INJection Analysis (NINJA) project is a collaborative effort
between members of the numerical relativity and gravitational-wave data
analysis communities. The purpose of NINJA is to study the sensitivity of
existing gravitational-wave search algorithms using numerically generated
waveforms and to foster closer collaboration between the numerical relativity
and data analysis communities. We describe the results of the first NINJA
analysis which focused on gravitational waveforms from binary black hole
coalescence. Ten numerical relativity groups contributed numerical data which
were used to generate a set of gravitational-wave signals. These signals were
injected into a simulated data set, designed to mimic the response of the
Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this
data using search and parameter-estimation pipelines. Matched filter
algorithms, un-modelled-burst searches and Bayesian parameter-estimation and
model-selection algorithms were applied to the data. We report the efficiency
of these search methods in detecting the numerical waveforms and measuring
their parameters. We describe preliminary comparisons between the different
search methods and suggest improvements for future NINJA analyses.Comment: 56 pages, 25 figures; various clarifications; accepted to CQ