Abstract

This paper investigates the positions and linear stability of an infinitesimal body around the equilibrium points in the framework of the Robe’s circular restricted three-body problem, with assumptions that the hydrostatic equilibrium figure of the first primary is an oblate spheroid and the second primary is an oblate body as well. It is found that equilibrium point exists near the centre of the first primary. Further, there can be one more equilibrium point on the line joining the centers of both primaries. Points on the circle within the first primary are also equilibrium points under certain conditions and the existence of two out-of-plane points is also observed. The linear stability of this configuration is examined and it is found that points near the center of the first primary are conditionally stable, while the circular and out of plane equilibrium points are unstable.

1. Introduction

Robe [1] considered a new kind of restricted three-body problem in which, one of the primaries of mass is a rigid spherical shell, filled with homogenous, incompressible fluid of density ; the second one is a point mass located outside the shell and moving around the mass in a Keplerian orbit; the infinitesimal mass is a small sphere of density , moving inside the shell and is subject to the attraction of and the buoyancy force due to the fluid of the first primary. Further, he discussed the linear stability of an equilibrium point obtained in two cases. In the first case, the orbit of around is circular and in the second case, the orbit is elliptic, but the shell is empty (there is no fluid inside it) or densities of and are equal. Since then various studies (e.g., [24]) under different assumptions have been carried out.

In his study, Robe [1] assumed that the pressure field of the fluid has a spherical symmetry around the center of the shell and he took into account only one out of the three components of the pressure field which is due to the own gravitational field of the fluid . He did not consider the other two components arising from the attraction of and the centrifugal force. Taking care of all these three components of the pressure field, A. R. Plastino and A. Plastino [5] reanalyzed the Robe’s. But in their study, they assumed the hydrostatic equilibrium figure of the first primary as Roche’s ellipsoid (see Figure 1). They found that when the density parameter is taken as zero, every point inside the fluid is an equilibrium point; otherwise the center of the ellipsoid is the only equilibrium point and it is linearly stable.

Hallan and Rana [3] investigated the existence of all equilibrium point and their stability in the Robe’s [1] restricted three-body problem. It was seen that the Robe’s elliptic restricted three-body problem has only one equilibrium point for all values of the density parameter and the mass parameter , while the Robe’s circular restricted three-body problem can have two, three, or infinite numbers of equilibrium points. As regards to the stability of these equilibria, they confirmed the stability result given by Robe [1] of the equilibrium point , whereas triangular and circular points are always unstable. The equilibrium point collinear with the center of the shell and the second primary was found to be stable under some conditions.

Hallan and Mangang [4] studied the Robe’s [1] restricted three-body problem by considering the full buoyancy force as in A. R. Plastino and A. Plastino [5] and assuming the hydrostatic equilibrium figure of the first primary as an oblate spheroid. They derived the pertinent equations of motion and discussed the existence of equilibrium point and their linear stability.

The participating bodies in the classical restricted three-body problem are strictly spherical in shape, but in actual situations several heavenly bodies, such as Saturn and Jupiter, are sufficiently oblate. The minor planets and meteoroids have irregular shape. The lack of sphericity, or the oblateness, of the planet causes large perturbations from a two-body orbit. The motions of artificial Earth satellites are examples of this. Global studies of problems with oblateness have been carried out by many researchers (e.g., [69]).

Therefore, our effort in this paper aims at investigating the equilibrium points and their stability in the Robe’s circular restricted three-body problem when the hydrostatic equilibrium figure of the fluid of the first primary is an oblate spheroid and the second one is an oblate spheroid as well. The model of this study can be used to study the small oscillation of the Earth’s inner core taking into account the Moon’s attraction.

This paper is organized as follows; Section 2 represents the equations of motion; the existence of the equilibrium points is mentioned in Section 3, while Section 4 investigates their linear stability; Section 5 discusses the results obtained; the conclusion is drawn in Section 6.

2. Equation of Motion

Let the first primary be a fluid of density in the shape of an oblate spheroid as assumed by Hallan and Mangang [4]; let the second primary be an oblate body too as Sharma and Subba Rao [6] assumed, which describes a circular orbit around .

We adopt a uniformly rotating coordinate system with origin at the center of mass , pointing towards , with being the orbital plane of coinciding with the equatorial plane of . Then, the equations of motion of the infinitesimal body of density in the coordinate system take the form [4, 6]: where Here is the potential that explains the combined forces upon the infinitesimal mass, denotes the potential due to the fluid mass of the first primary, stands for the potential due to the second primary, is the distance between the primaries, and is the gravitational constant. is the mean motion. and are the equatorial and polar radii of the first and second primary, respectively. stands for the polar moment of inertia, while    are the index symbols. and are the oblateness coefficients of the first and second primaries, respectively.

We choose the unit of mass such that the sum of the masses of the primaries is taken as unity, thus we take , . For the unit of length, we take the distance between the primaries as unity, that is, and the unit of time is also selected such that . With these units and substituting the expression for the potential due to the fluid in the first primary and the potential due to the second oblate primary, the equations of motion (2.1) are recast to the form: where

These above equations of motion of the infinitesimal mass under the framework of the Robe’s circular restricted three-body problem have been obtained by taking into account the shapes of the primaries, the full buoyancy force, the forces due to the gravitational attraction of the second primary, and the gravitational force exerted by the fluid of density . In the case when the second primary is not an oblate spheroid (i.e., ), the equations are the same as those of Hallan and Mangang [4].

3. Position of Equilibrium Points

The equilibrium points are the solutions of the equations: That is,

3.1. Equilibrium Points Near the Centre of the First Primary

The positions of the equilibrium points near the first primary are the solutions of (3.2) when , , , , and . The coordinate of the equilibrium points are then the roots of the equation:

We first determine the roots of (3.5) in the absence of oblateness, that is, the case when the primaries are spherical. In this case, the roots are [4] The term which appears in (3.2) and is due to the fluid mass affects these roots. Therefore, these roots will be real if the discriminant is nonnegative, that is if When , both roots are greater than unity and we reject them because they lie outside the first primary. Now, if , we have and . Further, we see that when . Thus, in the case when , the point lies within the first primary if . When , ; is a root of (3.5). Hence for , the only root is which lies outside the first primary and we neglect it. Hence, for , , is always a root of (3.5) and is also a root provided , .

Now, we find the roots of (3.5) when oblateness of both primaries is considered (i.e., , ).

Let the roots be such that Putting these values in (3.5), multiplying throughout by , expanding and neglecting second and higher powers of , as they are very small quantities, we have Similarly, putting in (3.5) and then simplifying it, we get Multiplying (3.8) by , simplifying and then using it in (3.10), we get A substitution of (3.11) in the second equation of (3.8) at once gives the position of the other equilibrium point near the center of the first primary.

3.1.1. Positions of Circular Points

The positions of the circular points are sought using the first two equations of system (3.1) with the conditions , , ; that is, they are the solutions of Solving the above equations and knowing that , we get We let Substituting (3.15) in (3.14), and simplifying, we get Now, we let Substituting (3.16) in (3.15), neglecting second and higher powers of , we get Therefore, (3.17) is now expressed as A substitution of (3.14) in (3.13) yields Therefore, when , the points on the circle given by (3.15) with and lying within the first primary are also equilibrium points. The general coordinates of these circular points are given by , where is a parameter. When , the circular points coalesce to those lying on the line joining the primaries.

3.1.2. Positions of Out-of-Plane Equilibrium Points

The out-of-plane points have no analogy in the classical restricted three-body problem. However the investigation concerning these points in the photogravitational restricted three-body problem was first carried out by Radzievskii [10]. Afterwards, other researchers, for instance Douskos and Markellos [8], Singh and Leke [11], and so forth, have worked on the out-of-plane points. In this section, we locate these points for our study, as it has remained an open problem to date.

The positions of the out-of-plane equilibrium points of the Robe’s problem with oblate primaries are the solutions of the first and last equations of (3.1) with , ; that is, From (3.22), since , we have Let Then, (3.23) and (3.21) may be written respectively: Now, from first equations (3.25), we get The use of (3.24) in second equation of (3.25) yields We use the software package Mathematica (Wolfram 2004) to compute the coordinates of the out-of-plane equilibrium points denoted by and starting with the initial values and in the case where we have kept up to first order terms in both the numerator and the denominator; we then get The location of the out-of-plane equilibrium points can be obtained by solving numerically equations (3.26) and (3.27) using (3.24).

Now, from the expression for the density parameter We assume that , then . In the case when the density parameter is positive, we have Hence, numerically we choose Then, Now, we perform a numerical exploration of computing the out-of-plane points in the case of the Earth-Moon system. To do this, we arbitrarily choose values for the   . We found that when and , the positions of the out-of-plane points : The abscissae of the out-of-plane point is outside the possible region of motion of the infinitesimal mass and so we neglect it. However, in the case when the    are chosen such that , the point , and are, respectively, and lies within the fluid.

4. Linear Stability of the Equilibrium Points

In order to study the linear stability of any equilibrium point of an infinitesimal body, we displace it to the position such that where are small displacements, and then linearize equation (2.3) to obtain the equations: where the partial derivatives are evaluated at the equilibrium points.

4.1. Equilibrium Points Near the Center of the First Primary

In order to consider the motion near any equilibrium point in the -plane, we let solutions of the first two equations of (4.2) be where , , and are constants.

Taking first and second derivatives of the above, substituting them into the first two equations of system (4.2) and has a non-zero solution when Expanding the determinant, we have Equation (4.5) is the characteristic equation corresponding to the variational equations (4.2) in the case when motion is considered in the -plane.

Now, the values of the second-order partial derivatives of the equilibrium point , where stands for the first equilibrium and for the second one, are given as Substituting these in (4.2), we at once have the variational equations: where the partial derivatives have been computed at each equilibrium point .

Now, (4.8) is independent of (4.7), the solution being a periodic function is bounded and therefore, the motion of the infinitesimal body in the direction is stable.

Now, the characteristic equation of the equilibrium points corresponding to the system (4.7) is where These equations have been obtained using binomial expansion and ignoring terms with second and higher power in , and their product.

Now, let and be the roots of (4.9), then, the equilibrium point is stable if both the roots are real and negative. This means that their sum must be negative and their product must be positive. Hence, the points will be stable if the following two conditions hold: Now, in the case of the first equilibrium point , if we suppose in (4.10) that then, since , , and when , we have .

Similarly, in (4.11), if we suppose then .

For the case , we will have when which is not possible, hence .

In the case , , and .

Also, if , we see that and . Thus, for the case , the equilibrium point is stable. For , the equilibrium point is stable if the condition (4.12) holds. When , the equilibrium point is unstable.

Next, for the other equilibrium point positioned at , when , then since and the equilibrium point is stable if the conditions (4.12) and (4.13) are satisfied. If then ; it makes , . Therefore, when , both the conditions (4.12) and (4.13) are fulfilled and the equilibrium point is stable.

4.2. Circular Points

At circular points , the values of the second partial derivatives with the use of (3.14) and neglecting the product are Substituting these values in the variational equations (4.2), we get Equation (4.16) is independent of (4.15), it shows that the motion of the infinitesimal mass along the -direction is stable.

Now, a substitution of these partial derivatives in the characteristic equation (4.5) yields Let in (4.17) then, we have Hence, either which implies that Therefore, (4.20) gives the roots of the characteristic equation (4.17). Hence we conclude that the circular points are unstable due to the presence of multiple roots.

4.3. Out-of-Plane Points

To determine the stability of the out-of-plane equilibrium points, we consider the following partial derivatives: Since , therefore the partial derivatives to be computed at the out-of-plane equilibrium points are Now, we let, Using (4.23), the variational equation can be recast in the form: In order to consider the motion of the out-of-plane points, we let solution of the system (4.24) be where , and are constants. , , and are the small displacements in the coordinates of the infinitesimal body.

Now, the characteristic equation corresponding to the variational equations (4.24) in the case of the out-of-plane point may be expressed as where the coefficients of the characteristic equation (4.26) are such that where .

These computations have been done using the software package Mathematica.

For the stability analysis of the out-of-plane equilibrium point, we compute numerically the partial derivatives calculated at the out-of-plane points with the use of (3.28) and the following numerical values:

Now, substituting the above values in the characteristic equation (4.27), we get Its roots are: The positive root and the positive real part of the complex roots induce instability at the out-of-plane point. Hence, the motion of the infinitesimal mass around the out-of-plane equilibrium points is unstable for the specific numerical example given here. However, fuller discussion of their stability remains a theme for future research.

5. Discussion

The equation of motion (2.3) is different from those of Hallan and Mangang [4] due to oblateness of the second primary. If we assume that the second primary is not oblate (i.e., ), then these equations will fully coincide with those of Hallan and Mangang [4].

Equation (3.9) gives the equilibrium position of the point near the center of the first primary and fully coincides with that of Hallan and Mangang [4]. It shows that the position of this equilibrium point does not depend on oblateness of the second primary, while the other equilibrium point given by (3.11) is different from that of Hallan and Mangang [4] due to the appearance of oblateness of the second primary. When , points on the circle , lying within the first primary are also equilibrium points. These points are affected by oblateness of both primaries. Equations (3.28) give the positions of the out-of-plane points when only linear terms in oblateness of the second primary are retained. We have been able to show that the oblateness of the primaries allows the existence of the out-of-plane equilibrium points in the -plane within the first primary. These points have no analogy in the previous studies of the Robe’s restricted three-body problem.

The linear stability analysis of the equilibrium solutions of the problem is investigated with the help of characteristic roots. The characteristic equation (4.9) in the case of the equilibrium point near the center of the first primary is the same as that of Hallan and Mangang [4], while that of the other point near the center differs from that of Hallan and Mangang [4] due to oblateness of the second primary. The characteristic equation of the circular case (4.17) also differs from that of Hallan and Mangang [4] due to oblateness of the second primary. The stability in the first approximation of this configuration shows that points near the centre of the first primary are conditionally stable; the circular points are unstable. This confirms the earlier results of Hallan and Rana [3], Hallan and Mangang [4]. A numerical exploration shows that the out-of-plane equilibrium points are also unstable. This outcome validates the earlier results of Douskos and Markellos [8] and Singh and Leke [11] that the points are unstable.

6. Conclusion

We have derived the equations of motion and established the positions of the equilibrium points of the infinitesimal body in the Robe’s [1] restricted three-body problem with oblateness. The term “oblateness” is used in the sense that both primaries are considered as oblate spheroids under the effects of the full buoyancy force exerted by the fluid on the infinitesimal mass.

We have obtained one equilibrium point near the centre of the first primary which will be on the left or right of the centre of the first primary accordingly as . This point is the same as that of Hallan and Mangang [4]. In addition to this, another equilibrium point is found within the first primary on the line joining the center of the primaries when and . When , points on the circle , lying within the first primary are also equilibrium points. We call them circular points. Finally, we have been able to show that the oblateness of the primaries allows the existence of the out-of-plane equilibrium points in the -plane within the first primary.

The result of this paper can be summarized as follows. The restricted three-body problem under the framework of the Robe’s [1] problem with oblate primaries has the equilibrium points of the type: points near the center of the first primary, points on the circle (circular points), and two out-of-plane points . It is seen that points near the first primary are conditionally stable, the circular points are unstable, while the out-of-plane equilibrium points are unstable for the specific numerical example given here. The effect of drag forces as considered by Giordano et al. [2] under the present context, particularly as regards the analysis of the properties of the equilibrium points located inside the first primary, will be interesting.