Related papers: Towards an exascale code for GRMHD on dynamical sp…

In this era of precision cosmology, a detailed physical understanding on the evolution of cosmic baryons is required. Cosmic magnetic fields, though still poorly understood, may represent an important component in the global cosmic energy…

We have recently constructed a numerical code that evolves a spherically symmetric spacetime using a hyperbolic formulation of Einstein's equations. For the case of a Schwarzschild black hole, this code works well at early times, but…

We evolve a scalar field in a fixed Kerr-Schild background geometry to test simple $(3+1)$-dimensional algorithms for singularity excision. We compare both centered and upwind schemes for handling the shift (advection) terms, as well as…

We present a new numerical tool for solving the special relativistic ideal MHD equations that is based on the combination of the following three key features: (i) a one-step ADER discontinuous Galerkin (DG) scheme that allows for an…

Numerical radiation-hydrodynamics (RHD) for non-relativistic flows is a challenging problem because it encompasses processes acting over a very broad range of timescales, and where the relative importance of these processes often varies by…

Recent years have seen a significant progress in the development of general relativistic codes for the numerical solution of the equations of magnetohydrodynamics in spacetimes with high and dynamical curvature. These codes are valuable…

With the forthcoming VLBI images of Sgr A* and M87, simulations of accretion flows onto black holes acquire a special importance to aid with the interpretation of the observations and to test the predictions of different accretion…

We present MHDuet, an open source evolution code for general relativistic magnetohydrodynamics with neutrino transport. The code solves the full set of Einstein equations coupled to a relativistic, magnetized fluid with an M1 neutrino…

Spherically symmetric (1D) black-hole spacetimes are considered as a test for numerical relativity. A finite difference code, based in the hyperbolic structure of Einstein's equations with the harmonic slicing condition is presented.…

General Relativistic Magnetohydrodynamics (GRMHD) simulations are an indispensable tool in studying accretion onto compact objects. The Event Horizon Telescope (EHT) frequently uses libraries of ideal GRMHD simulations to interpret…

We use rigorous techniques from numerical analysis of hyperbolic equations in bounded domains to construct stable finite-difference schemes for Numerical Relativity, in particular for their use in black hole excision. As an application, we…

This work presents a novel methodology for deriving stationary and axially symmetric solutions to Einstein field equations using the 1+3 tetrad formalism. This approach reformulates the Einstein equations into first order scalar equations,…

Fallback supernovae and the collapsar scenario for long-gamma ray burst and hypernovae have received considerable interest as pathways to black-hole formation and extreme transient events. Consistent simulations of these scenarios require a…

Based on the conformal energy theorem we prove the uniqueness theorem for static higher dimensional electrically and magnetically charged black holes being the solution of Einstein (n-2)-gauge forms equations of motion. Black hole spacetime…

In the harmonic description of general relativity, the principle part of Einstein's equations reduces to 10 curved space wave equations for the componenets of the space-time metric. We present theorems regarding the stability of several…

We present a novel method of magnetohydrodynamics (MHD) within the smoothed particle hydrodynamics scheme using the Geometric Density average force expression (GDSPH). GDSPH has recently been shown to reduce the leading order errors and…

We present a new python/SymPy based code generator for producing executable numerical expressions for partial differential equations in AMReX-based applications. We demonstrate the code generator capabilities for the case of $3+1$ ADM…

We present ErasureHead, a new approach for distributed gradient descent (GD) that mitigates system delays by employing approximate gradient coding. Gradient coded distributed GD uses redundancy to exactly recover the gradient at each…

We have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to IllinoisGRMHD, an open-source GRMHD code for dynamical spacetimes in the Einstein Toolkit. These…

We present a new pseudospectral code, bamps, for numerical relativity written with the evolution of collapsing gravitational waves in mind. We employ the first order generalized harmonic gauge formulation. The relevant theory is reviewed…