Conservation laws and symmetry analysis of a generalized Drinfeld-Sokolov system

: The generalized Drinfeld-Sokolov system is a widely-used model that describes wave phenomena in various contexts. Many properties of this system, such as Hamiltonian formulations and integrability, have been extensively studied and exact solutions have been derived for speciﬁc cases. In this paper we applied the direct method of multipliers to obtain all low-order local conservation laws of the system. These conservation laws correspond to physical quantities that remain constant over time, such as energy and momentum, and we provided a physical interpretation for each of them. Additionally, we investigated the Lie point symmetries and ﬁrst-order symmetries of the system. Through the point symmetries and constructing the optimal systems of one-dimensional subalgebras, we were able to reduce the system of partial di ﬀ erential equations to ordinary di ﬀ erential systems and obtain new solutions for the system.


Introduction
One of the most important, and surely one of the most studied, nonlinear evolution partial differential equations (PDEs) is the Korteweg-de Vries equation (KdV) where u has the physical meaning of the wave amplitude and also the wave speed.The KdV equation arises to model shallow water waves with weak nonlinearities and it has applications in a wide variety of physical phenomena, such as surface gravity waves [1], propagation of nonlinear acoustic waves in liquids with small volume concentrations of gas bubbles [2,3], magma flow and conduit waves [4] and many other applications [5].
Drinfeld and Sokolov [6,7] obtained all the simplest generalized KdV equations corresponding to the classical Kac-Moody algebras of ranks one and two.Among these generalizations is the system of equations given as follows: This system from Eq (1. 2) is an evolutionary system for two dependent variables that possesses a Lax pair and, after linear changes and scale transformations, yields a Hamiltonian form with the Hamiltonian given by H = u 2 + v 2 − 6uv. (1.3) In the present paper, we have considered a generalization of (1.2).It is called the generalized Drinfeld-Sokolov (gDS) system where α 1 , β 1 , α 2 , β 2 , γ 0 and δ −1 are arbitrary constants.There are several papers in which different kinds of generalizations of system (1.2) have been studied.However, in many of these papers, the gDS system is named in different ways.We summarized them, giving the background and the results obtained for some of these systems.
Nevertheless, system (1.5) is known as the Satsuma-Hirota system [9].In this paper, system (1.5) was transformed into the scalar sixth-order Drinfeld-Sokolov-Satsuma-Hirota equation by using a Miura-type transformation.Moreover, three different methods were applied to obtain solutions to the equation.
The same pair of equations from (1.5) were called Bogoyavlenskii coupled KdV equations in [10][11][12].Soliton solutions were found by means of the standard and non-standard truncation Painlevé expansion and a family of exact solutions was obtained by a Bäcklund transformation.
In [13], system (1.5) was considered but named as a Drinfeld-Sokolov-Satsuma-Hirota.Exact solutions were obtained using the (G'/G)-expansion method and the generalized tanh-coth method.
System (1.5) was also studied in [14].Their non-local symmetries were studied and exact solutions like solitons, cnoidal waves and Painlevé waves were derived through similarity reductions.
Through the Weiss method of truncated singular expansions, authors of [15] constructed an explicit Bäcklund transformation of the Drinfeld-Sokolov-Satsuma-Hirota system (1.5) into itself.They also obtained special solutions generated by this transformation.
Moreover, in [18] traveling wave solutions for system (1.6) were studied.In particular, the authors obtained solutions with physical interest like kink and anti-kink wave solutions and periodic wave solutions.
Therefore, the gDS system presented in this paper embraces all the previous equations and unifies all the cases under the study of the gDS system (1.4).Note that the named generalized Drinfeld-Sokolov (1.4) does not have an infinite number of conservation laws, which will be shown in Section 2, so this system is not completely integrable.
On the other hand, in order to simplify the number of arbitrary constants involved in Eq (1.4), we have applied a scaling equivalence transformation.It is a change of variables that acts on independent and dependent variables t = λ 1 t, (1.7) From (1.7), by taking β = β 1 β 2 and α = α 1 α 2 , system (1.4) can be invertibly reduced into a system involving four parameters instead of six given by where α, β 0, γ 2 = 1 and δ −1.Hereafter, without loss of generality, we can restrict the study of system (1.4) to system (1.8).Moreover, to avoid cases already studied, we consider the extra condition δ 0. This paper has two main objectives.First, it seeks to investigate the conservation laws of the gDS system (1.8) in order to determine physical properties that remain constant over time.These conservation laws are crucial in demonstrating whether quantities such as charge, momentum, angular momentum, and energy are conserved.Second, the paper aims to study the Lie point symmetries and first-order generalized symmetries of system (1.8), which can be used to obtain reductions and exact group-invariant solutions.
To achieve the first objective, we used the direct method of multipliers [19][20][21].This method updates the classical method in which Noether's theorem was applied to the variational symmetries of the equation [22].For the second objective, we determined all symmetries by using the Lie method [21,23,24], a widely used technique for deriving analytical solutions of nonlinear partial differential equations [25,26].
In a previous paper [27], Lie point symmetries were obtained for system (1.4) with different parameters and conditions than those of system (1.8).However, the present paper extends this work by updating the Lie point symmetry classification for the gDS system, studying the first-order generalized symmetries, constructing the optimal systems of one-dimensional subalgebras, and presenting all the reduced ordinary differential systems.
The paper is organized as follows: Section 2 presents the classification of low-order local conservation laws and interprets the associated conserved quantities.In Section 3, we update the Lie point symmetry classification of the gDS system (1.8) and study the first-order generalized symmetries.In Section 4, we construct the optimal systems of one-dimensional subalgebras using the Lie point symmetries.Finally, in Section 5, we show all the reduced ordinary differential systems as an application of the Lie point symmetries.

Conservation laws
For the generalized Drinfeld-Sokolov system (1.8), a local conservation law is a divergence expression of the form that holds for the whole set of solutions (u(t, x), v(t, x)), where D t and D x denote the total derivative operators with respect to t and x, respectively.T and X are functions depending on t, x, u, v and a finite number of derivatives of u and v. T is called conserved density, and X is called spatial flux.The expression (T, X) is the conserved current.
A conservation law is called locally trivial if there is a function where Φ = Ψ = 0 holds on all solutions (u(t, x), v(t, x)).Hence, conservation laws are considered to be locally equivalent if they differ by a locally trivial conservation law.Notice that the system is formed by a couple of evolution equations, so through use of the system (1.8) and its differential consequences, it is possible to eliminate t-derivatives of u and v in (T, X).
Every non-trivial local conservation law (2.1) for Eq (1.8) can be expressed in its characteristic form where the pair of functions is called a multiplier and δ δu and δ δv are the Euler operators with respect to the variables u and v, respectively [21,23,28].
Since the Euler operator has the property of annihilating total derivatives, all multipliers (Q u , Q v ) can be found by solving the system of determining equations that is given by holding for the set of solutions of Eq (1.8).
All non-trivial conservation laws arise from multipliers (2.3), and conservation laws of physical importance come from low-order multipliers [29][30][31] of the form (2.5) The term low-order is related to the dependencies of the multiplier function (2.5).In general, a loworder multiplier depends on the dependent variables u, v and those variables that can be differentiated up to the leading derivatives of the system.The leading derivatives of system (1.8)Similar to the determining equations for Lie symmetries, here the determining Eq (2.4) for loworder multipliers split with respect to u t , v t , u xxx , v xxx and differential consequences.
We have obtained a general classification for low-order multipliers (2.5), depending on the parameters that system (1.8) involves, by solving the determining Eq (2.4).For this purpose we have used Maple, "rifsimp" and "pdsolve" commands.
Proposition 2.1.The general classification of low-order multipliers (2.5) for the generalized Drinfeld-Sokolov system (1.8) is: (ii) There are additional multipliers for the following special α, β, δ: (2.9) Proof.The multipliers determining system (2.4) leads to an overdetermined linear system of 57 equations.By solving the easiest equations, we obtain that related by the following conditions: (2.11) The solutions of system (2.11) lead to the different cases listed in Proposition 2.1.
In the general case, where α, β 0, γ 2 = 1 and δ 0, −1, the admitted conservation law (2.12)yields the conserved quantity that is the mass density.Conservation law (2.13), in the case δ = 1, leads to the conserved quantity which is the momentum density.
In the case where δ = 1 and 2α = β, the conservation law (2.14)yields the conserved quantity which represents the total energy for solutions u(t, x), v(t, x).
For δ = 1 and α = −1, the conserved quantity associated to conservation law (2.15) is which is the Galilean momentum.
In the case where α = 1 2 and β = 1, the two admitted conservation laws (2.16) and (2.17) yield the conserved quantities ) The first one represents the angular momentum, and the second one the energy-momentum.

Lie point symmetries
Now, we consider a one-parameter Lie group of point transformations acting on independent and dependent variables t = t + ετ(t, x, u, v) where ε is the group parameter.The associated vector field has the form The Lie group of transformations (3.1) are the point symmetries admitted by the system if and only if it leaves invariant the solution space of system (1.8) .It is known as the infinitesimal criterion of invariance: when the gDS system (1.8) holds, where X (3) is the third-order prolongation of the vector field (3.2).Equation (3.3) is known as the determining equation.These equations can be split by considering u, v and its derivatives as independent variables.Consequently, it leads to an overdetermined linear system of equations for the infinitesimals τ, ξ, η u , η v together with the arbitrary parameters, α, β, γ, δ subject to the classification conditions α, β 0, γ 2 = 1, δ 0, −1.We obtain the following result.
Theorem 3.1.The classification of point symmetries admitted by the generalized Drinfeld-Sokolov system (1.8) is given by: (i) For the non-zero arbitrary parameters α, β, γ, δ, there are three infinitesimal generators (ii) For α = 1, an additional point symmetry is admitted Proof.The symmetry determining system (3.3)leads to an overdetermined linear system of 33 equations.By solving the easiest equations for the infinitesimals, we obtain that τ = τ(t), ξ = ξ(t, x), From the first two equations of system (3.6), we obtain that B(t, x, u) = E(t, x)u+F(t, x) and D(t, x, v) = G(t, x)v + H(t, x).Substituting the functions in the previous system, we get Now the functions involved in system (3.7)do not depend on u or v, so these equations split with respect to u, v and their derivatives.We derive and solve this overdetermined system with the help of Maple software.Using the commands "rifsimp" and "pdsolve", we split the system and solve the equations with results leading to the different cases listed in Theorem 3.1.
We interpret the symmetries by looking at the forms of the transformations.In case 1, transformation (3.8a) represents a time-translation, transformation (3.8b) represents a space-translation and transformation (3.8c) represents a scaling.For case 2, the additional symmetry (3.8d) is a Galilean boost.

First-order generalized symmetries
A first-order symmetry corresponds to a one-parameter Lie group of transformations in which the transformation of (t, x, u, v) necessarily depends on first-order derivatives of u and v. Hence, its transformations are given by where ε is the group parameter.Thus, the infinitesimal generator takes the form (3.10) Following a similar procedure to the one shown before for the point symmetries, applying the criterion of invariance and considering that ζ t , ζ x , σ t , σ x do not cancel at the same time, we have obtained the following result.
Proof.By splitting and solving the system of determining equations it is straightforward to prove that the solutions do not depend on first-order derivatives of the dependent variables.Therefore, all the first-order symmetries reduce to prolonged point symmetries of system (1.8).

Optimal systems
The construction of the optimal system is an efficient way of classifying all possible group-invariant solutions.From these group-invariant solutions, every other solution can be derived.
For one-dimensional subalgebras, finding an optimal system is the same as classifying the orbits of the adjoint representation.Essentially this problem is attacked by taking a general element X in the symmetry algebra and subjecting it to various adjoint transformations so as to simplify it as much as possible.
Let us consider a general element from the Lie algebra and reduce it to its simplest equivalent form by using the chosen adjoint transformations where ε is a real number, and X i , X j denotes the commutator operator defined by The commutator and adjoint tables of the generalized Drinfeld-Sokolov system (1.8) are given in Tables 1 and 2, corresponding to the general case and the α = 1 case, respectively.By using these tables, we simplify as much as possible the maximal Lie algebra (made up from a linear combination of all its infinitesimal generators) through the adjoint table in order to get the optimal systems.It is important to point out paper [34], where all three and four-dimensional indecomposable real Lie algebras are classified.
Theorem 4.1.For the generalized Drinfeld-Sokolov system (1.8), an optimal system of onedimensional subalgebras in each case is given by: (i) Case α, β, γ, δ non-zero arbitrary constants with α 1.An optimal system of one-dimensional subalgebras is given by {X (ii) Case α = 1 and β, γ, δ non-zero arbitrary constants.An optimal system of one-dimensional subalgebras is given by {X Proof.Consider a general element in the symmetry algebra The goal is to simplify as many of the a i coefficients as possible through the application of adjoint maps to X. Without loss of generality, suppose first that a 1 0. We can assume that a 1 = 1.Referring to Table 2, if we act on X by X 1 we can make the coefficient of X 2 vanish, obtaining a new X = X 1 +a 3 X 3 +a 4 X 4 .We continue acting on X to cancel the remaining coefficients until no further simplifications are possible.Following a similar procedure for the rest of cases, we get the results of Theorem 4.1.

Symmetry reductions
We reduce the generalized Drinfeld-Sokolov system (1.8) of PDEs to a system of two ordinary differential equations (ODEs) by using the symmetries shown in Theorem 3.1 and the optimal systems of one-dimensional subalgebras in Theorem 4.1.
These reductions are obtained by doing a change of variables with the solutions of the characteristic system dt τ From (5.1), we obtain the similarity variable z and the similarity solutions U(z), V(z), which are substituted into system (1.8) to get the reduced one.
5.1.Reduced system associated to X 1 + c 1 X 2 Considering α, β non-zero arbitrary parameters and by taking into account the generator X 1 + c 1 X 2 , we reduce the gDS system (1.8) into a nonlinear ODE system.
(5.8) 5.3.Reduced system associated to X 1 + c 2 X 4 In the case of α = 1, β a non-zero arbitrary parameter, and considering the generator X 1 + c 2 X 4 , we obtain the similarity variable and solutions (5.9) By substituting (5.9) into system (1.8), we reduce it to the third-order ODE system: (5.10)
are u t , u xxx , v t , and v xxx .Hence, u t , v t can be obtained by differentiation of u, v respectively, u xxx by differentiation of u, u x , u xx and v xxx by differentiation of v, v x , and v xx .