Free access article

A&A 499, 633-641 (2009)
DOI: 10.1051/0004-6361/20078822

Multiphase ISM simulations: comparing NIRVANA and ZEUS

R. A. Piontek, O. Gressel, and U. Ziegler

Astrophysikalisches Institut Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany
    e-mail: rpiontek@aip.de

Received 9 October 2007 / Accepted 12 February 2009

Abstract
Aims. Our aim is to compare the ZEUS and NIRVANA MHD codes in their application to multi-phase ISM simulations. A contrast in density of about 100 between the warm and cold ISM phases is typical, and discontinuities such as these can be difficult to handle numerically. It is plausible that the details of implementing a particular numerical method could affect the results of the simulation, particularly with regard to the mass and volume fractions of the warm, unstable, and cold gas. This work takes a step towards quantifying to what extent this may occur.
Methods. Using the NIRVANA code, we simulate a multi-phase ISM with turbulence driven by the magnetorotational instability. We make comparisons to our previously published ZEUS models and assess the impact of the numerical method on quantities such as the mass and volume fractions in the ISM phases, turbulent velocities, and probability distribution functions of density, temperature, pressure, and magnetic field strength. We also compare models with which turbulence is driven by an artificial forcing function.
Results. For the models that include turbulence driven by the magnetorotational instability, we find differences in the mass fractions that are typically on the order of 10%. This comparison, however, is complicated by the result that the saturation level of the MRI is different between the two codes. This suggests that, at the same saturation level, we would expect to find larger differences in the mass fractions, and this is indeed what we have found for our artificially forced models. For these runs we find differences in the cold mass fraction of approximately 30%, which is significant.

Key words: magnetohydrodynamics (MHD) -- ISM: kinematics and dynamics -- ISM: magnetic fields -- methods: numerical -- galaxies: ISM -- turbulence

© ESO 2009

The OpenURL standard is a protocol for transmission of metadata describing the resource that you wish to access. An OpenURL link contains article metadata and directs it to the OpenURL server of your choice. The OpenURL server can provide access to the resource and also offer complementary services (specific search engine, export of references...). The OpenURL link can be generated by different means.

• If your librarian has set up your subscription with an OpenURL resolver, OpenURL links appear automatically on the abstract pages.
• You can define your own OpenURL resolver with your EDPS Account. In this case your choice will be given priority over that of your library.
• You can use an add-on for your browser (Firefox or I.E.) to display OpenURL links on a page (see http://www.openly.com/openurlref/). You should disable this module if you wish to use the OpenURL server that you or your library have defined.