Evolution of planetesimals. I - Dynamics: Relaxation in a thin disk. II - Numerical simulations

The study examines the effects of density inhomogeneity and differential rotation as well as inelastic collisions on the dynamical evolution of planetesimals. Consideration is given to a three-step analysis: the dynamical evolution of the planetesimals, collisions and mass accumulation, and interaction with gas. It is shown that the velocity dispersion of a cold system of planetesimals increases rapidly due to elastic gravitational scattering. When the dispersion in the epicycle amplitude becomes comparable to the planetesimals' Roche radius, energy is transferred from the systematic Keplerian shear to the dispersive motion. With a numerical N-body scheme, gravitational scattering and physical collisions among a system of planetesimals is simulated. It is shown that dynamical equilibrium is attained with a velocity dispersion comparable to the surface escape velocity of those planetesimals which contribute most of the system mass.