A cosmological hydrodynamic code based on the piecewise parabolic method

Abstract

We present a hydrodynamical code for cosmological simulations which uses the piecewise parabolic method (PPM) to follow the dynamics of the gas component and an N-body particle-mesh algorithm for the evolution of the collisionless component. The gravitational interaction between the two components is regulated by the Poisson equation which is solved by a standard fast Fourier transform (FFT) procedure. In order to simulate cosmological flows we have introduced several modifications to the original PPM scheme which we describe in detail. Various tests of the code are presented including adiabatic expansion, single and multiple pancake formation and three-dimensional cosmological simulations with initial conditions based on the cold dark matter scenario.

hydrodynamics, methods: numerical, large-scale structure of Universe