Toy Stars in two dimensions

Abstract

Toy Stars are gas masses where the compressibility is treated without approximations but gravity is replaced by a force which, for any pair of masses, is along their line of centres and proportional to their separation. They provide an invaluable resource for testing the suitability of numerical codes for astrophysical gas dynamics. In this paper we derive the equations for both small amplitude oscillations and non linear solutions for rotating and pulsating Toy Stars in two dimensions, and show that the solutions can be reduced to a small number of ordinary differential equations. We compare the accurate solutions of these equations with Smoothed Particle Hydrodynamics (SPH) simulations. The two dimensional Toy Star solutions are found to provide an excellent benchmark for SPH algorithms, highlighting many of the strengths and also some weaknesses of the method.

16 pages, 16 figures, MNRAS accepted

References in corpus (3)

• Cosmological SPH simulations: The entropy equation
• SPH compressible turbulence
• Smoothed Particle Magnetohydrodynamics II. Variational principles and variable smoothing length terms

Cited by in corpus (6)

• Phantom: A smoothed particle hydrodynamics and magnetohydrodynamics code for astrophysics
• MAGMA: a 3D, Lagrangian magnetohydrodynamics code for merger applications
• On dust entrainment in photoevaporative winds
• A study of cool core resiliency and entropy mixing in simulations of galaxy cluster mergers
• Outcomes of Grazing Impacts Between Sub-Neptunes in Kepler Multis
• Lattice Melting and Rotation in Perpetually Pulsating Equilibria