We run hydrodynamical simulations of a 2D isothermal non self-gravitating inviscid gas flowing in a rigidly rotating externally imposed potential formed by only two components: a monopole and a quadrupole. We explore systematically the effects of var
ying the quadrupole while keeping fixed the monopole and discuss the consequences for the interpretation of longitude-velocity diagrams in the Milky Way. We find that the gas flow can constrain the quadrupole of the potential and the characteristics of the bar that generates it. The exponential scale length of the bar must be at least $1.5rm, kpc$. The strength of the bar is also constrained. Our global interpretation favours a pattern speed of $Omega=40,rm km s^{-1} {kpc}^{-1}$. We find that for most observational features, there exist a value of the parameters that matches each individual feature well, but is difficult to reproduce all the important features at once. Due to the intractably high number of parameters involved in the general problem, quantitative fitting methods that can run automatic searches in parameter space are necessary.
Spiral arms that emerge from the ends of a galactic bar are important in interpreting observations of our and external galaxies. It is therefore important to understand the physical mechanism that causes them. We find that these spiral arms can be un
derstood as kinematic density waves generated by librations around underlying ballistic closed orbits. This is even true in the case of a strong bar, provided the librations are around the appropriate closed orbits and not around the circular orbits that form the basis of the epicycle approximation. An important consequence is that it is a potentials orbital structure that determines whether a bar should be classified as weak or strong, and not crude estimates of the potentials deviation from axisymmetry.
We use a Cartesian grid to simulate the flow of gas in a barred Galactic potential and investigate the effects of varying the sound speed in the gas and the resolution of the grid. For all sound speeds and resolutions, streamlines closely follow clos
ed orbits at large and small radii. At intermediate radii shocks arise and the streamlines shift between two families of closed orbits. The point at which the shocks appear and the streamlines shift between orbit families depends strongly on sound speed and resolution. For sufficiently large values of these two parameters, the transfer happens at the cusped orbit as hypothesised by Binney et al. over two decades ago. For sufficiently high resolutions the flow downstream of the shocks becomes unsteady. If this unsteadiness is physical, as appears to be the case, it provides a promising explanation for the asymmetry in the observed distribution of CO.
We present a new method for fitting simple hydrodynamical models to the (l,v) distribution of atomic and molecular gas observed in the Milky Way. The method works by matching features found in models and observations. It is based on the assumption th
at the large-scale features seen in (l,v) plots, such as ridgelines and the terminal velocity curve, are influenced primarily by the underlying large-scale Galactic potential and are only weakly dependent on local ISM heating and cooling processes. In our scheme one first identifies by hand the features in the observations: this only has to be done once. We describe a procedure for automatically extracting similar features from simple hydrodynamical models and quantifying the distance between each models features and the observations. Application to models of the Galactic Bar region (|l|<30deg) shows that our feature-fitting method performs better than chi^2 or envelope distances at identifying the correct underlying galaxy model.
A re-investigation of the gravothermal catastrophe is presented. By means of a linear perturbation analysis, we study the dynamical stability of a spherical self-gravitating isothermal fluid of finite volume and find that the conditions for the onset
of the gravothermal catastrophe, under different external conditions, coincide with those obtained from thermodynamical arguments. This suggests that the gravothermal catastrophe may reduce to Jeans instability, rediscovered in an inhomogeneous framework. We find normal modes and frequencies for the fluid system and show that instability develops on the dynamical time scale. We then discuss several related issues. In particular: (1) For perturbations at constant total energy and constant volume, we introduce a simple heuristic term in the energy budget to mimic the role of binaries. (2) We outline the analysis of the two-component case and show how linear perturbation analysis can be carried out also in this more complex context in a relatively straightforward way. (3) We compare the behavior of the fluid model with that of the collisionless sphere. In the collisionless case the instability seems to disappear, which is at variance with the linear Jeans stability analysis in the homogeneous case; we argue that a key ingredient to understand the difference (a spherical stellar system is expected to undergo the gravothermal catastrophe only in the presence of some collisionality, which suggests that the instability is dissipative and not dynamical) lies in the role of the detailed angular momentum in a collisionless system. Finally, we briefly comment on the meaning of the Boltzmann entropy and its applicability to the study of the dynamics of self- gravitating inhomogeneous gaseous systems.
