Collisional particle disks

Yoram Lithwick*, Eugene Chiang

*Corresponding author for this work

Research output: Contribution to journalArticle

18 Scopus citations

Abstract

We present a new, simple, fast algorithm to numerically evolve disks of inelastically colliding particles surrounding a central star. Our algorithm adds negligible computational cost to the fastest existing collisionless N-body codes and can be used to simulate, for the first time, the interaction of planets with disks over many viscous times. Although the algorithm is implemented in two dimensions - i.e., the motions of bodies need only be tracked in a plane - it captures the behavior of fully three-dimensional disks in which collisions maintain inclinations that are comparable to random eccentricities. The method simulates vertically optically thin disks of identical collisional, massless, inelastic, indestructible test particles. We subject the algorithm to a battery of tests for the case of an isolated narrow circular ring. Numerical simulations agree with analytic theory with regard to how particles' random velocities equilibrate, how the ring viscously spreads, and how energy dissipation, angular momentum transport, and material transport are connected. We derive and measure the critical value of the coefficient of restitution, above which viscous stirring dominates inelastic damping and the particles' velocity dispersion runs away.

Original languageEnglish (US)
Pages (from-to)524-533
Number of pages10
JournalAstrophysical Journal
Volume656
Issue number1 I
DOIs
StatePublished - Feb 10 2007

    Fingerprint

Keywords

  • Accretion, accretion disks
  • Planets: rings

ASJC Scopus subject areas

  • Astronomy and Astrophysics
  • Space and Planetary Science

Cite this