Tree Code Simulations of Planetary Rings

A tree code method that incorporates a local shearing disc model and fourth-order integration algorithm is applied to the problem of planetary rings, with particular emphasis on the dynamics of Saturn's B ring. The new code, described in detail elsewhere, allows for particle self-gravity, a distribution of particle sizes, and surface friction (particle spin). Important changes made to the original code, to ensure an accurate treatment of collisions under severe high-density conditions, are described in detail. Comparison with work by Wisdom & Tremaine for the case of equal-size particles and mean self-gravity shows excellent agreement. Similar analysis is performed for the new regimes of selfgravity and particle size distributions, and it is shown that the condition for viscous instability is still not satisfied for these models. Particle spins lie generally within a rotational energy equipartition envelope, and are retrograde on average in the rotating (orbital) frame but exhibit a la..