Study accretion and ejection using a new GPU-accelerated GRMHD code

We study disks and jets in various accretion states (SANE and MAD) using novel, GPU-accelerated general-relativistic magneto-hydrodynamic (GR-MHD) code which we developed, based on HARM. This code, written in CUDA-c and uses OpenMP to parallelize multi-GPU setups, allows high resolution simulations of accretion disks and the formation and structure of jets without the need of multi-node supercomputer infrastructure. A 256³ simulation is well within the reach of an Nvidia DGX-V100 server, with the computation being a factor about 100 times faster if only the CPU was used. We use this code to examine several disk structures, wind and jet properties in the MAD and SANE states. In the MAD state, we find that the magnetic flux threading the horizon mostly depends on the spin of the BH. This implies that the jet structure and power are strong functions of the spin, with non-spinning BHs have the widest jets.