Assessment of the Exact Solutions of the Space and Time Fractional Benjamin-Bona-Mahony Equation via the G / G-Expansion Method , Modified Simple Equation Method , and Liu ’ s Theorem

Exact travelling wave solutions to the space and time fractional Benjamin-Bona-Mahony (BBM) equation defined in the sense of Jumarie’smodified Riemann-Liouville derivative via the (G/G) expansion and themodified simple equationmethods are presented in this paper. A fractional complex transformation was applied to turn the fractional BBM equation into an equivalent integer order ordinary differential equation. New complex type travelling wave solutions to the space and time fractional BBM equation were obtained with Liu’s theorem.Themodified simple equationmethod is not effective for constructing solutions to the fractional BBM equation.


Introduction
Nonlinear partial differential equations arise in a large number of physics, mathematics, and engineering problems.In the Soliton theory, the study of exact solutions to these nonlinear equations plays a very germane role, as they provide much information about the physical models they describe.In recent times, it has been found that many physical, chemical, and biological processes are governed by nonlinear partial differential equations of noninteger or fractional order [1][2][3][4].
In this paper, we apply the (  /) expansion method and the modified simple equation method to construct travelling wave solutions to the space and time fractional Benjamin-Bona-Mahony (BBM) equation in the sense of Jumarie's modified Riemann-Liouville derivative via a fractional complex transformation and we further get new complex type solutions to the equation by applying Liu's theorem [23].The Benjamin-Bona-Mahony equation is of the following form: (, )  +    (, )  +  (, )    (, )  −  2    (, )  2  = 0. ( This equation was introduced in [24] as an improvement of the Korteweg-de Vries equation (KdV equation) for modelling long waves of small amplitude in 1 + 1 dimensions.It is used in modelling surface waves of long wavelength in liquids, acoustic gravity waves in compressible fluids, and acoustic waves in anharmonic crystals.

ISRN Mathematical Physics
Jumarie's modified Riemann-Liouville derivative of order  with respect to  is defined as [25]      () = Some useful properties of the modified Riemann-Liouville derivative are listed below [25,26]: (3)

Description of the (𝐺 󸀠 /𝐺) Expansion Method
Here, we provide a brief explanation of the (  /) expansion method for finding travelling wave solutions of nonlinear fractional partial differential equations.Suppose that the nonlinear evolution equation is in the following form: where ℘ is a polynomial of (, ) and its derivatives (integer and fractional) with respect to  and . and  are parameters that describe the order of the space and time derivatives, respectively.

Theorem 1. Fractional Complex Transformation. To transform (4) into a nonlinear ordinary differential equation (ODE)
of integer order by applying a fractional complex transformation proposed by Li and He [27], where  is an arbitrary constant and (4) reduces to a nonlinear integer order ODE of the following form: The method assumes that the solution to (6) can be expressed as a polynomial in (  /): where  0 ,  1 , . . .,   are constants to be determined and  = () is a solution of the linear ordinary differential equation of the following form: where  and  are arbitrary constants.The general solution of (8) gives [28] the following: Equation ( 6) is integrated as long as all the terms contain derivatives, where integration constants are considered to be zero.
To determine , one considers the homogenous balance between the highest order derivative and the highest order nonlinear term(s).Substitute (7) with the determined value of  into (6), and collect all terms with the same order of (  /) together.If the coefficients of (  /) vanish separately, one has a set of algebraic equations in  0 ,  1 , . . .,   , ,  and  that is solved with the aid of Mathematica.
Theorem 2 (Liu's theorem [23]).If a nonlinear evolution equation has a kink-type solution in the form of where   is a polynomial of degree , then it has a certain kinkbell-type solution in the following form: where  is the imaginary number unit.

Description of the Modified Simple Equation Method
Here, we highlight the basic ideas of the modified simple equation method.Suppose that the fractional nonlinear evolution equation is in the following form: where ℘ is a polynomial of  and its fractional partial derivatives with respect to  and .Suppose a fractional complex transformation (see (5)) which allows us to reduce (12) to an ordinary differential equation (ODE) of the following form: (,   ,   , . ..) = 0.
We suppose that the solution to (13) can be expressed as a polynomial of (  /) in the following form: where  is a positive integer obtained by balancing the highest order derivatives and the highest order nonlinear terms in (13).  ( = 0, 1, 2, . . ., ) are arbitrary constants to be determined such that   ̸ = 0 and () is an unknown function to be determined and is not a solution of any predefined differential equation.
We substitute (14) and its derivatives into (13) taking into account the function ().As a result of this substitution we obtain a polynomial of  − () with the derivatives of ().Equating all the coefficients of  − () to zero yields a system of equations which can be solved to obtain   and ().Finally, substituting the values of   , () and its derivative   () into (14) gives the exact solution of (12).

Solutions for Space and Time
Fractional BBM Equation via (  /) Expansion Method.In this section, we apply the (  /) expansion method to construct travelling wave solution of the space and time fractional Benjamin-Bona-Mahony equation.Consider the space and time fractional BBM equation (1) which can be written in subscript notation as We make the transformation (, ) = (),  =   /Γ(1 + ) −   /Γ(1 + ).Equation ( 15) becomes Integrating ( 16) once with respect to  and setting the integration constant to zero yields To get , we consider the homogenous balance between the highest order derivative   and the highest order nonlinear term  2 : Then (7) becomes where  0 ,  1 , and  2 are constants to be determined later.From (19), Substituting (19) and its derivatives into (17) and collecting all terms with the same power of (  /) together yields a simultaneous set of nonlinear algebraic equations as follows: Solving this algebraic system of equations with the aid of Mathematica yields the following solution: Substituting the solution to the algebraic equation and the general solution to ( 8) into (19), we obtain three types of travelling wave solutions of the space and time fractional Benjamin-Bona-Mahony equation.

Solutions for Space and Time
Fractional BBM Equation via MSE Method.Now, we apply the modified simple equation method to construct travelling wave solutions of the space and time fractional Benjamin-Bona-Mahony equation.The procedure for transforming the nonlinear fractional PDE into an ODE in the MSE method is the same as that of the (  /) expansion method.Hence, from (17), we have We get  = 2 by considering the homogenous balance between the highest order derivative and the highest order nonlinear term.Then (14) becomes where  0 ,  1 , and  2 are constants to be determined such that  2 ̸ = 0 and () is an unidentified function to be determined.It is easy to find that Substituting these values into (17) and equating the coefficients of  0 ,  −1 ,  −2 ,  −3 , and  −4 , respectively, to zero, we obtain Solving ( 37) and ( 41), respectively, yields From ( 38), (39), and (40), we have Case 1.When  1 = 12  /  , we get  1 = ±12√( − 1).This yields an absurd solution and hence this case is discarded.
Case 2. When  1 = 0 and   = 0, we have where  1 is a constant of integration.Integrating (44) with respect to , we obtain where  2 is a constant of integration.Substituting the value of   and  from (44) and (45), respectively, into (35) yields when  0 = 0, we obtain the exact solution of (15) as when  0 = 2( − 1), we obtain the exact solution of (15) as Since  1 and  2 are arbitrary constants, therefore, if we set  1 = 1 and  2 = 0, we have the solutions ( 47) and (48) as Bulent and Erdal [29] used direct algebraic method to get new complex travelling wave solutions to the BBM equation ( 15) by employing a transformation given by (, ) = (),  = ( − ).The results reported by [29] are A generalized (  /) expansion method which uses the Klein-Gordon equation as the auxiliary equation was applied by Yanhong and Baodan [30] to a form of BBM equation given by The travelling wave solutions to the BBM equation obtained by [30] using a transformation given by (, ) = (),  =  −  are where (), (), (), (), (), and () are Jacobi elliptic functions and 0 <  < 1 gives the modulus of the Jacobi elliptic functions.These functions generate hyperbolic functions when  → 1 and trigonometric functions when  → 0 as follows: Remark 3. When the arbitrary constants  1 and  2 are taken to be zero separately in (23), (28), and (31), we get ( 24), ( 25), ( 29), (30), and (32) which are the travelling wave solutions to the space and time fractional BBM equation obtained via the (  /) expansion method.Equations (49) are the travelling wave solution to the space and time fractional BBM equation obtained through the modified simple equation method.We see that the modified simple equation is not very effective for constructing travelling wave solutions to the space and time fractional BBM equation because the integer  obtained from considering the homogenous balance between the highest nonlinear terms and the highest derivative is greater than 1.
Remark 5.All the travelling wave solutions to the space and time fractional Benjamin-Bona-Mahony equation obtained were checked by putting them back into (1) with the aid of Mathematica.
Remark 7. We also examine the solutions to the BBM equation from [30] for  = 1.For  =  = 1 and  = 1,  The shapes of ( 24), ( 25), (29), and (30) for selected values of , , and  at different  and  values are presented in Figures 1-4.The solution  1 (see (24)) is a singular kink shaped travelling wave solution of (15).Figure 1 shows the shape of (24) with  = 1,  = 3, and  = 2 for the interval 0 ≤ ,  ≤ 10 at different  and  values.From Figure 1 and (, ) = ( − ), we observe that, for  > 0, the propagation of the wave will be in the positive -direction.If we take  < 0, then the propagation of the wave will be in the negative -direction.The solution  2 (see (25)) is a bell-shaped Soliton solution of (15).Figure 2 shows the shape of (25) with  = 1,  = 3, and  = 2 for the interval 0 ≤ ,  ≤ 10 at different  and  values.From Figure 2 and (, ) = ( − ), we also observe that, for  > 0, the propagation of the wave will be in the positive -direction.If we take  < 0, then the propagation of the wave will be in the negative -direction.Figures 3 and 4 show, respectively, the shape of  3 (see (29)) and  4 (see (30)) with  = 1 and  =  = 2 for the interval 0 ≤ ,  ≤ 10 at different  and  values. 3 (see (29)) and  4 (see (30)) are exact periodic travelling wave solutions of the space and time fractional Benjamin-Bona-Mahony equation.

Conclusion
Hyperbolic, trigonometric, and rational function travelling wave solutions to the space and time fractional Benjamin-Bona-Mahony equation defined in the sense of Jumarie's modified Riemann-Liouville derivative have been obtained using the (  /) expansion method and modified simple equation method.Liu's theorem was also applied to obtain new complex travelling wave solutions to the space and time fractional BBM equation.The results were compared to those obtained by Bulent and Erdal [29] using the direct algebraic method for complex travelling wave solution and those by Yanhong and Baodan [30] using the generalized (  /) expansion method with the Klein-Gordon equation as an auxiliary equation.The (  /) expansion method is equivalent to the direct algebraic method and the generalized (  /) expansion method with some of their reported solutions found not to satisfy the BBM equation.Also, the modified simple equation method is not effective for constructing travelling wave solutions to the space and time fractional BBM equation because the index  obtained from considering the homogenous balance between the highest nonlinear term and the highest derivative is greater than 1.The improvement of the modified simple equation scheme for evolution equations with  > 1 is still an open problem.