We construct a simple theoretical model to investigate how entrainment gradually erodes high-speed FR II jets. This process is described by embedding a mixing-layer model developed originally to describe FR I objects in a self-similar model for the lobe structure of classical FR II sources. Following the classical FR II models, we assume that the lobe is dominated by the particles injected from the central jet. The entrainment produces a boundary shear layer which acts at the interface between the dense central jet and the less denser surrounding lobe, and the associated erosion of the jet places interesting limits on the maximum size of FR II sources. The model shows that this limit depends mainly on the initial bulk velocity of the relativistic jet triggered. The bulk velocities of FR IIs suggested by our model are in good agreement with that obtained from direct pc-scale observations on ordinary radio galaxies and quasars. Finally, we discuss how FR IIs may evolve into FR Is upon reaching their maximum, entrainment-limited sizes.

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. Abridged

It is generally accepted that the heating of gas in clusters of galaxies by active galactic nuclei (AGN) is a form of feedback. Feedback is required to ensure a long term, sustainable balance between heating and cooling. This work investigates the impact of proportional stochastic feedback on the energy balance in the intracluster medium. Using a generalised analytical model for a cluster atmosphere, it is shown that an energy equilibrium can be reached exponentially quickly. Applying the tools of stochastic calculus it is demonstrated that the result is robust with regard to the model parameters, even though they affect the amount of variability in the system.

We present an analytical model for jets in Fanaroff & Riley Class I (FRI) 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 FRI source 3C31 and show that there is a self-consistent solution. We derive the jet power, together with the variations of mass flux and 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 VLA observations. Our prediction for the shape of the laminar core can be tested with higher-resolution imaging.

There is growing evidence that the active galactic nuclei (AGN) associated with the central elliptical galaxy in clusters of galaxies are playing an important role in the evolution of the intracluster medium (ICM) and clusters themselves. We use high resolution three-dimensional simulations to study the interaction of the cavities created by AGN outflows (bubbles) with the ambient ICM. The gravitational potential of the cluster is modelled using the observed temperature and density profiles of the Virgo cluster. We demonstrate the importance of the hydrodynamical Kutta-Zhukovsky forces associated with the vortex ring structure of the bubbles, and discuss possible effects of diffusive processes on their evolution.

We present a phenomenological model of the dynamics of buoyant bubbles in the atmosphere of a cluster of galaxies. The derived equations describe velocity, size, mass, temperature and density of the buoyant bubbles as functions of time based on several simple approximations. The constructed model is then used to interpret results of a numerical experiment of heating of the cluster core with buoyant bubbles in a hydrodynamical approximation (i.e. in the absence of magnetic fields, viscosity, and thermal diffusion). Based on the model parameters we discuss possible limitations of the numerical treatment of the problem, and highlight the main physical processes that govern the dynamics of bubbles in the intracluster medium.

The results of hydrodynamic simulations of the Virgo and Perseus clusters suggest that thermal conduction is not responsible for the observed temperature and density profiles. As a result it seems that thermal conduction occurs at a much lower level than the Spitzer value. Comparing cavity enthalpies to the radiative losses within the cooling radius for seven clusters suggests that some clusters are probably heated by sporadic, but extremely powerful, AGN outflows interspersed between more frequent but lower power outflows.

A roughly constant temperature over a wide range of densities is maintained in molecular clouds through radiative heating and cooling. An isothermal equation of state is therefore frequently employed in molecular cloud simulations. However, the dynamical processes in molecular clouds include shock waves, expansion waves, cooling induced collapse and baroclinic vorticity, all incompatible with the assumption of a purely isothermal flow. Here, we incorporate an energy equation including all the important heating and cooling rates and a simple chemical network into simulations of three-dimensional, hydrodynamic, decaying turbulence. This allows us to test the accuracy of the isothermal assumption by directly comparing a model run with the modified energy equation to an isothermal model. We compute an extreme case in which the initial turbulence is sufficiently strong to dissociate much of the gas and alter the specific heat ratio. The molecules then reform as the turbulence weakens. We track the true specific heat ratio as well as its effective value. We analyse power spectra, vorticity and shock structures, and discuss scaling laws for decaying turbulence. We derive some limitations to the isothermal approximation for simulations of the interstellar medium using simple projection techniques. Overall, even given the extreme conditions, we find that an isothermal flow provides an adequate physical and observational description of many properties. The main exceptions revealed here concern behaviour directly related to the high temperature zones behind the shock waves.

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 motived 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 clasical 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 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 require duty cycles of roughly $10^{6}-10^{7}$ yrs 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.

Thermal conduction has been suggested as a possible mechanism by which sufficient energy is supplied to the central regions of galaxy clusters to balance the effect of radiative cooling. Here we present the results of a simulated, high-resolution, 3-d Virgo cluster for different values of thermal conductivity (1, 1/10, 1/100, 0 times the full Spitzer value). Starting from an initially isothermal cluster atmosphere we allow the cluster to evolve freely over timescales of roughly $ 1.3-4.7 \times 10^{9} $ yrs. Our results show that thermal conductivity at the Spitzer value can increase the central ICM radiative cooling time by a factor of roughly 3.6. In addition, for larger values of thermal conductvity the simulated temperature and density profiles match the observations significantly better than for the lower values. However, no physically meaningful value of thermal conductivity was able to postpone the cooling catastrophe (characterised by a rapid increase in the central density) for longer than a fraction of the Hubble time nor explain the absence of a strong cooling flow in the Virgo cluster today. We also calculate the effective adiabatic index of the cluster gas for both simulation and observational data and compare the values with theoretical expectations. Using this method it appears that the Virgo cluster is being heated in the cluster centre by a mechanism other than thermal conductivity. Based on this and our simulations it is also likely that the thermal conductvity is suppressed by a factor of at least 10 and probably more. Thus, we suggest that thermal conductvity, if present at all, has the effect of slowing down the evolution of the ICM, by radiative cooling, but only by a factor of a few.

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.

We present the results from three dimensional hydrodynamical simulations of decaying high-speed 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, to 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, due to the radiative cooling. The total thermal energy, initially raised by the introduction of turbulence, decays only a little slower 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.