Investigation of the effect of albedo and oblateness on the circular restricted four variable bodies problem

Open access

Abstract

The present paper investigates the motion of the variable infinitesimal body in circular restricted four variable bodies problem. We have constructed the equations of motion of the infinitesimal variable mass under the effect of source of radiation pressure due to which albedo effects are produced by another two primaries and one primary is considered as an oblate body which is placed at the triangular equilibrium point of the classical restricted three-body problem and also the variation of Jacobi Integral constant has been determined. We have studied numerically the equilibrium points, Poincaré surface of sections and basins of attraction in five cases (i. Third primary is placed at one of the triangular equilibrium points of the classical restricted three-body problem, ii. Variation of masses, iii. Solar radiation pressure, iv. Albedo effect, v. Oblateness effect.) by using Mathematica software. Finally, we have examined the stability of the equilibrium points and found that all the equilibrium points are unstable.

Abstract

The present paper investigates the motion of the variable infinitesimal body in circular restricted four variable bodies problem. We have constructed the equations of motion of the infinitesimal variable mass under the effect of source of radiation pressure due to which albedo effects are produced by another two primaries and one primary is considered as an oblate body which is placed at the triangular equilibrium point of the classical restricted three-body problem and also the variation of Jacobi Integral constant has been determined. We have studied numerically the equilibrium points, Poincaré surface of sections and basins of attraction in five cases (i. Third primary is placed at one of the triangular equilibrium points of the classical restricted three-body problem, ii. Variation of masses, iii. Solar radiation pressure, iv. Albedo effect, v. Oblateness effect.) by using Mathematica software. Finally, we have examined the stability of the equilibrium points and found that all the equilibrium points are unstable.

1 Introduction

Few-body problem attract many scientists for a long time in celestial mechanics and dynamical astronomy. The restricted three-body problem and four-body problem with many perturbations like different shapes of the primaries, resonance, variable mass of the primaries as well as infinitesimal body, Coriolis and centrifugal forces, Pointing-Robertson drag, solar radiation pressure and albedo effects etc., have been studied by many scientists. Simó [36] investigated the linear stability of the relative Lagrangian solutions in the four-body problem. Hadjidemetriou [17]studied the periodic orbits of fourth body in the restricted four-body problem with respect to rotating frame. Michalodmitrakis [27] generalized the restricted three-body problem to the restricted four-body problem and studied about the equilibrium points, regions of possible motion and periodic orbits. Kalvouridis et al. in the series of three papers [20,21 and 22] investigated the equilibrium points and stability in the restricted four-body problem under the effect of oblateness and radiation pressure. And also performed the zero-velocity surface and curves. Baltagiannis et al. [14] shown that equilibrium points depend on the mass of the primaries in the restricted four-body problem. Papadauris et al. [31,32] investigated the existence, locations, stability and periodic orbits of the equilibrium points in and out of the orbital plane in the photo-gravitational circular restricted four-body problem. Ansari [1] studied the periodic orbits in the restricted four-body problem around lagrangian points in three cases. In the first case, he has considered all three primaries as spherical in shape. In the second case, he has taken one of the three primaries as an oblate body. And in the third and last case, he has taken two of the primaries as oblate body and all the three primaries are source of radiation pressure. Falaye [16] investigated the stability of the equilibrium points in the restricted four-body problem under the effects of oblateness and solar radiation pressure and found that these equilibrium points are unstable. Arribas et al. [22] investigated the equilibria of the symmetric collinear restricted four-body problem where primaries are placed in a collinear central configuration with both masses and radiation pressure of the peripheral bodies are equal. Papadakis [30] performed the 21 families of simple 3D symmetric periodic orbits as well as the typical orbits of all symmetry type 3D orbits in the circular restricted four-body problem. After examined the stability, he illustrated the characteristic curves and stability diagrams of families of 3D periodic orbits. Asique et al. [10][11][12][13] studied the restricted four-body problem with different shapes of the primaries with solar radiation pressure. They have placed one of the primaries at the lagrangian points of the classical restricted three-body problem. They have illustrated the equilibrium points and zero-velocity curves for these models. Singh et al. [45,46] investigated in and out of plane equilibrium points in the circular restricted four-body problem with the effect of solar radiation pressure.

On the other-hand many scientists have studied about the albedo on these models. Anselmo et al. [8] performed the periodic perturbations of the satellite is the radiation pressure due to the sun-light reflected by the Earth. Rocco [34] evaluated the terrestrial albedo by using earth albedo model and orbital dynamics model and also calculated the irradiance at the satellite with radiant flux from each cells. Idrisi [18] investigated the existence and stability of the circular restricted three-body problem under the effect of albedo when smaller primary is an ellipsoid. Many scientists have investigated on these models with variable masses as Jeans [19], Meshcherskii [26], Shrivastava et al. [35], Lichtenegger [24], Singh et al. [37][38][39][40][41], Lukyanov [25], Zhang et al. [44], Abouelmagd and Mostafa [2], Mittal et al. [28], Ansari [3][4][5][6], etc. And also many scientists have explained the basins of attraction in these models as Douskos [15], Kumari and Kushvah [23], Paricio [33], Ansari [5,6], Zotos [45][46][47][48], etc.

This paper deals the motion of the infinitesimal body in the circular restricted four-body problem in which the masses of the primaries as well as the mass of the infinitesimal body vary and also one of the primaries is taken as source of radiation pressure due to which albedo effects are produced and another primary is considered as oblate body. We have studied this problem in various sections. In the statement of the problem and equations of motion section, we have constructed the equations of motion of the infinitesimal variable mass under the effect of these perturbations where the third primary is fixed at the triangular equilibrium point of the classical restricted three-body problem and also the variation of Jacobi Integral constant has been determined. In the computational work section, we have plotted numerically the equilibrium points, Poincaré surface of sections and basins of attraction in five cases by using Mathematica software. In the stability section, we have examined the stability of the equilibrium points under the effect of the perturbations. Finally, we have concluded the problem.

2 Statement of the problem and equations of motion

Let m1(t), m2(t), m3(t) and m(t) be four variable masses (m1(t) > m2(t) >> m3(t)). The configuration of the m1(t), m2(t) and m(t) is taken as the restricted three-body problem. The third primary m3(t) is considered as so small that it is not influencing the motion of the primaries m1(t) and m2(t). Thus the center of the rotation of the system remains the same as for the bodies m1(t) and m2(t). The body m3(t) is placed at one of the triangular equilibrium points of the R3BP. m(t) is moving in the space under the influence of these three primaries but not influencing them. In this way these four variable bodies formed the restricted circular four-body problem.

Here we have considered m3(t) as oblate body with oblateness factor σ (the shape of the oblate body remain same though the mass is variable) and m1(t) is a source of radiations(Fp) due to these radiations, m2(t) and m3(t) produce albedo FA and Fa respectively (i.e. Albedo = (radiation reflected back into the space)/(incident radiation)). The infinitesimal body m(t) is moving in the space under the gravitational forces F1, F2 and F3 of the primaries respectively. The total force on infinitesimal body will be F = F1(1−p1) + F2(1−p2) + F3(1−p3). Where 0≤p1=Fp/F1 < 1,0 ≤ p2 = FA/F2 < p1, 0 ≤ p3 = Fa/F3 < p1. The line joining m1(t) and m2(t) is taken as x-axis and are moving around their common center of mass which is taken as origin O and the line perpendicular to the x-axis and passing through the origin is considered as y-axis, the line passing through the origin and perpendicular to the plane of motion of the primaries is taken as z-axis. Let us consider the synodic coordinate system, initially coincide with the inertial coordinate system, and revolving with angular velocity ω(t) about z-axis. Let the coordinates of m1(t), m2(t), m3(t) and m(t) in the rotating frame be (x1, 0, 0), (x2, 0, 0), (x3, y3, 0) and (x, y, z) respectively (Fig. 1). Following the procedure given in Ansari [7], we can write the equations of motion of the infinitesimal variable mass in the rotating coordinate system when the variation of mass is non-isotropic with zero momentum as

Fig. 1

Download Figure

Fig. 1

Configuration of the problem in CRFVBP with the effect of oblateness and Albedo.

Citation: Applied Mathematics and Nonlinear Sciences 2, 2; 10.21042/AMNS.2017.2.00044

{m˙(t)m(t)(x˙ω(t)y)+(x¨ω˙(t)y2ω(t)y˙ω2(t)x)=μ1(t)(xx1)(1p1)r13μ2(t)(xx2)(1p2)r23μ3(t)(xx3)(1p3)r333μ3(t)(xx3)(1p3)σ2r35,m˙(t)m(t)(y˙+ω(t)x)+(y¨+ω˙(t)x+2ω(t)x˙ω2(t)y)=μ1(t)y(1p1)r13μ2(t)y(1p2)r23μ3(t)y(1p3)r333μ3(t)(yy3)(1p3)σ2r35,m˙(t)m(t)z˙+z¨=μ1(t)z(1p1)r13μ2(t)z(1p2)r23μ3(t)z(1p3)r333μ3(t)z(1p3)σ2r35.

where, ri2=(xxi)2+(yyi)2+z2 are the distances from the primaries to the infinitesimal body in the rotating coordinate system respectively, μi(t) = G mi(t), (i = 1, 2, 3).

Using Meshcherskii transformation

x=ξR(t),y=ηR(t),z=ζR(t),dtdτ=R2(t),ri=ρiR(t),ω(t)=ω0R2(t),xi=ξiR(t),yi=ηiR(t),μ(t)=μ1(t)+μ2(t)=μ0R(t),μi(t)=μi0R(t),m(t)=m0R(t),R(t)=αt2+2βt+γ,

where α, β, γ, μ0, μi0, m0 are constants for i = 1, 2, 3.

The system (1) becomes

{ξ2ω0η(αt+β)ξ=Wξ,η+2ω0ξ(αt+β)η=Wη,ζ(αt+β)ζ=Wζ.

where,

W=12((αt+β)2+ω02(αγβ2))(ξ2+η2)+12((αt+β)2(αγβ2))ζ2(αt+β)ξη+μ10(1p1)ρ1+μ20(1p2)ρ2+μ30(1p3)ρ3+μ30(1p3)σ2ρ33,ρi2=(ξξi)2+(ηηi)2+ζ2.

Prime (′) is the differentiation w.r.to τ. Taking unit of mass, distance and time at initial time t0 such that

μ0=1,ξ1+ξ2=1,G=1,ω0=1,αt0+β=α1(constant).

Introducing the new mass parameter as μ10 = 1−υ, μ20 = υ, μ30 = α2 υ, where α2 << 1 and υ is the ratio of the mass of the primary m2 to the total mass of the primaries m1 and m2. Finally, the system (2) becomes

{ξ2ηα1ξ=Uξ,η+2ξα1η=Uη,ζα1ζ=Uζ.

where, U=12(α12+k)(ξ2+η2+ζ2)12ζ2α1ξη+(1υ)(1p1)ρ1+υ(1p2)ρ2+α2υ(1p3)ρ3+α2υ(1p3)σ2ρ33, ρi2=(ξξi)2+(ηηi)2+ζ2,αγβ2=1k and (ξ1, η1) = (υ, 0),(ξ2, η2) = (1−υ, 0), (ξ3,η3)=(υ12,32).

If there are constant masses then there is a constant motion (i.e Moulton [29]), the Jacobi Integral Constant defined as

JIC=2U2((ξ)2+(η)2+(ζ)2),

Multiplying in the first equation of (3) by ξ′, in the second equation of (3) by η′ and in the third equation of (3) by ζ′ and add and using equation (4), we get the variation of the Jacobi Integral Constant as

dJICdτ=2α1((ξ)2+(η)2+(ζ)2),

Where JIC is the Jacobi Integral Constant.

3 Computational work

In this section, we have drawn the locations of equilibrium points, the Poincaré surfaces of section and the basins of attraction for five different cases by using Mathematica software (in all the calculations υ = 0.019, α2 = 0.01):

  • Third primary is placed at one of the Lagrangian points of the classical restricted three-body problem

    (i.e. k = 1, α1 = 0, p1 = 0, p2 = 0, p3 = 0, σ = 0),

  • Variation of masses

    (i.e. k = 0.4, α1 = 0.2, p1 = 0, p2 = 0, p3 = 0, σ = 0),

  • Solar radiation pressure

    (i.e. k = 0.4, α1 = 0.2, p1 = 0.5, p2 = 0, p3 = 0, σ = 0),

  • Albedo effect

    (i.e. k = 0.4, α1 = 0.2, p1 = 0.5, p2 = 0.3, p3 = 0.2, σ = 0),

  • Oblateness effect

    (i.e. k = 0.4, α1 = 0.2, p1 = 0.5, p2 = 0.3, p3 = 0.2, σ = 0.01).

3.1 Equilibrium points during in-plane and out of plane motions

The solutions of Uξ = 0, Uη = 0, and Uζ = 0, will be the location of equilibrium points but the solutions of these equations represent the location of equilibrium points during in-plane motions when (ξ≠0, η≠0, ζ = 0), (Fig. 2) and represent the location of equilibrium points of the out of plane when (ξ≠0, η = 0, ζ≠0), (Fig. 3) and (ξ = 0, η≠0, ζ≠0), (Fig. 4). At the equilibrium points, all the derivatives of the co-ordinates with respect to the time will be zero. Therefore the infinitesimal body will stay at these points with zero velocity. Hence the singularities are known as stationary points. Where

Fig. 2

Download Figure

Fig. 2

The locations of lagrangian points during in-plane motion (i.e. ξ≠0, η≠0, ζ = 0) in five cases: (a). Third primary is placed at one of the lagrangian points of the classical circular restricted three-body problem (points with blue color), (b). Variation of masses (points with green color), (c). Solar radiation pressure (points with red color), (d). Albedo effect (points with black color), (e). Oblateness effect (points with magenta color).

Citation: Applied Mathematics and Nonlinear Sciences 2, 2; 10.21042/AMNS.2017.2.00044

Uξ=(α12+k)ξα1η(1υ)(ξξ1)(1p1)ρ13υ(ξξ2)(1p2)ρ23υ3(ξξ3)(1p3)ρ333α2υ(ξξ3)(1p3)σ2ρ35,

Uη=(α12+k)ηα1ξ(1p1)(1υ)ηρ13(1p2)υηρ23υ3(ηη3)(1p3)ρ333α2υη(1p3)σ2ρ35,

Uζ=(α12+k1)ζ(1p1)(1υ)ζρ13(1p2)υζρ23υ3ζ(1p3)ρ333α2υζ(1p3)σ2ρ35.

3.1.1 During in-plane motion (i.e. ξ ≠ 0, η ≠ 0, ζ = 0) locations of lagrangian points

During in-plane motion, we have plotted graphs for the locations of the lagrangian points in five cases. We found six lagrangian points in which three points (L1, L2, L3) are collinear and rest three points (L4, L5, L6) are non-collinear when third primary is placed at one of the lagrangian points of the classical restricted circular three-body problem (Fig. 2(a) points with blue color). It is observed that the lagrangian points (L1, L4, L5, L6) and (L2, L3) lie left side and right side of the origin respectively. In the variable mass case, we found six lagrangian points (Fig. 2(b) points with green color) in which points (L2, L3, L6) and (L1, L4, L5) lie left side and right side of the origin respectively. Here (L3, L4) and (L5, L6) are symmetrical and rest points are non-symmetrical. In the rest three cases, we found eight lagrangian points with slight different locations (Fig. 2(c) points with red color, Fig. 2(d) points with black color, Fig. 2(e) points with magenta color). In all the figures black stars denote the locations of the primaries respectively.

3.1.2 During out of plane (i.e. ξ ≠ 0, η = 0, ζ ≠ 0 and ξ = 0, η ≠ 0, ζ ≠ 0) the locations of lagrangian points

During out of plane (i.e. ξ ≠ 0, η = 0, ζ ≠ 0), we found three lagrangian points on the ξ-axis in all the five cases (Fig. 3 (a) points in blue color, Fig. 3(b) points in green color, Fig. 3(c) points with red color, Fig. 3(d) points with black color, Fig. 3(e) points with magenta color). In all the cases, the points L1 and (L2, L3) lie left and right side of the origin with slight different locations respectively. The black stars denote the locations of the primaries m1 and m2.

Fig. 3

Download Figure

Fig. 3

The locations of lagrangian points during out of plane motion (i.e. ξ ≠ 0, η = 0, ζ ≠ 0) in five cases: (a). Third primary is placed at one of the lagrangian points of the classical circular restricted three-body problem (points with blue color), (b). Variation of masses (points with green color), (c). Solar radiation pressure (points with red color), (d). Albedo effect (points with black color), (e). Oblateness effect (points with magenta color).

Citation: Applied Mathematics and Nonlinear Sciences 2, 2; 10.21042/AMNS.2017.2.00044

On the other hand, during the out of plane (i.e. ξ = 0, η ≠ 0, ζ ≠ 0), we found three lagrangian points (Fig. 4) in all five cases. The points L1 and L3 lie left and right side of the origin respectively but the point L2 lies at the origin. It is also observed from the Figure 4 that lagrangian points are moving away from the origin from the first case to the variable mass case and then toward the origin in the other cases.

Fig. 4

Download Figure

Fig. 4

The locations of lagrangian points during out of plane motion (i.e.ξ = 0, η ≠ 0, ζ ≠ 0) in five cases: (a). Third primary is placed at one of the lagrangian points of the classical circular restricted three-body problem (points with blue color), (b). Variation of masses (points with green color), (c). Solar radiation pressure (points with red color), (d). Albedo effect (points with black color), (e). Oblateness effect (points with magenta color).

Citation: Applied Mathematics and Nonlinear Sciences 2, 2; 10.21042/AMNS.2017.2.00044

3.2 Poincaré Surface of Section

We also have drawn the Poincaré surface of sections for five cases in both the (ξξ′)−plane (Fig.5 (i)) and the (ηη′)−plane (Fig.5 (ii), Fig. 5 (iii)). It is observed that in the ξξ′-plane, the surfaces are shrinking and in the ηη′-plane, the surfaces first expanding and then shrinking.

Fig. 5

Download Figure

Fig. 5

Poincaré surface of sections for five cases: (a). Third primary is placed at one of the lagrangian points of the classical circular restricted three-body problem (with blue color), (b). Variation of masses (with green color), (c). Solar radiation pressure (with red color), (d). Albedo effect (with black color), (e). Oblateness effect (with magenta color). (i) Poincaré surface of sections in ξξ′-plane, (ii) Poincaré surface of sections in ηη′-plane, (iii) Zoomed part of figure (ii) near the origin.

Citation: Applied Mathematics and Nonlinear Sciences 2, 2; 10.21042/AMNS.2017.2.00044

3.3 Basins of attraction

In this section, we have drawn the basins of attraction for the circular restricted four variable bodies problem in which we have taken one primary as solar radiation pressure due to which other two primaries produced albedo and also one of the primaries as an oblate body, by using Newton-Raphson iterative method for the five cases (i. Third primary is placed at one of the lagrangian points of the restricted three-body problem, ii. Variation of mass case, iii. Solar radiation pressure effect, iv. Albedo effect, v. Oblateness effect). The algorithm of our problem is given by

{ξn+1=ξn(UξUηηUηUξηUξξUηηUξηUηξ)(ξn,ηn),ηn+1=ηn(UηUξξUξUηξUξξUηηUξηUηξ)(ξn,ηn).

Where ξn, ηn are the values of ξ and η coordinates of the nth step of the Newton-Raphson iterative process. If the initial point converges rapidly to one of the lagrangian points then this point (ξ, η) is a member of the basin of attraction of the root. This process stops when the successive approximation converges to an attractor. We used color code for the classification of the lagrangian points on the (ξ, η)−plane. In the first case (Fig. 6), L1, L2, L3 represent cyan color regions, L4, L5, L6 represents light blue color regions. The basins of attraction corresponding to the lagrangian points L1, L2, L3, L4, L5, L6 extend to infinity. In the variable mass case (Fig. 7), L1 represents cyan color region, L2, L3, L5 represent light blue color regions, L4 and L6 represent light green color regions. The basins of attraction corresponding to the lagrangian point L1 cover finite area but corresponding to the lagrangian points L2, L3, L4, L5, L6 extend to infinity. While in the rest three cases (Fig. 8, Fig. 9, Fig. 10), L1 represents cyan color region, L2, L3 represent light blue, L5 represents blue color region and L4, L6, L7, L8 represent light green color regions. The basins of attraction corresponding to the all lagrangian points L2, L3, L4, L5, L6, L7, L8 extend to infinity except L1 which covers finite area. In this way a complete view of the basin structures created by the attractors. We can observe in detail from the zoomed part of all the figures in Fig. 6(b), Fig. 7(b), Fig. 8(b), Fig. 9(b), Fig. 10(b). The black points and black stars denote the location of the lagrangian points and the primaries respectively.

Fig. 6

Download Figure

Fig. 6

(a): The basin of attraction for the case when third primary is placed at one of the lagrangian points of the classical circular restricted three- body problem. (b): Zoomed image of (a) near the lagrangian configuration.

Citation: Applied Mathematics and Nonlinear Sciences 2, 2; 10.21042/AMNS.2017.2.00044

4 Stability of the equilibrium points

We can examine the stability of the equilibrium points under the effect of Albedo and oblateness when all the masses are varying, by taking ξ = ξ0 + ξd, η = η0+ ηd, ζ = ζ0 + ζd in system (3), we get

{ξd2ηdα1ξd=ξd Uξξ0+ηd Uξη0+ζd Uξζ0ηd+2ξdα1ηd=ξd Uηξ0+ηd Uηη0+ζd Uηζ0ζdα1ζd=ξd Uζξ0+ηd Uζη0+ζd Uζζ0

Where ξd, ηd and ζd are the small displacements of the infinitesimal body from the equilibrium point. The superscript zero denotes the value at the equilibrium point.

To solve system (10), let ξd = C1eλτ, ηd = C2eλτ, ζd = C3eλτ, where C1, C2 and C3 are constant parameters. Substituting these values in system (10) and rearranging, we get:

{C1(λ2α1λUξξ0)C2(2λ+Uξη0)C3Uξζ0=0,C1(2λUηξ0)+C2(λ2α1λUηη0)C3Uηζ0=0,C1Uζξ0C2Uζη0+C3(λ2α1λUζζ0)=0.

The system (11), will have a non-trivial solution for C1, C2 and C3 if

|λ2α1λUξξ0(2λ+Uξη0)Uξζ02λUηξ0λ2α1λUηη0Uηζ0Uζξ0Uζη0λ2α1λUζζ0|=0,

which is equivalent to

λ63α1λ5+λ4(4+3α12Uξξ0Uηη0Uζζ0)+α1λ3(4α12+2Uξξ0+2Uηη0+2Uζζ0)+λ2((Uξη0)2(Uξζ0)2+Uξξ0Uηη0(Uηζ0)24Uζζ0+Uξξ0Uζζ0+Uζζ0Uηη0α12Uξξ0α12Uηη0α12Uζζ0)+α1λ((Uξη0)2+(Uξζ0)2Uξξ0Uηη0+(Uηζ0)2Uξξ0Uζζ0Uηη0Uζζ0)+(Uξζ0)2Uηη02Uξη0Uξζ0Uηζ0+Uξξ0(Uηζ0)2+(Uξη0)2Uζζ0Uξξ0Uηη0Uζζ0=0,

From the solution of the equation (12) we found that λ has complex values in all the cases and at least one of them has positive real value. Hence all the equilibrium points are unstable.

5 Conclusions

In this paper, we have investigated the effect of Albedo and oblateness of the primary in the circular restricted four-body problem with variable masses. We have determined the equations of motion, when the masses of the primaries as well as the infinitesimal body vary, which are different from the classical case by the variation parameters α1, k, oblateness factor σ and the radiations effect p1, p2, p3. Further the expression for the variation of Jacobi integral constant have been evaluated which is also depending on the variation parameter α1. We have plotted equilibrium points, Poincaré surface of sections and basins of attraction in five cases (i. Case when third primary is placed at one of the lagrangian points of the classical circular restricted three-body problem, ii. Variation of mass, iii. Solar radiation pressure effect, iv. Albedo effect, v. Oblateness effect) by using Mathematica software. The equilibrium points during in-plane motion (Fig. 2), we found six equilibrium points in which three are collinear and three are non-collinear equilibrium points in the first case, in the variable mass case, we got one is collinear and rest five are non-collinear, while in the cases solar radiation pressure, albedo and oblateness, we found eight equilibrium points. During out-of-plane motions (i.e. ξ ≠ 0, η = 0, ζ ≠ 0)(Fig. 3) and (i.e. ξ = 0, η ≠ 0, ζ ≠ 0)(Fig. 4), we found three equilibrium points in all cases. The poincaré surface of sections have been determined in two phase spaces (ξξFig. 5 (i) and ηηFig. 5 (ii)). In ξξ′-plane, the poincaré surface of sections are shrinking and ηη′-plane, the poincaré surface of sections are first expanding and then shrinking . The Newton-Raphson basins of attraction have been studied in all five cases (Fig. 6, Fig. 7, Fig. 8, Fig. 9, Fig. 10). We used color code for the classification of the equilibrium points on the (ξη)−plane. These points will be clearly visible in the zoomed part of all the figures near the lagrangian configuration (ie. Fig. 6 b, Fig. 7 b, Fig. 8 b, Fig. 9 b, Fig. 10 b). Finally, we have examined the stability of the equilibrium points in the circular restricted four variable bodies problem with oblateness and albedo effect and found that all the equilibrium points are unstable.

Fig. 7

Download Figure

Fig. 7

(a): The basin of attraction for the variable mass case. (b): Zoomed image of (a) near the lagrangian configuration.

Citation: Applied Mathematics and Nonlinear Sciences 2, 2; 10.21042/AMNS.2017.2.00044

Fig. 8

Download Figure

Fig. 8

(a): The basin of attraction for the solar radiation pressure case. (b): Zoomed image of (a) near the lagrangian configuration.

Citation: Applied Mathematics and Nonlinear Sciences 2, 2; 10.21042/AMNS.2017.2.00044

Fig. 9

Download Figure

Fig. 9

(a): The basin of attraction for the Albedo case. (b): Zoomed image of (a) near the lagrangian configuration.

Citation: Applied Mathematics and Nonlinear Sciences 2, 2; 10.21042/AMNS.2017.2.00044

Fig. 10

Download Figure

Fig. 10

(a): The basin of attraction for the oblateness case. (b): Zoomed image of (a) near the lagrangian configuration.

Citation: Applied Mathematics and Nonlinear Sciences 2, 2; 10.21042/AMNS.2017.2.00044

Communicated by Elbaz I. Abouelmagd

Acknowledgement

We are thankful to the Deanship of Scientific Research and the Department of Mathematics, College of Science Al-Zulfi, Majmaah University, Kingdom of Saudi Arabia for providing all the research facilities. We are also grateful to the referee for his valuable suggestions.

References

  • [1]

    Ansari, A. A, (2014),Periodic orbits around lagrangian points of the circular restricted four-body problem. Invertis Journal of Science and Technology. 7(1), 29-38.

  • [2]

    Abouelmagd, E. I., Mostafa, A.,(2015), Out of plane equilibrium points locations and the forbidden movement regions in the restricted three-body problem with variable mass. Astrophys. Space Sci 357, 58,

  • [3]

    Ansari, A. A., (2016(a)), Stability of the equilibrium points in the photogravitational circular restricted four body problem with the effect of perturbations and variable mass.Science International (Lahore). 28, 859- 866.

  • [4]

    Ansari, A. A., (2016(b)), Stability of the equilibrium points in the circular restricted four body problem with oblate primary and variable mass. International Journal of Advanced Astronomy. 4(1), 14-19.

  • [5]

    Ansari, A. A.(2016 (c)), The Photogravitational Circular Restricted Four-body Problem with Variable Masses. Journal of Engineering and Applied Sciences. 3(2), 30-38.

  • [6]

    Ansari, A. A.,(2017(a)), The circular restricted four body problem with variable masses. Nonlinear Sci. Lett. A. 8(3), 303-312.

  • [7]

    Ansari, A. A.,(2017(b)), Effect of Albedo on the motion of the infinitesimal body in circular restricted three-body problem with variable masses. Italian J. of Pure and Applied Mathematics. 38, 581-600.

  • [8]

    Anselmo, L., Farinella, P., Milani, A., Nobili, A. M.,(1983), Effects of the earth-reflected sunlight on the orbit of the LAGEOS Satellite. Astron. Astrophys. 117, 3-8.

  • [9]

    Arribas, M., Abad, A., Elipe, A., Palacios, M.,(2016), Out-of-plane equilibria in the symmetric collinear restricted four-body problem with radiation pressure.Astrophys. Space Sci 361, 270(8). 10.1007/s10509-016-2858-1

  • [10]

    Asique, M. C., et al.,(2015(a)),On the R4BP when third primary is an oblate spheroid. Astrophys. Space Sci 357, 82(1).

  • [11]

    Asique, M. C., et al.,(2015(b)),On the photogravitational R4BP when the third primary is an oblate/prolate spheroid. Astrophys. Space Sci 360, 13(1),

  • [12]

    Asique, M. C., et al., (2016), On the photogravitational R4BP when the third primary is a triaxial rigid body. Astrophys. Space Sci 361, 379.

  • [13]

    Asique, M. C., et al.,(2017), On the R4BP when Third Primary is an Ellipsoid. Journal of Astronaut. Sci 64, 231-250,

  • [14]

    Baltagiannis, A., Papadakis, K. E.,(2011), Equilibrium points and their stability in the restricted four body problem. International Journal of Bifurcation and Chaos. 21(8), 2179-2193,

  • [15]

    Douskos, C. N.,(2010), Collinear equilibrium points of Hill's problem with radiation pressure and oblateness and their fractal basins of attraction. Astrophys. Space Sci 326, 263-271.

  • [16]

    Falaye, B. J.,(2015), Effect of Oblateness, Radiation and a Circular Cluster of Material Points on the Stability of Equilibrium Points in the Restricted Four-Body Problem. Few body system. 56, 29-40,

  • [17]

    Hadjidemetriou, J. D.,(1980), The Restricted Planetary 4-Body Problem.Celest. Mech. 21, 63-71.

  • [18]

    Idrisi, M. J.(2017), A study of Libration Points in Modified CR3BP under Albedo Effect when smaller Primary is an Ellipsoid. J. of Astronaut. Sci,

  • [19]

    Jeans, J. H., (1928), Astronomy and Cosmogony, Cambridge University Press, Cambridge.

  • [20]

    Kalvouridis, T. J., Mavraganis, A. G.,(1995), Equilibria and stability of the restricted photogravitational problem of 2 +  2 bodies. Astrophys. Space Sci 226(1), 137-148.

  • [21]

    Kalvouridis, T. J.,(1996),The oblate spheroids version of the restricted photogravitational 2 + 2 body problem. Astrophys. Space Sci 246(2), 219-227.

  • [22]

    Kalvouridis, T. J., Arribas, M., Elipe, A.,(2006),Dynamical properties of the restricted four-body problem with radiation pressure. Mechanics Research Communications. 33, 811-817.

  • [23]

    Kumari, R., Kushvah, B. S.,(2014), Stability regions of equilibrium points in the restricted four body problem with oblateness effects. Astrophys. Space Sci 349, 693-704.

  • [24]

    Lichtenegger, H.,(1984),The dynamics of bodies with variable masses. Celest. Mech 34, 357-368.

  • [25]

    Lukyanov, L. G.,(2009), On the restricted circular conservative three-body problem with variable masses. Astronomy Letters. 35(5), 349-359.

  • [26]

    Meshcherskii, I. V.,(1952), Works on the Mechanics of Bodies of Variable Mass, GITTL, Moscow.

  • [27]

    Michalodmitrakis, M.,(1981), The circular restricted four-body problem. Astrophys. Space Sci 75, 289-310.

  • [28]

    Mittal, A., et. al.,(2016), Stability of libration points in the restricted four-body problem with variable mass. Astrophys. Space sci 361, 329,

  • [29]

    Moulton, F. R., (1914), An introduction to celestial mechanic, Second ed. Dover, New York.

  • [30]

    Papadakis, K. E.,(2016), Families of three-dimensional periodic solutions in the circular restricted four-body Problem. Astrophys. Space Sci 361, 129.

  • [31]

    Papadouris, J. P., Papadakis, K. E.,(2013), Equilibrium points in the photogravitational restricted four-body Problem. Astrophys. Space Sci 344, 21-38.

  • [32]

    Papadouris, J. P., Papadakis, K. E.,(2014), Periodic solutions in the photogravitational restricted four-body problem. MNRAS. 442, 1628-1639.

  • [33]

    Paricio, L. J. H.,(2016), Bivariate Newton-Raphson method and toroidal attraction basins. Numerical Algo 71, 349-381.

  • [34]

    Rocco, E. M.,(2009), Evaluation of the terestrial Albedo-Applied to some scientific missions. Space Sci. Rev 151, 135-147,

  • [35]

    Shrivastava, A. K., Ishwar, B.,(1983),Equations of motion of the restricted problem of three bodies with variable mass. Celest. Mech 30, 323-328.

  • [36]

    Simo, C.,(1978), Relative equilibrium solutions in the four body problem. Celest. Mech 18, 165-184.

  • [37]

    Singh, J., Ishwar, B.,(1984), Effect of perturbations on the location of equilibrium points in the restricted problem of three bodies with variable mass. Celest. Mech 32, 297-305.

  • [38]

    Singh, J., Ishwar, B.,(1985), Effect of perturbations on the stability of triangular points in the restricted problem of three bodies with variable mass. Celest. Mech 35, 201-207.

  • [39]

    Singh, J.,(2003),Photogravitational restricted three body problems with variable mass. Indian Journal of Pure and Applied Math 32 (2), 335-341.

  • [40]

    Singh, J., Leke, O.,(2010), Stability of photogravitational restricted three body problem with variable mass. Astrophys. Space Sci 326 (2), 305-314.

  • [41]

    Singh, J., Leke, O.,(2013), Existence and Stability of Equilibrium Points in the Robe's Restricted Three-Body Problem with Variable Masses. International Journal of Astronomy and Astrophysics. 3, 113-122.

  • [42]

    Singh, J., Vincent, A. E.,(2016(a)),Equilibrium points in the restricted four-body problem with radiation pressure. Few-Body Syst 57, 83-91.

  • [43]

    Singh, J.,Vincent, A. E.,(2016(b)),Out-of-plane equilibrium points in the photogravitational restricted four-Body problem with oblateness. British Journal of Mathematics and Computer Science. 19(5), 1-15.

  • [44]

    Zhang, M. J., Zhao, C. Y., Xiong, Y. Q.,(2012), On the triangular libration points in photo-gravitational restricted three body problem with variable mass. Astrophys. Space Sci 337, 107-113.

  • [45]

    Zotos, E. E.,(2016(a)), Fractal basins of attraction in the planar circular restricted three body problem with oblateness and radiation pressure. Astrophys. Space Sci 181(17).

  • [46]

    Zotos,E. E., (2016(b)), Fractal basins boundaries and escape dynamics in a multi-well potential. Nonlinear Dynamics 85, 1613- 1633.

  • [47]

    Zotos, E. E., (Oct 2016(c)), Investigating the Newton-Raphson basins of attraction in the restricted three body problem with modified Newtonian gravity. J. Appl. Math. Comput 1-19.

  • [48]

    Zotos, E. E.,(2017), Revealing the basins of convergence in the planar equilateral restricted four body problem. Astrophys. Space Sci 362(2).

[1]

Ansari, A. A, (2014),Periodic orbits around lagrangian points of the circular restricted four-body problem. Invertis Journal of Science and Technology. 7(1), 29-38.

[2]

Abouelmagd, E. I., Mostafa, A.,(2015), Out of plane equilibrium points locations and the forbidden movement regions in the restricted three-body problem with variable mass. Astrophys. Space Sci 357, 58,

[3]

Ansari, A. A., (2016(a)), Stability of the equilibrium points in the photogravitational circular restricted four body problem with the effect of perturbations and variable mass.Science International (Lahore). 28, 859- 866.

[4]

Ansari, A. A., (2016(b)), Stability of the equilibrium points in the circular restricted four body problem with oblate primary and variable mass. International Journal of Advanced Astronomy. 4(1), 14-19.

[5]

Ansari, A. A.(2016 (c)), The Photogravitational Circular Restricted Four-body Problem with Variable Masses. Journal of Engineering and Applied Sciences. 3(2), 30-38.

[6]

Ansari, A. A.,(2017(a)), The circular restricted four body problem with variable masses. Nonlinear Sci. Lett. A. 8(3), 303-312.

[7]

Ansari, A. A.,(2017(b)), Effect of Albedo on the motion of the infinitesimal body in circular restricted three-body problem with variable masses. Italian J. of Pure and Applied Mathematics. 38, 581-600.

[8]

Anselmo, L., Farinella, P., Milani, A., Nobili, A. M.,(1983), Effects of the earth-reflected sunlight on the orbit of the LAGEOS Satellite. Astron. Astrophys. 117, 3-8.

[9]

Arribas, M., Abad, A., Elipe, A., Palacios, M.,(2016), Out-of-plane equilibria in the symmetric collinear restricted four-body problem with radiation pressure.Astrophys. Space Sci 361, 270(8). 10.1007/s10509-016-2858-1

[10]

Asique, M. C., et al.,(2015(a)),On the R4BP when third primary is an oblate spheroid. Astrophys. Space Sci 357, 82(1).

[11]

Asique, M. C., et al.,(2015(b)),On the photogravitational R4BP when the third primary is an oblate/prolate spheroid. Astrophys. Space Sci 360, 13(1),

[12]

Asique, M. C., et al., (2016), On the photogravitational R4BP when the third primary is a triaxial rigid body. Astrophys. Space Sci 361, 379.

[13]

Asique, M. C., et al.,(2017), On the R4BP when Third Primary is an Ellipsoid. Journal of Astronaut. Sci 64, 231-250,

[14]

Baltagiannis, A., Papadakis, K. E.,(2011), Equilibrium points and their stability in the restricted four body problem. International Journal of Bifurcation and Chaos. 21(8), 2179-2193,

[15]

Douskos, C. N.,(2010), Collinear equilibrium points of Hill's problem with radiation pressure and oblateness and their fractal basins of attraction. Astrophys. Space Sci 326, 263-271.

[16]

Falaye, B. J.,(2015), Effect of Oblateness, Radiation and a Circular Cluster of Material Points on the Stability of Equilibrium Points in the Restricted Four-Body Problem. Few body system. 56, 29-40,

[17]

Hadjidemetriou, J. D.,(1980), The Restricted Planetary 4-Body Problem.Celest. Mech. 21, 63-71.

[18]

Idrisi, M. J.(2017), A study of Libration Points in Modified CR3BP under Albedo Effect when smaller Primary is an Ellipsoid. J. of Astronaut. Sci,

[19]

Jeans, J. H., (1928), Astronomy and Cosmogony, Cambridge University Press, Cambridge.

[20]

Kalvouridis, T. J., Mavraganis, A. G.,(1995), Equilibria and stability of the restricted photogravitational problem of 2 +  2 bodies. Astrophys. Space Sci 226(1), 137-148.

[21]

Kalvouridis, T. J.,(1996),The oblate spheroids version of the restricted photogravitational 2 + 2 body problem. Astrophys. Space Sci 246(2), 219-227.

[22]

Kalvouridis, T. J., Arribas, M., Elipe, A.,(2006),Dynamical properties of the restricted four-body problem with radiation pressure. Mechanics Research Communications. 33, 811-817.

[23]

Kumari, R., Kushvah, B. S.,(2014), Stability regions of equilibrium points in the restricted four body problem with oblateness effects. Astrophys. Space Sci 349, 693-704.

[24]

Lichtenegger, H.,(1984),The dynamics of bodies with variable masses. Celest. Mech 34, 357-368.

[25]

Lukyanov, L. G.,(2009), On the restricted circular conservative three-body problem with variable masses. Astronomy Letters. 35(5), 349-359.

[26]

Meshcherskii, I. V.,(1952), Works on the Mechanics of Bodies of Variable Mass, GITTL, Moscow.

[27]

Michalodmitrakis, M.,(1981), The circular restricted four-body problem. Astrophys. Space Sci 75, 289-310.

[28]

Mittal, A., et. al.,(2016), Stability of libration points in the restricted four-body problem with variable mass. Astrophys. Space sci 361, 329,

[29]

Moulton, F. R., (1914), An introduction to celestial mechanic, Second ed. Dover, New York.

[30]

Papadakis, K. E.,(2016), Families of three-dimensional periodic solutions in the circular restricted four-body Problem. Astrophys. Space Sci 361, 129.

[31]

Papadouris, J. P., Papadakis, K. E.,(2013), Equilibrium points in the photogravitational restricted four-body Problem. Astrophys. Space Sci 344, 21-38.

[32]

Papadouris, J. P., Papadakis, K. E.,(2014), Periodic solutions in the photogravitational restricted four-body problem. MNRAS. 442, 1628-1639.

[33]

Paricio, L. J. H.,(2016), Bivariate Newton-Raphson method and toroidal attraction basins. Numerical Algo 71, 349-381.

[34]

Rocco, E. M.,(2009), Evaluation of the terestrial Albedo-Applied to some scientific missions. Space Sci. Rev 151, 135-147,

[35]

Shrivastava, A. K., Ishwar, B.,(1983),Equations of motion of the restricted problem of three bodies with variable mass. Celest. Mech 30, 323-328.

[36]

Simo, C.,(1978), Relative equilibrium solutions in the four body problem. Celest. Mech 18, 165-184.

[37]

Singh, J., Ishwar, B.,(1984), Effect of perturbations on the location of equilibrium points in the restricted problem of three bodies with variable mass. Celest. Mech 32, 297-305.

[38]

Singh, J., Ishwar, B.,(1985), Effect of perturbations on the stability of triangular points in the restricted problem of three bodies with variable mass. Celest. Mech 35, 201-207.

[39]

Singh, J.,(2003),Photogravitational restricted three body problems with variable mass. Indian Journal of Pure and Applied Math 32 (2), 335-341.

[40]

Singh, J., Leke, O.,(2010), Stability of photogravitational restricted three body problem with variable mass. Astrophys. Space Sci 326 (2), 305-314.

[41]

Singh, J., Leke, O.,(2013), Existence and Stability of Equilibrium Points in the Robe's Restricted Three-Body Problem with Variable Masses. International Journal of Astronomy and Astrophysics. 3, 113-122.

[42]

Singh, J., Vincent, A. E.,(2016(a)),Equilibrium points in the restricted four-body problem with radiation pressure. Few-Body Syst 57, 83-91.

[43]

Singh, J.,Vincent, A. E.,(2016(b)),Out-of-plane equilibrium points in the photogravitational restricted four-Body problem with oblateness. British Journal of Mathematics and Computer Science. 19(5), 1-15.

[44]

Zhang, M. J., Zhao, C. Y., Xiong, Y. Q.,(2012), On the triangular libration points in photo-gravitational restricted three body problem with variable mass. Astrophys. Space Sci 337, 107-113.

[45]

Zotos, E. E.,(2016(a)), Fractal basins of attraction in the planar circular restricted three body problem with oblateness and radiation pressure. Astrophys. Space Sci 181(17).

[46]

Zotos,E. E., (2016(b)), Fractal basins boundaries and escape dynamics in a multi-well potential. Nonlinear Dynamics 85, 1613- 1633.

[47]

Zotos, E. E., (Oct 2016(c)), Investigating the Newton-Raphson basins of attraction in the restricted three body problem with modified Newtonian gravity. J. Appl. Math. Comput 1-19.

[48]

Zotos, E. E.,(2017), Revealing the basins of convergence in the planar equilateral restricted four body problem. Astrophys. Space Sci 362(2).

Journal Information

Figures

  • View in gallery

    Configuration of the problem in CRFVBP with the effect of oblateness and Albedo.

  • View in gallery

    The locations of lagrangian points during in-plane motion (i.e. ξ≠0, η≠0, ζ = 0) in five cases: (a). Third primary is placed at one of the lagrangian points of the classical circular restricted three-body problem (points with blue color), (b). Variation of masses (points with green color), (c). Solar radiation pressure (points with red color), (d). Albedo effect (points with black color), (e). Oblateness effect (points with magenta color).

  • View in gallery

    The locations of lagrangian points during out of plane motion (i.e. ξ ≠ 0, η = 0, ζ ≠ 0) in five cases: (a). Third primary is placed at one of the lagrangian points of the classical circular restricted three-body problem (points with blue color), (b). Variation of masses (points with green color), (c). Solar radiation pressure (points with red color), (d). Albedo effect (points with black color), (e). Oblateness effect (points with magenta color).

  • View in gallery

    The locations of lagrangian points during out of plane motion (i.e.ξ = 0, η ≠ 0, ζ ≠ 0) in five cases: (a). Third primary is placed at one of the lagrangian points of the classical circular restricted three-body problem (points with blue color), (b). Variation of masses (points with green color), (c). Solar radiation pressure (points with red color), (d). Albedo effect (points with black color), (e). Oblateness effect (points with magenta color).

  • View in gallery

    Poincaré surface of sections for five cases: (a). Third primary is placed at one of the lagrangian points of the classical circular restricted three-body problem (with blue color), (b). Variation of masses (with green color), (c). Solar radiation pressure (with red color), (d). Albedo effect (with black color), (e). Oblateness effect (with magenta color). (i) Poincaré surface of sections in ξξ′-plane, (ii) Poincaré surface of sections in ηη′-plane, (iii) Zoomed part of figure (ii) near the origin.

  • View in gallery

    (a): The basin of attraction for the case when third primary is placed at one of the lagrangian points of the classical circular restricted three- body problem. (b): Zoomed image of (a) near the lagrangian configuration.

  • View in gallery

    (a): The basin of attraction for the variable mass case. (b): Zoomed image of (a) near the lagrangian configuration.

  • View in gallery

    (a): The basin of attraction for the solar radiation pressure case. (b): Zoomed image of (a) near the lagrangian configuration.

  • View in gallery

    (a): The basin of attraction for the Albedo case. (b): Zoomed image of (a) near the lagrangian configuration.

  • View in gallery

    (a): The basin of attraction for the oblateness case. (b): Zoomed image of (a) near the lagrangian configuration.

Metrics

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 122 122 57
PDF Downloads 35 35 6