University of Oulu

A E Mondino-Llermanos, H Salo, An implementation of viscous pressure-force (‘soft-sphere’) model in REBOUND for local ring simulations, Monthly Notices of the Royal Astronomical Society, Volume 513, Issue 4, July 2022, Pages 4711–4728, https://doi.org/10.1093/mnras/stac1179

An implementation of viscous pressure-force (‘soft-sphere’) model in REBOUND for local ring simulations

Saved in:
Author: Mondino-Llermanos, A. E.1; Salo, H.1
Organizations: 1Space Physics and Astronomy Research Unit, University of Oulu, Pentti Kaiteran katu 1, FI-90014 Oulu, Finland
Format: article
Version: published version
Access: open
Online Access: PDF Full Text (PDF, 2.9 MB)
Persistent link: http://urn.fi/urn:nbn:fi-fe2023022328475
Language: English
Published: Oxford University Press, 2022
Publish Date: 2023-02-23
Description:

Abstract

Numerical simulations are an essential tool for studying dynamics of dense, self-gravitating planetary rings. Here, we present a new implementation of a soft-sphere collision model into the parallel gravitational N-body code REBOUND. Specifically, our module is designed for simulations done in the local shearing sheet approximation, and includes a simple viscoelastic model for the contact forces between colliding particles. The soft-sphere treatment is particularly important in regimes where collisions cannot be treated as independent binary impacts, as is done in the hard-sphere method already implemented in the code. The new module is fully parallelized with MPI and it can be used in conjunction with the available gravity octree method, also fully parallelized, to perform self-gravitating simulations. We check the validity of our soft-sphere model by successfully reproducing the main steady-state properties of non-gravitating and self-gravitating systems reported in earlier literature. Moreover, we show that it gives realistic results also in the high optical depth and strong self-gravity regimes, where the hard-sphere method fails due to particle overlaps. In addition, we test the CPU-scaling and efficiency of the MPI-parallelization. The algorithm proves to retain much of the original REBOUND speed without sacrificing the accuracy of the results. Our new tool can thus be applied with confidence to planetary ring studies, such as those aimed at understanding the onset and evolution of viscous overstability in dense self-gravitating rings.

see all

Series: Monthly notices of the Royal Astronomical Society
ISSN: 0035-8711
ISSN-E: 1365-8711
ISSN-L: 0035-8711
Volume: 513
Issue: 4
Pages: 4711 - 4728
DOI: 10.1093/mnras/stac1179
OADOI: https://oadoi.org/10.1093/mnras/stac1179
Type of Publication: A1 Journal article – refereed
Field of Science: 115 Astronomy and space science
Subjects:
Funding: We acknowledge funding by the Academy of Finland (HS) and Väisälä Foundation (AML).
Copyright information: © 2022 The Author(s). Published by Oxford University Press on behalf of Royal Astronomical Society. This is a pre-copyedited, author-produced version of an article accepted for publication in Monthly Notices of the Royal Astronomical Society following peer review. The version of record A E Mondino-Llermanos, H Salo, An implementation of viscous pressure-force (‘soft-sphere’) model in REBOUND for local ring simulations, Monthly Notices of the Royal Astronomical Society, Volume 513, Issue 4, July 2022, Pages 4711–4728, is available online at: https://doi.org/10.1093/mnras/stac1179.