Related papers: Towards an exascale code for GRMHD on dynamical sp…
Accreting supermassive black hole binaries are powerful multimessenger sources emitting both gravitational and EM radiation. Understanding the accretion dynamics of these systems and predicting their distinctive EM signals is crucial to…
This paper is concerned exclusively with axisymmetric spacetimes. We want to develop reductions of Einstein's equations which are suitable for numerical evolutions. We first make a Kaluza-Klein type dimensional reduction followed by an ADM…
We explore various aspects of dynamical black holes defined by a future outer trapping horizon in $n(\ge 5)$-dimensional Einstein-Gauss-Bonnet gravity. In the present paper, we assume that the spacetime has symmetries corresponding to the…
We present a new code for astrophysical magneto-hydrodynamics specifically designed and optimized for high performance and scaling on modern and future supercomputers. We describe a novel hybrid OpenMP/MPI programming model that emerged…
This thesis focuses on the application of numerical relativity methods to the solutions of problems in strong gravity. Our goal is the study of mergers of compact objects in the strong field regime where non-linear dynamics manifest and…
We derive the precise late-time asymptotics for solutions to the wave equation on extremal Reissner-Nordstr\"om black holes and explicitly express the leading-order coefficients in terms of the initial data. Our method is based on purely…
Non-continuous "jumps" of Apparent Horizons occur generically in 3+1 (binary) black hole evolutions. The dynamical trapping horizon framework suggests a spacetime picture in which these "Apparent Horizon jumps" are understood as spatial…
We present general relativistic radiation magnetohydrodynamics (GRRMHD) simulations of super-Eddington accretion flows around supermassive black holes (SMBHs) which may apply to tidal disruption events (TDEs). We perform long duration…
The next generation of Department of Energy supercomputers will be capable of exascale computation. For these machines, far more computation will be possible than that which can be saved to disk. As a result, users will be unable to rely on…
X-ray binaries (XRBs) display spectral state transitions that are accompanied by substantial changes in the hardness, luminosity, and structure of the accretion flow. We developed a GPU-accelerated cooling toolkit for general relativistic…
We propose two new alternative numerical schemes to solve the coupled Einstein-Euler equations in the Generalized Harmonic formulation. The first one is a finite difference (FD) Central Weighted Essentially Non-Oscillatory (CWENO) scheme on…
We present novel asymptotically global AdS$_4$ solutions, constructed by turning on a dipolar differential rotation at the conformal boundary. At fixed energy and boundary profile, we find two different geometries: a horizonless spacetime,…
Detailed observations of phenomena involving black holes, be it via gravitational waves or more traditional electromagnetic means, can probe the strong field regime of the gravitational interaction. The prediction of features in such…
We present a fourth-order accurate finite volume method for the solution of ideal magnetohydrodynamics (MHD). The numerical method combines high-order quadrature rules in the solution of semi-discrete formulations of hyperbolic conservation…
The architecture of Exascale computing facilities, which involves millions of heterogeneous processing units, will deeply impact on scientific applications. Future astrophysical HPC applications must be designed to make such computing…
We have developed a parallel cosmological hydrodynamic simulation code designed for the study of formation and evolution of cosmological structures. The gravitational force is calculated using the TreePM method and the hydrodynamics is…
Black hole (BH) accretion flows and jets are dynamic hot relativistic magnetized plasma flows whose radiative opacity can significantly affect flow structure and behavior. We describe a numerical scheme, tests, and an astrophysically…
In this paper we construct a class of hairy static black holes of higher dimensional Einstein-Skyrme theories with the cosmological constant $\Lambda \le 0$ whose scalar is an $SU(2)$ valued field. The spacetime is set to be conformal to $…
We present a new algorithm for evolving orbiting black-hole binaries that does not require excision or a corotating shift. Our algorithm is based on a novel technique to handle the singular puncture conformal factor. This system, based on…
We show explicitly that traversable wormholes requiring exotic matter in 4-dimensions nevertheless have acceptable stress-tensors obeying reasonable energy conditions in higher dimensions if the wormholes are regarded as being embedded in…