Advances in Mathematical Physics

Volume 2018, Article ID 4596506, 8 pages

https://doi.org/10.1155/2018/4596506

## Exact Solutions with Variable Coefficient Function Forms for Conformable Fractional Partial Differential Equations by an Auxiliary Equation Method

^{1}School of Mathematical Sciences, Qufu Normal University, Qufu, Shandong 273165, China^{2}School of Mathematics and Statistics, Shandong University of Technology, Zibo, Shandong 255049, China

Correspondence should be addressed to Qinghua Feng; moc.anis@auhqf

Received 22 March 2018; Accepted 29 July 2018; Published 5 August 2018

Academic Editor: Ziemowit Popowicz

Copyright © 2018 Fanwei Meng and Qinghua Feng. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

In this paper, an auxiliary equation method is introduced for seeking exact solutions expressed in variable coefficient function forms for fractional partial differential equations, where the concerned fractional derivative is defined by the conformable fractional derivative. By the use of certain fractional transformation, the fractional derivative in the equations can be converted into integer order case with respect to a new variable. As for applications, we apply this method to the time fractional two-dimensional Boussinesq equation and the space-time fractional (2+1)-dimensional breaking soliton equation. As a result, some exact solutions including variable coefficient function solutions as well as solitary wave solutions for the two equations are found.

#### 1. Introduction

Fractional differential equations are the generalizations of classical differential equations with integer order derivatives. It is well known that fractional partial differential equations have proved to be very useful in many research fields such as physics, mathematical biology, engineering, fluid mechanics, plasma physics, optical fibers, neural physics, solid state physics, viscoelasticity, electromagnetism, electrochemistry, signal processing, control theory, chaos, finance, and fractal dynamics. In the last few decades, research on various aspects for fractional partial differential equations has received more and more attention by many authors, such as the oscillation [1, 2] and existence and uniqueness [3–5]. In order to better understand the physical process described by fractional partial differential equations, one usually needs to obtain exact solutions or numerical solutions for fractional partial differential equations. And so far a lot of effective methods have been developed and used by many authors. For example, these methods include the finite difference method [6], the method [7–10], the Jacobi elliptic function method [11], the projective Riccati equation method [12], the modified Kudryashov method [13–19], the exp method [20–25], the ansatz method [26], the first integral method [27–29], and the subequation method [30–34]. Based on these methods, a lot of fractional partial differential equations have been investigated.

In this paper, we develop an auxiliary equation method for solving fractional partial differential equations, where the fractional derivative is defined in the sense of the conformable fractional derivative. Our aim is to seek exact solutions with variable coefficient function forms for some certain fractional partial differential equations.

The conformable fractional derivative is defined by [35]The following properties for the conformable fractional derivative are well known, which can be easily proved due to the definition of the conformable fractional derivative. , where is a constant.

The paper is organized as follows. In Section 2, we propose the description of the auxiliary equation method for seeking exact solutions with variable coefficient function forms for fractional partial differential equations. Then in Section 3, we apply the method to solve the time fractional two-dimensional Boussinesq equation and the space-time fractional (2+1)-dimensional breaking soliton equation. In Section 4, we present some concluding statements.

#### 2. Description of the Auxiliary Equation Method

In this section, we give the description of the auxiliary equation method for solving fractional partial differential equations.

Suppose that a fractional partial differential equation in the independent variables is given by where is an unknown function, the orders of the fractional derivatives are, for example, , and is a polynomial in and its various partial derivatives including fractional derivatives. Without loss of generality, next we may assume that the fractional partial derivatives are related to the variables , while the other variables are related to integer order derivatives.

*Step 1. *For those variables involving fractional derivatives, fulfil corresponding fractional transformations so that the fractional partial derivatives can be converted into integer order partial derivatives with respect to new variables.

Taking the expressions and , for example, one can use two fractional transformations and and denote . Then due to the properties and of the conformable fractional derivative, one can obtain that , Therefore, the original fractional partial differential equation can be converted into another partial differential equation of integer order as follows:

*Step 2. *Suppose that the solution of (3) can be expressed by a polynomial in as follows: where , are all unknown functions to be determined later with , and satisfies some certain auxiliary equation with the following form: whose solutions are known. Furthermore, as , then the degree of is usually one more than . The positive integer can be determined by considering the homogeneous balance of the degrees between the highest order derivatives and nonlinear terms appearing in (3).

*Step 3. *Substituting (4) into (3), using the relation between and derived from (5), and collecting all terms with the same order of together, the left-hand side of (3) is converted to another polynomial in . Equating each coefficient of this polynomial to zero yields a set of partial differential equations for , .

*Step 4. *Solving the equations yielded in Step 3 and using the solutions of (5), together with the fractional transformations introduced in Step 1, one can obtain exact solutions for (2).

*Remark 1. *In the present method, the expression of the transformation denoted by is underdetermined, and the coefficients in (4) are variable coefficient functions, which may contribute to the seeking of exact solutions with variable coefficient function forms. Furthermore, if (5) are selected for some different forms, such as the Riccati equation, Bernoulli equation, and Jacobi elliptic equation, then different exact solutions for (2) can be obtained correspondingly.

*Remark 2. *As the partial differential equations yielded in Step 3 are usually overdetermined, we may choose some special forms of as in the following.

#### 3. Applications of the Auxiliary Equation Method

##### 3.1. Time Fractional Two-Dimensional Boussinesq Equation

We consider the time fractional two-dimensional Boussinesq equation with the following form:

The fractional Boussinesq equation is used in the analysis of long waves in shallow water and also used in the analysis of many other physical applications, such as the percolation of water in a porous subsurface of a horizontal layer of material. Tasbozan et al. [36] solved (7) using the Jacobi elliptic function expansion method and obtained a series of exact solutions with Jacobi elliptic function forms.

Now we use the proposed auxiliary equation method to solve (6). First we let and . Then , and (6) can be converted into the following form: Suppose that the solution of (7) can be expressed by a polynomial in as follows: where , are underdetermined functions and satisfies (5). Balancing the degrees of and in (7), one can obtain , which means . Thus, one has

Next we will discuss two cases, in which satisfies two certain auxiliary equations.

*Case 1. * satisfies the following Riccati equation: where .

Substituting (9) into (7), using (10), collecting all the terms with the same power of together, and equating each coefficient to zero yield a set of underdetermined partial differential equations , and . Solving these equations yields the following families of results, where are arbitrary constants and , are arbitrary functions.

*Family 1*

*Family 2*

*Family 3*

*Family 4*

*Family 5*

*Remark 3. *It is obvious that if we let , in Family 1, then the transformation denoted by becomes , which has been used by many authors so far in the literature.

On the other hand, it is well known that the solution of (10) is denoted by where is an arbitrary constant, so In particular, when , one can obtain

By a combination of the results denoted in Families 1–5 and (17), together with the expression of , one can obtain a series of exact solutions for the time fractional two-dimensional Boussinesq equation as follows: where , and are arbitrary functions. where , and are arbitrary functions. where , is an arbitrary constant, and is an arbitrary function. where , and are arbitrary constants. where , are arbitrary functions.

Similarly, by a combination of the results in Families 1–5 and (18), one can obtain the following solitary wave solutions, in which ), , are arbitrary constants, and , are arbitrary functions. where where where where where

In Figure 1, the solitary wave solution in (24) with some special parameters is demonstrated.