** Journal of Applied Mathematics and Physics** Vol.2 No.4(2014), Article ID:43723,9 pages DOI:10.4236/jamp.2014.24004

Bessel Function and Damped Simple Harmonic Motion

Masoud Asadi-Zeydabadi

Department of Physics, University of Colorado Denver, Denver, USA

Email: Masoud.Asadi-Zeydabadi@ucdenver.edu

Copyright © 2014 by author and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

Received 5 February 2014; revised 5 March 2014; accepted 9 March 2014

ABSTRACT

A glance at Bessel functions shows they behave similar to the damped sinusoidal function. In this paper two physical examples (pendulum and spring-mass system with linearly increasing length and mass respectively) have been used as evidence for this observation. It is shown in this paper how Bessel functions can be approximated by the damped sinusoidal function. The numerical method that is introduced works very well in adiabatic condition (slow change) or in small time (independent variable) intervals. The results are also compared with the Lagrange polynomial.

**Keywords:**Bessel Function; Simple Harmonic Motion; Damped Sinusoidal Function; Lengthening Pendulum; Spring-Mass System; Lagrange Polynomial

1. Introduction

A brief view of the graphs of Bessel and sinusoidal functions shows they are very similar. Bessel functions look like damped sinusoidal functions. Sinusoidal functions are well known for all of us and we have seen the foot prints of them almost everywhere. We knew them from trigonometry but Bessel functions are new for college students and seem more complicated and the students get familiar with them usually in differential equation. Bessel and sinusoidal functions are solution of Bessel and harmonic differential equations. We know these differential equations belong to the family of Sturm-Liouville equation. Bessel and sinusoidal functions are orthogonal function and they appear in the solution of some partial differential equations. The type of orthogonal function that appears in the solution depends on the geometry, physics (the form of the differential equation) and the boundary conditions. In spite of the similarity between them still for the students dealing with Bessel function is more difficult than sinusoidal function.

The purpose of this paper is using a pedagogical method to show the similarity between them through two physical examples. Basically this is an attempt to understand the mathematics through physics. This method is based on the similarity between the form of Bessel and sinusoidal functions and their similarity has been interpreted by these examples. This paper is not aiming to discuss or prove fundamental similarities or differences between these two groups of functions. Basically our observation shows Bessel functions behave like sinusoidal functions with decreasing amplitude and varying period. In this paper two examples are given to understand the root of these behaviors. These examples are the lengthening pendulum and spring-mass system with variable mass. Both length and mass in the pendulum and the spring-mass system respectively increase linearly with time. Finally for these examples the results of the exact solution (Bessel function) are compared with the approximation method (damped sinusoidal function).

In addition to this pedagogical method (physical perspective of Bessel equation) the damped sinusoidal function is a good numerical approximation for Bessel function. It is compared with the Lagrange polynomial fitting; this method provides results better than the Lagrange polynomial fitting.

2. The Lengthening Pendulum

The lengthening pendulum which is known also as Lorentz’s pendulum is similar to a simple pendulum with increasing length at constant rate, with initial length of [1] -[8] . Its equation of motion is

(1)

where is the pendulum angle relative to the vertical axis and is gravitational acceleration. By changing variable from to and by using the small angle approximation (1) becomes:

(2)

The solution of (2) is given by Bessel function as following

(3)

where. If and at then the constants are given by and where is the initial value of. If and are related to each other such that is a zero of then the solution is

(4)

where

To understand the solution in terms of the sinusoidal function, Equation (1) for the small angle approximation can be written as

(5)

where and. Notice that since is time dependent, and are functions of time. If an adiabatic condition (slow change) is considered or if small time intervals are chosen then and are almost constant and the equation of motion is similar to damped simple pendulum. Then its solution for under damped condition is where frequency of motion is given by and it is function of time. For small time interval the solution at can be written in terms of the conditions of the pendulum at:

. This solution shows how Bessel functions can be related to the damped sinusoidal solution.

The equation of motion can also be written in terms of which is a dimensionless variable as following:

(6)

In this case the damping coefficient is. Again under adiabatic condition the solution at the neighborhood of is, where, and for we have

. Therefore for large the period approaches to like sinusoidal function. By comparing with the exact solution (Bessel function) at the neighborhood of we have: where

and are given by the condition of the problem at. This is an observation based on the solution of the lengthening pendulum and it is not a mathematical proof and depends on two constant (and) that should be determined by the condition at.

3. Spring-Mass System with Linearly Increasing Mass

In this case the mass is increased in steady rate: where is the initial mass and is the rate of change of the mass. The same treatment as previous case has been used. The equation of motion is:

(7)

where is the spring constant. By changing variable from to in (7) we have

(8)

The solution of (8) is given by Bessel function as following:

(9)

where. The initial conditions of and give and

where is the value of at. If and are adjusted such that

is a zero of then the solution is

(10)

where Equation (7) can be written as

, (11)

where and. Since is time dependent, again and are functions of time. The same as the pendulum case for adiabatic condition or for the small time interval and are almost constant then the equation of motion is similar to damped harmonic motion. Then its solution for under damped condition is where angular frequency of the motion is and it is function of time.

The equation of motion in terms of which is dimensionless is

(12)

In this case the damping coefficient is. Again under adiabatic condition the solution at the neighborhood of is, where, and for we have

. For large the period approaches to like sinusoidal function. By comparing with the exact solution (Bessel function) at the neighborhood of we have where and

are given by the conditions of the system at. This is physical evidence based on observation and it is not a mathematical proof.

4. The Quadratic Lagrange Polynomial for Numerical Comparison

In general the damped sinusoidal function provides a good approximation for Bessel functions. It can be compared with the quadratic Lagrange polynomial fitting which is given by [9] [10]

(13)

For the quadratic case and three points, and, are needed.

5. Results and Discussion

In this section some results are shown for both cases. The numerical values are used are not based on any physical reason and they are used just for comparison of these two methods.

In both of these problems there are two independent variables (for pendulum and for spring-mass system). Change in the momentum is due to both variables. Appearance of the first derivatives (and for pendulum and spring-mass system respectively) make them different from the simple harmonic motion. This is because of the change of the momentum due to the second independent variable (for the pendulum and for the spring-mass system). Mathematically this first derivative plays the role of the damping term that causes the amplitude decay. In reality there is no drag force like air resistance in these two cases but they are not conservative system their energies depend on time [1] . The total mechanical energy is the sum of the kinetic and the potential energy and for the pendulum case the energy density (energy per unit mass) is

. By substitution of (4) into this expression the energy density is given by. The last (Bessel’s) term is proportional to (like kinetic energy). This term decays down and the energy density approaches to. Since increases with time the energy decreases. The power density is given by

. Again the Bessel’s term decays and the power density converges to where is the string tension at this limit (at large value of time). Therefore the exact solution shows the energy of the system deceases as we expected from the approximation method based on the damped oscillation. Figure 1 shows the energy density and power density for the case of ν = 0.2048 m/s. It shows power density is negative and it decreases with time and approaches to W/kg.

The amplitude decays as but is inversely proportional to therefore the rate of the change of the amplitude decreases. The angular frequency increases as increases and approaches to 1, , therefore the period is inversely related to (it decreases as increases). This means for large the amplitude and frequency converge to almost constant values or the solution behaves like simple harmonic motion. We know by itself is function of time (and) and it increases by time and consequently the period increases by time. This can be clarified by giving the frequencies in terms of (and) which are functions of time. The frequencies in terms of for the pendulum and the spring-mass system are and respectively. For large value of time (large and) they approach to and. That means they get smaller and the periods become larger as time increases.

Figure 2 shows the results of two methods. For these results the small time intervals (less than the local period) have been used. This is a good criterion for the time step size in adaptive numerical method. Figure 2 shows the approximation method (damped sinusoidal method) agrees very well with the exact solution (Bessel function) for small time intervals. One can observe from this figure the period increases with time.

Suppose the pendulum problem has been solved exactly for some initial conditions. This solution is given in (4) by Bessel functions and at a given the exact conditions of the motion, (and), can be found and substituted into, , to find and. Figure 3 shows this solution for

(the first zero of) with initial conditions of and. The solution corresponding to has been shown by a dot on the graph. At this point the exact and the approximation solutions should

Figure 1. Energy density and power density for the lengthening pendulum: u_{0} = 30.571 m/s, and

Figure 2. The results of exact (Bessel) and the approximation (damped sinusoidal) solution for the lengthening pendulum (m/s, and). Since the small time interval has been used the results of both methods are quite agree with each other (no differences can be seen).

Figure 3. The results of the exact solution, Bessel function (solid line), and the approximation method, damped sinusoidal solution (dot line), for the lengthening pendulum: l_{0} = 1 m, and rad/s. The solution corresponding to has been shown by a dot on the graph.

match with each other through the condition at (i.e. and). Using these conditions from the exact solution and impose them on the approximation (damped sinusoidal) method gives us the approximation solution at the neighborhood of this point. Figure 3 shows the results agree with each other particularly near this point. One can see the deviation of the approximation method from the exact solution away from this point especially for small. In this case we are not using small time interval and the figure shows the results over several periods.

Figure 4 shows the same results as Figure 3 but for the spring-mass system, in this case the exact solution is given in terms of. Again the results of the exact solution and the approximation method are quite similar. The deviation of these two methods at small is shown in Figure 5. Again we can observe for the results are matched very well.

In the Lagrange polynomial fitting three points are needed but for the damped sinusoidal function only the value of Bessel function and its derivative at a point are needed. This is a big advantage of this method compare to the Lagrange polynomial fitting. The Bessel function is oscillatory therefore the order of polynomial in a large interval should be higher and more points are needed. Figure 6 shows the comparison of these two techniques with the exact Bessel function.

The results from Figure 6 shows the damped sinusoidal function is a better approximation compare to the Lagrange polynomial fitting. It is obvious for the larger interval the quadratic polynomial is not a reasonable approximation and higher order polynomial with more points are needed. The results of this paper shows in long interval the damped sinusoidal is good numerical approximation.

6. Conclusions

The graphs of Bessel functions are similar to the damped sinusoidal solution. In this paper this observation has been investigated by two physical examples. The Bessel’s equation has been compared with the equation of the damped simple harmonic motion. The solutions of these methods are compared at the neighborhood of an arbitrary point. The results are shown the approximation method works very well particularly when u is large (i.e. for large value of independent variable for example time in this paper for two examples). The damped sinusoidal

Figure 4. The results of the exact solution (Bessel function, solid line) and the approximation method (damped sinusoidal solution, dot line) for the mass-spring system: m_{0} = 1 kg, u_{0} = 3.832, x_{0} = 1 m, u_{n} = 100, x_{n} = −0.0496 m and m/s. The solution corresponding to u_{n} = 100 has been shown by a dot on the graph.

Figure 5. The results of the exact solution (Bessel function, solid line) and the approximation (damped sinusoidal, dot line) method for the mass-spring system: m_{0} = 1 kg, u_{0} = 3.832, x_{0} = 1 m, u_{n} = 10, x_{n} = 0.6106 m and m/s. The solution corresponding to u_{n} = 10, has been shown by a dot on the graph.

Figure 6. Bessel function, x = J_{0}(u), the damped sinusoidal function, and the quadratic Lagrange polynomial. The given points for Bessel function, J_{0}(0.5) = 0.9385, J_{0}(2.0) = 0.2239 and J_{0}(4.0) = −0.3971, are shown in figure.

function not only is a good way to interpret the property of the Bessel equation but also is a good numerical approximation. The comparison with the Lagrange polynomial shows the numerical advantage of the damped sinusoidal function.

In this paper we use two physical examples but this method can be generalized for any Bessel’s equation. The general form of Bessel’s equation is. Again by using the appropriate initial value for the final solution can be given in terms of either or. The damped sinusoidal solution at the neighborhood of a given point is given by () where and. The values of Bessel function and its derivative at are used to find and. This method can be considered as a good perturbation method for the problems that are involved Bessel function. This method can improve numerical techniques for some particular equations. In the same way we can look at other functions like Legendre polynomial. The Legendre differential equation is given by with . We can write the differential equation in the form damping harmonic oscillation with and . In this case the solution is growing up and the increases by. Notice that the polynomial is defined in region of .

References

- Werner, A. and Eliezer, J.C. (1969) The Lengthening Pendulum. Journal of Australian Mathematical Society, 9, 331- 336. http://dx.doi.org/10.1017/S1446788700007254
- Littlewood, J.E. (1963) Lorentz’s Pendulum Problem. Annals of Physics, 21, 233-249. http://dx.doi.org/10.1016/0003-4916(63)90107-6
- Littlewood, J.E. (1964) Adiabatic Invariance III. The Equation. Annals of Physics, 29, 1-12. http://dx.doi.org/10.1016/0003-4916(64)90188-5
- Littlewood, J.E. (1964) Adiabatic Invariance IV: Note on a New Method for Lorentz’s Pendulum Problem. Annals of Physics, 29, 13-18. http://dx.doi.org/10.1016/0003-4916(64)90189-7
- Littlewood, J.E. (1964) Adiabatic Invariance V. Multiple Periods. Annals of Physics, 30, 138-153. http://dx.doi.org/10.1016/0003-4916(64)90307-0
- Brearley, M.N. (1966) The Simple Pendulum with Uniformly Changing String Length. Proceedings of the Edinburgh Mathematical Society (Series 2), 15, 61-66.
- Sánchez-Soto, L.L. and Zoido, J. (2013) Variations on the Adiabatic Invariance: The Lorentz Pendulum. American Journal of Physics, 81, 57. http://dx.doi.org/10.1119/1.4763746
- Boas, M.L. (2006) Mathematical Methods in the Physical Science. 3rd Edition, Wiley, 598-599.
- Gil, A., Segura, J. and Temme, N. (2007) Numerical Methods for Special Functions. SIAM. http://dx.doi.org/10.1137/1.9780898717822
- Garcia, A.L. (2000) Method for Physics. 2nd Edition, Prentice-Hall, NJ.