Issue |
A&A
Volume 675, July 2023
|
|
---|---|---|
Article Number | A96 | |
Number of page(s) | 12 | |
Section | Numerical methods and codes | |
DOI | https://doi.org/10.1051/0004-6361/202346178 | |
Published online | 06 July 2023 |
Self-gravity in thin-disc simulations of protoplanetary discs: The smoothing length rectified and generalised to bi-fluids
1
Aix-Marseille Univ., CNRS, CNES, LAM,
Marseille,
France
e-mail: pierre.barge@lam.fr
2
Aix-Marseille Univ., CNRS, Centrale Marseille, IRPHE,
Marseille,
France
3
Leibniz-Institut für Astrophysik Potsdam (AIP),
An der Sternwarte 16,
14482
Potsdam,
Germany
e-mail: srendon@aip.de
Received:
17
February
2023
Accepted:
5
April
2023
Context. To mimic the evolution of protoplanetary discs (PPDs), 2D simulations that incorporate self-gravity must introduce a softening prescription of the gravitational potential. When the disc is only composed of gas, the smoothing length is proportional to the scale height of the gas. On the other hand, when a dust component is included, the question arises as to whether the smoothing length approach can still be used to not only quantify the dust self-gravity, but its gravitational interaction with gas as well.
Aims. We identified grey areas in the standard smoothing length formalism for computing self-gravity in PPDs entirely made up of gas. Our goal is to revisit the smoothing length approach, which can then be generalised to two phases, whereby the dust component may be considered as a pressureless fluid.
Methods. We used analytical developments to approximate the vertically averaged self-gravity when the smoothing length is not assumed to be constant, but by taking a spatial function instead.
Results. We obtained an analytical expression for the space-varying smoothing length, which strongly improves the accuracy of the self-gravity computation. For the first time, this method is generalised to address bi-fluid interactions in a PPD: two additional smoothing lengths are proposed for featuring an isolated dusty disc and gas-dust self-gravity interactions. On computational grounds, we prescribe the use of tapering functions for the purpose of avoiding numerical divergences. We also checked that our method continues to be compatible with standard fast Fourier transform algorithms and evaluated computational costs.
Conclusions. Our space-varying smoothing length allows us to: (i) solve the contradictions inherent in the constant smoothing length hypothesis; (ii) fit the 3D vertically averaged self-gravity with a high level of accuracy; and (iii) render it applicable to a bi-fluid description of PPDs with the use of two additional smoothing lengths. Such results are crucial to enable realistic 2D numerical simulations that account for self-gravity and are essential to improving our understanding of planetesimal formation and type I migration.
Key words: protoplanetary disks / planet–disk interactions / gravitation / hydrodynamics / methods: analytical / methods: numerical
© The Authors 2023
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.