Journal of Mathematics

Volume 2013 (2013), Article ID 936859, 8 pages

http://dx.doi.org/10.1155/2013/936859

Research Article

Pulsating Different Curves of Zero Velocity around Triangular Equilibrium Points in Elliptical Restricted Three-Body Problem

A. Narayan¹ and Amit Shrivastava²

¹Department of Mathematics, Bhilai Institute of Technology, Durg 491001, India

²Department of Mathematics, Rungta College of Engineering and Technology, Bhilai 490020, India

Received 30 December 2012; Revised 15 February 2013; Accepted 16 February 2013

Academic Editor: Baoding Liu

Copyright © 2013 A. Narayan and Amit Shrivastava. 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.

The oblateness and the photogravitational effects of both the primaries on the location and the stability of the triangular equilibrium points in the elliptical restricted three-body problem have been discussed. The stability of the triangular points under the photogravitational and oblateness effects of both the primaries around the binary systems Achird, Lyeten, Alpha Cen-AB, Kruger 60, and Xi-Bootis, has been studied using simulation techniques by drawing different curves of zero velocity.

1. Introduction

The present paper is devoted to the analysis of the photogravitational and the oblateness effects of both primaries on the stability of triangular equilibrium points of the planar elliptical restricted three-body problem. The elliptical restricted three-body problem describes the dynamical system more accurately on account that the realistic assumptions of the motion of the primaries are subjected to move along the elliptical orbit. We have attempted to investigate the stability of triangular equilibrium points under the photogravitational and oblateness effects of both the primaries. The bodies of the elliptical restricted three-body problem are generally considered to be spherical in shape, but in actual situations, we have observed that several heavenly bodies are either oblate spheroid or triaxial rigid bodies. The Earth, Jupiter, and Saturn are examples of the oblate spheroid. The lack of sphericity in heavenly bodies causes large perturbation. In addition to the oblateness of heavenly bodies, the triaxiality, the radiation forces of the bodies, the atmospheric drag, and the solar wind are also causes of perturbation.

This motivates studies of stability of triangular equilibrium points under the influence of oblateness and radiation of the primaries in the elliptical restricted three-body problem. The stability of the infinitesimal around the triangular equilibrium points in the elliptical restricted three-body problem described in considerable details is due to [1] and the problem was also studied [2–9]. The stability of motion of infinitesimal around one of the triangular equilibrium points () also depends on and . Nonlinear stability of the triangular equilibrium points of the elliptical restricted three-body problem with or without radiation pressure was studied [10–12]. Furthermore, the nonlinear stability of the infinitesimal in the orbits or the size of the stable region around was studied numerically by [11] and the parametric resonance stability around in the elliptical restricted three-body problem has been studied [10].

The existence of the libration points [13, 14] and their stability in the photogravitational elliptical restricted three-body problem has been studied. The different aspects of the problem in details of elliptical restricted three-body problem have been investigated [15–26].

The influence of the eccentricity of the orbits of the oblate primary bodies with one of the photogravitational effects on the location of collinear and triangular equilibrium points and their stability has been investigated [27–29]. The stability of triangular points in the elliptical restricted three-body problem under the radiating and oblate primaries has been recently discussed [30]. A similar problem has been discussed in detail by applying different techniques to investigate the stability of the system using simulation technique [31].

The present study aims to examine the motion of the infinitesimal body in the elliptical restricted three-body problem, when the primaries are oblate spheroid and are also a source of radiation. We have obtained the coordinate of the triangular equilibrium points of the problem. For the circular problem, the primaries are fixed with respect to uniformly rotating axes and hence the Hamiltonian does not involve time explicitly. But when the primaries move on elliptical orbits, the introduction of nonuniformly rotating and pulsating coordinate system results again in fixed location of the primaries. The elliptical restricted three-body problem generalizes the original circular restricted three-body problems, while some useful problem of circular model still can be satisfied by the elliptical case. The Hamiltonian, however, does not depend explicitly on independent variable in this case. The dimensionless variables are introduced by using the distance between primaries given by

where and are the semimajor axis and the eccentricity of the elliptical orbit of the primaries, moving along each other and is the true anomaly of . A co-ordinate system which rotates with the variable angular velocity is introduced. This angular velocity is given by

where is dimensionless time.

The equation follows from the principal of the conservation of angular momentum in the problem of two bodies formed by the primaries of masses and . This principle is expressed by

where , and and are the products of the universal gravitational constants with the masses of primaries.

The force of radiation is given by

where is the gravitational attraction force; is the radiation pressure; is the mass radiation factor. For simplicity of calculation, we have considered , where .

The present paper deals with the photogravitational and oblateness effects of both the primaries on the stability of infinitesimal by exploiting simulation techniques by drawing the different curves of zero velocity.

The present paper comprises three sections. The first section of the paper describes the variational equation of motion of the system. In the second section, we have described the triangular equilibrium points of the system, and in the third section the paper we have derived an expression which is suitable for tracing of different curves of zero velocity.

The curves of zero velocity around equilibrium points have been presented through simulation techniques, which shows the region of stability. The effect of the oblateness of primaries plays an important role in analyzing the stability of infinitesimal which is obvious from the curves of zero velocity traced.

2. Variational Equation of Motion

The differential equations of motion of infinitesimal in the elliptical restricted three-body problem under the oblate and radiating primaries in barycentric, pulsating, and nondimensional coordinates are represented as follows [31]:

where

where

where denotes the partial differentiation of with respect to , and denotes the differentiation of partially with respect to, where and are the oblateness parameters of the primaries. and are the mass radiation factors due to the source of radiation of the bigger primary and the smaller primary.

The co-ordinates of the triangular equilibrium points and are determined as follows [31]:

Thus by the coordinates of the triangular equilibrium points we obtain up to the first-order terms in the parameter , , , and , which is represented by (9), the location of equilibrium points is shown Figures 1, 2, 3, 4, and 5 for the infinitesimal moving around various binary systems.

Figure 1: Location of triangular points for Achird with , , , and , .

Figure 2: Location of triangular points for Luyten with , , , and , .

Figure 3: Location of triangular points for Achird with , and .

Figure 4: Location of triangular points for Kruger 60 with , , , and , .

Figure 5: Location of triangular points for Xi Bootis with , , , and , .

3. Different Curves of Zero Velocity

In order to discuss the different curves of zero velocity of the infinitesimal in the elliptical restricted three-body problem, when both the primaries are oblate spheroid and radiating, multiplying the first equation of (5) by and the second equation by and adding, we get;

we obtain

Since Ω does not contain the time (true anomaly) explicitly, (10) can be integrated to give

Due to the presence of in the denominator of (11), the equation is not possible to integrate to any defined form. Hence, in elliptical restricted three-body problem, it does not adjust the Jacobi integral of the classical circular problem at least in its usual sense.

The elliptical restricted three-body problem is different from the classical restricted problem in the sense that the Jacobi integral does not exist [16], and energy along any orbit is a time-dependent quantity. As we know no exact, complete, and general solution to the elliptical restricted three-body problem, Ω can be obtained unlike in classical restricted three-body problem, but this mathematical inconvenience is overcome along investigation of certain special cases of the problem based on simplifying the mathematical model under consideration [2]. Now, consider the potential function which is represented as follows:

Hence, depends not only on the position co-ordinate of the infinitesimal but also on an independent variable. We select the initial point and we consider only a part of the trajectory and , where δ is arbitrary sufficiently small time interval, during which the primaries describe the small one. We may define a Jacobi constant in elliptical case as follows:

Equation (13) describes different curves of zero velocity, at each given instant of time of elliptical restricted three-body problem. The zero velocity curves are now pulsating with frequency of the nominal elliptical motion. Therefore, in the planar elliptical restricted three-body problem, the zero velocity curves are obtained from the following equation:

with the help of (14), the different curves of zero velocity have been traced using the software MATLAB 7.1 of the infinitesimal around the binary system Achird, Luyten, Alpha Cen AB, Kruger 60 and Xi Bootis, taking into account various values of oblateness parameters and and critical mass parameters and from Table 1. We have traced different curves of zero velocity of the infinitesimal about triangular equilibrium points. Figures 6, 7, and 8 represent the same around binary system Achird, and likewise Figures 9, 10, and 11 for binary system Luyten, Figures 12, 13, and 14 for binary system Alpha Cen-AB, Figures 15, 16, and 17 for binary system Kruger 60, and Figures 18, 19, and 20 for binary system Xi Bootis.

Table 1: Location of triangular equilibrium points mentioned in the research paper of Jagdish and Aishetu [30].

Binary system
Mass ratio ()
Critical mass value
Oblateness
Location of triangular points

±

0
0
0.604
0.8654

Achird
0.3949
0.9971
0.9997
0.01
0.02
0.5891
0.8567

0.1
0.2
0.4544
0.7785

0
0
0.5
0.866

Luyten
0.5
0.999998
0.99999
0.01
0.02
0.485
0.8574

0.1
0.2
0.35
0.7794

0
0
0.5421
0.8626

Cen-AB
0.4519
0.9971
0.985
0.01
0.02
0.527
0.8538

0.1
0.2
0.3903
0.7744

0
0
0.6063
0.866

Kruger 60
0.3937
0.99992
0.99996
0.01
0.02
0.5913
0.8573

0.1
0.2
0.4563
0.7794

0
0
0.5764
0.8658

Xi Bootis
0.4231
0.9988
0.9998
0.01
0.02
0.5614
0.8571

0.1
0.2
0.4267
0.779

Figure 6: Curve of zero velocity (Archid-1) , , and , .

Figure 7: Curve of zero velocity (Archid-2) , , and , .

Figure 8: Curve of zero velocity (Archid-3) , , and , .

Figure 9: Curve of zero velocity (Luyten-1) , , and , .

Figure 10: Curve of zero velocity (Luyten-2) , , and , .

Figure 11: Curve of zero velocity (Luyten-3) , , and , .

Figure 12: Curve of zero velocity (α CenAB-1) , , and , .

Figure 13: Curve of zero velocity (α CenAB-2) , , and , .

Figure 14: Curve of zero velocity (α CenAB-3) , , and , .

Figure 15: Curve of zero velocity (Kruger60-1) , , and .

Figure 16: Curve of zero velocity (Kruger60-2) , , and , .

Figure 17: Curve of zero velocity (Kruger60-3) , , and , .

Figure 18: Curve of zero velocity (Xi Bootis-1) , , and , .

Figure 19: Curve of zero velocity (Xi Bootis-2) , , and , .

Figure 20: Curve of zero velocity (Xi Bootis-3) , , and , .

Hence, we observe typical behavior of the infinitesimal around the binary system, Achird, Luyten, Alpha Cen-AB, Kruger 60, and Xi Bootis.

4. Discussion and Conclusion

The oblateness and photogravitational effects of the primaries on the location and the stability of the triangular Lagrangian points in elliptical restricted three-body problem have been studied. The problem was studied under the assumption that the eccentricity of the orbit of the gravitating bodies is small. The oblateness of the more massive primary does not affect the motion of the smaller primary due to its large mass, whereas it affects the motion of infinitesimal body.

The differential equation governing the motion and stability of triangular equilibrium points of the elliptical restricted three-body problem under the oblate and radiating primaries has been analyzed, and configurations of the triangular equilibrium points are described. The stability of the triangular points under the photogravitational and oblateness effects of both the primaries around the binary systems Achird, Lyeten, Aipha Cen-AB, Kruger 60, and Xi Bootis have been studied using simulation technique by drawing different curves of zero velocity around triangular equilibrium point. It is observed that the region within the curves, the infinitesimal will remain stable.

References

J. M. A. Danby, “Stability of the triangular points in the elliptic restricted problem of three bodies,” The Astronomical Journal, vol. 69, pp. 165–172, 1964.

M. K. Ammar, “The effect of solar radiation pressure on the Lagrangian points in the elliptic restricted three-body problem,” Astrophysics and Space Science, vol. 313, no. 4, pp. 393–408, 2008.

A. Bennett, “Characteristic exponents of the five equilibrium solutions in the elliptically restricted problem,” Icarus, vol. 4, no. 2, pp. 177–187, 1965.

V. V. Markellos, E. Perdios, and P. Labropoulou, “Linear stability of the triangular equilibrium points in the photogravitational elliptic restricted problem, I,” Astrophysics and Space Science, vol. 194, no. 2, pp. 207–213, 1992.

R. Meire, “The stability of the triangular points in the elliptic restricted problem,” Celestial Mechanics, vol. 23, no. 1, pp. 89–95, 1981.

E. Rabe, “Two new classes of periodic Trojan liberation in the elliptical restricted problem and there stabilities,” in Periodic Orbits, Stability and Resonance, G. E. O. Giacaglia, Ed., pp. 33–44, Reidel Publication Company, Dordrecht, The Netherlands, 1970.

E. Rabe, “Elliptical restricted problem : fourth-order stability analysis of the triangular points,” in Recent Advances in Dynamical Astronomy, E. D. Toply and V. Szebehely, Eds., pp. 155–160, D. Reidel Publication Company, Dordrecht, The Netherlands, 1973.

G. E. Roberts, “Linear stability of the elliptic Lagrangian triangle solutions in the three-body problem,” Journal of Differential Equations, vol. 182, no. 1, pp. 191–218, 2002.

A. S. Zimovshchikov and V. N. Tkhai, “Instability of libration points and resonance phenomena in the photogravitational elliptic restricted three-body problem,” Solar System Research, vol. 38, no. 2, pp. 155–163, 2004.

B. Érdi, E. Forgács-Dajka, I. Nagy, and R. Rajnai, “A parametric study of stability and resonances around in the elliptic restricted three-body problem,” Celestial Mechanics & Dynamical Astronomy, vol. 104, no. 1-2, pp. 145–158, 2009.

J. Gyorgrey, “On the non-linear stability of motion's around L₅ in the elliptic restricted problem of the three bodies,” Celestial Mechanics, vol. 36, no. 3, pp. 281–285, 1985.

V. Kumar and R. K. Choudhry, “Nonlinear stability of the triangular libration points for the photogravitational elliptic restricted problem of three bodies,” Celestial Mechanics & Dynamical Astronomy, vol. 48, no. 4, pp. 299–317, 1990.

S. N. Khasan, “Three dimensional periodic solutions to the photogravitational Hill problem,” Cosmic Research, vol. 34, no. 5, pp. 299–317, 1990.

S. N. Khasan, “Liberation solutions to the photogravitational restricted three-body problem,” Cosmic Research, vol. 34, no. 2, pp. 146–151, 1996.

C. Pinyol, “Ejection-collision orbits with the more massive primary in the planar elliptic restricted three-body problem,” Celestial Mechanics & Dynamical Astronomy, vol. 61, no. 4, pp. 315–331, 1995.

L. Floria, “On an analytical solutions in the planar elliptic restricted three body problem,” Monografías del Seminario Matemático García de Galdeano, vol. 31, pp. 135–144, 2004.

P. P. Hallan and N. Rana, “The existence and stability of equilibrium points in the Robe's restricted three-body problem,” Celestial Mechanics & Dynamical Astronomy, vol. 79, no. 2, pp. 145–155, 2001.

W. H. Jefferys, “A new class of periodic solutions of the three dimensional restricted problem,” The Astronomical Journal, vol. 71, no. 2, pp. 99–102, 1966.

W. H. Jefferys, “Doubly symmetric periodic orbits in the three-dimensional restricted problem,” The Astronomical Journal, vol. 70, no. 6, pp. 393–394, 1965.

A. P. Markeev, “On one special case of parametric resonance in problems of celestial mechanics,” Astronomy Letters, vol. 31, no. 5, pp. 350–356, 2005.

A. P. Markeev, Libration Points in Celestial Mechanics and Cosmodynamics, Nauka, Moscow, Russia, 1978.

D. Selaru and C. Cucu-Dumitreacu, “An analysis asymptotical solution in the three body problem,” Romanian Astron. Journal, vol. 4, no. 1, pp. 59–67, 1994.

D. Şelaru and C. Cucu-Dumitrescu, “Infinitesimal orbits around Lagrange points in the elliptic, restricted three-body problem,” Celestial Mechanics & Dynamical Astronomy, vol. 61, no. 4, pp. 333–346, 1995.

T. Schauner, “Die Bevegung in der Nach der Dreieckspu mkte des elliptischen eingeschrinkten Dreikorpen problems,” Celestial Mechanics, vol. 3, no. 2, pp. 189–196, 1971.

Z. Sándor and B. Érdi, “Symplectic mapping for Trojan-type motion in the elliptic restricted three-body problem,” Celestial Mechanics & Dynamical Astronomy, vol. 86, no. 4, pp. 301–319, 2003.

Z. Sàndor, R. Balla, F. Téger, et al., “Short time lyapunov, indicator in the restricted three-body problem,” Celestial Mechanics and Dynamical Astronomy, vol. 79, no. 1, pp. 29–40, 2001.

A. Narayan and C. Ramesh Kumar, “Effects of photogravitational and oblantensss on the triangular Lagrangian points in the elliptical restricted three body problem,” International Journal of Pure and Applied Mathematics, vol. 68, no. 2, pp. 201–224, 2011.

A. Narayan and C. Ramesh, “Stability of triangular equilibrium points in elliptical restricted three body problem under the effects of photogravitational and oblateness of primaries,” International Journal of Pure and Applied Mathematics, vol. 70, no. 5, pp. 735–754, 2011.

A. Narayan and C. Ramesh, “Existence and stability of collinear equilibrium points in elliptical restricted three body problem under the effects of photogravitational and oblateness primaries,” International Journal of Pure and Applied Mathematics, vol. 80, no. 4, pp. 477–494, 2012.

S. Jagdish and U. Aishetu, “On the stability of triangular points in the elliptical R3BP under the radiating and oblate primaries,” Astrophysics and Space Science, vol. 341, no. 2, pp. 349–358, 2012.

A. Narayan and S. Amit, “Effects of oblateness and radiation of primaries on the equilibrium points in the ellipted restricted three body problem,” International Journal of Mathematical Science, vol. 32, no. 10, pp. 330–345, 2012.