RODEO: a new method for planet-disk interaction

¹ Leiden Observatory, Postbus 9513, 2300 RA Leiden, The Netherlands
² ASTRON, Postbus 2, 7990 AA Dwingeloo, The Netherlands e-mail: paardeko@strw.leidenuniv.nl;gmellema@astron.nl

Received: 5 July 2005
Accepted: 14 November 2005

Abstract

Aims.In this paper we describe a new method for studying the hydrodynamical problem of a planet embedded in a gaseous disk.

Methods.We use a finite volume method with an approximate Riemann solver (the Roe solver), together with a special way to integrate the source terms. This new source term integration scheme sheds new light on the Coriolis instability, and we show that our method does not suffer from this instability.

Results.The first results on flow structure and gap formation are presented, as well as accretion and migration rates. For M_p < 0.1  and M_p > 1.0 ( = Jupiter's mass) the accretion rates do not depend sensitively on numerical parameters, and we find that within the disk's lifetime a planet can grow to  . In between these two limits numerics play a major role, leading to differences of more than  for different numerical parameters. Migration rates are not affected by numerics at all as long as the mass inside the Roche lobe is not considered. We can reproduce the type I and type II migration for low-mass and high-mass planets, respectively, and the fastest moving planet of 0.1  has a migration time of only  yr.