Lie symmetry analysis , Lie-Bäcklund symmetries , explicit solutions , and conservation laws of Drinfeld-Sokolov-Wilson system

*Correspondence: zhangyfcumt@163.com 1College of Mathematics, China University of Mining and Technology, Xuzhou, 221116, P.R. China Full list of author information is available at the end of the article Abstract The symmetry analysis method is used to study the Drinfeld-Sokolov-Wilson system. The Lie point symmetries of this system are obtained. An optimal system of one-dimensional subalgebras is derived by using Ibragimov’s method. Based on the optimal system, similarity reductions and explicit solutions of the system are presented. The Lie-Bäcklund symmetry generators are also investigated. Furthermore, the method of constructing conservation laws of nonlinear partial differential equations with the aid of a new conservation theorem associated with Lie-Bäcklund symmetries is presented. Conservation laws of the Drinfeld-Sokolov-Wilson system are constructed by using this method.


Introduction
The Lie symmetry method was initiated by Lie [] in the second half of the nineteenth century. It has become one of the most powerful methods to study nonlinear partial differential equations (NLPDEs). The core idea of the Lie symmetry method is that one-parameter groups of transformations acting on the space of the independent and dependent variables leave the NLPDEs unchanged [-]. Most problems in science and engineering can be represented by NLPDEs [-]. Application of the Lie symmetry method for constructing the explicit solutions of the NLPDEs can always be regarded as one of the most important fields of mathematical physics. Many important properties of NLPDEs such as symmetry reductions, conservation laws, and explicit solutions by using symmetries can be considered successively [-].
Studying conservation laws is helpful in analyzing NLPDEs in the physical point of view []. Some studies have indicated that conservation laws play an important role in the numerical integration of PDEs []. Conservation laws are also helpful in solving equations by means of the double reduction theory [, ]. In order to construct conservation laws various methods have been developed, such as Noether's theorem [], the partial Noether approach [], the multiplier approach [], and a new conservation theorem [, ]. Noether's theorem and the partial Noether approach establish a relationship between symmetries and conservation laws for NLPDEs. Nevertheless, these methods do not work on a nonlinear equation without a Lagrangian. It is notable that a new conservation theorem was proposed by Ibragimov []. The conservation laws of a nonlinear equation and its adjoint equation can be constructed by using formal Lagrangian.
In the present paper, we shall consider the Drinfeld-Sokolov-Wilson system (DSWS) where a, b, and k are real valued constants. System () is a physical model for describing nonlinear surface gravity waves propagating over a horizontal sea bed. The DSWS has attracted the attention of many scholars. Some conservation laws of DSWS were obtained with the aid of the multiplier approach []. Then the double reduction analysis was employed to study the reductions of some conservation laws for DSWS []. This system has an infinite number of conservation laws, has a Lax representation, and is a member of KP hierarchy [, ], which indicates the integrality of this system. The scaling invariant Lax pairs of the DSWS were derived by Hickman and Hereman et al. []. The homotopy analysis method was applied to obtain the approximate solutions of DSWS []. Matjila et al. derived the exact solutions of system () by using the (G /G)-expansion function method. They also constructed conservation laws using Noether's approach []. We gave some symmetry reductions and conservation laws of this system [].
This paper is arranged as follows. In Section , we derive the Lie point symmetries of the DSWS using Lie group analysis and find the transformed solutions. In Section , a new optimal system of subalgebras of system () is constructed by using a concise method. The new optimal system contains five operators. Then in Section , based on the optimal system, the similarity reduced equations and the explicit solutions of system () are investigated systematically. In Section , the method of constructing conservation laws of NLPDEs with the aid of the new conservation theorem associated with Lie-Bäcklund symmetries is presented. In Section , the conservation laws of the DSWS are constructed. Finally, the conclusions are given in the last section.

Lie point symmetries
The infinitesimal generators will give rise to the Lie group of symmetries []. Then X should satisfy the following invariant surface conditions: where  = u t + vv x and  = v tav xxx + bu x v + kuv x . The invariant surface conditions give the overdetermined system of PDEs By solving (), we get where c  , c  , and c  are arbitrary constants. Therefore we obtain the three infinitesimal symmetry generators () The above three infinitesimal symmetry generators form a three-dimensional Lie algebra L  by commutation.
One can obtain the group transformation generated by the Lie point symmetry operator X i (i = , . . . , ) by solving the ordinary differential equations with initial conditions Then the one-parameter transformation groups G i of the DSWS are given as follows: Therefore, the following theorem is established. ()

Optimal system of subalgebras
Over the past few decades, the problem of classifying the subgroup was studied by many researchers. In order to demonstrate the classification of the group invariant solutions, Ovsiannikov put forward the concept of optimal systems of subalgebras of a . This method only relies on a commutator table. In this section, we shall use the Ibragimov method to construct an optimal system of one-dimensional subalgebra of the Lie algebra L  for the DSWS. An arbitrary infinitesimal operator X of Lie algebra L  can be written The transformations of the symmetry group with the Lie algebra L  provide the threeparameter group of linear transformations. The following generators are useful in finding these transformations, where c ij is determined by [X i , X j ] = c λ ij X λ . The commutator table of generators X  , X  , and X  is given in Table . Using equation () and commutator Table , E  , E  , and E  can be written For the generator E  , the Lie equations with the initial conditionl| a  = = l are written where a  is the parameter. The solutions of the above initial value problem give the transformatioñ Other Lie equations for E  and E  are where a  and a  are the parameters. In a similar way, we obtain the following transformations:l Theorem  The following five operators constitute an optimal system of one-dimensional subalgebras of the three-dimensional Lie algebra L  : Proof The process of constructing optimal systems is equivalent to the classification of the vector by means of the transformations (), (), and (). We focus on finding the simplest representatives of each class of similar vectors (). Case . l  = . First of all, we take a  = -l  l  in () and reduce the vector () to the form By taking a  = -l  l  in the transformations (), we makel  = . Hence the vector () is reduced to the form Thus this case provides the operator X  . Case . l  = .
(I) l  = . We assume l  =  and make l  = ± by the transformations (). In addition, taking into account the possibility l  = , we obtain the following representatives for the optimal system (II) Let l  = . If l  =  we can set l  =  and obtain the vector (, , ). This gives rise to the operator X  .

Symmetry reductions and explicit solutions of the DSWS
In Section , an optimal system of one-dimensional subalgebras of the DSWS has been obtained. We now make use of the symmetries of the optimal system () to reduce system () to a system of nonlinear ordinary differential equations in one variable. All the symmetry reductions are presented in Table . By solving the reduced equations, the group invariant solutions to the DSWS can be obtained.

Explicit solutions of system (1) using the simplest equation method
The method of the simplest equation was first proposed by Kudryashov where m and n are positive integers that can be determined by a homogeneous balance procedure and G(ξ ) satisfies the Bernoulli equation The Bernoulli equation () has the solution .
We get m = , n =  by balancing system (B). Hence the solution of (B) can be written () Substituting system () along with equation () into (B) and collecting the coefficient of the same power G i leads to a set of constraining equations for a i and b i . Solving these constraining equations, we obtain one solution Therefore, the solutions of system () are where ξ = xt and C is a constant of integration (see Figure ).

Explicit solutions of system (1) using generalized tanh method
The generalized tanh method was proposed by Fan [], who constructed a series of traveling wave solutions for some special types of equations [-]. We take the solution of reduced equations (B) of case () of Table  in the form where ϕ = ϕ(ξ ) is a solution of the Riccati equation The Riccati equation has a series of solutions () The solutions of (B) are of the form () Substituting system () along with the equation () into (B) and collecting the coefficient of the same power ϕ i , we can obtain a set of constraining equations for a i and b i . Solving these constraining equations, we obtain one solution Therefore we obtain five solutions, namely where ξ = xt (see Figure  and Figure ).

Explicit power series solutions of system (1)
We shall seek a solution to reduced equations (A) of case  in Table   Therefore, we get the following relationship between the coefficients r i and s i : The exact power series solution to system (A) can be written as follows: () Thus, we obtain the exact power series solution of system () as follows: Remark  The power series method is a useful approach to solve higher order variable coefficient ordinary differential equations. A large number of solutions for ordinary differ-

Construction of conservation laws using Lie-Bäcklund symmetries
In this section, we briefly present the notations and theorem which are useful for constructing conservation laws with the aid of Lie-Bäcklund symmetries. For more detailed information, the reader is referred to the literature [, -].
Consider a system of PDEs where u (k) are the kth order derivatives, x = (x  , x  , . . . , x n ), and u = (u  , u  , . . . , u m ). u α i denotes the derivatives of u α with respect to x i , which is explicitly expressed by u α D i is the total derivative operator. A Lie-Bäcklund operator is given by where ζ α i  i  ···i j are determined by the following relations: The Lie-Bäcklund operators () which have characteristic functions are where W α = η αξ j u α j , α = , , . . . , m are the Lie characteristic functions. The Euler-Lagrange operator is represented by A Noether operator associated with a Lie-Bäcklund operator X is defined as For example, the first order derivative of the Euler-Lagrange operator is given by The adjoint equations for the system of PDEs () are defined by The relationship between Lie-Bäcklund symmetries and conservation laws can be established by the following theorem.

Theorem  Every Lie-Bäcklund symmetry generator [, ]
of PDEs () can give rise to a conservation law for the system consisting of () and the adjoint system (). The components T i of the conserved vector T = (T  , . . . , T n ) are determined by the formula where the formal Lagrangian is Remark  For the formal Lagrangian L which has a third order derivative, the conservation laws equation () can be written (   )

Conservation laws of the DSWS
The Lie-Bäcklund transformation group [, ] can be regarded as a tangent transformation group. It is the extension of one-parameter continuous symmetry group transformations. Lie-Bäcklund symmetry group of () will be generated by the vector field of the form where U satisfies U () i | i = , i = , . By applying the third prolongation U () to system (), we obtain an overdetermined system of PDEs. The general solutions of the overdetermined system are where ϑ  , ϑ  , ϑ  are arbitrary constants. Thus, the third order Lie-Bäcklund symmetries of () are given by For system (), the formal Lagrangian is The adjoint equations of system () are given by where φ(x, t) and τ (x, t) are the new dependent variables. On the basis of equation (), one can derive components of the conservation vector for system () as follows: where W  = η u and W  = η v .
(I) For the Lie-Bäcklund symmetry generator U  = u x ∂ ∂u + v x ∂ ∂v , the Lie characteristic functions are W  = u x and W  = v x . Using equation (), we obtain the following components of the conserved vector: (II) The Lie-Bäcklund symmetry generator U  = vv x ∂ ∂u + (bu x v + kuv xav xxx ) ∂ ∂v has the Lie characteristic functions W  = vv x and W  = bu x v + kuv xav xxx . Using equation (), we obtain the following components of the conserved vector: aτ (bu xxx v + bu xx v x + bu x v xx + ku xx v x + ku x v xx + kuv xxxav xxxxx ). () (III) Finally, for the Lie-Bäcklund symmetry generator the Lie characteristic functions are W  = tvv xu xx -u and W  = btu x v + ktuv x -atv xxxv xx -v. Using equation (), we obtain the following components of the conserved vector: + (btu x v + ktuv x -atv xxxv xx -v)(kτ u + φvaτ xx ) + aτ x (btu xx v + btu x v x + ktu x v x + ktuv xx -atv xxxxv xxx -v x ) aτ (btu xxx v + btu xx v x + btu x v xx + ktu xx v x + ktu x v xx + ktuv xxx -atv xxxxxv xxxx -v xx ). ()

Conclusions
Lie symmetry analysis has been employed to investigate Lie point symmetries of the Drinfeld-Sokolov-Wilson system. The symmetries X  , X  , and X  form a threedimensional Lie algebra L  . By using Ibragimov's method, we have derived an optimal system of one-dimensional subalgebra. It is proved that the optimal system has five operators. Based on the optimal system, we have considered the symmetry reductions and group invariant solutions of the DSWS. To the best of our knowledge, very little work has been devoted to constructing conservation laws of NLPDEs by using Lie-Bäcklund symmetries. Lie-Bäcklund symmetries of the DSWS have been derived. The method of constructing conservation laws of NLPDEs with the aid of a new conservation theorem associated with Lie-Bäcklund symmetries has been presented. Conservation laws of the DSWS have been constructed by using this method.