Journal Article

On the convergence of the critical cooling time-scale for the fragmentation of self-gravitating discs

Farzana Meru and Matthew R. Bate

in Monthly Notices of the Royal Astronomical Society

Published on behalf of The Royal Astronomical Society

Volume 427, issue 3, pages 2022-2046
Published in print December 2012 | ISSN: 0035-8711
Published online December 2012 | e-ISSN: 1365-2966 | DOI: http://dx.doi.org/10.1111/j.1365-2966.2012.22035.x
On the convergence of the critical cooling time-scale for the fragmentation of self-gravitating discs

Show Summary Details

Preview

Abstract

We carry out simulations of gravitationally unstable discs using a smoothed particle hydrodynamics (SPH) code and a grid-based hydrodynamics code, fargo, to understand the previous non-convergent results reported by Meru & Bate. We obtain evidence that convergence with increasing resolution occurs with both SPH and fargo and in both cases we find that the critical cooling time-scale is larger than previously thought. We show that SPH has a first-order convergence rate, while fargo converges with a second-order rate. We show that the convergence of the critical cooling time-scale for fragmentation depends largely on the numerical viscosity employed in both SPH and fargo. With SPH, particle velocity dispersion may also play a role. We show that reducing the dissipation from the numerical viscosity leads to larger values of the critical cooling time at a given resolution. For SPH, we find that the effect of the dissipation due to the numerical viscosity is somewhat larger than had previously been appreciated. In particular, we show that using a quadratic term in the SPH artificial viscosity (βSPH) that is too low appears to lead to excess dissipation in gravitationally unstable discs, which may affect any results that sensitively depend on the thermodynamics, such as disc fragmentation. We show that the two codes converge to values of the critical cooling time-scale, βcrit > 20 (for a ratio of specific heats of γ = 5/3), and perhaps even as large as βcrit ≈ 30. These are approximately three to five times larger than has been found by most previous studies. This is equivalent to a maximum gravitational stress that a disc can withstand without fragmenting of αGI, crit ≈ 0.013 − 0.02, which is much smaller than the values typically used in the literature. It is therefore easier for self-gravitating discs to fragment than has been concluded from most past studies.

Keywords: accretion, accretion discs; gravitation; hydrodynamics; instabilities; planets and satellites: formation; protoplanetary discs

Journal Article.  16909 words.  Illustrated.

Subjects: Astronomy and Astrophysics

Full text: subscription required

How to subscribe Recommend to my Librarian

Users without a subscription are not able to see the full content. Please, subscribe or login to access all content.