`Mathematical Problems in EngineeringVolume 2011, Article ID 831647, 14 pageshttp://dx.doi.org/10.1155/2011/831647`
Research Article

## First Integrals for Two Linearly Coupled Nonlinear Duffing Oscillators

1Centre for Mathematics and Statistical Sciences, Lahore School of Economics, Lahore 53200, Pakistan
2Department of Mathematics, LUMS, School of Science and Engineering, Lahore Cantt 54792, Pakistan
3Centre for Differential Equations, Continuum Mechanics and Applications, University of the Witwatersrand, Johannesburg, Wits 2050, South Africa

Received 6 September 2010; Revised 1 December 2010; Accepted 3 January 2011

Copyright © 2011 R. Naz et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We investigate Noether and partial Noether operators of point type corresponding to a Lagrangian and a partial Lagrangian for a system of two linearly coupled nonlinear Duffing oscillators. Then, the first integrals with respect to Noether and partial Noether operators of point type are obtained explicitly by utilizing Noether and partial Noether theorems for the system under consideration. Moreover, if the partial Euler-Lagrange equations are independent of derivatives, then the partial Noether operators become Noether point symmetry generators for such equations. The difference arises in the gauge terms due to Lagrangians being different for respective approaches. This study points to new ways of constructing first integrals for nonlinear equations without regard to a Lagrangian. We have illustrated it here for nonlinear Duffing oscillators.

#### 1. Introduction

In this paper, we study a system of two linearly coupled nonlinear Duffing oscillators where and are nonlinearity linear coupling parameters, and prime denotes differentiation with respect to . The system of two second-order ordinary differential equations (ODEs) frequently arises in nonlinear oscillations, nonlinear dynamics, relativity, fluid mechanics, and so forth. These oscillators, in general, describe different mechanical systems of practical importance and have two types of characteristics (see, e.g., [1, 2]). The hardening characteristic approaches to linear for small amplitudes and it raises towards infinity when the amplitude tends to certain limits. The softening characteristic performs in a non-monotonic way, and for large amplitude the vibration shape changes to the rectangular shape.

There is great interest in studying the system of two linearly coupled nonlinear Duffing oscillators, especially for first integrals, which are important from a physical point of view and for reduction purposes as well. Of course, a Lagrangian exists for the system of two equations under study. Our purpose is two fold: one is how first integrals can be constructed for ODEs; secondly, we want to investigate the effectiveness of the partial Lagrangian approach.

The Noether symmetries and first integrals are important due to their wide range of applications for Euler-Lagrange equations (see, e.g., Noether  and later works [4, 5]). This theorem is based on the existence of a Lagrangian and variational symmetries. It gives the relationship between equivalence classes of point symmetries and first integrals. There are equations that arise in application that do not admit standard Lagrangians, for example, the simple system of two second-order equations ,; the curve family is nonextremal for this system. The scalar evolution equations also do not admit Lagrangians. Similarly, for the system of two coupled Van der pol oscillators with linear diffusive coupling where and are constants and is a small parameter and no variational problem exists which can be verified from Douglas . The interested reader is referred to the interesting paper cited above for the classification of Lagrangians in which Douglas  has provided the complete solution to the inverse problem for three-dimensional space (system of two second-order ODEs). Now, we raise the question as to how one can construct first integrals for equations which are not variational as mentioned above. It is well known that there are other methods which provide first integrals without making use of a Lagrangian. The most elementary method is the direct method [6, 7] which is often used for constructing first integrals without the variational principle. There are some other approaches as well (see, e.g., [4, 8, 9]) in which the equations can be expressed in characteristic form.

The Noether symmetries and the corresponding first integrals have been subjects of rigorous investigations and play an important role in the reduction of differential equations. Some important works have been done relating to first integrals. In [10, 11], the authors deduce a relation between “symmetries” and first integrals without a variational principle. The classification of Noether point symmetries for a three degrees of freedom Lagrangian system was given by Damianou and Sophocleous in , and the results for one and two degrees of freedom were also reported in their paper. The first integrals associated with Noether and partial Noether operators of point type for a linear system of two second-order ODEs with variable coefficients were constructed by the authors in . They showed that the first integrals resulting from Noether and partial Noether approaches are equivalent. The difference occurs in the gauge terms due to Lagrangians being different for the respective approaches. The classification of partial Noether operators and first integrals for a system with two degrees of freedom was also discussed by Naeem and Mahomed in . Furthermore, a new perturbation method based on integrating vectors and multiple scales for perturbed systems of ODEs is presented in . By using this perturbation method a strongly nonlinear forced oscillator based on integrating factors is studied which yields the approximate first integrals . For an account of this theory, the reader is referred to an interesting book .

In this paper, we show that the first integrals corresponding to Noether and partial Noether operators of point type for a system of two linearly coupled nonlinear Duffing oscillators are equivalent, and the algebras for both cases are isomorphic. We give an elegant way to construct first integrals for such equations with respect to point symmetries. Firstly, we obtain the partial Noether operators of point type and then first integrals are constructed by utilizing a partial Noether's theorem. This approach gives new ways of constructing first integrals for the equations without regard to a Lagrangian as partial Lagrangians do exist for such equations. It would definitely be of interest to obtain the first integrals of the Duffing oscillator system with respect to first-order Lie-Bäcklund operators. This will provide a more complete picture for such systems. Further work can then be done on reductions using these symmetries.

The outline of the work is as follows. In Section 2, the basic definitions are adapted from the literature. The Noether point symmetries and the corresponding first integrals are presented in Section 3. In Section 4, we construct the partial Noether operators of point type and the relating first integrals. Concluding remarks are summarized in Section 5.

#### 2. Preliminaries

Suppose is considered to be the independent variable and , the dependent variable with coordinates and . The derivatives of with respect to are denoted by in which is the total derivative operator with respect to . Note that the summation convention is assumed for repeated indices throughout. The collection of first derivative is denoted by , and all higher-order derivatives are denoted by .

The following definitions and results are easily adaptable from the literature (see, e.g., [1012, 1820]).(1)The Lie-Bäcklund operator is where are defined by (2)A Lie-Bäcklund operator defined in (2.3) is a Noether symmetry generator corresponding to a Lagrangian if there exists a function such that where is the total derivative operator given in (2.2). The function is known as the gauge term.(3)The Lie-Bäcklund operator is known to be a partial Noether operator corresponding to a partial Lagrangian if it is determined from for a suitable function where is the Euler operator.(4)If the Lie-Bäcklund operator is a Noether symmetry generator corresponding to a Lagrangian of Euler-Lagrange differential equations, then is known as the first integral of Euler-Lagrange differential equations.(5)If the Lie-Bäcklund operator is a partial Noether operator with respect to a partial Lagrangian of partial Euler-Lagrange equations, then the constant of motion or first integral can be constructed from (2.8).

#### 3. Noether's Approach

In this section, we derive the Noether point symmetry operators for the coupled nonlinear Duffing oscillators (1.1) corresponding to a standard Lagrangian and construct the first integrals by Noether's theorem.

##### 3.1. Noether Point Symmetry Generators for Coupled Nonlinear Duffing Oscillators

A Lagrangian for system (1.1) satisfying the Euler-Lagrange equations and is The Noether point symmetry operator corresponding to a Lagrangian for a system of two linearly coupled nonlinear Duffing oscillators (1.1) is calculated from the formula (2.5) with respect to some function . The Noether point symmetry determining equations are After some simple calculations, (3.2)–(3.4) give rise to where is an arbitrary constant and , , , and are arbitrary functions of . The insertion of (3.7) in (3.6) and separation with respect to powers of and reduce to the following system: Equation (3.14) results in where is a constant.

In order to solve system (3.8)–(3.13), the following cases arise.

Case 1. , .
The solution of system (3.8)–(3.13) with and yields where are constants.
Now, (3.7) together with (3.16) gives
Setting one of the constants equal to one and the rest to zero, we find the following Noether point symmetry operators and gauge terms:

Case 2. , .
The utilization of solution of (3.8)–(3.13) together with (3.7) results in where are constants and . If then is replaced with and with .
The Noether point symmetries for , are

Case 3. If , , then from system (3.8)–(3.13), we obtain where is a constant. Hence,
We get only one generator for this case as

Case 4. , .
System (3.8)–(3.13) with and yields the following results: where is a constant. Thus,
For this case, the symmetry generator is

The interpretation of the results in case of Noether point symmetries are given below.

Noether point symmetries
Case 1: in this case, we obtain an eight-dimensional Lie algebra.
Case 2: the Lie algebra is five dimensional.
Case 3: we deduce a one-dimensional Lie algebra.
Case 4: the Lie algebra for this case is also one dimensional.

##### 3.2. First Integrals

The first integrals for Noether point symmetry generators for system (1.1) corresponding to the Lagrangian (3.1) are (by (2.8)) where we have used and which hold for all the cases. The first integrals for each case computed with the help of (3.37) are summarized below.

Case 1. , .
If we invoke the Noether point symmetry generators and gauge terms from system (3.21) in (3.37), we obtain the following first integrals: For a system of two second-order ODEs, there are exactly four functionally independent first integrals. What we obtained are seven first integrals. Out of these, there are only four functionally independent first integrals.

Case 2. , .
The first integrals for this case are

Case 3. ,.
The simple calculations lead to the following integral:

Case 4. , .
Straightforward manipulations result in the following first integral: Note that the partial Lagrangian approach is similar to the Lagrangian approach; when one uses point-type operators one may not obtain all the first integrals. A familiar example is the Laplace-Runge-Lenz vector for the well-known Kepler plroblem which comes from a higher symmetry and not a geometric sysmmetric of point type. Here, too, we expect that the existence of an elliptic type integral will arise from a derivative dependent operator which will be of interest to look at in a future work. We only obtain the kinetic plus potential energy integrals here. Presumably, a higher derivative dependent operator could yield further integrals if they exist. These need to be further investigated.

#### 4. Partial Noether's Approach

In this section, we derive the partial Noether operators of point type associated with a partial Lagrangian of (1.1) and then construct the first integrals for (1.1).

##### 4.1. Partial Noether Operators of Point Type for Coupled Nonlinear Duffing Oscillators

A partial Lagrangian for system (1.1) is so that system (1.1) can be written as The partial Noether operator corresponding to a partial Lagrangian for the system (1.1) can be determined from (2.6) with respect to some function . The partial Noether operator of point-type determining equations are The solution of (4.3)–(4.6) yields where is an arbitrary constant and , , , and are arbitrary functions of . The substitution of (4.8) in (4.7) and separation with respect to powers of and reduce to the system of (3.8)–(3.13) that we have obtained in the case of Noether point symmetries and . In case of partial Noether operators of point type the same cases need to be considered as we have for the Noether point symmetry operators.

Case 1. , .
For this case , , and are the same as given in (3.17), (3.18), and (3.19) and, therefore, the partial Noether operators of point type are identical to the Noether point symmetry operators but the gauge terms are different. The gauge term for the partial Noether operators of point type is Notice from (3.20) and (4.9) that only the 's relating to , , and for the partial Noether's case are different from that for the Noether's approach. The 's relating to , , and for the partial Noether case are

Case 2. , .
For this case, ,, and are the same as given in (3.22), (3.23), and (3.24) and therefore the partial Noether operators of point type are similar to the case of Noether point symmetries (3.30) but the gauge terms are different. The gauge term for partial Noether case is Note that in (4.11) all the gauge terms are similar as in the case of Noether symmetries (see (3.25)) except one gauge term relating to is different which is given as follows:

Case 3. , .
The partial Noether operator of point type when , is which is identical to the case of Noether point symmetries but the difference occurs in the gauge term as

Case 4. , .
The operator for this case is also similar to the case of Noether point symmetries except in (3.36) is different as follows

##### 4.2. First Integrals

For the partial Lagrangian (4.1) and also hold and therefore the first integrals can be computed from (3.37). Moreover, since Noether and partial Noether operator of point type are the same for each case, only the term in formula (3.37) will be different from the Noether case.

Case 1. , .
For operators and to , we have and the gauge terms are the same as for the case of Noether point symmetries. Therefore, we obtain the same first integrals. The gauge term 's relating to , , and for the partial Noether case are different from that for the Noether's approach and for these operators; hence, formula (3.37) yields the same first integrals as we have constructed for the Noether's case.

Similarly, for other cases, the first integrals for the partial Noether operators of point type are identical to the case of Noether point symmetry generators.

The interpretation of the results in case of partial Noether operators of point type is as follows.

##### 4.3. Partial Noether Operators of Point Type

The partial Noether operators of point type for the coupled nonlinear Duffing oscillators are the same as the Noether point symmetry operators. The difference arises in the “gauge” terms due to the Lagrangians being different. The algebras for both cases are isomorphic due to the fact that and are independent of derivatives. The first integrals due to the partial Noether operators of point type associated with a partial Lagrangian are equivalent to the case of Noether point symmetries.

#### 5. Concluding Remarks

The Noether point symmetry operators corresponding to a Lagrangian of a system of two linearly coupled nonlinear Duffing oscillators have been constructed. We have computed the first integrals associated with Noether point operators by utilization of Noether's theorem. Then the partial Noether operators of point type and first integrals for the system under study are constructed corresponding to a partial Lagrangian. The partial Euler-Lagrange equations obtained herein are independent of derivatives, and therefore the partial Noether operators become point symmetry operators and the algebras are isomorphic for both the Noether and partial Noether cases. The first integrals corresponding to the Noether and partial Noether operators of point type are equivalent and these form one-, five-, and eight-dimensional Lie algebras. It would be of value to extend the results here to more general symmetries of Lie-Bäcklund type. Moreover, it is of interest to obtain reductions for the system using point or Lie-Bäcklund symmetries.

This study also points to new ways of constructing first integrals for systems without use of a Lagrangian. However, a partial Lagrangian can exist for such equations and one can invoke a formula for the construction of the first integrals. Since, in general, no variational problem exists for equations, one can utilize a partial Lagrangian approach which we have shown here to be effective for the coupled Duffing oscillators. The first integrals still can be constructed by using a Noether-like theorem.

#### References

1. M. F. Dimentberg and A. S. Bratus, “Bounded parametric control of random vibrations,” Proceedings of the Royal Society. A, vol. 456, no. 2002, pp. 2351–2363, 2000.
2. V. N. Pilipchuk, “Transitions from strongly to weakly-nonlinear dynamics in a class of exactly solvable oscillators and nonlinear beat phenomena,” Nonlinear Dynamics, vol. 52, no. 3, pp. 263–276, 2008.
3. E. Noether, “Invariant Variationsprobleme,” Nachrichten von der Königlicher Gesellschaft den Wissenschaft zu Göttingen, Mathematisch-Physikalische Klasse Heft, vol. 2, pp. 235–257, 1918, English translation in Transport Theory and Statical Physics, vol. 1, pp. 186–207, 1971.
4. P. J. Olver, Applications of Lie Groups to Differential Equations, vol. 107 of Graduate Texts in Mathematics, Springer, New York, NY, USA, 1986.
5. J. Douglas, “Solution of the inverse problem of the calculus of variations,” Transactions of the American Mathematical Society, vol. 50, pp. 71–128, 1941.
6. J. Hietarinta, Direct Methods for the Search of the Second Invariant, Report Series, Department of Physical Sciences, University of Turku, Turku, Finland, 1986.
7. H. R. Lewis and P. G. L. Leach, “A direct approach to finding exact invariants for one-dimensional time-dependent classical Hamiltonians,” Journal of Mathematical Physics, vol. 23, no. 12, pp. 2371–2374, 1982.
8. H. Steudel, “Über die Zuordnung zwischen Invarianzeigenschaften und Erhaltungssätzen,” Zeitschrift für Naturforschung A, vol. 17, pp. 129–132, 1962.
9. S. Moyo and P. G. L. Leach, “Symmetry properties of autonomous integrating factors,” SIGMA, vol. 1, 12 pages, 2005.
10. A. H. Kara and F. M. Mahomed, “Noether-type symmetries and conservation laws via partial Lagrangians,” Nonlinear Dynamics, vol. 45, no. 3-4, pp. 367–383, 2006.
11. A. H. Kara, F. M. Mahomed, I. Naeem, and C. Wafo Soh, “Partial Noether operators and first integrals via partial Lagrangians,” Mathematical Methods in the Applied Sciences, vol. 30, no. 16, pp. 2079–2089, 2007.
12. P. A. Damianou and C. Sophocleous, “Classification of Noether symmetries for Lagrangians with three degrees of freedom,” Nonlinear Dynamics, vol. 36, no. 1, pp. 3–18, 2004.
13. I. Naeem and F. M. Mahomed, “Noether, partial Noether operators and first integrals for a linear system,” Journal of Mathematical Analysis and Applications, vol. 342, no. 1, pp. 70–82, 2008.
14. I. Naeem and F. M. Mahomed, “Partial Noether operators and first integrals for a system with two degrees of freedom,” Journal of Nonlinear Mathematical Physics, vol. 15, supplement 1, pp. 165–178, 2008.
15. W. T. van Horssen, “A perturbation method based on integrating vectors and multiple scales,” SIAM Journal on Applied Mathematics, vol. 59, no. 4, pp. 1444–1467, 1999.
16. S. B. Waluya and W. T. van Horssen, “On approximations of first integrals for a strongly nonlinear forced oscillator,” Nonlinear Dynamics, vol. 33, no. 3, pp. 225–252, 2003.
17. F. Verhulst, Nonlinear Differential Equations and Dynamical Systems, Universitext, Springer, Berlin, Germany, 2nd edition, 1996.
18. V. M. Gorringe and P. G. L. Leach, “Lie point symmetries for systems of second order linear ordinary differential equations,” Quaestiones Mathematicae, vol. 11, no. 1, pp. 95–117, 1988.
19. C. Wafo Soh and F. M. Mahomed, “Linearization criteria for a system of second-order ordinary differential equations,” International Journal of Non-Linear Mechanics, vol. 36, no. 4, pp. 671–677, 2001.
20. C. Wafo Soh and F. M. Mahomed, “Canonical forms for systems of two second-order ordinary differential equations,” Journal of Physics A, vol. 34, no. 13, pp. 2883–2911, 2001.