Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical convergence in self-gravitating shearing sheet simulations and the stochastic nature of disc fragmentation
218
0
0.0
(
0
)
Added by
Sijme-Jan Paardekooper
Publication date
2012
fields
Physics
and research's language is
English
Authors
Sijme-Jan Paardekooper
Ask ChatGPT about the research
No Arabic abstract
We study numerical convergence in local two-dimensional hydrodynamical simulations of self-gravitating accretion discs with a simple cooling law. It is well-known that there exists a steady gravito-turbulent state, in which cooling is balanced by dissipation of weak shocks, with a net outward transport of angular momentum. Previous results indicated that if cooling is too fast (typical time scale 3/Omega, where Omega is the local angular velocity), this steady state can not be maintained and the disc will fragment into gravitationally bound clumps. We show that, in the two-dimensional local approximation, this result is in fact not converged with respect to numerical resolution and longer time integration. Irrespective of the cooling time scale, gravito-turbulence consists of density waves as well as transient clumps. These clumps will contract because of the imposed cooling, and collapse into bound objects if they can survive for long enough. Since heating by shocks is very local, the destruction of clumps is a stochastic process. High numerical resolution and long integration times are needed to capture this behaviour. We have observed fragmentation for cooling times up to 20/Omega, almost a factor 7 higher than in previous simulations. Fully three-dimensional simulations with a more realistic cooling prescription are necessary to determine the effects of the use of the two-dimensional approximation and a simple cooling law.
rate research
Read More
We study the numerical convergence of hydrodynamical simulations of self-gravitating accretion discs, in which a simple cooling law is balanced by shock heating. It is well-known that there exists a critical cooling time scale for which shock heating can no longer compensate for the energy losses, at which point the disc fragments. The numerical convergence of previous results of this critical cooling time scale was questioned recently using Smoothed Particle Hydrodynamics (SPH). We employ a two-dimensional grid-based code to study this problem, and find that for smooth initial conditions, fragmentation is possible for slower cooling as the resolution is increased, in agreement with recent SPH results. We show that this non-convergence is at least partly due to the creation of a special location in the disc, the boundary between the turbulent and the laminar region, when cooling towards a gravito-turbulent state. Converged results appear to be obtained in setups where no such sharp edges appear, and we then find a critical cooling time scale of ~ 4 $Omega^{-1}$, where $Omega$ is the local angular velocity.
We study the stability of gaps opened by a giant planet in a self-gravitating protoplanetary disc. We find a linear instability associated with both the self-gravity of the disc and local vortensity maxima which coincide with gap edges. For our models, these edge modes develop and extend to twice the orbital radius of a Saturn mass planet in discs with disc-to-star mass ratio >0.06, corresponding to a Toomre Q < 1.5 at the outer disc boundary. Unlike the local vortex-forming instabilities associated with gap edges in weakly or non-self-gravitating low viscosity discs, the edge modes are global and exist only in sufficiently massive discs, but for the typical viscosity values adopted for protoplanetary discs. Analytic modelling and linear calculations show edge modes may be interpreted as a localised disturbance associated with a gap edge inducing activity in the extended disc, through the launching of density waves excited at Lindblad resonances. Nonlinear hydrodynamic simulations are performed to investigate the evolution of edge modes in disc-planet systems. The form and growth rates of unstable modes are consistent with linear theory. Their dependence on viscosity and gravitational softening is also explored. We also performed a first study of the effect of edge modes on planetary migration. We found that if edge modes develop, then the average disc-on-planet torque becomes more positive with increasing disc mass. In simulations where the planet was allowed to migrate, although a fast type III migration could be seen that was similar to that seen in non-self-gravitating discs, we found that it was possible for the planet to interact gravitationally with the spiral arms associated with an edge mode and that this could result in the planet being scattered outwards. Thus orbital migration is likely to be complex and non monotonic in massive discs of the type we consider.
It is likely that most protostellar systems undergo a brief phase where the protostellar disc is self-gravitating. If these discs are prone to fragmentation, then they are able to rapidly form objects that are initially of several Jupiter masses and larger. The fate of these disc fragments (and the fate of planetary bodies formed afterwards via core accretion) depends sensitively not only on the fragments interaction with the disc, but with its neighbouring fragments. We return to and revise our population synthesis model of self-gravitating disc fragmentation and tidal downsizing. Amongst other improvements, the model now directly incorporates fragment-fragment interactions while the disc is still present. We find that fragment-fragment scattering dominates the orbital evolution, even when we enforce rapid migration and inefficient gap formation. Compared to our previous model, we see a small increase in the number of terrestrial-type objects being formed, although their survival under tidal evolution is at best unclear. We also see evidence for disrupted fragments with evolved grain populations - this is circumstantial evidence for the formation of planetesimal belts, a phenomenon not seen in runs where fragment-fragment interactions are ignored. In spite of intense dynamical evolution, our population is dominated by massive giant planets and brown dwarfs at large semimajor axis, which direct imaging surveys should, but only rarely, detect. Finally, disc fragmentation is shown to be an efficient manufacturer of free floating planetary mass objects, and the typical multiplicity of systems formed via gravitational instability will be low.
The gravitational interaction between a protoplanetary disc and planetary sized bodies that form within it leads to the exchange of angular momentum, resulting in migration of the planets and possible gap formation in the disc for more massive planets. In this article, we review the basic theory of disc-planet interactions, and discuss the results of recent numerical simulations of planets embedded in protoplanetary discs. We consider the migration of low mass planets and recent developments in our understanding of so-called type I migration when a fuller treatment of the disc thermodynamics is included. We discuss the runaway migration of intermediate mass planets (so-called type III migration), and the migration of giant planets (type II migration) and the associated gap formation in the disc. The availability of high performance computing facilities has enabled global simulations of magnetised, turbulent discs to be computed, and we discuss recent results for both low and high mass planets embedded in such discs.
Recently it has been suggested that the fragmentation boundary in Smoothed Particle Hydrodynamic (SPH) and FARGO simulations of self-gravitating accretion discs with beta-cooling do not converge as resolution is increased. Furthermore, this recent work suggests that by carefully optimising the artificial viscosity parameters in these codes it can be shown that fragmentation may occur for much longer cooling times than earlier work suggests. If correct, this result is intriguing as it suggests that gas giant planets could form, via direct gravitational collapse, reasonably close to their parent stars. This result is, however, slightly surprising and there have been a number of recent studies suggesting that the result is likely an indication of a numerical problem with the simulations. One suggestion, in particular, is that the SPH results are influenced by the manner in which the cooling is implemented. We extend this work here and show that if the cooling is implemented in a manner that removes a known numerical artefact in the shock regions, the fragmentation boundary converges to a value consistent with earlier work and that fragmentation is unlikely for the long cooling times suggested by this recent work. We also investigate the optimisation of the artificial viscosity parameters and show that the values that appear optimal are likely introducing numerical problems in both the SPH and FARGO simulations. We therefore conclude that earlier predictions for the cooling times required for fragmentation are likely correct and that, as suggested by this earlier work, fragmentation cannot occur in the inner parts (r < 50 au) of typical protostellar discs.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
