The exact solutions and approximate analytic solutions of the (2 + 1)-dimensional KP equation based on symmetry method.

In this paper, we successfully obtained the exact solutions and the approximate analytic solutions of the (2 + 1)-dimensional KP equation based on the Lie symmetry, the extended tanh method and the homotopy perturbation method. In first part, we obtained the symmetries of the (2 + 1)-dimensional KP equation based on the Wu-differential characteristic set algorithm and reduced it. In the second part, we constructed the abundant exact travelling wave solutions by using the extended tanh method. These solutions are expressed by the hyperbolic functions, the trigonometric functions and the rational functions respectively. It should be noted that when the parameters are taken as special values, some solitary wave solutions are derived from the hyperbolic function solutions. Finally, we apply the homotopy perturbation method to obtain the approximate analytic solutions based on four kinds of initial conditions.

Background

Recently, the nonlinear phenomenon has been extensively appeared in the fields of mathematical physics and engineering technology. A multitude of research focuses have been changed from linear problems to nonlinear ones. These problems can be ascribed to the research of nonlinear partial differential equations (NLPDE), as the complexity of equation, it becomes hard to get the exact solutions. Hence, the investigation of solving NLPDE has important theoretical and practical significance. In recent decades, a growing number of scholars established effective methods and these methods have obtained comprehensive applications. Such as the symmetry method (Lie 1881), the CK direct method (Clarkson and Kruskal 1989), the homogeneous balance method (Wang and Li 1996), the tanh-function method (Fan 2000; Akbar et al. 2013; Xie et al. 2005), the F-expansion method (Wang and Li 2005), sub-ODE method (Zhang et al. 2006), the simplest equation method (Kudryashov 2005; Sudao and Temuer 2010), the -expansion method (Wang et al. 2008; Alam and Akbar 2013), the homotopy perturbation method (He 1999; Narayanamoorthy and Sathiyapriya 2016; Filobello-Nino et al. 2016) and so on. However, a unified and systemic method, which can be applied to solve all kinds of equations, is still inexistence, and the above-mentioned methods all have individual range of applications. Therefore, summarizing and concluding, adopting the advantages and abandoning the disadvantages have been regarded as the effective approaches to investigate these existing methods. At the same time, it is worth obtaining more new solutions of NLPDE by using Lie symmetry and other methods.

As we all know, the symmetry method is the most universal method, and many traditional methods become its special cases. During the end of the nineteenth century, in order to unify and expand the methods used in solving the ordinary differential equations(ODE), Norwegian mathematician Sophus Lie (1842–1899) firstly proposed the symmetry theory of differential equations (Lie 1881). The investigations of the symmetry theory and approach have important theoretical and practical significance in modern mathematics, physics, mechanics and so on, at the same time, many successful applications have emerged in those fields (Bluman and Kumei 1989; Bluman et al. 2009; Noether 1918; Ma 1990; Clarkson and Kruskal 1989; Lou and Tang 2001; Ma and Chen 2009; Ma 2013). At present, using the symmetry method and others, such as the analytic solutions method, the approximate analytic solutions method and the numerical method with the aid of thorough considering mutual complementarity and availability to solve NLPDE are the new research subjects.

The premise of applying the symmetry method is to determine the all kinds of symmetries of the partial differential equations (PDEs). The main approach of determining the symmetries is the infinitesimal transform method which is proposed and constructed by Lie, called Lie’s algorithm. Lie’s algorithm, which is the major method with respect to determining symmetries, transforms the problem of determining symmetries into that of determining corresponding infinitesimal vectors whose infinitesimal functions are found as solutions of some over-determined system of PDEs, called the determining equations (Lie 1881). In determining symmetries, tedious, mechanical computations are involved and the order relation of unknown quantities have not been considered in conventional Lie’s algorithm, which result many problems, such as infinite loops on computers, a mass of work and so on. According to the investigations, differential form Wu’s method is one of effective methods to get rid of the defects of Lie’s algorithm. Therefore, Wu-differential characteristic set algorithm extended and constructed by Temuer Chaolu can partially solve the above-mentioned problems (Temuer 1999; Temuer and Bai 2010). This algorithm has been successfully applied to classical symmetries, nonclassical symmetries, high-order symmetries, approximate symmetries, potential symmetries, conservation laws and symmetry classification of PDEs, which has promoted the investigations of symmetry theory of PDEs (Bluman and Temuer 2006; Temuer et al. 2007; Temuer and Bai 2009; Temuer and Pang 2010; Sudao et al. 2014). Recently, we investigate the applications of the symmetry method in the boundary value problem of the nonlinear PDEs based on Wu-differential characteristic set algorithm and use the symmetry method and the homotopy analytic method to solve the boundary value problem (Sudao et al. 2014; Sudao 2011). Some other investigators use the symmetry method, the variational iterative method and the homotopy perturbation method to solve the boundary value problem based on Wu-differential characteristic set algorithm (Lu and Temuer 2011a, b; EerDun and Temuer 2012).

In this paper, we will construct the exact solutions and the approximate analytic solutions of the (2 + 1)-dimensional KP equation by using the Lie symmetry, the extended tanh method and the homotopy perturbation method. The Wu-differential characteristic set algorithm plays an important role in calculating the symmetries of the (2 + 1)-dimensional KP equation. This investigation will explore a new approaches of Lie symmetry in application of NLPDE. In addition, it will also effectively popularize the range of application and advance the efficiency of using method.

The symmetries and symmetry reduction of the (2 + 1)-dimensional KP equation

We consider the (2 + 1)-dimensional KP equation (Ding and Ji 2008) as follow:it is applied to describe the law of motion of water waves in (2 + 1)-dimensional spaces as well as plasmas in magnetic fields. Next, we will give the process of calculating the symmetry and reduction of Eq. (1).

The symmetries of the (2 + 1)-dimensional KP equation

The symmetry group of Eq. (1) will be generated by the vector field of the formwhere , , , are the infinitesimal generated functions of the symmetry. According to the Lie algorithm, we obtain the determining equations of symmetry (2), but it is too difficult to get its solutions. However, we can obtain the followig system of equations corresponding to the characteristic set which is equivalent to the determining equations by using Wu-differential characteristic set algorithm (Temuer 1999).

By solving the above PDEs, we get the infinitesimal functionswhere , , , are arbitrary constants, then the corresponding infinitesimal vector is the following form

Obviously, X has five one-parameter point symmetries, then the corresponding infinitesimal vectors are as follow:

The reduction of Eq. (1)

To facilitate solve the Eq. (1), we will reduce it by using the invariant form method. The resulting reduced PDE is fewer independent variable than Eq. (1).

Case 1

When  = , we obtain by solving the characteristic equation , where , are the invariants, then the reduction of Eq. (1) is

By the same token, we will get the following reductions.

Case 2

When  = , we reduce to Eq. (1) by using the invariant form method as follow:where , and , are the invariants.

Case 3

When  = , we reduce to Eq. (1) as follow:where and , are the invariants.

Case 4

When  = , we reduce to Eq. (1) as follow:where , and , are the invariants.

Case 5

When  = , we reduce to Eq. (1) as follow:where , and , are the invariants.

Case 6

When  = , we reduce to Eq. (1) as follow:where , and , are the invariants.

In all above, , , , such as , , and so on. From the above Eqs. (6) to (11), it is not difficult to find that Eq. (1) is reduced into the variable coefficient equations by using the symmetry and the constant coefficient equations by using the symmetries – respectively.

The exact travelling wave solutions of (8) based on the extended tanh method

Recently, as an effective approach, the extended tanh method is introduced to seek the exact solutions of the nonlinear evolution equations by Xie et al. (2005). This method is further improved by the generalized Riccati equation and introducing its twenty seven new solutions, these solutions are expressed by the hyperbolic functions, the trigonometric functions and the rational functions, respectively. When the parameters are taken as special values, some solitary wave solutions are derived from the hyperbolic function solutions.

Taking Eq. (8) for example from the symmetry reduction equations, we will get its exact travelling wave solutions by the extended tanh method and the process is composed of the following four steps.

Step 1

Doing the travelling wave transformations. In order to look for the travelling wave solutions of Eq. (8), we introduce the travelling wave transformation as follows:where k, c are constants and . Then we reduce Eq. (8) into ODE for , namely

Step 2

Choosing the expression of solution. By considering the homogeneous balance between the highest order derivatives and nonlinear terms appearing in Eq. (13), we choose the following expression of solution:where , , are undetermined coefficients. The function satisfies the second-order linear ODEwhere , , are constants. The ODE (15) has four cases of solutions as follows.

Case 1

When and (or ),where A and B are two nonzero constants and satisfies

Case 2

When and (or ),where A and B are two nonzero constants and satisfies for the above –, setting .

Case 3

When and where, is an arbitrary constant.

Case 4

When and where c is an arbitrary constant, and for the above –, setting

Step 3

Determining the coefficients. By substituting (14) into Eq. (13) and using ODE (15), collecting all terms with the same order of together, the left-hand side of Eq. (13) is converted into another polynomial in . Equating each coefficient of this different power terms to zero yields a set of nonlinear algebraic equations for and . With the aid of mathematica, we get the solutions as follows:

By analyzing (16), these solutions are suitable to all cases of the general solutions – to ODE (15).

Step 4

Acquiring the exact travelling wave solutions. By substituting (16) and the general solutions – of ODE (15) into (14) respectively, we obtain the exact travelling wave solutions as follows:the solutions (17) have 27 different cases, which are expressed by the hyperbolic functions, the trigonometric functions and the rational functions, respectively. The solitary wave solutions can be obtained (see Fig. 1) when the parameters are taken as special values.

Fig. 1

The solitary wave solutions of the exact solutions

The approximate analytic solutions of Eq. (8) based on the homotopy perturbation method

The homotopy perturbation method is proposed by He (1999), and it has successfully been applied to solve many types of linear and nonlinear functional equations. This method, which is a combination of homotopy in topology and classic perturbation techniques, provides us with a convenient way to obtain analytic or approximate solutions for a wide variety of problems arising indifferent fields. In recent years, the application of the homotopy perturbation method in nonlinear problems has been developed by scientists and engineers (He 2003, 2006; Olga 2011; Ebaid 2014; Najafi and Edalatpanah 2014).

Next, we construct the approximate analytic solutions of Eq. (8) by using the homotopy perturbation method. The is a function of in (17), and . Based on the solutions (17), we substitute into (17) and take , then the following initial conditions of the homotopy perturbation method can be obtained.

According to the homotopy perturbation method (He 1999), we construct the following homotopy

Equation (19) has the following form of solutionswhere p is an embedding parameter, and are undetermined. In order to be convenient for computing, we choose the following initial value approximation

By substituting (20) and (21) into Eq. (19) and collecting parameters with the aid of expansion as follows:where donates that takes derivative with respect to the variant. We choose the initial conditions as follows:and construct the approximate analytic solutions of Eq. (8) based on the following four cases.

When , satisfying the initial conditions as follows:the solutions can be obtained by (22) and (24) as follows:

Remark

Two variables , have are substituted in (25) and (26).

Then the second-order approximate solutions of Eq. (8) can be achieved by (25) and (26)

When , satisfying the initial conditions as follows:the solutions can be obtained by (22) and (28) as follows:

Then the second-order approximate solutions of Eq. (8) can also be achieved by (29) and (30)

When satisfying the initial conditions as follows:the solutions can be obtained by (22) and (32) as follows:

Then the second-order approximate solutions of Eq. (8) can also be achieved by (33) and (34)

When , satisfying the initial conditions as follows:the solutions can be obtained by (22) and (36) as follows:

Then the second-order approximate solutions of Eq. (8) can also be achieved by (37) and (38)

Figures 2 and 3 show the exact solutions (17) () and the second-order approximate solutions (31) based on the homotopy perturbation method of Eq. (8) respectively when the parameters are regarded as proper values. Table 1 shows the error comparison between the solutions (17) () and (27) when , , , , . According to the figure and table, the exact property of the homotopy perturbation method has been showed successfully.

Fig. 2

The exact solutions for , , , , ,

Fig. 3

The approximate analytic solutions for , , , , ,

Table 1

The error comparison between and at y=0.2

t	Error
	x
	0.01	0.02	0.03	0.04	0.05
0.05	2.70070 × 10−6	2.80067 × 10−6	2.90063 × 10−6	3.00058 × 10−6	3.10053 × 10−6
0.10	2.20078 × 10−6	2.30077 × 10−6	2.40076 × 10−6	2.50074 × 10−6	2.60072 × 10−6
0.15	1.70074 × 10−6	1.80076 × 10−6	1.90077 × 10−6	2.00078 × 10−6	2.10078 × 10−6
0.25	2.00189 × 10−7	3.00238 × 10−7	4.00286 × 10−7	5.00333 × 10−7	6.00378 × 10−7
0.30	1.60072 × 10−6	8.00465 × 10−7	9.00506 × 10−7	1.00054 × 10−6	1.10058 × 10−6

Conclusion

In this paper, we studied that construct the exact solutions and the approximate analytic solutions of NLPDE by using the Lie symmetry, the extended tanh method and the homotopy perturbation method. Specifically, we have constructed the abundant exact travelling wave solutions and approximate analytic solutions of the (2 + 1)-dimensional KP equation by using the above-mentioned three methods and obtained the high-precision approximate solutions by error analysis.

Lie symmetry, the extended tanh method and the homotopy perturbation method are effective methods which applied to solve PDEs. Hence, comprehensive use of them will advance their availability. The Wu-differential characteristic set algorithm is a key factor which influence the calculating the symmetry of PDEs. At present, combining the Wu-differential characteristic set algorithm, symmetry method and others to solve NLPDE has been regarded as a hot research topic and widened the application of symmetry and the Wu-differential characteristic set algorithm. This investigation is valuable in advanced research and development.