Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPH

Abstract

We implemented sink particles in the adaptive mesh refinement (AMR) hydrodynamics code FLASH. Sink particles are created in regions of local gravitational collapse, and their trajectories and accretion can be followed over many dynamical times. We perform a series of tests including the time integration of circular and elliptical orbits, the collapse of a Bonnor-Ebert sphere and a rotating, fragmenting cloud core. We compare the collapse of a highly unstable singular isothermal sphere to the theory by Shu (1977), and show that the sink particle accretion rate is in excellent agreement with the theoretical prediction. To model eccentric orbits and close encounters of sink particles accurately, we show that a very small timestep is often required, for which we implemented subcycling of the N-body system. We emphasize that a sole density threshold for sink particle creation is insufficient in supersonic flows, if the density threshold is below the opacity limit. In that case, the density can exceed the threshold in strong shocks that do not necessarily lead to local collapse. Additional checks for bound state, gravitational potential minimum, Jeans instability and converging flows are absolutely necessary for a meaningful creation of sink particles. We apply our new sink particle module for FLASH to the formation of a stellar cluster, and compare to a smoothed particle hydrodynamics (SPH) code with sink particles. Our comparison shows encouraging agreement of gas properties, indicated by column density distributions and radial profiles, and of sink particle formation times and positions. We find excellent agreement in the number of sink particles formed, and in their accretion and mass distributions.

30 pages, 17 figures, ApJ accepted, simulation movies available at http://www.ita.uni-heidelberg.de/~chfeder/videos.shtml?lang=en

References in corpus (15)

• Theory of Star Formation
• Analytical theory for the initial mass function: CO clumps and prestellar cores
• The Density Probability Distribution in Compressible Isothermal Turbulence: Solenoidal versus Compressive Forcing
• Radiation-Hydrodynamic Simulations of Collapse and Fragmentation in Massive Protostellar Cores
• Star formation through gravitational collapse and competitive accretion
• The Importance of Radiative Feedback for the Stellar Initial Mass Function
• Numerical simulations of compressively driven interstellar turbulence: I. Isothermal gas
• The Fractal Density Structure in Supersonic Isothermal Turbulence: Solenoidal versus Compressive Energy Injection
• The effect of magnetic fields on star cluster formation
• A test suite for quantitative comparison of hydrodynamics codes in astrophysics
• A Direct Multigrid Poisson Solver for Oct-Tree Adaptive Meshes
• The star formation efficiency and its relation to variations in the initial mass function
• Protostellar collapse: A comparison between SPH and AMR calculations
• Simulating hydromagnetic processes in star formation: introducing ambipolar diffusion into an adaptive mesh refinement code
• Numerical comparison of Riemann solvers for astrophysical hydrodynamics

Cited by in corpus (5)

• A new Jeans resolution criterion for (M)HD simulations of self-gravitating gas: Application to magnetic field amplification by gravity-driven turbulence
• Limiting Accretion onto Massive Stars by Fragmentation-Induced Starvation
• Understanding Spatial and Spectral Morphologies of Ultracompact H II Regions
• Metallicity and the Universality of the IMF
• The impact of X-rays on molecular cloud fragmentation and the IMF