Buoyancy response of a disk to an embedded planet: a cross-code comparison at high resolution
Abstract
In radiatively inefficient, laminar protoplanetary disks, embedded planets can excite a buoyancy response as gas gets deflected vertically near the planet. This results in vertical oscillations that drive a vortensity growth in the planet's corotating region, speeding up inward migration in the type-I regime. We present a comparison between PLUTO/IDEFIX and FARGO3D using 3D, inviscid, adiabatic numerical simulations of planet-disk interaction that feature the buoyancy response of the disk, and show that PLUTO/IDEFIX struggle to resolve higher-order modes of the buoyancy-related oscillations, weakening vortensity growth and the associated torque. We interpret this as a drawback of total-energy-conserving, finite-volume schemes. Our results indicate that a very high resolution or high-order scheme is required in shock-capturing codes in order to adequately capture this effect.
Cite
@article{arxiv.2308.14896,
title = {Buoyancy response of a disk to an embedded planet: a cross-code comparison at high resolution},
author = {Alexandros Ziampras and Sijme-Jan Paardekooper and Richard P. Nelson},
journal= {arXiv preprint arXiv:2308.14896},
year = {2023}
}
Comments
13 pages, 17 figures, 1 table; accepted by MNRAS