Advances in Astronomy

Advances in Astronomy / 2015 / Article

Research Article | Open Access

Volume 2015 |Article ID 850252 | 14 pages | https://doi.org/10.1155/2015/850252

Nonlinear Stability of the Triangular Libration Points for Radiating and Oblate Primaries in CR3BP in Nonresonance Condition

Academic Editor: Zdzislaw E. Musielak
Received01 Apr 2015
Revised17 Jun 2015
Accepted21 Jul 2015
Published24 Aug 2015

Abstract

This paper investigates the existence of resonance and nonlinear stability of the triangular equilibrium points when both oblate primaries are luminous. The study is carried out near the resonance frequency, satisfying the conditions , and in circular cases by the application of Kolmogorov-Arnold-Moser (KAM) theory. The study is carried out for the various values of radiation pressure and oblateness parameters in general. It is noticed that the system experiences resonance at for different values of radiation pressures and oblateness parameter. The case corresponds to the boundary region of the stability for the system. It is found that, except for some values of the radiation pressure, and oblateness parameters and for , the triangular equilibrium points are stable.

1. Introduction

The three-body problem, which describes three masses interacting through Newtonian gravity, has attracted many scientists for more than 300 years. The three-body problem is one of the most challenging problems in the history of science. In celestial mechanics the general three-body problem deals with gravitationally interacting astronomical bodies and intends to predict their motions. The restricted three-body problem (R3BP) is a particular case of the general three-body problem. The problem is restricted in the sense that the infinitesimal mass does not influence the motion of two gravitating primaries but is rather influenced by them. In R3BP, there are two possibilities; namely, the two bodies with dominant masses move around their common centre of mass along either circular or elliptic orbits, which lead to the respective circular or elliptic restricted three-body problems (ER3BP). The ER3BP describes the three-dimensional motion of small particle called the infinitesimal mass under the gravitational force of two finite bodies called primaries around their common centre of mass. The ER3BP generalizes the original circular restricted three-body problems, while some useful properties of circular model still can be satisfied in the elliptical case. The ER3BP describes the dynamical system more accurately as the primaries move along the elliptical orbit.

Modern applications of the three-body problem and restricted three-body problem have been extended to include the Earth, the Moon, and the artificial satellites as well as recently discovered exoplanets. The study of the stability of an elliptical or circular restricted three- body problem of the Hamiltonian system is generally performed by using Kolmogorov-Arnold-Moser (KAM) theorem. The KAM theory is very useful for studying the global stability in three-body problem. The KAM theorem is suitable for the motion undergoing small perturbations which preserves the features of the unperturbed motion.

The nonlinear stability for resonance as well as for the nonresonance cases of the triangular libration points, taking one of the bodies as radiating, was studied by Manju and Choudhry [1]. Kumar and Choudhry [2] investigated the stability of the triangular libration points for nonresonance as well as resonance case, taking both the bodies as radiating in circular restricted three-body problem in the presence of the third and fourth order resonance. Bhatnagar et al. [3] discussed the nonlinear stability of the triangular equilibrium points in circular restricted three bodies, considering bigger primary as a source of radiation. The nonlinear stability of the triangular Lagrangian points, considering the bigger primary as oblate spheroid in circular case, was examined by Markellos et al. [4]. Recently Narayan and Singh [5] studied the nonlinear stability of higher order for both radiating primaries and found that binary systems are stable.

The detailed description and the behavior of equilibrium points in ER3BP are touched upon by Danby [6], Bennett [7], Szebebely [8], Markeev [9], Selaru and Cucu-Dumitrescu [10], and Hallan and Rana [11]. The influence of the eccentricity of the orbits of the primaries with or without radiation pressure on the existence and stability of the equilibrium points was studied by Györgyey [12], Grebenikov [13], Kumar and Choudhry [2, 14], Markellos [4, 15], Sahoo and Ishwar [16], Roberts [17], Zimovshchikov and Tkhai [18], Ammar [19], Érdi et al. [20], Kumar and Ishwar [21], Singh and Umar [22, 23], Usha et al. [24], and Narayan and Singh [5, 25, 26].

It may be noted that the case when the frequencies are equal to zero or are equal to each other usually corresponds to the boundary of stability of the linear system, since, in the absence of oblateness parameters and radiation pressures the critical value of mass ratio denoted by is . Also, if it is considered that , then it is found that the inequality , provided ( are integers), is invalid for , which gives rise to the resonance cases.

The present paper investigates the existence of the resonance and the stability of the infinitesimal mass about the triangular equilibrium points by taking both primaries as radiating and oblate in nonresonance condition. The existence of the resonance and the stability of infinitesimal near the resonance frequency satisfying the conditions , , are studied in the circular cases. The study is carried out at various values of radiation pressures and oblateness parameters.

This paper has been organized in various sections.

Section 1 gives introduction, Section 2 describes the equations of motion of the problem, and Section 3 deals with the characteristics roots and first order stability of the triangular equilibrium points. The existence of resonance is discussed in Section 4, while Section 5 deals with normalization and higher order stability of the libration points. Finally Section 6 summarizes the discussion and conclusion of the paper.

2. Equations of Motion

The differential equations of the motion of the infinitesimal mass in elliptical restricted three-body problem under radiating primaries in pulsating system as given by Narayan and Shrivastava [27] arewhere the force function is defined aswhere    is the radiation pressure, is a true anomaly of the primaries, , are the oblateness parameter, is the eccentricity of the orbits, and is the mass ratio defined asThe coordinates of the triangular equilibrium points and as given by Narayan and Shrivastava [27] are

3. Characteristics Roots and First Order Stability of the Triangular Equilibrium Points

The stability of the elliptical restricted three-body problem is restricted to planar case only. Since the nature of stability about and is similar it is sufficient to study the stability only about . The Hamiltonian as described by Narayan and Shrivastava [27] is given bywhere , , , are the variations in the coordinates .

By substituting    we haveNow expand the Hamiltonian function aswhere is the sum of the terms of the th degree which is homogeneous in the variables , , , . is constant and . Consider where Restricting to alone, the characteristics equation can be given in the formwhere , , and are given by ((9), (10), and (11)). After further calculations the characteristics equation reduces to the formIf and are the frequencies then putting in (15), roots can be written asThe correlation between and and is shown in Figures 15. It is found that for fixed values of radiation pressure increases with an increase in whereas decreases.

4. Existence of Resonance in Circular Cases

In order to discuss the existence of resonance, we consider the following three cases.

Case 1. The first case is when ; that is, Solving the above equation we obtain that is,Solving for we obtain Since , the positive sign is inadmissible. Hence the region of stability in the first approximation can be written asThus, the value of responsible for stable equilibrium points is given by It is clear from (22) that in the absence of radiation pressures and oblateness parameters the critical value of mass ratio is when , which usually corresponds to a boundary of the region of stability of the system.

Case 2. It is when ; that is,that is,Solving the above equation for the resonance value is obtained as

Case 3. It is when ; that is,that is,Solving the above equation for the resonance value is obtained as Table 1 shows the values of corresponding to , , and . Figures 610 show correlation between and by varying , , and . Correlation between and is depicted in Figure 11 by varying and taking and . Figure 12 shows combined figures between and by taking , , and , 0.02, 0.03.



0.010.010.010.020.03670.0320.0129
0.010.010.020.030.03970.02500.0139
0.010.010.040.050.04730.02970.0165
0.010.020.010.020.03740.02360.0131
0.010.020.020.030.04050.02550.0142
0.010.020.040.050.04850.03050.0169
0.010.030.010.020.03670.02320.0129
0.010.030.020.030.04410.02780.0154
0.010.030.040.050.04980.03120.173
0.010.040.010.020.03900.02460.0137
0.010.040.020.030.04230.02660.0148
0.010.040.040.050.05110.03210.0178
0.010.050.010.020.03980.02510.0139
0.010.050.020.030.04330.02720.0151
0.010.050.040.050.05250.03290.0182
0.010.060.010.020.04060.02560.0142
0.010.060.020.030.04430.02790.0155
0.010.060.040.050.05400.03380.0187
0.010.070.010.020.04150.02620.0145
0.010.070.020.030.04530.02850.0158
0.010.070.040.050.05560.03480.0193
0.010.080.010.020.04250.02670.0149
0.010.080.020.030.04650.02920.0162
0.010.080.040.050.05730.03580.0198
0.010.090.010.020.04340.02730.0152
0.010.090.020.030.04760.02990.0166
0.010.090.040.050.05900.03690.0204
0.870.790.010.020.18350.10740.0572
0.870.790.020.030.31310.16480.0846
0.870.790.040.05
0.720.290.010.020.03510.02210.0123
0.720.290.020.030.03780.02380.0133
0.720.290.040.050.04090.02580.0143
0.620.390.010.020.04760.02990.0166
0.620.390.020.030.05270.03310.0183
0.620.390.040.050.06710.04180.231
0.520.490.010.020.07420.04610.0254
0.520.490.020.030.08750.05400.0296
0.520.490.040.050.13670.08230.0445
0.820.290.010.020.03250.02050.0115
0.820.290.020.030.03480.02200.0122
0.820.290.040.050.04050.02550.0142
0.920.290.010.020.03030.01920.0107
0.920.290.020.030.03230.02040.0114
0.920.290.040.050.03710.02340.0130

5. Normalization and Higher Order Stability of the Libration Points

In order to investigate the stability the Hamiltonian is normalized by Birkhoff’s method to the following form: where If is Hamiltonian of 2nd order, defined by (8), (9), (10), and (11) and is of a positive definite form, then the equilibrium position is stable by virtue of Lyapunov theorem [28]. Otherwise the problem of stability can be solved by KAM theorem as given by Arnold [29, 30]. To apply KAM theorem, linear canonical transformation of variations as given by Manju and Choudhry [1] is used which is given as where whereThe transformation (31) reduces the Hamiltonian to the following form:The Hamiltonian can be expanded in the following form:Similarly can be expanded asThe coefficient of third and fourth order terms of and given by Kumar and Choudhry [2] can be listed as Again the Hamiltonian is reduced to a more convenient form which is suitable for further investigation by using the following canonical transformation:Thus, the Hamiltonian (34) may be written asAs explained by Kumar and Choudhry [2], consider the following.

If then The other ten coefficients of third order terms of (39) are obtained by the formula