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π‘¦ξ…žξ…ž+Ξ©2𝑦=βˆ’π›½(π‘¦βˆ’π‘§)βˆ’π›Όπ‘¦3,π‘§ξ…žξ…ž+Ξ©2𝑧=βˆ’π›½(π‘§βˆ’π‘¦)βˆ’π›Όπ‘§3,(1.1) 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 [3] 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 π‘¦ξ…žξ…ž=𝑦2+𝑧2,π‘§ξ…žξ…ž=𝑦; 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π‘¦ξ…žξ…žξ€·π‘¦+πœ–2ξ€Έπ‘¦βˆ’1ξ…žξ€·π‘§+𝑦=𝐴(π‘§βˆ’π‘¦)+π΅ξ…žβˆ’π‘¦ξ…žξ€Έ,π‘§ξ…žξ…žξ€·π‘§+πœ–2ξ€Έπ‘§βˆ’1ξ…žξ€·π‘¦+𝑧=𝐴(π‘¦βˆ’π‘§)+π΅ξ…žβˆ’π‘§ξ…žξ€Έ,(1.2) where 𝐴 and 𝐡 are constants and πœ– is a small parameter and no variational problem exists which can be verified from Douglas [5]. The interested reader is referred to the interesting paper cited above for the classification of Lagrangians in which Douglas [5] 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 [12], 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 [13]. 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 [14]. Furthermore, a new perturbation method based on integrating vectors and multiple scales for perturbed systems of ODEs is presented in [15]. By using this perturbation method a strongly nonlinear forced oscillator based on integrating factors is studied which yields the approximate first integrals [16]. For an account of this theory, the reader is referred to an interesting book [17].

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 𝑒=(𝑒1,𝑒2)=(𝑦,𝑧), the dependent variable with coordinates 𝑦 and 𝑧. The derivatives of 𝑒𝛼 with respect to π‘₯ are denoted by𝑒𝛼π‘₯=𝐷π‘₯(𝑒𝛼),𝛼=1,2,(2.1) in which𝐷π‘₯=πœ•πœ•π‘₯+𝑒𝛼π‘₯πœ•πœ•π‘’π›Ό+𝑒𝛼π‘₯π‘₯πœ•πœ•π‘’π›Όπ‘₯+β‹―,(2.2) 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 π‘’ξ…ž=(𝑒1,𝑒2)=(π‘¦ξ…ž,π‘§ξ…ž), and all higher-order derivatives are denoted by π‘’ξ…ž,π‘’ξ…žξ…ž,π‘’ξ…žξ…žξ…ž,….

The following definitions and results are easily adaptable from the literature (see, e.g., [10–12, 18–20]).(1)The Lie-BΓ€cklund operator is πœ•π‘‹=πœ‰πœ•π‘₯+πœ‚π›Όπœ•πœ•π‘’π›Ό+πœπ›Όπ‘₯πœ•πœ•π‘’π›Όπ‘₯,(2.3) where πœπ›Όπ‘₯ are defined by πœπ›Όπ‘₯=𝐷π‘₯(πœ‚π›Ό)βˆ’π‘’π›Όπ‘₯𝐷π‘₯(πœ‰),𝛼=1,2.(2.4)(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 𝐷𝑋𝐿+π‘₯πœ‰ξ€ΈπΏ=𝐷π‘₯(𝐡),(2.5) 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 𝐷𝑋𝐿+π‘₯πœ‰ξ€Έξ€·πœ‚πΏ=π›Όβˆ’πœ‰π‘’π›Όπ‘₯𝛿𝐿𝛿𝑒𝛼+𝐷π‘₯(𝐡),(2.6) for a suitable function 𝐡(π‘₯,𝑒), where 𝛿𝛿𝑒𝛼=πœ•πœ•π‘’π›Όβˆ’π·π‘₯πœ•πœ•π‘’π›Ό+𝐷2π‘₯πœ•πœ•π‘’π›Όπ‘₯βˆ’π·3π‘₯πœ•πœ•π‘’π›Όπ‘₯π‘₯+β‹―,𝛼=1,2(2.7) 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 ξ‚Έξ€·πœ‚πΌ=π΅βˆ’πœ‰πΏ+π›Όβˆ’πœ‰π‘’π›Όπ‘₯ξ€Έπœ•πΏπœ•π‘’π›Όπ‘₯ξ‚Ή,𝛼=1,2(2.8) 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 𝛿𝐿/𝛿𝑦=0 and 𝛿𝐿/𝛿𝑧=0 is1𝐿=2π‘¦ξ…ž2+12π‘§ξ…ž2βˆ’Ξ©22𝑦2+𝑧2ξ€Έβˆ’π›½2𝑦2+𝑧2ξ€Έβˆ’π›Ό4𝑦4+𝑧4ξ€Έ+𝛽𝑦𝑧.(3.1) 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πœ‰π‘¦=0,πœ‰π‘§πœ‚=0,(3.2)1π‘¦βˆ’12πœ‰π‘₯=0,πœ‚2π‘§βˆ’12πœ‰π‘₯=0,πœ‚1𝑧+πœ‚2π‘¦πœ‚=0,(3.3)1π‘₯=𝐡𝑦,πœ‚2π‘₯=π΅π‘§πœ‚,(3.4)1ξ€·βˆ’Ξ©2π‘¦βˆ’π›½π‘¦βˆ’π›Όπ‘¦3ξ€Έ+𝛽𝑧+πœ‚2ξ€·βˆ’Ξ©2π‘§βˆ’π›½π‘§βˆ’π›Όπ‘§3ξ€Έ+𝛽𝑦(3.5)+πœ‰π‘₯ξ‚ƒβˆ’12𝑦2+𝑧2𝛽+Ξ©2ξ€Έβˆ’π›Ό4𝑦4+𝑧4ξ€Έξ‚„+𝛽𝑦𝑧=𝐡π‘₯.(3.6) After some simple calculations, (3.2)–(3.4) give rise toπœ‚πœ‰=𝑓(π‘₯),1=12π‘“ξ…žπ‘¦βˆ’π‘§π‘1πœ‚+𝑔(π‘₯),2=12π‘“ξ…žπ‘§+𝑦𝑏11+β„Ž(π‘₯),𝐡=4π‘“ξ…žξ…žξ€·π‘¦2+𝑧2ξ€Έ+π‘¦π‘”ξ…ž(π‘₯)+π‘§β„Žξ…ž(π‘₯)+𝑇(π‘₯),(3.7) where 𝑏1 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:𝛼𝑏11=0,𝛽𝑓′=0,(3.8)𝛼𝑓′(π‘₯)=0,𝛼𝑔(π‘₯)=0,π›Όβ„Ž(π‘₯)=0,(3.9)4π‘“ξ…žξ…žξ…ž+π›½π‘“ξ…ž+Ξ©2π‘“ξ…žβˆ’π›½π‘11=0,(3.10)4π‘“ξ…žξ…žξ…ž+π›½π‘“ξ…ž+Ξ©2π‘“ξ…ž+𝛽𝑏1𝑔=0,(3.11)ξ…žξ…žξ€·(π‘₯)+𝛽+Ξ©2ξ€Έβ„Žπ‘”(π‘₯)βˆ’π›½β„Ž(π‘₯)=0,(3.12)ξ…žξ…žξ€·(π‘₯)+𝛽+Ξ©2ξ€Έβ„Ž(π‘₯)βˆ’π›½π‘”(π‘₯)=0,(3.13)𝑇′(π‘₯)=0.(3.14) Equation (3.14) results in𝑇(π‘₯)=𝑐0,(3.15) where π‘π‘œ is a constant.

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

Case 1. 𝛼=0, 𝛽=0.
The solution of system (3.8)–(3.13) with 𝛼=0 and 𝛽=0 yields 𝑓(π‘₯)=𝑏2+𝑏3cos2Ξ©π‘₯+𝑏4𝑔sin2Ξ©π‘₯,(π‘₯)=𝑏5cosΞ©π‘₯+𝑏6sinΞ©π‘₯,β„Ž(π‘₯)=𝑏7cosΞ©π‘₯+𝑏8sinΞ©π‘₯,(3.16) where 𝑏1,…,𝑏8 are constants.
Now, (3.7) together with (3.16) gives πœ‰=𝑏2+𝑏3cos2Ξ©π‘₯+𝑏4πœ‚sin2Ξ©π‘₯,(3.17)1=ξ€·βˆ’π‘3sin2Ξ©π‘₯+𝑏4ξ€Έcos2Ξ©π‘₯Ξ©π‘¦βˆ’π‘1𝑧+𝑏5cosΞ©π‘₯+𝑏6πœ‚sinΞ©π‘₯,(3.18)2=ξ€·βˆ’π‘3sin2Ξ©π‘₯+𝑏4ξ€Έcos2Ξ©π‘₯Ω𝑧+𝑏1𝑦+𝑏7cosΞ©π‘₯+𝑏8sinΞ©π‘₯,(3.19)𝐡=βˆ’Ξ©2𝑏3cos2Ξ©π‘₯+𝑏4𝑦sin2Ξ©π‘₯ξ€»ξ€·2+𝑧2ξ€Έξ€Ί+π‘¦Ξ©βˆ’π‘5sinΞ©π‘₯+𝑏6ξ€»ξ€ΊcosΞ©π‘₯+π‘§Ξ©βˆ’π‘7sinΞ©π‘₯+𝑏8ξ€»cosΞ©π‘₯+π‘π‘œ.(3.20)
Setting one of the constants equal to one and the rest to zero, we find the following Noether point symmetry operators and gauge terms: 𝑋1πœ•=βˆ’π‘§πœ•πœ•π‘¦+π‘¦π‘‹πœ•π‘§,𝐡=0,2=πœ•π‘‹πœ•π‘₯,𝐡=0,3πœ•=cos2Ξ©π‘₯πœ•πœ•π‘₯βˆ’π‘¦Ξ©sin2Ξ©π‘₯πœ•πœ•π‘¦βˆ’π‘§Ξ©sin2Ξ©π‘₯πœ•π‘§,𝐡=βˆ’Ξ©2𝑦2+𝑧2𝑋cos2Ξ©π‘₯,4πœ•=sin2Ξ©π‘₯πœ•πœ•π‘₯+𝑦Ωcos2Ξ©π‘₯πœ•πœ•π‘¦+𝑧Ωcos2Ξ©π‘₯πœ•π‘§,𝐡=βˆ’Ξ©2𝑦2+𝑧2𝑋sin2Ξ©π‘₯,5πœ•=cosΞ©π‘₯π‘‹πœ•π‘¦,𝐡=βˆ’Ξ©π‘¦sinΞ©π‘₯,6πœ•=sinΞ©π‘₯π‘‹πœ•π‘¦,𝐡=Ω𝑦cosΞ©π‘₯,7πœ•=cosΞ©π‘₯π‘‹πœ•π‘§,𝐡=βˆ’Ξ©π‘§sinΞ©π‘₯,8πœ•=sinΞ©π‘₯πœ•π‘§,𝐡=Ω𝑧cosΞ©π‘₯.(3.21)

Case 2. 𝛼=0, 𝛽≠0.
The utilization of solution of (3.8)–(3.13) together with (3.7) results in 𝑏1=0,πœ‰=𝑓(π‘₯)=𝑐1,πœ‚(3.22)1𝑐=𝑔(π‘₯)=22𝑐cosΞ©π‘₯+32sinΞ©π‘₯+𝑐4cos2𝛽+Ξ©2ξ‚Άπ‘₯+𝑐5sin2𝛽+Ξ©2ξ‚Άπœ‚π‘₯,(3.23)2𝑐=β„Ž(π‘₯)=22𝑐cosΞ©π‘₯+32sinΞ©π‘₯βˆ’π‘4cos2𝛽+Ξ©2ξ‚Άπ‘₯βˆ’π‘5sin2𝛽+Ξ©2ξ‚ΆΞ©π‘₯,(3.24)𝐡=2ξ€·(𝑦+𝑧)βˆ’π‘2sinΞ©π‘₯+𝑐3ξ€ΈβˆšcosΞ©π‘₯+(π‘¦βˆ’π‘§)2𝛽+Ξ©2Γ—ξ‚ƒβˆ’π‘4ξ‚€βˆšsin2𝛽+Ξ©2π‘₯+𝑐5ξ‚€βˆšcos2𝛽+Ξ©2π‘₯ξ‚„+π‘π‘œ,(3.25) where 𝑐1,…,𝑐5 are constants and 2𝛽+Ξ©2>0. If 2𝛽+Ξ©2<0 then √cos(2𝛽+Ξ©2)π‘₯ is replaced with √cosh(2𝛽+Ξ©2)π‘₯ and √sin(2𝛽+Ξ©2)π‘₯ with √sinh(2𝛽+Ξ©2)π‘₯.
The Noether point symmetries for 𝛼=0, 𝛽≠0 are 𝑋1=πœ•π‘‹πœ•π‘₯,𝐡=0,(3.26)2=12ξ‚΅πœ•cosΞ©π‘₯+πœ•πœ•π‘¦ξ‚ΆΞ©πœ•π‘§,𝐡=βˆ’2𝑋(𝑦+𝑧)sinΞ©π‘₯,(3.27)3=12ξ‚΅πœ•sinΞ©π‘₯+πœ•πœ•π‘¦ξ‚ΆΞ©πœ•π‘§,𝐡=2𝑋(𝑦+𝑧)cosΞ©π‘₯,(3.28)4=cos2𝛽+Ξ©2ξ‚Άπ‘₯ξ‚΅πœ•βˆ’πœ•πœ•π‘¦ξ‚Άξ”πœ•π‘§,𝐡=βˆ’2𝛽+Ξ©2(π‘¦βˆ’π‘§)sin2𝛽+Ξ©2𝑋π‘₯,(3.29)5=sin2𝛽+Ξ©2ξ‚Άπ‘₯ξ‚΅πœ•βˆ’πœ•πœ•π‘¦ξ‚Άξ”πœ•π‘§,𝐡=2𝛽+Ξ©2(π‘¦βˆ’π‘§)cos2𝛽+Ξ©2ξ‚Άπ‘₯.(3.30)

Case 3. If 𝛼≠0, 𝛽=0, then from system (3.8)–(3.13), we obtain 𝑏1=0,𝑓(π‘₯)=𝑑1,β„Ž(π‘₯)=0,𝑔(π‘₯)=0,(3.31) where 𝑑1 is a constant. Hence, πœ‰=𝑑1,πœ‚1=0,πœ‚2=0,𝐡=0.(3.32)
We get only one generator for this case as 𝑋1=πœ•πœ•π‘₯,𝐡=0.(3.33)

Case 4. 𝛼≠0, 𝛽≠0.
System (3.8)–(3.13) with 𝛼≠0 and 𝛽≠0 yields the following results: 𝑏1=0,𝑓(π‘₯)=𝑒1,β„Ž(π‘₯)=0,𝑔(π‘₯)=0,(3.34) where 𝑒1 is a constant. Thus, πœ‰=𝑒1,πœ‚1=0,πœ‚2=0,𝐡=0.(3.35)
For this case, the symmetry generator is 𝑋1=πœ•πœ•π‘₯,𝐡=0.(3.36)

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))ξ€Ίξ€·πœ‚πΌ=π΅βˆ’πœ‰πΏ+1βˆ’πœ‰π‘¦ξ…žξ€Έπ‘¦ξ…ž+ξ€·πœ‚2βˆ’πœ‰π‘§ξ…žξ€Έπ‘§ξ…žξ€»,(3.37) 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. 𝛼=0, 𝛽=0.
If we invoke the Noether point symmetry generators and gauge terms from system (3.21) in (3.37), we obtain the following first integrals: 𝐼1=π‘¦ξ…žπ‘§βˆ’π‘¦π‘§ξ…ž,𝐼2=12ξ€·π‘¦ξ…ž2+π‘§ξ…ž2ξ€Έ+Ξ©22𝑦2+𝑧2ξ€Έ,𝐼3=βˆ’Ξ©2𝑦2+𝑧2ξ€Έ1cos2Ξ©π‘₯βˆ’2𝑦cos2Ξ©π‘₯ξ…ž2+π‘§ξ…ž2βˆ’Ξ©2𝑦2+𝑧2+ξ€·ξ€Έξ€»π‘¦π‘¦ξ…ž+π‘§π‘§ξ…žξ€Έξ‚€π‘¦Ξ©sin2Ξ©π‘₯+ξ…ž2+π‘§ξ…ž2𝐼cos2Ξ©π‘₯,4=βˆ’Ξ©2𝑦2+𝑧2ξ€Έ1sin2Ξ©π‘₯βˆ’2𝑦sin2Ξ©π‘₯ξ…ž2+π‘§ξ…ž2βˆ’Ξ©2𝑦2+𝑧2βˆ’ξ€·ξ€Έξ€»π‘¦π‘¦ξ…ž+π‘§π‘§ξ…žξ€Έξ€·π‘¦Ξ©cos2Ξ©π‘₯+ξ…ž2+π‘§ξ…ž2𝐼sin2Ξ©π‘₯,5=βˆ’Ξ©π‘¦sinΞ©π‘₯βˆ’π‘¦ξ…žπΌcosΞ©π‘₯,6=Ω𝑦cosΞ©π‘₯βˆ’π‘¦ξ…žπΌsinΞ©π‘₯,7=βˆ’Ξ©π‘§sinΞ©π‘₯βˆ’π‘§ξ…žπΌcosΞ©π‘₯,8=Ω𝑧cosΞ©π‘₯βˆ’π‘§ξ…žsinΞ©π‘₯.(3.38) 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. 𝛼=0, 𝛽≠0.
The first integrals for this case are 𝐼1=12𝑦′2+𝑧′2ξ€Έ+Ξ©22𝑦2+𝑧2ξ€Έ+𝛽2𝑦2+𝑧2ξ€ΈπΌβˆ’π›½π‘¦π‘§,2Ξ©=βˆ’2(1𝑦+𝑧)sinΞ©π‘₯βˆ’2𝐼cosΞ©π‘₯(𝑦′+𝑧′),3=Ξ©21(𝑦+𝑧)cosΞ©π‘₯βˆ’2𝐼sinΞ©π‘₯(𝑦′+𝑧′),4=βˆ’(π‘¦βˆ’π‘§)2𝛽+Ξ©2sin2𝛽+Ξ©2π‘₯βˆ’(π‘¦β€²βˆ’π‘§β€²)cos2𝛽+Ξ©2𝐼π‘₯,5=(π‘¦βˆ’π‘§)2𝛽+Ξ©2cos2𝛽+Ξ©2π‘₯βˆ’(π‘¦β€²βˆ’π‘§β€²)sin2𝛽+Ξ©2ξ‚Άπ‘₯.(3.39)

Case 3. 𝛼≠0,𝛽=0.
The simple calculations lead to the following integral: 𝐼1=12𝑦′2+𝑧′2ξ€Έ+Ξ©22𝑦2+𝑧2ξ€Έ+𝛼4𝑦4+𝑧4ξ€Έ.(3.40)

Case 4. 𝛼≠0, 𝛽≠0.
Straightforward manipulations result in the following first integral: 𝐼1=12ξ‚€π‘¦ξ…ž2+π‘§ξ…ž2+12𝑦2+𝑧2𝛽+Ξ©2ξ€Έ+𝛼4𝑦4+𝑧4ξ€Έβˆ’π›½π‘¦π‘§.(3.41) 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) is1𝐿=2π‘¦ξ…ž2+12π‘§ξ…ž2,(4.1) so that system (1.1) can be written as𝛿𝐿𝛿𝑦=Ξ©2𝑦+𝛽(π‘¦βˆ’π‘§)+𝛼𝑦3,𝛿𝐿𝛿𝑧=Ξ©2𝑧+𝛽(π‘§βˆ’π‘¦)+𝛼𝑧3.(4.2) 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πœ‰π‘¦=0,πœ‰π‘§πœ‚=0,(4.3)1π‘¦βˆ’12πœ‰π‘₯=0,πœ‚2π‘§βˆ’12πœ‰π‘₯=0,πœ‚1𝑧+πœ‚2π‘¦πœ‚=0,(4.4)1π‘₯=π΅π‘¦ξ€·Ξ©βˆ’πœ‰2𝑦+𝛽𝑦+𝛼𝑦3ξ€Έπœ‚βˆ’π›½π‘§,(4.5)2π‘₯=π΅π‘§ξ€·Ξ©βˆ’πœ‰2𝑧+𝛽𝑧+𝛼𝑧3ξ€Έπ΅βˆ’π›½π‘¦,(4.6)π‘₯+πœ‚1ξ€·Ξ©2𝑦+𝛽𝑦+𝛼𝑦3ξ€Έβˆ’π›½π‘§+πœ‚2ξ€·Ξ©2𝑧+𝛽𝑧+𝛼𝑧3ξ€Έβˆ’π›½π‘¦=0.(4.7) The solution of (4.3)–(4.6) yieldsπœ‚πœ‰=𝑓(π‘₯),1=12π‘“ξ…žπ‘¦βˆ’π‘§π‘1πœ‚+𝑔(π‘₯),2=12π‘“ξ…žπ‘§+𝑦𝑏11+β„Ž(π‘₯),𝐡=4π‘“ξ…žξ…žξ€·π‘¦2+𝑧2ξ€Έ+π‘¦π‘”ξ…ž(π‘₯)+π‘§β„Žξ…žξ‚ƒ1(π‘₯)+𝑓2ξ€·Ξ©2𝑦+𝛽2+𝑧2ξ€Έ+𝛼4𝑦4+𝑧4ξ€Έξ‚„βˆ’π‘¦π‘§π›½+𝑇(π‘₯),(4.8) where 𝑏1 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. 𝛼=0, 𝛽=0.
For this case πœ‰, πœ‚1, and πœ‚2 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 Ω𝐡=2𝑏22𝑦2+𝑧2ξ€Έβˆ’12Ξ©2𝑏3cos2Ξ©π‘₯+𝑏4𝑦sin2Ξ©π‘₯ξ€»ξ€·2+𝑧2ξ€Έξ€Ί+π‘¦Ξ©βˆ’π‘5sinΞ©π‘₯+𝑏6ξ€»ξ€ΊcosΞ©π‘₯+π‘§Ξ©βˆ’π‘7sinΞ©π‘₯+𝑏8ξ€»cosΞ©π‘₯+π‘π‘œ.(4.9) Notice from (3.20) and (4.9) that only the 𝐡's relating to 𝑋2, 𝑋3, and 𝑋4 for the partial Noether's case are different from that for the Noether's approach. The 𝐡's relating to 𝑋2, 𝑋3, and 𝑋4 for the partial Noether case are Ω𝐡=22𝑦2+𝑧2ξ€Έ,Ω𝐡=βˆ’22𝑦2+𝑧2ξ€ΈΞ©cos2Ξ©π‘₯,𝐡=βˆ’22𝑦2+𝑧2ξ€Έsin2Ξ©π‘₯.(4.10)

Case 2. 𝛼=0, 𝛽≠0.
For this case, πœ‰,πœ‚1, and πœ‚2 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 Ω𝐡=2ξ€·(𝑦+𝑧)βˆ’π‘2sinΞ©π‘₯+𝑐3cosΞ©π‘₯+(π‘¦βˆ’π‘§)2𝛽+Ξ©2ξ‚Έβˆ’π‘4sin2𝛽+Ξ©2ξ‚Άπ‘₯+𝑐5cos2𝛽+Ξ©2ξ‚Άπ‘₯ξ‚Ή+𝑐112ξ€·Ξ©2𝑦+𝛽2+𝑧2ξ€Έξ‚„βˆ’π‘¦π‘§π›½+π‘π‘œ.(4.11) 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 𝑋1 is different which is given as follows: 1𝐡=2ξ€·Ξ©2𝑦+𝛽2+𝑧2ξ€Έβˆ’π›½π‘¦π‘§.(4.12)

Case 3. 𝛼≠0, 𝛽=0.
The partial Noether operator of point type when 𝛼≠0, 𝛽=0 is 𝑋1=πœ•/πœ•π‘₯ which is identical to the case of Noether point symmetries but the difference occurs in the gauge term as Ω𝐡=22𝑦2+𝑧2ξ€Έ+𝛼4𝑦4+𝑧4ξ€Έ.(4.13)

Case 4. 𝛼≠0, 𝛽≠0.
The operator for this case is also similar to the case of Noether point symmetries except 𝐡 in (3.36) is different as follows 1𝐡=2ξ€·Ξ©2𝑦+𝛽2+𝑧2ξ€Έ+𝛼4𝑦4+𝑧4ξ€Έβˆ’π›½π‘¦π‘§.(4.14)

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. 𝛼=0, 𝛽=0.
For operators 𝑋1 and 𝑋5 to 𝑋8, we have πœ‰=0 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 𝑋2, 𝑋3, and 𝑋4 for the partial Noether case are different from that for the Noether's approach and πœ‰β‰ 0 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.