Consistent SPH Simulations of Protostellar Collapse and Fragmentation

Abstract

We study the consistency and convergence of smoothed particle hydrodynamics (SPH), as a function of the interpolation parameters, namely the number of particles $N$, the number of neighbors $n$, and the smoothing length $h$, using simulations of the collapse and fragmentation of protostellar rotating cores. The calculations are made using a modified version of the GADGET-2 code that employs an improved scheme for the artificial viscosity and power-law dependences of $n$ and $h$ on $N$, as was recently proposed by Zhu et al., which comply with the combined limit $N\to \infty$, $h\to 0$, and $n\to \infty$ with $n/N\to 0$ for full SPH consistency, as the domain resolution is increased. We apply this realization to the "standard isothermal test case" in the variant calculated by Burkert & Bodenheimer and the Gaussian cloud model of Boss to investigate the response of the method to adaptive smoothing lengths in the presence of large density and pressure gradients. The degree of consistency is measured by tracking how well the estimates of the consistency integral relations reproduce their continuous counterparts. In particular, $C^0$ and $C^1$ particle consistency is demonstrated, meaning that the calculations are close to second-order accuracy. As long as $n$ is increased with $N$, mass resolution also improves as the minimum resolvable mass $M_{\min }\sim n^{-1}$. This aspect allows proper calculation of small-scale structures in the flow associated with the formation and instability of protostellar disks around the growing fragments, which are seen to develop a spiral structure and fragment into close binary/multiple systems as supported by recent observations.

34 pages, 15 figures, ApJ in press

References in corpus (20)

• The cosmological simulation code GADGET-2
• Cosmological SPH simulations: The entropy equation
• Improving convergence in smoothed particle hydrodynamics simulations without pairing instability
• Smoothed Particle Hydrodynamics and Magnetohydrodynamics
• Inviscid SPH
• Hydra: An Adaptive--Mesh Implementation of PPPM--SPH
• A Triple Protostar System Formed via Fragmentation of a Gravitationally Unstable Disk
• Smoothed Particle Hydrodynamics with particle splitting, applied to self-gravitating collapse
• Resolving mixing in Smoothed Particle Hydrodynamics
• SPHGal: Smoothed Particle Hydrodynamics with improved accuracy for Galaxy simulations
• Numerical Convergence in Smoothed Particle Hydrodynamics
• Galaxy mergers on a moving mesh: a comparison with smoothed-particle hydrodynamics
• Resolution requirements for simulating gravitational fragmentation using SPH
• Protostellar collapse: A comparison between SPH and AMR calculations
• VINE -- A numerical code for simulating astrophysical systems using particles II: Implementation and performance characteristics
• Protostellar collapse and fragmentation using an MHD GADGET
• Particle Methods in Astrophysical Fluid Dynamics
• SPH Simulations of Accretion Disks and Narrow Rings
• Gravitational Collapse and Fragmentation in Molecular Clouds with Adaptive Mesh Refinement
• On the thermal sensitivity of binary formation in collapsing molecular clouds

Cited by in corpus (6)

• Strongly misaligned triple system in SR 24 revealed by ALMA
• Relativistic dust accretion onto a scale--dependent polytropic black hole
• From large-scale to protostellar disk fragmentation into close binary stars
• Impetus: Consistent SPH calculations of 3D spherical Bondi accretion onto a black hole
• Radiative non-isothermal Bondi accretion onto a massive black hole
• Multi-scale physics of cryogenic liquid helium-4: Inverse coarse-graining properties of smoothed particle hydrodynamics