435 research outputs found
Numerical relativity with characteristic evolution, using six angular patches
The characteristic approach to numerical relativity is a useful tool in
evolving gravitational systems. In the past this has been implemented using two
patches of stereographic angular coordinates. In other applications, a
six-patch angular coordinate system has proved effective. Here we investigate
the use of a six-patch system in characteristic numerical relativity, by
comparing an existing two-patch implementation (using second-order finite
differencing throughout) with a new six-patch implementation (using either
second- or fourth-order finite differencing for the angular derivatives). We
compare these different codes by monitoring the Einstein constraint equations,
numerically evaluated independently from the evolution. We find that, compared
to the (second-order) two-patch code at equivalent resolutions, the errors of
the second-order six-patch code are smaller by a factor of about 2, and the
errors of the fourth-order six-patch code are smaller by a factor of nearly 50.Comment: 12 pages, 5 figures, submitted to CQG (special NFNR issue
Study of multi black hole and ring singularity apparent horizons
We study critical black hole separations for the formation of a common
apparent horizon in systems of - black holes in a time symmetric
configuration. We study in detail the aligned equal mass cases for ,
and relate them to the unequal mass binary black hole case. We then study the
apparent horizon of the time symmetric initial geometry of a ring singularity
of different radii. The apparent horizon is used as indicative of the location
of the event horizon in an effort to predict a critical ring radius that would
generate an event horizon of toroidal topology. We found that a good estimate
for this ring critical radius is . We briefly discuss the
connection of this two cases through a discrete black hole 'necklace'
configuration.Comment: 31 pages, 21 figure
Black Hole--Scalar Field Interactions in Spherical Symmetry
We examine the interactions of a black hole with a massless scalar field
using a coordinate system which extends ingoing Eddington-Finkelstein
coordinates to dynamic spherically symmetric-spacetimes. We avoid problems with
the singularity by excising the region of the black hole interior to the
apparent horizon. We use a second-order finite difference scheme to solve the
equations. The resulting program is stable and convergent and will run forever
without problems. We are able to observe quasi-normal ringing and power-law
tails as well an interesting nonlinear feature.Comment: 16 pages, 26 figures, RevTex, to appear in Phys. Rev.
Initial Data and Coordinates for Multiple Black Hole Systems
We present here an alternative approach to data setting for spacetimes with
multiple moving black holes generalizing the Kerr-Schild form for rotating or
non-rotating single black holes to multiple moving holes. Because this scheme
preserves the Kerr-Schild form near the holes, it selects out the behaviour of
null rays near the holes, may simplify horizon tracking, and may prove useful
in computational applications. For computational evolution, a discussion of
coordinates (lapse function and shift vector) is given which preserves some of
the properties of the single-hole Kerr-Schild form
Are moving punctures equivalent to moving black holes?
When simulating the inspiral and coalescence of a binary black-hole system,
special care needs to be taken in handling the singularities. Two main
techniques are used in numerical-relativity simulations: A first and more
traditional one ``excises'' a spatial neighbourhood of the singularity from the
numerical grid on each spacelike hypersurface. A second and more recent one,
instead, begins with a ``puncture'' solution and then evolves the full
3-metric, including the singular point. In the continuum limit, excision is
justified by the light-cone structure of the Einstein equations and, in
practice, can give accurate numerical solutions when suitable discretizations
are used. However, because the field variables are non-differentiable at the
puncture, there is no proof that the moving-punctures technique is correct,
particularly in the discrete case. To investigate this question we use both
techniques to evolve a binary system of equal-mass non-spinning black holes. We
compare the evolution of two curvature 4-scalars with proper time along the
invariantly-defined worldline midway between the two black holes, using
Richardson extrapolation to reduce the influence of finite-difference
truncation errors. We find that the excision and moving-punctures evolutions
produce the same invariants along that worldline, and thus the same spacetimes
throughout that worldline's causal past. This provides convincing evidence that
moving-punctures are indeed equivalent to moving black holes.Comment: 4 pages, 3 eps color figures; v2 = major revisions to introduction &
conclusions based on referee comments, but no change in analysis or result
Revisiting Event Horizon Finders
Event horizons are the defining physical features of black hole spacetimes,
and are of considerable interest in studying black hole dynamics. Here, we
reconsider three techniques to localise event horizons in numerical spacetimes:
integrating geodesics, integrating a surface, and integrating a level-set of
surfaces over a volume. We implement the first two techniques and find that
straightforward integration of geodesics backward in time to be most robust. We
find that the exponential rate of approach of a null surface towards the event
horizon of a spinning black hole equals the surface gravity of the black hole.
In head-on mergers we are able to track quasi-normal ringing of the merged
black hole through seven oscillations, covering a dynamic range of about 10^5.
Both at late times (when the final black hole has settled down) and at early
times (before the merger), the apparent horizon is found to be an excellent
approximation of the event horizon. In the head-on binary black hole merger,
only {\em some} of the future null generators of the horizon are found to start
from past null infinity; the others approach the event horizons of the
individual black holes at times far before merger.Comment: 30 pages, 15 figures, revision
A template bank for gravitational waveforms from coalescing binary black holes: non-spinning binaries
Gravitational waveforms from the inspiral and ring-down stages of the binary
black hole coalescences can be modelled accurately by
approximation/perturbation techniques in general relativity. Recent progress in
numerical relativity has enabled us to model also the non-perturbative merger
phase of the binary black-hole coalescence problem. This enables us to
\emph{coherently} search for all three stages of the coalescence of
non-spinning binary black holes using a single template bank. Taking our
motivation from these results, we propose a family of template waveforms which
can model the inspiral, merger, and ring-down stages of the coalescence of
non-spinning binary black holes that follow quasi-circular inspiral. This
two-dimensional template family is explicitly parametrized by the physical
parameters of the binary. We show that the template family is not only
\emph{effectual} in detecting the signals from black hole coalescences, but
also \emph{faithful} in estimating the parameters of the binary. We compare the
sensitivity of a search (in the context of different ground-based
interferometers) using all three stages of the black hole coalescence with
other template-based searches which look for individual stages separately. We
find that the proposed search is significantly more sensitive than other
template-based searches for a substantial mass-range, potentially bringing
about remarkable improvement in the event-rate of ground-based interferometers.
As part of this work, we also prescribe a general procedure to construct
interpolated template banks using non-spinning black hole waveforms produced by
numerical relativity.Comment: A typo fixed in Eq.(B11
AMR, stability and higher accuracy
Efforts to achieve better accuracy in numerical relativity have so far
focused either on implementing second order accurate adaptive mesh refinement
or on defining higher order accurate differences and update schemes. Here, we
argue for the combination, that is a higher order accurate adaptive scheme.
This combines the power that adaptive gridding techniques provide to resolve
fine scales (in addition to a more efficient use of resources) together with
the higher accuracy furnished by higher order schemes when the solution is
adequately resolved. To define a convenient higher order adaptive mesh
refinement scheme, we discuss a few different modifications of the standard,
second order accurate approach of Berger and Oliger. Applying each of these
methods to a simple model problem, we find these options have unstable modes.
However, a novel approach to dealing with the grid boundaries introduced by the
adaptivity appears stable and quite promising for the use of high order
operators within an adaptive framework
Waveform propagation in black hole spacetimes: evaluating the quality of numerical solutions
We compute the propagation and scattering of linear gravitational waves off a
Schwarzschild black hole using a numerical code which solves a generalization
of the Zerilli equation to a three dimensional cartesian coordinate system.
Since the solution to this problem is well understood it represents a very good
testbed for evaluating our ability to perform three dimensional computations of
gravitational waves in spacetimes in which a black hole event horizon is
present.Comment: 13 pages, RevTeX, to appear in Phys. Rev.
Finding Apparent Horizons in Dynamic 3D Numerical Spacetimes
We have developed a general method for finding apparent horizons in 3D
numerical relativity. Instead of solving for the partial differential equation
describing the location of the apparent horizons, we expand the closed 2D
surfaces in terms of symmetric trace--free tensors and solve for the expansion
coefficients using a minimization procedure. Our method is applied to a number
of different spacetimes, including numerically constructed spacetimes
containing highly distorted axisymmetric black holes in spherical coordinates,
and 3D rotating, and colliding black holes in Cartesian coordinates.Comment: 19 pages, 13 figures, LaTex, to appear in Phys. Rev. D. Minor changes
mad
- …