Toy stars in one dimension

J.J. Monaghan, D.J. Price

School of Mathematical Sciences, Monash University, Clayton 3800, Australia
Institute of Astronomy, Madingley Rd, Cambridge, CB3 0HA, UK

MNRAS 350, 1449-1456, 2004


This paper describes a simple model of a star where compressibility is retained but the gravitational force is replaced by a force between pairs which is directed along their line of centres and proportional to their separation. This force was analysed by Newton. It may be considered the simplest many body force because the system reduces to a set of independent particles moving in a harmonic oscillator potential. We call such models Toy Stars. In this paper we simplify further and focus on the one dimensional problem. We show that a non linear solution exists where the velocity is a linear function of the coordinate and the density a quadratic function of the coordinates. We generalise this result to include a magnetic field. Our results provide a very useful bench mark for algorithms designed to simulate gravitating gas and we show that Smoothed Particle Hydrodynamics (SPH) simulations of Toy Stars with or without magnetic fields give results in good agreement with theory.

[ pdf offprint* ]

« Back to publications
« Daniel Price's home page

*The pdf version is an electronic version of an article published in Monthly Notices of the Royal Astronomical Society: complete citation information for the final version of the paper, as published in the print edition of Monthly Notices of the Royal Astronomical Society, is available on the Blackwell Science Synergy online delivery service, accessible via the journal's Website at: