Mean Free Path in Sinai's Billiards

Sinai's Billiards is a term referring to the study of the chaotic dynamical properties of hard elastic balls. (See the visualization page for details). Below we present a few computational results for the mean free path of a ray in a periodic cubic lattice of reflective spheres.

The following graph shows an experimental result, and a suitable fit to the data that was guessed at.



The raw data is here. It seems to be accurate to maybe three or four decimal places. The data seems to be fit with about 5% accuracy with the equation

mfp = exp (1-r8) / r4
where r is the sphere radius and the lattice is two units on a side. The green line in the graph (labelled 'fit') is a graph of this formula. I have no idea as to whether this is a well-known result or not; most technical papers on the web speak at a far more complex level, and introductory texts are not available on the web, or in my library. I don't know whether this formula has the right asymptotic behavior as r->0 or r >~ 1. (Because this is a three dimensional lattice, the ball radius can be as large as sqrt(3) before the space is filled, and as large as sqrt(2) before the space becomes disconnected. I don't know if there is a quantitatively different behavior for r > 1.)



Copyright (c) 2002 Linas Vepstas

Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.1; with no Invariant Sections, with no Front-Cover Texts, and with no Back-Cover Texts. A copy of the license is included at the URL http://www.linas.org/fdl.html, the web page titled "GNU Free Documentation License".

To contact Linas, see his Home Page.