Research

# Symmetry properties and explicit solutions of the nonlinear time fractional KdV equation

Gangwei Wang* and Tianzhou Xu

Author Affiliations

School of Mathematics, Beijing Institute of Technology, Beijing, 100081, P.R. China

For all author emails, please log on.

Boundary Value Problems 2013, 2013:232  doi:10.1186/1687-2770-2013-232

 Received: 8 September 2013 Accepted: 2 October 2013 Published: 8 November 2013

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

### Abstract

The time fractional KdV equation in the sense of the Riemann-Liouville derivatives is considered. The symmetry properties of the time fractional KdV equation is investigated by using the Lie group analysis method. On the basis of the point symmetry, the vector fields of the time fractional KdV equation are presented. And then, the symmetry reductions are constructed. By right of the obtained Lie point symmetries, it is shown that this equation could transform into a nonlinear ordinary differential equation of fractional order with the new independent variable ξ = x t α / 3 . The derivative is an Erdélyi-Kober derivative depending on a parameter α. At last, by means of the sub-equation method, some exact and explicit solutions of the time fractional KdV equation are constructed.

MSC: 22E70, 26A33.

##### Keywords:
fractional KdV equation; modified Riemann-Liouville derivative; Lie symmetry analysis; Erdélyi-Kober operators; exact solutions

### 1 Introduction

Lie point transformation groups have been extensively applied to the nonlinear partial differential equations arising in mathematics, physics and in many other scientific fields. Lie point transformation groups can deal with symmetry reductions, similarity solutions of nonlinear differential partial equations. It can be said that the method of Lie symmetry groups is the most important approach for constructing analytical solutions of nonlinear partial differential equations. A huge number of papers and many excellent textbooks (see, e.g., [1-16] and papers cited therein) have been devoted to the theory of Lie symmetry groups and their applications to differential equations. On the other hand, it is important to note that a very small number of them involve Lie group analysis to solve problems for fractional differential equations (FDEs) [17-23]. As is known to all, there is no existing general method for dealing with fractional differential equations. Lie group analysis, however, is a powerful, systematic and direct method for investigating fractional differential equations. It is a fact that the classical Lie theory applied to fractional differential equations (FDEs) for symmetry analysis should be modified since it fails in the case of fractional differential equations.

Since Leibniz’s note in his letter to L’Hospital, in which the meaning of the derivative of order one half was firstly discussed, the theory of derivatives of non-integer order has stimulated considerable interest in the areas of mathematics, physics and engineering [24-28]. Recently, the fractional calculus has been used successfully to describe many complex nonlinear phenomena and dynamic processes in physics, engineering, electromagnetics, acoustics, viscoelasticity, electrochemistry, material science, etc.[24-40]. In reality, the next state of a physical phenomenon might depend not only on its current state but also on its historical states (non-local property), which can be successfully modeled by using the theory of derivatives and integrals of fractional order [21,22]. Because of this, many effective analytic methods have been developed to derive exact, explicit and numerical solutions of nonlinear fractional partial differential equations (FPDEs), such as the Adomian decomposition method [29,30], the invariant subspace method [31], the transform method [32,33], the homotopy perturbation method [34], the variational iteration method [35], the sub-equation method [36-38], and so on.

In this paper, by means of the Lie symmetry group method, we consider the following nonlinear fractional KdV (FK) equation with time fractional derivatives of the form

α u t α + a u u x + b u x x x = 0 , (1)

where 0 < α 1 , a, c are parameters, α is a parameter describing the order of the fractional time derivative. The function u ( x , t ) is assumed to be a causal function of time. In physical applications, the special cases such as α = 1 2 and α = 1 3 may be more useful. In the case of α = 1 , the fractional equation reduces to the classical nonlinear KdV equation. In [21], invariant analysis of the time fractional generalized Korteweg-de Vries equation is investigated. However, the authors do not give explicit solutions.

The paper is organized as follows. In Section 2, a brief review of the main definitions and properties of fractional calculus and general concept of the invariance criterion for FDEs are given to provide a convenient reference. In Section 3, we perform Lie group classification on the time fractional KdV equation, the general similarity forms and symmetry reductions are established. In Section 4, some exact solutions are discussed. The conclusions are presented in Section 5.

### 2 Preliminaries

We state some notations and theorems. First, let us recall essentials of fractional calculus. The fractional calculus is a name for the theory of integrals and derivatives of arbitrary order. It unifies and generalizes the notions of integer-order differentiation and n-fold integration. Besides the R-L definition of fractional derivatives, there are several other different definitions, such as the modified R-L (mR-L) derivative, the Grünwald-Letnikov derivative (G-L) and Caputo’s fractional derivative, and so on. For different circumstances, they can be used for dealing with different properties of physical models. For example, Caputo’s fractional derivative is related to initial value problems; generally speaking, the mR-L derivative is used to investigate exact and explicit solutions of some FDEs [36-38]. The modified Riemann-Liouville derivative was defined by Jumarie [39,40].

Definition 2.1

D t α f ( t ) = { 1 Γ ( 1 α ) d d t 0 t ( t ξ ) α ( f ( ξ ) f ( 0 ) ) d ξ , 0 < α < 1 , [ f ( n ) ( t ) ] ( α n ) , n α < n + 1 , n 1 , (2)

where the Euler gamma function Γ ( z ) is defined by the integral

Γ ( z ) = 0 e t t z 1 d t , (3)

which converges in the right half of the complex plane Re ( z ) > 0 .

Properties 2.1

D t α t γ = Γ ( γ + 1 ) Γ ( γ + 1 α ) t γ α , γ > 0 , (4)

D t α [ u ( t ) v ( t ) ] = u ( t ) D t α v ( t ) + v ( t ) D t α u ( t ) , (5)

D t α [ f ( u ( t ) ) ] = f u [ u ( t ) ] D t α u ( t ) = D u α f [ u ( t ) ] ( u t ) α , (6)

which will be used in the following sections.

In what follows, considering the general case of a nonlinear system of partial differential equations of order n in p independent and q dependent variables [1,2,4-7]

Δ v ( x , u ( n ) ) = 0 , v = 1 , 2 , 3 , , l , (7)

here x = ( x 1 , x 2 , , x p ) , u = ( u 1 , u 2 , , u q ) , and the derivatives of u in reference to x up to n, where u ( n ) represents all the derivatives of u of all orders from 0 to n. Consider a one-parameter Lie group of infinitesimal transformations acting on system (7),

x ¯ i = x i + ϵ ξ i ( x , u ) + O ( ϵ 2 ) , i = 1 , , p , u ¯ j = u j + ϵ η j ( x , u ) + O ( ϵ 2 ) , j = 1 , , q , (8)

where ϵ is the parameter. The vector field V is associated with the above group of transformations as follows:

V = i = 1 p ξ i ( x , u ) x i + α = 1 q η α ( x , u ) u α . (9)

The invariance of system (7) under the infinitesimal transformations leads to the invariance conditions

pr ( n ) V [ Δ v ( x , u ( n ) ) ] = 0 , v = 1 , 2 , 3 , , l , (10)

whenever Δ v ( x , u ( n ) ) = 0 , where pr ( n ) V is called the nth order prolongation of the infinitesimal generator given by

pr ( n ) V = V + α = 1 q J η α J ( x , u ( n ) ) u j α , (11)

where J = ( j 1 , j 2 , , j k ) , with 1 j k p , 1 k n . The coefficient functions η α J of pr ( n ) V are given by the following formula:

η α J ( x , u ( n ) ) = D J ( η α i = 1 p ξ i u i α ) + i = 1 p ξ i u J , i α , (12)

where u i α = u α / x i and u J , i α = u J α / x i .

However, for a nonlinear system of fractional partial differential equations (FPDEs), the prolongation formula will be different. Now, we present below brief details of the Lie symmetry analysis to FPDEs with respect to two independent variables.

Consider a scalar time FPDE having the following form [21,22]:

α u t α = F ( x , t , u , u x , u x x , ) . (13)

If (13) is invariant under a one-parameter Lie group of point transformations

t = t + ϵ τ ( x , t , u ) + O ( ϵ 2 ) , x = x + ϵ ξ ( x , t , u ) + O ( ϵ 2 ) , u = u + ϵ η ( x , t , u ) + O ( ϵ 2 ) , α u ¯ t ¯ α = α u t α + ϵ η α 0 ( x , t , u ) + O ( ϵ 2 ) , u ¯ x ¯ = u x + ϵ η x ( x , t , u ) + O ( ϵ 2 ) , 2 u ¯ x ¯ 2 = 2 u x 2 + ϵ η x x ( x , t , u ) + O ( ϵ 2 ) , (14)

where

η α 0 = D t α ( η ) + ξ D t α ( u x ) D t α ( ξ u x ) + D t α ( D t ( τ ) u ) D t α + 1 ( τ u ) + τ D t α + 1 ( u ) = α η t α + ( η u α D t ( τ ) ) α u t α u α η u t α + μ + n = 1 [ ( a n ) α η u t α ( a n + 1 ) D t n + 1 ( τ ) ] D t α n ( u ) n = 1 ( a n ) D t n ( ξ ) D t α n ( u x ) , (15)

where

μ = n = 2 m = 2 n k = 2 m r = 0 k 1 ( a n ) ( n m ) ( k r ) 1 k ! t n α Γ ( n + 1 α ) [ u ] r m t m [ u k r ] n m + k η t n m u k (16)

and

η x = D x ( η ) u x D x ( ξ ) u t D x ( τ ) , η x x = D x ( η x ) u x t D x ( τ ) u x x D x ( ξ ) , η x x x = D x ( η x x ) u x x t D x ( τ ) u x x x D x ( ξ ) , (17)

Here, D x denotes the total derivative operator and it is defined by

D x = x + u x u + u x x u x + (18)

with the associated vector field of the form

V = τ ( x , t , u ) t + ξ ( x , t , u ) x + η ( x , t , u ) u , (19)

where the coefficient functions ξ ( x , t , u ) , τ ( x , t , u ) , and η ( x , t , u ) of the vector field are to be determined.

If vector field (19) generates a symmetry of (1), then V must satisfy Lie’s symmetry condition

pr ( n ) V ( Δ 1 ) | Δ 1 = 0 = 0 , (20)

where Δ 1 = α u t α F ( x , t , u , u x , u x x , ) .

### 3 Symmetry group analysis of the nonlinear time fractional KdV equation

In the preceding section, we have given some definitions and formulas about the Lie symmetry analysis method on the FPDEs. In the present section, through the above discussion and the Lie theory, we investigate the time-fractional KdV equation.

According to the Lie theory, applying the third prolongation pr ( 3 ) V to (1), we can get the following system of symmetry equations:

η α 0 + a u η x + a η u x + b η x x x = 0 . (21)

Substituting (14), (15) and (17) into (21), and equating the coefficients of the various monomials in partial derivatives with respect to x and various power of u, one can find the determining equations for the symmetry group of Eq. (1). Solving these equations, we obtain the following forms of the coefficient functions:

ξ = c 1 x + c 2 , τ = 3 c 1 t α , η = 2 c 1 u , (22)

where c 1 and c 2 are arbitrary constants. Therefore, we can obtain the corresponding vector fields

V = 3 c 1 t α t + ( c 1 x + c 2 ) x 2 c 1 u u . (23)

Thus, the Lie point symmetries are represented by the Lie group generators (vector fields)

V 1 = x , V 2 = 3 t α t + x x 2 u u . (24)

Moreover, it is easily checked that the vector fields are closed under the Lie bracket, respectively,

[ V 1 , V 2 ] = V 1 , [ V 2 , V 1 ] = V 1 . (25)

For the symmetry V 2 , we get the characteristic equation

d x x = α d t 3 t = d u 2 u , (26)

and the corresponding invariants are

ξ = x t α 3 , u = t 2 α 3 g ( ξ ) . (27)

From what has been discussed above, one can see that (1) can be reduced to a nonlinear ODE of fractional order with a new independent variable ξ = x t α / 3 . Consequently, one can get the following theorem.

TheoremTransformation (27) reduces (1) to the following nonlinear ordinary differential equation of fractional order:

( P 3 α 1 5 α 3 , α g ) ( ξ ) + a g g ξ + b g ξ ξ ξ = 0 , (28)

with the Erdélyi-Kober fractional differential operator P β τ , α of order[28]

( P β τ , α g ) : = j = 0 n 1 ( τ + j 1 β ξ d d ξ ) ( K β τ + α , n α g ) ( ξ ) , (29)

n = { [ α ] + 1 , α N , α , α N , (30)

where

( K β τ , α g ) ( ξ ) : = { 1 Γ ( α ) 1 ( u 1 ) α 1 u ( τ + α ) g ( ξ u 1 β ) d u , α > 0 , g ( ξ ) , α = 0 , (31)

is the Erdélyi-Kober fractional integral operator.

Proof Let n 1 < α < n , n = 1 , 2 , 3 ,  . Considering the Riemann-Liouville fractional derivative, one can get

α u t α = n t n [ 1 Γ ( n α ) 0 t ( t s ) n α 1 s α 3 g ( x s α 3 ) d s ] . (32)

Let v = t s , one can get d s = t v 2 d v , thus (32) can write

α u t α = n t n [ t n 5 α 3 1 Γ ( n α ) 1 ( v 1 ) n α 1 v ( n + 1 5 α 3 ) g ( ξ v α 3 ) d v ] . (33)

Based on the Erdélyi-Kober fractional integral operator (31), one can have

α u t α = n t n [ t n 5 α 3 ( K α 3 1 2 α 3 , n α g ) ( ξ ) ] . (34)

In view of the relation ( ξ = x t α 3 ), we can get

t t ϕ ( ξ ) = t x ( α 3 ) t α 3 1 ϕ ( ξ ) = α 3 ξ ξ ϕ ( ξ ) . (35)

Thus, one can obtain

n t n [ t n 5 α 3 ( K 3 α 1 2 α 3 , n α g ) ( ξ ) ] = n 1 t n 1 [ t ( t n 5 α 3 ( K 3 α 1 2 α 3 , n α g ) ( ξ ) ) ] = n 1 t n 1 [ t n 5 α 3 ( n 5 α 3 α 3 ξ ξ ( K 3 α 1 2 α 3 , n α g ) ( ξ ) ) ] . (36)

Repeating the similar procedure as above for n 1 times, one can get

n t n [ t n 5 α 3 ( K 3 α 1 2 α 3 , n α g ) ( ξ ) ] = n 1 t n 1 [ t ( t n 5 α 3 ( K 3 α 1 2 α 3 , n α g ) ( ξ ) ) ] = n 1 t n 1 [ t n 5 α 3 ( n 5 α 3 α 3 ξ ξ ( K 3 α 1 2 α 3 , n α g ) ( ξ ) ) ] = = t 5 α 3 j = 0 n 1 ( 1 5 α 3 + j α 3 ξ d d ξ ) ( K 3 α 1 2 α 3 , n α g ) ( ξ ) . (37)

Now using (29), we get

n t n [ t n 5 α 3 ( K 3 α 1 2 α 3 , n α g ) ( ξ ) ] = t 5 α 3 ( P 3 α 1 5 α 3 , n α g ) ( ξ ) . (38)

Substituting (38) into (34), one can get

α u t α = t 5 α 3 ( P 3 α 1 5 α 3 , n α g ) ( ξ ) . (39)

Therefore, the time fractional KdV equation can be reduced into an FODE (fractional ordinary differential equation)

( P 3 α 1 5 α 3 , α g ) ( ξ ) + a g g ξ + b g ξ ξ ξ = 0 . (40)

This completes the proof. □

In particular, if α = 1 2 , then the symmetries of the time fractional Korteweg-de Vries equation are as follows:

V = 6 c 1 t t + ( c 1 x + c 2 ) x 2 c 1 u u .

In this case, the time fractional KdV equation can be reduced into an FODE

( P 6 1 6 , 1 2 g ) ( ξ ) + a g g ξ + b g ξ ξ ξ = 0 .

Similarly, in view of (23), the symmetries can be obtained successively in terms of the specific parameters α, such as α = 1 3 , α = 1 4 and so on.

Remark 1 Through the above discussion, we can find that the point symmetries of the time fractional KdV equation are relatively fewer than those for the generalized KdV equation. The main reason is that the fractional order 0 < α 1 is an arbitrary parameter in our discussion. It is appropriate to mention here that the fractional order 0 < α 1 affects the properties of the equation greatly.

### 4 Exact and explicit solutions of the time fractional KdV equation

#### 4.1 Fractional sub-equation method

Now, we outline the main steps of the fractional sub-equation method for solving fractional differential equations.

For a given NFDE, consider two variables x and t,

P ( u , u t , u x , D t α u , D x α u , ) = 0 , 0 < α 1 , (41)

where D x α u and D t α u are the modified Riemann-Liouville derivatives of u with respect to t and x, respectively.

Step 1: By making use of the traveling wave transformation

u ( x , t ) = u ( ξ ) , ξ = x + c t , (42)

where c is a nonzero constant to be determined later, (41) can be reduced to a nonlinear fractional ordinary differential equation (NFODE)

P ( u , c u , u , c α D ξ α u , D ξ α u , ) = 0 , 0 < α 1 . (43)

Step 2: Suppose that Eq. (43) has the following solution:

u ( ξ ) = a 0 + i = 1 n a i ( ψ ( ξ ) ) i , (44)

where a i ( i = 1 , , n ) are constants to be determined later, positive integer n can be determined by balancing the highest order derivatives and nonlinear terms in Eq. (41) or Eq. (43). The function ψ ( ξ ) satisfies the following Bäklund transformation of the fractional Riccati equation [38]:

ψ ( ξ ) = σ B + D ϕ ( ξ ) D + B ϕ ( ξ ) , (45)

where B, D are arbitrary parameters, and B 0 . Meanwhile, ϕ ( ξ ) are decided by

D ξ α ϕ ( ξ ) = σ + ϕ 2 ( ξ ) , (46)

where σ is a constant. Eq. (46) has the following solutions:

ϕ ( ξ ) = { σ tanh ( σ ξ , α ) , σ < 0 , σ coth ( σ ξ , α ) , σ < 0 , σ tan ( σ ξ , α ) , σ > 0 , σ cot ( σ ξ , α ) , σ > 0 , Γ ( 1 + α ) ξ α + ω , ω  is constant , σ = 0 , (47)

with the generalized hyperbolic and trigonometric functions

sin α ( ξ ) = E α ( i ξ α ) E α ( i ξ α ) 2 i , cos α ( ξ ) = E α ( i ξ α ) + E α ( i ξ α ) 2 i , sinh α ( ξ ) = E α ( ξ α ) E α ( ξ α ) 2 , cosh α ( ξ ) = E α ( ξ α ) + E α ( ξ α ) 2 , tan α ( ξ ) = sin α ( ξ ) cos α ( ξ ) , cot α ( ξ ) = cos α ( ξ ) sin α ( ξ ) , tanh α ( ξ ) = sinh α ( ξ ) cosh α ( ξ ) , coth α ( ξ ) = cosh α ( ξ ) sinh α ( ξ ) , (48)

here E α ( ξ ) = k = 0 ξ k Γ ( 1 + k α ) ( α > 0 ) is the Mittag-Leffler function in one parameter.

Step 3: Substituting (44), (45) and (46) into (43) and setting the coefficients of the powers of ( ϕ ( ξ ) ) i to be zero, one can obtain an over-determined nonlinear algebraic system in a i ( i = 1 , , n ) and c.

Step 4: With the aid of Maple, solving the nonlinear algebraic system yields the explicit expressions of the parameters a i ( i = 1 , , n ) and c. Then substituting these constants and the solutions of Eq. (47) into Eq. (44), we can get the exact and explicit solutions of the nonlinear fractional partial differential equation (NFPDE) (41).

#### 4.2 Applications to the time fractional KdV equation

According to the above steps, firstly, we introduce the following transformations:

u ( x , t ) = u ( ξ ) , ξ = x + c t , (49)

where c is a constant. Substituting (49) into (1), then (1) can be reduced to the following nonlinear fractional ordinary differential equation (NFODE):

c α D ξ α u + a u u ξ + b u ξ ξ ξ = 0 . (50)

We suppose that Eq. (50) has the following solution:

u ( ξ ) = a 0 + i = 1 n a i ( ψ ( ξ ) ) i , (51)

where a i ( i = 1 , , n ) are constants to be determined later. Balancing the highest order derivative terms with nonlinear terms in Eq. (50), we get

u ( ξ ) = a 0 + a 1 ( ψ ) + a 2 ( ψ ) 2 . (52)

Substituting (52) along with (45) into (50) and then letting the coefficients of ( ϕ ) i to zero, one can get some algebraic equations about c, a 0 , a 1 and a 2 . Solving the algebraic equations by Maple, one can get the following.

Case 1:

B = B , D = D , a = a , b = b , c = ( a a 0 ) 1 α , α = α , σ = 0 , a 0 = a 0 , a 1 = 0 , a 2 = 12 b a . (53)

Case 2:

B = B , D = D , a = a , b = b , c = c , α = α , σ = D 2 B 2 , a 0 = a 0 , a 1 = a 1 , a 2 = a 2 . (54)

Case 3:

B = B , D = D , a = 12 b a 2 , b = b , c = ( 4 b ( 2 a 2 σ 3 a 0 ) a 2 ) 1 α , α = α , σ = σ , a 0 = a 0 , a 1 = 0 , a 2 = a 2 . (55)

Case 4:

B = B , D = 0 , a = a , b = 1 12 a a 2 , c = ( 2 3 a σ a 2 a a 0 ) 1 α , α = α , σ = σ , a 0 = a 0 , a 1 = 0 , a 2 = a 2 . (56)

In view of (54), one can get new types of explicit solutions of Eq. (1) as follows:

u = a 0 + a 1 ( σ B D σ tanh ( σ ξ , α ) D B σ tanh ( σ ξ , α ) ) + a 2 ( σ B D σ tanh ( σ ξ , α ) D B σ tanh ( σ ξ , α ) ) 2 , (57)

where σ < 0 , ξ = x + c t ,

u = a 0 + a 1 ( σ B D σ coth ( σ ξ , α ) D B σ coth ( σ ξ , α ) ) + a 2 ( σ B D σ coth ( σ ξ , α ) D B σ coth ( σ ξ , α ) ) 2 , (58)

where σ < 0 , ξ = x + c t ,

u = a 0 + a 1 ( σ B + D σ tan ( σ ξ , α ) D + B σ tan ( σ ξ , α ) ) + a 2 ( σ B + D σ tan ( σ ξ , α ) D + B σ tan ( σ ξ , α ) ) 2 , (59)

where σ > 0 , ξ = x + c t ,

u = a 0 + a 1 ( σ B D σ cot ( σ ξ , α ) D B σ cot ( σ ξ , α ) ) + a 2 ( σ B D σ cot ( σ ξ , α ) D B σ cot ( σ ξ , α ) ) 2 , (60)

where σ > 0 , ξ = x + c t .

If σ = 0 , one can get u = 0 .

Using (55), one can get new types of explicit solutions of Eq. (1) as follows:

u = a 0 + a 2 ( σ B D σ tanh ( σ ξ , α ) D B σ tanh ( σ ξ , α ) ) 2 , (61)

where σ < 0 , ξ = x + c t ,

u = a 0 + a 2 ( σ B D σ coth ( σ ξ , α ) D B σ coth ( σ ξ , α ) ) 2 , (62)

where σ < 0 , ξ = x + c t ,

u = a 0 + a 2 ( σ B + D σ tan ( σ ξ , α ) D + B σ tan ( σ ξ , α ) ) 2 , (63)

where σ > 0 , ξ = x + c t ,

u = a 0 + a 2 ( σ B D σ cot ( σ ξ , α ) D B σ cot ( σ ξ , α ) ) 2 , (64)

where σ > 0 , ξ = x + c t ,

u = a 0 + a 2 ( D Γ ( 1 + α ) D ( ξ α + ω ) + B Γ ( 1 + α ) ) 2 , (65)

where σ = 0 , ξ = x + c t .

Remark 2 Using (53) and (56), we can also get other exact solutions of (1). Here we do not list all of them.

Remark 3 To the best of our knowledge, the solutions obtained in this paper have not been reported in previous literature. Therefore, these solutions are new.

Remark 4 It is interesting to note that if α = 1 , FDEs (1) can be reduced to the conventional integer order KdV equation, and the obtained exact solutions can be reduced to the conventional hyperbolic and trigonometric functions.

### 5 Concluding remarks

In the present study, we investigated the symmetry properties, similarity reduction forms and explicit solutions of the time fractional KdV equation by using the Lie symmetry groups. Lie point symmetry groups of the governing equations were obtained. The obtained point transformation groups, however, for the time fractional KdV equation are fewer than those for the generalized KdV equation. The main reason is that the fractional order 0 < α 1 is an arbitrary parameter in our model. Based on the obtained Lie point symmetries, we have shown that this equation can transform into a nonlinear ordinary differential equation of fractional order with the new independent variable ξ = x t α / 3 . The derivative is an Erdélyi-Kober derivative depending on a parameter α. At last, some exact and explicit solutions of the equations are presented. These solutions include generalized hyperbolic function solutions, generalized trigonometric function solutions and rational function solutions and so on. These solutions may be useful to further investigate the complicated nonlinear physical phenomena.

### Competing interests

The authors declare that they have no competing interests.

### Authors’ contributions

All authors read and approved the final manuscript.

### Acknowledgements

The authors express their sincere thanks to the referees for their careful review of this manuscript and their useful suggestions. This work was completed with the support of the National Natural Science Foundation of China (NNSFC) (Grant No. 11171022).

### References

1. Olver, PJ: Application of Lie Group to Differential Equation, Springer, New York (1986)

2. Ovsiannikov, LV: Group Analysis of Differential Equations, Academic Press, New York (1982)

3. Lie, S: On integration of a class of linear partial differential equations by means of definite integrals. Arch. Math. Log.. 6(3), 328–368 (1881)

4. Bluman, GW, Kumei, S: Symmetries and Differential Equations, Springer, New York (1989)

5. Ibragimov NH (ed.): CRC Handbook of Lie Group Analysis of Differential Equations, CRC Press, Boca Raton (1994)

6. Ibragimov NH (ed.): CRC Handbook of Lie Group Analysis of Differential Equations, CRC Press, Boca Raton (1995)

7. Ibragimov NH (ed.): CRC Handbook of Lie Group Analysis of Differential Equations, CRC Press, Boca Raton (1996)

8. Liu, HZ, Geng, YX: Symmetry reductions and exact solutions to the systems of carbon nanotubes conveying fluid. J. Differ. Equ.. 254, 2289–2303 (2013). Publisher Full Text

9. Craddock, M, Lennox, K: Lie symmetry methods for multi-dimensional parabolic PDEs and diffusions. J. Differ. Equ.. 252, 56–90 (2012). Publisher Full Text

10. Kumar, S, Singh, K, Gupta, RK: Painlevé analysis, Lie symmetries and exact solutions for ( 2 + 1 ) -dimensional variable coefficients Broer-Kaup equations. Commun. Nonlinear Sci. Numer. Simul.. 17, 1529–1541 (2012). Publisher Full Text

11. Vaneeva, O: Lie symmetries and exact solutions of variable coefficient mKdV equations: an equivalence based approach. Commun. Nonlinear Sci. Numer. Simul.. 17, 611–618 (2012). Publisher Full Text

12. Naz, R, Khan, MD, Naeem, I: Conservation laws and exact solutions of a class of non linear regularized long wave equations via double reduction theory and Lie symmetries. Commun. Nonlinear Sci. Numer. Simul.. 18, 826–834 (2013). Publisher Full Text

13. Listopadovaa, V, Magdab, O, Pobyzhc, V: How to find solutions, Lie symmetries, and conservation laws of forced Korteweg-de Vries equations in optimal way. Nonlinear Anal., Real World Appl.. 14, 202–205 (2013). Publisher Full Text

14. Johnpillaia, AG, Karab, AH, Biswas, A: Symmetry reduction, exact group-invariant solutions and conservation laws of the Benjamin-Bona-Mahoney equation. Appl. Math. Lett.. 26, 376–381 (2013). Publisher Full Text

15. Jefferson, GF: On the second-order approximate symmetry classification and optimal systems of subalgebras for a forced Korteweg-de Vries equation. Commun. Nonlinear Sci. Numer. Simul.. 18, 2340–2358 (2013). Publisher Full Text

16. Cherniha, R, Pliukhin, O: New conditional symmetries and exact solutions of reaction-diffusion-convection equations with exponential nonlinearities. J. Math. Anal. Appl.. 403, 23–37 (2013). Publisher Full Text

17. Gazizov, RK, Kasatkin, AA, Lukashchuk, YS: Continuous transformation groups of fractional differential equations. Vestn. USATU. 9, 125–135 (in Russian) (2007)

18. Gazizov, RK, Kasatkin, AA, Lukashchuk, SY: Symmetry properties of fractional diffusion equations. Phys. Scr. T. 136, Article ID 014016 (2009)

19. Buckwar, E, Luchko, Y: Invariance of a partial differential equation of fractional order under the Lie group of scaling transformations. J. Math. Anal. Appl.. 227, 81–97 (1998). Publisher Full Text

20. Djordjevic, VD, Atanackovic, TM: Similarity solutions to nonlinear heat conduction and Burgers/Korteweg-de Vries fractional equations. J. Comput. Appl. Math.. 212, 701–714 (2008)

21. Sahadevan, R, Bakkyaraj, T: Invariant analysis of time fractional generalized Burgers and Korteweg-de Vries equations. J. Math. Anal. Appl.. 393, 341–347 (2012). Publisher Full Text

22. Wang, GW, Liu, XQ, Zhang, YY: Lie symmetry analysis to the time fractional generalized fifth-order KdV equation. Commun. Nonlinear Sci. Numer. Simul.. 18, 2321–2326 (2013). Publisher Full Text

23. Liu, HZ: Complete group classifications and symmetry reductions of the fractional fifth-order KdV types of equations. Stud. Appl. Math. (2013). Publisher Full Text

24. Diethelm, K: The Analysis of Fractional Differential Equations, Springer, Berlin (2010)

25. Miller, KS, Ross, B: An Introduction to the Fractional Calculus and Fractional Differential Equations, Wiley, New York (1993)

26. Podlubny, I: Fractional Differential Equations, Academic Press, San Diego (1999)

27. Oldham, KB, Spanier, J: The Fractional Calculus, Academic Press, San Diego (1974)

28. Kiryakova, V: Generalised Fractional Calculus and Applications (1994)

29. El-Sayed, AMA, Gaber, M: The Adomian decomposition method for solving partial differential equations of fractal order in finite domains. Phys. Lett. A. 359, 175–182 (2006). Publisher Full Text

30. Chen, Y, An, HL: Numerical solutions of coupled Burgers equations with time- and space-fractional derivatives. Appl. Math. Comput.. 200, 87–95 (2008). Publisher Full Text

31. Gazizov, RK, Kasatkin, AA: Construction of exact solutions for fractional order differential equations by the invariant subspace method. Comput. Math. Appl. (2013). Publisher Full Text

32. Odibat, Z, Momani, S: A generalized differential transform method for linear partial differential equations of fractional order. Appl. Math. Lett.. 21, 194–199 (2008). Publisher Full Text

33. Li, X, Chen, W: Analytical study on the fractional anomalous diffusion in a half-plane. J. Phys. A, Math. Theor.. 43(49), Article ID 495206 (2010)

34. He, JH: A coupling method of a homotopy technique and a perturbation technique for non-linear problems. Int. J. Non-Linear Mech.. 35, 37–43 (2000). Publisher Full Text

35. Wu, G, Lee, EWM: Fractional variational iteration method and its application. Phys. Lett. A. 374, 2506–2509 (2010). Publisher Full Text

36. Zhang, S, Zhang, HQ: Fractional sub-equation method and its applications to nonlinear fractional PDEs. Phys. Lett. A. 375, 1069–1073 (2011). Publisher Full Text

37. Guo, S, Mei, LQ, Li, Y, Sun, YF: The improved fractional sub-equation method and its applications to the space-time fractional differential equations in fluid mechanics. Phys. Lett. A. 376, 407–411 (2012). Publisher Full Text

38. Lu, B: Bäklund transformation of fractional Riccati equation and its applications to nonlinear fractional partial differential equations. Phys. Lett. A. 376, 2045–2048 (2012). Publisher Full Text

39. Jumarie, G: Modified Riemann-Liouville derivative and fractional Taylor series of nondifferentiable functions further results. Comput. Math. Appl.. 51, 1367–1376 (2006). Publisher Full Text

40. Jumarie, G: Cauchy’s integral formula via the modified Riemann-Liouville derivative for analytic functions of fractional order. Appl. Math. Lett.. 23, 1444–1450 (2010). Publisher Full Text