{"title":"Horseshoes and spiral waves: capturing the 3D flow induced by a low-mass planet analytically","authors":"Joshua J. Brown, Gordon I. Ogilvie","doi":"arxiv-2409.02687","DOIUrl":null,"url":null,"abstract":"The key difficulty faced by 2D models for planet-disc interaction is in\nappropriately accounting for the impact of the disc's vertical structure on the\ndynamics. 3D effects are often mimicked via softening of the planet's\npotential; however, the planet-induced flow and torques often depend strongly\non the choice of softening length. We show that for a linear adiabatic flow\nperturbing a vertically isothermal disc, there is a particular vertical average\nof the 3D equations of motion which exactly reproduces 2D fluid equations for\narbitrary adiabatic index. There is a strong connection here with the\nLubow-Pringle 2D mode of the disc. Correspondingly, we find a simple, general\nprescription for the consistent treatment of planetary potentials embedded\nwithin '2D' discs. The flow induced by a low-mass planet involves large-scale\nexcited spiral density waves which transport angular momentum radially away\nfrom the planet, and 'horseshoe streamlines' within the co-orbital region. We\nderive simple linear equations governing the flow which locally capture both\neffects faithfully simultaneously. We present an accurate co-orbital flow\nsolution allowing for inexpensive future study of corotation torques, and\npredict the vertical structure of the co-orbital flow and horseshoe region\nwidth for different values of adiabatic index, as well as the vertical\ndependence of the initial shock location. We find strong agreement with the\nflow computed in 3D numerical simulations, and with 3D one-sided Lindblad\ntorque estimates, which are a factor of 2 to 3 times lower than values from\nprevious 2D simulations.","PeriodicalId":501209,"journal":{"name":"arXiv - PHYS - Earth and Planetary Astrophysics","volume":"11 1","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2024-09-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - PHYS - Earth and Planetary Astrophysics","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2409.02687","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
The key difficulty faced by 2D models for planet-disc interaction is in
appropriately accounting for the impact of the disc's vertical structure on the
dynamics. 3D effects are often mimicked via softening of the planet's
potential; however, the planet-induced flow and torques often depend strongly
on the choice of softening length. We show that for a linear adiabatic flow
perturbing a vertically isothermal disc, there is a particular vertical average
of the 3D equations of motion which exactly reproduces 2D fluid equations for
arbitrary adiabatic index. There is a strong connection here with the
Lubow-Pringle 2D mode of the disc. Correspondingly, we find a simple, general
prescription for the consistent treatment of planetary potentials embedded
within '2D' discs. The flow induced by a low-mass planet involves large-scale
excited spiral density waves which transport angular momentum radially away
from the planet, and 'horseshoe streamlines' within the co-orbital region. We
derive simple linear equations governing the flow which locally capture both
effects faithfully simultaneously. We present an accurate co-orbital flow
solution allowing for inexpensive future study of corotation torques, and
predict the vertical structure of the co-orbital flow and horseshoe region
width for different values of adiabatic index, as well as the vertical
dependence of the initial shock location. We find strong agreement with the
flow computed in 3D numerical simulations, and with 3D one-sided Lindblad
torque estimates, which are a factor of 2 to 3 times lower than values from
previous 2D simulations.