Imagine the beginning of a game of pool, you have 16 balls, 15 of them in a triangle < and 1 of them being the cue ball off to the left of that triangle. Imagine that the rack (the 15 balls in a triangle) has every ball equally spaced apart and all balls touching all other appropriate balls. All balls are perfectly round. Now, imagine that the cue ball was hit along a friction free surface on the center axis for this triangle O< and hits the far left ball of the rack dead center on this axis. How would the rack react? I would imagine this would be an extension of newtons cradle and only the 5 balls on the far end would move at all. But in what way would they move? Thanks

2This seems like a nice question, but what does it have to do with the vector spaces. Did you try $3+1$ balls to begin with? – Alex Degtyarev Jan 31 '14 at 10:12

2For many particles in contact, conservation of energy and momentum do not determine the outcome, so I don't think your problem has a unique solution, unless only the far left ball moves. Physically, there are fractional forces between the balls and between the balls and surface, that depend on how the arrangement was set up. – user25199 Jan 31 '14 at 10:21

1I would appreciate it if someone else could give the correct tags for this problem. I'm not a complete genius at physics and I could not create the tags I thought were relevant, nor find any tags that I believed to be more relevant. Vector spaces was a guess as it had vector in it. I am writing a pool simulation program, and am looking for a qualitative answer to base my program around, rather than a specific quantitative result. – Phedg1 Jan 31 '14 at 11:53

2Hmm, I guess you are right. The problem seems underdetermined, with my solution being only one possibility. – Yoav Kallus Jan 31 '14 at 21:38

4@tacorama I think physically the roughness of the balls (and hence friction) is probably more important than quantum effects. The fact that the outcome depends on the chosen potential means that it is important to get the physics right to develop an accurate simulation tool. – user25199 Feb 6 '14 at 17:05
This question was crossposted on Math Stack Exchange. Here is a copy of my answer for it there.
This is it. The perfectly centered billiards break. Behold.
Setup
This break was computed in Mathematica using a numerical differential equations model. Here are a few details of the model:
 All balls are assumed to be perfectly elastic and almost perfectly rigid.
 Each ball has a mass of 1 unit and a radius of 1 unit.
 The cue ball has a initial speed of 10 units/sec.
 The force between two balls is given by the formula $$ F \;=\; \begin{cases}0 & \text{if }d \geq 2, \\ 10^{11}(2d)^{3/2} & \text{if }d<2,\end{cases} $$ where $d$ is the distance between the centers of the balls. Note that the balls overlap if and only if $d < 2$. The power of $3/2$ was suggested by Yoav Kallus in the comments, because it follows Hertz's theory of nonadhesive elastic contact.
The initial speed of the cue ball is immaterial  slowing down the cue ball is the same as slowing down time. The force constant $10^{11}$ has no real effect as long as it's large enough, although it does change the speed at which the initial collision takes place.
The Collision
For this model, the entire collision takes place in the first 0.2 milliseconds, and none of the balls overlap by more than 0.025% of their radius during the collision. (These figures are model dependent  real billiard balls may collide faster or slower than this.)
The following animation shows the forces between the balls during the collision, with the force proportional to the area of each yellow circle. Note that the balls themselves hardly move at all during the collision, although they do accelerate quite a bit.
The Trajectories
The following picture shows the trajectories of the billiard balls after the collision.
After the collision, some of the balls are travelling considerably faster than others. The following table shows the magnitude and direction of the velocity of each ball, where $0^\circ$ indicates straight up.
$$ \begin{array}{ccccccccccc} \hline \text{ball} & \text{cue} & 1 & 2,3 & 4,6 & 5 & 7,10 & 8,9 & 11,15 & 12,14 & 13 \\ \hline \text{angle} & 0^\circ & 0^\circ & 40.1^\circ & 43.9^\circ & 0^\circ & 82.1^\circ & 161.8^\circ & 150^\circ & 178.2^\circ & 180^\circ \\ \hline \text{speed} & 1.79 & 1.20 & 1.57 & 1.42 & 0.12 & 1.31 & 0.25 & 5.60 & 2.57 & 2.63 \\ \hline \end{array} $$
For comparison, remember that the initial speed of the cue ball was 10 units/sec. Thus, balls 11 and 15 (the back corner balls) shoot out at more than half the speed of the original cue ball, whereas ball 5 slowly rolls upwards at less than 2% of the speed of the original cue ball.
By the way, if you add up the sum of the squares of the speeds of the balls, you get 100, since kinetic energy is conserved.
Linear and Quadratic Responses
The results of this model are dependent on the power of $3/2$ in the force law  other force laws give other breaks. For example, we could try making the force a linear function of the overlap distance (in analogy with springs and Hooke's law), or we could try making the force proportional to the square of the overlap distance. The results are noticeably different
Stiff Response
Glenn the Udderboat points out that "stiff" balls might be best approximated by a force response involving a higher power of the distance (although this isn't the usual definition of "stiffness"). Unfortunately, the calculation time in Mathematica becomes longer when the power is increased, presumably because it needs to use a smaller time step to be sufficiently accurate.
Here is a simulation involving a reasonably "stiff" force law $$ F \;=\; \begin{cases}0 & \text{if }d \geq 2, \\ 10^{54}(2d)^{10} & \text{if }d<2.\end{cases} $$
As you can see, the result is very similar to my initial answer on Math Stack Exchange. This seems like good evidence that the behavior discussed in my initial answer is indeed the limiting behavior in the case where this notion of "stiffness" goes to infinity.
As you might expect, most of the energy in this case is transferred very quickly at the beginning of the collision. Almost all of the energy has moves to the back corner balls in the first 0.02 milliseconds. Here is an animation of the forces:
After that, the corner balls and the cue ball shoot out, and the remaining balls continue to collide gently for the next millisecond or so.
While the simplicity of this behavior is appealing, I would guess that "real" billard balls do not have such a force response. Of the models listed here, the intial Hertzbased model is probably the most accurate. Qualitatively, it certainly seems the closest to an "actual" break.
Note: I have now posted the Mathematica code on my web page.

11

2"Of the models listed here, I think the initial, linear model is probably the most accurate." Hertz's theory of nonadhesive elastic contact (en.wikipedia.org/wiki/Contact_mechanics) gives $F\sim(2d)^{3/2}$. – Yoav Kallus Feb 1 '14 at 21:36

1@YoavKallus Thanks for the suggestion! I have incorporated this into the model. – Jim Belk Feb 1 '14 at 22:16

2@user2376055 From a programming perspective, everything difficult here was done by Mathematica, which has a builtin numerical differential equations solver. I don't know what algorithm(s) it uses to solve a problem like this, but it's potentially very complicated, and I'm not sure I would have been capable of solving this problem in, say, Java (or at least, not without finding a good numerical mathematics library). If you like, I'd be happy to share my Mathematica code with you. – Jim Belk Feb 1 '14 at 23:11

1Using the trick suggested in this answer, you can probably figure out which numerical method Mathematica is using. – Federico Poloni Feb 2 '14 at 17:59
This is more of a comment inspired by Jim Belk's answer than an answer to the question in itself. It is also more physics than mathematics. However, I hope I can help readers see that this system is related to some very interesting science.
I wanted to say a little bit about the dynamics of the energy transfer in the billiard break. Very naïvely, you might think of a triangular lattice of interacting particles (the billiard rack) as acting like a discrete model for a crystalline solid. Solids support sound waves (just knock on your table and listen), and so again naïvely we might expect that the billiard break occurs via some strong compressional sound waves propagating through the system, akin to those vibrations induced when knocking on wood. However, this naïve intuition is wrong!
A system of billiards all just barely touching, like Newton's cradle, is an example of a sonic vacuum (see e.g. this expository physics article).
The term was coined by Vitali Nesterenko to describe mechanical systems where the speed of sound vanishes and the tiniest perturbation must lead to supersonic shocks  behavior very far from that of ordinary sound. The usual formula for the speed of sound in say a chain of identical masses joined by identical springs is $v_s=\sqrt{k/m}a$ where $k$ is the spring constant, $m$ is the mass, and $a$ is the length of the springs. In a system interacting via a nonlinear force law like Hertz's, the effective spring constant is $df/d\delta_{\delta=0}=0$ (where $\delta$ is the overlap between particles) which actually vanishes! Thus the speed of sound is zero. This is a very special situation which occurs because the particles are all just barely touching.
The consequence of all this is that energy cannot propagate via ordinary sound, and what appears instead is a shock wave, satisfying a nonlinear PDE. In Jim Belk's pictures, this is what we see as the "shock front" colored yellow traveling from the point of impact through the billiards.
Nesterenko wrote down a continuum limit for systems of Hertzian particles (the following is specialized to particles arranged on a triangular lattice):
$\frac{\partial^2\xi}{\partial t^2}=c^2\frac{\partial^2}{\partial x^2}\left[\xi^{3/2}+\frac{2R^2}{5}\xi^{1/4}\frac{\partial^2}{\partial x^2}(\xi^{5/4})\right]$
(If one wants to study shock waves in systems of particles interacting with other force laws, then some of the exponents here change). Also, none of the above discussion really depends on the "perfectly centered" initial condition. I think the paper I cite below may give you some idea of what can happen as we vary the initial angle.
It might be interesting to compare the solutions of this equation with suitable boundary conditions to simulations along the line of Jim Belk's for a very large pool ball rack (large enough at least so that one can follow the shock wave for an appreciable length of time).
The derivation of the above equation is apparently in Nesterenko's book "Dynamics of Heterogeneous Materials", though I learned it from this very relevant paper coauthored by some friends of mine: Transmission and reflection of strongly nonlinear solitary waves at granular interfaces, by A. M. Tichler, L. R. Gomez, N. Upadhyaya, X. Campman, V. F. Nesterenko, V. Vitelli. In this paper they study the propagation of shocks through interfaces between packings of billiard balls of two different masses. I think you'll agree it's quite close to the billiard break.
Here's figure 2 of Tichler et al's paper.
protected by user9072 Feb 6 '14 at 16:13
Thank you for your interest in this question.
Because it has attracted lowquality or spam answers that had to be removed, posting an answer now requires 10 reputation on this site (the association bonus does not count).
Would you like to answer one of these unanswered questions instead?