878 research outputs found
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
Electrostatic Patch Effect in Cylindrical Geometry. I. Potential and Energy between Slightly Non-Coaxial Cylinders
We study the effect of any uneven voltage distribution on two close
cylindrical conductors with parallel axes that are slightly shifted in the
radial and by any length in the axial direction. The investigation is
especially motivated by certain precision measurements, such as the Satellite
Test of the Equivalence Principle (STEP). By energy conservation, the force can
be found as the energy gradient in the vector of the shift, which requires
determining potential distribution and energy in the gap. The boundary value
problem for the potential is solved, and energy is thus found to the second
order in the small transverse shift, and to lowest order in the gap to cylinder
radius ratio. The energy consists of three parts: the usual capacitor part due
to the uniform potential difference, the one coming from the interaction
between the voltage patches and the uniform voltage difference, and the energy
of patch interaction, entirely independent of the uniform voltage. Patch effect
forces and torques in the cylindrical configuration are derived and analyzed in
the next two parts of this work.Comment: 26 pages, 1 Figure. Submitted to Classical and Quantum Gravit
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
Numerical Relativity Using a Generalized Harmonic Decomposition
A new numerical scheme to solve the Einstein field equations based upon the
generalized harmonic decomposition of the Ricci tensor is introduced. The
source functions driving the wave equations that define generalized harmonic
coordinates are treated as independent functions, and encode the coordinate
freedom of solutions. Techniques are discussed to impose particular gauge
conditions through a specification of the source functions. A 3D, free
evolution, finite difference code implementing this system of equations with a
scalar field matter source is described. The second-order-in-space-and-time
partial differential equations are discretized directly without the use first
order auxiliary terms, limiting the number of independent functions to
fifteen--ten metric quantities, four source functions and the scalar field.
This also limits the number of constraint equations, which can only be enforced
to within truncation error in a numerical free evolution, to four. The
coordinate system is compactified to spatial infinity in order to impose
physically motivated, constraint-preserving outer boundary conditions. A
variant of the Cartoon method for efficiently simulating axisymmetric
spacetimes with a Cartesian code is described that does not use interpolation,
and is easier to incorporate into existing adaptive mesh refinement packages.
Preliminary test simulations of vacuum black hole evolution and black hole
formation via scalar field collapse are described, suggesting that this method
may be useful for studying many spacetimes of interest.Comment: 18 pages, 6 figures; updated to coincide with journal version, which
includes some expanded discussions and a new appendix with a stability
analysis of a simplified problem using the same discretization scheme
described in the pape
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
Adaptive Mesh Refinement for Characteristic Grids
I consider techniques for Berger-Oliger adaptive mesh refinement (AMR) when
numerically solving partial differential equations with wave-like solutions,
using characteristic (double-null) grids. Such AMR algorithms are naturally
recursive, and the best-known past Berger-Oliger characteristic AMR algorithm,
that of Pretorius & Lehner (J. Comp. Phys. 198 (2004), 10), recurses on
individual "diamond" characteristic grid cells. This leads to the use of
fine-grained memory management, with individual grid cells kept in
2-dimensional linked lists at each refinement level. This complicates the
implementation and adds overhead in both space and time.
Here I describe a Berger-Oliger characteristic AMR algorithm which instead
recurses on null \emph{slices}. This algorithm is very similar to the usual
Cauchy Berger-Oliger algorithm, and uses relatively coarse-grained memory
management, allowing entire null slices to be stored in contiguous arrays in
memory. The algorithm is very efficient in both space and time.
I describe discretizations yielding both 2nd and 4th order global accuracy.
My code implementing the algorithm described here is included in the electronic
supplementary materials accompanying this paper, and is freely available to
other researchers under the terms of the GNU general public license.Comment: 37 pages, 15 figures (40 eps figure files, 8 of them color; all are
viewable ok in black-and-white), 1 mpeg movie, uses Springer-Verlag svjour3
document class, includes C++ source code. Changes from v1: revised in
response to referee comments: many references added, new figure added to
better explain the algorithm, other small changes, C++ code updated to latest
versio
New Upper Limit of Terrestrial Equivalence Principle Test for Rotating Extended Bodies
Improved terrestrial experiment to test the equivalence principle for
rotating extended bodies is presented, and a new upper limit for the violation
of the equivalence principle is obtained at the level of 1.6, which is limited by the friction of the rotating gyroscope. It
means the spin-gravity interaction between the extended bodies has not been
observed at this level.Comment: 4 page
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
Status of NINJA: the Numerical INJection Analysis project
The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise
Hyperboloidal foliations and scri-fixing
We discuss a gauge choice which allows us to avoid the introduction of
artificial timelike outer boundaries in numerical studies of test fields based
on a 3+1 decomposition of asymptotically flat background spacetimes. The main
idea is to include null infinity in the computational domain by conformally
compactifying the metric on hyperboloidal foliations and fixing the spatial
coordinate location of null infinity, i.e. scri-fixing. We construct such
coordinates explicitly on Minkowski, Schwarzschild and Kerr spacetimes.Comment: 14 pages, 14 figures. Published versio
- …