The Falling Body Problem in Quantum Calculus

1. Introduction

Basically, the regular calculus uses limits in calculating the derivatives of real functions. However, the calculus without limits is nowadays known as quantum calculus or q-calculus. Historically, in the eighteenth century, Euler obtained the basic formulae in q-calculus. However, Jackson [1] may have been the first to introduce the notion of the definite q-derivative and q-integral. Currently, there is a significant interest in implementing the q-calculus due to its applications in several areas, such as mathematics, number theory, and combinatorics [2]. Ernst [3, 4] pointed out that the majority of scientists who use q-calculus are physicists. Baxter [5] introduced the exact solutions of several models in Statistical Mechanics. Bettaibi and Mezlini [6] solved some q-heat and q-wave equations. Many interesting results in such area of research were also introduced by several authors in the literature [7–12].

In this paper, we aim to extend the applications of the q-calculus to study the falling body problem in a resisting medium. This problem and also the full projectile motion have been investigated by several authors [13–17] using various definitions in fractional calculus. However, the present paper may be the first to analyze the falling body problem in view of the q-calculus.

The basic formulae in q-calculus will be used to analyze the motion of a falling body in a resisting medium. Moreover, it will be shown that the exact solutions for the vertical velocity and distance reduce to the classical ones as q → 1. The paper is organized as follows. Section 2 presents the main aspects of the q-calculus. Sections 3 discusses the application of the q-calculus on the falling body problem. Section 4 includes an additional analysis. Finally, section 5 outlines the conclusions.

2. The Main Aspects of the q-Calculus

Let q ∈ ℝ and n ∈ ℕ, then [n]_q is defined as (first chapter in [18])

$\begin{array}{ll}{\left[n\right]}_q=\frac{1-q^n}{1-q},& (1)\end{array}$

and as q → 1, we have

$\begin{array}{ll}{\displaystyle \underset{q\to 1}{lim}}{\left[n\right]}_q=n.& (2)\end{array}$

The q-factorial [n]_q! of a positive integer n is given by

$\begin{array}{ll}{\left[n\right]}_q!={\left[1\right]}_q\times {\left[2\right]}_q\times {\left[3\right]}_q\times \cdots \times {\left[n\right]}_q.& (3)\end{array}$

The definition of q-differential is d_qf(t) = f(t) − f(qt) and the q-derivative of a function f(t) is defined by [18]

$\begin{array}{ll}{D}_{q}f(t):=\frac{d_{q}f(t)}{d_{q}t}=\frac{f(t)-f(qt)}{(1-q)t},t\ne 0,& (4)\end{array}$

such that

$\begin{array}{ll}{\displaystyle \underset{q\to 1}{lim}}{D}_{q}f(t)=f^{\prime }(t),& (5)\end{array}$

if f is differentiable at t, and we have at t = 0 that

$\begin{array}{ll}{D}_{q}f(0)={\displaystyle \underset{t\to 0}{lim}}{D}_{q}f(t).& (6)\end{array}$

According to (4) we have

$\begin{array}{ll}{D}_q{t}^n={\left[n\right]}_q{t}^{n-1}.& (7)\end{array}$

The small q-analog of the exponential function e^t denoted by e_q(t) (also called the small q-exponential function) is given as

$\begin{array}{ll}{e}_q(t)={\displaystyle \sum _{j=0}^{\infty }}\frac{t^j}{{\left[j\right]}_q!}.& (8)\end{array}$

The definite Jackson q-integral is defined by

$\begin{array}{ll}{\displaystyle {\int }_0^x}f(t)d_{q}t=(1-q)x{\displaystyle \sum _{j=0}^{\infty }}{q}^{j}f(q^{j}x),& (9)\end{array}$

and according to (4) and (9), we have

$\begin{array}{ll}{\displaystyle {\int }_0^x}{D}_{q}f(t)d_{q}t=f(x)-f(0).& (10)\end{array}$

The indefinite Jackson q-integral of the small q-exponential function e_q(αt) is given as [18]

$\begin{array}{ll}{\displaystyle \int }{e}_q(\alpha t)d_{q}t=\frac{1}{\alpha }{e}_q(\alpha t)+c,& (11)\end{array}$

where c is a real constant. The correctness of dimensionality of the physical quantities is actually guaranteed by the definition (4).

3. The Falling Body Problem

Consider the falling of an object of mass m in the Earth gravitational field through the air from a height h with initial velocity v₀. The classical equation of motion for the particle is given by [15, 16]

$\begin{array}{ll}m\frac{dv}{dt}=-mg-mkv,& (12)\end{array}$

where k is a positive constant and its dimensionality is the inverse of seconds, i.e., [k] = s⁻¹. The initial conditions are given as

$\begin{array}{ll}v(0)=v_0,z(0)=h,& (13)\end{array}$

where z(t) is the vertical distance of the particle at arbitrary time t and $\frac{dz(t)}{dt}=v(t)$. The equation of motion (12) in view of the quantum calculus becomes

$\begin{array}{ll}\frac{d_{q}v}{d_{q}t}:=-g-kv,q\in (0,1].& (14)\end{array}$

In order to solve Equation (14), we assume the solution in the series form:

$\begin{array}{ll}v(t)={\displaystyle \sum _{n=0}^{\infty }}{a}_n{t}^n,& (15)\end{array}$

and therefore

$\begin{array}{l}\frac{d_{q}v}{d_{q}t}={\displaystyle \sum _{n=0}^{\infty }}{\left[n\right]}_q{a}_n{t}^{n-1},\\ ={\displaystyle \sum _{n=1}^{\infty }}{\left[n\right]}_q{a}_n{t}^{n-1},\text{where }{\left[0\right]}_q=0,\\ ={\displaystyle \sum _{n=0}^{\infty }}{\left[n+1\right]}_q{a}_{n+1}{t}^n.& (16)\end{array}$

Substituting (15) and (16) into (14), yields

$\begin{array}{ll}{\displaystyle \sum _{n=0}^{\infty }}{\left[n+1\right]}_q{a}_{n+1}{t}^n=-g-k{\displaystyle \sum _{n=0}^{\infty }}{a}_n{t}^n,& (17)\end{array}$

or

$\begin{array}{ll}{\left[1\right]}_q{a}_1+{\displaystyle \sum _{n=1}^{\infty }}{\left[n+1\right]}_q{a}_{n+1}{t}^n=-g-k{a}_0-k{\displaystyle \sum _{n=1}^{\infty }}{a}_n{t}^n,& (18)\end{array}$

which gives

$\begin{array}{l}{a}_1=\frac{-g-k{a}_0}{{\left[1\right]}_q},\\ a_{n+1}=\frac{-k{a}_n}{{\left[n+1\right]}_q},n\ge 1,& (19)\end{array}$

From (19), we have

$\begin{array}{l}{a}_2=\frac{-k{a}_1}{{\left[2\right]}_q}=\frac{{(-1)}^{2}kg+{(-k)}^2{a}_0}{{\left[1\right]}_q{\left[2\right]}_q},\\ a_3=\frac{-k{a}_2}{{\left[3\right]}_q}=\frac{{(-1)}^3{k}^{2}g+{(-k)}^3{a}_0}{{\left[1\right]}_q{\left[2\right]}_q{\left[3\right]}_q},\\ a_4=\frac{-k{a}_3}{{\left[4\right]}_q}=\frac{{(-1)}^4{k}^{3}g+{(-k)}^4{a}_0}{{\left[1\right]}_q{\left[2\right]}_q{\left[3\right]}_q{\left[4\right]}_q},\\ .\\ .\\ a_n=\frac{{(-1)}^n{k}^{n-1}g+{(-k)}^n{a}_0}{{\left[1\right]}_q{\left[2\right]}_q{\left[3\right]}_q\dots {\left[n\right]}_q},n\ge 1.& (20)\end{array}$

This n-term coefficient can expressed in terms of the q-factorial [n]_q! as

$\begin{array}{ll}{a}_n=\frac{{(-1)}^n{k}^{n-1}g+{(-k)}^n{a}_0}{{\left[n\right]}_q!},n\ge 1.& (21)\end{array}$

The instantaneous velocity is obtained as

$\begin{array}{l}v(t)=a_0+{\displaystyle \sum _{n=1}^{\infty }}{a}_n{t}^n,\\ =a_0+{\displaystyle \sum _{n=1}^{\infty }}\left[\frac{{(-1)}^n{k}^{n-1}g+{(-k)}^n{a}_0}{{\left[n\right]}_q!}\right]t^n.\\ =a_0+{\displaystyle \sum _{n=1}^{\infty }}\left[\frac{(g/k){(-kt)}^n+{(-kt)}^n{a}_0}{{\left[n\right]}_q!}\right],& (22)\end{array}$

which can be written as

$\begin{array}{ll}v(t)=a_0+\left(\frac{g}{k}+a_0\right){\displaystyle \sum _{n=1}^{\infty }}\frac{{(-kt)}^n}{{\left[n\right]}_q!}.& (23)\end{array}$

In terms of the small exponential function e_q(−kt), we have

$\begin{array}{ll}v(t)=a_0+\left(\frac{g}{k}+a_0\right)\left[e_q(-kt)-1\right].& (24)\end{array}$

Applying the first initial condition in (13) on (24), we obtain a₀ = v₀ and therefore v(t) becomes

$\begin{array}{ll}v(t)=v_0+\left(\frac{g}{k}+v_0\right)\left[e_q(-kt)-1\right],& (25)\end{array}$

which can be simplified as

$\begin{array}{ll}v(t)=-\frac{g}{k}+\left(\frac{g}{k}+v_0\right)e_q(-kt).& (26)\end{array}$

The vertical distance z(t) in quantum calculus is governed by,

$\begin{array}{ll}{D}_{q}z(t)=-\frac{g}{k}+\left(\frac{g}{k}+v_0\right)e_q(-kt),& (27)\end{array}$

where v(t) = D_qz(t). Integrating (27), it then follows;

$\begin{array}{ll}{\displaystyle {\int }_0^t}{D}_{q}z(\tau )d_q\tau ={\displaystyle {\int }_0^t}\left(-\frac{g}{k}\right)d_q\tau +\left(\frac{g}{k}+v_0\right){\displaystyle {\int }_0^t}{e}_q(-k\tau )d_q\tau ,& (28)\end{array}$

and hence,

$\begin{array}{ll}z(t)-z(0)=-\frac{g}{k}{\left[\frac{\tau }{{\left[1\right]}_q}\right]}_0^t+\left(\frac{g}{k}+v_0\right){\left[-\frac{e_q(-k\tau )}{k}\right]}_0^t,& (29)\end{array}$

or

$\begin{array}{ll}z(t)=h-\frac{g}{k}\left(\frac{t}{{\left[1\right]}_q}\right)+\left(\frac{g}{k}+v_0\right)\left(-\frac{e_q(-kt)}{k}+\frac{1}{k}\right),& (30)\end{array}$

i.e.,

$\begin{array}{ll}z(t)=h-\frac{gt}{k}+\frac{1}{k}\left(\frac{g}{k}+v_0\right)\left(1-e_q(-kt)\right),& (31)\end{array}$

where [1]_q = 1. The exact solutions (26) and (31) should be reduced to the corresponding solutions in classical Newtonian mechanics when q → 1. In addition, if the acceleration due to gravity is measured in ms⁻², then the vertical velocity in (26) must has dimension ms⁻¹ and the vertical distance in (31) must has dimension m. These issues are addressed in the following section.

4. Analysis and Applications

First of all, we investigate the solutions (26) and (31) when q → 1. In this case, the small exponential function e_q(−kt) reduces to the standard exponential function e^−kt in classical calculus. Hence, (26) becomes

$\begin{array}{ll}v(t)=-\frac{g}{k}+\left(v_0+\frac{g}{k}\right)e^{-kt},& (32)\end{array}$

which is the analytic expression for velocity in the case of the classical Newtonian mechanics (see Equation 16 in reference [15]). Besides, the vertical distance in (31) reduces to

$\begin{array}{ll}z(t)=h-\frac{gt}{k}+\frac{1}{k}\left(\frac{g}{k}+v_0\right)\left(1-e^{-kt}\right),& (33)\end{array}$

which is also the analytic expression for the vertical distance in the classical Newtonian mechanics (see Equation 17 in reference [15]).

In addition, in the case of no air resistance, i.e., the parameter k vanishes, we obtain from (32) that

$\begin{array}{l}v(t){|}_{k\to 0}={\displaystyle \underset{k\to 0}{lim}}\left[v_0{e}^{-kt}+g\left(\frac{e^{-kt}-1}{k}\right)\right],\\ =v_0+g{\displaystyle \underset{k\to 0}{lim}}\left(\frac{e^{-kt}-1}{k}\right),\\ =v_0+g{\displaystyle \underset{k\to 0}{lim}}\left(\frac{-t{e}^{-kt}}{1}\right),\\ =v_0-gt.& (34)\end{array}$

Also, the vertical distance in (33) in the absence of air resistance becomes

$\begin{array}{l}z(t){|}_{k\to 0}=h+{\displaystyle \underset{k\to 0}{lim}}\left[\frac{-gtk+(g+k{v}_0)(1-e^{-kt})}{k^2}\right],\\ =h+{\displaystyle \underset{k\to 0}{lim}}\left[\frac{v_0-gt+\left[(g+k{v}_0)t-v_0\right]e^{-kt}}{2k}\right],\\ =h+{\displaystyle \underset{k\to 0}{lim}}\left[\frac{-\left[(g+k{v}_0)t-v_0\right]t{e}^{-kt}+v_{0}t{e}^{-kt}}{2}\right],\\ =h+\left(\frac{-(gt-v_0)t+v_{0}t}{2}\right),\\ =h+v_{0}t-\frac{1}{2}g{t}^2.& (35)\end{array}$

Here, it should be noted that L'H$\widehat{\mathrm{\text{o}}}$pital's rule was applied to calculate the above limits. The Equations (34) and (35) are the same of the corresponding equations for the vertical velocity and vertical distance in Newtonian mechanics in the absence of air resistance.

Regarding the dimensions of the q-forms of v(t) and z(t) in (26) and (31), respectively, it should be first to specify the dimensions of the quantities e_q(−kt) and (1 − e_q(−kt)) as indicated below:

$\begin{array}{l}\left[kt\right]=\left[k\right]\times \left[s\right]=s^{-1}\times s=\text{Scalar},\\ \left[e_q(-kt)\right]=\text{Scalar},\\ \left[1-e_q(-kt)\right]=\text{Scalar}.& (36)\end{array}$

By this, e_q(−kt) and (1 − e_q(−kt)) are dimensionless quantities, i.e., e_q(−kt) and (1 − e_q(−kt)) are scalar quantities. Accordingly, v(t) in (26) always has dimension ms⁻¹ for all values of the quantum parameter q. Also z(t) in (31) always has dimension m ∀q ∈ (0, 1]. The correctness of dimensions of the q-vertical velocity and the q-height was actually guaranteed by the definition (4) without any need to involve an auxiliary parameter as in the literature [15, 16].

Although the present model of the falling body problem seems simple, the authors believe that the current work is worthy of exploration. This is because the present solution was provided to the first time for the falling problem in view of q-calculus. In addition, it was shown in this paper the way of obtaining the solutions in exact forms and also how to check the dimensions of the physical quantities in terms of q-parameter. Furthermore, the obtained solutions can be verified by direct substitutions into the governing equations. Therefore, the present work is a first step for further studies in future to explore various physical models in applied mathematics implementing the q-calculus.

5. Conclusion

In this paper, the quantum calculus was applied to solve the falling body problem. The exact solutions for the q-vertical velocity and the q-distance have been obtained. The obtained exact solutions were expressed in terms of the small q-exponential function. The correctness of dimensionality of the obtained formulae of the velocity and the distance was proved. Moreover, The present exact solutions reduced to the corresponding solutions in classical Newtonian mechanics when the quantum parameter q tends to one. The present work can be further extended to explore the physical properties of the projectile motion in two and three dimensions in view of the q-calculus.

Data Availability Statement

All datasets generated for this study are included in the article/supplementary material.

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Funding

The authors extend their appreciation to the Deanship of Scientific Research at University of Tabuk for funding this work through Research Group no. RGP-0207-1440.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

1. Jackson FH. On a q-definite integrals. Q J Pure Appl Math. (1910) 41:193–203.

Google Scholar

2. Andrews GE. q-Series: Their Development and Applications in Analysis, Number Theory, Combinatorics, Physics and Computer Algebra. CBMS Regional Conference Series in Mathematics, Vol. 66. Providence, RI: American Mathematical Society (1986).

Google Scholar

3. Ernst T. The history of q-calculus and a new method (Licentiate thesis), U.U.D.M. Report 2000: 16. Available online at: https://pdfs.semanticscholar.org/895b/4f8b059dc7ad18807f5aabdf09f2158d0c66.pdf

Google Scholar

4. Ernst T. A method for q-calculus. J Nonlin Math Phys A Method. (2003) 10:487–525. doi: 10.2991/jnmp.2003.10.4.5

CrossRef Full Text | Google Scholar

5. Baxter R. Exact Solved Models in Statistical Mechanics. New York, NY: Academic Press (1982).

Google Scholar

6. Bettaibi N, Mezlini K. On the use Of the q-Mellin transform to solve some q-heat and q-wave equations. Int J Math Arch. (2012) 3:446–55.

7. Li YQ, Sheng ZM. A deformation of quantum mechanics. J Phys A Math Gen. (1992) 25:6779–88. doi: 10.1088/0305-4470/25/24/028

CrossRef Full Text | Google Scholar

8. Annaby MH, Mansour ZS. q-Fractional Calculus and Equations. Heidelberg; New York, NY: Springer (2012).

Google Scholar

9. Koca I, Demirci E. On local asymptotic stability of q-fractional nonlinear dynamical systems. Appl Appl Math. (2016) 11:174–83. Available online at: http://pvamu.edu/aam

Google Scholar

10. Rangaig NA, Pada CT, Convicto VC. On the existence of the solution for q-Caputo fractional boundary value problem. Appl Math Phys. (2017) 5:99–102. doi: 10.12691/amp-5-3-4

CrossRef Full Text

11. Tang Y, Zhang T. A remark on the q-fractional order differential equations. Appl Math Comput. (2019) 350:198–208. doi: 10.1016/j.amc.2019.01.008

CrossRef Full Text | Google Scholar

12. Chanchlani L, Alha S, Gupta J. Generalization of Taylor's formula and differential transform method for composite fractional q-derivative. Ramanujan J. (2019) 48:21–32. doi: 10.1007/s11139-018-9997-7

CrossRef Full Text | Google Scholar

13. Kwok SF. A falling body problem through the air in view of the fractional derivative approach. Phys A. (2005) 350:199–206. doi: 10.1016/j.physa.2004.11.041

CrossRef Full Text | Google Scholar

14. Ebaid A. Analysis of projectile motion in view of the fractional calculus. Appl Math Model. (2011) 35:1231–9. doi: 10.1016/j.apm.2010.08.010

CrossRef Full Text | Google Scholar

15. Garcia JJR, Calderon MG, Ortiz JM, Baleanu D. Motion of a particle in a resisting medium using fractional calculus approach. Proc Roman Acad Ser A. (2013) 14:42–7.

Google Scholar

16. Ebaid A, Masaedeh B, El-Zahar E. A new fractional model for the falling body problem. Chin Phys Lett. (2017) 34:020201. doi: 10.1088/0256-307X/34/2/020201

CrossRef Full Text | Google Scholar

17. Alharbi FM, Baleanu D, Ebaid A. Physical properties of the projectile motion using the conformable derivative. Chin J Phys. (2019) 58:18–28. doi: 10.1016/j.cjph.2018.12.010

CrossRef Full Text | Google Scholar

18. Kac VG, Cheung P. Quantum Calculus. New York, NY: Springer-Verlag (2002).

Google Scholar