Related papers: IllinoisGRMHD: An Open-Source, User-Friendly GRMHD…
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…
In this work we present a proof of concept of CUDA-capable, resistive, multi-fluid models of relativistic magnetohydrodynamics (RMHD). Resistive and multi-fluid codes for simulating models of RMHD suffer from stiff source terms, so it is…
The powerful high-energy phenomena typically encountered in astrophysics invariably involve physical engines, like neutron stars and black hole accretion disks, characterized by a combination of highly magnetized plasmas, strong…
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…
High-energy phenomena in astrophysics involve quite generally a combination of relativistic motions and strong gravity. The simultaneous solution of Einstein equations and General Relativistic MHD equations is thus necessary to model with…
We present a new numerical code designed to solve the Einstein field equations for axisymmetric spacetimes. The long term goal of this project is to construct a code that will be capable of studying many problems of interest in axisymmetry,…
Here we extend the exploration of significantly super-Chandrasekhar magnetised white dwarfs by numerically computing axisymmetric stationary equilibria of differentially rotating magnetised polytropic compact stars in general relativity…
The numerical simulation of turbulence and flows in almost ideal, large-Reynolds-number astrophysical plasmas motivates the implementation of almost conservative MHD computer codes. They should efficiently calculate, use highly parallelized…
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…
The advent of high-fidelity imaging of supermassive black holes calls for efficient and robust data-analysis methods. In this work, we use $\texttt{Jipole}$, a differentiable, $\texttt{ipole}$-based radiative transfer code, to enable…
Current state-of-the-art simulations of accretion flows onto black holes require a significant level of numerical sophistication, in order to allow the three-dimensional modelling of relativistic magnetised plasma in a regime of strong…
Merging binary neutron stars (BNSs) represent the ultimate targets for multimessenger astronomy, being among the most promising sources of gravitational waves (GWs), and, at the same time, likely accompanied by a variety of electromagnetic…
Accretion of magnetized gas on compact astrophysical objects such as black holes has been successfully modeled using general relativistic magnetohydrodynamic (GRMHD) simulations. These simulations have largely been performed in the Kerr…
According to the theory of general relativity, the relative acceleration of masses generates gravitational radiation. Although gravitational radiation has not yet been detected, it is believed that extremely violent cosmic events, such as…
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…
Given the compact binary evolution problem of numerical relativity, in the finite-difference, block-based, adaptive mesh refinement context, choices must be made on how evolved fields are to be discretized. In GR-Athena++, the space-time…
Next-generation gravitational wave detectors such as Cosmic Explorer, the Einstein Telescope, and LISA, demand highly accurate and extensive gravitational wave (GW) catalogs to faithfully extract physical parameters from observed signals.…
This paper presents a fast, economical particle-multiple-mesh N-body code optimized for large-N modelling of collisionless dynamical processes, such as black-hole wandering or bar-halo interactions, occurring within isolated galaxies. The…
We present version 1.0 of our Lagrangian Numerical Relativity code SPHINCS_BSSN. This code evolves the full set of Einstein equations, but contrary to other Numerical Relativity codes, it evolves the matter fluid via Lagrangian particles in…
The applications of numerical relativity to cosmology are on the rise, contributing insight into such cosmological problems as structure formation, primordial phase transitions, gravitational-wave generation, and inflation. In this paper, I…