Next Article in Journal
Some Remarks on Fuzzy sb-Metric Spaces
Previous Article in Journal
Fuzzy Multicriteria Decision Mapping to Evaluate Implant Design for Maxillofacial Reconstruction
Previous Article in Special Issue
Fractional Diffusion–Wave Equation with Application in Electrodynamics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ordinary Differential Equation with Left and Right Fractional Derivatives and Modeling of Oscillatory Systems

Institute of Applied Mathematics and Automation, Kabardino-Balkarian Scientific Center of Russian Academy of Sciences, 89-A Shortanov Street, 360000 Nalchik, Russia
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(12), 2122; https://doi.org/10.3390/math8122122
Submission received: 30 October 2020 / Revised: 19 November 2020 / Accepted: 26 November 2020 / Published: 27 November 2020
(This article belongs to the Special Issue Mathematical Modeling of Hereditarity Oscillatory Systems)

Abstract

:
We consider the principle of least action in the context of fractional calculus. Namely, we derive the fractional Euler–Lagrange equation and the general equation of motion with the composition of the left and right fractional derivatives defined on infinite intervals. In addition, we construct an explicit representation of solutions to a model fractional oscillator equation containing the left and right Gerasimov–Caputo fractional derivatives with origins at plus and minus infinity. We derive a representation for the composition of the left and right derivatives with origins at plus and minus infinity in terms of the Riesz potential, and introduce special functions with which we give solutions to the model fractional oscillator equation with a complex coefficient. This approach can be useful for describing dissipative dynamical systems with the property of heredity.

1. Introduction

The basic framework for Lagrangian and Hamiltonian mechanics is the principle of least action. In problems of mechanics, the principle of stationary action is most commonly used, and is the most important among all extreme principles. The equations of motion obtained using these principles are called the Euler–Lagrange equations. The Euler–Lagrange equation with the Lagrange function can describe any mechanical motion.
In recent years, fractional integro-differentiation has been widely used to model the properties of dynamic systems [1,2,3]. Fractional calculus enables one to naturally take into account the heredity, time non-locality and memory effects in dynamic processes [2].
The principle of least action in the context of the fractional calculus is studied in [4,5,6,7,8,9,10,11]. In [4,5], it was shown that applying the Lagrangian with fractional derivatives allows one to describe the dynamics of nonconservative systems. In [6], two fractional variational problems were considered. A generalized Euler–Lagrange equation with fractional differentiation was derived. The classical Euler–Lagrange equation is a special case of the latter. In addition, practical examples of the application of the results obtained were considered. It was shown that this approach would make it possible to simulate dynamic systems with greater accuracy. In [8], the Lagrange formalism with fractional derivatives with infinite limits was considered for the first time. In [9], the Hamiltonian and Lagrangian formulations with fractional differentiation were used to describe constrained systems. It was shown that both treatments for systems with linear velocities are equivalent. In [7,12,13] the principle resulted in very exotic differential equations of fractional order, with composition of the Riemann–Liouville or Caputo operators with different origins. These equations and their possible application to modeling dynamic processes can be found in [12,13,14,15,16,17].
We point out the major challenges related to the equations with different origins. First, methods for obtaining solutions to these equations, both analytically and numerically, are not yet well developed. Next, physical interpretation of the fractional Euler–Lagrange equations is currently not well understood. The first-order Euler–Lagrange equation does not contain the limits of integration in the action functional, unlike the fractional Euler–Lagrange equation. Therefore, the following question arises: is it possible to use a fractional Euler–Lagrange equation if the time variable falls outside the interval within which the action functional is considered? Any value can be chosen as a time interval, including a small physical time (for example, the Planck time, which is approximately 5.4 × 10 44 s). In this case, studying any physical process in this interval is believed to be incorrect.
This ambiguity in a sense can be eliminated by using the Gerasimov–Caputo fractional derivatives (fractional operators defined on an infinite interval [18,19]) in the Euler–Lagrange equation, and we offer this in the present article.
In the present paper, we derive the fractional Euler–Lagrange equation and the general equation of motion with the composition of the left and right fractional derivatives defined on infinite intervals, and construct an explicit representation of solutions to a model fractional oscillator equation containing the left and right Gerasimov–Caputo fractional derivatives with origins at plus and minus infinity.
The rest of the paper is structured as follows: Section 2 contains the definitions of fractional integrals and derivatives; in Section 3, we derive the Euler–Lagrange equation and the equation of motion with fractional derivatives defined on infinite intervals; and, in Section 4, we construct solutions of a model one-dimensional equation for the fractional oscillator; Section 5 presents the main findings of this work.

2. Fractional Differentiation

The fractional derivatives of order α ( 0 < α < 1 ) with respect to t, having origins at t = a ( a ), in the Riemann–Liouville and Gerasimov–Caputo senses, are given by ([18,19,20,21])
D a t α g ( t ) = sign ( t a ) d d t D a t α 1 g ( t ) and a t α g ( t ) = sign ( t a ) D a t α 1 d d t g ( t ) ,
respectively. Here, D a t α 1 stands for the Riemann–Liouville fractional integral:
D a t α 1 g ( t ) = sign ( t a ) a t g ( η ) | t η | α Γ ( 1 α ) d η .

3. Euler–Lagrange Equation

Introduce a velocity vector for a material point motion in a certain medium:
v ( t ) = 1 τ Γ ( 1 α ) t ( t η ) α d d η r ( η ) d η = 1 τ t α r ( t ) ( 0 < α 1 ) ,
where r ( t ) is the radius vector, t is the dimensionless (related to τ ) time and τ is some characteristic time for the process. On account of Equation (1), the generalized action functional can be written as
S = L [ t , r ( t ) , v ( t ) ] d t .
In Equation (2) infinite limits for integration are established. To determine the velocity vector of Equation (1), it is necessary to take into account the entire time axis. Assume r 0 (t) is a solution for the extremum problem. Then by Equation (2) and the conditions
r ( t ) = r 0 ( t ) + ε h ( t ) ,
h ( ) = h ( ) = 0
( ε is a small parameter) we get
S = L t , r 0 ( t ) + ε h ( t ) , 1 τ t α r 0 ( t ) + ε τ t α h ( t ) d t .
The extremum of the functional satisfies the condition
d S d ε = d L d ε d t = h ( t ) r + 1 τ D t α v L [ t , r ( t ) , v ( t ) ] d t = 0 .
Using Equation (3) we obtain the Euler–Lagrange equation
r + 1 τ D t α v L [ t , r ( t ) , v ( t ) ] = 0 .
It should be noted that the same problem for classical fields was first considered in [8]. In our case, Equation (4) defines the action function for a mechanical system. Moreover, the other techniques can be applied to get Equation (4), for example, using the results of [22], obtained by exploiting the Riesz fractional derivative.
Next, assume that the space of the considered motion is homogeneous and isotropic. Then the scalar function L should not be explicitly dependent on time. The function L should undergo power series expansion that involves only even powers v ( t ) . Up to the second term of the expansion, we have
L = m v 2 ( t ) 2 U [ r ( t ) , t ] .
Expansion coefficients m and U [ r ( t ) , t ] are independent of the derivative fractional order α . Therefore, there is no indication that these coefficients are different from the mass and potential energy of an object. Substituting Equation (5) into Equation (4), we obtain the equation of motion
m τ 2 D t α t α r ( t ) = F [ r ( t ) , t ] ,
where F [ r ( t ) , t ] = ( / r ) U [ r ( t ) , t ] . The right-hand side of Equation (6) does not explicitly depend on α , and, consequently, F [ r ( t ) , t ] determines the conservative force acting on a material point. The parameter α , according to [23,24] takes into account energy dissipation. Therefore, the equation of motion in the form of Equation (6) is apparently obtained for the first time.
As α = 1 , by Equation (6) we have the classical equation of motion
m τ 2 d 2 r ( t ) d t 2 = F [ r ( t ) , t ] .

4. Solution of Fractional Oscillator Equation

Here, inspired by Equation (6), we are going to solve the equation
D t α t α x ( t ) + λ x ( t ) = f ( t ) ( λ C ) ,
which is a model one-dimensional equation for the fractional oscillator.
Definition 1.
By a solution of Equation (7), we mean a function x ( t ) such that x ( t ) L 1 ( R ) , D t α 1 D t α 1 x ( t ) C 1 ( R ) and x ( t ) satisfies Equation (7).

4.1. Riesz Potentials

Now, we will rewrite Equation (7) in terms of the Riesz potential.
Definition 2.
Let μ ( 0 , 1 ) . The operator K μ is called the Riesz potential, and, for an integrable function x ( t ) , it is defined by
K μ x ( t ) = K μ x ( t ) = c μ | t s | μ 1 x ( s ) d s , c μ = Γ ( 1 μ ) Γ μ 2 Γ 1 μ 2 .
Lemma 1.
Let μ ( 0 , 1 / 2 ) . If x ( t ) L 1 ( R ) , then
D t μ D t μ x ( t ) = K 2 μ x ( t ) .
Proof. 
Let M = max { s , t } and m = min { s , t } . It is easy to check that
Γ ( μ ) 2 D t μ D t μ x ( t ) = t ( ξ t ) μ 1 ξ ( ξ s ) μ 1 x ( s ) d s d ξ
= x ( s ) m ( M ξ ) μ 1 ( m ξ ) μ 1 x ( s ) d ξ d s
= x ( s ) 0 ( M m + η ) μ 1 η μ 1 d η d s
= x ( s ) ( M m ) 2 μ 1 d s · 0 ( 1 + z ) μ 1 z μ 1 d z .
Combining this with the relations
M m = | t s | and 0 ( 1 + z ) μ 1 z μ 1 d z = Γ ( μ ) Γ ( 1 2 μ ) Γ ( 1 μ ) ,
we get Equation (8). □
By Equation (8) we have
D t α t α x ( t ) = d d t D t α 1 D t α 1 d d t x ( t ) = d d t K 2 2 α d d t x ( t ) ,
and, consequently, Equation (7) can be rewritten in the form
d d t K 2 2 α d d t x ( t ) + λ x ( t ) = f ( t ) .

4.2. Special Functions

Consider the functions
C α , λ ( t ) = 1 π 0 cos ξ t ξ 2 α + λ d ξ
and
S α , λ μ ( t ) = 1 π 0 ξ μ sin ξ t ξ 2 α + λ d ξ
for t R , λ C \ { 0 } , α ( 1 / 2 , 1 ) and μ ( 2 , 2 α ) .
Lemma 2.
Let λ C \ { 0 } and α ( 1 / 2 , 1 ) . Then
d d t C α , λ ( t ) = S α , λ 1 ( t ) , d d t S α , λ 1 ( t ) = C α , λ ( t ) ,
and
K 2 2 α S α , λ 1 ( t ) = sign t 2 λ S α , λ 1 ( t ) .
Moreover, if | arg λ | < ( 1 α ) π ε for some ε > 0 , then
lim | t | t 2 C α , λ ( t ) < a n d lim | t | t 3 S α , λ 1 ( t ) < .
Proof. 
The formulas in Equation (9) are derived by direct differentiation. In the same way, an easy computation gives
K 2 2 α S α , λ 1 ( t ) = S α , λ 2 α 1 ( t )
= 1 π 0 ξ 2 α 1 sin ξ t ξ 2 α + λ d ξ = 1 π 0 sin ξ t ξ d ξ λ S α , λ 1 ( t ) .
This proves Equation (10).
Using the properties of the Mittag–Leffler function (see e.g., [25]), we can write
C α , λ ( t ) = 1 π 0 cos ξ t 0 e s ξ s 2 α 1 E 2 α , 2 α λ s 2 α d s d ξ
= 1 π 0 s 2 α 1 E 2 α , 2 α λ s 2 α 0 e s ξ cos ξ t d ξ d s = 1 π 0 E 2 α , 2 α λ s 2 α s 2 α s 2 + t 2 d s .
This gives Equation (11). □

4.3. Solution Representation

Theorem 1.
Let α ( 1 / 2 , 1 ) , | λ | > 0 , and | arg λ | < ( 1 α ) π ε for some ε > 0 . If f ( t ) C ( R ) L 1 ( R ) , then the function
x ( t ) = f ( s ) C α , λ ( t s ) d s
is a regular solution of Equation (7).
Proof. 
Lemmas 1 and 2 allow us to write
D t α t α x ( t ) = d d t K 2 2 α f ( s ) S α , λ 1 ( t s ) d s
= d d t f ( s ) sign ( t s ) 2 λ S α , λ 1 ( t s ) d s = f ( t ) λ f ( s ) C α , λ ( t s ) d s .

5. Conclusions

The principle of least action is often used to describe dissipative systems, in particular in the case of linear oscillators with viscous friction [26]. In this paper, we used this principle in the context of fractional calculus. In particular, we considered the concept of effective velocity for a moving material point, which was defined in terms of the Gerasimov–Caputo fractional derivative. This allowed us to take into account the heredity and nonlocality in dynamic processes of a material point with respect to the time. Using the fractional Euler–Lagrange equation, we derived the general equation of motion with Gerasimov–Caputo fractional derivatives. The equation of motion of such a form (see Equation (6)), has apparently been obtained for the first time in this paper.
In addition, we expressed the composition of the left and right fractional derivatives with origins at plus and minus infinities in terms of the Riesz potential, and constructed an explicit representation of solutions to a model equation of motion for a fractional oscillator.
It is of interest to give a detailed physical interpretation for Equations (4) and (6), and apply these equations to real physical processes. Undoubtedly, it is relevant to develop algorithms for the numerical solution of equations in the form of Equation (7). These issues should be the subjects of separate works.

Author Contributions

Conceptualization, L.E., A.P., and S.R.; methodology, L.E., A.P., and S.R.; validation, L.E., A.P., and S.R.; formal analysis, L.E., A.P., and S.R.; investigation, L.E., A.P., and S.R.; writing–original draft preparation, L.E. (Section 1, Section 2, Section 3, Section 4 and Section 5), A.P. (Section 1, Section 2, Section 4, Section 5), and S.R. (Section 1, Section 2 and Section 3, Section 5); writing–review and editing, L.E., A.P., and S.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tarasov, V.E. Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  2. Uchaikin, V.V. Fractional Derivatives for Physicists and Engineers; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  3. Herrmann, R. Fractional Calculus—An Introduction for Physicists; World Scientific: Singapore, 2018. [Google Scholar]
  4. Riewe, F. Non-conservative Lagrangian and Hamiltonian mechanics. Phys. Rev. E 1996, 53, 1890–1899. [Google Scholar] [CrossRef] [PubMed]
  5. Riewe, F. Mechanics with fractional derivatives. Phys. Rev. E 1997, 55, 3581–3592. [Google Scholar] [CrossRef]
  6. Agrawal, O.P. Formulation of Euler-Lagrange equations for fractionalvariational problems. J. Math. Anal. Appl. 2002, 72, 368–379. [Google Scholar] [CrossRef] [Green Version]
  7. Rekhviashvili, S.S. The Lagrange formalism with fractional derivatives in problems of mechanics. Tech. Phys. Lett. 2004, 30, 33–37. [Google Scholar] [CrossRef]
  8. Baleanu, D.; Muslih, S.I. Lagrangian formulation of classical fields within Riemann-Liouville fractional derivatives. Phys. Scr. 2005, 72, 119–121. [Google Scholar] [CrossRef] [Green Version]
  9. Muslih, S.I.; Baleanu, D. Hamiltonian formulation of systems withlinear velocities within Riemann-Liouville fractional derivatives. J. Math. Anal. Appl. 2005, 304, 599–606. [Google Scholar] [CrossRef] [Green Version]
  10. Rabei, E.M.; Nawafleh, K.I.; Hijjawi, R.S.; Muslih, S.I.; Baleanu, D. The Hamilton formalism with fractional derivatives. J. Math. Anal. Appl. 2007, 327, 891–897. [Google Scholar] [CrossRef]
  11. Atanacković, T.M.; Konjik, S.; Pilipovic, S. Variational Problems with Fractional Derivatives: Euler-Lagrange Equations. J. Phys. Math. Theor. 2008, 41, 095201. [Google Scholar] [CrossRef]
  12. Atanacković, T.M.; Stanković, B. On a differential equation with left and right fractional derivatives. Fract. Calc. Appl. Anal. 2007, 10, 139–150. [Google Scholar] [CrossRef]
  13. Atanacković, T.M.; Stanković, B. On a class of differential equations with left and right fractional derivatives. ZAMM Z. Angew. Math. Mech. 2007, 87, 537–564. [Google Scholar] [CrossRef]
  14. Eneeva, L.M. Boundary value problem for differential equation with fractional order derivatives with different origins. Vestn. KRAUNC. Fiz. Mat. Nauk. 2015, 2, 39–44. [Google Scholar]
  15. Eneeva, L.M. On Neumann problem for equation with fractional derivatives with different starting points. Vestn. KRAUNC. Fiz. Mat. Nauk. 2018, 4, 61–65. [Google Scholar]
  16. Eneeva, L.M. Lyapunov inequality for an equation with fractional derivatives with different origins. Vestn. KRAUNC. Fiz. Mat. Nauk. 2019, 3, 32–39. [Google Scholar]
  17. Eneeva, L.M. A priori estimate for an equation with fractional derivatives with different origins. Vestn. KRAUNC. Fiz. Mat. Nauk. 2019, 4, 41–47. [Google Scholar]
  18. Gerasimov, A.N. A generalization of linear laws of deformation and its application to internal friction problem. Prikl. Mat. Mekh. 1948, 12, 251–260. [Google Scholar]
  19. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Taylor & Francis Books: Didcot, UK, 1993. [Google Scholar]
  20. Nakhushev, A.M. Fractional Calculus and Its Applications; Fizmatlit: Moscow, Russia, 2003. [Google Scholar]
  21. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; North-Holland Math. Stud.; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  22. Agrawal, O.P. Fractional variational calculus in terms of Riesz fractional derivatives. J. Phys. A: Math. Theor. 2007, 40, 6287–6303. [Google Scholar] [CrossRef]
  23. Rekhviashvili, S.S.; Pskhu, A.V.; Agarwal, P.; Jain, S. Application of the fractional oscillator model to describe damped vibrations. Turk. J. Phys. 2019, 43, 236–242. [Google Scholar] [CrossRef]
  24. Rekhviashvili, S.S.; Pskhu, A.V. New Method for Describing Damped Vibrations of a Beam with a Built-in End. Tech. Phys. 2019, 64, 1237–1241. [Google Scholar] [CrossRef]
  25. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer: Berlin/Heidelberg, Germany, 2014. [Google Scholar]
  26. Shalashov, A.G. Can we refer to Hamilton equations for an oscillator with friction? Phys. Usp. 2018, 61, 1082–1088. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Eneeva, L.; Pskhu, A.; Rekhviashvili, S. Ordinary Differential Equation with Left and Right Fractional Derivatives and Modeling of Oscillatory Systems. Mathematics 2020, 8, 2122. https://doi.org/10.3390/math8122122

AMA Style

Eneeva L, Pskhu A, Rekhviashvili S. Ordinary Differential Equation with Left and Right Fractional Derivatives and Modeling of Oscillatory Systems. Mathematics. 2020; 8(12):2122. https://doi.org/10.3390/math8122122

Chicago/Turabian Style

Eneeva, Liana, Arsen Pskhu, and Sergo Rekhviashvili. 2020. "Ordinary Differential Equation with Left and Right Fractional Derivatives and Modeling of Oscillatory Systems" Mathematics 8, no. 12: 2122. https://doi.org/10.3390/math8122122

APA Style

Eneeva, L., Pskhu, A., & Rekhviashvili, S. (2020). Ordinary Differential Equation with Left and Right Fractional Derivatives and Modeling of Oscillatory Systems. Mathematics, 8(12), 2122. https://doi.org/10.3390/math8122122

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop