We describe and test an updated version of radiation-hydrodynamics (RHD) in the RAMSES code, that includes three new features: i) radiation pressure on gas, ii) accurate treatment of radiation diffusion in an unresolved optically thick medium, and ii
i) relativistic corrections that account for Doppler effects and work done by the radiation to first order in v/c. We validate the implementation in a series of tests, which include a morphological assessment of the M1 closure for the Eddington tensor in an astronomically relevant setting, dust absorption in a optically semi-thick medium, direct pressure on gas from ionising radiation, convergence of our radiation diffusion scheme towards resolved optical depths, correct diffusion of a radiation flash and a constant luminosity radiation, and finally, an experiment from Davis et al. of the competition between gravity and radiation pressure in a dusty atmosphere, and the formation of radiative Rayleigh-Taylor instabilities. With the new features, RAMSES-RT can be used for state-of-the-art simulations of radiation feedback from first principles, on galactic and cosmological scales, including not only direct radiation pressure from ionising photons, but also indirect pressure via dust from multi-scattered IR photons reprocessed from higher-energy radiation, both in the optically thin and thick limits.
In a recent study (Martizzi et al. 2012), we used cosmological simulations to show that active galactic nuclei (AGN) feedback on the gas distribution in clusters of galaxies can be important in determining the spatial distribution of stars and dark m
atter in the central regions of these systems. The hierarchical assembly of dark matter, baryons and black holes obscures the physical mechanism behind the restructuring process. Here we use idealized simulations to follow the response of a massive dark matter halo as we feed the central black hole with a controlled supply of cold gas. This removes most of the complexity taking place in the cosmological simulations that may have biased our previous study. We confirm our previous results: gas heated and expelled from the central regions of the halo by AGN feedback can return after cooling; repeated cycles generate gravitational potential fluctuations responsible for irreversible modifications of the dark matter mass profile. The main result is the expulsion of large amounts of baryons and dark matter from the central regions of the halo. According to the work presented here, outflow induced fluctuations represent the only mechanism able to efficiently create dark matter cores in clusters of galaxies.
We have implemented non-ideal Magneto-Hydrodynamics (MHD) effects in the Adaptive Mesh Refinement (AMR) code RAMSES, namely ambipolar diffusion and Ohmic dissipation, as additional source terms in the ideal MHD equations. We describe in details how w
e have discretized these terms using the adaptive Cartesian mesh, and how the time step is diminished with respect to the ideal case, in order to perform a stable time integration. We have performed a large suite of test runs, featuring the Barenblatt diffusion test, the Ohmic diffusion test, the C-shock test and the Alfven wave test. For the latter, we have performed a careful truncation error analysis to estimate the magnitude of the numerical diffusion induced by our Godunov scheme, allowing us to estimate the spatial resolution that is required to address non-ideal MHD effects reliably. We show that our scheme is second-order accurate, and is therefore ideally suited to study non-ideal MHD effects in the context of star formation and molecular cloud dynamics.
