Abstract

The fractional derivatives in the sense of the modified Riemann-Liouville derivative and Feng’s first integral method are employed to obtain the exact solutions of the nonlinear space-time fractional ZKBBM equation and the nonlinear space-time fractional generalized Fisher equation. The power of this manageable method is presented by applying it to the above equations. Our approach provides first integrals in polynomial form with high accuracy. Exact analytical solutions are obtained through establishing first integrals. The present method is efficient and reliable, and it can be used as an alternative to establish new solutions of different types of fractional differential equations applied in mathematical physics.

1. Introduction

Fractional differential equations have been the focus of many studies due to their frequent appearance in various applications in physics, biology, engineering, signal processing, systems identification, control theory, finance, and fractional dynamics [13]. Recently, a large amount of literature has been provided to construct the solutions of fractional ordinary differential equations, fractional partial differential equations (PDEs), and integral equations of physical interest. Several powerful methods have been proposed to obtain approximate and exact solutions of fractional differential equations, such as the Adomian decomposition method [4, 5], the variational iteration method [68], the homotopy analysis method [912], the homotopy perturbation method [1315], the Lagrange characteristic method [16], and the fractional subequation method [17].

In [18], Jumarie proposed a modified Riemann-Liouville derivative. With this kind of fractional derivative and some useful formulas, we can convert fractional differential equations into integer-order differential equations by variable transformation.

Feng [19] has introduced a reliable and effective method called Feng’s first integral method to look for travelling wave solutions of nonlinear PDEs. The basic idea of this method is to construct a first integral with polynomial coefficients of an explicit form to an equivalent autonomous planar system by using division theorem. This method in comparison with other methods has many advantages; it avoids a great deal of complicated and tedious calculation and provides exact and explicit travelling solutions with high accuracy.

Feng’s first integral method [2026] can be used to construct the exact solutions for some time fractional differential equations.

Among the nonlinear PDEs there are two important examples of fundamental interest in mathematical-physical models: the nonlinear ZKBBM equation that describes the wave water phenomena in the large wavelength limit and one special case of nonlinear reaction-diffusion equation better known as the generalized Fisher equation. Let us review some general aspects about these two types of nonlinear PDEs.

The well-known Zakharov-Kuznetsov-Benjamin-Bona-Mahony equation (ZKBBM equation) also referred to as the regularized long wave equation (see [27, 28]) was originally derived for water waves under the assumption of small wave amplitude and large wavelength. The ZKBBM equation describes approximately the unidirectional propagation of long waves in certain nonlinear dispersive systems and has been introduced in order to give a satisfactory description for the small wave amplitude in the large wavelength regime in many other physical systems. It has been used to account adequately for observable phenomena such as the interaction of solitary waves and dissipationless undular shocks; see [27].

The details of the derivation of the ZKBBM equation differ, in different examples, but in general the essentials are as follows. The first property of the systems in question is that dispersive effects on infinitesimal waves vanish in the limit as wavelength becomes infinite, and the limiting phase speed is a constant (). Thus waves of extreme length, along the propagation direction, are described by whose general solution is an arbitrary differentiable function of , where and . In certain applications (2) already has some validity as an approximation for real waves of sufficiently small amplitude and large wavelength but is not a valid approximation for very large times, during which nonlinear and frequency dispersive effects are accumulated to a significant level. These effects require the introduction of at least a first order approximation model to extend the range of validity of the model (1).

For the linearized equation (1) all values of are propagated with the same velocity . When the nonlinear effects are taken into account, the characteristic velocity becomes dependent linearly on for the waves of finite but small amplitude; therefore in the first order approximation we obtain where is a constant. After the substitution of (3) into (2) the following equation is obtained: Additionally the dispersive effects as suffered by waves of finite but large wavelength in the model (2) are ignored. In the first order approximation the dispersive effects transform the model (2) into the following equation: Although for a self-consistent theory the two kinds of small effects have to be considered simultaneously, the outcome of such theories can generally be anticipated by considering the two effects separately. That is, in a first approximation, the above corrections are added separately into (2), and higher order terms appear negligible in this approximation.

When the physical parameters and the scaling factors presented in a particular example are appropriately absorbed into the definition of the dependent variable and the independent variables and , which are, respectively, proportional to distance in the physical system and to time, the ZKBBM equation is given by: Equation (6) can be extended in order to consider fractional order derivatives obtaining the ZKBBM space-time fractional equation [29]: where and denotes the fractional derivative of order with respect to the independent variables and .

The second nonlinear PDE mentioned before, the generalized Fisher equation [3032], can be obtained as a special case of the reaction-diffusion equation. The reaction-diffusion equation plays a fundamental role in a great number of various models of heat and reaction-diffusion processes, mathematical biology, chemistry, genetics, and many other fields.

We start with the reaction-diffusion equation with power nonlinearity (see, for example, [33]): here is a sufficiently smooth polynomial function satisfying the relations .

One particular expression for the function in (8) which satisfies these conditions is [33] where are constants that satisfy the condition imposed over and is an arbitrary constant. Equation (9) defines a special class of the nonlinear heat equation, which includes many important applications in some of the areas mentioned before. If we consider the special case in (9) where , , and , we obtain the generalized Fisher equation: where in the original mathematical biology model denotes the population density and represents the population supply due to births and deaths.

The above equation has been introduced as a model for the propagation of a mutant gene in mathematical biology [3032] and can be extended to fractional order derivatives obtaining the generalized Fisher space-time fractional equation: The present work investigates the applicability and effectiveness of Feng’s first integral method to obtain new exact analytical solutions for the nonlinear space-time fractional Zakharov-Kuznetsov-Benjamin-Bona-Mahony (ZKBBM) partial differential equation and the nonlinear space-time fractional generalized Fisher equation.

We will show that, for the analytical solutions of the ZKBBM fractional partial differential equation, one of the advantages of Feng’s first integral method is the evaluation of the constants involved in the analytical solutions that have not been evaluated within the subequation method [29]. Furthermore for the analytical solutions of the space-time fractional generalized Fisher equation, by applying Feng’s first integral method, we will obtain three new analytical solutions that have not been obtained in previous works [30, 31, 33].

2. The Modified Riemann-Liouville Derivative and Feng’s First Integral Method

In this section we present the main ideas of Feng’s first integral method. This method considers the Jumarie modified Riemann-Liouville fractional derivative of order ; we first give some definitions and properties of the modified Riemann-Liouville derivative which are used further in this paper.

Assume that , denotes a continuous (but not necessarily differentiable) function. The Jumarie modified Riemann-Liouville derivative of order is defined by the expression [18]

Some properties of the fractional modified Riemann-Liouville derivative are

Now in order to introduce Feng’s first integral method [19], let us consider the space-time fractional differential equation with independent variables and some function : Using the variable transformation where and are constants to be determined later, the fractional differential equation (14) is reduced to a nonlinear ordinary differential equation where .

We assume that (16) has a solution in the form and we introduce a new independent variable , which leads to the following new system of equations: Now, let us introduce the central idea of Feng’s first integral method. By using the division theorem for two variables in the complex domain which is based on the Hilbert-Nullstellensatz theorem [34], we can obtain one first integral to (18) which can reduce (16) to a first order integrable ordinary differential equation. An exact solution to (14) is then obtained by solving this equation directly.

Division Theorem. Suppose that and are polynomials in , and is irreducible in . If vanishes at all zero points of , then there exists a polynomial in such that

3. Applications

In this section, we illustrate the applicability of Feng’s first integral method to solve the nonlinear fractional partial differential equations (7) and (11).

Example 1. As a first example we consider the space-time fractional ZKBBM equation (7): where and are arbitrary constants. By considering the travelling wave transformation where and are constants, substituting (21) into (20), we can reduce (20) into an ordinary differential equation (ODE): or Integrating (23) and being the integration constant, we have and using (17) and (18), (24) is equivalent to the two-dimensional autonomous system According to Feng’s first integral method, we suppose that and are nontrivial solutions of (25), and is an irreducible polynomial in the complex domain such that where are polynomials of and . Due to the division theorem, there exists a polynomial in the complex domain such that Suppose that ; by equating the coefficients of on both sides of (27), we havewhere since are polynomials, then from (28a) we deduce that is constant and . For simplicity, we take . Balancing the degrees of and , we conclude that where is an arbitrary constant. Substituting , , , and into (28c) and setting all the coefficients of powers of to be zero, we obtain a system of nonlinear algebraic equations and by solving it, we obtain where From the above conditions, four solutions can be obtained for the space-time fractional ZKBBM partial differential equation.
Case  1. When , , , and , Now from the condition , we obtain and therefore taking into account that and from (33) it follows that Combining this first integral with the two-dimensional autonomous system of (18), the second order differential equation (24) can be reduced to the following first order differential equation: We consider the following transformation: introduced by Wang [35]. With this transformation (36) can be written as or where and .
By solving (39), the following solutions are obtained: or if we consider the case and , and taking into account (29) and (37), the solutions (40) and (41) reduce to or with , where and .
Case  2. When , , , and , Taking into account that and from (44) it follows that Combining this first integral with the two-dimensional autonomous system given by (18), the second order differential equation (24) can be reduced to the following first order differential equation: By using the transformation defined in (37), (46) can be written as where .
By solving (47), the following solutions are obtained: or if we consider the case and , and taking into account (29) and (37), the solutions (48) and (49) reduce to or with , where and .
Case  3. When , , , and , Taking into account that and from (52) it follows that Combining this first integral with the two-dimensional autonomous system given by (18), the second order differential equation (24) can be reduced to the following first order differential equation: By using the transformation defined in (37), (54) can be written as where and .
By solving (55), the following solutions are obtained: or if we consider the case and , and taking into account (29) and (37), the solutions (56) and (57) reduce to or with , where and .
Case  4. When , , , and , Taking into account that and from (60) it follows that Combining this first integral with the two-dimensional autonomous system given by (18), the second order differential equation (24) can be reduced to the following first order differential equation: By using the transformation defined in (37), (62) can be written as where .
By solving (63), the following solutions are obtained: or if we consider the case and , and taking into account (29) and (37), the solutions (64) and (65) reduce to or with , where and .
For the nonlinear space-time fractional ZKBBM equation, we have successfully recovered the previously known solutions, (42), (43), (48), (49), (58), (59), (66), and (67), that have been found in [29]; furthermore, Feng’s first integral method has evaluated the constant which has not been evaluated in this previous work.

Example 2. Next we consider the nonlinear fractional generalized Fisher equation (11): By considering the following transformation, introduced by Wang [35], (68) can be written as By considering the travelling wave transformation, where and are constants, (70) can be reduced to the following nonlinear ordinary differential equation: We assume that (72) has a solution in the form , and using (18), (72) is equivalent to the two-dimensional autonomous system: According to Feng’s first integral method, we suppose that and are nontrivial solutions of (73), and is an irreducible polynomial in the complex domain such that where are polynomials of and . Due to the division theorem, there exists a polynomial in the complex domain such that Suppose that ; by equating the coefficients of on both sides of (75), we haveSince are polynomials, then from (76a) we deduce that is constant and . For simplicity, we take . Balancing the degrees of and , we conclude that Substituting , , , and into (76c) and setting all the coefficients of powers of to be zero, we obtain a system of nonlinear algebraic equations and by solving it, we obtain therefore we obtain the following set of solutions: From the above conditions, four solutions can be obtained.
Case  1. When and , Now from the condition we obtain and therefore taking into account that and from (81) it follows that If we consider the definition given in (18) one finally obtains By solving the above equation the solution to (72) is obtained: if we consider the case , the solution (84) reduces to with . Taking into account (69), the solution to the Fisher equation (68) is given by with
Case  2. When and , Now from the condition we obtain and therefore taking into account that and from (88), it follows that If we consider the definition given in (18), one finally obtains By solving the above equation, the solution to (72) is obtained: if we consider the case , the solution (92) reduces to with . Taking into account (69), the solution to the Fisher equation (68) is given by with
Case  3. When and , Now from the condition , we obtain and therefore taking into account that and from (96), it follows that If we consider the definition given in (18), one finally obtains By solving the above equation, the solution to (72) is obtained: if we consider the case , the solution (100) reduces to with . Taking into account (69), the solution to the Fisher equation (68) is given by with
Case  4. When and , Now from the condition , we obtain and therefore taking into account that and from (104), it follows that If we consider the definition given in (18), one finally obtains By solving the above equation, the solution to (72) is obtained: if we consider the case , the solution (108) reduces to with . Taking into account (69), the solution to the Fisher equation (68) is given by with For the nonlinear space-time generalized fractional Fisher equation, we have successfully recovered the previously known solution (110) that has been found in [31], but to the best of our knowledge the solutions (86), (94), and (102) have not been obtained previously in the literature.

4. Conclusions

Feng’s first integral method was applied successfully to obtain new exact analytical solutions of the nonlinear space-time fractional ZKBBM equation and the nonlinear space-time fractional generalized Fisher equation. These solutions can be very useful as a starting point of comparison when some approximate methods are applied to these nonlinear space-time fractional equations. The performance of Feng’s first integral method is reliable and effective to obtain new solutions. This method has more advantages: it is direct and concise. Thus, the proposed method can be extended to solve many systems of nonlinear fractional partial differential equations in mathematical and physical sciences.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors gratefully acknowledge the Universidad Autónoma de la Ciudad de Mexico for supporting and facilitating this research work.