89 research outputs found
IllinoisGRMHD: An Open-Source, User-Friendly GRMHD Code for Dynamical Spacetimes
In the extreme violence of merger and mass accretion, compact objects like
black holes and neutron stars are thought to launch some of the most luminous
outbursts of electromagnetic and gravitational wave energy in the Universe.
Modeling these systems realistically is a central problem in theoretical
astrophysics, but has proven extremely challenging, requiring the development
of numerical relativity codes that solve Einstein's equations for the
spacetime, coupled to the equations of general relativistic (ideal)
magnetohydrodynamics (GRMHD) for the magnetized fluids. Over the past decade,
the Illinois Numerical Relativity (ILNR) Group's dynamical spacetime GRMHD code
has proven itself as a robust and reliable tool for theoretical modeling of
such GRMHD phenomena. However, the code was written "by experts and for
experts" of the code, with a steep learning curve that would severely hinder
community adoption if it were open-sourced. Here we present IllinoisGRMHD,
which is an open-source, highly-extensible rewrite of the original
closed-source GRMHD code of the ILNR Group. Reducing the learning curve was the
primary focus of this rewrite, with the goal of facilitating community
involvement in the code's use and development, as well as the minimization of
human effort in generating new science. IllinoisGRMHD also saves computer time,
generating roundoff-precision identical output to the original code on
adaptive-mesh grids, but nearly twice as fast at scales of hundreds to
thousands of cores.Comment: 37 pages, 6 figures, single column. Matches published versio
Prompt Electromagnetic Transients from Binary Black Hole Mergers
Binary black hole (BBH) mergers provide a prime source for current and future
interferometric GW observatories. Massive BBH mergers may often take place in
plasma-rich environments, leading to the exciting possibility of a concurrent
electromagnetic (EM) signal observable by traditional astronomical facilities.
However, many critical questions about the generation of such counterparts
remain unanswered. We explore mechanisms that may drive EM counterparts with
magnetohydrodynamic simulations treating a range of scenarios involving
equal-mass black-hole binaries immersed in an initially homogeneous fluid with
uniform, orbitally aligned magnetic fields. We find that the time development
of Poynting luminosity, which may drive jet-like emissions, is relatively
insensitive to aspects of the initial configuration. In particular, over a
significant range of initial values, the central magnetic field strength is
effectively regulated by the gas flow to yield a Poynting luminosity of
, with BBH mass
scaled to and ambient density . We also calculate the
direct plasma synchrotron emissions processed through geodesic ray-tracing.
Despite lensing effects and dynamics, we find the observed synchrotron flux
varies little leading up to merger.Comment: 22 pages, 21 figures; additional reference + clarifying text added to
match published versio
Accretion disks around binary black holes of unequal mass: GRMHD simulations near decoupling
We report on simulations in general relativity of magnetized disks onto black
hole binaries. We vary the binary mass ratio from 1:1 to 1:10 and evolve the
systems when they orbit near the binary-disk decoupling radius. We compare
(surface) density profiles, accretion rates (relative to a single, non-spinning
black hole), variability, effective -stress levels and luminosities as
functions of the mass ratio. We treat the disks in two limiting regimes: rapid
radiative cooling and no radiative cooling. The magnetic field lines clearly
reveal jets emerging from both black hole horizons and merging into one common
jet at large distances. The magnetic fields give rise to much stronger shock
heating than the pure hydrodynamic flows, completely alter the disk structure,
and boost accretion rates and luminosities. Accretion streams near the horizons
are among the densest structures; in fact, the 1:10 no-cooling evolution
results in a refilling of the cavity. The typical effective temperature in the
bulk of the disk is yielding characteristic thermal frequencies . These systems are
thus promising targets for many extragalactic optical surveys, such as LSST,
WFIRST, and PanSTARRS.Comment: 29 pages, 23 captioned figures, 3 tables, submitted to PR
GiRaFFE: An Open-Source General Relativistic Force-Free Electrodynamics Code
We present GiRaFFE, the first open-source general relativistic force-free
electrodynamics (GRFFE) code for dynamical, numerical-relativity generated
spacetimes. GiRaFFE adopts the strategy pioneered by McKinney and modified by
Paschalidis and Shapiro to convert a GR magnetohydrodynamic (GRMHD) code into a
GRFFE code. In short, GiRaFFE exists as a modification of IllinoisGRMHD, a
user-friendly, open-source, dynamical-spacetime GRMHD code. Both GiRaFFE and
IllinoisGRMHD leverage the Einstein Toolkit's highly-scalable infrastructure to
make possible large-scale simulations of magnetized plasmas in strong,
dynamical spacetimes on adaptive-mesh refinement (AMR) grids. We demonstrate
that GiRaFFE passes a large suite of both flat and curved-spacetime code tests
passed by a number of other state-of-the-art GRFFE codes, and is thus ready for
production-scale simulations of GRFFE phenomena of key interest to relativistic
astrophysics.Comment: 23 pages, 4 figures. Consistent with published versio
Improved Moving Puncture Gauge Conditions for Compact Binary Evolutions
Robust gauge conditions are critically important to the stability and
accuracy of numerical relativity (NR) simulations involving compact objects.
Most of the NR community use the highly robust---though
decade-old---moving-puncture (MP) gauge conditions for such simulations. It has
been argued that in binary black hole (BBH) evolutions adopting this gauge,
noise generated near adaptive-mesh-refinement (AMR) boundaries does not
converge away cleanly with increasing resolution, severely limiting
gravitational waveform accuracy at computationally feasible resolutions. We
link this noise to a sharp (short-wavelength), initial outgoing gauge wave
crossing into progressively lower resolution AMR grids, and present
improvements to the standard MP gauge conditions that focus on stretching,
smoothing, and more rapidly settling this outgoing wave. Our best gauge choice
greatly reduces gravitational waveform noise during inspiral, yielding less
fluctuation in convergence order and lower waveform phase and
amplitude errors at typical resolutions. Noise in other physical quantities of
interest is also reduced, and constraint violations drop by more than an order
of magnitude. We expect these improvements will carry over to simulations of
all types of compact binary systems, as well as other +1 formulations of
gravity for which MP-like gauge conditions can be chosen.Comment: 25 pages, 16 figures, 2 tables. Matches published versio
- …