اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدNumerical simulations of the type III migration:I. Disc model and convergence tests
93
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We investigate the fast (type III) migration regime of high-mass protoplanets orbiting in protoplanetary disks. This type of migration is dominated by corotational torques. We study the details of flow structure in the planets vicinity, the dependence of migration rate on the adopted disc model, and the numerical convergence of models (independence of certain numerical parameters such as gravitational softening). We use two-dimensional hydrodynamical simulations with adaptive mesh refinement,based on the FLASH code with improved time-stepping scheme. We perform global disk simulations with sufficient resolution close to the planet, which is allowed to freely move throughout the grid. We employ a new type of equation of state in which the gas temperature depends on both the distance to the star and planet, and a simplified correction for self-gravity of the circumplanetary gas. We find that the migration rate in the type III migration regime depends strongly on the gas dynamics inside the Hill sphere (Roche lobe of the planet) which, in turn, is sensitive to the aspect ratio of the circumplanetary disc. Furthermore, corrections due to the gas self-gravity are necessary to reduce numerical artifacts that act against rapid planet migration. Reliable numerical studies of Type III migration thus require consideration of both the thermal andthe self-gravity corrections, as well as a sufficient spatial resolution and the calculation of disk-planet attraction both inside and outside the Hill sphere. With this proviso, we find Type III migration to be a robust mode of migration, astrophysically promising because of a speed much faster than in the previously studied modes of migration.
قيم البحث
اقرأ أيضاً
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 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 dis
sipation 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 present a numerical study of rapid, so called type III migration for Jupiter-sized planets embedded in a protoplanetary disc. We limit ourselves to the case of outward migration, and study in detail its evolution and physics, concentrating on the
structure of the co-rotation and circumplanetary regions, and processes for stopping migration. We also consider the dependence of the migration behaviour on several key parameters. We perform this study using global, two-dimensional hydrodynamical simulations with adaptive mesh refinement. We find that the outward directed type III migration can be started if the initial conditions support $Z > 1$, that corresponds to initial value $M_rmn{Delta} ga 1.5$. Unlike the inward directed migration, in the outward migration the migration rate increases due to the growing of the volume of the co-orbital region. We find the migration to be strongly dependent on the rate of the mass accumulation in the circumplanetary disc, leading to two possible regimes of migration, fast and slow. The structure of the co-orbital region and the stopping mechanism differ between these two regimes.
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 planet
s. 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.
We present a numerical study of rapid, so called type III migration for Jupitersized planets embedded in a protoplanetary disc. We limit ourselves to the case of inward migration, and study in detail its evolution and physics, concentrating on the st
ructure of the corotation and circumplanetary regions, and processes for stopping migration. We also consider the dependence of the migration behaviour on several key parameters. We perform this study using the results of global, two-dimensional hydrodynamical simulations with adaptive mesh refinement. The initial conditions are chosen to satisfy the condition for rapid inward migration. We find that type III migration can be divided into two regimes, fast and slow. The structure of the coorbital region, mass accumulation rate, and migration behaviour differ between these two regimes. All our simulations show a transition from the fast to the slow regime, ending type III migration well before reaching the star. The stopping radius is found to be larger for more massive planets and less massive discs. A sharp density drop is also found to be an efficient stopping mechanism. In the fast migration limit the migration rate and induced eccentricity are lower for less massive discs, but almost do not depend on planet mass. Eccentricity is damped on the migration time scale.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
