Lie Algebra Classification, Conservation Laws, and Invariant Solutions for a Generalization of the Levinson–Smith Equation

We obtain the optimal system’s generating operators associated with a generalized Levinson–Smith equation; this one is related to the Liénard equation which is important for physical, mathematical, and engineering points of view.)e underlying equation has applications in mechanics and nonlinear dynamics as well. )is equation has been widely studied in the qualitative scheme. Here, we treat the equation by using the Lie group method, and we obtain certain operators; using those operators, we characterized all invariants solutions associated with the generalized equation of Levinson Smith considered in this paper. Finally, we classify the Lie algebra associated with the given equation.


Introduction
Lie group symmetry method is a powerful tool employed to study ODEs, PDEs, FPDEs, FODEs, and so on. is theory was introduced in the 19th century by Sophus Lie [1], following the idea of Galois theory in algebra. Lie group method applied to differential equations has received great interest among researchers in different fields of science such as mathematics, theoretical, and applied physics, due to the physical interpretations of the underlying equations that are studied. As a consequence, this method leads to construct, for example, conservation laws, using the well known Noether's theorem [2], even more applying Ibragimov's approach [3]. In the same way, it is possible to build similarity solutions which, in the traditional methods, are not possible.
Furthermore, this method contributes to establish schemes and the usefulness of some numerical methods; here, many packages are being built in different environments of computations, e.g., [4,5]. In general, taking into account the importance of the equations' study (such as ODEs, PDEs, and others), this method can be interesting to different researchers. A vast reference in Lie group method can be found in the literature, e.g., [6][7][8][9]. Recently, the Lie group method approach has been applied to solve and analyze different problems in many scientific fields, e.g., in [10], the authors applied the Lie symmetry method to investigate a fourth-order 1 + 2 evolutionary partial differential equation which has been proposed for the image processing noise reduction. References in the latest progress in symmetry analysis can be found in [11][12][13][14][15][16][17][18] and therein.
In [19], Kamke proposes the following differential equation: where f, g, h are the arbitrary functions, for this equation presents the transformation y(x) � η(ξ(x)) which reduces this equation to a system of two first-order equations. Equation (1) can be written as y xx + ϕ x, y, y x y x � c(x, y), where ϕ(x, y, y x ) � − ((f y (y)/f(y))y x + g(x)) and c(x, y) � − h(x)f(y).
is means that the coefficient of friction, i.e., ϕ is a function that depends on x, y, and y x , and it will almost always be a nonlinear function, and on the other hand, function c, which is known as function of disturbance, is also nonlinear. It is worth remembering that this type of equation (2) is known as generalizations of the Levinson-Smith equation. Also, equation (2), which is a particular case of the generalized Levinson-Simth equations, is related to the Liénard-type second-order nonlinear differential equation [20]. e underlying equations describe several phenomena in different areas such as electronics, biology, mechanics, seismology, chemistry, physics, and cosmology, for example, an important model in physical and biological sciences is the van der Pol equation, which describes a nonconservative oscillator with nonlinear damping. Levinson and Smith in [21] studied a general equation for relaxation oscillations.
In [8], the equation, is presented; note that equation (3) is a particular case of equation [8], Cantwell states that equation (3) has a group of symmetries 5− dimensional; but it does not exhibit the development of said statement; they affirm the Lie group of symmetries of (3), using a ODEtools Maple package. In fact, the goal of this work is (i) to calculate the 5− dimensional Lie symmetry group in detail, (ii) to present the optimal system (optimal algebra) for (3), (iii) to make use of all elements of the optimal algebra to propose invariant solutions for (3), then (iv) to construct the Lagrangian with which we could determine the variational symmetries using Noether's theorem and thus to present conservation laws associated, and (iv) also to build some nontrivial conservation laws using Ibragimov's method, and finally (v) to classify the Lie algebra associated to (3), corresponding to the symmetry group.

Continuous Group of Lie Symmetries
In this section, we study the Lie symmetry group for (3). e main result of this section can be presented as follows.

Proposition 1.
e Lie symmetry group for equation (3) is generated by the following vector fields: Proof. A general form of the one-parameter Lie group admitted by (3) is given by where ϵ is the group parameter. e vector field associated with the group of transformations shown above can be written as Γ � ξ(x, y)(z/zx) + η(x, y)(z/zy), where ξ, η are differentiable functions in R 2 . Applying its second prolongation, to equation (3), we must find the infinitesimals ξ, η, satisfying the symmetry condition, associated with (3). Here, η [x] , η [xx] are the coefficients in Γ (2) given by Being D x as the total derivative operator, D x � z x + y x z y + y xx z y x + · · ·. Replacing (8) into (7) and using (3), we obtain 2 International Journal of Differential Equations From (9), canceling the coefficients of the monomial variables in derivatives 1, y 3 x , y 2 x , and y x , we obtain the determining equations for the symmetry group of (3). at is, (9d) Solving the system of equations (9a)-(9d) for ξ and η, we get η � c 1 y + c 2 y 2 + c 3 xy 2 + c 4 x 2 y 2 + c 5 − 2x 2 y 2 + x 3 y 2 .

Optimal System
Taking into account [22][23][24][25], we present in this section the optimal system associated to the symmetry group of (3), which shows a systematic way to classify the invariant solutions. To obtain the optimal system, we should first calculate the corresponding commutator table, which can be obtained from the operator where i � 1, 2, with α, β � 1, . . . , 5, and ξ i α , ξ i β are the corresponding coefficients of the infinitesimal operators Π α , Π β . After applying the operator (11) to the symmetry group of (3), we obtain the operators that are shown in Table 1. Now, the next thing is to calculate the adjoint action representation of the symmetries of (3), and to do that, we use Table 1 Making use of this operator, we can construct Table 2, which shows the adjoint representation for each Π i .

Proposition 2.
e optimal system associated to equation (3) is given by the vector fields

International Journal of Differential Equations
Proof. To calculate the optimal system, we start with the generators of symmetries (4) and a generic nonzero vector. Let e objective is to simplify as many coefficients, a i , as possible, through maps adjoint to G, using Table 2. (1) Assuming a 5 � 1 in (14), we have that Applying the adjoint operator to (Π 1 , G) and (Π 3 , G), we do not have any reduction; on the other hand, applying the adjoint operator to (Π 2 , G), we get (1.1.A.A 1 ) Case a 1 + a 3 ≠ 0. Using λ 3 � (1/2(a 1 + a 3 ))), with a 1 + a 3 ≠ 0, in (16), Π 5 is eliminated, therefore G 3 � a 1 Π 1 + a 3 Π 3 . en, we have the first element of the optimal system.
en, we have the other element of the optimal system.
with a 3 ≠ 0. is is how the other reduction of the generic element (14) ends.
Applying the adjoint operator to (Π 5 , G 2 ), we have en, we have the other element of the optimal system.
en, we have the other element of the optimal system.
is is how the other reduction of the generic element (14) ends.
It is clear that we do not have any reduction.
Applying the adjoint operator to (Π 5 , G 16 ), we have Using 4 . en, we have the other element of the optimal system.
It is clear that we do not have any reduction.
en, we have the other element of the optimal system.
(2) Assuming a 5 � 0 and a 4 � 1 in (14), we have that G � a 1 Π 1 + a 2 Π 2 + a 3 Π 3 + Π 4 . Applying the adjoint operator to (Π 1 , G) and (Π 3 , G), we do not have any reduction; on the other hand, applying the adjoint operator to (Π 2 , G), we get Applying the adjoint operator to (Π 5 , G 5 ), we get It is clear that we do not have any reduction.
(2.1.A.A 2 ) Case a 1 + a 3 � 0. We get G 6 � − a 3 Π 1 + a 3 Π 3 , and then we have the other element of the optimal system, with a 1 , a 3 ≠ 0 and a 1 + 2a is is how the other reduction of the generic element (14) ends.
It is clear that we do not have any reduction; it is also clear that a 1 ≠ 0 and then a 3 ≠ 0; then, substituting λ 23 � b 13 , we have the other element of the optimal system with a 3 ≠ 0 y a 1 + 2a 3 � 0. is is how the other reduction of the generic element (14) ends.
(2.2.A.1) Case a 2 ≠ 0. en, substituting λ 21 � (+b 11 /2a 2 ) with a 2 ≠ 0, we have other element of the optimal system with a 2 , a 3 ≠ 0. is is how the other reduction of the generic element (14) is is how the other reduction of the generic element (14) ends.
(2.2.B.2) Case a 2 � 0. We get G 22 � Π 4 , and then we have the other element of the optimal system is is how the other reduction of the generic element (14) ends.
(3) Following a procedure analogous to the previous one and analyzing the respective cases for a 4 � a 5 � 0, a 3 � 1 in (14); a 3 � a 4 � a 5 � 0, a 2 � 1 in (14) and a 2 � a 3 � a 4 � a 5 � 0, a 1 � 1 in (14); we can reduce and obtain all the elements presented for the optimal system. □

Invariant Solutions by Some Generators of the Optimal System
In this section, we characterize invariant solutions taking into account all operators that generate the optimal system presented in Proposition 2. For this purpose, we use the method of invariant curve condition [23] (presented in Section 4.3), which is given by the following equation: Using the element Π 5 from Proposition 2, under the condition (39), we obtain that Q � η 5 − y x ξ 5 � 0, which implies (− 2x 2 y + x 3 y 2 ) − y x (x 3 ) � 0; then, solving this ODE, we have y(x) � (1/x(c 1 x + 1)), where c is an arbitrary constant, which is an invariant solution for (3); using an analogous procedure with all of the elements of the optimal system (Proposition 2), we obtain both implicit and explicit invariant solutions that are shown in Table 3, with c being a constant.

Variational Symmetries and Conserved Quantities
In this section, we present the variational symmetries of (3) and we are going to use them to define conservation laws via Noether's theorem [26]. First of all, we are going to determine the Lagrangian using the Jacobi last multiplier method, presented by Nucci in [27], and for this reason, we are urged to calculate the inverse of the determinant Δ, Δ � x y x y xx where Π 2,x , Π 2,y , Π 4,x , and Π 4,y are the components of the symmetries Π 2 , Π 4 shown in the Proposition 4 and Π (1) 2 , Π (1) 4 are its first prolongations. en, we get Δ � 2xy 4 which implies that M � (1/Δ) � (x − 1 y − 4 /2). Now, from [27], we know that M can also be written as M � L y x y x which means that L y x y x � (x − 1 y − 4 /2); then, integrating twice with respect to y x , we obtain the Lagrangian, where f 1 , f 2 are arbitrary functions. From the preceding expression, we can consider f 1 � f 2 � 0. It is possible to find more Lagrangians for (3) by considering other vector fields given in the Proposition 4. We then calculate using (41) and (8). us, we get From the preceding expression, rearranging and associating terms with respect to 1, y x , y 2 x and y 3 x , we obtain the following determinant equations: Solving the preceding system for ξ, η and f, we obtain the infinitesimal generators of Noether's symmetries η � − 2a 1 yx 2 − a 2 y 2 x 2 + a 4 y 2 + a 5 y, with a 1 , a 2 , a 3 , a 4 , and a 5 arbitrary constants. en, the Noether symmetry group or variational symmetries are International Journal of Differential Equations
Definition 1. Let F be a differential function and ] � ](x) be the new dependent variable, known as the adjoint variable or nonlocal variable [31]. e formal Lagrangian for F � 0 is the differential function defined by Definition 2. Let F be a differential function and for the differential equation (49), denoted by F[y] � 0, we define the adjoint differential function to F by and the adjoint differential equation by where the Euler operator and D x i is the total derivative operator with respect to x i defined by Definition 3. e differential equation (40) is said to be nonlinearly self-adjoint if there exists a substitution such that for some undetermined coefficients λ � λ(x, y, . . .). If ] � ϕ(y) in the two previous expressions, equation (49) is called quasi-self-adjoint. If ] � y, we say that equation (49) is strictly self-adjoint. Now, we shall obtain the adjoint equation to equation (3). For this purpose, we write (3) in the form (49), where en, the corresponding formal Lagrangian (50) is given by and the Euler operator (53) assumes the following form: We calculate explicitly the Euler operator previously applied to L determined by (58). In this way, we obtain the adjoint equation (52) to (3): (60) e main result in this section can be stated as follows.

Proposition 3. Equation (3) is nonlinearly self-adjoint, with the substitution given by
where k 1 , k 2 are the arbitrary constants.

International Journal of Differential Equations
Proof. Substituting in (60), and then in (56), ] � ϕ(x, y) and its respective derivatives, and comparing the corresponding coefficients, we get the following five equations: We observe that equation (62c) is obtained from equation (62b) by differentiation with respect to x. erefore, we have to study only equations (62b), (62d), and (62e). Solving for ϕ in (62b), we obtain where c 1 (x) is the arbitrary function. Using (63) in (62e), we get − (c 1 (x)/x 2 ) + c 1,xx (x) � 0, thus solving for c 1 (x), we have then, substituting in (63), the statement in the theorem is obtained.

Conservation Laws
In this section, we shall establish some conservation laws for equation (3) using the conservation theorem of N [31]. Since equation (3) is of second order, the formal Lagrangian contains derivatives up to order two. us, the general formula in [31] for the component of the conserved vector is reduced to where where j � 1, . . . , 5 is the formal Lagrangian (58), and η j , ξ j are the infinitesimals of a Lie point symmetry admitted by equation (3), given in (4). Using (3), (4), and (3) in (64), we obtain the following conservation vectors for each symmetry stated in (4). where

Classification of Lie Algebra
Generically, a finite-dimensional Lie algebra in a field of characteristic 0 is classified by the Levi's theorem, which states that any finite-dimensional Lie algebra can be written as a semidirect product of a semisimple Lie algebra and a solvable Lie algebra; the solvable Lie algebra is the radical of that algebra. In other words, there exist two important classes of Lie algebras, the solvable and the semisimple. In each class mentioned above, there are some particular classes that have other classifications, for example, in the solvable one, we have the nilpotent Lie algebra.
According the Lie group symmetry of generators given in Table 1, we have a five-dimensional Lie algebra. First of all, we remember some basic criteria to classify a Lie algebra. In the case of solvable and semisimple Lie algebra, we will denote K(·, ·) as the Cartan-Killing form. e next propositions can be found in [32]. Proposition 6. Let g 1 and g 2 be two Lie algebras of dimension n. Suppose each has a basis with respect to which the structure constants are the same. en, g 1 and g 2 are isomorphic.
Let g be the Lie algebra related to the symmetry group of infinitesimal generators of equation (1) as stated by the table of the commutators; it is enough to consider the next relations: Using that, we calculate Cartan-Killing form K as follows: where the determinant vanishes, and thus by Cartan criterion, it is not semisimple (see Proposition 4). Since a nilpotent Lie algebra has a Cartan-Killing form, that is, identically zero, we conclude, using the counter-reciprocal of the last claim, that the Lie algebra g is not nilpotent. We verify that the Lie algebra is solvable using the Cartan criteria to solvability (Proposition 5), and then we have a solvable non-nilpotent Lie algebra. e nilradical of the Lie algebra g, M is generated by Π 2 , Π 4 , Π 5 , and it is isomorphic to h 3 , the Heisenberg Lie algebra, and so we have a solvable Lie algebra with three-dimensional nilradical. Let m be the dimension of the nilradical M of a solvable Lie algebra. In this case, in fifth-dimensional Lie algebra, we have that 3 ≤ m ≤ 5.
Mubarakzyanov in [33] classified the 5-dimensional solvable non-nilpotent Lie algebras, in particular the solvable non-nilpotent Lie algebra with three-dimensional nilradical. en, by Proposition 6 and consequently, we establish an isomorphism of Lie algebras with g and the Lie algebra g 5,35 . In summary, we have the next proposition.

Proposition 7.
e 5-dimensional Lie algebra g related to the symmetry group of equation (1) is a solvable non-nilpotent Lie algebra with three-dimensional nilradical; this nilradical is isomorphic to h 3 , the Heisenberg Lie algebra. Besides that, Lie algebra is isomorphic with g 5,35 in the Mubarakzyanov's classification.

Conclusion
For a generalized Levinson-Smith equation (3), we obtained the optimal system's generating operators (see Proposition 2); using those operators, it was possible to characterize all invariant solutions as it is shown in Table 3; these invariant solutions do not appear in the literature known until today.
It has been shown in the variational symmetries for (3), as it was shown in (46) with its corresponding conservation laws (48) and all these were using Noether's theorem, but nontrivial conservation laws were also calculated using the Ibragimov's method as it is shown in (67) using the nonlinearly self-adjoint of equation (3) as mentioned in Proposition 3. e results obtained in this study are new, and according to the phenomena that govern this equation, which reaches several fields of science, for instance, the nonlinear oscillators, it may be of significant importance for several researchers. erefore, the goal initially proposed was achieved.
e Lie algebra associated to equation (3) is a solvable non-nilpotent Lie algebra with three-dimensional nilradical, and it is isomorphic with g 5,35 in the Mubarakzyanov's classification; therefore, the goal initially proposed was achieved.
For future works, equivalence group theory could be also considered to obtain preliminary classifications associated to a complete classification of (3).

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.