Convergence of SPH simulations of self-gravitating accretion discs: Sensitivity to the implementation of radiative cooling

Details Convergence of SPH simulations of self-gravitating accretion discs: Sensitivity to the implementation of radiative cooling Convergence of SPH simulations of self-gravitating accretion discs: Sensitivity to the implementation of radiative cooling

2011-11-14

arxiv.org/abs/1111.3147v1

Astronomy and Astrophysics

Astrophysics

Earth and Planetary Astrophysics

9 pages, 11 figures, MNRAS in press

Scientific paper

Recent simulations of self-gravitating accretion discs, carried out using a three-dimensional Smoothed Particle Hydrodynamics (SPH) code by Meru and Bate, have been interpreted as implying that three-dimensional global discs fragment much more easily than would be expected from a two-dimensional local model. Subsequently, global and local two-dimensional models have been shown to display similar fragmentation properties, leaving it unclear whether the three-dimensional results reflect a physical effect or a numerical problem associated with the treatment of cooling or artificial viscosity in SPH. Here, we study how fragmentation of self-gravitating disc flows in SPH depends upon the implementation of cooling. We run disc simulations that compare a simple cooling scheme, in which each particle loses energy based upon its internal energy per unit mass, with a method in which the cooling is derived from a smoothed internal energy density field. For the simple per particle cooling scheme, we find a significant increase in the minimum cooling time scale for fragmentation with increasing resolution, matching previous results. Switching to smoothed cooling, however, results in lower critical cooling time scales, and tentative evidence for convergence at the highest spatial resolution tested. We conclude that precision studies of fragmentation using SPH require careful consideration of how cooling (and, probably, artificial viscosity) is implemented, and that the apparent non-convergence of the fragmentation boundary seen in prior simulations is likely a numerical effect. In real discs, where cooling is physically smoothed by radiative transfer effects, the fragmentation boundary is probably displaced from the two-dimensional value by a factor that is only of the order of unity.

Affiliated with

Also associated with

No associations

Rating

Convergence of SPH simulations of self-gravitating accretion discs: Sensitivity to the implementation of radiative cooling does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Convergence of SPH simulations of self-gravitating accretion discs: Sensitivity to the implementation of radiative cooling, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Convergence of SPH simulations of self-gravitating accretion discs: Sensitivity to the implementation of radiative cooling will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-217788