Next Article in Journal
Fixed Point of Interpolative Rus–Reich–Ćirić Contraction Mapping on Rectangular Quasi-Partial b-Metric Space
Next Article in Special Issue
Monotone Iterative Technique for the Periodic Solutions of High-Order Delayed Differential Equations in Abstract Spaces
Previous Article in Journal
Generalized Attracting Horseshoe in the Rössler Attractor
Previous Article in Special Issue
Global Stabilization of a Reaction Wheel Pendulum: A Discrete-Inverse Optimal Formulation Approach via A Control Lyapunov Function
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Note on the Periodic Solutions for a Class of Third Order Differential Equations

1
Department of Mathematics and Computer Science, Larbi Tebessi University, 12002 Tebessa, Algeria
2
Departamento de Matemática Aplicada y Estadística, Universidad Politécnica de Cartagena, Hospital de Marina, 30203 Cartagena, Región de Murcia, Spain
3
Centro Universitario de la Defensa, Academia General del Aire, Universidad Politécnica de Cartagena, 30720 Santiago de la Ribera, Región de Murcia, Spain
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(1), 31; https://doi.org/10.3390/sym13010031
Submission received: 9 November 2020 / Revised: 17 December 2020 / Accepted: 23 December 2020 / Published: 27 December 2020
(This article belongs to the Special Issue Advances in Nonlinear, Discrete, Continuous and Hamiltonian Systems)

Abstract

:
The aim of the present work is to study the necessary and sufficient conditions for the existence of periodic solutions for a class of third order differential equations by using the averaging theory. Moreover, we use the symmetry of the Monodromy matrix to study the stability of these solutions.
MSC:
37C80; 37G15; 37C30

1. Introduction and Statement of the Main Results

Third-order differential equations
x + g ( x , x ˙ , x ¨ ) = 0
are associated with so-called jerk dynamics. In the literature, the third derivative of the position is called the jerk, see [1] and references therein. Many applications of the jerk dynamics to multidisciplinary sciences are been development, see for instance [2].
From a purely mathematical view, various authors have done interesting research about the periodic structure, chaos and other modern topics on nonlinear analysis related to the jerk dynamics; see for example [3,4,5,6].
It is well known that the analytical study of the periodic orbits of differential equations is one of the most important problems in qualitative theory. It is in terms of number, existence, non-existence, uniqueness and stability of such periodic orbits. Motivated by it, it will make sense to study the perturbed jerk dynamics, i.e, the third-order differential equations,
x + g ( x , x ˙ , x ¨ ) = ε f ( t , x , x ˙ , x ¨ , ε ) ,
with ε is a small parameter and f , g are smooth functions.
In what follows, we endeavor to study the periodic structure of this type of problem.
There are different methods to face such a study. In the literature, there are some papers analyzing the periodic orbits for this class of third order differential equations using the Leray-Schauder method; see [7,8]. There are also some works that use the Leray-Schauder method in the field of partial differential equations, for example, studying the existence of weak solutions to the nonlinear elliptic problem [9] and studying the existence of traveling waves for delay reaction-diffusion equations, see [10].
In the present paper, we shall study the periodic orbits of this class of third-order differential equations, by using the averaging theory of dynamical systems, which is a completely different approach to the ones used in the mentioned papers. In [11,12,13] averaging theory was applied to order three systems too.
One of the main problems in the qualitative study of a system of nonlinear differential equations is the investigation of the periodic structure of the problem. Many of the current methods used by researchers who study these problems are born from the ideas and considerations expressed by Henri Poincaré in his treatise [14].
Among the analytical methods considered, the averaging method is a clear candidate to address to the problems of the following type of vectorial differential equations
d x d t + f ( x ) = ε g ( x , t ) .
In the literature on this matter, we can go back to Laplace where he introduces the original foundations of the averaging method [15]. Already in the previous century, the averaging theory was theoretically-based. In numerous works, different theoretical and applied results of the development of the theory were collected.
Please note that in modern language, the averaging theory relies on the so-called displacement map of the system and by mean of this method allows obtaining an asymptotic expansion in the small parameter ε of this map. The zeros of the map which are hyperbolic provide an asymptotic expansion of the initial condition of the periodic orbit. For more details, see [16] and its references.
In [17], Garcia and Llibre studied the periodic structure of second order differential equations of the form
x ¨ + g ( x ) = ε f ( t , x , x ˙ , , ε ) ,
and they chose a perturbation family of a Duffing equation and forced pendulum in the form:
x ¨ + g ( x ) = μ 2 n + 1 p ( t ) + μ 4 n + 1 q ( t , x , x ˙ , μ ) ,
where n N , μ is a small parameter, and
g ( x ) = x + x 2 n + 1 ( b + x h ( x ) ) ,
where h ( x ) is smooth, b 0 , p ( t ) and q ( t , x , y , μ ) are smooth and periodic with period 2 π in the variable t.
Therefore, inspired by this study we decided to study the generalization of this system considering the jerk variable in the form
x + g ( x , x ˙ , x ¨ ) = μ 2 n + 1 p ( t ) + μ 4 n + 1 q ( t , x , x ˙ , x ¨ , μ ) ,
where n N , μ is a small parameter, and
g ( x , x ˙ , x ¨ ) = x ¨ + x ˙ x + x x ˙ 2 n + 1 ( b + x h ( x ) ) ,
and q ( t , x , x ˙ , x ¨ , μ ) , p ( t ) are smooth and periodic with period 2 π in t, h ( x ) is even smooth function with b 0 .
We introduce the following variables y , z = x ˙ , x ¨ , then we write (2) as the differential system
x ˙ = y , y ˙ = z , z ˙ = x y + z x y 2 n + 1 ( b + x h ( x ) ) + μ 2 n + 1 p ( t ) + μ 4 n + 1 q ( t , x , y , z , μ ) .
Let α 1 and β 1 be the first Fourier coefficients of p ( t ) , i.e.,
α 1 = 1 π 0 2 π p ( t ) cos t d t , β 1 = 1 π 0 2 π p ( t ) sin t d t ,
and let Γ ( x ) the gamma function.
In the following definition, we are concerned with the limit cycle and its properties:
Definition 1.
A periodic orbit in the plane is called a limit cycle if it is an isolated trajectory in the set of all periodic orbits of a differential equation.
The limit cycle is topologically distinguished with neighboring trajectories that are not closed.
The statement of our main result on the periodic structure of system (2) is the following theorem:
Theorem 1.
If α 1 β 1 b 0 , α 1 β 1 then for μ 0 sufficiently small the differential Equation (2) has a periodic solution x ( t , μ ) of period 2 π satisfying the initial condition
x ( 0 , μ ) = μ π 1 4 n + 2 Γ ( n + 2 ) 2 b Γ n + 3 2 1 2 n + 1 α 1 2 + β 1 2 n 1 + 2 n α 1 β 1 2 + O μ 2 n + 1 .

2. Some Results on the Other Averaging Theory

We consider the differential system of the form
x ˙ = F 0 ( t , x ) + ε F 1 ( t , x ) + ε 2 R ( t , x , ε ) ,
with ε 0 is small parameter, F 0 , F 1 : R × D R n and R : R × D × ε 0 , ε 0 R n are C 2 functions, 2 π –periodic in the first variable, and D is an open subset of R n . We suppose that the differential system
x ˙ = F 0 ( t , x ) ,
has a submanifold of periodic solutions. The system (5) is called the unperturbed system.
In what follows, we denote by x ( t , z ) the solution of the unperturbed system (5) satisfying x ( 0 , z ) = z .
We consider the differential system
y ˙ = D x F 0 ( t , x ( t , z ) ) y ,
this system represents the first order variational equation on the periodic solution x ( t , z ) .
Let M z ( t ) be the fundamental matrix of the linear differential system (6) such that M z ( 0 ) be the identity matrix and let the projection ξ : R k × R n k R k defined by ξ x 1 , , x n = x 1 , , x k .
In what follows, we denote by C l ( V ) the closure of V. The closure of a subset V is the smallest closed subset containing V.
Theorem 2.
Assume that V R k be an open and bounded subset and β 0 : C l ( V ) R n k be a C 2 function. We consider the following statements:
(i) 
Z = z α = α , β 0 ( α ) : α C l ( V ) D and that for every z α Z the solution x t , z α of (5) is periodic of period T.
(ii) 
for every z α Z there exists a fundamental matrix M z α ( t ) of (6) such that the matrix M z α 1 ( 0 ) M z α 1 ( T ) has in the upper right corner the k × ( n k ) zero matrix, and in the lower right corner a ( n k ) × ( n k ) matrix Δ α with d e t Δ α 0 .
Let us define the function F : C l ( V ) R k as
F ( α ) = ξ 0 T M z α 1 ( t ) F 1 t , x t , z α d t .
If there exists a V with F ( a ) = 0 and
d e t ( ( d F / d α ) ( a ) ) 0 ,
then for ε 0 sufficiently small there exists a periodic solution x ( t , ε ) of period 2 π for differential system (4) such that x ( 0 , ε ) z a as ε 0 .
For a proof on the existence of the periodic solution of Theorem 2 see Malkin [18] and Roseau [19], or [20] for a proof in a shorter way. By using an idea similar to the one used to proved the stability of the periodic solutions of Theorem 11.6 in [21] we can prove the stability of periodic orbits.

3. Proof of Theorem 1

To apply the averaging theory summarized at Section 2 to the differential system (3) we do the following rescaling x μ x , y μ y and z μ z . Hence the system (3) writes as follows
x ˙ = y , y ˙ = z , z ˙ = x y + z + μ 2 n b x y 2 n + 1 + p ( t ) μ 2 n + 1 x x y 2 n + 1 h ( μ x ) + μ 4 n q ( t , μ x , μ y , μ z , μ ) ,
Suppose that μ 2 n = ε , O ε = μ 2 n + 1 x x y 2 n + 1 h ( μ x ) + μ 4 n q ( t , μ x , μ y , μ z , μ ) and F t , x , y , z = b x y 2 n + 1 + p ( t ) the system (9) becomes
x ˙ = y , y ˙ = z , z ˙ = x y + z + ε b x y 2 n + 1 + p ( t ) + O ε 2 , = x y + z + ε F t , x , y , z + O ε 2 .
The origin is the unique equilibrium point of the unperturbed system (10) ε = 0 , the eigenvalues of this unperturbed system at the origin are i , i , 1 . We do the following linear change
X , Y , Z T = C ( x , y , z ) T ,
where the matrix
C = 1 1 0 0 1 1 1 0 1 .
System (10) is transformed into the next system by using the previous linear change
X = Y , Y = X + ε F ˜ ( t , X , Y , Z ) + O ε 2 , Z = Z + ε F ˜ ( t , X , Y , Z ) + O ε 2 ,
where the linear part of this system in the real Jordan canonical normal form and
F ˜ ( t , X , Y , Z ) = F t , X Y + Z 2 , X Y + Z 2 , X + Y + Z 2 .
Returning to (4), we consider x = ( X , Y , Z ) , F 0 ( t , x ) = ( Y , X , Z ) , F 1 ( t , x ) = ( 0 , F ˜ , F ˜ ) and F 2 ( t , x ) = 0 , O ε 2 , O ε 2 . Let x t , X 0 , Y 0 , Z 0 , ε be the solution to system (11) with x 0 , X 0 , Y 0 , Z 0 , ε = X 0 , Y 0 , Z 0 . Please note that the unperturbed system (11) ε = 0 has a linear center at the origin in the ( X , Y ) plane, and the periodic solution x t , X 0 , Y 0 , 0 , 0 = ( X ( t ) , Y ( t ) , Z ( t ) ) is
X ( t ) = X 0 cos t Y 0 sin t , Y ( t ) = Y 0 cos t + X 0 sin t , Z ( t ) = 0 .
Please note that ( X ( t ) , Y ( t ) , Z ( t ) ) has period 2 π . For our system, V and α from Theorem 2 are
V = ( X , Y , 0 ) : 0 < X 2 + Y 2 < ρ ,
where ρ > 0 is arbitrary and α = X 0 , Y 0 V .
The fundamental matrix M z ( t ) = M ( t ) of the first order variational equation (6) satisfying the system (11) ε = 0 is
M ( t ) = cos t sin t 0 sin t cos t 0 0 0 e t ,
We note that M ( t ) is independent from the initial condition ( X 0 , Y 0 , 0 ) , and
M 1 ( 0 ) M 1 ( 2 π ) = 0 0 0 0 0 0 0 0 1 e 2 π .
Very briefly, we confirmed that the statements of Theorem 2 hold. Now we study the solutions α = X 0 , Y 0 V of the function F α given in (7). In other words, we have F α = ( F 1 ( α ) , F 2 ( α ) ) where
F ( α ) = 0 2 π M 1 ( t ) F 1 ( t , x ( t , z α , 0 ) ) d t = F 1 X 0 , Y 0 F 2 X 0 , Y 0 = 0 2 π sin t F ˜ t , x t , X 0 , Y 0 , 0 , 0 d t 0 2 π cos t F ˜ t , x t , X 0 , Y 0 , 0 , 0 d t = 0 2 π sin t F t , X ( t ) Y ( t ) 2 , X ( t ) + Y ( t ) 2 , X ( t ) + Y ( t ) 2 d t 0 2 π cos t F t , X ( t ) Y ( t ) 2 , X ( t ) + Y ( t ) 2 , X ( t ) + Y ( t ) 2 d t = 0 2 π sin t b X 0 cos t Y 0 sin t 2 n + 1 + p ( t ) d t 0 2 π cos t b X 0 cos t Y 0 sin t 2 n + 1 + p ( t ) d t ,
Using simple calculations, it can be shown that
0 2 π sin t X 0 cos t Y 0 sin t 2 n + 1 d t = 2 π Γ 3 2 + n Γ ( 2 + n ) Y 0 X 0 2 + Y 0 2 n , 0 2 π cos t X 0 cos t Y 0 sin t 2 n + 1 d t = 2 π Γ 3 2 + n Γ ( 2 + n ) X 0 X 0 2 + Y 0 2 n .
We solve the following system
F 1 X 0 , Y 0 F 2 X 0 , Y 0 = b 2 π Γ 3 2 + n Γ ( 2 + n ) Y 0 X 0 2 + Y 0 2 n + π β 1 b 2 π Γ 3 2 + n Γ ( 2 + n ) X 0 X 0 2 + Y 0 2 n + π α 1 = 0 0 .
We obtain a unique solution
X 0 * Y 0 * = π 1 4 n + 2 Γ ( n + 2 ) 2 b Γ n + 3 2 1 2 n + 1 α 1 2 + β 1 2 n 1 + 2 n α 1 β 1 .
Returning to (8), the d e t D F X 0 * , Y 0 * is
1 + 2 n π 4 n + 1 2 n + 1 Γ ( n + 2 ) 2 b Γ n + 3 2 2 1 + 2 n α 1 2 + β 1 2 2 n 1 + 2 n ,
in order because α 1 β 1 b 0 , the statements of Theorem 2 hold. Using the following relationship
x ( t ) = X ( t ) Y ( t ) + Z ( t ) 2
where X ( t ) , Y ( t ) , Z ( t ) ) are given by (12), we can prove that
x ( 0 , μ ) = μ X 0 * Y 0 * 2 + O μ 2 n + 1 = μ π 1 4 n + 2 Γ ( n + 2 ) 2 b Γ n + 3 2 1 2 n + 1 α 1 2 + β 1 2 n 1 + 2 n α 1 β 1 2 + O μ 2 n + 1 .
This completes the proof of Theorem 1. □
Remark 1.
The eigenvalues of the Jacobian matrix D F X 0 * , Y 0 * are
λ 1 , 2 = ± i 2 n + 1 1 2 π 4 n + 1 4 n + 2 2 b Γ 3 2 + n Γ ( 2 + n ) 4 n 2 + 2 n + 1 2 n + 1 α 1 2 + β 1 2 n 1 + 2 n .
In consequence, we cannot use the averaging theory found in Section 2 to study the stability of the periodic solution of the differential system (10). To solve this problem we use the symmetry of the Monodromy matrix.
The variational equation of this periodic orbit is a given by
y · = A ( t ) y ,
where y = x , y , z R 3 and A ( t ) a 2 π –periodic matrix given by
A ( t ) = A 0 + ε A 1 ( t ) .
with
A 0 = 0 1 0 0 0 1 1 1 1 , A 1 ( t ) = 0 0 0 0 0 0 b ( 2 n + 1 ) ( x ( t ) y ( t ) ) 2 n b ( 2 n + 1 ) ( x ( t ) y ( t ) ) 2 n 0
The multipliers of the monodromy matrix of the variational equation are 1 + O ( ε 2 ) , 1 + O ( ε 2 ) and e 2 π + O ( ε 2 ) . From it, we conclude that the periodic solution of the differential system (10) is unstable because the multiplier e 2 π + O ( ε 2 ) > 1 .

4. Conclusions

Appropriately introducing the small parameter enables us to apply the averaging theory to prove analytically the existence of periodic solutions for a class of third-order differential equations. The application of this theory to the third-order differential equations gave interesting results on the stability of the periodic orbits. Our future work will consist of applying this theory to differential equations of order n, with n 4 with particular physical applications.

Author Contributions

Investigation, Z.D., J.L.G.G. and J.A.V. All authors have read and agreed to the published version of the manuscript.

Funding

Foundation Seneca: 20783/PI/18 Ministry of Science, Innovation and Universities: PGC2018-097198-B-I00.

Data Availability Statement

This paper has not data involved.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schot, S.H. Jerk: The time rate of change of acceleration. Am. J. Phys. 1978, 46, 1090. [Google Scholar] [CrossRef]
  2. Zeithamer, T. Methodology of Theoretical Physics in Economics: Jerk Dynamics. Procedia Econ. Financ. 2015, 30, 967–973. [Google Scholar] [CrossRef] [Green Version]
  3. Gottlieb, H.P.W. Simple nonlinear jerk functions with periodic solutions. Am. Phys. 1998, 66, 903–906. [Google Scholar] [CrossRef]
  4. Patidar, V.; Sud, K.K. Bifurcation and chaos in simple jerk dynamical systems. Pramana 2005, 64, 75–93. [Google Scholar] [CrossRef]
  5. Sprott, J.C. Some simple chaotic jerk functions. Am. J. Phys. 1997, 65, 537–543. [Google Scholar] [CrossRef]
  6. Spprott, J.C. Elegant Chaos. Algebraically Simple Chaotic Flows; World Scientific: Singapore, 2010. [Google Scholar]
  7. Reissig, R. Periodic solutions of a third order nonlinear differential equation. Ann. Mat. Pura Appl. 1972, 92, 193–198. [Google Scholar] [CrossRef]
  8. Xiankai, H. On the existence of periodic solutions for the third-order nonlinear ordinary differential equations. Appl. Math.-A J. Chin. Univ. 1999, 34, 125–130. [Google Scholar] [CrossRef]
  9. Jung, T.; Choi, Q.H. Existence of four solutions for the elliptic problem with nonlinearity crossing one eigenvalue. J. Inequal. Appl. 2013, 2013, 187. [Google Scholar] [CrossRef] [Green Version]
  10. Volpert, V. Existence of waves for a bistable reaction-diffusion system with delay. J. Dyn. Differ. Equ. 2020, 32, 615–629. [Google Scholar] [CrossRef]
  11. Llibre, J.; Roberto, L. On the periodic orbits of the third-order differential equation x′′′μx′′ + xμx = εF(x,x,x′′). Appl. Math. Lett. 2014, 26, 425–430. [Google Scholar] [CrossRef] [Green Version]
  12. Llibre, J.; Yu, J.; Zhang, X. Limit Cycles for a Class of Third-Order Differential Equations. Rocky Mount. J. Math. 2010, 40, 581–594. [Google Scholar] [CrossRef]
  13. Toniol Cardin, P.; Llibre, J. Transcritical and zero-Hopf bifurcations in the Genesio system. Nonlinear Dyn. 2017, 88, 547–553. [Google Scholar] [CrossRef] [Green Version]
  14. Poincaré, H. Les m éthodes Nouvelles de la Mécanique Céleste; Gauthier-Villars: Paris, France, 1989; Volume 3. [Google Scholar]
  15. Marquis de Laplace, P.S. Traité de Mécanique Céleste; Bachelier, 1829; Volume 1. [Google Scholar]
  16. Llibre, J.; Novaes, D.D.; Teixeira, M.A. Higher order averaging theory for finding periodic solutions via Brouwer degree. Nonlinearity 2014, 27, 563–583. [Google Scholar] [CrossRef] [Green Version]
  17. Garcia, A.; Llibre, J. Existence of Periodic Solutions for a Class of Second Order Ordinary Differential Equations. Acta Appl. Math. Vol. 2020, 169, 193–197. [Google Scholar] [CrossRef]
  18. Malkin, I.G. Some Problems of the Theory of Nonlinear Oscillations; State Pub. House of Technical and Theoretical Literature Moscow: Moscow, Russia, 1956. (In Russian) [Google Scholar]
  19. Roseau, M. Vibrations Non Linéaires et Théorie de la Stabilité; Springer Tracts in Natural Philosophy; Springer: New York, NY, USA, 1976; Volume 8. [Google Scholar]
  20. Buică, A.; Françoise, J.P.; Llibre, J. Periodic solutions of nonlinear periodic differential systems with a small parameter. Commun. Pure Appl. Anal. 2006, 6, 103–111. [Google Scholar] [CrossRef]
  21. Verhulst, F. Nonlinear Differential Equations and Dynamical Systems, Universitext; Springer: Berlin, Germany, 1996. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Diab, Z.; Guirao, J.L.G.; Vera, J.A. A Note on the Periodic Solutions for a Class of Third Order Differential Equations. Symmetry 2021, 13, 31. https://doi.org/10.3390/sym13010031

AMA Style

Diab Z, Guirao JLG, Vera JA. A Note on the Periodic Solutions for a Class of Third Order Differential Equations. Symmetry. 2021; 13(1):31. https://doi.org/10.3390/sym13010031

Chicago/Turabian Style

Diab, Zouhair, Juan L. G. Guirao, and Juan A. Vera. 2021. "A Note on the Periodic Solutions for a Class of Third Order Differential Equations" Symmetry 13, no. 1: 31. https://doi.org/10.3390/sym13010031

APA Style

Diab, Z., Guirao, J. L. G., & Vera, J. A. (2021). A Note on the Periodic Solutions for a Class of Third Order Differential Equations. Symmetry, 13(1), 31. https://doi.org/10.3390/sym13010031

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