Systems of pinned billiard balls serve as simplified models of collisions, where all particles remain fixed in their positions while their (pseudo-)velocities evolve in accordance with the laws of conservation of energy and momentum. For some families of ball configurations, Athreya, Burdzy, and Duarte have established the maximum upper bound for the number of pseudo-collisions, thereby demonstrating that the number of collisions is finite. The result has been extended to all ball configurations. In this project, we do extensive simulations to study two specific configurations. First, we consider balls arranged in a half-space and assign a single ball an inward (pseudo-) velocity. Simulations suggest that in the long run, most of the energy is concentrated near the boundary. Second, when the balls are arranged on a flat torus, we find that in the stationary regime, the distributions of the velocity components are i.i.d. normal. Additionally, we find that the components of the velocities in the direction of impact between two touching balls are uncorrelated.