A new version of the Athena++ AMR MHD code implemented using the Kokkos library for performance-portability will be described. The code is being used for a variety of problems related to accretion onto (and merger of) compact objects. This talk will focus on results from general relativistic radiation MHD models of accretion in luminous systems such as X-ray binaries and AGN. The models allow observationally important quantities such as the radiative efficiency, variability, and kinetic feedback in winds to be measured over a range of luminosities, from sub-Eddington to highly super-Eddington. A new “cyclic-zoom” approach for modeling such flows over very long time scales will also be briefly introduced.
