Skip to main content

ORIGINAL RESEARCH article

Front. Phys., 06 March 2020
Sec. Statistical and Computational Physics
This article is part of the Research Topic Analytical and Numerical Methods for Differential Equations and Applications View all 11 articles

The Falling Body Problem in Quantum Calculus

  • Department of Mathematics, Faculty of Sciences, University of Tabuk, Tabuk, Saudi Arabia

The quantum calculus, q-calculus, is a relatively new branch in which the derivative of a real function can be calculated without limits. In this paper, the falling body problem in a resisting medium is revisited in view of the q-calculus to the first time. The q-differential equations describing the vertical velocity and distance of the body are obtained. Accordingly, exact expressions for the vertical velocity and the vertical distance are provided. The solutions are expressed in terms of the small q-exponential function which is an elementary function in the q-calculus. The dimensionality of the obtained formulae of the velocity and the distance are also analyzed. In addition, the present exact solutions reduce to the corresponding solutions in classical Newtonian mechanics when the quantum parameter q tends to one.

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 [712].

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 [1317] 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])

[n]q=1-qn1-q,    (1)

and as q → 1, we have

limq1 [n]q=n.    (2)

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

[n]q!=[1]q×[2]q×[3]q××[n]q.    (3)

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

Dqf(t):=dqf(t)dqt=f(t)-f(qt)(1-q)t,   t0,    (4)

such that

limq1Dqf(t)=f(t),    (5)

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

Dqf(0)=limt0Dqf(t).    (6)

According to (4) we have

Dqtn=[n]qtn-1.    (7)

The small q-analog of the exponential function et denoted by eq(t) (also called the small q-exponential function) is given as

eq(t)=j=0tj[j]q!.    (8)

The definite Jackson q-integral is defined by

0xf(t)dqt=(1-q)xj=0qjf(qjx),    (9)

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

0xDqf(t)dqt=f(x)-f(0).    (10)

The indefinite Jackson q-integral of the small q-exponential function eqt) is given as [18]

eq(αt)dqt=1αeq(αt)+c,    (11)

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 v0. The classical equation of motion for the particle is given by [15, 16]

mdvdt=-mg-mkv,    (12)

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

v(0)=v0,   z(0)=h,    (13)

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

dqvdqt:=-g-kv,      q(0,1].    (14)

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

v(t)=n=0antn,    (15)

and therefore

dqvdqt=n=0[n]qantn-1,         =n=1[n]qantn-1,   where  [0]q=0,         =n=0[n+1]qan+1tn.    (16)

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

n=0[n+1]qan+1tn=-g-kn=0antn,    (17)

or

[1]qa1+n=1[n+1]qan+1tn=-g-ka0-kn=1antn,    (18)

which gives

        a1=-g-ka0[1]q,an+1=-kan[n+1]q,   n1,    (19)

From (19), we have

a2=-ka1[2]q=(-1)2kg+(-k)2a0[1]q[2]q,a3=-ka2[3]q=(-1)3k2g+(-k)3a0[1]q[2]q[3]q,a4=-ka3[4]q=(-1)4k3g+(-k)4a0[1]q[2]q[3]q[4]q,..an=(-1)nkn-1g+(-k)na0[1]q[2]q[3]q[n]q,   n1.    (20)

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

an=(-1)nkn-1g+(-k)na0[n]q!,   n1.    (21)

The instantaneous velocity is obtained as

v(t)=a0+n=1antn,        =a0+n=1[(-1)nkn-1g+(-k)na0[n]q!]tn.        =a0+n=1[(g/k)(-kt)n+(-kt)na0[n]q!],    (22)

which can be written as

v(t)=a0+(gk+a0)n=1(-kt)n[n]q!.    (23)

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

v(t)=a0+(gk+a0)[eq(-kt)-1].    (24)

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

v(t)=v0+(gk+v0)[eq(-kt)-1],    (25)

which can be simplified as

v(t)=-gk+(gk+v0)eq(-kt).    (26)

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

Dqz(t)=-gk+(gk+v0)eq(-kt),    (27)

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

0tDqz(τ)dqτ=0t(-gk)dqτ+(gk+v0)0teq(-kτ)dqτ,    (28)

and hence,

z(t)-z(0)=-gk[τ[1]q]0t+(gk+v0)[-eq(-kτ)k]0t,    (29)

or

z(t)=h-gk(t[1]q)+(gk+v0)(-eq(-kt)k+1k),    (30)

i.e.,

z(t)=h-gtk+1k(gk+v0)(1-eq(-kt)),    (31)

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−2, then the vertical velocity in (26) must has dimension ms−1 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 eq(−kt) reduces to the standard exponential function ekt in classical calculus. Hence, (26) becomes

v(t)=-gk+(v0+gk)e-kt,    (32)

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

z(t)=h-gtk+1k(gk+v0)(1-e-kt),    (33)

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

v(t)|k0=limk0[v0e-kt+g(e-kt-1k)],                     =v0+glimk0(e-kt-1k),                     =v0+glimk0(-te-kt1),                     =v0-gt.    (34)

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

z(t)|k0=h+limk0[-gtk+(g+kv0)(1-e-kt)k2],                     =h+limk0[v0-gt+[(g+kv0)t-v0]e-kt2k],                     =h+limk0[-[(g+kv0)t-v0]te-kt+v0te-kt2],                     =h+(-(gt-v0)t+v0t2),                     =h+v0t-12gt2.    (35)

Here, it should be noted that L'Ho^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 eq(−kt) and (1 − eq(−kt)) as indicated below:

[kt]=[k]×[s]=s-1×s=Scalar,[eq(-kt)]=Scalar,[1-eq(-kt)]=Scalar.    (36)

By this, eq(−kt) and (1 − eq(−kt)) are dimensionless quantities, i.e., eq(−kt) and (1 − eq(−kt)) are scalar quantities. Accordingly, v(t) in (26) always has dimension ms−1 for all values of the quantum parameter q. Also z(t) in (31) always has dimension mq ∈ (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

Keywords: q-calculus, q-series, exact solution, falling body problem, q-differential equation

2010 Mathematics Subject Classification:35A24, 49K15

Citation: Alanazi AM, Ebaid A, Alhawiti WM and Muhiuddin G (2020) The Falling Body Problem in Quantum Calculus. Front. Phys. 8:43. doi: 10.3389/fphy.2020.00043

Received: 02 November 2019; Accepted: 13 February 2020;
Published: 06 March 2020.

Edited by:

Jesus Martin-Vaquero, University of Salamanca, Spain

Reviewed by:

Praveen Agarwal, Anand International College of Engineering, India
Ding-fang Zeng, Beijing University of Technology, China

Copyright © 2020 Alanazi, Ebaid, Alhawiti and Muhiuddin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Abdelhalim Ebaid, aebaid@ut.edu.sa

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.