Lattice Points on Circles and Discrete Velocity Models for the Boltzmann   Equation
  The construction of discrete velocity models or numerical methods for the Boltzmann equation, may lead to the necessity of computing the collision operator as a sum over lattice points. The collision operator involves an integral over a sphere, which corresponds to the conservation of energy and momentum. In dimension two there are difficulties even in proving the convergence of such an approximation since many circles contain very few lattice points, and some circles contain many badly distributed lattice points. However, by showing that lattice points on most circles are equidistributed we find that the collision operator can indeed be approximated as a sum over lattice points in the two-dimensional case. For higher dimensions, this result has already been obtained by A. Bobylev et. al. (SIAM J. Numerical Analysis 34 no 5 p. 1865-1883 (1997)) 
L. Fainsilber, P. Kurlberg and B. Wennberg