Application of the Gillespie algorithm to a granular intruder particle
arXiv:cond-mat/0608158 · doi:10.1088/0305-4470/39/35/001
Abstract
We show how the Gillespie algorithm, originally developed to describe coupled chemical reactions, can be used to perform numerical simulations of a granular intruder particle colliding with thermalized bath particles. The algorithm generates a sequence of collision ``events'' separated by variable time intervals. As input, it requires the position-dependent flux of bath particles at each point on the surface of the intruder particle. We validate the method by applying it to a one-dimensional system for which the exact solution of the homogeneous Boltzmann equation is known and investigate the case where the bath particle velocity distribution has algebraic tails. We also present an application to a granular needle in bath of point particles where we demonstrate the presence of correlations between the translational and rotational degrees of freedom of the intruder particle. The relationship between the Gillespie algorithm and the commonly used Direct Simulation Monte Carlo (DSMC) method is also discussed.
13 pages, 8 figures, to be published in J. Phys. A Math. Gen
References in corpus (9)
- Breakdown of Energy Equipartition in a 2D Binary Vibrated Granular Gas
- Detecting a small perturbation through its non-Gaussianity
- Computer simulation of uniformly heated granular fluids
- The Inelastic Maxwell Model
- Energy partition and segregation for an intruder in a vibrated granular system under gravity
- Dynamics of a tracer granular particle as a non-equilibrium Markov process
- Stationary states and energy cascades in inelastic gases
- Power-law velocity distributions in granular gases
- Thermalization of an anisotropic granular particle