Related papers: grim: A Flexible, Conservative Scheme for Relativi…
Performing accurate large eddy simulations in compressible, turbulent magnetohydrodynamics is more challenging than in non-magnetized fluids due to the complex interplay between kinetic, magnetic and internal energy at different scales.…
Based on a 14-moment closure for non-resistive (general-) relativistic viscous plasmas, we describe a new numerical scheme that is able to handle all first-order dissipative effects (heat conduction, bulk and shear viscosities), as well the…
We perform three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations of a near-maximally spinning black hole (spin parameter, a = 0.998) with varying initial magnetic field geometries, systematically exploring the…
We present 3D general relativistic magnetohydrodynamic(GRMHD) simulations of zero angular momentum accretion around a rapidly rotating black hole, modified by the presence of initially uniform magnetic fields. We consider serveral angles…
The general relativistic magnetohydrodynamic (GRMHD) simulations are widely used to study accretion disk and jet dynamics around a black hole. Despite strong observational evidences for intrinsically nonlinear behavior, the interpretations…
This paper describes the development and testing of a general relativistic magnetohydrodynamic (GRMHD) code to study ideal MHD in the fixed background of a Kerr black hole. The code is a direct extension of the hydrodynamic code of Hawley,…
A generalization of implicit conservative numerics to multiple dimensions requires advanced concepts of tensor analysis and differential geometry and hence a more thorough dedication to mathematical fundamentals than maybe expected at first…
This paper investigates the application of mesh adaptation techniques in the Non-Ideal Compressible Fluid Dynamic (NICFD) regime, a region near the vapor-liquid saturation curve where the flow behavior significantly departs from the ideal…
MHD turbulence is likely to play an important role in several astrophysical scenarios where the magnetic Reynolds is very large. Numerically, these cases can be studied efficiently by means of Large Eddy Simulations, in which the…
A widely accepted picture of an accretion flow in the luminous soft spectral state of X-ray binary systems is a geometrically thin disc structure much like the classic analytic solution of Shakura \& Sunyaev. Although the analytic models…
General-relativistic magnetohydrodynamic (GRMHD) simulations are an important tool to study a variety of astrophysical systems such as neutron star mergers, core-collapse supernovae, and accretion onto compact objects. A conservative GRMHD…
In this proceeding, we provide a novel approach to study the General Relativistic Magnetohydrodynamic (GRMHD) accretion flows around rotating black holes (BHs). In doing so, we choose a sub-Keplerian distribution of angular momentum of the…
We introduce a new GPU-accelerated general-relativistic magneto-hydrodynamic (GR-MHD) code based on HARM which we call cuHARM. The code is written in CUDA-C and uses OpenMP to parallelize multi-GPU setups. Our code allows us to run high…
We present a self-contained overview of GR-Athena++, a general-relativistic magnetohydrodynamics (GRMHD) code, that incorporates treatment of dynamical space-time, based on the recent work of (Daszuta+, 2021)[49] and (Cook+, 2023)[45].…
Observational efforts to image the immediate environment of a black hole at the scale of the event horizon benefit from the development of efficient imaging codes that are capable of producing synthetic data, which may be compared with…
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…
Recent observations by the Event Horizon Telescope (EHT) of supermassive black holes M87* and Sgr A* offer valuable insights into their spacetime properties and astrophysical conditions. Utilizing a library of model images (~2 million for…
We present an adaptive-order positivity-preserving conservative finite-difference scheme that allows a high-order solution away from shocks and discontinuities while guaranteeing positivity and robustness at discontinuities. This is…
Observations are providing increasingly detailed quantitative information about the accretion flows that power such high energy systems as X-ray binaries and active galactic nuclei. Analytic models of such systems must rely on assumptions…
We present the combination of dynamical accretion model based on 3D GRMHD simulations and general relativistic (GR) polarized radiative transfer. We write down the formalism of and perform the GR ray-tracing of cyclo-synchrotron radiation…