![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
Abstract
The use of improved generalized Riccati equation mapping method has been demonstrated to find some new exact travelling wave solutions to space-time fractional non-liner double dispersive equation (DDE). The equation is used for modelling wave propagation in elastic inhomogeneous Murnaghan's rod. We have used Caputo's fractional derivative to achieve the fractional version of Murnaghan's rod equation. Improved generalized Riccati equation mapping method proves to be very effective tool to find a variety of soliton solutions. As a result, we found many new and general solutions including dark, combined dark-bright, singular periodic wave, combined singular periodic wave solutions and rational solutions. We have simulated the solitons, to check their types, with the help of graphs and all the solutions obtained in this article have been verified by back substitution in original equation by using Maple 17.
1. Introduction
The use of fractional calculus to model certain real-life phenomena is getting a great attention nowadays. Nonlinear fractional differential equations (NLFDEs) appear as a direct result of this attention. Nonlinear fractional partial differential equations (NLFPDEs) cover a major share of those NLFDEs and they are used to model such phenomena where the dependent variable is reliant on more than one independent variable. NLFPDEs are generalization of nonlinear partial differential equations (NPDEs) in which the orders of derivatives involved are fractional. These equations have numerous applications in different fields of engineering and physical sciences such as in fluid mechanic, fractional dynamics, wave propagation and so on [Citation1]. It is very important not only to formulate the governing nonlinear PDE of a certain phenomenon but also to find out its exact solutions. Solutions of an equation, governing a certain real-life phenomenon, give us very useful details of the phenomenon itself and can be used to understand and predict the variations in the depended variable (and the quantities driven by it). Many scientist and researchers have extensively studied nonlinear PDEs and their methods to construct exact solutions. Some of the important works are done in [Citation2–13].
In this study, we are interested in a special type of exact solutions of NLFPDEs known as solitary wave solutions. Since solitons have been proved to be the exact solutions of a large class of NLPDEs, their complete understanding would lead us to a broad understanding of the real-life phenomena themselves. Some of the methods that are already being used to find solution of fractional order nonlinear partial differential equations are Homotopy perturbation method (HPM) [Citation14], Lie algebra method [Citation15,Citation16], Variational iteration method (VIM) [Citation17,Citation18], F-expansion method [Citation19], Exp-function method [Citation20,Citation21], Fan sub-equation method [Citation22], -expansion method [Citation23], Improved tan (
)-expansion method [Citation24], Exp
method [Citation25], Kudryashov method [Citation26], etc. Some of these methods provide exact solutions to NLFPDEs (like Exp-function method, Fan sub-equation method,
-expansion method, etc.) while the others provide series solution (like VIM and HPM). Nowadays mathematicians are trying to extend conventional methods to make them capable of solving fractional order partial differential equations. These extended methods would enable scientist working on fractional models to deal with them more effectively. Finding exact solutions of NLFPDEs used to be a herculean task, however, with modern symbolic computation tools have make the task relatively easier. In a result of these computational tools, the efforts to extend the methods used to solve integer order NLPDEs to their fractional counterparts, and apply them to solve real-life fractional models, have gain a tremendous popularity.
Improved generalized Riccati equation method is one of the methods to the get exact traveling wave solutions to the PDEs having both steepening and spreading effects. It is a straight-forward and easy-to-use method that, by symbolic computation, can generate many different types of exact travelling wave solutions. Zhu [Citation27] introduced this method with the extended tanh-function method to solve (2 + 1) dimensional Boiti–Leon–Pempinelle equation. Bekir and Cevikel [Citation28] used Riccati equation combined with tanh–coth method to solve nonlinear coupled equation in mathematical physics. Li et al. [Citation29] used this method to find exact solutions of (3 + 1)-dimensional Jimbo–Miwa equation. Tala-Tebue et al. [Citation30] used this method to solve discrete nonlinear electrical transmission lines in (2 + 1) dimension. Salathiel et al. [Citation31] utilized generalized Riccati equation mapping method to construct soliton and travelling wave solutions for discrete electrical lattice. Koonprasert et al. [Citation32] implemented this method to find more explicit solitary solutions to the space–time fractional fifth-order nonlinear Sawada–Kotera equation. Most recently, Bibi et al. [Citation33] have used this method on Caudrey–Dodd–Gibsson equation. Their work shows that the improved generalized Riccati equation method has a great protentional for solving fractional order partial differential equations. In this context, it seems appropriate to apply this method to some other real-life fractional models and obtain their unexplored exact travelling wave solutions.
The doubly dispersive equation (DDE) is an important nonlinear physical model describing the nonlinear wave propagation in the elastic inhomogeneous circular cylinder Murnaghan’s rod. It is an important model to study the wave propagation in nonlinear elastic solids. Solitary strain waves have great importance in the study of seismology, acoustics, introscopy, examination of sudden destruction, long distance energy transfer, vibro-impact treatments of hard materials also in studies for the development of non-destructive testing techniques especially for pipelines, and to understand the physical properties and internal structure of solids like brass, steel, glass and polymers. This relatively new area of study might increase interest in researchers to study nonlinear physics of solids, nonlinear mechanics and condensed matter physics [Citation34]. The global existence and blow-up of solutions for doubly dispersive equation were discussed by Harby et al. [Citation35]. Cattani et al. [Citation36] had used extended Sinh-Gordon equation expansion method (ShGEEM) and the modified exp(−ϕ(ζ))-expansion function method, to find different types of solitary wave solutions. Moreover, Baskonus et al. [Citation37] solved inhomogeneous Murnaghan’s rod by F-expansion method and obtained Jacobi elliptic function solutions including bright and dark solitons, topological, non-topological, singular, periodic, their combinations and compound solitons.
From the above references, it is clear that the solutions to DD equation (even of integer order) are a very recent topic and it still needs a lot exploration. Up to the best of our knowledge, no one has ever studied the fractional model for DD equation. To fill this gap, in this paper, we have used generalized Riccati equation method and able to construct abundant exact solutions (including solitary wave solutions, periodic wave solutions and rational solutions) to nonlinear DD equation (with space–time fractional derivatives) which describes wave propagation in nonlinear elastic inhomogeneous Murnaghan’s rod.
This paper is organized as follows: Section 1 includes introduction, Section 2 provides definition and properties of Caputo fractional derivative, in Section 3, improved generalized Riccati equation mapping method has been described briefly, Section 4 includes application of method described in Section 3 to obtain new solitary wave solutions of Murnaghan’s rod equation, Section 5 consists results and discussion on some of obtained results, Section 6 includes conclusion.
2. Definitions and properties of Caputo fractional derivative
Let be a smallest integer that is greater than
, the Caputo time fractional derivative operator of order
of the function
is defined as follows [Citation38]:
(1)
(1)
Important characteristics of Caputo fractional derivative can be seen in [Citation38] and keeping in view the length of the manuscript we see fit not to rewrite them here.
3. Improved generalized Riccati equation mapping method
Let us consider the fractional order differential equation with independent variables and some dependent function
(2)
(2)
where
and
are known as Caputo fractional derivatives of
and
is a polynomial in
with its various orders of nonlinear partial fractional derivatives.
Step 1. Let
(3)
(3)
where
(4)
(4)
is a fractional complex transformation which can convert nonlinear fractional differential equation (2) into nonlinear ordinary differential equation, where
is a constant which is to be determined, this fractional complex transform is an easy transform to convert fractional differential equation into ordinary differential equation. As the classical chain rule does not hold for fractional derivatives, [Citation39] introduced new chain rule in terms of
index and suggested that
(5)
(5)
where
,
are called sigma indexes [Citation40], without loss of generality we can take
Where
is constant. By using the definition of Caputo derivative and applying the chain rule mentioned in Equation (5) along with complex transformation (4) into Equation (2) we get the following nonlinear ordinary differential equation.
(6)
(6)
where
indicates derivative in term of
. We integrate Equation (6) as many times as we get at least one term without derivative.
Step 2. We suppose that the following series expansion is the solution of Equation (6).
(7)
(7)
where
being constants, which are to be determined provided
The function
satisfies the Riccati differential equation
(8)
(8)
Step 3. Positive integer in Equation (7) can be found by using homogeneous balance between the derivatives of highest order and the nonlinear terms in Equation (6).
Step 4. Substituting Equation (7) along with Equation (8) into Equation (6) followed by collecting all the same order terms together, we get the polynomial equation in
and
, where
Equalizing coefficients of the resulting polynomial to zero, we get over-determined system of algebraic equations for
where
Step 5. With the help of Maple, we solve the system described in step 4, and obtain , where
. We substitute these values in Equation (7) coupled with solutions of Equation (8) and applying the transformation in Equation (6) we construct several exact solutions of Equation (2), establishing four families [Citation32].
4. Implementation of proposed method to Murnaghan’s rod
In this section, we apply improved generalized Riccati equation mapping method on space–time fractional double dispersive equation which is given as
(9)
(9)
where
is the strain wave function,
are combinations of the constant scale factors,
is the small parameter,
is the nonlinearity,
is the density,
is the Poisson coefficient, for details see [Citation34–37]. Parameter
is the order of fractional time and space derivatives. Where
and
are the Caputo fractional derivative of
with respect to
and
respectively. Now, consider the following nonlinear fractional order wave transformation:
where
is constant. Substituting the above-mentioned transformation (incorporating Caputo’s fractional derivative of power function) and Equation (5) into Equation (9) we get the following nonlinear ODE:
(10)
(10)
by integrating Equation (10) and using homogeneous balance principle between the highest order derivative and nonlinearity yields
Therefore, Equation (10) has a solution
(11)
(11)
now substituting Equation (11) along with Equation (8) into Equation (10) after collecting all terms with the same order in
and
, where
and equating each coefficient to zero, we obtain a system of nonlinear algebraic equations. Solving these equations by using Maple 17, we get following non-trivial solutions:
Case 1:
(12)
(12)
(13)
(13)
Case 2:
(14)
(14)
(15)
(15)
Case 3:
(16)
(16)
(17)
(17)
Case 4:
(18)
(18)
(19)
(19)
According to the method the improved generalized Riccati equation mapping technique yield 27 solutions of the Riccati equation in Equation (8) incorporating four different families [Citation32]. For the case 1, substituting the values from Equation (12) into Equation (13) along with the Riccati equations solutions, we can get following different types of solutions. Please note, in all the following solutions, these substitutions have been made to make the results more elegant:
Family 1:
When and
or
, the hyperbolic function solutions of Equation (9) are as follows:
(20)
(20)
(21)
(21)
(22)
(22)
(23)
(23)
(24)
(24)
(25)
(25)
(26)
(26)
where two non-zero real constants A and B satisfies
(27)
(27)
(28)
(28)
(29)
(29)
(30)
(30)
(31)
(31)
Family 2:
If and
, we have the following trigonometric solutions:
(32)
(32)
(33)
(33)
(34)
(34)
(35)
(35)
(36)
(36)
(37)
(37)
(38)
(38)
where two non-zero real constants A and B satisfies
(39)
(39)
(40)
(40)
(41)
(41)
(42)
(42)
(43)
(43)
For case 2, we have following solutions.
Family 1:
The hyperbolic function solutions of Equation (9) (when and
or
) are:
(44)
(44)
(45)
(45)
(46)
(46)
(47)
(47)
(48)
(48)
(49)
(49)
(50)
(50)
Where two non-zero real constants and
satisfy
(51)
(51)
(52)
(52)
(53)
(53)
(54)
(54)
(55)
(55)
Family 2:
The trigonometric solutions of Equation (9) (if and
) are:
(56)
(56)
(57)
(57)
(58)
(58)
(58a)
(58a)
(59)
(59)
(60)
(60)
(61)
(61)
where two non-zero real constants
and
satisfies
(62)
(62)
(63)
(63)
(64)
(64)
(65)
(65)
(66)
(66)
For case 3, we have
Family 1:
When and
or
, we have
(67)
(67)
(68)
(68)
(69)
(69)
(70)
(70)
(71)
(71)
(72)
(72)
(73)
(73)
where two non-zero real constants A and B satisfies
(74)
(74)
(75)
(75)
(76)
(76)
(77)
(77)
(78)
(78)
Family 2:
If and
, we have the following solutions:
(79)
(79)
(80)
(80)
(81)
(81)
(82)
(82)
(83)
(83)
(84)
(84)
(85)
(85)
where two non-zero real constants A and B satisfies
(86)
(86)
(87)
(87)
(88)
(88)
(89)
(89)
(90)
(90)
Family 3:
When and
we get soliton-like solutions
(91)
(91)
(92)
(92)
where
is constant.
Family 4:
When and
, we have following rational solution:
(93)
(93)
where C is an arbitrary constant.
In the case 4, we get
Family 1:
When and
or
, the hyperbolic function solutions are as follows:
(94)
(94)
(95)
(95)
(96)
(96)
(97)
(97)
(98)
(98)
(99)
(99)
(100)
(100)
where two non-zero real constants A and B satisfies
(101)
(101)
(102)
(102)
(103)
(103)
(104)
(104)
(105)
(105)
Family 2:
When and
or
, the trigonometric solutions of Equation (9) are
(106)
(106)
(107)
(107)
(108)
(108)
(109)
(109)
(110)
(110)
(111)
(111)
(112)
(112)
where two non-zero real constants A and B satisfies
(113)
(113)
(114)
(114)
(115)
(115)
(116)
(116)
(117)
(117)
Family 3:
When and
the hyperbolic function solutions are
(118)
(118)
(119)
(119)
where
is constant.
5. Results and discussion
In this section, we compare our results with already present results in the literature, obtained by different mathematical techniques. Here we have used improved generalized Riccati equation mapping method on Equation (9) and constructed different types of exact solutions including dark, combined dark-bright, singular periodic wave, combined singular periodic wave solutions and one rational solution. In [Citation36], authors construct topological, non-topological, singular, compound topological-non-topological bell-type and compound singular, soliton-like, singular periodic wave and exponential function solution and in [Citation37] authors calculated singular, periodic, bright, dark, their combinations and compound solitons. We have compared these and found out that our results are different and new from the results obtained by the authors in [Citation36,Citation37].
It is worth mentioning that the solutions obtained in this study represent certain real-life situations. For example, the tan-hyperbolic solutions are useful in calculating the magnetic moment and rapidity of special relativity, cos-hyperbolic solutions represent the shape of hanging cable, cot-hyperbolic solutions appear in the Langevin function which arise in magnetic polarization, sec-hyperbolic solutions represent the laminar jet profile [Citation41]. Similarly, exact solutions with the periodic functions exhibit periodic wave phenomena. It is significant to mention here that a lot of new solutions have been produced for Murnaghan’s rod, and for the first time this equation has been solved for space–time fractional order. The reason of using fractional differential equation is that it is naturally related to physical phenomena with memory. Many well-known equations can be solved by space–time fractional differential equations to get variety of new solutions. For the better understanding of the solitary wave phenomenon graphs of some obtained solutions has been discussed here with the aid of mathematical software Maple 17. Figure represents 3D graphs of which is dark soliton generated by taking fractional order
, with some given parameters
. Figure represents 3D graphs of solutions
with fractional order
, exhibits combined singular periodic wave solution by taking parameters
. Figure represents 3D graphs dark singular solitons of
with fractional order
by choosing parameters
Figure represents 3D graphs of combined dark-bright soliton generated by
with fractional order
, by taking
. Figure represents 3D graphs of solitary wave solution
with fractional order
with parameters
6. Conclusion
In this paper, improved generalized Riccati equation mapping method has been successfully applied to secure new and more general solutions to the space–time fractional Murnaghan’s rod equation. As a result, some totally new solutions have been obtained including several solitary wave solutions: dark, combined dark-bright, singular periodic wave, combined singular periodic wave solutions and one rational solution. These new solutions might be very useful in the study of seismology, physical properties of brass, steel and new elastic materials like polymers. The obtained results show that this mathematical technique is very effective and to find more general solutions of many NPDEs arising in plasma physics, mathematical physics, chemistry, many engineering disciplines and in natural sciences. The physical interpretation of these solutions has been highlighted with the help of 3D graphs and all the solutions obtained in this article have been verified by back substitution in original equation by using Maple 17.
Disclosure statement
No potential conflict of interest was reported by the authors.
References
- K. Diethelm, The analysis of fractional differential equations: An application-oriented exposition using differential operators of caputo type. Berlin: Springer-Verlag, 2010.
- Iqbal M, seadawy AR, Lu D. Construction of solitary wave solutions to the nonlinear modified Kortewege-de Vries dynamical equation in unmagnetized plasma via mathematical methods. Mod Phys Lett A. 2018;33(32):1850183.
- Lu D, Seadawy AR, Iqbal M. Mathematical methods via construction of traveling and solitary wave solutions of three coupled system of nonlinear partial differential equations and their applications. Results Phys. 2018;11:1161–1171.
- Iqbal M, Seadawy AR, Lu D, et al. Construction of a weakly nonlinear dispersion solitary wave solution for the Zakharov–Kuznetsov-modified equal width dynamical equation. Indian J Phys. 2020;94:1465–1474.
- Seadawy AR, Iqbal M, Lu D. Propagation of kink and anti-kink wave solitons for the nonlinear damped modified Korteweg–de Vries equation arising in ion-acoustic wave in an unmagnetized collisional dusty plasma. Physica A. 2020;544:123560.
- Iqbal M, Seadawy AR, Lu D. Dispersive solitary wave solutions of nonlinear further modified Korteweg–de Vries dynamical equation in an unmagnetized dusty plasma. Mod Phys Lett A. 2018;33(37):1850217.
- Iqbal M, Seadawy AR, Lu D. Applications of nonlinear longitudinal wave equation in a magneto-electro-elastic circular rod and new solitary wave solutions. Mod Phys Lett B. 2019;33(18):1950210.
- Seadawy AR, Iqbal M, Lu D. Applications of propagation of long-wave with dissipation and dispersion in nonlinear media via solitary wave solutions of generalized Kadomtsev–Petviashvili modified equal width dynamical equation. Comput Math Appl. 2019;78(11):3620–3632.
- Seadawy AR, Iqbal M, Lu D. Propagation of long-wave with dissipation and dispersion in nonlinear media via generalized Kadomtsive–Petviashvili modified equal width-Burgers equation. Indian J Phys. 2020;94(5):675–687.
- Iqbal M, Seadawy AR, Lu D, et al. Construction of bright–dark solitons and ion-acoustic solitary wave solutions of dynamical system of nonlinear wave propagation. Mod Phys Lett A. 2019;34(37):1950309.
- Seadawy AR, Iqbal M, Lu D. Construction of soliton solutions of the modify unstable nonlinear Schrödinger dynamical equation in fiber optics. Indian J Phys. 2020;94:823–832.
- Seadawy AR, Iqbal M, Lu D. Nonlinear wave solutions of the Kudryashov–Sinelshchikov dynamical equation in mixtures liquid-gas bubbles under the consideration of heat transfer and viscosity. J Taibah Univ Sci. 2019;13(1):1060–1072.
- Iqbal M, Seadawy AR, Khalil OH, et al. Propagation of long internal waves in density stratified ocean for the (2+1)-dimensional nonlinear Nizhnik–Novikov–Vesselov dynamical equation. Results hys. 2020;16:102838.
- Mohyud-Din ST, Noor MA. Homotopy perturbation method for solving partial differential equations. Z Naturforsch A. 2009;64(3–4):157–170.
- Shang Y. A lie algebra approach to susceptible-infected-susceptible epidemics. Electron J Differ Equ. 2012;2012:1–7.
- Shang Y. Lie algebraic discussion for affinity based information diffusion in social networks. Open Phys. 2017;15(1):705–711.
- Sakar MG, Ergören H. Alternative variational iteration method for solving the time-fractional Fornberg–whitham equation. Appl Math Model. 2015;39(14):3972–3979.
- Abdou M, Soliman A. New applications of variational iteration method. Physica D. 2005;211(1):1–8.
- Zhang S, Xia T. An improved generalized F-expansion method and its application to the (2+1)-dimensional KdV equations. Commun Nonlinear Sci Numer Simul. 2008;13(7):1294–1301.
- Navickas Z, Ragulskis M, Listopadskis N, et al. Comments on ‘soliton solutions to fractional-order nonlinear differential equations based on the exp-function method’. Optik. 2017;132:223–231.
- Navickas Z, Telksnys T, Ragulskis M. Comments on ‘The exp-function method and generalized solitary solutions’. Comput Math Appl. 2015;69(8):798–803.
- Zhang S, Xia TC. A further improved extended Fan sub-equation method and its application to the (3+1)-dimensional Kadomstev–Petviashvili equation. Phys Lett A. 2006;356(2):119–123.
- Guner O, Atik H, Kayyrzhanovich AA. New exact solution for space–time fractional differential equations via (G′/G)-expansion method. Optik. 2017;130:696–701.
- Manafian J, Lakestani M. Application of tan(ϕ/2)-expansion method for solving the Biswas–Milovic equation for Kerr law nonlinearity. Optik. 2016;127(4):2040–2054.
- Akbulut A, Kaplan M, Tascan F. The investigation of exact solutions of nonlinear partial differential equations by using exp(−Φ(ξ)) method. Optik. 2017;132:382–387.
- Kudryashov NA. One method for finding exact solutions of nonlinear differential equations. Commun Nonlinear Sci Numer Simul. 2012;17(6):2248–2253.
- Zhu S. The generalizing Riccati equation mapping method in non-linear evolution equation: application to (2+1)-dimensional Boiti–Leon–Pempinelle equation. Chaos Solitons Fractals. 2008;37(5):1335–1342.
- Bekir A, Cevikel AC. The tanh–coth method combined with the Riccati equation for solving nonlinear coupled equation in mathematical physics. J King Saud Univ - Sci. 2011;23(2):127–132.
- Li B, Chen Y, Xuan H, et al. Generalized Riccati equation expansion method and its application to the (3+1)-dimensional Jumbo–Miwa equation. Appl Math Comput. 2004;152(2):581–595.
- Tala-Tebue E, Tsobgni-Fozap DC, Kenfack-Jiotsa A, et al. Envelope periodic solutions for a discrete network with the Jacobi elliptic functions and the alternative (G′/G)-expansion method including the generalized Riccati equation. Eur Phys J Plus. 2014;129(6):136.
- Salathiel Y, Amadou Y, Betchewe G, et al. Soliton solutions and traveling wave solutions for a discrete electrical lattice with nonlinear dispersion through the generalized Riccati equation mapping method. Nonlinear Dyn. 2017;87(4):2435–2443.
- Koonprasert S, Sirisubtawee S, Ampun S. More explicit solitary solutions of the space-time fractional fifth order nonlinear Sawada–Kotera equation via the improved generalized Riccati equation mapping method. 2017.
- Bibi S, Ahmed N, Faisal I, et al. Some new solutions of the Caudrey–Dodd–Gibbon (CDG) equation using the conformable derivative. Adv Differ Equ. 2019;2019(1):89.
- Samsonov A. Strain solitons in solids and how to construct them. Boca Raton (FL): Chapman and Hall/CRC; Dec 18, 2020.
- Erbay HA, Erbay S, Erkip A. Thresholds for global existence and blow-up in a general class of doubly dispersive nonlocal wave equations. Nonlinear Anal: Theory Methods Appl. 2014;95:313–322.
- Cattani C, Sulaiman TA, Baskonus HM, et al. Solitons in an inhomogeneous Murnaghan’s rod. Eur Phys J Plus. 2018;133(6):228.
- Silambarasan R, Baskonus HM, Bulut H. Jacobi elliptic function solutions of the double dispersive equation in the Murnaghan’s rod. Eur Phys J Plus. 2019;134(3):125.
- I. Podlubny, Fractional differential equations: An introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications. Slovak Republic: Academic Press, 1999.
- He J-H, Elagan SK, Li ZB. Geometrical explanation of the fractional complex transform and derivative chain rule for fractional calculus. Phys Lett A. 2012;376(4):257–259.
- Guner O, Bekir A. On the concept of exact solution for nonlinear differential equations of fractional-order. Math Methods Appl Sci. 2016;39(14):4035–4043.
- Abbott S, Weisstein EW. CRC concise encyclopedia of mathematics. Math Gazette. 2007;84(501):549.