Related papers: GRHayL: a modern, infrastructure-agnostic, extensi…
Apparent horizon (AH) finders are essential for characterizing black holes and excising their interiors in numerical relativity (NR) simulations. However, open-source AH finders to date are tightly coupled to individual NR codes. We…
We introduce an extension to the AthenaK code for general-relativistic magnetohydrodynamics (GRMHD) in dynamical spacetimes using a 3+1 conservative Eulerian formulation. Like the fixed-spacetime GRMHD solver, we use standard finite-volume…
We model the structure and evolution of black hole accretion disks, and their neighboring regions, using numerical simulations. The numerics is governed by the equations of general relativistic magneto-hydrodynamics (GRMHD). In particular,…
Astrophysical plasmas in relativistic spacetimes, such as black hole accretion flows, are often weakly collisional and require kinetic modeling to capture non-local transport and particle acceleration. However, the extreme scale separation…
A new code and methodology are introduced for solving the general relativistic magnetohydrodynamic (GRMHD) equations in fixed background spacetimes using time-explicit, finite-volume discretization. The code has options for solving the…
We present a new numerical code, X-ECHO, for general relativistic magnetohydrodynamics (GRMHD) in dynamical spacetimes. This is aimed at studying astrophysical situations where strong gravity and magnetic fields are both supposed to play an…
There is great interest in numerical relativity simulations involving matter due to the likelihood that binary compact objects involving neutron stars will be detected by gravitational wave observatories in the coming years, as well as to…
The force-free limit of magnetohydrodynamics (MHD) is often a reasonable approximation to model black hole and neutron star magnetospheres. We describe a general relativistic force-free (GRFFE) formulation that allows general relativistic…
We study the evolution of two fiducial configurations for binary neutron stars using two different general relativistic hydrodynamics(GRHD), distributed adaptive mesh codes. One code, HAD, has for many years been used to study mergers of…
Hot, diffuse, relativistic plasmas such as sub-Eddington black hole accretion flows are expected to be collisionless, yet are commonly modeled as a fluid using ideal general relativistic magnetohydrodynamics (GRMHD). Dissipative effects…
HII regions powered by ionizing radiation from massive stars drive the dynamical evolution of the interstellar medium. Fast radiative transfer methods for incorporating photoionization effects are thus essential in astrophysical…
Black hole-neutron star mergers, together with binary neutron star mergers, are key laboratories for neutron star physics. They enable us to probe merger dynamics imprinted in gravitational waves and potential electromagnetic counterparts.…
We investigate the long-term evolution of black hole accretion disks formed in neutron star mergers. These disks expel matter that contributes to an $r$-process kilonova, and can produce relativistic jets powering short gamma-ray bursts.…
Global, general relativistic magnetohydrodynamic (GRMHD) simulations of nonradiative, magnetized disks are widely used to model accreting black holes. We have performed a convergence study of GRMHD models computed with HARM3D. The models…
We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code, RAISHIN, using a conservative, high resolution shock-capturing scheme. The numerical fluxes are calculated using the Harten, Lax, & van Leer…
Almost a hundred compact binary mergers have been detected via gravitational waves by the LIGO-Virgo-KAGRA collaboration in the past few years providing us with a significant amount of new information on black holes and neutron stars. In…
The integration of distributed renewable energy sources and the multi-domain behaviours inside the cyber-physical energy system (smart grids) draws up major challenges. Their validation and roll out requires careful assessment, in term of…
An implicit algorithm for solving the equations of general relativistic hydrodynamics in conservative form in three-dimensional axi-symmetry is presented. This algorithm is a direct extension of the pseudo-Newtonian implicit radiative…
One of the main challenges in the numerical modeling of binary neutron-star (BNS) mergers is long-term simulations of the post-merger remnant over timescales of the order of seconds. When this modeling includes all the aspects of complex…
This work summarizes a program intended to unify three burgeoning branches of the high-energy astrophysics of relativistic jets: general relativistic magnetohydrodynamic (GRMHD) simulations of ever-increasing dynamical range, the…