14 research outputs found
Characteristic extraction tool for gravitational waveforms
We develop and calibrate a characteristic waveform extraction tool whose major improvements and corrections of prior versions allow satisfaction of the accuracy standards required for advanced LIGO data analysis. The extraction tool uses a characteristic evolution code to propagate numerical data on an inner worldtube supplied by a 3+1 Cauchy evolution to obtain the gravitational waveform at null infinity. With the new extraction tool, high accuracy and convergence of the numerical error can be demonstrated for an inspiral and merger of mass M binary black holes even for an extraction worldtube radius as small as R=20M. The tool provides a means for unambiguous comparison between waveforms generated by evolution codes based upon different formulations of the Einstein equations and based upon different numerical approximations
Strategies for the characteristic extraction of gravitational waveforms
We develop, test, and compare new numerical and geometrical methods for improving the accuracy of extracting waveforms using characteristic evolution. The new numerical method involves use of circular boundaries to the stereographic grid patches which cover the spherical cross sections of the outgoing null cones. We show how an angular version of numerical dissipation can be introduced into the characteristic code to damp the high frequency error arising form the irregular way the circular patch boundary cuts through the grid. The new geometric method involves use of the Weyl tensor component Psi4 to extract the waveform as opposed to the original approach via the Bondi news function. We develop the necessary analytic and computational formula to compute the O(1/r) radiative part of Psi4 in terms of a conformally compactified treatment of null infinity. These methods are compared and calibrated in test problems based upon linearized waves
Binary Black Hole Waveform Extraction at Null Infinity
In this work, we present a work in progress towards an efficient and
economical computational module which interfaces between Cauchy and
characteristic evolution codes. Our goal is to provide a standardized waveform
extraction tool for the numerical relativity community which will allow CCE to
be readily applied to a generic Cauchy code. The tool provides a means of
unambiguous comparison between the waveforms generated by evolution codes based
upon different formulations of the Einstein equations and different numerical
approximation.Comment: 11 pages, 7 figure
Finite difference schemes for second order systems describing black holes
In the harmonic description of general relativity, the principle part of
Einstein's equations reduces to 10 curved space wave equations for the
componenets of the space-time metric. We present theorems regarding the
stability of several evolution-boundary algorithms for such equations when
treated in second order differential form. The theorems apply to a model black
hole space-time consisting of a spacelike inner boundary excising the
singularity, a timelike outer boundary and a horizon in between. These
algorithms are implemented as stable, convergent numerical codes and their
performance is compared in a 2-dimensional excision problem.Comment: 19 pages, 9 figure
Testing numerical relativity with the shifted gauge wave
Computational methods are essential to provide waveforms from coalescing
black holes, which are expected to produce strong signals for the gravitational
wave observatories being developed. Although partial simulations of the
coalescence have been reported, scientifically useful waveforms have so far not
been delivered. The goal of the AppleswithApples (AwA) Alliance is to design,
coordinate and document standardized code tests for comparing numerical
relativity codes. The first round of AwA tests have now being completed and the
results are being analyzed. These initial tests are based upon periodic
boundary conditions designed to isolate performance of the main evolution code.
Here we describe and carry out an additional test with periodic boundary
conditions which deals with an essential feature of the black hole excision
problem, namely a non-vanishing shift. The test is a shifted version of the
existing AwA gauge wave test. We show how a shift introduces an exponentially
growing instability which violates the constraints of a standard harmonic
formulation of Einstein's equations. We analyze the Cauchy problem in a
harmonic gauge and discuss particular options for suppressing instabilities in
the gauge wave tests. We implement these techniques in a finite difference
evolution algorithm and present test results. Although our application here is
limited to a model problem, the techniques should benefit the simulation of
black holes using harmonic evolution codes.Comment: Submitted to special numerical relativity issue of Classical and
Quantum Gravit
Implementation of standard testbeds for numerical relativity
We discuss results that have been obtained from the implementation of the
initial round of testbeds for numerical relativity which was proposed in the
first paper of the Apples with Apples Alliance. We present benchmark results
for various codes which provide templates for analyzing the testbeds and to
draw conclusions about various features of the codes. This allows us to sharpen
the initial test specifications, design a new test and add theoretical insight.Comment: Corrected versio
Characteristic Evolution and Matching
I review the development of numerical evolution codes for general relativity
based upon the characteristic initial value problem. Progress in characteristic
evolution is traced from the early stage of 1D feasibility studies to 2D
axisymmetric codes that accurately simulate the oscillations and gravitational
collapse of relativistic stars and to current 3D codes that provide pieces of a
binary black hole spacetime. Cauchy codes have now been successful at
simulating all aspects of the binary black hole problem inside an artificially
constructed outer boundary. A prime application of characteristic evolution is
to extend such simulations to null infinity where the waveform from the binary
inspiral and merger can be unambiguously computed. This has now been
accomplished by Cauchy-characteristic extraction, where data for the
characteristic evolution is supplied by Cauchy data on an extraction worldtube
inside the artificial outer boundary. The ultimate application of
characteristic evolution is to eliminate the role of this outer boundary by
constructing a global solution via Cauchy-characteristic matching. Progress in
this direction is discussed.Comment: New version to appear in Living Reviews 2012. arXiv admin note:
updated version of arXiv:gr-qc/050809
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