Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] adbe_v1_0.gz(49 Kbytes)
Manuscript Title: Reduced event-list on an array for many-body simulation.
Authors: K. Shida, S. Yamada
Program title: N_3dim
Catalogue identifier: ADBE_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 86(1995)289
Programming language: C.
Computer: SPARCstation IPX.
Operating system: SUN-OS4.1.3, HPUX9.05.
RAM: 110K words
Word size: 32
Keywords: N-body simulation, Event-list, Hard-core particle, Solid state physics, Other, Many-body system, Simulation algorithm, Data reduction, Collisional system.
Classification: 7.7.

Nature of problem:
Inelastic colliding N-body simulations without mutual attractive force, such as some kind of molecular dynamics [1] or dynamics of planetary ringlets [2].

Solution method:
Alder and Wainwright's method [3] with our event-reduction mechanism.

Restrictions:
No external field. This method is essentially applicable for the cases of central gravitational field, f = k r**-2, harmonic oscillation, f = k r, and other problems in which particle trajectories are predictable without numerical integration.

Unusual features:
This program runs even if the number of potential events exceeds the capacity of the event-list. Typical optimul capacity is only N.

Running time:
Several ten seconds for the test data (1600 particles).

References:
[1] K. Shida and T. Kawai, Cluster formation by inelastically colliding particles in one-dimensional space, Physica A 162 (1989) 145.
[2] K. Shida, K. Suzuki and T. Kawai, Formation of sharp ringlets by inelastic collisions, J. Phys. Soc. Japan 60 (1991) 3953.
[3] B.J. Alder and T.E. Wainwright, Studies in molecular dynamics. I. General method, J. Chem. Phys. 11 (1959) 459.