Skip to content
BY 4.0 license Open Access Published by De Gruyter Open Access September 24, 2020

Closed-form wave structures of the space-time fractional Hirota–Satsuma coupled KdV equation with nonlinear physical phenomena

  • Md Nur Alam EMAIL logo , Aly R. Seadawy EMAIL logo and Dumitru Baleanu
From the journal Open Physics

Abstract

The present paper applies the variation of ( G / G ) -expansion method on the space-time fractional Hirota–Satsuma coupled KdV equation with applications in physics. We employ the new approach to receive some closed form wave solutions for any nonlinear fractional ordinary differential equations. First, the fractional derivatives in this research are manifested in terms of Riemann–Liouville derivative. A complex fractional transformation is applied to transform the fractional-order ordinary and partial differential equation into the integer order ordinary differential equation. The reduced equations are then solved by the method. Some novel and more comprehensive solutions of these equations are successfully constructed. Besides, the intended approach is simplistic, conventional, and able to significantly reduce the size of computational work associated with other existing methods.

1 Introduction

Nonlinear dynamical systems play an important role in physics. New ideas in physics including astrophysics, computational physics, fluid dynamic, mathematical physics, biological physics and medical physics often explain the fundamental mechanisms studied by other sciences and suggest new avenues of research in academic disciplines such as mathematics. Nature and research problems and most of the life and compound phenomena are explained and invented through fractional wave models. For example, the fluid-dynamic traffic model [1] could succeed the deficiency based upon the hypothesis of continuum traffic flow, the principles of fractional calculus could model the dynamical methods in fluids and porous structures [2] and fractional derivatives can form the nonlinear oscillation of earthquake [3]. For all the attention, it is an indispensable task to attain the closed-form wave structures of the fractional differential equations (FDEs). Various meaningful and truthful approaches have been interpolated for achieving the closed-form wave structures of FDEs, including the exp ( ϕ ( ξ ) ) -expansion approach [4,5], extended Jacobi elliptic function expansion method [6], improved sub-equation scheme [7], modified fractional reduced differential transform method [8], sub-equation method [9], singular manifold method [10], fractional homotopy method [11], fractional reduced differential transform method [12], modified ( G / G ) -expansion approach [13], extended modified mapping method [14], Sine–Gordon expansion method [15], extended trial equation method [16], iterative method [17], simplest equation method [18], ansatz scheme [19], F-expansion method [20], modified Kudryashov method [21], extended mapping method [22], homo separation analysis method [23], modified simple equation method [24], reduced differential transform scheme [25], modified extended mapping method [26], functional variable method [27], extended direct algebraic method [28], Darcy’s law rule [29], function transformation method [30], the variation of ( G / G ) -expansion method [31], differential transform method [32], unified scheme, ( G / G ) -expansion approach [33,34], Kudryashov method [35], new extended Kudryashov process [36], Sine-cosine approach [37], auxiliary equation scheme [38] and many other techniques [39,40,41,42,43,44,45,46,47,48,49,50].

Here, we introduce a new method [31] for nonlinear FDEs based on the homogenous balancing method employing wave transformation. In this research, through the transformation ξ = k x β Γ ( 1 + β ) + λ t β Γ ( 1 + α ) for ordinary differential equation and ξ = x + y + z V t α α for partial differential equation, a given fractional ordinary and a partial differential equation turn into a fractional ordinary differential equation. Suppose the solutions have the form V ( ξ ) = i = 0 M A i L i + i = 1 M B i L i 1 H , where L = G G , H = N N , where G = G ( ξ ) and N = N ( ξ ) represent the solution of the coupled Riccati equation G ( ξ ) = G ( ξ ) N ( ξ ) and N ( ξ ) = 1 N ( ξ ) 2 . This coupled Riccati equation gives us four types of hyperbolic function solutions, such as sech, tanh, csch and coth. We have manipulated the unique method for obtaining more general closed form wave structures of the nonlinear FDEs, such as the space-time fractional Hirota–Satsuma coupled KdV equation [36]. The improvement of this process over the existing rule is that it presents a few novel closed-form wave solutions. Apart from the physical importance, the closed-form wave structures of the nonlinear FDEs might be helpful to the analytical solvers to compare the exactness of their outcomes and also develop the stability analysis of the nonlinear FDEs. In Section 2, some definitions, characters, properties, theorems and fundamental facts of fractional derivatives are introduced. Section 3 shows that the current scheme of the closed-form wave solutions of the proposed models is obtained in Section 4. Discussion and future works are given in Section 5.

2 Preliminaries

Herein, we acquaint a few definitions, characters, properties, theorems and fundamental facts, which are implemented throughout this article.

2.1 The modified Riemann–Liouville derivative

Fractional calculus possesses several ways to generalize the concept of differential derivatives to fractional derivatives [51,52].

Definition 2.1

A real function g ( t ) , t > 0 is called in the space C k where k if there exists a real number p > k , for example, g ( t ) = t p g 1 ( t ) , where g 1 ( t ) C ( 0 , ) and called in the space C k n if g n C k , where n is a positive integer.

Definition 2.2

Suppose that g : x g ( x ) . It denotes a continuous but not significantly differentiable function, then the fractional derivatives of order α are represented through the following expression [51,52]:

(1) D x α = 1 Γ ( α ) 0 x ( x ξ ) α 1 [ G ( ξ ) G ( 0 ) ] d ξ , α < 0 , 1 Γ ( α ) d d x 0 x ( x ξ ) α [ G ( ξ ) G ( 0 ) ] d ξ , 0 < α < 1 , ( G ( n ) ( x ) ) α m , m α m + 1 , m 1 ,

in which Γ ( . ) is the Gamma function illustrated by [53].

(2) Γ ( α ) = lim m m ! m α α ( α + 1 ) ( α + 2 ) ( α + m )

or

(3) Γ ( x ) = 0 e t t x 1 d x

and

(4) D x α = lim h 0 h α k = 0 ( 1 ) k G [ x + ( α k ) h ] .

Definition 2.3

The Mittag–Leffler function including two parameters is represented as follows [54]

(5) E α , β ( x ) = i = 0 x i Γ ( α i + β ) , Re ( α ) > 0 , β , x C .

Some notable characteristics for the fractional derivative are as follows:

  1. D x α x γ = Γ ( 1 + γ ) Γ ( 1 + γ α ) x γ α , γ > 0 .

  2. D x α ( c G ( x ) ) = c D x α ( G ( x ) ) .

  3. D x α ( a G ( x ) + b H ( x ) ) = a D x α G ( x ) + b D x α H ( x ) .

2.2 The fractional complex transformation

This section performs the complex fractional transformation for the fractional-order ordinary differential equation (ODE). First, we explore the nonlinear fractional ODE:

(6) P ( u , D t α u , D x β u , D t α D t α u , D t α D x β u , D x β D x β u , ) = 0 ,

where 0 < α 0 , 0 < β 0 , D t α u and D x β u are the fractional derivatives of u with respect to t and x.

We investigate equation (6) through the transformation u = u ( x , t ) = u ( ξ ) , ξ = k x β Γ ( 1 + β ) + λ t β Γ ( 1 + α ) , where k and λ are nonzero arbitrary constants, then equation (6) becomes

(7) Q = u , u ξ , 2 u ξ 2 , 3 u ξ 3 , = 0 .

3 Glimpse of the method

At this moment, the general scheme of the variation of ( G / G ) -expansion method [31] is shortened as follows:

  • Step 1: Calculate N through rule of the homogeneous balance in equation (7).

  • Step 2: Considering the method can be expressed as the form:

    (8) V ( ξ ) = i = 0 M A i L i + i = 1 M B i L i 1 H ,

    where L = G G , H = N N and G = G ( ξ ) and N = N ( ξ ) represent the solution of the coupled Riccati equations

    (9) G ( ξ ) = G ( ξ ) N ( ξ ) ,

    (10) N ( ξ ) = 1 N ( ξ ) 2 .

    These coupled Riccati equations give us four types of hyperbolic function solutions including sech, tanh, csch and coth such as

    (11) G ( ξ ) = ± sech ( ξ ) , N ( ξ ) = tanh ( ξ ) ,

    (12) G ( ξ ) = ± csch ( ξ ) , N ( ξ ) = coth ( ξ ) .

  • Step 3: A polynomial in L or N is accomplished plugging equation (8) into equation (7). Determining the coefficients of the equivalent power of L or N produces a system of algebraic equations, which can be determined to construct the values of A i and B i using MAPLE. Turning the over measured values of A i and B i in 12, the general solutions of equation (16) and (17) complete the calculation of the result of equation (6).

4 Implementation of the method

To illustrate the idea of the new method, we implement the process on the space-time fractional Hirota–Satsuma coupled KdV equation. We implement the new method to the considered model, which models the intercommunication between two long waves that have well-defined dispersion connection:

(13) D t α u = 1 4 u x x x + 3 u u x + 3 ( v 2 + w ) x D t α v = 1 2 v x x x 3 u v x D t α w = 1 2 w x x x 3 u w x ,

where u = u ( x , t ) , u = u ( x , t ) and u = u ( x , t ) , t > 0 and 0 < α 0 . For our goal, we utilize the transformations

u ( x , t ) = 1 λ u ( ξ ) 2 , v ( x , t ) = λ + u ( ξ ) , w ( x , t ) = 2 λ 2 2 λ u ( ξ ) ,

where ξ = x λ t β Γ ( 1 + α ) . Then equation (13) reduced to the ODE

(14) λ 2 u ξ 2 + 2 u 3 2 λ 2 u .

Applying the homogeneous balance rule on equation (14) yields M = 2 . Therefore,

(15) u ( ξ ) = ( a 0 b 2 ) + b 1 N ( a 1 + b 1 ) N + ( a 2 + b 2 ) N 2 ,

where a 0 , a 1 , a 2 , b 1 and b 2 are free parameters. Plugging equation (15) into equation (14) and calculating each coefficient of L or N to zeros, we have:

  • Case 1: Let λ = 1 , a 0 = b 2 , a 1 = ± 1 , a 2 = b 2 , b 1 = 0 , and insert the values of Case 1 into equation (15), we find:

    (16) u 1 ( x , t ) = 1 + sech 2 2 x + 2 t α Γ ( 1 + α ) ,

    (17) u 2 ( x , t ) = 1 csch 2 2 x + 2 t α Γ ( 1 + α ) ,

    (18) v 1 ( x , t ) = 1 tanh 2 x + 2 t α Γ ( 1 + α ) ,

    (19) v 2 ( x , t ) = 1 coth 2 x + 2 t α Γ ( 1 + α ) ,

    (20) w 1 ( x , t ) = 2 2 tanh 2 x + 2 t α Γ ( 1 + α ) ,

    (21) w 2 ( x , t ) = 2 2 coth 2 x + 2 t α Γ ( 1 + α ) .

  • Case 2: Let λ = 4 , a 0 = b 2 , a 1 = 4 , a 2 = b 2 , b 1 = 2 and put the values of Case 2 into equation (15), we get:

    (22) u 3 ( x , t ) = tanh 2 2 x + 2 t α Γ ( 1 + α ) + 1 2 tanh 2 2 x + 2 t α Γ ( 1 + α ) ,

    (23) u 4 ( x , t ) = coth 2 2 x + 2 t α Γ ( 1 + α ) + 1 2 coth 2 2 x + 2 t α Γ ( 1 + α ) ,

    (24) v 3 ( x , t ) = 4 + 2 tanh x + 4 t α Γ ( 1 + α ) + 2 tanh x + 4 t α Γ ( 1 + α ) ,

    (25) v 4 ( x , t ) = 4 + 2 coth x + 4 t α Γ ( 1 + α ) + 2 coth x + 4 t α Γ ( 1 + α ) ,

    (26) w 3 ( x , t ) = 32 + 16 tanh x + 4 t α Γ ( 1 + α ) + 16 tanh x + 4 t α Γ ( 1 + α ) ,

    (27) w 4 ( x , t ) = 32 + 16 coth x + 4 t α Γ ( 1 + α ) + 16 coth x + 4 t α Γ ( 1 + α ) .

  • Case 3: Let λ = 1 , a 0 = b 2 , a 1 = 1 , a 2 = b 2 , b 1 = ± 1 and substitute the values of Case 3 into equation (15), we achieve:

    (28) u 5 ( x , t ) = coth 2 2 x + 2 t α Γ ( 1 + α ) ,

    (29) u 6 ( x , t ) = tanh 2 2 x + 2 t α Γ ( 1 + α ) ,

    (30) v 5 ( x , t ) = 1 1 tanh 2 x + 2 t α Γ ( 1 + α ) ,

    (31) v 6 ( x , t ) = 1 1 coth 2 x + 2 t α Γ ( 1 + α ) ,

    (32) w 5 ( x , t ) = 2 2 tanh 2 x + 2 t α Γ ( 1 + α ) ,

    (33) w 6 ( x , t ) = 2 2 coth 2 x + 2 t α Γ ( 1 + α ) .

  • Case 4: Let λ = 4 , a 0 = b 2 , a 1 = 4 , a 2 = b 2 , b 1 = 2 and use the values of Case 4 into equation (15), we find:

    (34) u 7 ( x , t ) = 4 + sech 2 2 x + 2 t α Γ ( 1 + α ) csch 2 2 x + 2 t α Γ ( 1 + α ) ,

    (35) u 8 ( x , t ) = 4 + sech 2 2 x + 2 t α Γ ( 1 + α ) csch 2 2 x + 2 t α Γ ( 1 + α ) ,

    (36) v 7 ( x , t ) = 4 2 tanh x + 4 t α Γ ( 1 + α ) 2 tanh x + 4 t α Γ ( 1 + α ) ,

    (37) v 8 ( x , t ) = 4 2 coth x + 4 t α Γ ( 1 + α ) 2 coth x + 4 t α Γ ( 1 + α ) ,

    (38) w 7 ( x , t ) = 32 16 tanh x + 4 t α Γ ( 1 + α ) 16 tanh x + 4 t α Γ ( 1 + α ) ,

    (39) w 8 ( x , t ) = 32 16 coth x + 4 t α Γ ( 1 + α ) 16 coth x + 4 t α Γ ( 1 + α ) .

  • Case 5: Let λ = 2 , a 0 = b 2 , a 1 = 0 , a 2 = b 2 , b 1 = ± 2 with b 2 being a free parameter and plug the values of Case 5 into equation (15), we obtain:

(40) u 9 ( x , t ) = tanh 2 2 x + 2 t α Γ ( 1 + α ) coth 2 2 x + 2 t α Γ ( 1 + α ) ,

(41) u 10 ( x , t ) = coth 2 2 x + 2 t α Γ ( 1 + α ) tanh 2 2 x + 2 t α Γ ( 1 + α ) ,

(42) v 9 ( x , t ) = 2 ± 2 1 tanh x 2 t α Γ ( 1 + α ) tanh x 2 t α Γ ( 1 + α ) ,

(43) v 10 ( x , t ) = 2 ± 2 1 coth x 2 t α Γ ( 1 + α ) coth x 2 t α Γ ( 1 + α ) ,

(44) w 9 ( x , t ) = 8 4 2 tanh x 2 t α Γ ( 1 + α ) ¯ tanh x 2 t α Γ ( 1 + α )

(45) w 10 ( x , t ) = 8 4 2 ( coth x 2 t α Γ ( 1 + α ) coth x 2 t α Γ ( 1 + α ) ) .

A picture is an indispensable tool for conversation and to demonstrate the results to the difficulties lucidly. When performing the computation in daily life, we require the fundamental understanding of constructing the use of pictures. Therefore, the graphical presentations of few obtained results are depicted in Figures 1–6.

Figure 1 
               The picture of the result in 
                     
                        
                        
                           
                              
                                 u
                              
                              
                                 1
                              
                           
                           (
                           x
                           ,
                           t
                           )
                        
                        {u}_{1}(x,t)
                     
                   with the values 
                     
                        
                        
                           α
                           =
                           0.35
                        
                        \alpha =0.35
                     
                  , 
                     
                        
                        
                           
                              
                                 b
                              
                              
                                 2
                              
                           
                           =
                           0.5
                        
                        {b}_{2}=0.5
                     
                   and 
                     
                        
                        
                           t
                           =
                           0.01
                        
                        t=0.01
                     
                   for 
                     
                        
                        
                           2
                           D
                        
                        2\text{D}
                     
                   graphics. (a) Real 3D surface, (b) complex 3D surface, (c) 2D shape, (d) real contour shape, (e) complex contour shape.
Figure 1

The picture of the result in u 1 ( x , t ) with the values α = 0.35 , b 2 = 0.5 and t = 0.01 for 2 D graphics. (a) Real 3D surface, (b) complex 3D surface, (c) 2D shape, (d) real contour shape, (e) complex contour shape.

Figure 2 
               The picture of the result in 
                     
                        
                        
                           
                              
                                 u
                              
                              
                                 2
                              
                           
                           (
                           x
                           ,
                           t
                           )
                        
                        {u}_{2}(x,t)
                     
                   with the values 
                     
                        
                        
                           α
                           =
                           0.35
                        
                        \alpha =0.35
                     
                  , 
                     
                        
                        
                           
                              
                                 b
                              
                              
                                 2
                              
                           
                           =
                           0.5
                        
                        {b}_{2}=0.5
                     
                   and 
                     
                        
                        
                           t
                           =
                           0.01
                        
                        t=0.01
                     
                   for 
                     
                        
                        
                           2
                           D
                        
                        2\text{D}
                     
                   graphics. (a) Real 3D surface, (b) complex 3D surface, (c) 2D shape, (d) real contour shape, (e) complex contour shape.
Figure 2

The picture of the result in u 2 ( x , t ) with the values α = 0.35 , b 2 = 0.5 and t = 0.01 for 2 D graphics. (a) Real 3D surface, (b) complex 3D surface, (c) 2D shape, (d) real contour shape, (e) complex contour shape.

Figure 3 
               The picture of the result in 
                     
                        
                        
                           
                              
                                 v
                              
                              
                                 1
                              
                           
                           (
                           x
                           ,
                           t
                           )
                        
                        {v}_{1}(x,t)
                     
                   with the values 
                     
                        
                        
                           α
                           =
                           0.35
                        
                        \alpha =0.35
                     
                  , 
                     
                        
                        
                           
                              
                                 b
                              
                              
                                 2
                              
                           
                           =
                           0.5
                        
                        {b}_{2}=0.5
                     
                   and 
                     
                        
                        
                           t
                           =
                           0.01
                        
                        t=0.01
                     
                   for 
                     
                        
                        
                           2
                           D
                        
                        2\text{D}
                     
                   graphics. (a) Real 3D surface, (b) complex 3D surface, (c) 2D shape, (d) real contour shape, (e) complex contour shape.
Figure 3

The picture of the result in v 1 ( x , t ) with the values α = 0.35 , b 2 = 0.5 and t = 0.01 for 2 D graphics. (a) Real 3D surface, (b) complex 3D surface, (c) 2D shape, (d) real contour shape, (e) complex contour shape.

Figure 4 
               The picture of the result in 
                     
                        
                        
                           
                              
                                 v
                              
                              
                                 2
                              
                           
                           (
                           x
                           ,
                           t
                           )
                        
                        {v}_{2}(x,t)
                     
                   with the unknown parameter values 
                     
                        
                        
                           α
                           =
                           0.35
                        
                        \alpha =0.35
                     
                  , 
                     
                        
                        
                           
                              
                                 b
                              
                              
                                 2
                              
                           
                           =
                           0.5
                        
                        {b}_{2}=0.5
                     
                   and 
                     
                        
                        
                           t
                           =
                           0.01
                        
                        t=0.01
                     
                   for 
                     
                        
                        
                           2
                           D
                        
                        2\text{D}
                     
                   graphics. (a) Real 3D surface, (b) complex 3D surface, (c) 2D shape, (d) real contour shape, (e) complex contour shape.
Figure 4

The picture of the result in v 2 ( x , t ) with the unknown parameter values α = 0.35 , b 2 = 0.5 and t = 0.01 for 2 D graphics. (a) Real 3D surface, (b) complex 3D surface, (c) 2D shape, (d) real contour shape, (e) complex contour shape.

Figure 5 
               The picture of the result in 
                     
                        
                        
                           
                              
                                 v
                              
                              
                                 6
                              
                           
                           (
                           x
                           ,
                           t
                           )
                        
                        {v}_{6}(x,t)
                     
                   with the unknown parameter values 
                     
                        
                        
                           α
                           =
                           0.35
                        
                        \alpha =0.35
                     
                   and 
                     
                        
                        
                           
                              
                                 b
                              
                              
                                 2
                              
                           
                           =
                           0.5
                        
                        {b}_{2}=0.5
                     
                   and 
                     
                        
                        
                           t
                           =
                           0.01
                        
                        t=0.01
                     
                   for 
                     
                        
                        
                           2
                           D
                        
                        2\text{D}
                     
                   graphics. (a) Real 3D surface, (b) complex 3D surface, (c) 2D shape, (d) real contour shape, (e) complex contour shape.
Figure 5

The picture of the result in v 6 ( x , t ) with the unknown parameter values α = 0.35 and b 2 = 0.5 and t = 0.01 for 2 D graphics. (a) Real 3D surface, (b) complex 3D surface, (c) 2D shape, (d) real contour shape, (e) complex contour shape.

Figure 6 
               The picture of the result in 
                     
                        
                        
                           
                              
                                 v
                              
                              
                                 10
                              
                           
                           (
                           x
                           ,
                           t
                           )
                        
                        {v}_{10}(x,t)
                     
                   with the unknown parameter values 
                     
                        
                        
                           α
                           =
                           0.35
                        
                        \alpha =0.35
                     
                  , 
                     
                        
                        
                           
                              
                                 b
                              
                              
                                 2
                              
                           
                           =
                           0.5
                        
                        {b}_{2}=0.5
                     
                   and 
                     
                        
                        
                           t
                           =
                           0.01
                        
                        t=0.01
                     
                   for 
                     
                        
                        
                           2
                           D
                        
                        2\text{D}
                     
                   graphics. (a) Real 3D surface, (b) complex 3D surface, (c) 2D shape, (d) real contour shape, (e) complex contour shape.
Figure 6

The picture of the result in v 10 ( x , t ) with the unknown parameter values α = 0.35 , b 2 = 0.5 and t = 0.01 for 2 D graphics. (a) Real 3D surface, (b) complex 3D surface, (c) 2D shape, (d) real contour shape, (e) complex contour shape.

5 Discussions and future work

In this research, the suggested method has been triumphantly displayed to attain new closed-form wave solutions of equation (13). It has been determined that the complex fractional transformation and the advanced method are an essential and significant mathematical device in analyzing closed-form wave structures of a whole class of fractional FDEs. Consequently, a few new closed-form wave answers of the equation are determined.

The compensations and legality of the recommended approach over the extended Kudryashov approach are as follows. The critical compensation of the proposed method over the extended Kudryashov approach is that the proposed method implements numerous general and plentiful new closed-form wave structures. The closed-form wave solutions of nonlinear FDE have its vital importance to show the complicated physical aspects. For example, [36] extends the Kudryashov method to solve equation Q ξ ( ξ ) = Q 3 ( ξ ) Q ( ξ ) as an auxiliary equation and the closed-form solutions presented by u ( ξ ) = i = 0 n α i Q i ( ξ ) , where Q ( ξ ) = ± 1 1 ± e 2 ξ . Our solutions u 1 , v 1 , w 1 , u 2 , v 2 , w 2 , u 5 and u 6 are similar to solutions of [36] and solutions u 3 , v 3 , w 3 , u 4 , v 4 , w 4 , v 5 , w 5 , v 6 , w 6 , u 7 , v 7 , w 7 , u 8 , v 8 , w 8 , u 9 , v 9 , w 9 , u 10 , v 10 and w 10 are all new exact solutions obtained in this paper, which validates our proposed methods. Although we have presented diverse closed-form wave solutions by performing the recommended method, our obtained solutions show that our proposed methods are more helpful and extraordinary in contributing much more general and many new closed-form wave solutions at the same time. To state the effectiveness and present the insight of processes and the comparison among our proposed method, the new ansatz method is much more proficient than the new extended Kudryashov method. We can underline from our understanding that the method can be performed in other nonlinear FDEs and can reduce the amount of computational work. Therefore, the investigation of closed-form wave structures of other nonlinear FDEs deserves further analysis.

Acknowledgements

The authors would like to acknowledge CAS-TWAS president’s fellowship program. The author thanks the referees for their suggestions and comments.

  1. Competing interests: This research received no specific grant from any funding agency in the public, commercial or not-for-profit sectors. The authors did not have any competing interests in this research.

  2. Author contributions: All parts contained in the research carried out by the authors through hard work and a review of the various references and contributions in the field of mathematics and Applied physics.

  3. Availability of data and materials: Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.

  4. Funding: No funding for this article.

References

[1] Podlubny I . Fractional Differential Equations. New York: Academic Press; 1999.Search in Google Scholar

[2] Miller KS , Ross B . An Introduction to the Fractional Calculus and Fractional Differential Equations. New York: John Wiley; 1993.Search in Google Scholar

[3] Mainardi F . Fractional Calculus Some Basic Problems in Continum and Statistical Mechanics. New York: Springer-Verlag; 1997.10.1007/978-3-7091-2664-6Search in Google Scholar

[4] Alam M , Belgacem F . Microtubules nonlinear models dynamics investigations through the exp(−ϕ(ξ))-expansion method implementation. Mathematics. 2016;4:6.10.3390/math4010006Search in Google Scholar

[5] Alam M , Tunc C . An analytical method for solving exact solutions of the nonlinear bogoyavlenskii equation and the nonlinear diffusive predator-prey system. Alex Eng J. 2016;55:1855–65.10.1016/j.aej.2016.04.024Search in Google Scholar

[6] Sonmezoglu A . Exact solutions for some fractional differential equations. Adv Math Phys. 2015;2015:567842 10.1155/2015/567842Search in Google Scholar

[7] Karaagac B . New exact solutions for some fractional order differential equations via improved sub-equation method. Discret Cont Dyn Syst – S. 2019;12:447–54.10.3934/dcdss.2019029Search in Google Scholar

[8] 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.10.1142/S0217732318501833Search in Google Scholar

[9] Martineza HY , Aguilarb JFG , Baleanu D . Beta-derivative and sub-equation method applied to the optical solitons in medium with parabolic law nonlinearity and higher order dispersion. Optik. 2018;155:357–65.10.1016/j.ijleo.2017.10.104Search in Google Scholar

[10] Saleh R , Kassem M , Mabrouk SM . Exact solutions of nonlinear fractional order partial differential equations via singular manifold method. Chin J Phys. 2019;61:290–300.10.1016/j.cjph.2019.09.005Search in Google Scholar

[11] Delgado VFM , Aguilar JFG , Torres L , Jimenez RFE , Hernandez MAT . Exact solutions for the liénard type model via fractional homotopy methods. Fract Derivatives Mittag–Leffler Kernel. 2019;194:269–91.10.1007/978-3-030-11662-0_16Search in Google Scholar

[12] Abuasad S , Moaddy K , Hashim I . Analytical treatment of two-dimensional fractional helmholtz equations. J King Saud Univ – Sci. 2019;31:659–66.10.1016/j.jksus.2018.02.002Search in Google Scholar

[13] Alam MN , Tunc C . Constructions of the optical solitons and others soliton to the conformable fractional zakharov–kuznetsov equation with power law nonlinearity. J Taibah Univ Med Sci. 2020;15:263–72. 10.1080/16583655.2019.1708542.Search in Google Scholar

[14] Seadawy AR . Three-dimensional weakly nonlinear shallow water waves regime and its travelling wave solutions. Int J Comput Methods. 2018;15:1850017.10.1142/S0219876218500172Search in Google Scholar

[15] Kadkhoda N , Jafari H . An analytical approach to obtain exact solutions of some space-time conformable fractional differential equations. Adv Differ Equ. 2019;2019;428.10.1186/s13662-019-2349-0Search in Google Scholar

[16] Seadawy AR , Manafian J . New soliton solution to the longitudinal wave equation in a magneto-electro-elastic circular rod. Results Phys. 2018;8:1158–67.10.1016/j.rinp.2018.01.062Search in Google Scholar

[17] Lu D , Seadawy AR , Ali A . Dispersive traveling wave solutions of the equal-width and modified equal-width equations via mathematical methods and its applications. Results Phys. 2018;9:313–20.10.1016/j.rinp.2018.02.036Search in Google Scholar

[18] Helal MA , Seadawy AR , Zekry MH . Stability analysis of solitary wave solutions for the fourth-order nonlinear Boussinesq water wave equation. Appl Math Comput. 2014;232:1094–103.10.1016/j.amc.2014.01.066Search in Google Scholar

[19] Shi D , Zhang Y . Diversity of exact solutions to the conformable space-time fractional mew equation. Appl Math Lett. 2020;99:105994.10.1016/j.aml.2019.07.025Search in Google Scholar

[20] Seadawy AR , El-Rashidy K . Dispersive solitary wave solutions of kadomtsev–petviashivili and modified kadomtsev–petviashivili dynamical equations in unmagnetized dust plasma. Results Phys. 2018;8:1216–22.10.1016/j.rinp.2018.01.053Search in Google Scholar

[21] Lu D , Seadawy AR , Ali A . Applications of exact traveling wave solutions of modified Liouville and the symmetric regularized long wave equations via two new techniques. Results Phys. 2018;9:1403–10.10.1016/j.rinp.2018.04.039Search in Google Scholar

[22] Arshad M , Seadawy A , Lu D . Elliptic function and solitary wave solutions of the higher-order nonlinear Schrodinger dynamical equation with fourth-order dispersion and cubic-quintic nonlinearity and its stability. Eur Phys J Plus. 2017;132:371.10.1140/epjp/i2017-11655-9Search in Google Scholar

[23] Zuriqat M . Exact solution for the fractional partial differential equation by homo separation analysis method. Afr Matematika. 2019;30:1133–43.10.1007/s13370-019-00707-xSearch in Google Scholar

[24] Arnous AH , Seadawy AR , Alqahtani RT , Biswas A . Optical solitons with complex ginzburg-landau equation by modified simple equation method. Opt – Int J Light Electron Opt. 2017;144:475–80.10.1016/j.ijleo.2017.07.013Search in Google Scholar

[25] Owyed S , Abdou MA , Abdel-Atyb AH , Nekhili WAR . Numerical and approximate solutions for coupled time fractional nonlinear evolutions equations via reduced differential transform method. Chaos, Solitons Fractals. 2019;217:109474.10.1016/j.chaos.2019.109474Search in Google Scholar

[26] Abdullah , Seadawy AR , Wang J . Mathematical methods and solitary wave solutions of three-dimensional zakharov–kuznetsov–burgers equation in dusty plasma and its applications. Results Phys. 2017;7:4269–77.10.1016/j.rinp.2017.10.045Search in Google Scholar

[27] Cenesiz Y , Tasbozan O , Kurt A . Functional variable method for conformable fractional modified kdv-zkequation and maccari system. Tbilisi Math J. 2017;10:117–25.10.1515/tmj-2017-0010Search in Google Scholar

[28] Seadawy AR , Lu D , Yue C . Travelling wave solutions of the generalized nonlinear fifth-order kdv water wave equations and its stability. J Taibah Univ Sci. 2017;11:623–33.10.1016/j.jtusci.2016.06.002Search in Google Scholar

[29] Sheikholeslami M . Cuo-water nanofluid free convection in a porous cavity considering darcy law. Eur Phys J Plus. 2017;132:132.10.1140/epjp/i2017-11330-3Search in Google Scholar

[30] Khater AH , Callebaut DK , Seadawy AR . General soliton solutions for nonlinear dispersive waves in convective type instabilities. Phys Scr. 2006;74:384–93.10.1088/0031-8949/74/3/015Search in Google Scholar

[31] Shehata AR , Abu-Amra SSM . Geometrical properties and exact solutions of three (3 + 1) -dimensional nonlinear evolution equations in mathematical physics using different expansion methods. J Adv Math Comput Sci. 2019;33:1–19.10.9734/jamcs/2019/v32i430149Search in Google Scholar

[32] Unal E , Gokdogan A . Solution of conformable fractional ordinary differential equations via differential transform method. Opt – Int J Light Electron Opt. 2017;128:264–73.10.1016/j.ijleo.2016.10.031Search in Google Scholar

[33] Alam MN , Li X . Exact traveling wave solutions to higher order nonlinear equations. J Ocean Eng Sci. 2019;4:276–88.10.1016/j.joes.2019.05.003Search in Google Scholar

[34] Arshad M , Seadawy A , Lu D . Modulation stability and optical soliton solutions of nonlinear Schrodinger equation with higher order dispersion and nonlinear terms and its applications. Superlattices Microstructures. 2017;112:422–34.10.1016/j.spmi.2017.09.054Search in Google Scholar

[35] Eslami M . Exact traveling wave solutions to the fractional coupled nonlinear schrodinger equations. Appl Math Comput. 2016;285:141–8.10.1016/j.amc.2016.03.032Search in Google Scholar

[36] Ege S , Misirli E . A new method for solving nonlinear fractional differential equations. N Trends Math Sci. 2017;5:225–33.10.20852/ntmsci.2017.141Search in Google Scholar

[37] Wazwaz A . Exact solutions with solitons and periodic structures for the zakharov–kuznetsov equation and its modified form. Commun Nonlinear Sci Numer Simul. 2005;10:597–606.10.1016/j.cnsns.2004.03.001Search in Google Scholar

[38] Akbulut MKA , Bekir A . Auxiliary equation method for fractional differential equations with modified riemann–liouville derivative. Int J Nonlin Sci Num. 2016;17:413–20.10.1515/ijnsns-2016-0023Search in Google Scholar

[39] Aly Seadawy , El-Rashidy K . Dispersive solitary wave solutions of Kadomtsev–Petviashivili and modified Kadomtsev–Petviashivili dynamical equations in unmagnetized dust plasma. Results Phys. 2018;8:1216–22.10.1016/j.rinp.2018.01.053Search in Google Scholar

[40] Kumar S , Kumar R , Singh J , Nisar KS , Kumard D . An efficient numerical scheme for fractional model of HIV-1 infection of CD4 + T-cells with the effect of antiviral drug therapy. Alex Eng J. 2020;216:634–42. 10.1016/j.aej.2019.12.046 Search in Google Scholar

[41] Kumar S , Kumar A , Momani S , Aldhaifallah M , Nisar KS . Numerical solutions of nonlinear fractional model arising in the appearance of the stripe patterns in two-dimensional systems. Adv Differ Equ. 2019;2019(1):413.10.1186/s13662-019-2334-7Search in Google Scholar

[42] Kumar S , Kumar A , Abbas S , Qurashi MA , Baleanu D . A modified analytical approach with existence and uniqueness for fractional Cauchy reaction-diffusion equations. Adv Differ Equ. 2020;2020(1):1–8.10.1186/s13662-019-2488-3Search in Google Scholar

[43] Sharma B , Kumar S , Cattani C , Baleanu D . Nonlinear dynamics of Cattaneo–Christov heat flux model for third-grade power-law fluid. J Comput Nonlinear Dyn. 2020 Jan 1;15(1):927–41.10.1115/1.4045406Search in Google Scholar

[44] Kumar S , Ghosh S , Samet B , Goufo EFD . An analysis for heat equations arises in diffusion process using new Yang–AbdelAty–Cattani fractional operator. Math Methods Appl Sci. 2020;162:572–84.Search in Google Scholar

[45] Tassaddiq A , Khan I , Nisar KS . Heat transfer analysis in sodium alginate based nanofluid using MoS2 nanoparticles: Atangana–Baleanu fractional model. Chaos, Solitons Fractals. 2020;130:109445.10.1016/j.chaos.2019.109445Search in Google Scholar

[46] Gao W , Senel M , Yel G , Baskonus HM , Senel B . New complex wave patterns to the electrical transmission line model arising in network system. Aims Math. 2020;5(3):1881–92.10.3934/math.2020125Search in Google Scholar

[47] Gao W , Ismael HF , Husien AM , Bulut H , Baskonus HM . Optical soliton solutions of the nonlinear Schrödinger and resonant nonlinear Schrödinger equation with Parabolic law. Appl Sci. 2020;10(1):219 , 1–20.10.3390/app10010219Search in Google Scholar

[48] Gao W , Yel G , Baskonus HM , Cattani C . Complex solitons in the conformable (2 + 1)-dimensional Ablowitz–Kaup–Newell–Segur equation. Aims Math. 2020;5(1):507–21.10.3934/math.2020034Search in Google Scholar

[49] Gao W , Rezazadeh H , Pinar Z , Baskonus HM , Sarwar S , Yel G . Novel explicit solutions for the nonlinear zoomeron equation by using newly extended direct algebraic technique. Opt Quant Electron. 2020;52(52):1–13.10.1007/s11082-019-2162-8Search in Google Scholar

[50] Cordero A , Jaiswal JP , Torregrosa JR . Stability analysis of fourth-order iterative method for finding multiple roots of non-linear equations. Appl Math Nonlinear Sci. 2019;4(1):43–56.10.2478/AMNS.2019.1.00005Search in Google Scholar

[51] Jumarie G . Table of some basic fractional calculus formulae derived from a modified riemann–liouville derivative for non-differentiable functions. Appl Math Lett. 2009;22:378–85.10.1016/j.aml.2008.06.003Search in Google Scholar

[52] Guner O . Singular and non-topological soliton solutions for nonlinear fractional differential equations. Chin Phys B. 2015;24:100201.10.1088/1674-1056/24/10/100201Search in Google Scholar

[53] Hartley T , Lorenzo C . Dynamics and control of initialized fractional-order systems. Nonlinear Dyn. 2002;29:201–33.10.1023/A:1016534921583Search in Google Scholar

[54] Kilbas A , Srivastava H , Trujillo J . Theory and applications of fractional differential equations. Amsterdam: Elsevier; 2006.Search in Google Scholar

Received: 2020-05-27
Revised: 2020-07-29
Accepted: 2020-07-29
Published Online: 2020-09-24

© 2020 Md Nur Alam et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 25.4.2024 from https://www.degruyter.com/document/doi/10.1515/phys-2020-0179/html
Scroll to top button