We analyzed the feedback of AGN jets on cooling flow clusters using three-dimensional AMR hydrodynamic simulations. We studied the interaction of the jet with the intracluster medium and creation of low X-ray emission cavities (Bubbles) in cluster plasma. The distribution of energy input by the jet into the system was quantified in its different forms, i.e. internal, kinetic and potential. We find that the energy associated with the bubbles, (pV + gamma pV/(gamma-1)), accounts for less than 10 percent of the jet energy.

We examine the distribution of injected energy in three-dimensional, adaptive-grid simulations of the heating of cooling flows. We show that less than 10 percent of the injected energy goes into bubbles. Consequently, the energy input from the nucleus is underestimated by a factor of order 6 when it is taken to be given by PVgamma/(gamma-1), where P and V are the pressure and volume of the bubble, and gamma the ratio of principal specific heats.

Three-dimensional hydrodynamical simulations are used to investigate the structural stability of cooling flows that are episodically heated by jets from a central AGN. The radial profile of energy deposition is controlled by (a) the power of the jets, and (b) the pre-outburst density profile. A delay in the ignition of the jets causes more powerful jets to impact on a more centrally concentrated medium. The net effect is a sufficient increase in the central concentration of energy deposition to cause the post-outburst density profile to be less centrally concentrated than that of an identical cluster in which the outburst happened earlier and was weaker. These results suggest that the density profiles of cooling flows oscillate around an attracting profile, thus explaining why cooling flows are observed to have similar density profiles. The possibility is raised that powerful FR II systems are ones in which this feedback mechanism has broken down and a runaway growth of the source parameters has occurred.

Active galactic nuclei are clearly heating gas in `cooling flows'. The effectiveness and spatial distribution of the heating are controversial. We use three-dimensional simulations on adaptive grids to study the impact on a cooling flow of weak, subrelativistic jets. The simulations show cavities and vortex rings as in the observations. The cavities are fast-expanding dynamical objects rather than buoyant bubbles as previously modelled, but shocks still remain extremely hard to detect with X-rays. At late times the cavities turn into overdensities that strongly excite the cluster's g-modes. These modes damp on a long timescale. Radial mixing is shown to be an important phenomenon, but the jets weaken the metallicity gradient only very near the centre. The central entropy density is modestly increased by the jets. We use a novel algorithm to impose the jets on the simulations.