Next Article in Journal
Numerical Study for Magnetohydrodynamic Flow of Nanofluid Due to a Rotating Disk with Binary Chemical Reaction and Arrhenius Activation Energy
Previous Article in Journal
Convolutional Neural Network for Copy-Move Forgery Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lie Symmetry Analysis of the Time Fractional Generalized KdV Equations with Variable Coefficients

1
School of Science, Xi’an University of Posts and Telecommunications, Xi’an 710121, China
2
School of Mathematics and Statistics, Xi’an Jiaotong University, Xi’an 710049, China
3
School of Electronic Engineering, Xi’an University of Posts and Telecommunications, Xi’an 710121, China
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(10), 1281; https://doi.org/10.3390/sym11101281
Submission received: 9 August 2019 / Revised: 3 October 2019 / Accepted: 9 October 2019 / Published: 14 October 2019

Abstract

:
The group classification of a class of time fractional generalized KdV equations with variable coefficient is presented. The Lie symmetry analysis method is extended to the certain subclasses of time fractional generalized KdV equations with initial and boundary values. Under the corresponding similarity transformation with similarity invariants, KdV equations with initial and boundary values have been transformed into fractional ordinary differential equations with initial value. Then we use the power series method to obtain the exact solution of the reduced equation with the Erdélyi-Kober fractional differential operator.

1. Introduction

In recent years, the fractional partial differential equations (FPDEs) have been widely studied. They can represent numerous intricate phenomena in many fields, such as engineering, mechanics, mechatronics, physics, control theory, and other fields of science [1,2,3,4,5] well. Consequently, some methods have been proposed to find exact solutions to FPDEs, such as the variational method [6,7], the Homotopy analysis method [8], the sub-equation method [9], the Lie symmetry analysis method [10,11,12,13], and so on.
The Korteweg-de Vries (KdV) equation is a very important mathematical model that was firstly derived by Korteweg and de Vries in 1895 and it has been used to describe the evolution and interaction of nonlinear waves [14,15]. Many scholars are also devoted to the study of the time fractional KdV equations. In this article, we are concerned with the time fractional generalized KdV equations of the form
α u t α + a u m u x + g ( t ) u x x x = 0 , x > 0 , t > 0 ,
where 0 < α < 1 , a is a nonzero constant, m ( m 1 ) is a positive integer, and g ( t ) is an arbitrary smooth non-vanishing function. In Equation (1), α u t α is the Riemann-Liouville fractional derivative [2,3] with respect to t, which is defined as
α u t α = 1 Γ ( n α ) n t n 0 t u ( s , x ) ( t s ) α + 1 n d s , ( 0 n 1 < α < n , n N ) , n u t n , ( α = n N ) ,
where Γ ( z ) = 0 + e t t z 1 d t is Gamma function.
When α = 1 , Equation (1) can be reduced to variable coefficient generalized KdV equations, which is studied by way of the Lie group analysis method in detail in Reference [16]. The exhaustive group classification was presented. O.O. Vaneeva et al. applied the Lie symmetry method, which was proposed by Bluman [17] and the generalized KdV equations with the initial and boundary value was reduced to nonlinear third-order ODEs with initial value.
This paper is organized as follows. In Section 2, Lie symmetry analysis method of the fractional partial differential equation is presented. In Section 3, the group classification of the time fractional generalized KdV equations with variable coefficients is carried out by using the classical Lie symmetry analysis method. Then in Section 4, we first try to extend the Lie symmetry analysis method to the time fractional generalized KdV equations with initial and boundary value and use the power series method to obtain the exact solution. Finally, the conclusions are given.

2. Lie Symmetry Analysis of the Time Fractional Partial Differential Equation

Consider the following time fractional partial differential equation
α u t α F ( t , x , u , u x , u x x , ) = 0
with two independent variables t , x and one dependent variable u = u ( t , x ) . The one-parameter( ε ) Lie group of transformations is given by
t = t + ε · τ ( t , x , u ) + o ( ε 2 ) , x = x + ε · ς ( t , x , u ) + o ( ε 2 ) , u = u + ε · η ( t , x , u ) + o ( ε 2 ) , α u t α = α u t α + ε · η α 0 ( t , x , u ) + o ( ε 2 ) , u x = u x + ε · η x ( t , x , u ) + o ( ε 2 ) , 2 u x 2 = u x x + ε · η x x ( t , x , u ) + o ( ε 2 ) ,
where
η x = η x + ( η u ς x ) u x τ x u t ς u ( u x ) 2 τ u u x u t , η x x = η x x + ( 2 η x u ς x x ) u x τ x x u t ( η u u 2 ς x u ) ( u x ) 2 2 τ x u u x u t ς u u ( u x ) 3 τ u u ( u x ) 2 u t + ( η u 2 ς x ) u x x 2 τ x u x t 3 ς u u x x u x τ u u x x u t 2 τ u u x t u x , η x x x = D x ( η x x ) u x x t D x ( τ ) u x x x D x ( ς ) ,
η α 0 = α η t α + η u α D t ( τ ) α u t α u α η u t α + μ + n = 1 [ n α n η u t n n + 1 α D t n + 1 ( τ ) ] D t α n ( u ) n = 1 n α D t n ( ς ) D t α n ( u x ) ,
here
μ = n = 2 m = 2 n k = 2 m r = 0 k 1 n α m n r k 1 k ! t n α Γ ( n + 1 α ) ( u ) r m t m ( u k r ) n m + k η t n m u k .
(see [12,13,14,15,18,19,20,21,22]).
Definition 1.
Under transformations (3), Equation (2) is invariant (or Equation(2) admits transformations (3)) if Equation (2) has the same form with respect to new variables t , x , u , i.e.,
α u t α F ( t , x , u , u x , u x x , ) = 0 .
On the basis of the Lie theory [17,23,24,25], we can obtain
Theorem 1.
Equation (2) is invariant under a one-parameter (ε) Lie group of infinitesimal transformations if and only if
P r V ( k ) ( α u t α F ( t , x , u , u x , u x x , ) ) α u t α F ( t , x , u , u x , u x x , ) = 0 = 0 .
where P r V ( k ) is k order prolongation infinitesimal operator of
V = ς ( t , x , u ) x + τ ( t , x , u ) t + η ( t , x , u ) u .
This theorem is called criterion of invariance.

3. Lie Symmetry Analysis of the Time Fractional Generalized KdV Equations with Variable Coefficients

In this section, we carry out the group classification of Equation (1). The third prolongation of the infinitesimal operator P r V ( 3 ) is as follows
P r V ( 3 ) = V + η α 0 ( α u t α ) + η x u x + η x x x u x x x .
Applying P r V ( 3 ) to Equation (1), one can get
τ ( g ( t ) ) t u x x x + a m η u m 1 u x + a η x u m + η x x x g ( t ) + η α 0 = 0 .
Substituting the expressions of η α 0 , η x , η x x x into Equation (4), we sort it out in the form of power of partial derivatives of u. Then equating the coefficients to zero, we get a system of linear PDEs and FPDEs, called the determining equations, as follows
τ x = τ u = ς u = ς t = η u u = 0 , 3 g ( t ) ς x = α ( τ ( t ) ) t g ( t ) + τ ( t ) ( g ( t ) ) t , g ( t ) ( η x u ς x x ) = 0 , a ( α τ t ς x ) u m + a m η u m 1 + g ( t ) ( 3 η x x u ς x x x ) = 0 , a u m η x + g ( t ) η x x x + α η t α u α η u t α = 0 , n α n η u t n n + 1 α D t ( n + 1 ) ( τ ) = 0 , n = 1 , 2 , .
From the first five equations of the determining equations, we can obtain
τ = τ ( t ) , ς = ς ( x ) , η = η 1 ( t , x ) u + η 0 ( t , x ) .
By analyzing the second equation of the determining equations, we have
ς = a 1 x + a 2 ,
where a 1 , a 2 are arbitrary constants. Then the second equation is reduced to
3 g ( t ) a 1 = α ( τ ( t ) ) t g ( t ) + τ ( t ) ( g ( t ) ) t .
To solve the above differential equation, τ ( t ) is given as
τ ( t ) = g ( t ) 1 α 3 a 1 α g ( t ) 1 α d t + a 3 ,
where a 3 is a constant.
By analyzing the third equation of the determining equations, we can get η x = 0 , then the fourth equation is reduced to
( α τ t ς x ) u + m η = 0 .
Namely,
η = a 1 α ( τ ( t ) ) t m u .
According to n α = ( 1 ) n 1 α Γ ( n α ) Γ ( 1 α ) Γ ( n + 1 ) , Γ ( n + 1 ) = n Γ ( n ) , the last equation of the determining equations is written as
( n + 1 ) n η u t n + ( n α ) D t ( n + 1 ) ( τ ) = 0 , n = 1 , 2 , .
When n = 1 , equation is
2 η u t + ( 1 α ) D t 2 ( τ ) = 0 ,
which gives
2 α m ( τ ( t ) ) t t t + ( 1 α ) ( τ ( t ) ) t t = 0 .
The invariance condition τ ( t ) t = 0 = 0 under the Lie group transformations (3), because of the conservative property of the Riemann-Liouville fractional derivative operator.
Next, we will carry out group classification according to different forms of function g ( t ) as follows.
(1)
If g ( t ) is arbitrary, we obtain the explicit form of infinitesimals
τ = 0 , ς = a 1 , η = 0
where a 1 is arbitrary constant. Hence the infinitesimal operator becomes
V = a 1 x .
(2)
If g ( t ) is a nonzero constant, the infinitesimals are
τ = 3 a 1 α t , ς = a 1 x + a 2 , η = 2 a 1 m u ,
where a 1 , a 2 are arbitrary constants. Hence, the general form of the infinitesimal operator is
V = 3 a 1 α t t + ( a 1 x + a 2 ) x 2 a 1 m u u .
A two-dimensional Lie algebra is spanned by
V 1 = 3 α t t + x x 2 m u u , V 2 = x .
(3)
If g ( t ) = t p ( p 0 , p α ), the expressions of infinitesimals are given as
τ = 3 a 1 p + α t , ς = a 1 x + a 2 , η = p 2 α m ( p + α ) a 1 u
where a 1 and a 2 are arbitrary constants. The infinitesimal operator takes the form
V = 3 a 1 p + α t t + ( a 1 x + a 2 ) x + p 2 α m ( p + α ) a 1 u u
We obtain a two-dimensional Lie algebra, spanned by
V 1 = 3 t p + α t + x x + p 2 α m ( p + α ) u u , V 2 = x .
When p = 2 α , the infinitesimals are given as
τ = a 1 α t , ς = a 1 x + a 2 , η = 0
where a 1 and a 2 are arbitrary constants. The corresponding infinitesimal generators are
V 1 = t α t + x x , V 2 = x .
When p = α , the infinitesimals are given as
τ = a 1 t , ς = a 2 , η = α a 1 m u
where a 1 and a 2 are arbitrary constants. The corresponding infinitesimal generators are
V 1 = t t α u m u , V 2 = x .
(4)
If g ( t ) = e p t (p is a nonzero real number), we obtain the following expressions of infinitesimals
τ = 3 a 1 p , ς = a 1 x + a 2 , η = a 1 m u ,
where a 1 and a 2 are arbitrary constants. Then we have
V = 3 a 1 p t + ( a 1 x + a 2 ) x + a 1 m u u .
We have a two-dimensional Lie algebra, which is spanned by
V 1 = 3 p t + x x + u m u , V 2 = x .

4. Lie Symmetry Analysis of Time Fractional Generalized KdV Equations with Initial and Boundary Value

4.1. Lie Symmetry Analysis of Time Fractional Partial Differential Equations with Initial and Boundary Value

In applications, the partial differential equation with the initial and boundary value has been paid more attention. In Reference [17], Bluman successfully applied the Lie symmetry analysis method to deal with the partial differential equations with the initial and boundary value. In this paper we will attempt this method to solve the following time fractional partial differential equations with initial and boundary value
α u t α F ( t , x , u , u x , u x x , ) = 0 ,
u | S = h ( x ) ,
where S is the manifold.
Assuming that Equation (5) admits l one-parameter symmetries V i , i = 1 , 2 , , l , we say that Equations (5) and (6) are invariant under a symmetry of the form
V = i = 1 l a i V i
for some constants a i . The following two conditions also need to be satisfied
(i)
under symmetry V, the manifold S is invariant;
(ii)
under symmetry V restricted to manifold S, the initial (boundary) condition u | S = h ( x ) is invariant.

4.2. Reduction for Time Fractional Generalized KdV Equations with Initial and Boundary Value

We first apply this method to the time fractional generalized KdV equation with initial and boundary value as follows
α u t α + a u m u x + t p u x x x = 0 , x > 0 , t > 0 ,
u ( 0 , x ) = 0 , x > 0 , u ( t , 0 ) = q ( t ) , t > 0 , u x ( t , 0 ) = 0 , t > 0 , u x x ( t , 0 ) = 0 , t > 0 ,
where 0 < α < 1 , a and p α are arbitrary nonzero constants, m ( m 1 ) is a positive integer, q ( t ) is a non-vanishing smooth function with respect to t.
Then we search for the infinitesimal generators, which leave initial and boundary conditions (8) of Equation (7) invariant.
Theorem 2.
The infinitesimal generator
V = 3 t p + α t + x x + p 2 α m ( p + α ) u u
leaves Equations (7) and (8) invariant, and boundary condition u ( t , 0 ) should satisfy
q ( t ) = λ t p 2 α 3 m , t > 0 ,
where λ > 0 is a constant.
Proof. 
In order to find infinitesimal generator admitted by Equation (7) with initial and boundary conditions (8), we construct a special linear combination,
V = a 1 V 1 + a 2 V 2 ,
where a 1 and a 2 are constants to be determined, which cannot be zero at the same time. Namely,
V = a 1 ( 3 t p + α t + x x + p 2 α m ( p + α ) u u ) + a 2 x .
Applying (9) to first boundary condition x = 0 , u ( t , 0 ) = q ( t ) , and initial condition u ( 0 , x ) = 0 , we have
V ( x ) | x = 0 = 0 ,
V ( u q ( t ) ) | x = 0 , u ( t , 0 ) = q ( t ) = 0 ,
V ( t ) | t = 0 = 0 ,
V ( u ) | t = 0 , u ( 0 , x ) = 0 = 0 .
Equation (10) gives a 1 × 0 + a 2 × 1 = 0 , namely a 2 = 0 .
From Equation (11), we obtain
a 1 × ( 3 t p + α q t + p 2 α m ( p + α ) q ) + a 2 × 0 = 0 .
For nonzero a 1 , we have
3 t p + α q t + p 2 α m ( p + α ) q = 0 ,
which gives
q ( t ) = λ t p 2 α 3 m ,
where λ > 0 is an integral constant.
It is obvious that the initial conditions (12) and (13) are invariant under V with a 2 = 0 .
Using the second extension of V with a 2 = 0 , that is,
P r V ( 2 ) = V 1 + ( 1 m ) p ( m + 2 ) α 3 m u x u x + ( 1 2 m ) p 2 ( m + 1 ) α 3 m u x x u x x ,
it can be shown that P r V ( 2 ) leaves the remaining boundary conditions u x ( t , 0 ) = 0 , u x x ( t , 0 ) = 0 invariant.
Namely,
P r ( 2 ) ( u x ) | u x ( t , 0 ) = 0 = 0 , P r ( 2 ) ( u x x ) | u x x ( t , 0 ) = 0 = 0 ,
which are obviously established.
Hence the V 1 is the general form of infinitesimal generator admitted by Equation (7) with Equation (8). □
Remark 1.
If boundary condition q ( t ) λ t p 2 α 3 m , then Equation (7) with (8) is not invariant under V 1 . Hence, q ( t ) = λ t p 2 α 3 m is an indispensable condition.
Remark 2.
To obtain one-parameter ( ε ) transformation by V 1 , one solves the system of first order ODEs,
d t d ε = 3 t p + α , d x d ε = x , d u d ε = p 2 α m ( p + α ) u ,
with t = t , x = x , u = u at ε = 0 . Then the corresponding one-parameter Lie group of transformations is
t = t e 3 p + α , x = x e ε , u = u e p 2 α m ( p + α ) .
Under transformations (14), the equations are obtained as following
α u ( t ) α + a ( u ) m u x + ( t ) p u x x x = 0 , x > 0 , t > 0 ,
u ( 0 , x ) = 0 , x > 0 , u ( t , 0 ) = λ ( t ) p 2 α 3 m , t > 0 , u x ( t , 0 ) = 0 , t > 0 , u x x ( t , 0 ) = 0 , t > 0 ,
which we have checked.
Hence the admitted infinitesimal generator can be used to reduce Equation (7) with the initial and boundary conditions (8) to an ordinary differential equation. The characteristic equations corresponding to infinitesimal generator V = V 1 are
d t 3 t p + α = d x x = d u p 2 α m ( p + α ) u .
The similarity variable and similarity transformation take the following form
z = x t p + α 3 , u = ϕ ( z ) t p 2 α 3 m ,
respectively.
According to the similarity transformation u = ϕ ( z ) t p 2 α 3 m , z = x t p + α 3 , the initial and boundary values are reduced to
ϕ ( 0 ) = λ , ϕ ( 0 ) = 0 , ϕ ( 0 ) = 0 .
Theorem 3.
Under the similarity transformation u = ϕ ( z ) t p 2 α 3 m with z = x t p + α 3 , time fractional generalized KdV equations(7) with initial and boundary value(8) is reduced into fractional ODEs of the form
( P 3 p + α 1 + p 2 α 3 m α , α ϕ ) ( z ) + a ( ϕ ( z ) ) m ϕ ( z ) + ϕ ( z ) = 0 , z [ 0 , + ) , ϕ ( 0 ) = λ , ϕ ( 0 ) = 0 , ϕ ( 0 ) = 0 ,
with the Erdélyi-Kober fractional differentia operator [10,12]
( P β τ , α ϕ ) ( z ) : = j = 0 n 1 ( τ + j 1 β z d d z ) ( K β τ + α , n α ϕ ) ( z ) , z > 0 , β > 0 , α > 0 ,
n = [ α ] + 1 , ( α N ) , α , ( α = N ) ,
where
( K β τ , α ϕ ) ( z ) = 1 Γ ( α ) 1 [ ( u 1 ) α 1 u ( τ + α ) ϕ ( z u 1 β ) ] d u , ( α > 0 ) , ϕ ( z ) , ( α = 0 ) ,
is the Erd e ´ lyi-Kober fractional integral operator and λ > 0 is a constant, the prime denotes differentiation with respect to z.
Proof. 
Suppose n 1 < α < n , n = 1 , 2 , 3 , . Under the similarity transformation (15), Riemann-Liouville fractional derivative becomes
α u t α = n t n [ 1 Γ ( n α ) 0 t ( t s ) n α 1 s p 2 α 3 m ϕ ( x s p + α 3 ) d s ] .
Letting v = t s , we have d s = t v 2 d v . Then the above equation can be written as
α u t α = n t n 1 Γ ( n α ) 0 t s n α 1 ( t s 1 ) n α 1 s p 2 α 3 m ϕ ( x t α + p 3 s α + p 3 t ( α + p ) 3 ) d s = n t n 1 Γ ( n α ) 0 t s n α 1 + p 2 α 3 m t n α 1 + p 2 α 3 m t n α 1 + p 2 α 3 m ( t s 1 ) n α 1 ϕ ( x t α + p 3 ( t s ) α + p 3 ) d s = n t n t n α 1 + p 2 α 3 m Γ ( n α ) 0 t ( s t ) n α 1 + p 2 α 3 m ( t s 1 ) n α 1 ϕ ( x t α + p 3 ( t s ) α + p 3 ) d s = n t n t n α 1 + p 2 α 3 m Γ ( n α ) + 1 ( 1 v ) n α 1 + p 2 α 3 m ( v 1 ) n α 1 ϕ ( z v α + p 3 ) ( t v 2 ) d v = n t n t n α + p 2 α 3 m Γ ( n α ) 1 v ( n α + 1 + p 2 α 3 m ) ( v 1 ) n α 1 ϕ ( z v p + α 3 ) ) d v = n t n t n α + p 2 α 3 m ( K 3 p + α 1 + p 2 α 3 m , n α ϕ ) ( z ) .
Using the relation ( z = x t ( p + α ) 3 , φ C 1 ( 0 , + ) ) , we have
t t φ ( z ) = p + α 3 z d d z φ ( z ) .
Then, Equation (16) is simplified as follows
n t n t n α + p 2 α 3 m ( K 3 p + α 1 + p 2 α 3 m , n α ϕ ) ( z ) = n 1 t n 1 t ( t n α + p 2 α 3 m ( K 3 p + α 1 + p 2 α 3 m , n α ϕ ) ( z ) ) = n 1 t n 1 t n α 1 + p 2 α 3 m ( n α + p 2 α 3 m p + α 3 z d d z ) ( K 3 p + α 1 + p 2 α 3 m , n α ϕ ) ( z ) .
Repeating the similar procedure for n 1 times, we obtain
n t n t n α + p 2 α 3 m ( K 3 p + α 1 + p 2 α 3 m , n α ϕ ) ( z ) = = t α + p 2 α 3 m j = 0 n 1 ( 1 α + p 2 α 3 m + j p + α 3 z d d z ) ( K 3 p + α 1 + p 2 α 3 m , n α ϕ ) ( z ) .
By the definition of the Erdélyi-Kober fractional differential operator given in Theorem 3, the above equation can be written as
α u t α = t α + p 2 α 3 m ( P 3 p + α 1 + p 2 α 3 m + α , α ϕ ) ( z ) .
According to u = ϕ ( z ) t p 2 α 3 m , z = x t p + α 3 , we obtain
u x = t p + α 3 + p 2 α 3 m ϕ ( z ) , u x x = t 2 p + 2 α 3 + p 2 α 3 m ϕ ( z ) , u x x x = t ( p + α ) + p 2 α 3 m ϕ ( z ) .
Substituting (17) and (18) into (7), the time fractional generalized KdV equations (7) with initial and boundary value (8) are reduced into the fractional ordinary differential equation with the form
( P 3 p + α 1 + p 2 α 3 m α , α ϕ ) ( z ) + a ( ϕ ( z ) ) m ϕ ( z ) + ϕ ( z ) = 0 , z [ 0 , + ) ,
ϕ ( 0 ) = λ , ϕ ( 0 ) = 0 , ϕ ( 0 ) = 0 ,
where the prime denotes differentiation with respect to z. □

4.3. Explicit Power Series Solutions of Time Fractional Generalized KdV Equations with Initial and Boundary Value

Then, we investigate explicit power series solutions of Equation (19) with conditions (20) applying the power series method in Reference [12,20,22].
Suppose that the power series solutions for Equation (19) is of the form
ϕ ( z ) = n = 0 a n z n .
From Equation (21), we can obtain
ϕ ( z ) = n = 0 ( n + 1 ) a n + 1 z n , ϕ ( z ) = n = 0 ( n + 1 ) ( n + 2 ) a n + 2 z n , ϕ ( z ) = n = 0 ( n + 1 ) ( n + 2 ) ( n + 3 ) a n + 3 z n .
Here, we only consider the following cases:
Case 1. When m = 1 in Equation (19), Equation (19) is written as
( P 3 p + α 1 + p 2 α 3 α , α ϕ ) ( z ) + a ( ϕ ( z ) ) ϕ ( z ) + ϕ ( z ) = 0 .
After substituting (21) and (22) into Equation (23), we obtain
n = 0 Γ ( 1 α + p 2 α 3 n ( p + α ) 3 ) Γ ( 1 + p 2 α 3 n ( p + α ) 3 ) a n z n + a n = 0 a n z n n = 0 ( n + 1 ) a n + 1 z n + n = 0 ( n + 1 ) ( n + 2 ) ( n + 3 ) a n + 3 z n = 0 .
When n = 0 in Equation (24), we have
a 3 = 1 6 Γ ( 1 α + p 2 α 3 ) Γ ( 1 + p 2 α 3 ) a 0 + a a 0 a 1 .
When n 1 in Equation (24), comparing coefficient of z n we have
a n + 3 = ( n + 3 ) ! n ! Γ ( 1 α + p 2 α 3 n ( p + α ) 3 ) Γ ( 1 + p 2 α 3 n ( p + α ) 3 ) a n + a k = 0 n ( n + 1 k ) a k a n + 1 k .
Therefore, according to (25) and (26) the power series solution for Equation (23) is of form as:
ϕ ( z ) = a 0 + a 1 z + a 2 z 2 1 6 Γ ( 1 α + p 2 α 3 ) Γ ( 1 + p 2 α 3 ) a 0 + a a 0 a 1 z 3 n = 1 ( n + 3 ) ! n ! Γ ( 1 α + p 2 α 3 n ( p + α ) 3 ) Γ ( 1 + p 2 α 3 n ( p + α ) 3 ) a n + a k = 0 n ( n + 1 k ) a k a n + 1 k z n + 3 .
Then substituting the initial conditions (20) into Equation (27), we obtain
a 0 = λ , a 1 = a 2 = 0 , a 3 = 1 6 Γ ( 1 α + p 2 α 3 ) Γ ( 1 + p 2 α 3 ) λ .
Hence, when m = 1 , applying the transformation of the Lie group, the explicit power series solutions of the time fractional generalized KdV equations with initial and boundary (Equations (7) and (8)) can be obtained as follows
u ( x , t ) = 1 1 6 Γ ( 1 α + p 2 α 3 ) Γ ( 1 + p 2 α 3 ) x 3 t ( p + α ) λ t p 2 α 3 n = 1 ( n + 3 ) ! n ! Γ ( 1 α + p 2 α 3 n ( p + α ) 3 ) Γ ( 1 + p 2 α 3 n ( p + α ) 3 ) a n + a k = 0 n ( n + 1 k ) a k a n + 1 k x n + 3 t p 2 α ( p + α ) ( n + 3 ) 3 .
Case 2. When m = 2 in Equation (19), Equation (19) is written as
( P 3 p + α 1 + p 2 α 6 α , α ϕ ) ( z ) + a ( ϕ ( z ) ) 2 ϕ ( z ) + ϕ ( z ) = 0 .
After substituting (21) and (22) into equation Equation (29), we obtain
n = 0 Γ ( 1 α + p 2 α 6 n ( p + α ) 3 ) Γ ( 1 + p 2 α 6 n ( p + α ) 3 ) a n z n + a n = 0 a n z n n = 0 a n z n n = 0 ( n + 1 ) a n + 1 z n + n = 0 ( n + 1 ) ( n + 2 ) ( n + 3 ) a n + 3 z n = 0 .
When n = 0 in Equation (30), we have
a 3 = 1 6 Γ ( 1 α + p 2 α 6 ) Γ ( 1 + p 2 α 6 ) a 0 + a a 0 2 a 1 .
When n 1 in Equation (30), comparing coefficient of z n we have
a n + 3 = ( n + 3 ) ! n ! Γ ( 1 α + p 2 α 6 n ( p + α ) 3 ) Γ ( 1 + p 2 α 6 n ( p + α ) 3 ) a n + a k = 0 n j = 0 k ( n + 1 k ) a j a k j a n + 1 k .
Therefore, according to (31) and (32), the power series solution for Equation (29) is of the form:
ϕ ( z ) = a 0 + a 1 z + a 2 z 2 1 6 Γ ( 1 α + p 2 α 6 ) Γ ( 1 + p 2 α 6 ) a 0 + a a 0 2 a 1 z 3 n = 1 ( n + 3 ) ! n ! Γ ( 1 α + p 2 α 6 n ( p + α ) 3 ) Γ ( 1 + p 2 α 6 n ( p + α ) 3 ) a n + a k = 0 n j = 0 k ( n + 1 k ) a j a k j a n + 1 k z n + 3 .
Then, substituting the initial conditions (20) into Equation (33), we obtain
a 0 = λ , a 1 = a 2 = 0 , a 3 = 1 6 Γ ( 1 α + p 2 α 6 ) Γ ( 1 + p 2 α 6 ) λ .
Hence, when m = 2 , applying the transformation of the Lie group, the explicit power series solutions of the time fractional generalized KdV equations with initial and boundary (Equations (7) and (8)) can be obtained as follows
u ( x , t ) = 1 1 6 Γ ( 1 α + p 2 α 6 ) Γ ( 1 + p 2 α 6 ) x 3 t ( p + α ) λ t p 2 α 3 n = 1 ( n + 3 ) ! n ! Γ ( 1 α + p 2 α 6 n ( p + α ) 3 ) Γ ( 1 + p 2 α 6 n ( p + α ) 3 ) a n + a k = 0 n j = 0 k ( n + 1 k ) a j a k j a n + 1 k x n + 3 t p 2 α ( p + α ) ( n + 3 ) 3 .
Case 3. When m = 3 in Equation (19), Equation (19) is written as
( P 3 p + α 1 + p 2 α 9 α , α ϕ ) ( z ) + a ( ϕ ( z ) ) 3 ϕ ( z ) + ϕ ( z ) = 0 .
After substituting (21) and (22) into equation Equation (35), we obtain
n = 0 Γ ( 1 α + p 2 α 9 n ( p + α ) 3 ) Γ ( 1 + p 2 α 9 n ( p + α ) 3 ) a n z n + a n = 0 a n z n n = 0 a n z n n = 0 a n z n n = 0 ( n + 1 ) a n + 1 z n + n = 0 ( n + 1 ) ( n + 2 ) ( n + 3 ) a n + 3 z n = 0 .
When n = 0 in Equation (36), we have
a 3 = 1 6 Γ ( 1 α + p 2 α 9 ) Γ ( 1 + p 2 α 9 ) a 0 + a a 0 3 a 1 .
When n 1 in Equation (36), comparing coefficient of z n , we have
a n + 3 = ( n + 3 ) ! n ! Γ ( 1 α + p 2 α 9 n ( p + α ) 3 ) Γ ( 1 + p 2 α 9 n ( p + α ) 3 ) a n + a k = 0 n j = 0 k l = 0 j ( n + 1 k ) a l a j l k a k j a n + 1 k .
Therefore, according to (37) and (38), the power series solution for Equation (35) is of form as:
ϕ ( z ) = a 0 + a 1 z + a 2 z 2 1 6 Γ ( 1 α + p 2 α 3 ) Γ ( 1 + p 2 α 3 ) a 0 + a a 0 a 1 z 3 n = 1 ( n + 3 ) ! n ! Γ ( 1 α + p 2 α 9 n ( p + α ) 3 ) Γ ( 1 + p 2 α 9 n ( p + α ) 3 ) a n + a k = 0 n j = 0 k l = 0 j ( n + 1 k ) a l a j l k a k j a n + 1 k z n + 3 .
Then substituting the initial conditions (20) into Equation (39), we obtain
a 0 = λ , a 1 = a 2 = 0 , a 3 = 1 6 Γ ( 1 α + p 2 α 9 ) Γ ( 1 + p 2 α 9 ) λ .
Hence, when m = 3 , applying the transformation of the Lie group, the explicit power series solutions of the time fractional generalized KdV equations with initial and boundary (Equations (7) and (8)) can be obtained as follows
u ( x , t ) = 1 1 6 Γ ( 1 α + p 2 α 9 ) Γ ( 1 + p 2 α 9 ) x 3 t ( p + α ) λ t p 2 α 3 n = 1 ( n + 3 ) ! n ! ( Γ ( 1 α + p 2 α 9 n ( p + α ) 3 ) Γ ( 1 + p 2 α 9 n ( p + α ) 3 ) a n + a k = 0 n j = 0 k l = 0 j ( n + 1 k ) a l a j l k a k j a n + 1 k ) x n + 3 t p 2 α ( p + α ) ( n + 3 ) 3 .

5. Conclusions

In this paper, the method of classical Lie symmetry analysis can be successfully extended to the time fractional partial differential equation with initial and boundary values. Then time fractional generalized KdV equations have been transformed into fractional ordinary differential equations with initial value. Then we utilize the power series method to find exact solutions of the reduced equation with the Erdélyi-Kober fractional differential operator. The exact solutions of the time fractional partial differential equation with initial and boundary value are given. Hence, it confirms that our work is feasible.

Author Contributions

These authors contributed equally to this work.

Funding

This work was supported by the Natural Science Foundation of China (NSFC) under grant (11871393), Shaanxi Provincial Key Research and Development Program (2018KW-004), International Science and Technology Cooperation Program of Shaanxi Key Research & Development Plan under grant 2019KWZ-08.

Acknowledgments

The authors are very grateful to the anonymous referees for their helpful comments and valuable suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bagley, R.L.; Torvik, P.J. A theoretical basis for the application of fractional calculus to viscoelasticity. J. Rheol. 1983, 27, 201–210. [Google Scholar] [CrossRef]
  2. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and applications of fractional differential equations. N.-Holl. Math. Stud. 2006, 204, 2453–2461. [Google Scholar]
  3. Miller, K.S.; Ross, B. An introduction to the fractional calculus and fractional order differential equations. Kybernetes 1993, 3, 66–67. [Google Scholar]
  4. Carpinteri, A.; Mainardi, F. (Eds.) Fractals and Fractional Calculus in Continuum Mechanics; Springer: Belin, Gemany, 1997; pp. 223–276. [Google Scholar]
  5. Diethelm, K. The analysis of fractional differential equations. An application-oriented exposition using differential operators of caputo type. Lect. Notes Math. Vol. 2010, 265, 229–248. [Google Scholar]
  6. El-Ajou, A.; Arqub, O.A.; Momani, S. Approximate analytical solution of the nonlinear fractional KdV–Burgers equation: A new iterative algorithm. J. Comput. Phys. 2015, 293, 81–95. [Google Scholar] [CrossRef]
  7. Su, W.H.; Baleanu, D.; Yang, X.J.; Jafari, H. Damped wave equation and dissipative wave equation in fractal strings within the local fractional variational iteration method. Fixed Point Theory Appl. 2013, 2013, 89. [Google Scholar] [CrossRef] [Green Version]
  8. Sakar, M.G.; Uludag, F.; Erdogan, F. Numerical solution of time-fractional nonlinear PDEs with proportional delays by homotopy perturbation method. Appl. Math. Model. 2016, 40, 6639–6649. [Google Scholar] [CrossRef]
  9. Ghany, H.A.; Hyder, A.A. Exact solutions for the wick-type stochastic time-fractional kdv equations. Kuwait J. Sci. 2014, 41, 75–84. [Google Scholar]
  10. Sahadevan, R.; Bakkyaraj, T. Invariant analysis of time fractional generalized Burgers and Korteweg-de Vries equations. J. Math. Anal. Appl. 2012, 393, 341–347. [Google Scholar] [CrossRef]
  11. Sahoo, S.; Garai, G.; Ray, S.S. Lie symmetry analysis for similarity reduction and exact solutions of modified KdV–Zakharov–Kuznetsov equation. Nonlinear Dyn. 2017, 87, 1995–2000. [Google Scholar] [CrossRef]
  12. Inc, M.; Yusuf, A.; Aliyu, A.I.; Baleanu, D. Time-fractional Cahn–Allen and time-fractional Klein–Gordon equations: Lie symmetry analysis, explicit solutions and convergence analysis. Phys. A Stat. Mech. Appl. 2018, 493, 94–106. [Google Scholar] [CrossRef]
  13. Chen, C.; Jiang, Y.L. Lie Group Analysis and Invariant Solutions for Nonlinear Time-Fractional Diffusion-Convection Equations. Commun. Theor. Phys. 2017, 68, 295–300. [Google Scholar] [CrossRef]
  14. Fung, M.K. KdV equation as an Euler-Poincare equation. Chin. J. Phys. 1997, 35, 789–796. [Google Scholar]
  15. Korteweg, D.D.J.; de Vries, D.G. XLI, on the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves. Philos. Mag. 2009, 39, 422–443. [Google Scholar] [CrossRef]
  16. Vaneeva, O.O.; Papanicolaou, N.C.; Christou, M.A.; Sophocleous, C. Numerical solutions of boundary value problems for variable coefficient generalized KdV equations using Lie symmetries. Commun. Nonlinear Sci. Numer. Simul. 2013, 19, 3074–3085. [Google Scholar] [CrossRef]
  17. Bluman, G.; Kumei, S. Symmetries and Differential Equations; Springer Science and Business Media: Belin, Germany, 2013. [Google Scholar]
  18. Gazizov, R.K.; Kasatkin, A.A.; Lukashchuk, S.Y. Symmetry properties of fractional diffusion equations. Phys. Scr. 2009, 2009, 014016. [Google Scholar] [CrossRef]
  19. Gazizov, R.K.; Kasatkin, A.A.; Lukashchuk, S.Y. Group-invariant solutions of fractional differential equations. Nonlinear Sci. Complex. 2011, 51–59. [Google Scholar] [CrossRef]
  20. Baleanu, D.; Inc, M.; Yusuf, A.; Aliyu, A.I. Lie symmetry analysis, exact solutions and conservation laws for the time fractional Caudrey-Dodd-Gibbon-Sawada-Kotera equation. Commun. Nonlinear Sci. Numer. Simul. 2018, 59, 222–234. [Google Scholar] [CrossRef]
  21. Chen, C.; Jiang, Y.L. Lie group analysis method for two classes of fractional partial differential equations. Commun. Nonlinear Sci. Numer. Simul. 2015, 26, 24–35. [Google Scholar] [CrossRef]
  22. Baleanu, D.; Inc, M.; Yusuf, A.; Aliyu, A.I. Lie symmetry analysis, conservation laws and explicit solutions for the time fractional Rosenau-Haynam equation. Waves Random Complex Media 2017, 27, 308–324. [Google Scholar]
  23. Ovsiannikov, L.V. Group Analysis of Differential Equations; Acadmic Press: New York, NY, USA, 1982. [Google Scholar]
  24. Ibragimov, N.H. CRC Handbook of Lie Group Analysis of Differential Equations; CRC Press: Boca Raton, FL, USA, 1995. [Google Scholar]
  25. Olver, P.J. Applications of Lie Groups to Differential Equations; Springer Science and Business Media: Berlin, Germany, 2000. [Google Scholar]

Share and Cite

MDPI and ACS Style

Chen, C.; Jiang, Y.-L.; Wang, X.-T. Lie Symmetry Analysis of the Time Fractional Generalized KdV Equations with Variable Coefficients. Symmetry 2019, 11, 1281. https://doi.org/10.3390/sym11101281

AMA Style

Chen C, Jiang Y-L, Wang X-T. Lie Symmetry Analysis of the Time Fractional Generalized KdV Equations with Variable Coefficients. Symmetry. 2019; 11(10):1281. https://doi.org/10.3390/sym11101281

Chicago/Turabian Style

Chen, Cheng, Yao-Lin Jiang, and Xiao-Tian Wang. 2019. "Lie Symmetry Analysis of the Time Fractional Generalized KdV Equations with Variable Coefficients" Symmetry 11, no. 10: 1281. https://doi.org/10.3390/sym11101281

APA Style

Chen, C., Jiang, Y. -L., & Wang, X. -T. (2019). Lie Symmetry Analysis of the Time Fractional Generalized KdV Equations with Variable Coefficients. Symmetry, 11(10), 1281. https://doi.org/10.3390/sym11101281

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

Article Metrics

Back to TopTop