Treating gravity in thin-disk simulations

¹ Institut für Astronomie & Astrophysik, Universität Tübingen, Auf der Morgenstelle 10, 72076 Tübingen, Germany
e-mail: Tobias_Mueller@twam.info
² Institute for Astronomy, ETH Zürich, Wolfgang-Pauli-Strasse 27, 8093 Zürich, Switzerland

Received: 23 December 2011
Accepted: 6 March 2012

Abstract

Context. In 2D-simulations of thin gaseous disks with embedded planets or self-gravity the gravitational potential needs to be smoothed to avoid singularities in the numerical evaluation of the gravitational potential or force. The softening prescription used in 2D needs to be adjusted properly to correctly resemble the realistic case of vertically extended 3D disks.

Aims. We analyze the embedded planet and the self-gravity case and provide a method to evaluate the required smoothing in 2D simulations of thin disks.

Methods. Starting from the averaged hydrodynamic equations and using a vertically isothermal disk model, we calculated the force to be used in 2D simulations. We compared our results to the often used Plummer form of the potential, which runs as  ∝ 1/(r² + ϵ²)^1/2. For that purpose we computed the required smoothing length ϵ as a function of distance r to the planet or to a disk element within a self-gravitating disk.

Results. We find that for longer distances ϵ is determined solely by the vertical disk thickness H. For the planet case we find that outside r ≈ H a value of ϵ = 0.7H describes the averaged force very well, while in the self-gravitating disk the value needs to be higher, ϵ = 1.2H. For shorter distances the smoothing needs to be reduced significantly. Comparing torque densities of 3D and 2D simulations we show that the modification to the vertical density stratification as induced by an embedded planet needs to be taken into account to obtain agreeing results.

Conclusions. It is very important to use the correct value of ϵ in 2D simulations to obtain a realistic outcome. In disk fragmentation simulations the choice of ϵ can determine whether a disk will fragment or not. Because a wrong smoothing length can change even the direction of migration, it is very important to include the effect of the planet on the local scale height in 2D planet-disk simulations. We provide an approximate and fast method for this purpose that agrees very well with full 3D simulations.