34 research outputs found
The [N II] Kinematics of R Aquarii
We report a kinematic study of the symbiotic star system R Aqr derived from [N H]lambda 6584 emission observations with a Fabry-Perot imaging spectrometer. The [N II] spatial structure of the R Aqr jet, first observed circa 1977, and surrounding hourglass-shaped nebulosity, due to an explosion approximately 660 years ago, are derived from 41 velocity planes spaced at approximately 12 km/s intervals. Fabry-Perot imagery shows the elliptical nebulosity comprising the waist of the hourglass shell is consistent with a circular ring expanding radially at 55 km/s as seen at an inclination angle, i approximately 70 deg. Fabry-Perot imagery shows the two-sided R Aqr jet is collimated flow in opposite directions. The intensity-velocity structure of the strong NE jet component is shown in contrast to the amorphous SW jet component. We offer a idealized schematic model for the R Aqr jet motion which results in a small-scale helical structure forming around a larger-scale helical path. The implications of such a jet model are discussed. We present a movie showing a side-by-side comparison of the spatial structure of the model and the data as a function of the 41 velocity planes
New, efficient, and accurate high order derivative and dissipation operators satisfying summation by parts, and applications in three-dimensional multi-block evolutions
We construct new, efficient, and accurate high-order finite differencing
operators which satisfy summation by parts. Since these operators are not
uniquely defined, we consider several optimization criteria: minimizing the
bandwidth, the truncation error on the boundary points, the spectral radius, or
a combination of these. We examine in detail a set of operators that are up to
tenth order accurate in the interior, and we surprisingly find that a
combination of these optimizations can improve the operators' spectral radius
and accuracy by orders of magnitude in certain cases. We also construct
high-order dissipation operators that are compatible with these new finite
difference operators and which are semi-definite with respect to the
appropriate summation by parts scalar product. We test the stability and
accuracy of these new difference and dissipation operators by evolving a
three-dimensional scalar wave equation on a spherical domain consisting of
seven blocks, each discretized with a structured grid, and connected through
penalty boundary conditions.Comment: 16 pages, 9 figures. The files with the coefficients for the
derivative and dissipation operators can be accessed by downloading the
source code for the document. The files are located in the "coeffs"
subdirector
Characteristic extraction in numerical relativity: binary black hole merger waveforms at null infinity
The accurate modeling of gravitational radiation is a key issue for
gravitational wave astronomy. As simulation codes reach higher accuracy,
systematic errors inherent in current numerical relativity wave-extraction
methods become evident, and may lead to a wrong astrophysical interpretation of
the data. In this paper, we give a detailed description of the
Cauchy-characteristic extraction technique applied to binary black hole
inspiral and merger evolutions to obtain gravitational waveforms that are
defined unambiguously, that is, at future null infinity. By this method we
remove finite-radius approximations and the need to extrapolate data from the
near zone. Further, we demonstrate that the method is free of gauge effects and
thus is affected only by numerical error. Various consistency checks reveal
that energy and angular momentum are conserved to high precision and agree very
well with extrapolated data. In addition, we revisit the computation of the
gravitational recoil and find that finite radius extrapolation very well
approximates the result at \scri. However, the (non-convergent) systematic
differences to extrapolated data are of the same order of magnitude as the
(convergent) discretisation error of the Cauchy evolution hence highlighting
the need for correct wave-extraction.Comment: 41 pages, 8 figures, 2 tables, added references, fixed typos. Version
matches published version
Matching post-Newtonian and numerical relativity waveforms: systematic errors and a new phenomenological model for non-precessing black hole binaries
We present a new phenomenological gravitational waveform model for the
inspiral and coalescence of non-precessing spinning black hole binaries. Our
approach is based on a frequency domain matching of post-Newtonian inspiral
waveforms with numerical relativity based binary black hole coalescence
waveforms. We quantify the various possible sources of systematic errors that
arise in matching post-Newtonian and numerical relativity waveforms, and we use
a matching criteria based on minimizing these errors; we find that the dominant
source of errors are those in the post-Newtonian waveforms near the merger. An
analytical formula for the dominant mode of the gravitational radiation of
non-precessing black hole binaries is presented that captures the phenomenology
of the hybrid waveforms. Its implementation in the current searches for
gravitational waves should allow cross-checks of other inspiral-merger-ringdown
waveform families and improve the reach of gravitational wave searches.Comment: 22 pages, 11 figure
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
Initial data transients in binary black hole evolutions
We describe a method for initializing characteristic evolutions of the
Einstein equations using a linearized solution corresponding to purely outgoing
radiation. This allows for a more consistent application of the characteristic
(null cone) techniques for invariantly determining the gravitational radiation
content of numerical simulations. In addition, we are able to identify the {\em
ingoing} radiation contained in the characteristic initial data, as well as in
the initial data of the 3+1 simulation. We find that each component leads to a
small but long lasting (several hundred mass scales) transient in the measured
outgoing gravitational waves.Comment: 18 pages, 4 figure
Exact boundary conditions in numerical relativity using multiple grids: scalar field tests
Cauchy-Characteristic Matching (CCM), the combination of a central 3+1 Cauchy
code with an exterior characteristic code connected across a time-like
interface, is a promising technique for the generation and extraction of
gravitational waves. While it provides a tool for the exact specification of
boundary conditions for the Cauchy evolution, it also allows to follow
gravitational radiation all the way to infinity, where it is unambiguously
defined.
We present a new fourth order accurate finite difference CCM scheme for a
first order reduction of the wave equation around a Schwarzschild black hole in
axisymmetry. The matching at the interface between the Cauchy and the
characteristic regions is done by transfering appropriate characteristic/null
variables. Numerical experiments indicate that the algorithm is fourth order
convergent. As an application we reproduce the expected late-time tail decay
for the scalar field.Comment: 14 pages, 5 figures. Included changes suggested by referee
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
The Samurai Project: verifying the consistency of black-hole-binary waveforms for gravitational-wave detection
We quantify the consistency of numerical-relativity black-hole-binary
waveforms for use in gravitational-wave (GW) searches with current and planned
ground-based detectors. We compare previously published results for the
mode of the gravitational waves from an equal-mass
nonspinning binary, calculated by five numerical codes. We focus on the 1000M
(about six orbits, or 12 GW cycles) before the peak of the GW amplitude and the
subsequent ringdown. We find that the phase and amplitude agree within each
code's uncertainty estimates. The mismatch between the modes
is better than for binary masses above with respect to
the Enhanced LIGO detector noise curve, and for masses above
with respect to Advanced LIGO, Virgo and Advanced Virgo. Between the waveforms
with the best agreement, the mismatch is below . We find that
the waveforms would be indistinguishable in all ground-based detectors (and for
the masses we consider) if detected with a signal-to-noise ratio of less than
, or less than in the best cases.Comment: 17 pages, 9 figures. Version accepted by PR
A pilgrimage to gravity on GPUs
In this short review we present the developments over the last 5 decades that
have led to the use of Graphics Processing Units (GPUs) for astrophysical
simulations. Since the introduction of NVIDIA's Compute Unified Device
Architecture (CUDA) in 2007 the GPU has become a valuable tool for N-body
simulations and is so popular these days that almost all papers about high
precision N-body simulations use methods that are accelerated by GPUs. With the
GPU hardware becoming more advanced and being used for more advanced algorithms
like gravitational tree-codes we see a bright future for GPU like hardware in
computational astrophysics.Comment: To appear in: European Physical Journal "Special Topics" : "Computer
Simulations on Graphics Processing Units" . 18 pages, 8 figure