Open Access
 Issue Emergent Scientist Volume 1, 2017 International Physicists' Tournament 2016 6 6 https://doi.org/10.1051/emsci/2017006 03 October 2017 This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## 1 Introduction

The Gaussian cannon is a popular device for teaching since it is easy to construct and affordable. It is constructed as seen in Figure 1. A line of steel balls are attached to a strong magnet. If an additional ball collides with the magnet the outermost ball can be ejected at high speed. In the presence of a magnetic field the colliding ball, or trigger ball, loses magnetic potential energy and gains kinetic energy. By putting more balls on the firing side of the magnet the potential energy lost by the trigger ball is made smaller than the potential energy gained by the ejected ball. Hence the ejected ball can reach higher speeds than the trigger ball. However, due to energy losses in collisions the ejection speed will not increase indefinitely. Thus the Gaussian cannon is a good example of a system which can be optimized.

Though the Gaussian cannon is easily constructed, a thorough analysis is far from trivial due to the sheer amount of parameters. This report focuses on determining the speed of the ejected ball as a function of system parameters and how to maximize this speed. Adjustable parameters include the number of balls on the firing side of the cannon, the size of the magnet, the radius of the balls, and the incoming speed of the trigger ball. Earlier works have measured the force on the steel balls as a function of distance to compute the potential energies of the trigger ball and the ejected ball, thereby estimating the gained kinetic energy of the system [1,2]. In this report we present a more detailed analysis of the firing process and a novel approach to determine coefficients of restitution. The structure is IMRaD plus model, dead end, and conclusion. Dead end contains discarded experiments and lines of thought that was not useful. Fig. 1 The firing process of a Gaussian cannon. (a) A trigger ball moves towards a magnet and a line of steel balls. (b) The trigger ball collides with the magnet at speed vin. The collision energy is transferred through the magnet and the balls in line. (c) The rightmost ball is ejected with an initial speed vout. (d) The ejected ball reaches a speed of v1.

## 2 Model

The firing process is modeled as follows: the balls are subjected to magnetic force from the magnet, gravitational force from the earth and frictional and normal force from the track. Figure 1 presents the four different steps we consider in our model. (a) The trigger ball is attracted by the magnet, and rolls towards the bottom of the potential well with an initial velocity v0. Assuming rolling without slipping, the sum of the translational and rotational energy is and the magnetic potential energy is 0. (b) As the target ball hits the magnet, its magnetic potential energy is E0 and its velocity is vin. If the ball is still rolling without slipping, its rotational energy is . However, the ball can also start sliding during its approach, in which case the rotational energy would be small than this value. Assuming that the rotation can not decrease as compared to the initial motion, a lower bound of the rotational energy is given by . (c) As the kinetic energy of the target ball is transferred to the firing ball by collisions, the final ball of the line, with magnetic potential energy En, is ejected. Its velocity vout is expressed as a function of the total restitution coefficient vout = εtot vin. As no torque is applied to it, the eject ball has initially no rotation speed. (d) The ejected ball leaves the area of influence of the magnet, and is forced to roll by friction forces on the rail. We will assume that the energy lost in the process is negligible.

In summary, we obtain the following system of equations from conservation of energy and the definition of the coefficient of restitution: (1) (2) (3)

with describing the unknown rotational energy of the trigger ball. corresponds to a no-slip rotation, i.e. vin = in. The lowest possible value for α is (no angular acceleration). From equation (1) and ω0 < ωin, we can prove Therefore, if then v0 ≈ vin and . This will also be proven experimentally. Until then, the composite formula for v1 (v0 ) is: equations (1), (2), and (3) together give: (4)

This relation suggests a linear relation between the squared velocities, and depends on three parameters describing the kinetic energy of the trigger ball α, the energy loss occurring during the energy transfer ε and the potential energy experienced by the balls En. In the following, we will estimate experimentally these terms to gain understanding on the problem, and verify the linear relationship between and .

## 3 Methods

### 3.1 Velocity measurement

The experimental setup for velocity measurement consisted of a Gaussian cannon on a wooden track (of species fir). In order to measure speed, lasers were directed across the track and an oscilloscope registered whenever a ball interrupted a laser beam, and output time between interruptions. As the speed was expected to decrease due to friction, several lasers were used. Figure 2 presents an annotated picture of the setup used. Two laser beams crossed the path at five different points in total, the closest being 10 mm from the to-be ejected ball, i.e. at x = 10 mm. The speed was obtained from v = (x − x)/(t − t) where x − x denotes the distance between two nearby beams and t and t denotes the corresponding times as the beams were disrupted by a fired ball. As x − x is constant and t − t is measured with nanosecond precision, the measurement error is dominated by the random nature of the firing. How early the ball touches the track appears to matter, and we could not fully control this. To combat this, all experiments with controllable parameters were repeated five times. For such measurements, the standard error of mean is displayed as error bars.

The magnet was made of several connected cylindrical neodymium magnets of radius 10 mm and height 2.5 mm, magnetized along the symmetry axis . Both big and small balls were used, with ball diameter 12.0 mm and 9.5 mm (errors being less than 1 μm) and mass 7.10 g and 3.51 g respectively (errors being less than 0.01 g). The wooden track was U-shaped with slightly larger radius than the balls. Hot melt adhesive was used to fasten the magnet. Fig. 2 An experimental setup used to measure speed. The Gaussian cannon is resting on a wooden track. A trigger ball is approaching the cannon from the right and on the left lasers are arranged in order to measure the firing velocity.

### 3.2 Force measurement and energy determination

To determine the energy – En required to fully dislodge a ball from the system, consisting of a magnet and n other balls stuck to it, the attractive force F(d) was measured for varying distances d. En could then be calculated by integrating Fn (d) from 2rn to . The actual measurement was done in two slightly different ways.

For large distances, d > 3 mm, the ball was stuck to a nonmagnetic horizontal plate, whose height y could be adjusted with precision Δy < 0.01 mm. Beneath the plate was the magnet and the other balls, anchored to a weight, making the total weight M. This system rested upon a scale, showing a value m with precision Δm < 0.01 g. The attractive force F was calculated as F = (M − m) g. This method has good precision for both d and F, but for d < 3 mm F becomes so great that the plate bends, and therefore the precision of y and d becomes too low.

For small distances, d < 3 mm, tape strips with a total thickness of d were pasted onto the cannon and the ball was initially resting on top of the stacked strips. The cannon was pressed against the scale using heavy weights. The ball was then pulled by elevating the same plate as above as shown in Figure 3. As the plate was elevated, the displayed mass m gradually lowered, until it hit a certain mc, at which the ball snapped loose, causing m to rapidly increase. F was calculated as the maximum decrease in weight, F = (M − mc) g. This method has great precision for d, but a larger error for F, ΔF < 0.1 N. Fig. 3 The measurement setup for the force applied by magnets on a single steel ball. The ball is suspended above the magnets, which are attached to weights. The reciprocal force working on the magnets reduce the weight displayed on the scale, allowing for calculation of the magnetic force as a function of the difference in displayed weight. Themethod by which the ball was suspended in this image is the one we used for very small distances (d < 3 mm) between theball and the magnets.

### 3.3 COR determination

The COR εtot is defined as These velocities are hard to measure, since the magnetic field accelerates the balls rapidly over a short distance very close to the magnet or close to the line of balls. Therefore we measured v0 and v1 and used equations (1) and (3) to obtain vin and vout. This was done for number of balls n = 1, …, 7. The number of cylinder magnets was also varied, in the interval 3–9.

## 4 Results

### 4.1 Force and energy

We generally define En as the magnetic potential energy between the magnet and a firing ball having n balls between itself and the magnet. For Figure 1, n = 1. Note that we define the magnetic potential energy terms E0, E1, … so that they are always negative.

The force Fn(d) was measured for small and big balls. For brevity, only the plot for big balls (Fig. 4) is shown. The plots were qualitatively similar. The energies En in Table 1 was determined by integrating the cubic interpolant of Fn(d). Higher n results in less energy being needed to escape the cannon, which increases v1. Fig. 4 Pull forces F0, F1, and F2 for bigballs, as a function of distance d between themagnet and the border of the ball. Note that d is definedfrom the magnet itself, which is why Fn is defined only for d ≥ n ⋅ 12 mm. En was obtained by integrating the cubic interpolant (solid line).
Table 1

Values of En in mJ for big and small balls.

### 4.2 Velocity

When measuring the final velocity v1, the following parameters were varied:

• number of magnets;

• number of balls;

• ball size;

• trigger ball velocity v0.

Figure 5 shows the two first parameters being varied. The 13 magnets data set shows that ejection velocity does not increase monotonously with the number of magnets. The optimal number of magnets was found to be 6, adding more causes too much energy dissipation to be worth the extra magnetic flux. All following cannons will therefore consist of 6 magnets, in an effort to maximize velocities.

The influence of ball size can be seen in Figure 6. A smaller ball size allows higher v1, and increase the optimal n. Small balls were used for most following experiments.

Finally, the parameter v0 was varied, for virtually every data point in Figure 6. The result was qualitatively similar for every case, and therefore only one such graph is presented (Fig. 7). Fig. 5 Ejection velocities v1 for small balls with a v0 as small as possible, as a function of the number of balls and magnets. Fig. 6 Ejection velocities v1 for balls of various sizes when v0 = 0, as a function of the number of balls in-between the magnet and the ejected ball n. n = 0 is left out, as no ejection occurs. For small balls, n = 5 is optimal. For big balls, its n = 1, the smallest possible n.

### 4.3 COR

log  (εtot ) as a function of n is shown in Figure 8, with six magnets. The exponential relation suggests that , where εbm is the COR between ball and magnet, and εbb is the COR between balls. Table 2 shows εbb and εbm for different number of magnets. The value of α affects how vin is calculated from v0, and in turn εtot. Higher n results in more energy being dissipated, which lowers v1. Fig. 8 The logarithm of COR values are shown on the y-axis and the number of balls are shown on the x-axis. The ball diameter was 12 mm, the number of magnets was 6 and α = 0.7. The linear regression, y = −0.12 n − 0.08, is shown as a solid red line.
Table 2

COR values obtained from linearization of Figure 8, under the assumption .

## 5 Discussion

### 5.1 n dependence

The general qualitative dependence on n was as predicted by theory. For small n, the ejected ball transfers a significant amount of energy while escaping the magnet. For large n, a significant amount of energy is dissipated during the collisions between balls.

### 5.2 Linearity

Another strength of the model was the predicted (assuming constant α) and measured linearity between and . Equations (1), (2), and (3) together gives a linear relation between v1 and v0, explicitly stated in equation (4). The linear relation between the squared velocities was tested for big and small balls with n = 3, 4, …, 7. Figure 7 shows the test for small balls with n = 5. The linear relation hold Fig. 7 Ejection velocity squared as a function of trigger velocity , for small balls and n = 5. Each circle corresponds to one ejected ball. The linear relation is a general feature, regardless of choice of ball size and n. Chained cannon: When the cannons are coupled together, the speed squared follows a convergent series, represented by the green dashed line. With a large number of cannons, the final speed corresponds to the purple intersection cross.

### 5.3 Chained cannon

When two cannons are coupled together such that v1 of the first one becomes v0 of the second, the result is a better cannon. With an arbitrarily long chain of cannons, what is the final velocity?

It is known that for a single cannon, is a linear function of , in particular equation (4). In general, the sequence xn+1 = Axn + B converges to if |A| < 1. For our system, . This is positive and smaller than 1, because εtot < 1 and 2α > 1. Figure 7 shows the convergence of the series To test this quantitatively, we compare the velocities from intersection points between linear regressions and v0 = v1 (seen as the purple dot in Fig. 7) to calculated values, . The comparison can be seen in Figure 9. It not only gives us the expected value α = 0.7, but can also be seen as a test of equation (4), if α = 0.7 is assumed.

Figure 9 suggests that α = 0.7 is a good approximation, which would allow equations (1) and (4) to be simplified. Analysis of the optimal α values show that they have mean 0.68 and estimated SD 0.028. Thus 0.7 is 0.74 SD from the α distribution, which corresponds to failing to discard null hypothesis α = 0.7 with p < 0.5. In summary, we can accept the approximation α = 0.7. An interesting observation of Figure 9 is that the optimal n is lower than if a single cannon were used. This is due to the increased importance of energy dissipation in a chained system. Fig. 9 Comparison between intersection points from linear regressions (yellow line in Fig. 7), and intersection points from our model (fix points from Eq. (4)). Since equation (4) depends on α, which is not yet known, several plots are made depending on the choice of α.

## 6 Dead end

Early on, we considered comparing measured pull force to calculated force from B-field equations. This approach was discarded, due to the increased complexity and the magnetic dipole moment being unknown.

We tried filming the ejected ball to learn about its rotation, but the time resolution was not good enough for quantitative analysis. Sometimes balls rotated around the vertical axis due to contact with the side of the firing track, which made it impossible to determine angular velocity.

## 7 Conclusion

The Gaussian canon is a simple device that constitutes an insightful optimization problem. The ejection of the firing ball is due to the large kinetic energy transferred by collisions from the trigger ball, and to the limitation of the magnetic potential energy by the spacing balls. While a large number of spacing balls decreases the depth of the potential well from which the firing ball has to escape, it also increases the energy loss occurring during the energy transfer. This two competing behaviors result in the existence of an optimal number of spacing balls.

We presented a simple model based on energy conservation to relate the velocity of the firing ball to that of the trigger ball, and found a linear relation between this two quantities squared. Our model relies on three parameters, namely a rotational kinetic term α, the coefficient of restitution εtot and the magnetic potential energy En, of which the latter two were investigated experimentally. For α, we tested all possible values, and found in each case that with α = 0.7, the model agreed with experiments. That α = 0.7 also brought better understanding to the physics of the ejection process. The linear behavior was verified by the experiments, and shed light on the terminal velocity reached by a succession of Gaussian canon. This constitutes the main result of this work.

In addition to the optimization study presented here, the Gaussian canon exhibits numerous features that deserve attention. For instance, the symmetry between the number of trigger and firing ball (one and one in this work) can be seen as a generalization of Newton’s cradle, but also studied as the propagation of a soliton wave through the spacing balls . The inelastic collisions leading to finite losses in the coefficient of restitution also constitutes a problem of interest. We therefore believe that the Gaussian canon is a powerful tool not only to fire metallic balls, but also to teach and understand physics.

## References

1. J.A. Rabchuk, The gauss rifle and magnetic energy, Phys. Teach. 41, 158 (2003) [CrossRef] [Google Scholar]
2. D. Kagan, Energy and momentum in the gauss accelerator, Phys. Teach. 42, 24 (2004) [CrossRef] [Google Scholar]
3. Neodymmagnet, product number 31-1467, http://www.clasohlson.com [Google Scholar]
4. The sine-gordon model and its applications (2014), http://www.springer.com/in/book/9783319067216 [Google Scholar]

Cite this article as: Åke Andersson, Carl-Joar Karlsson, Hampus Lane, The Gaussian cannon, Emergent Scientist 1, 6 (2017)

## All Tables

Table 1

Values of En in mJ for big and small balls.

Table 2

COR values obtained from linearization of Figure 8, under the assumption .

## All Figures Fig. 1 The firing process of a Gaussian cannon. (a) A trigger ball moves towards a magnet and a line of steel balls. (b) The trigger ball collides with the magnet at speed vin. The collision energy is transferred through the magnet and the balls in line. (c) The rightmost ball is ejected with an initial speed vout. (d) The ejected ball reaches a speed of v1. In the text Fig. 2 An experimental setup used to measure speed. The Gaussian cannon is resting on a wooden track. A trigger ball is approaching the cannon from the right and on the left lasers are arranged in order to measure the firing velocity. In the text Fig. 3 The measurement setup for the force applied by magnets on a single steel ball. The ball is suspended above the magnets, which are attached to weights. The reciprocal force working on the magnets reduce the weight displayed on the scale, allowing for calculation of the magnetic force as a function of the difference in displayed weight. Themethod by which the ball was suspended in this image is the one we used for very small distances (d < 3 mm) between theball and the magnets. In the text Fig. 4 Pull forces F0, F1, and F2 for bigballs, as a function of distance d between themagnet and the border of the ball. Note that d is definedfrom the magnet itself, which is why Fn is defined only for d ≥ n ⋅ 12 mm. En was obtained by integrating the cubic interpolant (solid line). In the text Fig. 5 Ejection velocities v1 for small balls with a v0 as small as possible, as a function of the number of balls and magnets. In the text Fig. 6 Ejection velocities v1 for balls of various sizes when v0 = 0, as a function of the number of balls in-between the magnet and the ejected ball n. n = 0 is left out, as no ejection occurs. For small balls, n = 5 is optimal. For big balls, its n = 1, the smallest possible n. In the text Fig. 8 The logarithm of COR values are shown on the y-axis and the number of balls are shown on the x-axis. The ball diameter was 12 mm, the number of magnets was 6 and α = 0.7. The linear regression, y = −0.12 n − 0.08, is shown as a solid red line. In the text Fig. 7 Ejection velocity squared as a function of trigger velocity , for small balls and n = 5. Each circle corresponds to one ejected ball. The linear relation is a general feature, regardless of choice of ball size and n. Chained cannon: When the cannons are coupled together, the speed squared follows a convergent series, represented by the green dashed line. With a large number of cannons, the final speed corresponds to the purple intersection cross. In the text Fig. 9 Comparison between intersection points from linear regressions (yellow line in Fig. 7), and intersection points from our model (fix points from Eq. (4)). Since equation (4) depends on α, which is not yet known, several plots are made depending on the choice of α. In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.