Spectral fits to M87 core data from radio to hard x-ray are generated via a specially selected software suite, comprised of the HARM GRMHD accretion disk model and a 2D Monte Carlo radiation transport code. By determining appropriate parameter changes necessary to fit x-ray quiescent and flaring behavior of M87's core, we assess the reasonableness of various flaring mechanisms. This shows that an accretion disk model of M87's core out to 28 GM/c^2 can describe the inner emissions. High spin rates show GRMHD-driven polar outflow generation, without citing an external jet model. Our results favor accretion rate changes as the dominant mechanism of x-ray flux and index changes, with variations in density of approximately 20% necessary to scale between the average x-ray spectrum and flaring or quiescent spectra. The best fit parameters are black hole spin a/M > 0.8 and maximum accretion flow density n <= 3x10^7 cm^-3, equivalent to horizon accretion rates between m_dot = M_dot/M_dot_Edd ~ 2x10^-6 and 1x10^-5 (with M_dot_Edd defined assuming a radiative efficiency eta = 0.1). These results demonstrate that the immediate surroundings of M87's core are appropriate to explain observed x-ray variability.

We present results of models of the physical space and parameters of the accretion disk of Sagittarius A*, as well as simulations of its emergent spectrum. This begins with HARM, a 2D general relativistic magneto-hydrodynamic (GRMHD) model, specifically set up to evolve the space around a black hole. Data from HARM are then fed into a 2D Monte-Carlo (MC) code which generates and tracks emitted photons, allowing for absorption and scattering before they escape the volume.

We present results of simulations of the spectrum of the accretion flow onto the supermassive black hole in our Galactic Centre, Sagittarius A*, generated with a coupling of Monte-Carlo (MC) radiation and general relativistic magnetohydrodynamic (GRMHD) codes. In our modeling, we use the 2D HARM GRMHD code to first model the physical parameters of the disk, then feed its results into our 2D MC photon transport code. We will discuss results obtained which fit radio, IR, and Chandra-obtained flaring or quiescent x-ray data points, as well as the validity of the amount of scaling of input parameters (density, temperature, and magnetic field) required to fit these points. HARM output will be used to suggest whether the scaling is within reasonable limits.