A Generalized Simplest Equation Method and Its Application to the Boussinesq-Burgers Equation

In this paper, a generalized simplest equation method is proposed to seek exact solutions of nonlinear evolution equations (NLEEs). In the method, we chose a solution expression with a variable coefficient and a variable coefficient ordinary differential auxiliary equation. This method can yield a Bäcklund transformation between NLEEs and a related constraint equation. By dealing with the constraint equation, we can derive infinite number of exact solutions for NLEEs. These solutions include the traveling wave solutions, non-traveling wave solutions, multi-soliton solutions, rational solutions, and other types of solutions. As applications, we obtained wide classes of exact solutions for the Boussinesq-Burgers equation by using the generalized simplest equation method.


Introduction
It is important to seek more exact solutions of nonlinear evolution equations (NLEEs) in mathematical physics. In past decades, many powerful methods have been presented such as the inverse scattering method [1], the Darboux transformation [2], the Hirota bilinear method [3], the Painlevé expansion method [4], the Bäcklund transformation method [5,6], the multilinear variable separation method [7], the homogeneous balance method [8], the Jacobi elliptic function expansion method [9], the tanh-function method [10,11], the F-expansion method [12], the auxiliary equation method [13], the sub-ODE method [14], the Exp-function method [15], the (G 0 /G)-expansion method [16,17], the simplest equation method [18,19]. Thousands of examples have shown that these methods are powerful for obtaining exact solutions of NLEEs, especially, for traveling wave solutions. Up to now, a unified method that can be used to deal with all types of NLEEs has not been discovered. Hence, developing new method and finding more general exact solutions of NLEEs have drawn a lot of interests of a diverse group of scientists.
Particularly, N.A. Kudryashov first proposed the simplest equation method and showed that it is powerful for finding analytic solutions of NLEEs [18,19]. Two basic ideas are at the heart of the approach. The first idea is to apply the simplest nonlinear differential equations (the Riccati equation, the equation for the Jacobi elliptic faction, the equation for the Weierstrass ellipic function and so on) that have lesser order than the equation studied. The second idea is to use possible singularities of general solution for the equation studied. Recently, there are many applications and generalizations of the method [20][21][22][23][24][25].
In the present paper, based on a solution expression with a variable coefficient and variable coefficient auxiliary equation (a special simplest nonlinear differential equation), we proposed a generalized simplest equation method to seek exact solutions of NLEEs. This method can yield infinite number of exact solutions of NLEEs. For illustration, we apply this method to the Boussinesq-Burgers equation and successfully find many different types of exact solutions. These solutions include traveling wave solutions, non-traveling wave solutions, multi-soliton solutions, rational solutions and other types of solutions.
The rest of the paper is organized as follows. In Section 2, we describe the generalized simplest equation method to look for exact solutions of NLEEs. In Section 3, we firstly construct the most general form exact solutions and the exact traveling wave solutions of the Boussinesq-Burgers equation by using our method. Secondly, we obtain the Bäcklund transformation between the Boussinesq-Burgers equation and the constraint equation. Thirdly, we give some concrete new exact solutions of the Boussinesq-Burgers equation. In Section 4, we give some discussions and conclusion remarks.

Description of the generalized simplest equation method
We consider a nonlinear evolution equation (NLEE) where u = u(t, x) is an unknown function, t, x are two independent variables, P is a polynomial in u and its various partial derivatives. We consider a generalized simplest equation method for obtaining exact solutions to the given NLEE Eq (1), which is described in the following five steps.
Step 3. By substituting Eq (2) into Eq (1) and using Eq (3), collecting all terms with the same order of ϕ together, the left-hand side of Eq (1) is converted into another polynomial in ϕ. Equating each coefficient of powers of ϕ to zero, yields a set of over-determined partial differential equations for α i (t, x)(i = 0, 1, Á Á Á, M) and ξ(t, x).
Remark 2: It is the first time that Eq (3) is used as an auxiliary equation for obtaining exact solutions of NLEEs.
Remark 3: It is emphasized that in our designs Eqs (2) and (3), the parameters δ, ν will play the role of adjusters in solving the constraint equation and in obtaining a transformation between Eq (1) and linear heat equation (see Section 3).

Application of the Method
The exact solutions to the Boussinesq-Burgers equation In this section, we are aimed to first give the most general form exact solutions, then we will determine the exact traveling wave and non-traveling wave solutions for the Boussinesq-Burgers equation [26][27][28][29][30][31][32] The general form exact solutions. In the subsection, we solve the most general form exact solutions of Eqs (8) and (9).
Step 2. Considering the homogeneous balance between the highest order derivatives v x (u xxx ) and nonlinear terms uu x (uv x ) in Eqs (8) and (9), we get balance numbers M = 1 for u and N = 2 for v. Thus, we can denote the solutions of Eqs (8) and (9)  where functions a 0 (t, x), a 1 (t, x), b 0 (t, x), b 1 (t, x), b 2 (t, x) and ξ = ξ(t, x) are to be determined later and ϕ = ϕ(ξ) satisfies Eq (3).
Step 3. We substitute Eqs (10) and (11) along with the auxiliary Eq (3) into Eqs (8) and (9) and collect all terms with the same order of ϕ. As a result, the left-hand sides of Eqs (8) and (9) are converted into another polynomials in ϕ. Equating each coefficient of powers of ϕ to zero, we obtain a set of over-determined partial differential equations for a 0 (t, x), a 1 Remark 4: Unlike pure algebraic equations used in the standard auxiliary method, here we obtain differential equations for a i (t, x), b j (t, x) and ξ(t, x). This provides us with more general types of exact solutions to Eqs (8) and (9).
Step 4. Solving the differential system obtained in Step 3 by Mathematica, we obtain a 0 (t, x), x) and δ, ν as follows: where ξ = ξ(t, x) satisfies the equation Eq (13) is called the constraint equation of ξ.
Here for obtaining more general types of solutions, we consider the case of ξ x 6 ¼ 0. Any solution of the constraint Eq (13) leads to a group of coefficients Eq (12) which together with Eqs (4)-(6) result in three classes of exact solutions Eqs (10) and (11) of Eqs (8) and (9) as follows: ( in which i = 1 for s > 1, i = 2 for s < 1 and i = 3 for s = 1. Therefore, the solution expressions Eq (14) have established a Bäcklund transformation between Eqs (8) and (9) and constraint Eq (13). Because of the Bäcklund transformation, the method can give infinite number of exact solutions for the considering NLEEs immediately.
Seen from the formulae Eq (14), if ξ(t, x) is a solution of Eq (13) with the form ξ = ξ(x−Vt), then it yields exact traveling wave solutions of Eqs (8) and (9). While if ξ 6 ¼ ξ(x−Vt), then it yields exact non-traveling wave solutions of Eqs (8) and (9). Hence the formulae Eq (14) provide us with abundance of general form exact solutions to Eqs (8) and (9) once the solutions of Eq (13) are given.
Step 5. In the following, we determine the exact traveling wave and non-traveling wave solutions of Eqs (8) and (9) through exactly solving Eq (13).
Case 3: s = 1. In this case, there are four kinds of solutions of Eq (13) shown in the following subcases.
Substituting the obtained solutions ξ = ξ(x−Vt) in all of above cases Eqs (18)- (29) into Eqs (4)-(6) and (12) and assembling them in Eq (14), we obtain twelve kinds of exact traveling wave solutions of Eqs (8) and (9). Here we omit the duplicate and long expressions of these solutions.
Remark 5: A feature of these solutions is that they are multi-composite solutions, i.e., they are constructed by compounding several elementary functions. For example, in the case s > 1, the solutions Eqs (18)- (20) are the compounding of five elementary functions exp, arctan, tanh, tan, arccos etc. These solutions can not be obtained by tanh-function method, the (G 0 /G)expansion method and other method given [6] etc. The solution Eqs (21), (22), (26) and (27) have included soliton solutions when the parameters are taken properly. The solutions Eqs (24) and (25) are rational or irrational solutions at some values of these parameters appeared in these solutions. These already significantly expand the set of exact solutions for Eqs (8) and (9).

A transformation between Eqs (8) and (9) and a linear heat equation
It is observed that if the t-derivative is equal to second order x-derivative, then each monomial term in Eq (13) has the same sum of x-derivative orders. This hints a relationship between Eq (13) and a kind of heat equation. In fact, Eq (13) can be written as where P = ξ t −τξ xx −aξ x for arbitrary constant a, and In order to get most general form exact solutions, we consider the following two cases. Case I. For xx xx À x 2 x ¼ 0 and s = 2δ−δ 2 +4ν. In this case, by solving Eq (13) we have ξ(t, x) = d 4 exp[(d 2 x+d 3 )/(t+d 1 )], which results in exact solutions of Eqs (8) and (9) given by Eq (14) with Eq (12) as follows in which i = 1 for s > 1, i = 2 for s < 1 and i = 3 for s = 1. ϕ 1 (ξ), ϕ 2 (ξ) and ϕ 3 (ξ) are given in Eqs (4)- (6). d 1 , d 2 , d 3 , d 4 are arbitrary constants.
Case II. For s = 2δ−δ 2 +4ν = 0. In this case, Eq (13) or Eq (30) admits the solutions ξ(t, x) of the linear heat equation with a conductive term aξ x for arbitrary constant a.
It is interesting that a transformation in the form of Eq (14) between Eqs (8) and (9) and the linear heat Eq (35) has been found. Using the transformation, we will obtain infinite number of exact solutions of Eqs (8) and (9). That means each solution of the heat Eq (35) (such as those given in the following Eq (34)) yields a set of exact solutions of Eqs (8) and (9). Moreover, the connection provides us with an insight into integrability of Eqs (8) and (9). For example, using the transformation, we can investigate the multi-soliton solutions, rational solutions and other types of solutions to the Boussinesq-Burgers equation.
As we know, the heat Eq (35) has infinite number of specific known exact solutions. Some representatives of them are listed as follows: xðt; xÞ ¼ A exp ðÀtk 2 tÞ cos ðkðx þ atÞ þ BÞ þ C; ðt; xÞ 2 R 2 ; xðt; xÞ ¼ A exp ðÀkðx þ atÞÞ cos ðkx þ kða À 2tkÞt þ BÞ þ C; ðt; xÞ 2 R 2 ; where A, B, C, c 0 and k are arbitrary constants, ω(θ) is any integrable function, erfc(y) is error function and c n , n = 1, 2, Á Á Á are constants given by the initial or boundary conditions of Eq (35). τ = ±1/2 is given in Eq (31). It is noticed that the last solution in Eq (34) yields an exact solution to the Boussinesq-Burgers equation with containing an arbitrary function, which may give more freedom to solve related problems of the equations [7].
constraint equation, we can derive infinite number of exact solutions for NLEEs. These solutions include the traveling wave solutions, non-traveling wave solutions, multi-soliton solutions, rational solutions and other types of solutions for NLEEs. According to the auxiliary equation and above two features of our method, we call the method as the generalized simplest equation method. As application of our method, a great number of new exact solutions of the Boussinesq-Burgers equation are obtained. Essentially, the constraint equation is related to linear heat equation which result in new exact solutions of the Boussinesq-Burgers equation. These solutions are different from those solutins given in literatures [26][27][28][29][30][31][32]. It shows that our method is more flexible in finding more general form exact solutions and the method can be used for many other NLEEs in mathematical physics. As a result, our method effectively enhances the existing auxiliary equation methods, such as those given in [10,16], which demonstrates that the proposed method is effective and prospective.
We have also tested other several commonly used auxiliary equations, such as Riccati equation [10,11], the auxiliary equation of F-expansion method [12], to get exact solutions of the Boussinesq-Burgers equation in the frame of our method. Unfortunately, we can not get the results presented in the paper. What's more, the relationship between the Boussinesq-Burgers equation and the linear heat equation can not be found. Hence the auxiliary Eq (3) has its own advantages which are not admitted by other auxiliary equations.
In addition, the true idea of our method makes one further consider finding more general exact solutions of other NLEEs by developing the existing methods, such as Tanh-method, the simplest equation method, etc, or designing new method, in which a connection between the constraint equation and other easily solved equation may be set up. It is an excellent topic for further research.