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.

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.

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.