Next Article in Journal
A Proposal about the Meaning of Scale, Scope and Resolution in the Context of the Information Interpretation Process
Next Article in Special Issue
Lyapunov Functions to Caputo Fractional Neural Networks with Time-Varying Delays
Previous Article in Journal
Incentive Compatible Decision Making: Real Options with Adverse Incentives
Previous Article in Special Issue
Mild Solutions to the Cauchy Problem for Some Fractional Differential Equations with Delay
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exact Solutions to the Fractional Differential Equations with Mixed Partial Derivatives

1
School of Science, Wuhan University of Science and Technology, Wuhan 430065, Hubei, China
2
Hubei Province Key Laboratory of Systems Science in Metallurgical Process, Wuhan 430065, Hubei, China
*
Author to whom correspondence should be addressed.
Axioms 2018, 7(1), 10; https://doi.org/10.3390/axioms7010010
Submission received: 19 December 2017 / Revised: 4 February 2018 / Accepted: 8 February 2018 / Published: 11 February 2018
(This article belongs to the Special Issue Fractional Differential Equations)

Abstract

:
In this paper, the solvability of nonlinear fractional partial differential equations (FPDEs) with mixed partial derivatives is considered. The invariant subspace method is generalized and is then used to derive exact solutions to the nonlinear FPDEs. Some examples are solved to illustrate the effectiveness and applicability of the method.

1. Introduction

In recent years, fractional order calculus has been one of the most rapidly developing areas of mathematical analysis. In fact, a natural phenomenon may depend not only on the time instant but also on the previous time history, which can be successfully modeled by fractional calculus. Fractional-order differential equations are naturally related to systems with memory, as fractional derivatives are usually nonlocal operators. Thus fractional differential equations (FDEs) play an important role because of their application in various fields of science, such as mathematics, physics, chemistry, optimal control theory, finance, biology, engineering and so on [1,2,3,4,5,6,7].
It is of importance to find efficient methods for solving FDEs. More recently, much attention has been paid to the solutions of FDEs using various methods, such as the Adomian decomposition method (2005) [8], the first integral method (2014) [9], the Lie group theory method (2012, 2015) [10,11], the homotopy analysis method (2016) [12], the inverse differential operational method (2016) [13,14,15], the F-expansion method (2017) [16], M-Wright transforms (2017) [17], exponential differential operators (2017, 2018) [18,19], and so on. In reality, the finding of exact solutions of the FDEs is hard work and remains a problem.
Recently, investigations have shown that a new method based on the invariant subspace provides an effective tool to find the exact solution of FDEs. This method was initially proposed by Galaktionov and Svirshchevskii (1995, 1996, 2007) [20,21,22]. The invariant subspace method was developed by Later Gazizov and Kasatkin (2013) [23], Harris and Garra (2013, 2014) [24,25], Sahadevan and Bakkyaraj (2015) [26], and Ouhadan and El Kinani (2015) [27].
In 2016, R. Sahadevan and P. Prakash [28] showed how the invariant subspace method could be extended to time fractional partial differential equations (FPDEs) and could construct their exact solutions.
α u t α = F ^ [ u ] , α > 0
where α t α ( · ) is a fractional time derivative in the Caputo sense, and F ^ [ u ] is a nonlinear differential operator of order k.
In 2016, S. Choudhary and V. Daftardar-Gejji [29] developed the invariant subspace method for deriving exact solutions of partial differential equations with fractional space and time derivatives.
( λ 0 α t α + λ 1 α + 1 t α + 1 + + λ m α + m t α + m ) f ( x , t ) = N ( x , f , β f x β , β + 1 f x β + 1 , , β + n f x β + n )
where N [ f ] is the linear/nonlinear differential operator; α + j f t α + j , j = 0 , 1 , , m and β + i f x β + i , i = 0 , 1 , , n are Caputo time derivatives and Caputo space derivatives, respectively; 0 < α , β 1 and λ i R .
In 2017, K.V. Zhukovsky [30] used the inverse differential operational method to obtain solutions for differential equations with mixed derivatives of physical problems.
Motivated by the above results, in this paper, we develop the invariant subspace method for finding exact solutions to some nonlinear partial differential equations with fractional-order mixed partial derivatives (including both fractional space derivatives and time derivatives).
( λ 0 α t α + λ 1 α + 1 t α + 1 + + λ m α + m 1 t α + m 1 ) f ( x , t ) = N ( x , f , β f x β , β + 1 f x β + 1 , , β + m 2 f x β + m 2 ) + μ α t α ( β f x β )
where f = f ( x , t ) , N [ f ] is a linear/nonlinear differential operator; α + j f t α + j , j = 0 , 1 , , m 1 , m 1 N and β + i f x β + i , i = 0 , 1 , , m 2 , m 2 N are Caputo time derivatives and Caputo space derivatives, respectively; α t α ( β f x β ) is the Caputo mixed partial derivative of space and time; k 1 < α k 1 + 1 , k 2 < β k 2 + 1 , k 1 , k 2 N and λ i , μ R .
Using the invariant subspace method, the FPDEs are reduced to the systems of FDEs that can be solved by familiar analytical methods.
The rest of this paper is organized as follows. In Section 2, the preliminaries and notations are given. In Section 3, we develop the invariant subspace method for solving fractional space and time derivative nonlinear partial differential equations with fractional-order mixed derivatives. In Section 4, illustrative examples are given to explain the applicability of the method. Initial value problems are considered. Finally in Section 5, we give conclusions.

2. Preliminaries and Notation

In this section, we recall some standard definitions and notation.
Definition 1.
(See [7]) The Riemann–Liouville fractional integral of order α and function f is defined as
I α f ( t ) = 1 Γ ( α ) 0 t f ( x ) ( t x ) 1 α d x , t > 0
Definition 2.
(See [7]) The Caputo fractional derivative of order α and function f is defined as
d α f ( t ) d t α = I n α D n f ( t ) = 1 Γ ( n α ) 0 t f ( n ) ( x ) ( t x ) α n + 1 d x , n 1 < α < n f ( n ) ( t ) , α = n , n N
The Riemann–Liouville fractional integral and the Caputo fractional derivative satisfy the following properties [3]:
I α t β = Γ ( β + 1 ) Γ ( β + α + 1 ) t β + α , α > 0 , β > 1 , t > 0
d α t β d t α = 0 , α = n , β { 0 , 1 , 2 , , n 1 } Γ ( β + 1 ) t β α Γ ( β α + 1 ) , α = n , β N , a n d   β n ; o r   β N , a n d   β > n 1
I α ( d α f ( t ) d t α ) = f ( t ) k = 0 n 1 d k f ( 0 ) d t k t k k ! , n 1 < α < n , t > 0
Definition 3.
(See [7]) A two-parametric Mittag–Leffler function is defined as
E α , β ( z ) = k = 0 z k Γ ( k α + β ) , α , β C , R ( α ) , R ( β ) > 0
noting that E α , 1 ( z ) = E α ( z ) .
Derivatives of the Mittag–Leffler function are given as
E α , β ( n ) ( z ) = d n d z n E α , β ( z ) = k = 0 ( k + n ) ! z k k ! Γ ( α k + α n + β ) , n = 0 , 1 , 2 ,
d γ d t γ ( t β 1 E α , β ( a t α ) ) = t β γ 1 E α , β γ ( a t α ) , α , γ > 0 , a R
d α d t α ( E α ( a t α ) ) = a E α ( a t α ) , α > 0 , a R
The Laplace transform of the α th order Caputo derivative is
T { d α f ( t ) d t α ; s } = s α f ˜ ( s ) k = 0 n 1 s α k 1 f ( k ) ( 0 ) , n 1 < α < n , n N , R ( s ) > 0
where
f ˜ ( s ) = T { f ( t ) ; s } = 0 e s t f ( t ) d t , s R
The Laplace transform of the function t α n + β 1 E α , β ( n ) ( ± a t α ) is as follows [9]:
T { t α n + β 1 E α , β ( n ) ( ± a t α ) ; s } = n ! s α β ( s α a ) n + 1 , R ( s ) > | a | 1 α , n = 0 , 1 , 2 ,
We let I n be the n-dimensional linear space over R. It is spanned by n linearly independent functions φ 0 ( x ) , φ 1 ( x ) , , φ n 1 ( x ) :
I n = L { φ 0 ( x ) , φ 1 ( x ) , , φ n 1 ( x ) } = { i = 0 n 1 k i φ i ( x ) | k i R , i = 0 , 1 , , n 1 }
We let M be a differential operator; if M [ f ] I n , f I n , then a finite-dimensional linear space I n is invariant with respect to a differential operator M.

3. Invariant Subspace Method; Fractional Partial Differential Equations with Fractional-Order Mixed Partial Derivative

The FPDE with fractional-order mixed partial derivative is as follows:
( λ 0 α t α + λ 1 α + 1 t α + 1 + + λ m 1 α + m 1 t α + m 1 ) f = N [ f ] + μ α t α ( β x β f )
where
f = f ( x , t ) , N [ f ] = N ( x , f , β x β f , β + 1 x β + 1 f , , β + m 2 x β + m 2 f )
Here, α + j f t α + j , j = 0 , 1 , , m 1 ; m 1 N and β + i f x β + i , i = 0 , 1 , , m 2 ; m 2 N are Caputo time derivatives and Caputo space derivatives respectively; α t α ( β x β f ) is the Caputo mixed partial derivative of space and time; k 1 < α k 1 + 1 , k 2 < β k 2 + 1 , k 1 , k 2 N and λ i , μ R .
Theorem 1.
Suppose I n + 1 = L { φ 0 ( x ) , φ 1 ( x ) , , φ n ( x ) } is a finite-dimensional linear space, and it is invariant with respect to the operators N [ f ] and β x β f ; then FPDE (1) has an exact solution as follows:
f ( x , t ) = i = 0 n k i ( t ) φ i ( x )
where { k i ( t ) } satisfies the following system of FDEs:
j = 0 m 1 λ j d α + j k i ( t ) d t α + j μ d α ψ n + 1 + i ( k 0 ( t ) , k 1 ( t ) , , k n ( t ) ) d t α = ψ i ( k 0 ( t ) , k 1 ( t ) , , k n ( t ) )
Here i = 0 , 1 , , n , { ψ 0 , ψ 1 , , ψ n } are the expansion coefficients of N [ f ] with respect to { φ 0 ( x ) , φ 1 ( x ) , , φ n ( x ) } ; { ψ n + 1 , ψ n + 2 , , ψ 2 n + 1 } are the expansion coefficients of β x β f with respect to { φ 0 ( x ) , φ 1 ( x ) , , φ n ( x ) } .
Proof. 
Using Equation (2) and the linearity of Caputo fractional derivatives, we obtain
j = 0 m 1 λ j α + j f ( x , t ) t α + j = j = 0 m 1 λ j α + j t α + j ( i = 0 n k i ( t ) φ i ( x ) ) = i = 0 n ( j = 0 m 1 λ j d α + j k i ( t ) d t α + j ) φ i ( x )
Further, as I n + 1 is an invariant space under the operator N [ f ] and β x β f , there exist 2 n + 2 functions ψ 0 , ψ 1 , , ψ n ; ψ n + 1 , ψ n + 2 , , ψ 2 n + 1 such that
N ( i = 0 n k i ( t ) φ i ( x ) ) = i = 0 n ψ i ( k 0 ( t ) , k 1 ( t ) , , k n ( t ) ) φ i ( x )
β x β f ( x , t ) = i = 0 n ψ n + 1 + i ( k 0 ( t ) , k 1 ( t ) , , k n ( t ) ) φ i ( x )
where { ψ 0 , ψ 1 , , ψ n } are the expansion coefficients of N [ f ] with respect to { φ 0 ( x ) , φ 1 ( x ) , , φ n ( x ) } ; { ψ n + 1 , ψ n + 2 , , ψ 2 n + 1 } are the expansion coefficients of β x β f with respect to { φ 0 ( x ) , φ 1 ( x ) , , φ n ( x ) } .
In view of Equations (2), (5) and (6),
N [ f ( x , t ) ] + μ α t α ( β x β f ( x , t ) ) = i = 0 n ψ i ( k 0 ( t ) , k 1 ( t ) , , k n ( t ) ) φ i ( x ) + μ α t α ( i = 0 n ψ n + 1 + i ( k 0 ( t ) , k 1 ( t ) , , k n ( t ) ) φ i ( x ) ) = i = 0 n ψ i ( k 0 ( t ) , k 1 ( t ) , , k n ( t ) ) φ i ( x ) + μ ( i = 0 n d α ψ n + 1 + i ( k 0 ( t ) , k 1 ( t ) , , k n ( t ) ) d t α φ i ( x ) ) = i = 0 n ( ψ i ( k 0 ( t ) , k 1 ( t ) , , k n ( t ) ) φ i ( x ) + μ d α ψ n + 1 + i ( k 0 ( t ) , k 1 ( t ) , , k n ( t ) ) d t α φ i ( x ) ) )
Equations (4) and (7) are substituted in Equation (1) to obtain
i = 0 n ( j = 0 m 1 λ j d α + j k i ( t ) d t α + j ψ i ( k 0 ( t ) , k 1 ( t ) , , k n ( t ) ) μ d α ψ n + 1 + i ( k 0 ( t ) , k 1 ( t ) , , k n ( t ) ) d t α ) φ i ( x ) = 0
Using Equation (8) and the fact that φ 0 ( x ) , φ 1 ( x ) , , φ n ( x ) are linearly independent, we have the system of FDEs that follows:
j = 0 m 1 λ j d α + j k i ( t ) d t α + j μ d α ψ n + 1 + i ( k 0 ( t ) , k 1 ( t ) , , k n ( t ) ) d t α = ψ i ( k 0 ( t ) , k 1 ( t ) , , k n ( t ) )
where i = 0 , 1 , , n . ☐
If FPDE (1) satisfies the conditions of Theorem 1, then FPDE (1) has a particular solution given by Equation (2).
We consider the following FPDE:
( λ 1 α t α + λ 2 2 α t 2 α + + λ m 1 m 1 α t m 1 α ) f = N [ f ] + μ α t α ( β x β f ) = N [ x , f , β x β f , 2 β x 2 β f , , m 2 β x m 2 β f ] + μ α t α ( β x β f )
where f = f ( x , t ) , N [ f ] is a linear/nonlinear differential operator; j α f t j α , j = 1 , 2 , , m 1 ; m 1 N and i β f x i β , i = 1 , 2 , , m 2 ; m 2 N are Caputo time derivatives and Caputo space derivatives, respectively; k 1 < α k 1 + 1 , k 2 < β k 2 + 1 , k 1 , k 2 N and λ i , μ R .
Theorem 2.
Suppose I n = L { φ 1 ( x ) , φ 2 ( x ) , , φ n ( x ) } is a finite-dimensional linear space, and it is invariant with respect to the operator N [ f ] and β x β f ; then FPDE (10) has an exact solution as follows:
f ( x , t ) = i = 1 n k i ( t ) φ i ( x )
where { k i ( t ) } satisfy the following system of FDEs:
j = 1 m 1 λ j d j α k i ( t ) d t j α μ d α ψ n + i ( k 1 ( t ) , k 2 ( t ) , , k n ( t ) ) d t α = ψ i ( k 1 ( t ) , k 2 ( t ) , , k n ( t ) )
Here, i = 1 , 2 , , n , { ψ 1 , ψ 2 , , ψ n } are the expansion coefficients of N [ f ] with respect to { φ 1 ( x ) , φ 2 ( x ) , , φ n ( x ) } ; { ψ n + 1 , ψ n + 2 , , ψ 2 n } are the expansion coefficients of β x β f with respect to { φ 1 ( x ) , φ 2 ( x ) , , φ n ( x ) } .
Proof. 
Using Equation (11) and the linearity of Caputo fractional derivatives, we obtain
j = 1 m 1 λ j j α f ( x , t ) t j α = j = 1 m 1 λ j j α t j α ( i = 1 n k i ( t ) φ i ( x ) ) = i = 1 n ( j = 1 m 1 λ j d j α k i ( t ) d t j α ) φ i ( x )
Further, as I n is an invariant space under the operator N [ f ] and β x β f , there exist 2 n functions ψ 1 , ψ 2 , , ψ n ; ψ n + 1 , ψ n + 2 , , ψ 2 n such that
N ( i = 1 n k i ( t ) φ i ( x ) ) = i = 1 n ψ i ( k 1 ( t ) , k 2 ( t ) , , k n ( t ) ) φ i ( x )
β x β f ( x , t ) = i = 1 n ψ n + i ( k 1 ( t ) , k 2 ( t ) , , k n ( t ) ) φ i ( x )
where { ψ 1 , ψ 2 , , ψ n } are the expansion coefficients of N [ f ] with respect to { φ 1 ( x ) , φ 2 ( x ) , , φ n ( x ) } ; { ψ n + 1 , ψ n + 2 , , ψ 2 n } are the expansion coefficients of β x β f with respect to { φ 1 ( x ) , φ 2 ( x ) , , φ n ( x ) } .
In view of Equations (11), (14) and (15),
N [ f ( x , t ) ] + μ α t α ( β x β f ( x , t ) ) = i = 1 n ( ψ i ( k 1 ( t ) , k 2 ( t ) , , k n ( t ) ) φ i ( x ) + μ d α ψ n + i ( k 1 ( t ) , k 2 ( t ) , , k n ( t ) ) d t α φ i ( x ) )
Equations (13) and (16) are substituted into Equation (10), to obtain
i = 1 n ( j = 1 m 1 λ j d j α k i ( t ) d t j α ψ i ( k 1 ( t ) , k 2 ( t ) , , k n ( t ) ) μ d α ψ n + i ( k 1 ( t ) , k 2 ( t ) , , k n ( t ) ) d t α ) φ i ( x ) = 0
Using Equation (17) and the fact that φ 1 ( x ) , φ 2 ( x ) , , φ n ( x ) are linearly independent, we have the system of FDEs as follows:
j = 1 m 1 λ j d j α k i ( t ) d t j α μ d α ψ n + i ( k 1 ( t ) , k 2 ( t ) , , k n ( t ) ) d t α = ψ i ( k 1 ( t ) , k 2 ( t ) , , k n ( t ) )
here i = 0 , 1 , , n . ☐
Remark: Theorems 1 and 2 in [27] are special cases of our results for μ = 0 .

4. Illustrative Examples

In this section, we give several examples to illustrate Theorems 1 and 2.
Example 1.
The fractional diffusion equation is as follows:
α f t α = C β f x β + μ α t α ( β x β f )
where C = constant.
Diffusion is a process in which molecules move around until they are evenly spread out in the area. For α > 1 , the phenomenon is referred to as super-diffusion, and for α = 1 , it is called normal diffusion, whereas α < 1 describes subdiffusion.
We consider two cases of Equation (19): case 1: α ( 0 , 1 ] , β ( 1 , 2 ] ; case 2: α ( 1 , 2 ] , β ( 1 , 2 ] .
Case 1: α ( 0 , 1 ] , β ( 1 , 2 ] .
The subspace I 2 = L { 1 , x β } is invariant under N [ f ] = C β f x β and β x β f as
N [ k 0 + k 1 x β ] = C k 1 Γ ( β + 1 ) I 2
β ( k 0 + k 1 x β ) x β = k 1 Γ ( β + 1 ) I 2
It follows from Theorem 1 applied to Equation (19) that Equation (19) has the exact solution that follows:
f ( x , t ) = k 0 ( t ) + k 1 ( t ) x β
where k 0 ( t ) and k 1 ( t ) satisfy the system of FDEs as follows:
d α k 0 ( t ) d t α d α k 1 ( t ) d t α μ C Γ ( β + 1 ) = k 1 ( t ) C Γ ( β + 1 )
d α k 1 ( t ) d t α = 0
Solving the above FDE (22), we obtain
k 1 ( t ) = b
Substituting Equation (23) into Equation (21), we obtain
d α k 0 ( t ) d t α = b C Γ ( β + 1 )
Then
k 0 ( t ) = a + b C Γ ( β + 1 ) Γ ( α + 1 ) t α
Substituting Equations (23) and (25) into Equation (19), we obtain Equation (19) with the solution as follows:
f ( x , t ) = a + b C Γ ( β + 1 ) Γ ( α + 1 ) t α + b x β
where a and b are arbitrary constants.
It is clearly verified that the subspace I 3 = L { 1 , x β , x 2 β } is invariant under N [ f ] = C β f x β and β x β f as
N [ k 0 + k 1 x β + k 2 x 2 β ] = C k 1 Γ ( β + 1 ) + C k 2 Γ ( 2 β + 1 ) Γ ( β + 1 ) x β I 3
β ( k 0 + k 1 x β + k 2 x 2 β ) x β = k 1 Γ ( β + 1 ) + k 2 Γ ( 2 β + 1 ) Γ ( β + 1 ) x β I 3
We let Equation (19) have the exact solution that follows:
f ( x , t ) = k 0 ( t ) + k 1 ( t ) x β + k 2 ( t ) x 2 β
where k 0 ( t ) , k 1 ( t ) and k 2 ( t ) satisfy the system of FDEs as follows:
d α k 0 ( t ) d t α μ d α k 1 ( t ) d t α Γ ( β + 1 ) = C Γ ( β + 1 ) k 1 ( t )
d α k 1 ( t ) d t α μ d α k 2 ( t ) d t α Γ ( 2 β + 1 ) Γ ( β + 1 ) = C Γ ( 2 β + 1 ) Γ ( β + 1 ) k 2 ( t )
d α k 2 ( t ) d t α = 0
Equation (30) implies that k 2 ( t ) = a 2 . Thus Equation (29) takes the form
d α k 1 ( t ) d t α = a 2 C Γ ( 2 β + 1 ) Γ ( β + 1 )
which has the following solution:
k 1 ( t ) = a 1 + a 2 C Γ ( 2 β + 1 ) Γ ( β + 1 ) t α
Similarly, Equation (28) yields
k 0 ( t ) = a 0 + a 1 C Γ ( β + 1 ) 2 + a 2 μ C Γ ( 2 β + 1 ) Γ ( α + 1 ) Γ ( β + 1 ) t α + a 2 C 2 Γ ( 2 β + 1 ) Γ ( 2 α + 1 ) t 2 α
Thus, Equation (19) has the following solution:
f ( x , t ) = ( a 0 + a 1 C Γ ( β + 1 ) 2 + a 2 μ C Γ ( 2 β + 1 ) Γ ( α + 1 ) Γ ( β + 1 ) t α + a 2 C 2 Γ ( 2 β + 1 ) Γ ( 2 α + 1 ) t 2 α ) + ( a 1 + a 2 C Γ ( 2 β + 1 ) Γ ( β + 1 ) t α ) x β + a 2 x 2 β
where a 0 , a 1 and a 2 are arbitrary constants.
It can be easily verified that I 2 = L { 1 , E β ( x β ) } is also an invariant subspace with respect to N [ f ] = C β f x β and β x β f , as
N [ k 0 + k 1 E β ( x β ) ] = C β x β ( k 0 + k 1 E β ( x β ) ) = C k 1 E β ( x β ) I 2
β ( k 0 + k 1 E β ( x β ) ) x β = k 1 E β ( x β ) I 2
We consider the exact solution of the form
f ( x , t ) = k 0 ( t ) + k 1 ( t ) E β ( x β )
where k 0 ( t ) and k 1 ( t ) satisfy the following system of FDEs:
d α k 0 ( t ) d t α = 0
( 1 μ C ) d α k 1 ( t ) d t α = C k 1 ( t )
Clearly, k 0 ( t ) = a . Solving Equation (33) with the Laplace transform method, we obtain the following:
If μ C 1 ,
k 1 ( t ) = b E α ( C 1 μ C t α )
Thus Equation (19) has the exact solution that follows:
f ( x , t ) = a + b E α ( C 1 μ C t α ) E β ( x β )
where a and b are arbitrary constants.
We find that Equations (26), (31) and (34) are distinct particular solutions of Equation (19) under distinct invariant subspaces. Subspace I n + 1 = L { 1 , x β , x 2 β , , x n β } , n N is invariant under N [ f ] = C β f x β and β x β f , as
N [ k 0 + k 1 x β + + k n x n β ] = C k 1 Γ ( β + 1 ) + C k 2 Γ ( 2 β + 1 ) Γ ( β + 1 ) x β + + C k n Γ ( n β + 1 ) Γ ( ( n 1 ) β + 1 ) x ( n 1 ) β I n + 1
Thus we obtain infinitely many invariant subspaces for Equation (19), which in turn yield infinitely many particular solutions.
Case 2: α ( 1 , 2 ] , β ( 1 , 2 ] .
Clearly, subspace I 3 = L { 1 , x β , x 2 β } is an invariant subspace under N [ f ] = C β f x β and β x β f , as
N [ k 0 + k 1 x β + k 2 x 2 β ] = C k 1 Γ ( β + 1 ) + C k 2 Γ ( 2 β + 1 ) Γ ( β + 1 ) x β I 3
β ( k 0 + k 1 x β + k 2 x 2 β ) x β = k 1 Γ ( β + 1 ) + k 2 Γ ( 2 β + 1 ) Γ ( β + 1 ) x β I 3
We look for the exact solution that follows:
f ( x , t ) = k 0 ( t ) + k 1 ( t ) x β + k 2 ( t ) x 2 β
where k 0 ( t ) , k 1 ( t ) and k 2 ( t ) are unknown functions to be determined; k 0 ( t ) , k 1 ( t ) and k 2 ( t ) satisfy the system of FDEs as follows:
d α k 0 ( t ) d t α μ d α k 1 ( t ) d t α Γ ( β + 1 ) = C Γ ( β + 1 ) k 1 ( t )
d α k 1 ( t ) d t α μ d α k 2 ( t ) d t α Γ ( 2 β + 1 ) Γ ( β + 1 ) = C Γ ( 2 β + 1 ) Γ ( β + 1 ) k 2 ( t )
d α k 2 ( t ) d t α = 0
Solving Equations (35)–(37), we obtain
k 2 ( t ) = d 1 + d 2 t
k 1 ( t ) = b 1 + b 2 t + C d 1 Γ ( 2 β + 1 ) Γ ( β + 1 ) Γ ( α + 1 ) t α + C d 2 Γ ( 2 β + 1 ) Γ ( β + 1 ) Γ ( α + 2 ) t α + 1
k 0 ( t ) = a 1 + a 2 t + C b 1 Γ ( β + 1 ) + μ C d 2 Γ ( 2 β + 1 ) Γ ( α + 1 ) t α + C b 2 Γ ( β + 1 ) + μ C d 2 Γ ( 2 β + 1 ) Γ ( α + 2 ) t α + 1 + C 2 d 1 Γ ( 2 β + 1 ) Γ ( 2 α + 1 ) t 2 α + C 2 d 2 Γ ( 2 β + 1 ) Γ ( 2 α + 2 ) t 2 α + 1
Then, we obatin the exact solution of Equation (19) as
f ( x , t ) = ( a 1 + a 2 t + C b 1 Γ ( β + 1 ) + μ C d 2 Γ ( 2 β + 1 ) Γ ( α + 1 ) t α + C b 2 Γ ( β + 1 ) + μ C d 2 Γ ( 2 β + 1 ) Γ ( α + 2 ) t α + 1 + C 2 d 1 Γ ( 2 β + 1 ) Γ ( 2 α + 1 ) t 2 α + C 2 d 2 Γ ( 2 β + 1 ) Γ ( 2 α + 2 ) t 2 α + 1 ) + ( b 1 + b 2 t + C d 1 Γ ( 2 β + 1 ) Γ ( β + 1 ) Γ ( α + 1 ) t α + C d 2 Γ ( 2 β + 1 ) Γ ( β + 1 ) Γ ( α + 2 ) t α + 1 ) x β + ( d 1 + d 2 t ) x 2 β
where a 1 , a 2 , b 1 , b 2 , d 1 and d 2 are arbitrary constants.
When α and β are other numbers, we can similarly obtain the exact solution of Equation (19).
Next, we find the closed-form solutions of FPDEs satisfying initial conditions using the invariant subspace method.
Example 2.
We have the following FPDE with the initial condition as follows:
α f t α = ( β f x β ) 2 f ( β f x β ) + μ α t α ( β x β f ) , α , β ( 0 , 1 ]
f ( x , 0 ) = 3 + 5 2 E β ( x β )
The subspace I 2 = L { 1 , E β ( x β ) } is invariant under N [ f ] = ( β f x β ) 2 f ( β f x β ) and β x β f , as
N [ k 0 + k 1 E β ( x β ) ] = ( k 1 E β ( x β ) ) 2 ( k 0 + k 1 E β ( x β ) ) k 1 E β ( x β ) = k 0 k 1 E β ( x β ) I 2
β ( k 0 + k 1 E β ( x β ) ) x β = k 1 E β ( x β ) I 2
We consider the exact solution that follows:
f ( x , t ) = k 0 ( t ) + k 1 ( t ) E β ( x β )
where k 0 ( t ) and k 1 ( t ) are unknown functions to be determined.
By substituting Equation (40) into Equation (38) and equating coefficients of different powers of x, we obtain the following system of FDEs:
d α k 0 ( t ) d t α = 0
( 1 μ ) d α k 1 ( t ) d t α = k 0 ( t ) k 1 ( t )
We obtain k 0 ( t ) = a , and Equation (42) takes the following form:
If μ 1 ,
d α k 1 ( t ) d t α = a μ 1 k 1 ( t )
Then using the Laplace transform technique, we obtain
s α k ˜ 1 ( s ) s α 1 k 1 ( 0 ) = a μ 1 k ˜ 1 ( s )
k ˜ 1 ( s ) = s α 1 s α a μ 1 k 1 ( 0 )
Using the inverse Laplace transform, we obtain
k 1 ( t ) = k 1 ( 0 ) E α ( a μ 1 t α )
which leads to the exact solution of Equation (38) that follows:
f ( x , t ) = a + b E α ( a μ 1 t α ) E β ( x β )
where a and b are arbitrary constants.
Thus the exact solution of Equation (38) along with the initial condition of Equation (39) is
f ( x , t ) = 3 + 5 2 E α ( 3 μ 1 t α ) E β ( x β )
Example 3.
The fractional wave equation is used as an example to model the propagation of diffusive waves in viscoelastic solids. We considered the fractional wave equation with a constant absorption term as follows:
2 α f t 2 α = β x β ( f β f x β ) 1 + μ α t α ( β x β f ) , α , β ( 0 , 1 ]
f ( x , 0 ) = e + x β Γ ( β + 1 ) , f t ( x , 0 ) = 1 x β
Clearly, the subspace I 2 = L { 1 , x β } is invariant under N [ f ] = β x β ( f β f x β ) 1 and β x β f , as
N [ k 0 + k 1 x β ] = β x β ( ( k 0 + k 1 x β ) k 1 Γ ( β + 1 ) ) 1 = k 1 2 Γ 2 ( β + 1 ) 1 I 2
β ( k 0 + k 1 x β ) x β = k 1 Γ ( β + 1 ) I 2
By an application of Theorem 2, we know that Equation (43) has the exact solution as follows:
f ( x , t ) = k 0 ( t ) + k 1 ( t ) x β
where k 0 ( t ) and k 1 ( t ) satisfy the system of FDEs as follows:
d 2 α k 0 ( t ) d t 2 α μ d α k 1 ( t ) d t α Γ ( β + 1 ) = k 1 2 ( t ) Γ 2 ( β + 1 ) 1
d 2 α k 1 ( t ) d t 2 α = 0
Solving Equations (45) and (46) we obtain the following:
Case 1: when 0 < α 1 2 :
k 1 ( t ) = b 1
k 0 ( t ) = a 1 + b 1 2 Γ 2 ( β + 1 ) 1 Γ ( 2 α + 1 ) t 2 α
Thus Equation (43) has the exact solution that follows:
f ( x , t ) = ( a 1 + b 1 2 Γ 2 ( β + 1 ) 1 Γ ( 2 α + 1 ) t 2 α ) + b 1 x β
where a 1 and b 1 are arbitrary constants.
By the initial conditions of Equation (44), we obtain a 1 = e and b 1 = 1 Γ ( β + 1 ) .
Hence the exact solution of Equations (40) and (41) is
f ( x , t ) = e + x β Γ ( β + 1 )
Case 2: when 1 2 < α 1 :
k 1 ( t ) = b 1 + b 2 t
k 0 ( t ) = a 1 + a 2 t + μ b 2 Γ ( β + 1 ) Γ ( α + 2 ) t α + 1 + b 1 2 Γ 2 ( β + 1 ) 1 Γ ( 2 α + 1 ) t 2 α + 2 b 1 b 2 Γ 2 ( β + 1 ) Γ ( 2 α + 2 ) t 2 α + 1 + 2 b 2 2 Γ 2 ( β + 1 ) Γ ( 2 α + 3 ) t 2 α + 2
Thus Equation (43) has the exact solution that follows:
f ( x , t ) = ( a 1 + a 2 t + μ b 2 Γ ( β + 1 ) Γ ( α + 2 ) t α + 1 + b 1 2 Γ 2 ( β + 1 ) 1 Γ ( 2 α + 1 ) t 2 α + 2 b 1 b 2 Γ 2 ( β + 1 ) Γ ( 2 α + 2 ) t 2 α + 1 + 2 b 2 2 Γ 2 ( β + 1 ) Γ ( 2 α + 3 ) t 2 α + 2 ) + ( b 1 + b 2 t ) x β
where a 1 , a 2 , b 1 and b 2 are arbitrary constants.
Substituting the initial conditions of Equation (44), we obtain a 1 = e , a 2 = 1 , b 1 = 1 Γ ( β + 1 ) and b 2 = 1 .
Thus the exact solution of Equations (43) and (44) is
f ( x , t ) = ( e + t μ Γ ( β + 1 ) Γ ( α + 2 ) t α + 1 2 Γ ( β + 1 ) Γ ( 2 α + 2 ) t 2 α + 1 + 2 Γ 2 ( β + 1 ) Γ ( 2 α + 3 ) t 2 α + 2 ) + ( 1 Γ ( β + 1 ) t ) x β
We consider the following fractional generalization of the wave equation with a constant absorption term:
α + 1 f t α + 1 = β x β ( f β f x β ) 1 + μ α t α ( β x β f ) , α , β ( 0 , 1 ]
We know that the subspace I 2 = L { 1 , x β } is invariant from the above. In view of Theorem 1, Equation (47) has the exact solution that follows:
f ( x , t ) = k 0 ( t ) + k 1 ( t ) x β
where k 0 ( t ) and k 1 ( t ) satisfy the system of FDEs as follows:
d α + 1 k 0 ( t ) d t α + 1 μ d α k 1 ( t ) d t α Γ ( β + 1 ) = k 1 2 ( t ) Γ 2 ( β + 1 ) 1
d α + 1 k 1 ( t ) d t α + 1 = 0
Solving the system of FDEs (48) and (49), we obtain
k 1 ( t ) = b 1 + b 2 t
k 0 ( t ) = a 1 + a 2 t + b 1 2 Γ 2 ( β + 1 ) 1 Γ ( α + 2 ) t α + 1 + μ b 2 Γ ( β + 1 ) 2 t 2 + 2 b 1 b 2 Γ 2 ( β + 1 ) Γ ( α + 3 ) t α + 2 + 2 b 2 2 Γ 2 ( β + 1 ) Γ ( α + 4 ) t α + 3
Therefore Equation (47) has the exact solution that follows:
f ( x , t ) = ( a 1 + a 2 t + b 1 2 Γ 2 ( β + 1 ) 1 Γ ( α + 2 ) t α + 1 + μ b 2 Γ ( β + 1 ) 2 t 2 + 2 b 1 b 2 Γ 2 ( β + 1 ) Γ ( α + 3 ) t α + 2 + 2 b 2 2 Γ 2 ( β + 1 ) Γ ( α + 4 ) t α + 3 ) + ( b 1 + b 2 t ) x β
where a 1 , a 2 , b 1 and b 2 are arbitrary constants.
Example 4.
The Korteweg–de Vries (KdV) equation describes the evolution in time of long, unidirectional, nonlinear shallow water waves. We considered the fractional KdV equation that follows:
α f t α = β f x β ( β x β ( β x β ( f 2 2 ) ) ) + μ α t α ( β x β f ) , α , β ( 0 , 1 ]
I 3 = L { 1 , x β , x 2 β } is an invariant subspace under N [ f ] = β f x β ( β x β ( β x β ( f 2 2 ) ) ) and β x β f , as
N [ k 0 + k 1 x β + k 2 x 2 β ] = k 1 k 2 Γ ( 3 β + 1 ) + k 2 2 Γ ( 4 β + 1 ) 2 Γ ( β + 1 ) x β I 3
β ( k 0 + k 1 x β + k 2 x 2 β ) x β = k 1 Γ ( β + 1 ) + k 2 Γ ( 2 β + 1 ) Γ ( β + 1 ) x β I 3
We consider an exact solution that follows:
f ( x , t ) = k 0 ( t ) + k 1 ( t ) x β + k 2 ( t ) x 2 β
where k 0 ( t ) , k 1 ( t ) and k 2 ( t ) are unknown functions to be determined. It follows from Theorem 1 applied to Equation (47) that k 0 ( t ) , k 1 ( t ) and k 2 ( t ) satisfy the FDEs as follows:
d α k 0 ( t ) d t α μ d α k 1 ( t ) d t α Γ ( β + 1 ) = Γ ( 3 β + 1 ) k 1 ( t ) k 2 ( t )
d α k 1 ( t ) d t α μ d α k 2 ( t ) d t α Γ ( 2 β + 1 ) Γ ( β + 1 ) = Γ ( 4 β + 1 ) 2 Γ ( β + 1 ) k 2 2 ( t )
d α k 2 ( t ) d t α = 0
Solving Equation (53), we obtain k 2 ( t ) = c .
Hence Equation (52) has the form
d α k 1 ( t ) d t α = c 2 Γ ( 4 β + 1 ) 2 Γ ( β + 1 )
We obtain
k 1 ( t ) = b + c 2 Γ ( 4 β + 1 ) 2 Γ ( β + 1 ) Γ ( α + 1 ) t α
Similarly, we obtain
k 0 ( t ) = a + 2 b c Γ ( 3 β + 1 ) + μ c 2 Γ ( 4 β + 1 ) 2 Γ ( α + 1 ) t α + c 3 Γ ( 4 β + 1 ) Γ ( 3 β + 1 ) 2 Γ ( β + 1 ) Γ ( 2 α + 1 ) t 2 α
Thus the exact solution of Equation (50) is
f ( x , t ) = ( a + 2 b c Γ ( 3 β + 1 ) + μ c 2 Γ ( 4 β + 1 ) 2 Γ ( α + 1 ) t α + c 3 Γ ( 4 β + 1 ) Γ ( 3 β + 1 ) 2 Γ ( β + 1 ) Γ ( 2 α + 1 ) t 2 α ) + ( b + c 2 Γ ( 4 β + 1 ) 2 Γ ( β + 1 ) Γ ( α + 1 ) t α ) x β + c x 2 β
where a , b and c are arbitrary constants.
Example 5.
The fractional version of the nonlinear heat equation is as follows:
α f t α = β f x β ( f β x β ) + μ α t α ( β x β f ) , α , β ( 0 , 1 ]
Clearly, the subspace I 2 = L { 1 , x β } is invariant under N [ f ] = β f x β ( f β x β ) and β x β f , as
N [ k 0 + k 1 x β ] = β x β ( ( k 0 + k 1 x β ) k 1 Γ ( β + 1 ) ) = k 1 2 Γ 2 ( β + 1 ) I 2
β ( k 0 + k 1 x β ) x β = k 1 Γ ( β + 1 ) I 2
It follows from Theorem 1 that we consider the exact solution of Equation (54) as follows:
f ( x , t ) = k 0 ( t ) + k 1 ( t ) x β
such that
d α k 0 ( t ) d t α μ d α k 1 ( t ) d t α Γ ( β + 1 ) = k 1 2 ( t ) Γ 2 ( β + 1 )
d α k 1 ( t ) d t α = 0
Solving Equations (55) and (56), we obtain
k 1 ( t ) = b
k 0 ( t ) = a + b 2 Γ 2 ( β + 1 ) Γ ( α + 1 ) t α
We obtain an exact solution as follows:
f ( x , t ) = ( a + b 2 Γ 2 ( β + 1 ) Γ ( α + 1 ) t α ) + b x β
where a and b are arbitrary constants.
Next, we consider the integer-order differential equations in [30]. We can obtain some new different solutions using the invariant subspace method.
Example 6.
The modified hyperbolic heat conduction equation with the mixed derivative term is as follows ([30]):
2 f t 2 = 2 f x 2 + κ 2 f ε 2 t x f
f ( x , 0 ) = g ( x ) , f ( x , ) <
where ε , κ = const.
Clearly, the subspace I 2 = L { 1 , x } is an invariant subspace under N [ f ] = 2 f x 2 + κ 2 f and x f , as
N [ k 0 + k 1 x ] = κ 2 k 0 + κ 2 k 1 x I 2
( k 0 + k 1 x ) x = k 1 I 2
We let the exact solution be as follows:
f ( x , t ) = k 0 ( t ) + k 1 ( t ) x
where k 0 ( t ) , k 1 ( t ) are unknown functions to be determined, and k 0 ( t ) and k 1 ( t ) satisfy the system of differential equations as follows:
d 2 k 0 ( t ) d t 2 + ε d k 1 ( t ) d t = κ 2 k 0 ( t )
d 2 k 1 ( t ) d t 2 = κ 2 k 1 ( t )
Solving Equation (60), we obtain
k 1 ( t ) = b 1 e κ t + b 2 e κ t
Hence Equation (59) has the form
d 2 k 0 ( t ) d t 2 κ 2 k 0 ( t ) = b 1 ε κ e κ t + b 2 ε κ e κ t
We obtain
k 0 ( t ) = a 1 e κ t + a 2 e κ t a 1 ε 2 t e κ t a 2 ε 2 t e κ t
Then, we obtain the exact solution of Equation (57) as
f ( x , t ) = a 1 e κ t + a 2 e κ t a 1 ε 2 t e κ t a 2 ε 2 t e κ t + ( b 1 e κ t + b 2 e κ t ) x
where a 1 , a 2 , b 1 and b 2 are arbitrary constants.
Substituting the conditions of Equation (58) into Equation (61), we obtain a 2 = b 2 = 0 and a 1 + b 1 x = g ( x ) .
When g ( x ) has linear dependence on x, Equations (57) and (58) have the partial solution
f ( x , t ) = a 1 e κ t a 1 ε 2 t e κ t + b 1 x e κ t
where a 1 + b 1 x = g ( x ) .
When g ( x ) is not linearly dependent on x, Equations (57) and (58) do not have the form of the solution given by Equation (61).
The subspace I n + 1 = L { 1 , x , x 2 , , x n } , n N is invariant under N [ f ] = 2 f x 2 + κ 2 f and x f , as
N [ k 0 + k 1 x + + k n x n ] = ( κ 2 k 0 + 2 k 2 ) + + ( κ 2 k n 2 + n ( n 1 ) k n ) x n 2 + κ 2 k n 1 x n 1 + κ 2 k n x n I n + 1
( k 0 + k 1 x + + k n x n ) x = k 1 + k 2 x + + k n 1 x n 1 I n + 1
Thus we obtain infinitely many invariant subspaces for Equation (57), which in turn yield infinitely many solutions. If g ( x ) is a polynomial, we can obtain an exact solution of Equations (57) and (58).
Example 7.
The Fokker–Planck equation is the following ([30]):
2 f t 2 = α 2 f x 2 + β x f x ε 2 t x f
where α , β , κ = const.
Clearly, the subspace I 2 = L { 1 , x } is an invariant subspace under N [ f ] = α 2 f x 2 + β x f x and x f , as
N [ k 0 + k 1 x ] = β k 1 x I 2
( k 0 + k 1 x ) x = k 1 I 2
We suppose the exact solution that follows:
f ( x , t ) = k 0 ( t ) + k 1 ( t ) x
where k 0 ( t ) and k 1 ( t ) are unknown functions to be determined; k 0 ( t ) and k 1 ( t ) satisfy the system of differential equations as follows:
d 2 k 1 ( t ) d t 2 = β k 1 ( t )
d 2 k 0 ( t ) d t 2 = ε d k 1 ( t ) d t
Case 1: when β > 0 :
k 1 ( t ) = b 1 e β t + b 2 e β t
k 0 ( t ) = a 1 + a 2 t + ε b 1 β e β t ε b 2 β e β t
Thus Equation (62) has the exact solution that follows:
f ( x , t ) = ( a 1 + a 2 t + ε b 1 β e β t ε b 2 β e β t ) + ( b 1 e β t + b 2 e β t ) x
where a 1 , a 2 , b 1 and b 2 are arbitrary constants.
Case 2: when β = 0 :
k 1 ( t ) = b 1 + b 2 t
k 0 ( t ) = ε b 2 2 t 2 + a 1 + a 2 t
Thus Equation (62) has the exact solution that follows:
f ( x , t ) = ( ε b 2 2 t 2 + a 1 + a 2 t ) + ( b 1 + b 2 t ) x
where a 1 , a 2 , b 1 and b 2 are arbitrary constants.
Case 3: when β < 0 :
k 1 ( t ) = b 1 cos β t + b 2 sin β t
k 0 ( t ) = a 1 + a 2 t ε b 1 β sin β t + ε b 2 β cos β t
Thus Equation (62) has the exact solution that follows:
f ( x , t ) = ( a 1 + a 2 t ε b 1 β sin β t + ε b 2 β cos β t ) + ( b 1 cos β t + b 2 sin β t ) x
where a 1 , a 2 , b 1 and b 2 are arbitrary constants.
The subspace I n + 1 = L { 1 , x , x 2 , , x n } , n N is invariant under N [ f ] = α 2 f x 2 + β x f x and x f , as
N [ k 0 + k 1 x + + k n x n ] = α k 2 + ( α k 3 + β k 1 ) x + + ( β k n 1 ( n 1 ) ) x n 1 + β k n n x n I n + 1
( k 0 + k 1 x + + k n x n ) x = k 1 + k 2 x + + k n 1 x n 1 I n + 1
Thus we obtain infinitely many invariant subspaces for Equation (62), which in turn yield infinitely many solutions.

5. Conclusions

The present article develops the invariant subspace method for solving certain fractional space and time derivative nonlinear partial differential equations with fractional-order mixed partial derivatives. Using the invariant subspace method, FPDEs are reduced to systems of FDEs; then they are solved by known analytic methods. In general, FPDEs admit more than one invariant subspace, each of which that has the exact solution. In fact, FPDEs admit infinitely many invariant subspaces. The invariant subspace method is used to derive closed-form solutions of fractional space and time derivative nonlinear partial differential equations with fractional-order mixed partial derivatives along with certain kinds of initial conditions. Thus, the invariant subspace method represents an effective and powerful tool for exact solutions of a wide class of linear/nonlinear FPDEs.
The bases of invariant subspaces usually are orthogonal polynomials, Mittag–Leffler functions, trigonometric functions, and so on. What kinds of spaces are the invariant subspaces of one FPDE? At present we can only try one by one. Although we have found some invariant subspaces of the equations examples above, are there any more invariant subspaces of the equations? We hope to find a simple discriminant method for finding the correct invariant subspaces for FPDEs.

Acknowledgments

This research is partially supported by the Doctoral Fund of Education Ministry of China (20134219120003) and the National Natural Science Foundation of China (61473338).

Author Contributions

All authors made equal contribution to this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Oldham, K.B.; Spanier, J. The Fractional Calculus; Academic Press: New York, NY, USA, 1974. [Google Scholar]
  2. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; Wiley: New York, NY, USA, 1993. [Google Scholar]
  3. Diethelm, K. The Analysis of Fractional Differential Equations; Springer: New York, NY, USA, 2010. [Google Scholar]
  4. Kilbas, A.A.; Trujillo, J.J.; Srivastava, H.M. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  5. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  6. Ansari, A.; Askari, H. On fractional calculus of A2n+1(x) function. Appl. Math. Comput. 2014, 232, 487–497. [Google Scholar]
  7. Podlubny, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1998. [Google Scholar]
  8. Daftardar-Gejji, V.; Jafari, H. Adomian decomposition: A tool for solving a system of fractional differential equations. J. Math. Anal. Appl. 2005, 301, 508–518. [Google Scholar] [CrossRef]
  9. Eslami, M.; Vajargah, B.F.; Mirzazadeh, M.; Biswas, A. Applications of first integral method to fractional partial differential equations. Indian J. Phys. 2014, 88, 177–184. [Google Scholar] [CrossRef]
  10. Sahadevan, R.; Bakkyaraj, T. Invariant analysis of time fractional generalized Burgers and Korteweg-de Vries equations. J. Math. Anal. Appl. 2012, 2, 341–347. [Google Scholar] [CrossRef]
  11. Ouhadan, A.; El Kinani, E.H. Lie symmetry analysis of some time fractional partial differential equations. Int. J. Mod. Phys. Conf. Ser. 2015, 38, 1560075. [Google Scholar]
  12. Bakkyaraj, T.; Sahadevan, R. Approximate analytical solution of two coupled time fractional nonlinear Schrodinger equations. Int. J. Appl. Comput. Math. 2016, 2, 113–135. [Google Scholar] [CrossRef]
  13. Zhukovsky, K.V.; Srivastava, H.M. Operational solution of non-integer ordinary and evolution-type partial differential equations. Axioms 2016, 5, 29. [Google Scholar] [CrossRef]
  14. Zhukovsky, K.V. The Operational Solution of Fractional-Order Differential Equations as well as Black-Scholes and Heat-Conduction Equations. Mosc. Univ. Phys. Bull. 2016, 71, 237–244. [Google Scholar] [CrossRef]
  15. Zhukovsky, K.V. Operational method of solution of linear non-integer ordinary and partial differential equations. SpringerPlus 2016, 5, 119. [Google Scholar] [CrossRef] [PubMed]
  16. Pandir, Y.; Huseyin Duzgun, H. New Exact Solutions of Time Fractional Gardner Equation by Using New Version of F-Expansion Method. Commun. Theor. Phys. 2017, 67, 9–14. [Google Scholar] [CrossRef]
  17. Moslehi, L.; Ansari, A. On M-Wright transforms and time-fractional diffusion equations. Integral Transforms Spec. Funct. 2017, 28, 113–129. [Google Scholar] [CrossRef]
  18. Aghili, A.; Aghili, J. Exponential differential operators for singular integral equations and space fractional Fokker-Planck equation. Bol. Soc. Parana. Mat. 2018, 36, 223–233. [Google Scholar] [CrossRef]
  19. Aghili, A. Fractional Black-Scholes equation. Int. J. Financ. Eng. 2017, 4, 1750004. [Google Scholar] [CrossRef]
  20. Galaktionov, V.A. Invariant subspaces and new explicit solutions to evolution equations with quadratic nonlinearities. Proc. R. Soc. Edinb. Sect. A Math. 1995, 125, 225–246. [Google Scholar] [CrossRef]
  21. Svirshchevskii, S.R. Invariant linear spaces and exact solutions of nonlinear evolution equations. J. Nonlinear Math. Phys. 1996, 3, 164–169. [Google Scholar] [CrossRef]
  22. Galaktionov, V.A.; Svirshchevskii, S.R. Exact Solutions and Invariant Subspaces of Nonlinear Partial Differential Equations in Mechanics and Physics; Chapman and Hall/CRC: London, UK, 2007. [Google Scholar]
  23. Gazizov, R.K.; Kasatkin, A.A. Construction of exact solutions for fractional order differential equations by invariant subspace method. Comput. Math. Appl. 2013, 66, 576–584. [Google Scholar] [CrossRef]
  24. Harris, P.A.; Garra, R. Analytic solution of nonlinear fractional Burgers-type equation by invariant subspace method. Nonlinear Stud. 2013, 20, 471–481. [Google Scholar]
  25. Harris, P.A.; Garra, R. Nonlinear time-fractional dispersive equations. arXiv, 2014; arXiv:1410.8085v1. [Google Scholar]
  26. Sahadevan, R.; Bakkyaraj, T. Invariant subspace method and exact solutions of certain nonlinear time fractional partial differential equations. Fract. Calc. Appl. Anal. 2015, 18, 146–162. [Google Scholar] [CrossRef]
  27. Ouhadan, A.; El Kinani, E.H. Invariant subspace method and fractional modified Kuramoto-Sivashinsky equation. arXiv, 2015; arXiv:1503.08789v1. [Google Scholar]
  28. Sahadevan, R.; Prakash, P. Exact solution of certain time fractional nonlinear partial differential equations. Nonlinear Dyn. 2016, 85, 659–673. [Google Scholar] [CrossRef]
  29. Choudhary, S.; Daftardar-Gejji, V. Invariant Subspace Method: A tool for solving fractional partial differential equations. arXiv, 2016; arXiv:1609.04209v1. [Google Scholar]
  30. Zhukovsky, K.V. Operational solution for some types of second order differential equations and for relevant physical problems. J. Math. Anal. Appl. 2016, 446, 628–647. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Jiang, J.; Feng, Y.; Li, S. Exact Solutions to the Fractional Differential Equations with Mixed Partial Derivatives. Axioms 2018, 7, 10. https://doi.org/10.3390/axioms7010010

AMA Style

Jiang J, Feng Y, Li S. Exact Solutions to the Fractional Differential Equations with Mixed Partial Derivatives. Axioms. 2018; 7(1):10. https://doi.org/10.3390/axioms7010010

Chicago/Turabian Style

Jiang, Jun, Yuqiang Feng, and Shougui Li. 2018. "Exact Solutions to the Fractional Differential Equations with Mixed Partial Derivatives" Axioms 7, no. 1: 10. https://doi.org/10.3390/axioms7010010

APA Style

Jiang, J., Feng, Y., & Li, S. (2018). Exact Solutions to the Fractional Differential Equations with Mixed Partial Derivatives. Axioms, 7(1), 10. https://doi.org/10.3390/axioms7010010

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