Issue |
A&A
Volume 678, October 2023
|
|
---|---|---|
Article Number | A134 | |
Number of page(s) | 11 | |
Section | Numerical methods and codes | |
DOI | https://doi.org/10.1051/0004-6361/202245601 | |
Published online | 17 October 2023 |
Measuring the numerical viscosity in simulations of protoplanetary disks in Cartesian grids
The viscously spreading ring revisited
1
Institut für Astronomie und Astrophysik, Universität Tübingen,
Auf der Morgenstelle 10,
72076
Tübingen, Germany
2
Astronomy Unit, Department of Physics and Astronomy, Queen Mary University of London,
London
E1 4NS, UK
e-mail: a.ziampras@qmul.ac.uk; r.p.nelson@qmul.ac.uk
3
Leibniz-Institut für Astrophysik Potsdam (AIP),
An der Sternwarte 16,
14482
Potsdam, Germany
e-mail: jjoseph@aip.de
Received:
1
December
2022
Accepted:
24
July
2023
Context. Hydrodynamical simulations solve the governing equations on a discrete grid of space and time. This discretization causes numerical diffusion similar to a physical viscous diffusion, the magnitude of which is often unknown or poorly constrained. With the current trend of simulating accretion disks with no or very low prescribed physical viscosity, it has become essential to understand and quantify this inherent numerical diffusion in the form of a numerical viscosity.
Aims. We study the behavior of the viscous spreading ring and the spiral instability that develops in it. We aim to use this setup to quantify the numerical viscosity in Cartesian grids and study its properties.
Methods. We simulated the viscous spreading ring and the related instability on a two-dimensional polar grid using PLUTO as well as FARGO, ensuring the convergence of our results with a resolution study. We then repeated our models on a Cartesian grid and measured the numerical viscosity by comparing results to the known analytical solution using PLUTO and Athena++.
Results. We find that the numerical viscosity in a Cartesian grid scales with resolution as approximately vnum ∝ Δx2 and is equivalent to an effective α ~ 10−4 for a common numerical setup. We also showed that the spiral instability manifests as a single leading spiral throughout the whole domain on polar grids. This is contrary to previous results and indicates that sufficient resolution is necessary in order to correctly resolve the instability.
Conclusions. Our results are relevant in the context of models where the origin should be included in the computational domain, or when polar grids cannot be used. Examples of such cases include models of disk accretion onto a central binary and, inherently, Cartesian codes.
Key words: hydrodynamics / accretion, accretion disks / protoplanetary disks / 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.