Large-scale N-body simulations of the viscous overstability in Saturn's rings

Abstract

We present results from large-scale particle simulations of the viscous overstability in Saturn's rings. The overstability generates a variety of structure on scales covering a few hundred metres to several kilometres, including axisymmetric wavetrains and larger-scale modulations. Such patterns have been observed in Saturn's rings by the Cassini spacecraft. Our simulations model the collisional evolution of particles in a co-rotating patch of the disk. These are the largest N-body simulations of the viscous overstability yet performed. The radial box size is five orders of magnitude larger than a typical particle radius, and so describes a 20-50 km radial portion of the rings. Its evolution is tracked for more than 10,000 orbits. In agreement with hydrodynamics, our N-body simulations reveal that the viscous overstability exhibits a rich set of dynamics characterised by nonlinear travelling waves with wavelengths of a few hundred meters. In addition, wave defects, such as sources and shocks, punctuate this bed of waves and break them up into large-scale divisions of radial width ~5 km. We find that the wavelength of the travelling waves is positively correlated with the mean optical depth. In order to assess the role of the numerical boundary conditions and also background ring structure, we include simulations of broad spreading rings and simulations with a gradient in the background surface density. Overall, our numerical results and approach provide a tool with which to interpret Cassini occultation observations of microstructure in Saturn's rings. We present an example of such a synthetic occultation observation and discuss what features to expect. We also make the entire source code freely available.

15 pages, 12 figures, accepted for publication by MNRAS

References in corpus (5)

• REBOUND: An open-source multi-purpose N-body code for collisional dynamics
• Viscous overstability and eccentricity evolution in three-dimensional gaseous discs
• The linear stability of dilute particulate rings
• Dense planetary rings and the viscous overstability
• The ballistic transport instability in Saturn's rings I: formalism and linear theory

Cited by in corpus (13)

• Mapping spiral waves and other radial features in Saturn's rings
• Numerical Simulations of Saturn's B Ring: Granular Friction as a Mediator between Self-Gravity and Viscous Overstability
• Exploring overstabilities in Saturn's A ring using two stellar occultations
• Probing the inner boundaries of Saturn's A ring with the Iapetus -1:0 nodal bending wave
• Viscous Driving of Global Oscillations in Accretion Discs Around Black Holes
• First observational evidence of a relation between globular clusters' internal rotation and stellar masses
• Density Waves and the Viscous Overstability in Saturn's Rings
• Viscous Overstability in Saturn's Rings: Influence of Collective Self-gravity
• Filament formation due to diffusive instabilities in dusty protoplanetary disks
• On the tidal environment of an outwardly migrated F-ring
• Planetary Ring Dynamics -- The Streamline Formalism -- 2. Theory of Narrow Rings and Sharp Edges
• Implementation and Performance of Barnes-Hut N-body algorithm on Extreme-scale Heterogeneous Many-core Architectures
• Thermal hysteresis and front propagation in dense planetary rings