We investigate the effect of three important processes by which AGN-blown bubbles transport material: drift, wake transport and entrainment. The first of these, drift, occurs because a buoyant bubble pushes aside the adjacent material, giving rise to a net upward displacement of the fluid behind the bubble. For a spherical bubble, the mass of upwardly displaced material is roughly equal to half the mass displaced by the bubble, and should be ~ 10^{7-9} solar masses depending on the local ICM and bubble parameters. We show that in classical cool core clusters, the upward displacement by drift may be a key process in explaining the presence of filaments behind bubbles. A bubble also carries a parcel of material in a region at its rear, known as the wake. The mass of the wake is comparable to the drift mass and increases the average density of the bubble, trapping it closer to the cluster centre and reducing the amount of heating it can do during its ascent. Moreover, material dropping out of the wake will also contribute to the trailing filaments. Mass transport by the bubble wake can effectively prevent the build-up of cool material in the central galaxy, even if AGN heating does not balance ICM cooling. Finally, we consider entrainment, the process by which ambient material is incorporated into the bubble. AbridgedComment: Accepted for publication in MNRAS. 17 pages, 4 figures, 2 tables. Formatted for letter paper and adjusted author affiliations

The buoyant rise of hot plasma bubbles inflated by AGN outflows in galaxy clusters can heat the cluster gas and thereby compensate radiative energy losses of this material. Numerical simulations of this effect often show the complete disruption of the bubbles followed by the mixing of the bubble material with the surrounding cluster gas due to fluid instabilities on the bubble surface. This prediction is inconsistent with the observations of apparently coherent bubble structures in clusters. We derive a general description in the linear regime of the growth of instabilities on the surface between two fluids under the influence of a gravitational field, viscosity, surface tension provided by a magnetic field and relative motion of the two fluids with respect to each other. We demonstrate that Kelvin-Helmholtz instabilities are always suppressed, if the fluids are viscous. They are also suppressed in the inviscid case for fluids of very different mass densities. We show that the effects of shear viscosity as well as a magnetic fields in the cluster gas can prevent the growth of Rayleigh-Taylor instabilities on relevant scale lengths. R-T instabilities on pc-scales are suppressed even if the kinematic viscosity of the cluster gas is reduced by two orders of magnitude compared to the value given by Spitzer for a fully ionised, unmagnetised gas. Similarly, magnetic fields exceeding a few microG result in an effective surface tension preventing the disruption of bubbles. For more massive clusters, instabilities on the bubble surface grow faster. This may explain the absence of thermal gas in the north-west bubble observed in the Perseus cluster compared to the apparently more disrupted bubbles in the Virgo cluster.Comment: 12 pages, 5 figures, accepted by MNRA

In recent years, evidence has accumulated suggesting that the gas in galaxy clusters is heated by non‐gravitational processes. Here, we calculate the heating rates required to maintain a physically motivated mass flow rate, in a sample of seven galaxy clusters. We employ the spectroscopic mass deposition rates as an observational input along with temperature and density data for each cluster. On energetic grounds, we find that thermal conduction could provide the necessary heating for A2199, Perseus, A1795 and A478. However, the suppression factor of the classical Spitzer value is a different function of radius for each cluster. Based on the observations of plasma bubbles, we also calculate the duty cycles for each active galactic nucleus (AGN), in the absence of thermal conduction, which can provide the required energy input. With the exception of Hydra‐A, it appears that each of the other AGNs in our sample requires duty cycles of roughly 106–107 yr to provide their steady‐state heating requirements. If these duty cycles are unrealistic, this may imply that many galaxy clusters must be heated by very powerful Hydra‐A type events interspersed between more frequent smaller scale outbursts. The suppression factors for the thermal conductivity required for combined heating by AGN and thermal conduction are generally acceptable. However, these suppression factors still require ‘fine‐tuning’ of the thermal conductivity as a function of radius. As a consequence of this work, we present the AGN duty cycle as a cooling flow diagnostic.

We present an analytical model for jets in Fanaroff & Riley Class I (FR I) radio galaxies, in which an initially laminar, relativistic flow is surrounded by a shear layer. We apply the appropriate conservation laws to constrain the jet parameters, starting the model where the radio emission is observed to brighten abruptly. We assume that the laminar flow fills the jet there and that pressure balance with the surroundings is maintained from that point outwards. Entrainment continuously injects new material into the jet and forms a shear layer, which contains material from both the environment and the laminar core. The shear layer expands rapidly with distance until finally the core disappears, and all of the material is mixed into the shear layer. Beyond this point, the shear layer expands in a cone and decelerates smoothly. We apply our model to the well‐observed FR I source 3C 31 and show that there is a self‐consistent solution. We derive the jet power, together with the variations of mass flux and entrainment rate with distance from the nucleus. The predicted variation of bulk velocity with distance in the outer parts of the jets is in good agreement with model fits to Very Large Array observations. Our prediction for the shape of the laminar core can be tested with higher‐resolution imaging.

We present the results from three-dimensional hydrodynamical simulations of decaying highspeed turbulence in dense molecular clouds. We compare our results, which include a detailed cooling function, molecular hydrogen chemistry and a limited C and O chemistry, with those previously obtained for decaying isothermal turbulence.After an initial phase of shock formation, power-law decay regimes are uncovered, as in the isothermal case. We find that the turbulence decays faster than in the isothermal case because the average Mach number remains higher, owing to the radiative cooling. The total thermal energy, initially raised by the introduction of turbulence, decays only a little more slowly than the kinetic energy.We discover that molecule reformation, as the fast turbulence decays, is several times faster than that predicted for a non-turbulent medium. This is caused by moderate speed shocks which sweep through a large fraction of the volume, compressing the gas and dust. Through reformation, the molecular density and molecular column appear as complex patterns of filaments, clumps and some diffuse structure. In contrast, the molecular fraction has a wider distribution of highly distorted clumps and copious diffuse structure, so that density and molecular density are almost identically distributed during the reformation phase. We conclude that molecules form in swept-up clumps but effectively mix throughout via subsequent expansions and compressions.