1 Introduction

Magnet falling through a conducting pipe is generally demonstrated in almost every classroom to amuse the students about the effect of eddy current. Here an experiment is described to analyze the damping fall of a magnet through a metallic pipe using the ExpEYES17 data logger. The fall of a magnet inside a conducting pipe is damped by the opposing force due to the eddy current generated in the conducting pipe. This phenomenon is demonstrated in many places and many people have described the mathematical model of the fall of the magnet through the conducting pipe [1,2,3,4,5,6,7,8,9,10,11,12]. Further, this magnetic damping is exploited in a diversity of important applications including advanced car suspension systems that use a colloidal suspension, the smooth braking. Here we are presenting a case study by verifying the theoretical models by finding the drag coefficient and the conductivity of a conducting pipe using the previously derived equations.

2 Theory

The basic theory is explained in many papers [1,2,3,4,5,6,7,8,9,10,11,12], though a short description is given here for understanding the problem. When a permanent magnet is released through a metallic pipe, the magnetic flux through the metallic surface changes as it falls. According to Faraday’s Law, the flux change induces an electromotive force and then an electric current through the surface of the metallic pipe. This electric current produces a magnetic field to oppose the action that induced the emf. This repulsion yields a damping force on the magnet, so the magnet drifts through the pipe. The action of the forces on the magnet is shown in Fig. 1. A generally accepted model of velocity-dependent resistive forces can be written as Eq. (1) [13].

$${\text{F}} = {\text{k}}_{1} {\text{v}} + {\text{k}}_{2} {\text{v}}^{2}$$
(1)
Fig. 1
figure 1

Free body diagram of the falling magnet in a metallic pipe

In our case, the velocity is assumed to be small, therefore the term ‘v2’ will be negligible and it is ignored here. Further, the air resistance is also neglected. Therefore, two forces are acting on the magnet, one is the gravity which pulls the magnet down and the other is the induced magnetic force which tries to slow down the falling magnet.

The net force on the magnet can be calculated by summing the downward force and the dragging force. When ‘m’ and ‘v’ are the mass and velocity of the magnet respectively, the governing equation of motion can be written as;

$${\text{m}}\frac{{\rm dv}}{{\rm dt}} = {{\rm mg}} - {{\rm kv}}$$
(2)

where ‘g’ is the acceleration due to gravity (value of g is 9.8 m s−2) and ‘k’ is the damping coefficient. When ‘v0’ is the initial velocity then the velocity will be found by integrating Eq. (2) as;

$${\text{v}}\left( {\text{t}} \right) = \left( {{\text{v}}_{0} - \frac{{\rm mg}}{{\rm k}} } \right){{\rm e}}^{{ - \frac{{\rm k}}{{\rm m}}{{\rm t}}}} + \frac{{\rm mg}}{{\rm k}} = \frac{{{{\rm dx}}\left( {{\rm t}} \right)}}{{\rm dt}}$$
(3)

Then integrate the Eq. (3) to get x(t),

$${\text{x}}\left( {\text{t}} \right) = \left( { \frac{{{{\rm kv}}_{0} {{\rm m}} - {{\rm m}}^{2 } {{\rm g}}}}{{{{\rm k}}^{2} }}} \right) - \left( { \frac{{{{\rm kv}}_{0} {{\rm m}} - {{\rm m}}^{2 } {{\rm g}}}}{{{{\rm k}}^{2} }}} \right){{\rm e}}^{{ - \left( {\frac{{\rm k}}{{\rm m}}{{\rm t}}} \right)}} + \frac{{\rm mg}}{{\rm k}}{{\rm t}}$$
(4)

By rearranging the equation, we get

$$\frac{{\rm x}}{{{{\rm m}}^{2 } {{\rm g}}}}{{\rm k}}^{2 } - \left( {\frac{{{{\rm v}}_{0} }}{{\rm mg}} + \frac{{\rm t}}{{\rm m}}} \right) {\text{k}} + 1 = \left( {1 - \frac{{{{\rm v}}_{0} {{\rm k}}}}{{\rm mg}}} \right){{\rm e}}^{{ - \left( {\frac{{\rm t}}{{\rm m}}} \right) {{\rm k}}}}$$
(5)

Damping coefficient, ‘k’ can be found by solving the above Eq. (5) by substituting the corresponding values of the distance ‘x’ (length of the pipe), mass of the magnet ‘m’, initial velocity ‘v0’, the value of ‘g’ and the time taken by the magnet to pass through the pipe ‘t’. The time ‘t’ to pass the length of the pipe will be obtained using ExpEYES17 data logger by connecting two coils fixed at each end of the pipe (for details http://expeyes.in). Equation (5) solved by using www.wolframalpha.com to find k values.

For finding the conductivity of metallic pipe, Iniguez et al. [7] developed an equation by connecting the conductivity and speed of the falling object as given below,

$$\upsigma = \frac{{\sqrt 2 \uppi^{2} {{\rm k}}}}{{2 {\varphi }_{0}^{2} }} \sqrt {\frac{{{{\rm c}}^{2} - {{\rm b}}^{2} }}{{{ \ln }^{3} \left( {\frac{{\rm c}}{{\rm b}}} \right)}} }$$
(6)

where ‘b’ and ‘c’ are the inner and the outer radius of the pipe respectively, and ‘k’ is the damping coefficient. Where ‘φ0’ is the maximum value of magnetic flux ‘φ’, and will be obtained by numerically integrating (trapezoidal rule) the Eq. (7) using Microsoft Excel sheet. Expeyes17 was used to collect the emf versus time data.

$$\epsilon = - \frac{{{{\rm d}}{\varphi }}}{{\rm dt}}$$
(7)

The initial velocity of the magnet was found by the time of flight of the magnet through the PVC of length 8.4 cm fixed above the experimenting pipes. The whole idea is depicted in the Fig. 2. The time of flight can be calculated theoretically using the Eq. (8) and experimental time was found by fixing two coils at each end of the pipe as given in the Fig. 3,

$${\text{t}}_{0} = \sqrt {\frac{{2{{\rm L}}_{1} }}{{\rm g}}}$$
(8)
Fig. 2
figure 2

Description of time and velocity

Fig. 3
figure 3

Schematic diagram of the experimental setup

Then the initial velocity of the magnet when the magnet entering the conductive pipe will be

$${\text{V}}_{0} = {\text{gt}}_{0}$$
(9)

The effect of air in the downward motion of the magnet when it is passing through a pipe was found by dropping the test magnet through a non-conducting PVC pipe of known length. Then the theoretical time of flight through the larger PVC, 25.3 cm long, can be determined by solving the quadratic equation

$${\text{L}}_{2} = {\text{vt}} + \frac{1}{2}{{\rm gt}}^{2}$$
(10)

3 Experimental setup

Damping coefficient was found by solving the Eq. (5) after substituting the corresponding values of ‘x’, ‘m’, ‘g’, ‘v0’ and ‘t’. For our experiments, we had used fixed length of metallic pipes (as obtained from market) and magnets of known weights as given in the Tables 1 and 2. So, finding the time covered by the test magnet to cross the metallic pipes helps us to find the damping coefficient in that pipe. Recording of time in these small pipes are tedious and difficult using ordinary stopwatches and other conventional methods. Accurate time recording is only possible by sensors controlled by microcontrollers. So, here a computer interfaced data acquisition tool ExpEYES17 was used to record accurate time along with a coil as a sensor. ExpEYES17 is a hardware and software framework for real-time measurement and data analysis. The sensor used in our experiment was a coil that senses the change of magnetic flux. The sensor coils were fixed at each end of the metallic pipe and on PVC, as the center of the coil touches exactly at the ends of the pipes. Then the coil wires were connected to the ExpEYES17 as given in the Induction Coil experiments section of ExpEYES17 software. Schematically the experimental setup is shown in Fig. 3 and the photograph of the experimental setup with pipes and magnets are shown in Figs. 4 and 5.

Table 1 Characteristics of the pipes
Table 2 Characteristics of the magnets
Fig. 4
figure 4

Photograph of the experimental setup

Fig. 5
figure 5

Pipes and magnets

Magnets of different weights and sizes were dropped through a plastic non-conducting pipe from 8.4 cm above from the end of the metallic pipes as shown in Fig. 3 by keeping the ExpEYES17 data logger in ON state. When the magnet passes through the coils fixed at each end of the metallic pipes, the time changing flux induces some emf. This emf versus time data were recorded using the ExpEYES17 data logger option. ExpEYES17 software has an option to plot this recorded data, however we have used Gnuplot software for finding exact time in between the zero crossing of the emf versus time plots generated by the two coils. The zero crossing or the midpoint of the resultant plots shows that the magnet reached the ends of the metallic pipe. For this experiment, pipes and magnets of different sizes and weights were used and their diameters were measured using a screw gauge. The measured values for different pipes are given in Table 1. Magnets length and mass are given in Table 2. Magnetic moments were found using deflection magneto meter setup and tabled in Table 2.

4 Results and discussion

The initial velocity of the magnet was found by fixing two coils at each end of the PVC pipe of length 8.4 cm. Observed time and initial velocities are given in Table 3. The time values were substituted in Eq. (9) for finding initial velocity ‘vo’. Before performing the experiment with conducting pipes, we found the time to pass the magnets through a non-conducting pipe (here PVC) of length 25.3 cm and 2.5 cm of internal diameter. The obtained time values from ExpEYES17 data were compared with theoretical time. Theoretical time was found out as 0.1313 s by using Eq. (8). The time data shows that there is a difference in theoretical and experimental time that may be due to the air drag in pipes.

Table 3 The initial velocity of the magnet

Theoretical and experimental time for PVC for the test magnets are given in Table 4. The calculated data given in Table 4 shows that the magnet feels some upward thrust due to the air present in the PVC pipe. From this information, we deduce that the test magnet feels both the magnetic drag and the air resistance when the magnet moves through metallic pipe. The experiment was repeated for different conducting pipes and the results are tabulated in Table 5.

Table 4 Initial velocity, time and drag coefficient of PVC (theoretical and experimental)
Table 5 Drag coefficient k of different pipes

The data is also plotted in the Fig. 6 with the moment per mass of magnet versus drag coefficient. The plot shows that the copper has the highest drag coefficient than the other tested materials. When comparing the ‘k’ values of copper, the plot shows that the ‘k’ value increases with the thickness of the pipe as well as the strength of the magnet. But, in the case of SS, the magnetic drag force is not considerable and it approached the air drag as seen in PVC when the mass of the magnet increased. This plot also shows that the arrangement of materials based on their electrical conductivity in decreasing order from top to bottom.

Fig. 6
figure 6

Variation of damping coefficient k with the ratio of the moment to the mass of the magnets

Here in our case, the strength of the magnet is also increasing with increasing the mass. From the Fig. 6, it is clear that when the strength of the magnet increases the damping coefficient ‘k’ increases exponentially and tending to a constant value in all the pipes we used except in the case of stainless steel pipe. From the data, it is deduced that the drag coefficient is not a constant value for a particular material, but it depends on many other quantities such as magnetic moment, conductivity of the material, inner and outer diameter of the pipe.

Terminal velocity of the magnets when it moves in a conducting pipe can also be found by these experiments. Terminal velocity was found by plotting velocity versus time data. Velocity versus time plot for Cu pipe (P1) is given in Fig. 7, other plots have not shown here to save space but the extracted data is given in Table 5. It is seen that the magnets attains the terminal velocity within 0.1 s in highly conducting pipes than the other. Further, from the plots we can say that the terminal velocity is almost independent of the magnet weight in the case of higher conducting pipes.

Fig. 7
figure 7

Velocity versus time plot for Cu pipe (P1)

Further it is observed from the velocity–time plot of different conducting pipes, that when the magnet falls into the pipe, their velocity is exponentially decreases by the action of opposing force due to eddy current and it attains the terminal velocity very quickly. It is clear that the time taken for attaining the terminal velocity is dependent on the conductivity of the material of the pipe.

The conductivity of the pipes is calculated using the Eq. (6) and tabulated in Table 7. For finding the conductivity using the equation, maximum flux (φ0) is required. Therefore, maximum flux (φ0) was found by integrating the induced emf values using trapezoidal rule. The induced emf and the integrated values were plotted for all magnets and found the φ0 value. Plots of two and four magnets are only shown in Fig. 8. Extracted values of maximum flux for each magnet from their corresponding plots are given in Table 6.

Fig. 8
figure 8

emf versus time and flux versus time plots for two and four magnets

Table 6 Maximum value of flux for each magnet

Table 7 shows the conductivity of the materials of the pipe by magnet damped fall method and they are compared with standard values. The calculated values are almost matching with the standard value of the tested materials. Though, it is understood from this experiment that this is not an efficient method to find the conductivity of materials.

Table 7 Calculated and standard value of conductivity

5 Conclusions

Experiments have been performed to test the previously reported theoretical model of the retarded fall of a cylindrical magnet inside a vertical conductive pipe and the results are discussed here. We have used previously derived equations to find the damping coefficient ‘k’ during the retarded fall of the magnet through a conducting pipe. Further, the conductivity of seven different pipes was also calculated using ‘k’ values and maximum flux through the pipes. From this experiment, we observed that the drag force depends on the conductivity of the pipe, its inner and outer radii and the magnetic moment of the magnet used. Further, it is seen that the velocity of the magnet exponentially decreases and attains the terminal velocity within seconds in electrically conducting pipes.