Mixing of Primordial Elements in a Supernova-Driven ISM

Abstract

In order to understand the time scales for complete mixing in a supernova-driven ISM we used the model of Avillez (2000) to follow the mixing and dispersal of inhomogeneities with different length scales. We included tracer fields to follow the inhomogeneities in kpc-scale simulations with different supernova rates. We find that an initially inhomogeneous ISM still shows local variations and a clumpy distribution on time scales of several hundred Myr. This clumpy structure does not disapear over a time of a few tens of Myr, even if the rate of SNe is increased to ten times the Galactic rate. These local variations occur because, even in the presence of numerical diffusion, gas does not mix quickly between hot and cold regions. The lower limit that we are able to place on the mixing time at the kpc scale is already longer than the time for chemistry to reach local equilibrium in the gas. The simulations also show that the time scale to erase inhomogeneities in the ISM is nearly independent of their length scale over the range 25-500 pc, contrary to what is expected from classical mixing length theory.