Long-term and large-scale viscous evolution of dense planetary rings

Planetary rings are common in the outer Solar System but their origin and long-term evolution is still a matter of debate. It is well known that viscous spreading is a major evolutionary process for rings, as it globally redistributes the disk’s mass and angular momentum, and can lead to the disk’s loosing mass by infall onto the planet or through the Roche limit. However, describing this process is highly dependent on the model used for the viscosity. In this paper we investigate the global and long-term viscous evolution of a circumplanetary disk. We have developed a simple 1D numerical code, but we use a physically realistic viscosity model derived from N-body simulations (Daisaka et al., 2001), and dependent on the disk’s local properties (surface mass density, particle size, distance to the planet). Particularly, we include the effects of gravitational instabilities (wakes) that importantly enhance the disk’s viscosity. This method allows to study the global evolution of the disk over the age of the Solar System.Common estimates of the disk’s spreading time-scales with constant viscosity significantly underestimate the rings’ lifetime. We show that, with a realistic viscosity model, an initially narrow ring undergoes two successive evolutionary stages: (1) a transient rapid spreading when the disk is self-gravitating, with the formation of a density peak inward and an outer region marginally gravitationally stable, and with an emptying time-scale proportional to (where M₀ is the disk’s initial mass), (2) an asymptotic regime where the spreading rate continuously slows down as larger parts of the disk become non-self-gravitating due to the decrease of the surface density, until the disk becomes completely non-self-gravitating. At this point its evolution dramatically slows down, with an emptying time-scale proportional to 1/M₀, which significantly increases the disk’s lifetime compared to the case with constant viscosity. We show also that the disk’s width scales like t^1/4 with the realistic viscosity model, while it scales like t^1/2 in the case of constant viscosity, resulting in much larger evolutionary time-scales in our model. We find however that the present shape of Saturn’s rings looks like a 100 million-years old disk in our simulations. Concerning Jupiter’s, Uranus’ and Neptune’s rings that are faint today, it is not likely that they were much more massive in the past and lost most of their mass due to viscous spreading alone.     

A numerical model of cohesion in planetary rings

We present a numerical method that incorporates particle sticking in simulations using the N-body code pkdgrav to study motions in a local rotating frame, such as a patch of a planetary ring. Particles stick to form non-deformable but breakable aggregates that obey the (Eulerian) equations of rigid-body motion. Applications include local simulations of planetary ring dynamics and planet formation, which typically feature hundreds of thousands or more colliding bodies. Bonding and breaking thresholds are tunable parameters that can approximately mimic, for example, van der Waals forces or interlocking of surface frost layers. The bonding and breaking model does not incorporate a rigorous treatment of internal fracture; rather the method serves as motivation for first-order investigation of how semi-rigid bonding affects the evolution of particle assemblies in high-density environments.We apply the method to Saturn’s A ring, for which laboratory experiments suggest that interpenetration of thin, frost-coated surface layers may lead to weak cohesive bonding. These experiments show that frost-coated icy bodies can bond at the low impact speeds characteristic of the rings. Our investigation is further motivated by recent simulations that suggest a very low coefficient of restitution is needed to explain the amplitude of the azimuthal brightness asymmetry in Saturn’s A ring, and the hypothesis that fine structure in Saturn’s B ring may in part be caused by large-scale cohesion.This work presents the full implementation of our model in pkdgrav, as well as results from initial tests with a limited set of parameters explored. We find a combination of parameters that yields aggregate size distribution and maximum radius values in agreement with Voyager data for ring particles in Saturn’s outer A ring. We also find that the bonding and breaking parameters define two strength regimes in which fragmentation is dominated either by collisions or other stresses, such as tides. We conclude our study with a discussion of future applications of and refinements to our model.     

Adhesion and collisional release of particles in dense planetary rings

We propose a simple theoretical model for aggregative and fragmentative collisions in Saturn’s dense rings. In this model the ring matter consists of a bimodal size distribution: large (meter sized) boulders and a population of smaller particles (tens of centimeters down to dust). The small particles can adhesively stick to the boulders and can be released as debris in binary collisions of their carriers. To quantify the adhesion force we use the JKR theory (Johnson, K., Kendall, K., Roberts, A. [1971]. Proc. R. Soc. Lond. A 324, 301–313). The rates of release and adsorption of particles are calculated, depending on material parameters, sizes, and plausible velocity dispersions of carriers and debris particles. In steady state we obtain an expression for the amount of free debris relative to the fraction still attached to the carriers. In terms of this conceptually simple model a paucity of subcentimeter particles in Saturn’s rings (French, R.G., Nicholson, P.D. [2000]. Icarus 145, 502–523; Marouf, E. et al. [2008]. Abstracts for “Saturn after Cassini–Huygens” Symposium, Imperial College London, UK, July 28 to August 1, p. 113) can be understood as a consequence of the increasing strength of adhesion (relative to inertial forces) for decreasing particle size. In this case particles smaller than a certain critical radius remain tightly attached to the surfaces of larger boulders, even when the boulders collide at their typical speed. Furthermore, we find that already a mildly increased velocity dispersion of the carrier-particles may significantly enhance the fraction of free debris particles, in this way increasing the optical depth of the system.     

N-body simulations of cohesion in dense planetary rings: A study of cohesion parameters

We present results from a large suite of simulations of Saturn’s dense A and B rings using a new model of particle sticking in local simulations (Perrine, R.P., Richardson, D.C., Scheeres, D.J. [2011]. Icarus 212, 719–735). In this model, colliding particles can be incorporated into or help fragment rigid aggregations on the basis of certain user-specified parameters that can represent van der Waals forces or interlocking surface frost layers.Our investigation is motivated by laboratory results that show that interpenetration of surface layers can allow impacting frost-covered ice spheres to stick together. In these experiments, cohesion only occurs below specific impact speeds, which happen to be characteristic of impact speeds in Saturn’s rings. Our goal is to determine if weak bonding is consistent with ring observations, to constrain cohesion parameters in light of existing ring observations, to make predictions about particle populations throughout the rings, and to discover other diagnostics that may constrain bonding parameters.We considered the effects of five parameters on the equilibrium characteristics of our ring simulations: speed-based merge and fragmentation limits, bond strength, ring surface density, and patch orbital distance (i.e., the A or B ring), some with both monodisperse and polydisperse comparison cases. In total, we present data from 95 simulations.We find that weak cohesion is consistent with observations of the A and B rings (e.g., French, R.G., Nicholson, P.D. [2000]. Icarus 145, 502–523), and we present a range of simulation parameters that reproduce the observed size distribution and maximum particle size. It turns out that the parameters that match observations differ between the A and B rings, and we discuss the potential implications of this result. We also comment on other observable consequences of cohesion for the rings, such as optical depth and scale height effects, and discuss whether very large objects (e.g., “propeller” source objects) are grown bottom-up from cohesion of smaller ring particles.     

The viscous overstability, nonlinear wavetrains, and finescale structure in dense planetary rings

This paper addresses the fine-scale axisymmetric structure exhibited in Saturn's A and B-rings. We aim to explain both the periodic microstructure on 150-220 m, revealed by the Cassini UVIS and RSS instruments, and the irregular variations in brightness on 1-10 km, reported by the Cassini ISS. We propose that the former structures correspond to the peaks and troughs of the nonlinear wavetrains that form naturally in a viscously overstable disk. The latter variations on longer scales may correspond to modulations and defects in the wavetrains' amplitudes and wavelength. We explore these ideas using a simple hydrodynamical model which captures the correct qualitative behaviour of a disk of inelastically colliding particles, while also permitting us to make progress with analytic and semi-analytic techniques. Specifically, we calculate a family of travelling nonlinear density waves and determine their stability properties. Detailed numerical simulations that confirm our basic results will appear in a following paper.     

Gravitational scattering in planetary rings

The gravitational influence of a large body (moonlet; satellite) on the radial structure of planetary rings has been calculated numerically. A drastical change of the surface mass density is obtained even after a single scattering process of the ring-particles on a moonlet (satellite). The final surface density shows a significant radial structure, which has been used to estimate radius and mass of satellites embedded in rings of low optical depth (E-ring, Cassini-division, C-ring of Saturn).     

A planetary dynamo model

Using the Lagrangian approach, the author considered the temporal evolution of an ensemble of interacting magnetohydrodynamic cyclones, obeying equations of the Langevin type, in a rotating medium. The problem is topical for fast-rotating convective objects: cores of planets and a number of stars, where the Rossby numbers are much less than unity and the geostrophic balance of forces is observed. In this work, results of simulation are given both for the two-dimensional case, when axes of cyclones can rotate only in the vertical plane, and for the three-dimensional case when the axes are rotating by two angles. It is shown that a change in the heat flux on the shell boundary impacts the frequency of reversals of the mean dipole magnetic field, which agrees with results of simulation in three-dimensional models of a planetary dynamo. Applications of the model for the giant planets are considered, and an explanation of certain episodes of the geomagnetic field in the past is offered.     

A new formulation of the viscosity in planetary rings

We present a new formulation of the viscosity in planetary rings, where particles interact through their gravitational forces and direct collisions. In the previous studies on the viscosity in self-gravitating rings, the viscosity consists of three components, which are defined separately in different ways. The complex definitions make it difficult to evaluate the viscosity in N-body simulation of rings. In our new formulation, the viscosity is expressed in terms of changes in orbital elements of particles due to particle interactions. This makes the expression of the viscosity simple. The new formulation gives a simple way to evaluate the viscosity in N-body simulation. We find that for practical evaluation of the viscosity of planetary rings, only energy dissipation at direct inelastic collisions is needed.For tenuous particle disks (i.e., optically thin disks), we further derive a formula of the viscosity. The formula requires only a numerical coefficient that can be obtained from three-body calculation. Since planetesimal disks are also tenuous, the viscosity in planetesimal disks can be also obtained from this formula. In a subsequent paper, we will evaluate this coefficient through three-body calculation and obtain the viscosity for a wide range of parameters such as the restitution coefficient and the radial location in rings.     

Inelastic collisions in narrow planetary rings

It is shown that a particle ring with energy dissipation has an extremum in its energy when all the particles are in the same circular orbit. This extremum is a relative maximum in radial directions, indicating possible radial expansion; but it is a relative minimum in the particle velocity components, indicating a tendency for the velocity distribution to collapse. An N-body model of ring evolution incorporating two-body dynamics, oblateness perturbations, inelastic collisions, and phase averaging is described. By local analysis of impact statistics, it is shown that the velocity distribution of the ring will collapse if the coefficient of restitution ∈ ≲ 0.7. The collapse of the velocity distribution stabilizes a ring of point particles against radial dispersion. Furthermore, for ∈ ≲ 0.25, the semimajor axis distribution tends to collapse toward its local mean value, leading to radial collapse. A survey of ring evolution is presented for different values of coefficient of restitution and initial velocity dispersion. As has been predicted uy Goldreich and Tremaine, rings with equilibrium velocity distributions are unstable and expand, while rings where the velocity distribution collapses are shown to undergo massive, pervasive fragmentation into a myriad of ringlets. It is proposed that such fragmented rings are stable in their own right, and that the observational test to discriminate between the two cases is simply the presence of a smooth, featureless surface or the presence of intricate radial structure in the ring.     

Nonlinear density waves in planetary rings

We discuss the steady-state structure of the nonlinear density waves generated in a planetary ring at the Lindblad resonances of a satellite. We show that strong density waves lead to an enhancement of the background surface density in the wave zone.     

The disruption of planetary satellites and the creation of planetary rings

The Voyager spacecraft discovered that small moons orbit within all four observed ring systems coincident with the discovery of narrow and dusty rings around Jupiter, Saturn, Uranus and Neptune. These moons can provide the source for new rings if they are catastrophically disrupted by a comet or large meteoroid impact. This hypothesis for ring origins provides a natural mechanism for the ongoing creation of planetary rings. While it relieves somewhat the problem of explaining the continued existence of rings with apparently short evolutionary lifetimes, it raises the problem of explaining the continued existence of small moons, and the coexistence of moons and rings at comparable locations within the Roche zones of the giant planets. This problem has been studied in some detail recently, and the present work is a review of our current understanding of the processes in satellite disruption that pertain to the creation of planetary rings and the collisional cascade of circumplanetary bodies. Significant progress has been made. Narrow rings are produced by disruption of small moons in numerical simulations, and a self-consistent model of the collisional cascade can explain present-day moon populations. Absolute timescales and initial moon populations remain uncertain due to our poor knowledge of the impactor population and uncertainties in the strength of planetary satellites. More pressing are the qualitative issues that remain to be resolved including the nature of reaccretion of the debris and the origin of Saturn's rings.     

Dense planetary rings and the viscous overstability

This paper examines the onset of the viscous overstability in dense particulate rings. First, we formulate a dense gas kinetic theory that is applicable to the saturnian system. Our model is essentially that of Araki and Tremaine [Araki, S., Tremaine, S., 1986. Icarus 65, 83-109], which we show can be both simplified and generalised. Second, we put this model to work computing the equilibrium properties of dense planetary rings, which we subsequently compare with the results of N-body simulations, namely those of Salo [Salo, H., 1991. Icarus 90, 254-270]. Finally, we present the linear stability analyses of these equilibrium states, and derive criteria for the onset of viscous overstability in the self-gravitating and non-self-gravitating cases. These are framed in terms of particle size, orbital frequency, optical depth, and the parameters of the collision law. Our results compare favourably with the simulations of Salo et al. [Salo, H., Schmidt, J., Spahn, F., 2001. Icarus 153, 295-315]. The accuracy and practicality of the continuum model we develop encourages its general use in future investigations of nonlinear phenomena.     

Simulation of evolutionary mechanism of planetary rings

We devised a numerical model of planetary ring in which the inelastic collision and the gravitation between ring particles are considered. We adopt Hill's equations and the differential algorithm of circular mesh for computation of the particle orbits. The evolutional processes are presented for different coefficient of restitution and dynamical optical depth . The results show that the semi-major axis and eccentricity of the ring particles are changed with and . We compute the average energies transferred and loss in inelastic collisions for various values of the parameters. The dynamical equilibrium properties are discussed in the different cases.     

A Roche model for uniformly rotating rings

A Roche model for describing uniformly rotating rings is presented, and the results are compared with the numerical solutions to the full problem for polytropic rings. In the thin ring limit, the surfaces of constant pressure including the surface of the ring itself are given in analytical terms, even in the mass-shedding case.     

Size distribution of particles in planetary rings

Harris (Icarus24, 190–192) has suggested that the maximum size of particles in a planetary ring is controlled by collisional fragmentation rather than by tidal stress. While this conclusion is probably true, estimated radius limits must be revised upward from Harris' values of a few kilometers by at least an order of magnitude. Accretion of particles within Roche's limit is also possible. These considerations affect theories concerning the evolution of Saturn's rings, of the Moon, and of possible former satellites of Mercury and Venus. In the case of Saturn's rings, comparison of various theoretical scenarios with available observational evidence suggests that the rings formed from the breakup of larger particles rather than from original condensation as small particles. This process implies a distribution of particle sizes in Saturn's rings possibly ranging up to ～100 km but with most cross-section in cm-scale particles.     

Description and behavior of streamlines in planetary rings

The dynamical behavior of planetary rings is often described in terms of streamline deformations and motions. We argue that the standard formula giving the shape of the streamlines involves geometric elements rather than orbital elements, when the oblateness of the planet is taken into account. We relate the geometric elements to the orbital elements and derive the differential equations which govern the variations with time of the geometric elements.     

N-body simulations of viscous instability of planetary rings

We study viscous instability of planetary rings in terms of N-body simulations. We show that for rings composed of fairly elastic particles (e.g. as in Hatzes et al. [Hatzes, A., Bridges, F.G., Lin, D.N.C., 1988. Collisional properties of ice spheres at low impact velocities. Mon. Not. R. Astron. Soc. 231, 1091-1115]) the instability may lead to the spontaneous formation of dense ringlets in a background of lower density. In most parts of Saturn’s rings the particle collisions are probably much more dissipative, as suggested by the presence of self-gravity wakes, and classic viscous instability should be suppressed. However, our results demonstrate that the mechanism of viscous instability itself is valid. The dynamical effects of size-dependent elasticity in a system with a size distribution have never been studied before. We show that this may in principle lead to a size-selective viscous instability, small particles concentrating on ringlets against the more uniform background of large particles.     

Electrostatic forces and shadow effects in planetary dust rings

The vertical dust distribution of dust clouds around planets, resulting from electrostatic forces, is calculated as a function of dust and plasma parameters. Photoelectric charging is included and differences between clouds on the illuminated side and in the shadow zone are examined. We compute ring structures for conditions which may apply in the spoke-forming regions and study at what dust and plasma conditions the shadow has a significant effect on the vertical dust cloud structure.