Abstract

Differential invariants and their corresponding canonical forms for systems of three 2nd-order ODEs possessing three-dimensional Lie algebras are constructed. Their extension up to th-order system of three 2nd-order ODEs is presented. Furthermore singularity in invariant structure for the canonical forms is investigated. In addition integrability of these canonical forms is discussed. Illustrative physical examples from mechanics of system of particles are provided.

1. Introduction

There are several physical phenomena whose mathematical modeling is associated with system of 2nd-order ODEs and hence analysis of various aspects of these systems of ordinary differential equations plays a vital role in the applied sciences. Due to this contributing role and importance, these systems have been extensively studied over the years. Different methods and approaches have been introduced and constructed for various aspects of the analysis of these systems.

One of the prominent approaches is Lie symmetry method, initiated by Sophus Lie 1842–1899. This is a general method in which exact solutions and various other aspects of linear and nonlinear differential equations are analysed. Lie theory is based on Lie groups and the corresponding Lie algebras. Lie proved that the integrability of a differential equation is contingent on the properties of the Lie algebra it admits; for details see [15].

In Lie theory, there are two approaches for the integrability of differential equations: direct and indirect approach. In the direct approach classical reduction is performed via differential invariants and canonical variables and so forth, while in the indirect approach equations are classified via their admitted Lie algebras in prescribed number of underlying variables with corresponding set of structure constants. Then algebraic realizations are used for integration of these classified forms. The indirect approach is also known as Lie algebraic approach. The direct approach is useful for scalar ordinary differential equations only but does not work for system of ordinary differential equations [6]. The indirect approach is equally valid for both scalar ODEs as well as system of ordinary differential equations.

In this paper, we are interested in Lie algebraic approach (Realization approach) for system of 2nd-order ordinary differential equations.

To apply the algebraic approach, equations have to be classified according to their admitted symmetry algebras. Moreover this classification scheme depends not only upon the abstract Lie algebras (and their subalgebras) but also the realizations of the admitted Lie algebras in prescribed number of variables.

The initial work on classification and realizations of Lie algebras was done by Lie [7]. After Lie, several researchers have been involved in analysis of different aspects of this work; see [1, 6, 812]. Although the algebraic approach is more involved, it generates new cases and more insight to mathematical properties whose base is algebra (see [8, 9, 13, 14]). Thus the classification of differential equations via algebraic approach is more efficient than other classical approaches. Unluckily, very limited literature is available on this approach for system of 2nd-order ordinary differential equations.

Like systems of two 2nd-order ODEs, systems of three 2nd-order ODEs have vital role in different mathematical modeling of physical situations such as small oscillation problems, wave propagation problems, and problems of mechanics. There is very limited and restricted work on classification of system of three 2nd-order ODEs according to their symmetry algebras (see [15, 16]). Thus in this paper, we are interested in investigating classification of system of three 2nd-order ODEs via the algebraic approach for dimension 3.

The outline of the rest of the paper is as follows. In Section 2 invariant structure of canonical forms and their Lie algebraic properties are constructed. Section 3 is about the integrability of underlying canonical forms, whereas Section 4 presents illustrative examples of the results obtained in this paper. Paper ends with a brief conclusion.

2. Differential Invariants and Canonical Forms

By fixing number of prescribed variables for a given set of structure constants associated with an underlying Lie algebra one can obtain the most general classes of system of differential equations which admit the investigated Lie algebra in prescribed number of variables.

Here group classification of system of three 2nd order ODEs, admitting three-dimensional Lie algebra, has been investigated via algebraic approach. The classification results of [11] are utilized but in the format of [8].

2.1. Invariant Construction

Following result given in [8] is fundamental in construction of invariants and their corresponding canonical forms.

Proposition 1. If and are invariants of a Lie algebra admitted by any system of ordinary differential equations, then is also its invariant.

For the construction of canonical forms the approach developed in [8, 9, 13] is employed. By utilizing realizations of underlying 3-dimensional Lie algebras admitted by system of three 2nd-order ordinary differential equations we find differential invariants associated with them. These invariants are then used to construct invariant system of differential equations, called canonical forms of system.

Notations. The following notions are used in rest of the paper: denotes the th Lie algebra of dimension whereas superscripts indicate parameters on which the algebra possibly depends; the column in Table 1 gives the algebra realizations; the realization is referred to by a superscript ; as usual . denotes the elements of a basis of a given Lie algebra; here is less than or equal to the dimension of the underlying real Lie algebra. The rank of the associated realization is denoted by and is the realization of the corresponding algebra.

The realizations given in [8] are used to obtain the differential invariants and their corresponding invariant representation for the system of three 2nd-order ordinary differential equations. There are various methods and approaches for finding differential invariants but we follow the peculiar approach developed by Ayub et al. in [8, 9, 13]. In this approach invariants are interrelated via basic invariant. We explain the algorithm by the following example.

Example 2. Consider the second realization from [8]. We find the second prolongation of generators and then adopting the procedure of construction of differential invariants described in [8, 9, 13], we obtain the following differential invariants:This gives the following invariant representation and associated canonical form for the system of three 2nd-order ODEs admitting Lie algebra : Here and in the subsequent work , , , , , , and denote invariants whereas , , and are arbitrary functions.

Example 3. Consider the algebra given in [8]. Proceeding as in above example we arrive at the following invariants:Using these invariants we obtain the following invariant representation and corresponding canonical form for the system of three 2nd-order ODEs possessing Lie algebra :

In classification problem related to system of differential equations, invariant interrelation plays a key role in different types of analysis, namely, integrability, Linearization, and singularity as well as some other aspects especially related to system of three 2nd-order ODEs (see [8, 9, 13]). By careful observations here we have obtained some relations between invariants. On the basis of these relations we categorize the regular differential invariants in the following main two types.

(i) Invariant Representation in First-Order System of ODEs. In this system, for underlying Lie algebra with peculiar realization, we obtain seven linearly independent invariants whose invariant representation is in the form of system of first-order ordinary differential equations. These invariants contain four first-order differential invariants and three second-order differential invariants. The four first-order differential invariants include one invariant that plays a role of basic invariant (). The second-order differential invariants are derived via differentiation property of first-order differential invariants with respect to the basic invariant, for example, .

(ii) Invariant Representation in Mixed Type System of 2nd-Order ODE. In this system, for underlying Lie algebra with peculiar realization, we obtain seven linearly independent invariants whose invariant representation is in the form of mixed type system of 2nd-order ordinary differential equations. These invariants contain four first-order differential invariants and three 2nd-order differential invariants. The four first-order differential invariants include one invariant that plays a role of basic invariant () and remaining three; one of the invariant is derived by the differentiation of another invariant from the remaining two with respect to the basic invariant. The 2nd-order differential invariants contain one independent invariant () and the remaining two are obtained via differentiation of first-order differential invariants with respect to the basic invariant, for example, .

Both types of invariants are presented in Table 1. These observations led us to the following.

Remark 4. The number of functionally independent invariants for the regular system of three 2nd-order ODEs admitting three-dimensional Lie algebra is .

2.2. Singularity in Invariant Structure

If the rank of solution manifold is less than the rank of generic manifold, then invariant differential equation is singular [13]. This singularity plays an important role in the integrability of such systems.

For example, consider the case of Lie algebra . By employing the same procedure as used for the invariant construction of and in the previous subsection, the following set of differential invariants is found:The corresponding invariant representation and canonical form for the system of three 2nd-order ordinary differential equations having the Lie algebra , respectively, areBy employing the same procedure as given in [13], it is deduced that the condition of singularity is satisfied for the set of invariants given in (5) and hence possesses the singularity structure.

This type of singularity in invariant structure occurs in presented in Table 1.

If we consider the Lie algebra , we obtain the following set of differential invariants:whereas and do not form the system of three 2nd-order ODEs.

2.3. Extension in System of Three kth-Order ODEs

Utilizing the interrelationship between invariants of three second-order ODEs, three th-order ordinary differential equations admitting three-dimensional algebras can be constructed, as done in [8].

For this purpose we utilize the basic invariant that is introduced in the invariant representation for the classification. Consider the case of that belongs to the invariant representation in mixed type system of 2nd-order ODEs. We have the following invariant structure for this case:Differentiating the set of second-order differential invariants with respect to basic invariant (), we find the following 3rd-order differential invariants:Thus the set of functionally independent invariants for the system of three 3rd-order ordinary differential equations possessing iswhere are as mentioned in (9). Consequently, the number of functionally independent invariants for the system of three 3rd-order ODEs admitting is . We arrive at the following invariant representation as well as corresponding canonical form for the system of three 3rd-order ODEs possessing Lie algebra :Proceeding in similar manner, we arrive at the following invariant representation for the system of three th-order ODEs possessing Lie algebra :By careful observation we can classify the invariant representation for all cases of the three-dimensional Lie algebras with corresponding realization admitted by the system of three th-order ordinary differential equations admitting three-dimensional Lie algebra.

Theorem 5. If a regular system of three th-order ordinary differential equations admits three-dimensional Lie algebra then its canonical form can be represented by a set of basic functionally independent invariants in one of the following forms:where (a) corresponds to these algebras(b) corresponds toand (c) corresponds to

Remarks 6. The classification results of (a) and (b) mentioned in Theorem 5 are also valid for .
There are 31 canonical forms for the case of a system of three 3rd- and higher-order ordinary differential equations admitting three-dimensional Lie algebras, and these can be constructed by invariant differentiation as described above.
Note that and are admitted by a regular system of two 3rd- and higher-order ordinary differential equations whereas they do not form a system of three 2nd-order ordinary differential equations.

3. Integrability

It is well known that symmetries can be used to integrate ODEs either by successive reduction or by using the canonical forms. Though the later method requires classification according to their admitted algebras before integration but this method is applicable to all kinds of systems of ODEs whereas the former applies only to scalar ODEs. Here we provide an integration algorithm for canonical forms (given in Table 1). The algorithm depends on the subdivision scheme of the canonical forms that is based on general observations.

3.1. Integration Algorithm

The canonical forms for systems of three 2nd-order ODEs presented in Table 1 are subdivided mainly in two classes: Type-I and Type-II.

Type-I: Case of Singular Invariants. This consists of two cases and .

(i) . The corresponding canonical form and invariant representation for this symmetry Lie algebra arewhere and are taken from Table 1. By using (21), we arrive atThus utilizing the values of (23) in (20) and (22), we deduce thatInverting back into original variables, (23) and (24), we get the following form:The solution of (27) gives the value of , and then by utilizing this value in (26), one can get the solution of system (17), (18), and (19).

Remark 7. The solution of (27) can be obtained by employing different symmetry approaches including Lie canonical form method.
(ii) . The associated canonical form possessing symmetry Lie algebra iswhere and are taken from Table 1. By using (28), we deduce thatThe system of three second-order ODEs admitting symmetry Lie algebra is integrable by quadrature provided that system of two second-order ODEs formed by (29) and (30) is integrable by quadrature.

On the basis of reductions discussed above, we arrive at the following results.

Proposition 8. If a singular system of three 2nd-order ordinary differential equations admits a three-dimensional solvable symmetry algebra, then the general solution of the singular system can be obtained by quadratures from the general solution of the invariant representation of the singular system given by the solution of (a)a 2nd-order differential equation in the case of ;(b)system of two 2nd-order ordinary differential equations in the case of .

Remark 9. Proposition 8 is related to algebraic approach for analysis of singular system of ODEs. Proposition 8 leads to the solution of underlying singular system of three 2nd-order ODEs by either of two ways: case (a) when corresponding reduced 2nd-order differential equation has integrable class or case (b) when associated system of two 2nd-order ordinary differential equations is integrable. In both of the above cases ((a) and (b)), for integrability one can consult the Lie canonical form method of scalar second-order ODEs and other symmetry approaches as well in the former case and in latter case one can utilize approaches discussed in [6, 8, 9, 13].

Proposition 10. If a regular system of three th-order ordinary differential equations possesses a three-dimensional solvable symmetry algebra, then the general solution of the system can be obtained by quadratures from the general solution of the invariant representation given by one of the following systems:

The proof follows from Proposition 8 and the discussion preceding Proposition 8.

Remark 11. Proposition 10 is also related to algebraic approach for analysis of regular system of ODEs. Proposition 10 leads to the solution of underlying regular systems of three th-order ODEs admitting a three-dimensional solvable symmetry algebra when corresponding canonical forms of these systems given in cases (a), (b), and (c) are integrable.

Type-II: Case of Regular Invariants. This type consists of all the cases of regular invariants mentioned in Table 1 and their integrability depends on the peculiar structure of associated system of three 2nd-order ordinary differential equations with admitting Lie algebras.

3.1.1. Examples

Let us consider two examples of physical systems which describe the integration algorithm that depends on the peculiar structure of associated system of three 2nd-order ODEs and the admitted Lie algebra.

() Consider the motion of a system of three particles defined by in a small region; such that for any small number , where There is a time varying applied force which produces the motion after neglecting the effects of gravity due to small region. System (33) represents the motion of system of interacting particles in which the first particle faces the resistance proportional to the velocity of third particle, while second particle faces the drag resistance and third particle is facing resistance equal to the product of the velocities of second and third. This system admits the symmetry Lie algebra and belongs to the class of canonical form associated with this algebra with and . The corresponding invariants and invariant representation for are Plugging in the values of , and we arrive at which gives Some manipulations giveSimilarly which givesWe deduce thatFrom , by utilizing the value of , we obtain We deduce thatThus (37), (40), and (42) are solutions of system (33).

() Consider the oscillatory motion of a system of three particles governed by in a region; such that for any small number , where and ; are arbitrary constants whereas is a drag constant. There is a time varying applied force which produces the motion after neglecting the effects of gravity due to small region. System (43) admits the Lie algebra spanned by symmetry generatorsBy using the invertible transformationsthis Lie algebra is transformed into the canonical form of the Lie algebra mentioned in [8]. System (43) is transformed into System (46) is matched with corresponding canonical form by taking and in it, plugging with mentioned in Table 1. The associated differential invariants and canonical representation are Using invariant representations, inverting back into original variables and dropping the bars, we get These are solutions of system (43).

4. Conclusion

In symmetry analysis of differential equations, integrability is an important part of the analysis. In the literature different approaches have been constructed by employing admitted symmetries of differential equations. The classical approach of reductions is not valid in usual way for systems of ordinary differential equations as used for scalar ordinary differential equations (see, e.g., [5, 6, 12]). In this research work, our focus is on systems of three th-order () ordinary differential equations that admit three-dimensional Lie algebras.

By utilizing the approach developed by Ayub et al. in [8, 9, 13], the bases of differential invariants for system of three 2nd-order ODEs admitting three-dimensional Lie algebras which are constructed. This is Lie algebraic approach with special type of invariant representation, obtained by exploiting the interrelationships of the invariants. Moreover canonical forms associated with these invariants for system of three 2nd-order ODEs are obtained and listed in Table 1. Furthermore singularity in invariant structures has also been investigated. In addition an algorithm for construction of systems of three th-order () ordinary differential equations possessing three-dimensional Lie algebras is provided. Integrability of canonical forms presented in Table 1 is investigated in detailed manner. Illustrative examples from mechanics of system of particles are presented to explain the integrability of canonical forms. This classification will provide a framework for classification of system of three 2nd-order ODEs admitting higher dimensional Lie algebras.

Remarks 12. (i) For Table 1 notations are mentioned in Section 2.1 whereas , and are arbitrary functions.
(ii) Both of algebras and possess the same set of invariants plugging with a single equation (are not associated with systems of three 2nd-order ODEs); the canonical forms are not obtained for them. For example corresponds to the following invariants and the related single equation, namely, . The invariant representations of these cases are not listed in Table 1.
(iii) The algebras possess singular invariants structure.

Disclosure

M. N. Qureshi is presently at Allama Iqbal Open University (AIOU), Islamabad, Pakistan.

Conflicts of Interest

The authors declare that they have no conflicts of interest.