cuHARM: A New GPU-accelerated GRMHD Code and Its Application to ADAF Disks

We introduce a new GPU-accelerated general relativistic magnetohydrodynamic 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-resolution simulations of accretion disks and the formation and structure of jets without a need for multinode supercomputer infrastructure. A 256³ simulation is well within the reach of an Nvidia DGX-V100 server, with the computation being about 10 times faster than if only the CPU is used. We use this code to examine several disk structures all in the SANE state. We find that (i) increasing the magnetic field in the SANE state does not affect the mass accretion rate; (ii) simultaneously increasing the disk size and the magnetic field, while keeping the ratio of energies fixed, leads to the destruction of the jet once the magnetic flux through the horizon decreases below a certain limit (this demonstrates that the existence of a jet is a function of the initial mass distribution, and not of the initial intensity of the magnetic field, since the magnetorotational instability dictates the evolution of the magnetic field); and (iii) the structure of the jet is a weak function of the adiabatic index of the gas, with relativistic gas tending to have a wider jet.