Glimm-Godunov's Method for Cosmic-ray-hydrodynamics

A numerical method for integrating the equations describing a dynamically coupled system made of a fluid and cosmic-rays is developed. In smooth flows the effect of CR pressure is accounted for by modification of the characteristic equations and the energy exchange between cosmic-rays and the fluid, due to diffusive processes in configuration and momentum space, is modeled with a flux conserving method. Provided the shock acceleration efficiency as a function of the upstream conditions and shock Mach number, we show that the Riemann solver can be modified to take into account the cosmic-ray mediation without having to resolve the cosmic-ray induced substructure. Shocks are advanced with Glimm's method which preserves their discontinuous character without any smearing, thus allowing to maintain self-consistency in the shock solutions. In smooth flows either Glimm's or a higher order Godunov's method can be applied, with the latter producing better results when approximations are introduced in the Riemann solver.Comment: 32 pages, 4 figs, JCP in press, improved description of boundary conditions at high momenta, references updated, version matching the one accepted for publicatio

Simulating molecular cloud regulated star formation in Galaxies

This thesis is primarily concerned with understanding the process of galaxy formation via the simulation of the interstellar medium, star formation and supernova feedback. In order to probe galaxy formation it is necessary that we first obtain a basic knowledge of the cosmological framework in which we are working. Therefore in chapter 1 we provide a brief overview of the salient features of the current cosmological paradigm in addition to discussing in some detail the physics of the interstellar medium. In chapter 2 we focus on the numerical methods necessary to perform accurate cosmological simulations. We begin by providing an overview of the different simulation techniques currently in use in the field before performing comparisons of two simulation codes that work via two completely different methods. We then introduce a code for generating high-resolution initial conditions for the simulation of individual objects and investigate the numerical effects of mass resolution in cosmological simulation. In chapter 3 we describe a statistical model of the interstellar medium, in which the conversion of cooling gas to stars in the multiphase interstellar medium is governed by the rate at which molecular clouds are formed and destroyed. In the model, clouds form from thermally unstable ambient gas and get destroyed by star formation, feedback and thermal conduction. In chapter 4 this statistical model is applied to the simulation of isolated disk galaxies. We show that it naturally produces a multiphase medium with cold clouds, a warm disk and hot supernova bubbles. We illustrate this by evolving an isolated Milky Way like galaxy. Many observed properties of disk galaxies are reproduced well, including the molecular cloud mass spectrum, the molecular fraction as a function of radius, the Schmidt law, the stellar density profile and the appearance of a galactic fountain. Finally in chapter 5 we perform an investigation into more dynamic situations, namely the evolution of gravitationally interacting disk galaxies and the formation of a galaxy in a fully cosmological simulation. It is found that the sticky particle model does a good job of reproducing many of the observed properties of interacting galaxies, including the properties of the ISM in the resulting elliptical galaxy