ABSTRACT

In the present article, we construct the exact traveling wave solutions of nonlinear PDEs in mathematical physics via the variant Boussinesq equations and the coupled KdV equations by using the extended mapping method and auxiliary equation method. This method is more powerful and will be used in further works to establish more entirely new solutions for other kinds of nonlinear partial differential equations arising in mathematical physics.

1. Introduction

The effort in finding exact solutions to nonlinear equations is important for the understanding of most nonlinear physical phenomena. For instance, the nonlinear wave phenomena observed in fluid dynamics, plasma and optical fibers are often modeled by the bell shaped sech solutions and the kink shaped tanh solutions. Many effective methods have been presented, such as inverse scattering transform method [1], Bäcklund transformation [2], Darboux transformation [3], Hirota bilinear method [4], variable separation approach [5], various tanh methods [6-9], homogeneous balance method [10], similarity reductions method [11,12], -expansion method [13], the reduction mKdV equation method [14], the tri-function method [15,16], the projective Riccati equation method [17], the Weierstrass elliptic function method [18], the Sine-Cosine method [19,20], the Jacobi elliptic function expansion [21,22], the complex hyperbolic function method [23], the truncated Painleve’ expansion [24], the F-expansion method [25], the rank analysis method [26], the ansatz method [27,28], the exp-function expansion method [29], the sub-ODE method [30], and so on.

The main objective of this paper is used the extended mapping method and auxiliary equation method to construct the exact solutions for nonlinear evolution equations in the mathematical physics via the variant Boussinesq equations and the coupled KdV equations.

2. Description of the Extended Mapping Method

Suppose we have the following nonlinear PDE:

(1)

where is an unknown function, F is a polynomial in and its various partial derivatives in which the highest order derivatives and nonlinear terms are involved. In the following we give the main steps of a deformation method:

Step 1. The traveling wave variable

(2)

where k and are the wave number and the wave speed, respectively. Under the transformation (2), Equation (1) becomes an ordinary differential equation (ODE) as

(3)

where

Step 2. Suppose that the solution Equation (3) has the following form:

(4)

where and are constants to be determined later, while satisfies the auxiliary ordinary differential equation:

(5)

where p, q and r are arbitrary constants.

Step 3. The positive integer “n” can be determined by considering the homogeneous balance between the highest derivative term and the nonlinear terms appearing in Equation (3). Therefore, we can get the value of n in Equation (4).

Step 4. Substituting Equation (4) into Equation (3) with the condition (5), we obtain polynomial in ,. Setting each coefficient of this polynomial to be zero, yields a set of algebraic equations for and k.

Step 5. Solving the algebraic equations by use of Maple or Mathematica, we have and k expressed by p, q, r.

Step 6. Since the general solutions of auxiliary Equation (5) have been well known for us (see Appendix A), then substituting the obtained coefficients and the general solution of Equation (5) into Equation (4), we have the travelling wave solutions of the nonlinear PDE (1).

3. Applications of the Method

In this section, we apply the extended mapping method to construct the exact solutions for the variant Boussinesq equations and the coupled KdV equations. Which are very important nonlinear evolution equations in mathematical physics and have been paid attention by many researchers.

3.1. Example 1. The Variant Boussinesq Equations

We start with the variant Boussinesq equations [31] in the following form:

(6)

(7)

where is a constant. As models for water waves, u is the velocity and v is the total depth. Wang [32] obtained their solitary wave solutions by using homogeneous balance method while Fan et al. [33] got a series of new traveling wave solutions of this system by using an algebraic method. The traveling wave variables below

(8)

permit us converting the Equations (6) and (7) into ODEs for and as follows:

(9)

and

(10)

On integrating Equations (9) and (10) with respect to once, we get

(11)

and

(12)

where L₁ and L₂ are integration constants.

Suppose that the solutions of Equations (11) and (12) can be expressed by

(13)

(14)

where and are constants to be determined later.

Considering the homogeneous balance between the highest order derivatives and the nonlinear terms in Equations (11) and (12), we get n = m = 2, hence the exact solutions of Equations (11) and (12) can be rewritten as:

(15)

and

(16)

where and are constants to be determined later. Substituting Equations (15) and (16) with the condition (5) into Equations (11) and (12) and collecting all terms with the same power of, . Setting each coefficients of this polynomial to be zero, we get a system of algebraic equations which can be solved by Maple or Mathematica. Thus we obtain the two sets of solutions as:

Case 1:

(17)

Case 2:

(18)

Note that, there are other cases which are omitted here. Since the solutions obtained here are so many, we just list some of the exact solutions corresponding to Case 2 to illustrate the effectiveness of the extended mapping method.

Substituting (18) into Equation (15) and Equation (16) yields

(19)

and

(20)

where

(21)

According to the Appendix A, we have the following families of exact solutions:

Family 1. If, , , , then we get

(22)

and

(23)

where

(24)

Family 2. If, , , , then we get

(25)

and

(26)

where

(27)

Family 3. If, , , then we get

(28)

and

(29)

where

(30)

Family 4. If, , , , then we get

(31)

and

(32)

where

(33)

Family 5. If, , , then we get

(34)

and

(35)

where

(36)

Family 6. If, , ,

, then we get

(37)

and

(38)

where

(39)

Similarly, we can write down the other families of exact solutions of Equations (6) and (7) which are omitted for convenience.

3.2. Example 2. The Coupled KdV Equations

In this subsection, we consider the coupled KdV Equations [34]:

(40)

(41)

where is a constant. The traveling wave variables (8) permit us converting the Equations (40) and (41) into ODEs in the forms:

(42)

(43)

Suppose that the solutions of Equations (42) and (43) can be expressed by

(44)

(45)

where and are constants to be determined later.

Considering the homogeneous balance between the highest order derivatives and the nonlinear terms in Equation (63) and Equation (64), we get n = m = 2, so Equations (44) and (45) can be rewritten as:

(46)

(47)

where and are constants to be determined later. Substituting Equation (46) and Equation (47) with the condition (5) into Equation (42) and Equation (43) and collecting all terms with the same power of, . Setting each coefficients of this polynomial to be zero, we get a system of algebraic equations which can be solved by Maple or Mathematica to get the following sets of solutions:

Case 1:

(48)

Case 2:

(49)

Note that, there are other cases which are omitted here. Since the solutions obtained here are so many, we just list some of the exact solutions corresponding to Case 2 to illustrate the effectiveness of the extended mapping method.

Substituting (49) into Equation (46) and Equation (47) yields

(50)

and

(51)

where

(52)

According to the Appendix A, we have the following families of exact solutions:

Family 1. If, , , , then we get

(53)

and

(54)

where

(55)

Family 2. If, , , , then we get

(56)

and

(57)

where

(58)

Family 3. If, , ,

, then we get

(59)

and

(60)

where

(61)

Family 4. If, , , , then we get

(62)

and

(63)

where

(64)

Family 5. If, , , , then we get

(65)

and

(66)

where

(67)

Family 6. If, , , , then we get

(68)

and

(69)

where

(70)

4. Conclusion and Discussion

In this article, we have found the exact solutions of the variant Boussinesq equations and the coupled KdV equations by using the extended mapping method and the auxiliary equation method. Also, we conclude according to the Appendix B that our results in terms of Jacobi elliptic functions generate into hyperbolic functions when m ® 1 and generate into trigonometric functions when m ® 0. This method provides a powerful mathematical tool to obtain more general exact solutions of a great many nonlinear PDEs in mathematical physics.

REFERENCES

Appendix A

The general solutions of the auxiliary Equation (5) and its derivatives [35] are listed as follows:

where 0 < m < 1 is the modulus of the Jacobi elliptic functions and.

Appendix B

The Jacobi elliptic functions generate into hyperbolic functions when as follows:

and into trigonometric functions when as follows:

Appendix C