5 research outputs found
Locating Boosted Kerr and Schwarzschild Apparent Horizons
We describe a finite-difference method for locating apparent horizons and
illustrate its capabilities on boosted Kerr and Schwarzschild black holes. Our
model spacetime is given by the Kerr-Schild metric. We apply a Lorentz boost to
this spacetime metric and then carry out a 3+1 decomposition. The result is a
slicing of Kerr/Schwarzschild in which the black hole is propagated and Lorentz
contracted. We show that our method can locate distorted apparent horizons
efficiently and accurately.Comment: Submitted to Physical Review D. 12 pages and 22 figure
Extended Lifetime in Computational Evolution of Isolated Black Holes
Solving the 4-d Einstein equations as evolution in time requires solving
equations of two types: the four elliptic initial data (constraint) equations,
followed by the six second order evolution equations. Analytically the
constraint equations remain solved under the action of the evolution, and one
approach is to simply monitor them ({\it unconstrained} evolution).
The problem of the 3-d computational simulation of even a single isolated
vacuum black hole has proven to be remarkably difficult. Recently, we have
become aware of two publications that describe very long term evolution, at
least for single isolated black holes. An essential feature in each of these
results is {\it constraint subtraction}. Additionally, each of these approaches
is based on what we call "modern," hyperbolic formulations of the Einstein
equations. It is generally assumed, based on computational experience, that the
use of such modern formulations is essential for long-term black hole
stability. We report here on comparable lifetime results based on the much
simpler ("traditional") - formulation.
We have also carried out a series of {\it constrained} 3-d evolutions of
single isolated black holes. We find that constraint solution can produce
substantially stabilized long-term single hole evolutions. However, we have
found that for large domains, neither constraint-subtracted nor constrained
- evolutions carried out in Cartesian coordinates admit
arbitrarily long-lived simulations. The failure appears to arise from features
at the inner excision boundary; the behavior does generally improve with
resolution.Comment: 20 pages, 6 figure
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
Generic Tracking of Multiple Apparent Horizons with Level Flow
We report the development of the first apparent horizon locator capable of
finding multiple apparent horizons in a ``generic'' numerical black hole
spacetime. We use a level-flow method which, starting from a single arbitrary
initial trial surface, can undergo topology changes as it flows towards
disjoint apparent horizons if they are present. The level flow method has two
advantages: 1) The solution is independent of changes in the initial guess and
2) The solution can have multiple components. We illustrate our method of
locating apparent horizons by tracking horizon components in a short
Kerr-Schild binary black hole grazing collision.Comment: 13 pages including figures, submitted to Phys Rev
Three-dimensional adaptive evolution of gravitational waves in numerical relativity
Adaptive techniques are crucial for successful numerical modeling of
gravitational waves from astrophysical sources such as coalescing compact
binaries, since the radiation typically has wavelengths much larger than the
scale of the sources. We have carried out an important step toward this goal,
the evolution of weak gravitational waves using adaptive mesh refinement in the
Einstein equations. The 2-level adaptive simulation is compared with unigrid
runs at coarse and fine resolution, and is shown to track closely the features
of the fine grid run.Comment: REVTeX, 7 pages, including three figures; submitted to Physical
Review