Symmetry reduction and exact solutions of two higher-dimensional nonlinear evolution equations.

In this paper, symmetries and symmetry reduction of two higher-dimensional nonlinear evolution equations (NLEEs) are obtained by Lie group method. These NLEEs play an important role in nonlinear sciences. We derive exact solutions to these NLEEs via the [Formula: see text]-expansion method and complex method. Five types of explicit function solutions are constructed, which are rational, exponential, trigonometric, hyperbolic and elliptic function solutions of the variables in the considered equations.

Introduction

In 1998, Yu et al. [1] extended the Bogoyavlenskii Schiff equation to the -dimensional NLEE

Setting , equation (2) is changed into the -dimensional potential Yu-Toda-Sasa-Fukuyama (YTSF) equation

The generalized -dimensional Zakharov-Kuznetsov (gZK) equation is given by Here () are arbitrary constants.

We note that equation (4) includes many famous NLEEs as its special cases. For instance, if , then equation (4) is the Korteweg-de Vries equation [2, 3]. If , then equation (4) is the dimensional ZK-MEW equation [4]. If , then equation (4) is the Gardner equation [5]. If , then equation (4) is the modified Zakharov-Kuznetsov equation [6].

In recent years, it has aroused widespread interest in the study of NLEEs [7-13]. Equations (3) and (4) are very meaningful higher-dimensional NLEEs which can describe many dynamic processes and important phenomena in engineering and physics. The YTSF equation is a mostly used model for investigating the dynamics of solitons and nonlinear waves in fluid dynamics, plasma physics and weakly dispersive media [13]. Zakharov and Kuznetsov [14] proposed the ZK equation to describe nonlinear ion-acoustic waves in a plasma comprised of cold ions and hot isothermal electrons in the presence of a uniform magnetic field. Many physical phenomena, in the purely dispersive limit, are governed by this type of equation, such as the long waves on a thin liquid film [15], the Rossby waves in a rotating atmosphere [16], and the isolated vortex of drift waves in a three-dimensional plasma [17]. The gZK equation is of a generalized setting of ZK equation. Seeking exact solutions of NLEEs is an interesting and significant subject. Over the past few years, many powerful methods for constructing the solutions of NLEEs have been used, for instance, the Bäcklund transform method [18], direct algebraic method [19], modified simple equation method [20], Lie group method [21, 22], -expansion method [8, 9, 23, 24], and so on. Recently, Yuan et al. [25-27] introduced the complex method to find the exact solutions of NLEEs in mathematical physics. In this paper, we study symmetries, symmetry reduction of the two higher-dimensional NLEEs, and then we obtain their exact solutions via the -expansion method and complex method.

Description of the methods

Description of the -expansion method

Suppose that a nonlinear partial differential equation (PDE) is given by where P is a polynomial of an unknown function and its derivatives in which nonlinear terms and highest order derivatives are involved. The main steps of this method are given in the following.

Step 1. Substituting the traveling wave transform into equation (5) converts it to the following ordinary differential equation (ODE): in which F is a polynomial of and its derivatives, while .

Step 2. Assume that equation (6) has the following traveling wave solution: where () are constants to be determined, such that and satisfies the ODE as follows: Equation (8) has the following solutions.

When , , When , , When , , , When , , , When , , , Here , δ, μ are constants that will be determined later and c is an arbitrary constant. We take the homogeneous balance between nonlinear terms and highest order derivatives of equation (6) to determine the positive integer n.

Step 3. Substituting equation (7) into equation (6) and accounting the function , we obtain a polynomial of . Equating all the coefficients of the same power of to zero yields a set of algebraic equations. By solving the algebraic equations, we get the values of , δ, μ, and then we substitute them into equation (7) along with equations (9)-(15) to complete the determination of the solutions of equation (5).

Description of the complex method

Let , , , , and then is the degree of . Let the differential polynomial be defined by where J is a finite index set, and are constants, then is the degree of .

Consider the following differential equation: where , , d are constants.

Set , and the meromorphic solutions w of equation (16) have at least one pole. If equation (16) has exactly p distinct meromorphic solutions, and their multiplicity of the pole at is q, then equation (16) is said to satisfy the condition. It might not be easy to show that the condition of equation (16) holds, so we need the weak condition as follows.

Inserting the Laurent series into equation (16), we can determine exactly p different Laurent singular parts: then equation (16) is said to satisfy the weak condition.

Given two complex numbers , such that , and let L be the discrete subset , and L is isomorphic to . Let the discriminant and

A meromorphic function with double periods , , which satisfies the following equation: in which , , and , is called the Weierstrass elliptic function and satisfies an addition formula [28] as follows:

If a meromorphic function g is a rational function of z, or a rational function of , , or an elliptic function, then we say that g belongs to the class W.

In 2009, Eremenko et al. [29] studied the mth-order Briot-Bouquet equation (BBEq) where are constant coefficients polynomials, . For the mth-order BBEq, we have the following lemma.

Lemma 2.1

([28, 30, 31])

Let , , and a mth-order BBEq satisfies the weak condition, then the meromorphic solutions . Supposing for some values of the parameters the solution w exists, then any other meromorphic solutions will be one parameter family , . In addition, every elliptic solution w with a pole at is expressed as where are determined by (17), and .

Every rational function solution is expressed as which has h (≤p) distinct poles of multiplicity q.

Every simply periodic solution is a rational function of (), and is expressed as which has h (≤p) distinct poles of multiplicity q.

By the above definitions and lemma, we now present the complex method.

Step 1. Insert the transformation defined by into a given PDE to yield a nonlinear ODE.

Step 2. Insert (17) into the ODE to determine whether the weak condition holds.

Step 3. Insert the indeterminate solutions introduced in Lemma 2.1 into the ODE, and then get meromorphic solutions of the ODE with a pole at .

Step 4. Obtain meromorphic solutions by Lemma 2.1 and the addition formula.

Step 5. Inserting the inverse transformation into the meromorphic solutions, we get the exact solutions for the original PDE.

Symmetries and symmetry reduction

Symmetries

In order to find the symmetry of equation (4), we set where u is the solution of equation (4), a, b, c, d, e, f are unknown functions of real variables x, y, s, t. According to Lie group analysis [21, 22], σ satisfies Substituting equation (21) into equation (22), we have a new differential equation, where By equation (21), equation (22) and equation (23), we have where () are real constants. Substituting equations (24) into equation (21), we achieve the symmetry of the gZK equation,

In order to find the symmetry of equation (3), we set Here u is the solution of equation (3), a, b, c, d, e, f are unknown functions of real variables x, y, s, t. According to Lie group analysis, σ satisfies Substituting equation (26) into equation (27), we have a new differential equation, where By equation (26), equation (27) and equation (28), we have where () are real constants, , , are arbitrary real functions of t. Substituting equations (29) into equation (26), we achieved the symmetry of YTSF equation

Symmetry reduction

By solving the characteristic equation (25) of σ we find different symmetry reduced equations. Without loss of generality, we have two reduced equations as follows.

Setting , , we have the first similarity solution of equation (4) where , . Substituting equation (32) into equation (4), we have the first symmetry reduced equation of equation (4)

Setting , , solving , we have the second similarity solution of equation (4) where , . Substituting equation (34) into equation (4), we have the second symmetry reduced equation of equation (4)

By solving the characteristic equation (30) of σ we obtain symmetry reduction of equation (3). Without loss of generality, we have two reduced equations as follows.

Setting , , , solving , we have the first similarity solution of equation (3) where , . Substituting equation (36) into equation (3), we have the first symmetry reduced equation of equation (3)

Setting , , , solving , we have the second similarity solution of equation (3) Substituting equation (38) into equation (3), we have the second symmetry reduced equation of equation (3)

Exact solutions

Exact solutions of gZK equation via the -expansion method

Substituting the traveling wave transform into equation (33), then integrating it with respect to z, we obtain where γ is the integration constant which can be determined later.

Taking the homogeneous balance between and in equation (39) yields where , are constants to be determined, whereas δ and μ are arbitrary constants.

Substitute w, , , into equation (39) and equate the coefficients of to zero, then Solving the above algebraic equations, we obtain where μ and δ are arbitrary constants.

Substituting equations (41) into equation (40) yields We apply equation (9) to equation (15) into equation (42), respectively, to get traveling wave solutions of the gZK equation as follows.

When , ,

When , ,

When , , ,

When , , ,

When , , ,

Exact solutions of gZK equation via the complex method

Inserting (17) into equation (39) we have , , , , , and is an arbitrary constant.

Therefore, equation (39) is a second order BBEq and satisfies the weak condition. Hence, by Lemma 2.1, we see that meromorphic solutions of equation (39) belong to W. We will show meromorphic solutions of equation (39) in the following.

By (19), we infer that the indeterminate rational solutions of equation (39) are with a pole at .

Substituting into equation (39), we have where and ; where and .

So the rational solutions of equation (39) are and where , . , in the former case, or , in the latter case.

To obtain simply periodic solutions, let , and substitute into equation (39), then

Substituting into equation (43), we obtain and where , in the former case, or , in the latter case.

Inserting into equation (44) and equation (45), we can get simply periodic solutions to equation (39) with a pole at , where , in the former case, or , in the latter case.

So simply periodic solutions of equation (39) are and where , . , in the former case, or , in the latter case.

From (18), we have the indeterminate relations for the elliptic solutions of equation (39) with a pole at , where . Considering the results obtained above, we infer that , , . So we obtain where .

Thus, the elliptic function solutions of equation (39) are where , , is arbitrary. Applying the addition formula, we can rewrite it as where , , E and are arbitrary.

Exact solutions of YTSF equation via the -expansion method

Substituting the traveling wave transform into equation (37), then integrating it with respect to z, we obtain where γ is the integration constant which can be determine later.

Setting , equation (46) becomes

Taking the homogeneous balance between and in equation (47) yields where , () are constants to be determined, whereas δ and μ are arbitrary constants.

Substitute w, , into equation (47) and equate the coefficients of to zero, then

Solving the above algebraic equations, we obtain where μ and δ are arbitrary constants.

Substituting equations (49) into equation (48), yields We apply equation (9) to equation (15) into equation (50), respectively, to get traveling wave solutions of the YTSF equation as follows.

When , ,

When , ,

When , , ,

When , , ,

When , , ,

Exact solutions of YTSF equation via the complex method

Inserting (17) into equation (47) we have , , , , , , , and is an arbitrary constant.

Therefore, equation (47) is a second order BBEq and satisfies the weak condition. Hence, by Lemma 2.1, we see that meromorphic solutions of equation (47) belong to W. We will show meromorphic solutions of equation (47) in the following.

By (19), we deduce the indeterminacy rational solutions of equation (47) are with a pole at .

Substituting into equation (47), we get the following form: where .

So the rational solutions of equation (47) are where , .

To obtain simply periodic solutions, let , and substitute into equation (47), then we get Substituting into equation (51), we obtain where . Substituting into equation (52), we can obtain simply periodic solutions of equation (47), with a pole at .

Therefore the simply periodic solutions of equation (47) are where , .

From (18), we can express the elliptic solutions of equation (47) as with a pole at .

Substituting into equation (47), we obtain where , is arbitrary.

Therefore, the elliptic solutions of equation (47) are in which . Applying the addition formula, we can rewrite it as where , , is arbitrary.

Comparison

Implementing the -expansion method, we found seven solutions for the gZK and YSFT equation, respectively. Using the complex method, we found five solutions for the gZK equation and three solutions for the YSFT equation. Rational solutions and are obtained via the -expansion method, and and are obtained via the complex method. If we let , then is equivalent to , and is equivalent to . For getting rational solutions, these two methods are in good agreement. Rational solutions and simply periodic solutions and are new and cannot be degenerated successively through elliptic function solutions. From the results, we can find more solutions by the -expansion method, whereas we can obtain elliptic function solutions just by the complex method. These two methods are very useful tools in finding the exact solutions of NLEEs.

Computer simulations

In this section, we illustrate some results by the computer simulations. We carry out further analysis to the properties of simply periodic solutions and as in Figures 1 and 2.

Figure 1

The solution of the gZK equation corresponding to . (a) , (b) , (c) .

Figure 2

The solution of the YSTF equation corresponding to . (a) , (b) , (c) .

By employing the complex method, we are capable to obtain simply periodic solutions and of the gZK equation. The solutions and come from hyperbolic function. Figure 1 shows the shape of solutions for , , , , , , , and within the interval . Note that they have two distinct generation poles which are showed by Figure 1.

By using the complex method, we achieve to obtain simply periodic solutions of the YSTF equation. The solutions are in terms of the hyperbolic function solution. The solutions in Figure 2 of the YSTF equation are represented the singular soliton solution for the parameters , , , and within the interval .

Conclusions

In this article, we utilize Lie group analysis to obtain symmetries and symmetry reduction for two higher-dimensional NLEEs. In this way, we can reduce the dimension of the NLEEs, which is relevant in the fields of mathematical physics and engineering. Five types of explicit function solutions are constructed by the -expansion method and complex method. It demonstrates these methods are very efficient and powerful to seek the exact solutions of NLEEs. We can apply the idea of this study to other NLEEs.