Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical convergence in self-gravitating disc simulations: initial conditions and edge effects
91
0
0.0
(
0
)
Added by
Sijme-Jan Paardekooper
Publication date
2011
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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.
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.
Generating pre-initial conditions (or particle loads) is the very first step to set up a cosmological N-body simulation. In this work, we revisit the numerical convergence of pre-initial conditions on dark matter halo properties using a set of simulations which only differs in initial particle loads, i.e. grid, glass, and the newly introduced capacity constrained Voronoi tessellation (CCVT). We find that the median halo properties agree fairly well (i.e. within a convergence level of a few per cent) among simulations running from different initial loads. We also notice that for some individual haloes cross-matched among different simulations, the relative difference of their properties sometimes can be several tens of per cent. By looking at the evolution history of these poorly converged haloes, we find that they are usually merging haloes or haloes have experienced recent merger events, and their merging processes in different simulations are out-of-sync, making the convergence of halo properties become poor temporarily. We show that, comparing to the simulation starting with an anisotropic grid load, the simulation with an isotropic CCVT load converges slightly better to the simulation with a glass load, which is also isotropic. Among simulations with different pre-initial conditions, haloes in higher density environments tend to have their properties converged slightly better. Our results confirm that CCVT loads behave as well as the widely used grid and glass loads at small scales, and for the first time we quantify the convergence of two independent isotropic particle loads (i.e. glass and CCVT) on halo properties.
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
