SciELO - Scientific Electronic Library Online

 
vol.57 número2Estimating the propagation of a uniformly accelerated jetA study of the secular variation of the high-amplitude Delta Scuti star AD CMi índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de astronomía y astrofísica

versión impresa ISSN 0185-1101

Rev. mex. astron. astrofis vol.57 no.2 Ciudad de México oct. 2021  Epub 15-Feb-2022

https://doi.org/10.22201/ia.01851101p.2021.57.02.05 

Artículos

Location and stability of the triangular points in the triaxial elliptic restricted three-body problem

M. Radwan1 

Nihad. S. Abd El Motelp2 

1 Astronomy; Space Science and Metreology Department, Faculty of Science, Cairo University, 12613 Giza, Egypt.

2 Astronomy Department, National Research Institute of Astronomy and Geophysics (NRIAG), 11421 Helwan, Cairo, Egypt.


ABSTRACT

The main goal of the present paper is to evaluate the perturbed locations and investigate the linear stability of the triangular points. We studied the problem in the elliptic restricted three body problem frame of work. The problem is generalized in the sense that the two primaries are considered as triaxial bodies. It was found that the locations of these points are affected by the triaxiality coefficients of the primaries and the eccentricity of orbits. Also, the stability regions depend on the involved perturbations. We also studied the periodic orbits in the vicinity of the triangular points.

Key Words: celestial mechanics; gravitation; methods; analytical; methods; numerical; planets and satellites; dynamical evolution and stability

RESUMEN

Nuestro objetvo principal es determinar las ubicaciones perturbadas e investigar la estabilidad lineal de los puntos triangulares. Realizamos nuestro estudio en el marco del problema elíptico restringido de tres cuerpos. Generalizamos el problema considerando a los dos cuerpos primarios como cuerpos triaxiales. Encontramos que la ubicación de los puntos triangulares se ve afectada por los coeficientes de triaxialidad de los cuerpos primarios y por la excentricidad de las órbitas. También notamos que las regiones de estabilidad dependen de las perturbaciones inducidas y estudiamos las órbitas periódicas en la vecindad de los puntos triangulares.

1. INTRODUCTION

The restricted three-body problem (R3BP) is the most significant problem in celestial mechanics. This due to its wide applications in space dynamics and solar system dynamics. The knowledge of the R3BP is crucial in almost all space applications. We can define this problem as a dynamical system in which an infinitesimal body moves in the gravitational field of two massive bodies (primaries). The two primaries move about their common centre of mass in either circular or elliptic orbits, which leads to the famous circular or elliptic restricted problems, respectively (Musielak & Quarles 2014). In the elliptic restricted three-body problem, ER3BP, it is assumed that the smaller primary moves around the more massive one in an elliptical orbit. The influence of the infinitesimal mass on the motion of the primaries is negligible. This significantly simplified model has more applications in space dynamics than the more general problem.

Despite the fact that the RTBP is not integrable, five specified solutions in the rotating frame are known. These solutions correspond to equilibrium positions. Three of these equilibrium points, known as collinear points, are located on the line joining the two massive primaries in the rotating reference frame. The remaining two points, called triangular points, form equilateral triangles with the massive primaries. The determination of the location of these libration points is of great interest for different space applications (Abd El-Salam 2019). In general the equilateral equilibrium points are conditionally stable, while the collinear libration points are usually unstable (Singh & Tyokyaa 2016).

Most of the celestial bodies are of irregular nature and cannot be assumed to be spherical in the RTBP, because their stability of movement is influenced by their shape. The planets of the solar system and their satellites are assumed as extended bodies and cannot be regarded as spheres. Furthermore, the effect of eccentricity on the orbits is significant. Hence, in order to achieve an acceptable accuracy, we have to take into account that the orbits of most celestial bodies are elliptic rather than circular (Singh & Tyokyaa 2016).

From the theory of small oscillations around the triangular points two types of periodic orbits can be obtained; short-period orbits with periods approximately the same as the orbital period of the primaries, and long-period orbits with periods many times that of the primaries. The importance of investigating these periodic orbits lies in their frequent appearance in nature; they provide us with important information about orbital resonances, spin orbits and are also useful to approximate quasi-periodic trajectories (Abouelmagd (2013), Mittal et al. (2009), Kumar & Sharma (2020), Ahmed et al. (2020).

A number of studies have been performed on the equilateral triangular points of the restricted threebody problem. Singh & Begha (2011) investigated the existence of periodic orbits about the triangular points, namely L 4 and L 5, in the RTBP when the primaries are considered as triaxial and oblate spheroid bodies, besides perturbations due to Coriolis and centrifugal forces. The authors found that long and short periodic orbits exist around these points, and that their eccentricities and periods are influenced by the involved perturbations. Singh & Umar (2012) studied the stability of triangular points in the ER3BP, when both oblate primaries emit light energy. The locations of L 4 and L 5 are shifted away from the line joining the two primaries, as compared to the classical case. The stability region decreases or increases with variations in eccentricity, oblateness and radiation pressure. Naryan & Usha (2014) investigated the stability of infinitesimal motions around the equilateral points. They considered the more massive primary as a radiant source and the less massive one as a triaxial body. Zahra et al. (2016), Zahra et al. (2017) studied the locations and linear stability of the triangular points under the effects of oblateness and triaxiality of the primaries, plus small relativistic perturbations. They found that the locations of the triangular points are affected by the considered perturbations. Numerical explorations were carried out to show the effect of perturbations on the position and stability. Singh & Tyokyaa (2016) investigated the positions and stability of the libration points when both primaries are taken as oblate spheroids with oblateness up to the fourth zonal harmonic J 4. They found that both harmonic coefficients J 2 and J 4, eccentricity and semimajor axis have destabilizing tendencies and as a result the size of the region of stability decreases with an increase in the parameters involved.

Recently, Hadia et al. (2019) studied analytically the existence and the linear stability of the libration points assuming the primaries as triaxial bodies and the Euler angles as specified values. The authors proved that the positions and the stability of the triangular points change according to the effect of the triaxiality of the primaries. Furthermore, they presented the solution of long and short periodic orbits for a stable motion. Cárcamo-Díaz et al. (2020) revisited the circular spatial restricted problem. The authors presented some new results in the light of the concept of Lie stability. This paper studies the positions and stability of the well-known triangular points, in the elliptic restricted problem frame of work. The triaxiality of the primary bodies is taken into consideration. In addition, we study the periodic orbits around the triangular points. To examine the present problem we carry out several numerical explorations.

2. DYNAMICAL MODEL

Let m 1, m 2 and m 3 be the masses of the primary, secondary, and the third infinitesimal body, respectively. The third body moves under the effect of the triaxial primaries, but it does not affect their motion. The two primaries move around their common centre of mass in elliptic orbits. We consider a rotating coordinate system (x,y,z) with origin at the centre of mass of the primaries. The primaries are permanently located on the x-axis, which is the line joining them. When investigating the restricted three-body problem it is better to choose a system of normalized units. So, we assume that the distance between the primaries is unity and we choose also the unit of time to make the gravitational constant unity. The total mass of the primaries is considered as 1, i.e. m 1 + m 2 = 1. In this normalized system we adopt μ=m2m1+m2 as a dimensionless mass parameter. According to the above considerations, the equations of motion of the third infinitesimal body moving under the effect of the triaxial primaries located at (µ,0), (1−µ,0) are given by (Ahmed et al. 2020)

x¨-2ny˙=Ux,y¨+2nx˙=Uy, (1)

with

U=n22[(1-μ)r12+μr22]+(1-μ)r1+μr2+(1-μ)(2σ1-σ2)2r13-3(1-μ)(σ1-σ2)y22r15+μ(2γ1-γ2)2r23-3μ(γ1-γ2)y22r25. (2)

r1=(x+μ)2+y2,r2=(x+μ-1)2+y2, (3)

n=1a(1-e2){1+32[(2σ1-σ2)+(2γ1-γ2)]}. (4)

Here, U is the restricted three-body potential; n is the mean motion of the primaries; r 1 and r 2 are the distances of the infinitesimal mass from the primaries; σ 1, σ 2 and γ 1, γ 2 are the triaxiality coefficients of the bigger and smaller primaries, respectively.

3. LOCATIONS OF THE TRIANGULAR POINTS

The libration points represent stationary solutions of the restricted three-body problem. The positions of these points can be obtained by setting all components of the relative velocity and acceleration equal to zero, i.e. we found them by setting = = 0 = = ÿ. Consequently, the locations of these points can be found by solving simultaneously the nonlinear equations U x = U y = 0, and y 6= 0. Setting

Aγ=γ1-γ2,Aσ=σ1-σ2,Aα=2γ1-γ2,Aβ=2σ1-σ2,

we have

Ux=-μ(-1+x+μ)[1r23+32r25Aα-15y22r27Aγ]+xa(1-e2)+(1-μ)(x+μ)[1r13-32r15Aβ+32a(1-e2)(Aα+Aβ)-15y22r17Aσ], (5)

Uy=-y{μr23+1-μr13+3μ2r25Aα+3(1-μ)2r15Aβ-1a(1-e2)[1+32(Aα+Aβ)]-15y22[μAγr27+(1-μ)r17Aσ]-3μr25Aγ+3(1-μ)r15Aσ}, (6)

Ignoring the triaxiality perturbations due to the primary bodies yields the equilateral solution of the classical restricted three-body problem i.e. r 1 = r 2 = 1. Then, it may be reasonable to consider that the locations of these points are the same as given by classical problem, but perturbed by small terms, i.e.

ri=1+δi,δi<<1,(i=1,2). (7)

Substituting this assumption (7) into (3), solving for x and y and retaining terms up to the first order in the small quantities δ i we get:

x=12(2δ1-2δ2-2μ+1),y=±12(3+4(δ1+δ2)). (8)

Substituting the values of r 1, r 2, x and y into equations (5) and (6) and ignoring the terms higher than order one in δ i we obtain the following two simultaneous equations in δ 1 and δ 2.

A1+A2δ1+A3δ2=0,B1+B2δ1+B3δ2=0. (9)

The corresponding solution is:

δ1=-A3B1-A1B3A3B2-A2B3,A3B2-A2B30, (10)

δ2=--A2B1+A1B2A3B2-A2B3,A3B2-A2B30, (11)

where:

A1=-12+μ+34(μAα-Aβ+μAβ)+34a(1-e2)(Aα+Aβ)+32a(1-e2)(1-μ(2+3Aα+3Aβ))+4516(-μ(Aγ+Aσ)+Aσ), (12)

A2=12-1a(1-e2)(1+32(Aα+Aβ))+32(1-μAα+32(1-μ)Aβ-558(1-μ)Aσ+54μAγ), (13)

A3=1+32Aβ+1a(1-e2)(1+32(Aα+Aβ))-158Aσ-32μ(1+32Aα+Aβ-558-54Aσ), (14)

B1=-32{1-1a(1-e2)[1+32(Aα+Aβ)]-32[Aβ-74Aσ+μ(Aα-Aβ-74Aγ+74Aσ)]}, (15)

B2=723+13a(1-e2)[1+32(Aα+Aβ)]+32[132Aβ-1218Aσ-μ(3+Aα+132Aβ-374Aγ-1218Aσ)], (16)

B3=-13(1+1a(-1+e2))+332[374Aσ-Aβ3+μ(1+132Aα+Aβ-1218Aγ-374Aσ)-(Aα+Aβ)a(-1+e2)]. (17)

4. GRAPHICAL REPRESENTATIONS (LOCATION)

In Figures 1-2 we plot the locations r L4,5 of the triangular points in the elliptic restricted three body problem, taking into consideration triaxiality and eccentricity effects. Figure 1 represents the variation of r 4,5 against the mass ratio µ. The calculations are carried out for several dynamical models with different coefficients γ i and σ i ; and constant eccentricity of the orbit of the primary e = 0.2. As is clear from the figure, the locations begin to decrease sharply until we reach the value µ ≈ 0.15, then they decrease slowly. We also observe that the variation in the locations r L4,5 is small, because the triaxiality coefficients are very small. Figure 2 depicts the variation of r L4,5 against the mass ratio µ with constant triaxiality coefficients. We applied our case study to a range of eccentricities; the eccentricity e takes the values 0.01,0.04,0.06,0.08 and 0.1. We notice that the change in the curves is nearly the same for all considered cases. This is because the chosen values of e are small. The dynamics is approximately the same but with different sizes of the perturbations. We expect that the larger the eccentricity the larger will be the perturbation in the positions of the triangular points.

Fig. 1 The locations of L 4,5 with different triaxiality coefficients and eccentricity e = 0.2. The color figure can be viewed online.  

Fig. 2 The locations of L 4,5 with different eccentricities and triaxiality coefficients σ 1 = 0.0002, σ 2 = 0.0005, γ 1 = 0.003 and γ 2 = 0.002. The color figure can be viewed online.  

5. STABILITY OF THE TRIANGULAR POINTS

The dynamical models which describe the restricted problem are very difficult. The great complexity of these models leads us to focus our attention on linearized dynamics, because in this case we obtain simplified mathematical expressions that could be handled easily. The location of the third body is displaced little from the equilibrium point due to the considered perturbations. If the resultant motion of this body is a rapid departure from this point, we call such a position of equilibrium point an unstable one; if the body just oscillates about the point, it is said to be a stable position (Abd El-Salam 2019). To examine the possible motion of the infinitesimal body in the proximity of the triangular points, the equations of motion should be linearized. Let the position of the triangular points be (x ,y ) and let the third body be displaced to the point (x +ξ, y +η), where ξ and η are small displacement in (x ,y ), so that the linearized equations can be written as:

ξ¨-2nη˙=1a(1-e2)(UxxL4,5ξ+UxyL4,5η),η¨-2nξ˙=1a(1-e2)(UyyL4,5η+UxyL4,5ξ). (18)

Here UxxL4,5 denotes the second derivative of U with respect to x computed at the stationary solution and the rest of the derivatives are defined similarly. The characteristic equation corresponding to equation (18) may be written as

λ4+(4n2-UxxL4,5-UyyL4,5)λ2+UxxL4,5UyyL4,5-(UxyL4,5)2=0, (19)

which may be written in the form

λ4+ULiλ2+νLi=0, (20)

with

ULi=(4n2-UxxL4,5-UyyL4,5),νLi=UxxL4,5UyyL4,5-(UxyL4,5)2. (21)

The roots of equation (21) are

λ1,2=±-ULi-ULi2-4νLi2, (22)

and

λ3,4=±-ULi+ULi2-4νLi2. (23)

The possible solutions for equation (19) will determine the stability of the equilibrium points. We have three possibilities: asymptotically stable, when all the solutions are real and negative; unstable, when at least one of the solutions is positive and real; and finally, the stable case when all solutions are purely imaginary (Szebehely 1967). We will focus our attention only on the case with oscillatory stable solutions about the triangular points. To investigate the stability/instability we first compute the partial derivatives required for equations (22) and (23) as

Uxx=1a(1-e2)[1+32(Aα+Aβ)]+μ{3(μ+x-1)2r25-1r23+32Aα[5(μ+x-1)2r27-1r25]-152y2Aγ[7r29(μ+x-1)2-1r27]}+(1-μ){32Aβ[5(μ+x)2r17-1r15]+[3(μ+x)2r15-1r13]-152y2Aσ[7(μ+x)2r19-1r17]}, (24)

Uyy=1a(1-e2)[1+32(Aα+Aβ)]+μ{3Aγ[10y2r27-5y22(7y2r29-1r27)-1r25]}+(1-μ){3y2r15-1r13+32(5y2r17-1r15)Aβ-3Aσ[1r15+5y2(-2r17+7y22r19-12r17)]+3Aα2(5y2r27-1r25)+3y2r25-1r23}, (25)

Uxy=y{μ(μ+x-1)[-105y2Aγ2r29+15Aα2r27+15Aγr27+3r25]+y(1-μ)(μ+x)[3r15-105y2Aσ2r19+15Aβ2r17+15Aσr17]}. (26)

6. GRAPHICAL REPRESENTATIONS (STABILITY)

Figures 3-5 illustrate the stability regions taking into consideration the triaxiality coefficients and the eccentricity of the orbits. Figure 3 represents the classical case i.e. when σ i = 0, γ i = 0 and e = 0 as well. Figure 4 depicts the stability regions for dynamical systems with different triaxiality coefficients, while Figure 5 shows these regions for different eccentricities and constant coefficients. We can see from Figure 4 that each stability region is characterized by different values of the triaxiality coefficients. The size of the regions is not much different from that of the classical case; this because the perturbation due to the triaxiality coefficients is small. In Figure 5 the effect of different eccentricities on the size of the regions is obvious. Comparing with previous work e.g. Abd El-Salam (2015), we found that the size of the regions is slightly different from that found in the quoted paper due to the perturbations considered in each case. We expect that as the eccentricity increases the stability regions will be destroyed. Table 1 gives the values of the critical mass for some selected dynamical models. These values are computed under the effect of the included perturbations. We notice from the table that when the triaxiality coefficients (σ i and γ i ) are equal to zero, µ c reduces to the critical mass value of the classical restricted problem (see Figure 3). We also notice from Figures 4, 5 that the triaxiality coefficients of the primaries change the range of stability.

Fig. 3 The imaginary part of the eigenvalues versus the mass parameter in the classical case, σ i = γ i = 0. The color figure can be viewed online.  

Fig. 4 The imaginary part of the eigenvalues versus the mass parameter with eccentricity e = 0.01 and different triaxial parameters. The color figure can be viewed online. 

Fig. 5 The imaginary part of the eigenvalues versus the mass parameter with different eccentricities and triaxial parameters, σ 1 = 0.001, σ 2 = 0.004, γ 1 = 0.003 and γ 2 = 0.002. The color figure can be viewed online. 

TABLE 1 VALUES OF THE CRITICAL MASS RATIO FOR DIFFERENT DYNAMICAL MODELS 

σ1 σ2 γ1 γ2 e a μc
0 0 0 0 0 1 0.0385209
0.006 0.0004 0.03 0.002 0.2 0.92 0.0443621
0.006 0.0004 0.03 0.002 0.4 0.92 0.0401412
0.006 0.0004 0.03 0.002 0.6 0.92 0.0132093
0.003 0.0002 0.007 0.0004 0.02 0.92 0.0456447
0.003 0.0002 0.007 0.0004 0.04 0.92 0.0424109
0.003 0.0002 0.007 0.0004 0.07 0.92 0.00909276
0.009 0.002 0.008 0.002 0.02 0.92 0.042956
0.009 0.004 0.008 0.002 0.02 0.92 0.0425331
0.009 0.008 0.008 0.002 0.02 0.92 0.0416678
0.006 0.001 0.008 0.001 0.08 0.92 0.0466464
0.006 0.001 0.008 0.005 0.08 0.92 0.0431538
0.006 0.001 0.008 0.008 0.08 0.92 0.0427982

In addition, we can also investigate the behavior of the system by plotting the critical mass ratio against the eccentricity for different values of the triaxial parameters σ i and γ i . We plot different combinations of triaxial parameters and a stable region is observed. Figure 6 shows that the increase in the eccentricity reduces the region of stability of the infinitesimal body about the triangular points. Figures 7 and 8 represent the variation of the critical mass with the triaxial parameters when the eccentricity is constant. We can see from the two figures that increasing the triaxial parameters results in a destabilizing effect in the system.

Fig. 6 The critical mass ratio versus eccentricity for different values of the triaxial parameters. The color figure can be viewed online. 

Fig. 7 The critical mass ratio versus triaxial parameter σ 2, for different values of the eccentricity. The color figure can be viewed online. 

Fig. 8 The critical mass ratio versus triaxial parameter γ 2, for different values of the eccentricity. The color figure can be viewed online. 

7. MOTION IN THE PROXIMITY OF THE TRIANGULAR POINTS

It is known that periodic orbits represent the backbone in the study of dynamical systems in celestial mechanics. Since in the range 0 ≤ µµ c , the roots of the characteristic polynomial of the present system are purely imaginary, we conclude that the motion around the triangular points L 4 (or L 5) is stable and composed of two harmonic motions. The solutions for the perturbed motion about the triangular point L 4 (in the synodic reference frame) may be written as

X[t]=j=14αjeλjt, (27)

Y[t]=j=14βjeλjt, (28)

where the terms with coefficients α j and β j represent the short and long periodic motion. The constants α j and β j are not independent. Applying the initial conditions to the solutions of the system these solutions may be written in the form

X[t]=5.80204*10-5sin(0.241582t)-1.31958*10-5sin(1.06221t)+3.19255*10-5cos(0.241582t)-2.19255*10-5cos(1.06221t). (29)

Y[t]=8.06288*10-5cos(0.241582t)-7.062876*10-5cos(1.06221t)-3.054769*10-5sin(0.241582t)+6.9475842*10-5sin(1.06221t). (30)

We may note from equations (29) and (30) that the solution is of the oscillatory type with fundamental periods of 0.241582 2 π and 2π1.06221 . times the orbital periods of the orbiting body. Hence the body will remain in the proximity of the triangular point L 4 and the motion is stable. These two different types of frequencies in the solution exist due to the fact that the resulting motion of the small body is composed of two types of motion. Firstly, a short-period motion with a period close to the orbital period of the less massive primary, and secondly a superimposed longperiod one known as a libration around the point L 4. The two solutions are depicted in Figures 9 and 10. The graph of the trajectory in the vicinity of equilibrium point L 4 is shown in Figures 11 and 12 with magenta color. The graphs contain both the short and the long-period motions. Figure 11 shows a comparison between the long-period motions in both the classical case (blue curve) and the present dynamical system (magenta curve), while Figure 12 shows a comparison between the two systems in the shortperiod case. It is observed that the trajectories (magenta curve) in both figures are shifted from those of the classical case. This is due to the included perturbations.

Fig. 9 x[t] as a function of time with triaxiality coefficients σ 1 = 0.007, σ 2 = 0.002, γ 1 = 0.005, γ 2 = 0.002, µ = 0.01, e = 0.01 and a = 0.92. The color figure can be viewed online. 

Fig. 10 y[t] as a function of time with triaxiality coefficients σ 1 = 0.007, σ 2 = 0.002,γ 1 = 0.005, γ 2 = 0.002, µ = 0.01, e = 0.01 and a = 0.92. The color figure can be viewed online.  

Fig. 11 Comparison between the long period motion in both the classical case (blue curve) and the present dynamical system (magenta curve). The color figure can be viewed online.  

Fig. 12 Comparison between the short period motion in both the classical case (blue) and the present dynamical system (magenta). The color figure can be viewed online.  

8. CONCLUSIONS

We have studied the existence and linear stability of the triangular points in the elliptic restricted threebody problem when the primaries are triaxial bodies. We found that the locations and the linear stability of the triangular points are affected by the triaxiality coefficients of the primaries and the eccentricity of the orbit. We also noticed that the size of the stability regions depends on the triaxiality coefficients and the eccentricity of the orbits. In addition, we studied the periodic orbits in the vicinity of the triangular points. We showed that these orbits are elliptical. Further, we observed that the shape of the periodic orbits changed due to the triaxiality of the primaries and showed a deviation from the classical case.

REFERENCES

Abd El-Salam, F. A., 2015, Ap&SS, 357, 15 [ Links ]

______. 2019, NewA, 68, 25 [ Links ]

Ahmed, A. A., Selim, H. H., Kamel, A. K. G., et al. 2020, Applied Mathematics and Nonlinear Sciences, 5, 157 [ Links ]

Abouelmagd, Elbaz I. 2013, EM&P, 110, 143 [ Links ]

Cárcamo-Díaz, D., Palacián, J. F., Vidal, C., & Yanguas, P. 2020, RCD, 25, 131 [ Links ]

Kumar, P. & Sharma, R. K. 2020, Indian Journal of Science and Technology, 13, 32 [ Links ]

Mittal, A., Ahmad, I., & Bhatnagar, K. B., 2009, Ap&SS , 323, 65 [ Links ]

Musielak, Z. E. & Quarles, B. 2014, RPPh, 77, 6, 065901 [ Links ]

Narayan, A. & Usha, T. 2014, Ap&SS , 351, 135 [ Links ]

Selim, H. H., Guirao, J. L. G., & Abouelmagd, E. I. 2019, Discrete & Continuous Dynamical Systems-S, 12, 703 [ Links ]

Singh, J. & Begha, J. M. 2011, Ap&SS , 332, 319 [ Links ]

Singh, J. & Umar, A. 2012, AJ, 143, 109 [ Links ]

Singh, J. & Tyokyaa, R. K. 2016, EPJP, 131, 365 [ Links ]

Szebehely, V. G. 1967, Theory of orbits. The restricted problem of three bodies (New York, NY: Academic Press) [ Links ]

Zahra, K., Radwan, M., & Awad, M. E. 2016, Nonlinear Analysis and Differential Equations, 4, 609 [ Links ]

Zahra, K., Awad, Z., Dwidar, H. R., & Radwan, M. 2017, SerAJ, 195, 47 [ Links ]

Received: February 04, 2021; Accepted: April 12, 2021

Nihad S. Abd El Motelp: Astronomy; Astronomy Department, National Research Institute of Astronomy and Geophysics (NRIAG), 11421 Helwan, Cairo, Egypt (nihad.saad@nriag.sci.eg).

M. Radwan: Astronomy; Space Science and Metreology Department, Faculty of Science, Cairo University, 12613 Giza, Egypt (mradwan@sci.cu.edu.eg).

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License