[This page uses sraightforward Html and JavaScript. You can also get an XHtml version using MathML and SVG]

As if we didn’t need one, a new type of weapon is being designed which transfers the energy from a large moving mass to a much smaller and faster projectile.

The large mass is fired at the smaller projectile at rest in a tube which is closed at the far end. The projectile can bounce off the far end and come back to hit the input mass again, or the far end can be opened letting the projectile out. All collisions are elastic, there is no energy loss.

For each number of hits between the input mass and the projectile, there is a size for the large input mass such that all its energy is transferred to the projectile, which is finally fired at a higher speed.

The first few optimal input masses and resultant output velocities are as follows:

N m _{1}/m_{2}v _{2}/v_{1}1 1 1 2 3+2√2 1+√2 3 7+4√3 2+√3

N | Number of hits between the input mass and the projectile |
---|---|

m_{1} | Mass fired in |

m_{2} | Mass of the projectile |

v_{1} | Initial velocity of input mass |

v_{2} | Final velocity of projectile |

A good guess from a quick inspection might be that the
velocity ratio v_{2}/v_{1} goes up as
(N−1)+√N.
However it isn't that easy I’m afraid!

So the problem is get a simple formula for the N^{th}
velocity ratio v_{2}/v_{1}.

The N^{th} optimal mass ratio m_{1}/m_{2}
is this squared as the energy ∑½mv^{2} is conserved.

Firstly we can consider a bounce of the projectile off the end and then a hit with the input mass as a unit as they always come in pairs, we can treat that the iniial state as having its zero velocity reversed!

To cut through the straightforward dynamics here is how the
velocities of the input mass and projectile are changed when
the projectle bounces off the end and then collides with the
input mass (i.e. v_{2} → −v_{2} first):

v_{1}′ = ((m_{1}−m_{2})v_{1}− 2m_{2}v_{2}) / (m_{1}+m_{2})

v_{2}′ = (2m_{1}v_{1}+ (m_{1}−m_{2})v_{2}) / (m_{1}+m_{2})

These are got from the conservation of momentum ∑mv and the
conservation of energy ∑½mv^{2}.

Thus the velocities are changed by the matrix:

m _{1}−m_{2}m _{1}+m_{2}

−2m _{2}m _{1}+m_{2}

2m _{1}m _{1}+m_{2}

m _{1}−m_{2}m _{1}+m_{2}

We want the result of applying the matrix N times to [v_{1}, 0]
to result in [0, v_{2}] i.e. v_{1}′ is zero.

In fact this determinant of this matrix is 1 (It had to be because of conservation of energy but it's good to see that it actually is!). Also it can be written as:

d ^{-1}0 0 d

c -s s c

d 0 0 d ^{-1}

where:

d =

^{4}√

m _{1}m _{2}(Yes that’s right, the fourth root)

c =

m _{1}−m_{2}m _{1}+m_{2}

s =

−2√(m _{1}m_{2})m _{1}+m_{2}

So the matrix can be decomposed as:

D^{-1}RD

Multiplying the matrix N times gives D^{-1}R^{N}D.

R is a rotation matrix as c^{2}+s^{2} is 1. So N of
them are equivalent to rotating by N times the amount. To get the final
v_{1}′ equal to zero we need the cosine of the rotation to be
zero.

The rotation is by:

φ = cos ^{-1}(

m _{1}−m_{2}m _{1}+m_{2})

So we need cos(Nφ) = 0. To do this we need (higher values of φ won't occur in the problem):

φ = π / (2N)

(m

_{1}+m_{2}) cos φ = m_{1}−m_{2}m

_{1}(1−cos φ) = m_{2}(1+cos φ)m

_{1}/ m_{2}= (1+cos φ) / (1−cos φ)

Multiply top and bottom by (1+cos φ) and you’ll see that
the square root of that last expression,
which has to be the required result v_{2} / v_{1}
because of energy conservation, is:

v

_{2}/ v_{1}= (1+cos φ) / sin φwhere φ = π / (2N)

The limit as N→∞ is 4N/π.