A note on improved F-expansion method combined with Riccati equation applied to nonlinear evolution equations

The purpose of this article is to present an analytical method, namely the improved F-expansion method combined with the Riccati equation, for finding exact solutions of nonlinear evolution equations. The present method is capable of calculating all branches of solutions simultaneously, even if multiple solutions are very close and thus difficult to distinguish with numerical techniques. To verify the computational efficiency, we consider the modified Benjamin–Bona–Mahony equation and the modified Korteweg-de Vries equation. Our results reveal that the method is a very effective and straightforward way of formulating the exact travelling wave solutions of nonlinear wave equations arising in mathematical physics and engineering.


Introduction
Nonlinear evolution equations (NLEEs) are widely used to describe many important phenomena and dynamical processes in mathematical physics and engineering. The investigation of exact solutions of NLEEs plays an important role in the study of these physical phenomena. In this paper, we highlight an analytical method, namely the improved F-expansion method, for finding exact solutions of NLEEs. Exact solutions provide a means to describe the salient features in various science, technology and engineering applications and can serve as a basis for perfecting and testing computer algebra software packages for solving NLEEs. It is significant that many equations of physics, chemistry and biology contain empirical parameters or empirical functions. Exact solutions allow researchers to design and run experiments by creating appropriate natural conditions in order to determine these parameters or functions.
In each of these aforementioned works, a variety of ansätz have been proposed for seeking travelling wave solutions of nonlinear differential equations. The choice of an appropriate ansätz is of great importance when using these analytical methods. Among those approaches, the proposed improved F-expansion method is a powerful tool to reveal more general solitons of NLEEs in mathematical physics and engineering. The main idea of this method is to express the travelling wave solutions of NLEEs in terms of functions that satisfy the Riccati equation F (ξ ) = k + F 2 (ξ ).
The major advantage of the improved F-expansion method over the existing other methods is that it provides more new exact travelling wave solutions. To demonstrate the efficiency and reliability of this proposed method, the mBBM equation and the mKdV equation have been solved in which new solutions are found. It is important to note that analysis of convergence and stability for the numerical methods is required, a distinct disadvantage when compared with analytical methods that do not require such an analysis. Apart from the physical relevance, the closed-form solutions of NLEEs can assist the numerical solvers to measure up to the accuracy of their results and thus aid in the convergence analysis.
The rest of the article has been prepared as follows. In §2, the improved F-expansion method is discussed in detail. In §3, we apply this method to obtain solutions to the NLEES mentioned above. In §4, we provide graphical representations of some of the obtained solutions. Section 5 contains the conclusion of our study.

Algorithm of improved F-expansion method
In this section, we describe the main steps of the improved F-expansion method for finding exact travelling wave solutions of NLEEs. To begin, consider the well-known Riccati equation: where F = F(ξ ) and the prime stands for derivatives with respect to ξ . We now present the three cases of the general solutions of the Riccati equation (3.1).

Case 1.
When k < 0, the general solutions are Case 2. When k > 0, the general solutions are Case 3. When k = 0, the general solution is where k is the real parameter. where u = u (x, y, t) is an unknown function, P is a polynomial in u (x, y, t) and its partial derivatives in which the highest order partial derivatives and the nonlinear terms are involved and the subscripts stand for the partial derivatives. The foremost steps of the method are given below.
Step 1. We introduce the travelling wave transformation, where λ is the speed of the travelling wave and substitute this into equation (3.2), yielding the ordinary differential equation (ODE): where Q is a polynomial of u and its derivatives and the superscripts specify the ordinary derivatives with respect to ξ .
Step 2. In many instances, equation (3.4) can be integrated term by term one or more times, yielding constants of integration, which can be set equal to zero for simplicity.
Step 3. We assume the travelling wave solution of equation (3.4) can be expressed by a polynomial in F(ξ ) as follows: N), λ and m are constants to be determined later.
Step 4. The positive integer N can be determined by using homogeneous balance between the highest order derivatives and the nonlinear terms appearing in the ODE (3.4).
then the degree of the other expressions will be given by Therefore, we can find the value of N in equation (3.5), using equation (3.6).
Step 5. Substituting equation (3.5) into equation (3.4) together with the value of N obtained in step 4, we obtain polynomials in F(ξ ). We set each coefficient of the resulting polynomial to zero, yielding an over-determined set of algebraic equations for α N , β N , m and λ.
Step 6. We suppose the values of the constants α N , β N , m and λ can be determined by solving the algebraic equations obtained in step 5. As the general solution of equation (  The Benjamin-Bona-Mahony (BBM) equation u t + u x + au n u x + bu xxt = 0 is a well-known NLEE that models long waves in a nonlinear dispersive system. The solution of the BBM equation exhibits solitonlike behaviour. The BBM equation is used in the analysis of the surface waves of long wavelength in liquids, hydromagnetic waves in cold plasma, acoustic-gravity waves in compressible fluids and acoustic waves in harmonic crystals. When n = 2, the BBM equation is called the modified BBM equation [3,38] and is given by

Applications
Now integrating equation (4.2) with respect to ξ once and setting the constant of integration to zero, we obtain bλu + a 3 Taking the homogeneous balance between the highest order nonlinear term u 3 and the derivative term u from equation (4.3), yields 3N = N + 2, i.e. N = 1.
Now substituting equation (4.4) into equation (4.3), we obtain a polynomial in F(ξ ). Setting the coefficients of the powers of F(ξ ) to zero, we obtain the following system of algebraic equations: Solving the above system of equations for α 0 , α 1 , β 1 , m and λ, we obtain the following values: Set-1: Set-2: Set-3: ab(8bk + 1) .

Family-4:
Family-5: Family-6: Case 2. When k > 0, we get the following trigonometric solutions. Family-7: Family-8: Family-9: Family-10: Family-11: Family-12: Case 3. When k = 0, we get the following solutions. Family-13: Family-14: Family-15: Family-16: Family-17: Remark. All of these solutions have been verified with MAPLE by substituting them into the original equations. In this section, we consider the mKdV equation given by where δ is a non-zero constant. The mKdV equation is similar to the KdV equation in that both are completely integrable and each has infinitely many conserved quantities. The mKdV equation appears in the study of electric circuits and multi-component plasmas [39,40]. We substitute the travelling wave transformation u(x, t) = u(ξ ), ξ = x + λt into equation (4.5) and obtain the ODE: Now integrating equation (4.6) with respect to ξ once and setting the constant of integration to zero, we obtain Taking the homogeneous balance between the highest order nonlinear term u 3 and the derivative term u from equation ( Hence for N = 1 equation (3.5) reduces to u(ξ ) = α 0 + α 1 (m + F(ξ )) + β 1 (m + F(ξ )) −1 . (4.8) Now substituting equation (4.8) into equation (4.7), we get a polynomial in F(ξ ). Setting the coefficients of the powers of F(ξ ) to zero, we obtain the following system of algebraic equations: Solving the above system of equations for α 0 , α 1 , β 1 , m and λ, we get the following values: Set-1: Set-2: Case 1. When k < 0, we get the following hyperbolic trigonometric solutions. Family-1: Family-2: Family-3: u 9,10 (ξ ) = ∓ √ 6δ −k tanh −kξ and u 11,12  Family-4: Case 2. When k > 0, we get the following trigonometric solutions.   Family-8: Case 3. When k = 0, we get the following solutions. Family-9: Family-10:   where ξ = x.
Remark. All of these solutions have been verified with MAPLE by substituting them into the original equations.