Hydrodynamics and stability of galactic cooling flows

A. Kritsuk,H. Boehringer,E. Mueller

DOI: https://doi.org/10.1046/j.1365-8711.1998.01897.x

1998-02-07

Abstract:Using numerical techniques we studied the global stability of cooling flows in giant elliptical galaxies. As an initial equilibrium state we choose the hydrostatic gas recycling model (Kritsuk 1996). Non-equilibrium radiative cooling, stellar mass loss, heating by type Ia supernovae, distributed mass deposition, and thermal conductivity are included. Although the recycling model reproduces the basic X-ray observables, it appears to be unstable with respect to the development of inflow or outflow. In spherically symmetry the inflows are subject to a central cooling catastrophe, while the outflows saturate in a form of a subsonic galactic wind. Two-dimensional axisymmetric random velocity perturbations of the equilibrium model trigger the onset of a cooling catastrophe, which develops in an essentially non-spherical way. The simulations show a patchy pattern of mass deposition and the formation of hollow gas jets, which penetrate through the outflow down to the galaxy core. The X-ray observables of such a hybrid gas flow mimic those of the equilibrium recycling model, but the gas temperature exhibits a central depression. The mass deposition rate M_dot consists of two contributions of similar size: (i) a hydrostatic one resembling that of the equilibrium model, and (ii) a dynamical one which is related to the jets and is more concentrated to the centre. For a model galaxy, like NGC 4472, our 2D simulations predict M_dot = 2 M_sun/yr within the cooling radius for the advanced non-linear stage of the instability. We discuss the implications of these results to H_alpha nebulae and star formation in cooling flow galaxies and emphasize the need for high-resolution 3D simulations.

Astrophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is about the global stability of galactic cooling flows. Specifically, the author uses numerical techniques to study the global stability of cooling flows in giant elliptical galaxies. The initial equilibrium state is chosen as the hydrostatic gas - circulation model (Kritsuk 1996). Factors such as non - equilibrium radiative cooling, stellar mass loss, Type Ia supernova heating, distributed mass deposition, and thermal conduction are considered in the study. Although the circulation model can reproduce the basic X - ray observations, it seems to be unstable in developing inflows or outflows.

### Main problems
1. **Stability of cooling flows**:
- The paper explores the global stability of the hydrostatic equilibrium gas configuration in the absence of effective thermal conduction. The research shows that in the absence of effective thermal conduction, the equilibrium state is unstable. For example, for a 10% positive density perturbation, an inflow will be generated on the central cooling time scale \( t_c\approx0.5t_s \) (see Figure 1), while the same negative density perturbation will trigger an outflow (see Figure 2).

2. **Development of inflows and outflows**:
- The positive gas density perturbation first breaks the balance between cooling and heating in the core region. Gravity immediately begins to accelerate the gas inflow. As a result, a dense approximately isothermal gas core with a temperature of about \( 10^6 \) K is formed, and this core is condensation - unstable. The further development of instability depends on the detailed structure of the inflow velocity profile. Condensation first occurs at the radius where \( \nabla\cdot\mathbf{v} \) has a local minimum. In the example shown in Figure 1, this occurs almost simultaneously at the inner boundary and at \( r\approx0.5 \) kpc. There are two distinct minima in the temperature profile (with a temperature of about \( 10^4 \) K), and there are also obvious dips in the density distribution due to high condensation efficiency.
- For \( \eta\in[0, 0.1] \), during this transition stage, the mass loss rate \( \dot{M}(r) \) is 2 to 4 \( M_\odot \) yr\(^{-1}\) within the cooling radius \( r_{\text{cool}} \). Most of the condensate is deposited within a radius of 1 to 2 kpc. The surface brightness in the core region is 10 times higher than the equilibrium value.

3. **Effect of two - dimensional axisymmetric random velocity perturbations**:
- Two - dimensional axisymmetric random velocity perturbations trigger the occurrence of a cooling catastrophe, which develops in a non - spherical manner. The simulations show a patchy pattern of mass deposition and hollow gas jets that penetrate the outflow and reach the galactic core. The X - ray observable properties of these mixed gas flows are similar to those of the equilibrium circulation model, but the gas temperature drops in the central region. The mass deposition rate \( \dot{M} \) consists of two contributions of similar magnitudes: (i) the hydrostatic part, which is similar to the equilibrium model; (ii) the dynamic part, which is related to the jets and is more concentrated in the center.

### Conclusions
- This study emphasizes the need for high - resolution three - dimensional simulations to more comprehensively understand the complex physical processes of cooling flows and their impact on Hα nebulae and star formation in cooling - flow galaxies.

Through these studies, the author hopes to better understand the physical processes in cooling - flow galaxies and provide a theoretical basis for future research.