Aims. Our goal is to find a converged cooling limit for fragmentation in self-gravitating disks. This is especially interesting for the formation of planets, brown dwarfs or stars and the growth of black holes. While investigating the limit, we want to give a clear criterion for the state of convergence. Methods. We run two-dimensional shearingsheet simulations with the hydrodynamic package Fosite at high resolutions. Thereby resolution and limiters are altered. Subsequently, we investigate the spectra of important physical quantities at the length scales where fragmentation occurs. In order to avoid prompt fragmentation at high resolutions we start these simulations with a fully developed gravitoturbulent state obtained at a lower resolution. Results. We show nearly converged results for fragmentation with a critical cooling timescale $t_{\mathrm{crit}} \sim 10\,\Omega^{-1}$ . We can backtrace this claim by investigating the spectra of relevant physical variables at length scales around and below the pressure scale height. We argue that well behaved results cannot be expected if counteracting quantities are varying too much on these critical length scales, either by change of resolution or numerical method. A comparison of fragmentation behaviour with the related spectra reveals that simulations behave similar, if the spectra are converged to the length scales where self-gravity leads to instabilities. Observable deviations in the results obtained with different numerical setup are confined to scales below these critical length scales.

Context. It is still an open issue whether a self-gravitating accretion disk fragments. There are many different physical and numerical explanations for fragmentation, but simulations often show a non-convergent behavior for ever better resolution. Aims. We investigate the influence of different numerical limiters in Godunov type schemes on the fragmentation boundary in self- gravitating disks. Methods. We compare the linear and non-linear outcome in two-dimensional shearingsheet simulations using the VANLEER and the SUPERBEE limiter. Results. We show that choosing inappropriate limiting functions to handle shock-capturing in Godunov type schemes can lead to an overestimation of the surface density in regions with shallow density gradients. The effect amplifies itself on timescales comparable to the dynamical timescale even at high resolutions. This is exactly the environment, where clumps are expected to form. The effect is present without, but scaled up by, self-gravity and also does not depend on cooling. Moreover it can be backtracked to a well known effect called oversteepening. If the effect is also observed in the linear case, the fragmentation limit is shifted to larger values of the critical cooling timescale.