Simulating planet migration in globally evolving disks

A. Crida (1) , A. Morbidelli (2) , F. Masset (3, 4)

1 OCA - Observatoire de la Côte d'Azur
2 CASSIOPEE - Laboratoire de Cosmologie, Astrophysique Stellaire & Solaire, de Planétologie et de Mécanique des Fluides
3 DAPNIA - Département d'Astrophysique, de physique des Particules, de physique Nucléaire et de l'Instrumentation Associée
4 UNAM - Universidad Nacional Autónoma de México = National Autonomous University of Mexico

Context.Numerical simulations of planet-disk interactions are usually performed with hydro-codes that – because they consider only an annulus of the disk, over a 2D grid – cannot take into account the global evolution of the disk. However, the global evolution governs type II planetary migration, so that the accuracy of the planetary evolution can be questioned. Aims.To develop an algorithm that models the local planet-disk interactions together with the global viscous evolution of the disk. Methods.We surround the usual 2D grid with a 1D grid ranging over the real extension of the disk. The 1D and 2D grids are coupled at their common boundaries via ghost rings, paying particular attention to the fluxes at the interface, especially the flux of angular momentum carried by waves. The computation is done in the frame centered on the center of mass to ensure angular momentum conservation. Results.The global evolution of the disk and the local planet-disk interactions are both well described and the feedback of one on the other can be studied with this algorithm, for a negligible additional computing cost with respect to the usual algorithms.

Mots clés: methods: numerical; solar system: formation; accretion disks

Domaines: Planétologie et astrophysique de la terre [astro-ph.EP]

A. Crida, A. Morbidelli, F. Masset. Simulating planet migration in globally evolving disks. Astronomy and Astrophysics - A&A, 2007, 461 (3), pp.1173-1183. ⟨10.1051/0004-6361:20065870⟩. ⟨hal-00388110⟩