This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions,
and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the codes parallel performance, and discuss the Enzo collaborations code development methodology.
Galaxies moving through the intracluster medium (ICM) of a cluster of galaxies can lose gas via ram pressure stripping. This stripped gas forms a tail behind the galaxy which is potentially observable. In this paper, we carry out hydrodynamical simul
ations of a galaxy undergoing stripping with a focus on the gas properties in the wake and their observational signatures. We include radiative cooling in an adaptive hydrocode in order to investigate the impact of a clumpy, multi-phase interstellar medium. We find that including cooling results in a very different morphology for the gas in the tail, with a much wider range of temperatures and densities. The tail is significantly narrower in runs with radiative cooling, in agreement with observed wakes. In addition, we make detailed predictions of H I, Halpha and X-ray emission for the wake, showing that we generally expect detectable H I and Halpha signatures, but no observable X-ray emission (at least for our chosen ram-pressure strength and ICM conditions). We find that the relative strength of the Halpha diagnostic depends somewhat on our adopted minimum temperature floor (below which we set cooling to zero to mimic physics processes not included in the simulation).
UV radiation from early astrophysical sources could have a large impact on subsequent star formation in nearby protogalaxies. Here we study the radiative feedback from the first, short-lived stars using hydrodynamical simulations with transient UV ba
ckgrounds (UVBs) and persistent Lyman-Werner backgrounds (LWBs) of varying intensity. We extend our prior work in Mesinger et al. (2006), by studying a more typical region whose proto-galaxies form at lower redshifts, z~13-20, in the epoch likely preceding the bulk of reionization. We confirm our previous results that feedback in the relic HII regions resulting from such transient radiation, is itself transient. Feedback effects dwindle away after ~30% of the Hubble time, and the same critical specific intensity of J_UV~0.1 x 10^{-21} ergs/s/cm^2/Hz/sr separates positive and negative feedback regimes. Additionally, we discover a second episode of eventual positive feedback in halos which have not yet collapsed when their progenitor regions were exposed to the transient UVB. This eventual positive feedback appears in all runs, regardless of the strength of the UVB. However, this feedback regime is very sensitive to the presence of Lyman-Werner radiation, and notable effects disappear under fairly modest background intensities of J_LW>10^{-3} x 10^{-21} ergs/s/cm^2/Hz/sr. We conclude that UV radiative feedback in relic HII regions, although a complicated process, seems unlikely to have a major impact on the progress of cosmological reionization, provided that present estimates of the lifetime and luminosity of a PopIII star are accurate. More likely is that the build-up of the LWB ultimately governs the feedback strength until a persistent UV background can be established. [abridged]
Modelling global disc galaxies is a difficult task which has previously resulted in the small scale physics of the interstellar medium being greatly simplified. In this talk, I compare simulations of galaxies with different ISM properties to determin
e the importance of the ISM structure in the star formation properties of the disc.
We perform a set of non--radiative cosmological simulations of a preheated intracluster medium in which the entropy of the gas was uniformly boosted at high redshift. The results of these simulations are used first to test the current analytic techni
ques of preheating via entropy input in the smooth accretion limit. When the unmodified profile is taken directly from simulations, we find that this model is in excellent agreement with the results of our simulations. This suggests that preheated efficiently smoothes the accreted gas, and therefore a shift in the unmodified profile is a good approximation even with a realistic accretion history. When we examine the simulation results in detail, we do not find strong evidence for entropy amplification, at least for the high-redshift preheating model adopted here. In the second section of the paper, we compare the results of the preheating simulations to recent observations. We show -- in agreement with previous work -- that for a reasonable amount of preheating, a satisfactory match can be found to the mass-temperature and luminosity-temperature relations. However -- as noted by previous authors -- we find that the entropy profiles of the simulated groups are much too flat compared to observations. In particular, while rich clusters converge on the adiabatic self--similar scaling at large radius, no single value of the entropy input during preheating can simultaneously reproduce both the core and outer entropy levels. As a result, we confirm that the simple preheating scenario for galaxy cluster formation, in which entropy is injected universally at high redshift, is inconsistent with observations.
