Abstract

The restricted three body problem was outlined. The acceleration due to planetary magnetic field, in terms of space craft’s orbital elements, was analysed. The conditions for calculating the liberation points including the mutual gravitational attraction and the effect of Lorentz acceleration were derived for the case of circular planer restricted three bodies. The stability of the solution for the artificial Lorentz triangular liberation points was studied. Finally, numerical investigation for the case of Sun-Jupiter system was calculated as case study. The results show the ability of changing the position of the triangular liberation points by an order from 10−7 to 10−6 for the dimensionless x, y coordinates and distance r from Jupiter. This is equivalent to about hundreds of Kilometers which is considerable.

1. Introduction

The history of the RTBP begins in 1772 by Euler and Lagrange, continues in 1836 by Hill [1], and is followed by Poincare [2], Levi-Civita [3], Birkhoff [4], and then by Szebhley [5]. Till now there are great names and important contributions to the RTBP.

Most perturbation techniques are applied to solve the equations of motion of the RTBP. Delva [6], used a Lie series to determine the orbits of massless bodies in the elliptic RTBP. Also Sandor et al. [7], applied the method of the short time Lyapunov indicators to the RTBP in order to study the structure of the phase space in some selected regions.

The idea of the effect of geomagnetic field on the rotation and motion of Earth’s artificial satellite was initiated by Fain and Greer [8]. Westerman, [9], used their idea to study the perturbation on circular orbit in the plane of the geomagnetic equator. Sehnal [10], concluded, at that time, no precise information available to evaluation of the effects of the Earth’s magnetic field on the orbit of a charged satellite.

The concept of generating AEP in the RTBP was initiated by Dusek [11]. The stability and location of the AEP were investigated since that time, from different aspects, by several authors [1216]. On the other side, McInnes and others produced an extensive subsequently researches to find the equilibrium positions using solar sail propulsion such as [1721]. Such work lead to find infinite equilibrium surfaces depending on the magnitude of the propulsive acceleration.

In the direction to finding the AEP using low-thrust systems Morimoto, et al. [22], proved the existence of infinite equilibrium surfaces depending on the magnitude of the propulsive acceleration.

However, only a subset of the potentially achievable AEPs turns out to be stable and, as thus, could not be exploited by a spacecraft without the use of a suitable control system. The topology of such subset of stable AEPs is strictly dependent on the propulsion system type employed by the spacecraft. In fact, as was recently pointed out by Bombardelli and Peláez [23], if the available propulsive acceleration is low, the stable AEPs are confined to a very restricted region around the classical Lagrange points.

In the last decade, several important articles worked on the effects of Lorentz force to control both orbit motion and attitude of a spacecraft whether alone or combined with other forces, for example, in 2011 Pollock et al. [24] proposed analytical solutions for the relative motion of a Lorentz spacecraft, while in the year, 2020, Mostafa et al. [25], studied the use of Lorentz force to control the perturbations due to solar radiation pressure, and Sun et al. [26], studied the coupled effect of Lorentz force with aerodynamics on controlling a spacecraft. Huang et al. [27], introduced a hybrid control system that contains both the Lorentz force and torque to control the orbit and attitude of a satellite.

In this article we study the effect of Lorentz force on the position and stability of the triangular liberation points in the restricted three body problem with an application given to the Sun-Jupiter system.

2. Formulation of the Problem

Consider the problem of three-bodies with the masses of the primaries and the mass of space craft (SC) in which is negligibly small compared to and . Assume that the primaries revolve in circular orbits around their center of mass with mean motion and constant separation L. Then, with is the Gaussian constant and .

In Figure 1, let the plane be the plane of motion of and , with X axis is the line joining and , Y axis is normal to X in their orbital plane while Z is perpendicular to the XY plane. Let be the coordinates of the SC in this coordinate system (the synodic coordinate system), then the equations of motion of SC in the framework of the circular restricted three body problem (CRTBP), using dimensionless units, become.

If we transform the coordinate system to the center of the second primary , with the coordinate x, y, and z, the equations of motion will be:

With,

Equations (2) and (3) give the base to calculate the liberation points (Lagrange points). To change the position of such points, additional acceleration, artificially, must be added to the equations. Let such acceleration be denoted by:with is adopted such as:

3. Lorentz Acceleration

The Lorentz acceleration for charged SC with mass m experienced by a charge q (Coulombs) moving through a planet magnetic field B, in plant center coordinate system and , is given by Brett Jordan Streetman [28]:where and are unit vectors in the equatorial plane of the planet while in the direction of rotational axis. With is the charge-to-mass ratio of the SC. in Coulombs per kilogram (C/kg). is the inertial spacecraft velocity while is the velocity of the magnetic field which is rotate with the spin speed , of the planet and can written as .

The general vector model of a dipole magnetic field is [29]:where a unit is vector along the north magnetic pole and is the strength of the field in Weber-meters. The unit vectors included in Eq. (12) are described in Figure 2.

For a tilted magnetic dipole, the angle α between and , in the Earth case is 11.7° and is independent of time.

If we neglect the tilt angle of the magnetic field (i.e. α=0), the unit vector will be coincidence with the unit vector . Then, we can write the last equation in the form:

Noting that from the spherical triangle ZbL, is related to the orbital elements, and equals and we will denote it “Γ”. Substituting Eq. (13) into Eq. (8), (9) and (10) to describe Lorentz acceleration as:

Completing the required vector operations, we get:where, and .

4. Controlling Equations

To control the position of liberation points using Lorentz force, we must add the Lorentz acceleration to the equations of motion and applying the stationary conditions.

If we neglect the inclination of the equatorial plane of the planet on its orbit, and choose the x-axis to be the line joining the Sun with planet, y-axis normal to it in the planet's orbital plane and z-axis normal to the orbital plane, then Eqns. (3), (4) and (15) will be described in the same coordinate system.

Substituting Eqns. (15) into Eqns. (3), (4) the equations of motion are derived in the form:

Equations (16), (17) and (18) are the governing equations for the circular restricted three body problem with the magnetic field of the primary considered, and the small mass is controlled through a charge q.

5. Dimensionless Form of Controlling Equations

In what follows, we put equations (16), (17) and (18) in the usual dimensionless form of the restricted three body problem so that it is easy to compare the case of artificial points by the addition of the Lorentz force with the natural points when q is set equal zero.

The equations are diveded by (L·n2) which is clearly the dimensions of acceleration to have it in dimensionless form, we further define the quantities:

t=n·t dimensionless time, and dimensionless distances. Also let , and finally define dimensionless angular velocity. We get

We note that the quantity is dimensionless since the units of is: CM−1 (MC−1T−1l3) which is after simplifying T−1l3 where we use C for unit of electric charge, M for unit of mass, T for time and l for length. So, we also define the dimensionless quantity Q=. Next, we substitute , then introduce the well-known parameter in the three-body problem, . The equations take the simpler dimensionless form,

Finally, we omit the stars for simplicity of writing, keeping in mind that all involved quantities are now dimensionless,

Equations (21), (22) and (23) are the dimensionless controlling equations of the motion of a small body with an electric charge, under the effect of two big bodies one of them has a significant magnetic field. The equations for planar problem are:

6. Artificial Equilibrium Points

The liberation points are found by letting . We also take L as the unit of distance = 1, and the unit of time, hence n= 1, thus Eqns (24) and (25) will give,

Equations (26) and (27) represent the controlling equations depending on the position of the space craft’s liberation pointes and the charge (per unit mass) required to generate Lorenz acceleration.

From Eq. (27), we get:, is the solution of collinear liberation points while . This leads to the equilateral triangle solution.

In the equilateral triangle solution, solving for r, we get:

In other side, solving Eq. (28) for Q we get:

We can express another useful relation since Q must satisfy both Eqns. (26) and (27), then:

Simplifying yields,

Which leads to conclude that the artificial liberation points due to Lorentz force generated from the small primary exist on a circular arc with center at the big primary with radius the unit distance, with angles where is the small deviation from the original point of equilateral triangular case, and the positive sign is for L4 while the negative sign is for L5, as shown in Figure 3.

The charge per unit mass , on the SC can be expressed from Eq. (30), as:

Applying the constraint (32), with , leads to:

When the body is not charged (Q = 0), we get the classical case of the Lagrangian triangular points L4 and L5 with r = 1.

7. Triangular Liberation Points Coordinates

We proved in the previous section that the liberation points exist. In this section we will calculate their position as a function of the charge, of the SC, per unit mass.

Using Eqns. (5) and (6), with , L = 1, a = µ, b = 1−µ, yields:

Equations (35) and (36) give the coordinates of the artificial triangular points in terms of the parameter µ, the magnetic field strength B0 of the primary at m1, the dimensionless angular velocity of the primary m1, and the charge per unit mass .

8. Stability of Artificial Lorentz Triangular Liberation Points

Recalling the equations of motion in planar case:

These equations can be written aswhere,

Then,

Now, we linearize the equations about the equilibrium points (x0, y0) by letting,where.

Then, we expand the equations keeping only first order terms. This gives,where,

The characteristic equation for the system (40) is then,

This can be written as:where,

The condition for to be imaginary is that is a real negative value. This is guaranteed by the satisfaction of the conditions:

Equations (47), (48) and (49) are the necessary conditions for the artificial liberation points to be stable.

9. Numerical Investigation for the Case of the Sun-Jupiter System

The magnetic field of Jupiter has the values Tesla meter3. The parameter µ= 0.000954, the mean motion n= 0.00006042 hour−1, the sun-Jupiter distance L= 779 × 106 km., and the angular velocity hour−1 [30]. The dimensionless quantity Q will thus have the value . Also, in this case . If the distance from the original triangular liberation points to the small primary (Jupiter in this case) be r = 1, the x and y coordinates X4,5 = −0.5 and Y4,5 =  . These values will differ by a small amount due to the existence of the Lorentz force. We plot the deviations of r, X4,5, and Y4,5 denoted by δr, δX4 and δY4 respectively against the charge per unit mass in Figure 4–Figure 6, where we note that X4 = X5 and Y4 = −Y5, so we plot only the graphs for L4. The results show a variation of an order from 10−7 to 10−6 for the dimensionless x, y coordinates and distance r from Jupiter. This is equivalent to about hundreds of kilometers which is considerable. The range of is from −0.1 to 0.1 C/kg.

Concerning the stability, the graphs of F, H, and F2-4H are plotted in the same range of the specific charge −0.1 to 0.1 to check conditions (47), (48), and (49). The graphs show values greater than zero in the whole range which guarantees the stability of the artificial triangular points as the case of the natural points. The results are given in Figure 7−Figure 9.

Finally, the Eigen values of equation (44) for different values of are given in Table 1.

10. Conclusion

The artificial triangular liberation points L4,5 in the circular restricted three body problem in the case when the small primary has a significant magnetic field are studied. In this work the control for the motion is taken the charge per unit mass of the charged space craft under study. The origin of the coordinates was taken at the small primary which has the center of the magnetic field. A numerical application is taken for the case of Sun-Jupiter system. The choice of this system is because of the small eccentricity of Jupiter’s orbit around the sun and the weak effects of other perturbing bodies on the motion, and most important for this study is because of the strong magnetic field of Jupiter.

The results show the ability of variating the triangular liberation points by an order from 10−7 to 10−6 for the dimensionless x, y coordinates and distance r from Jupiter. This is equivalent to about hundreds of kilometers which is considerable. The artificial points lie on a circular arc with center at the big primary with radius the unit distance, with angles where is the small deviation from the original point of equilateral triangular case, and the positive sign is for L4 while the negative sign is for L5. It should be mentioned that the results are made for a range of charge per unit mass up to order 0.1 C/kg., where the existing values till now for the specific charge are 0.001–0.01 C/kg. However, it is predicted to overcome the technical problems of this issue to reach values greater than 0.1 C/kg in the future [31] (Peck et al.).

The results also show the stability of the Lorentz artificial L4,5, in the range chosen for the specific charge, as the original points for Sun-Jupiter system.

With the improvement of the studies with artificial liberation points, future studies may generalize the work either to stronger kinds of control or to work with the 5 body problem or even the n-body problem [32, 33].

Abbreviations

SC:Space craft
TBP:Three body problem
RTBP:Restricted three body problem
CRTBP:Circular restricted three body problem
AEP:Artificial equilibrium points
LA:Lorentz acceleration
LF:Lorentz force
(X, Y, Z):The center of mass coordinates system (synodic coordinate system)
(x, y, z):The planet centered coordinates system
x, y, z and r:Dimensionless distances
:Unit vector along the north magnetic pole of the planet
k:Unit in the direction of rotational axis of the planet
:Lorentz acceleration vector
:Charge-to-mass ratio of the SC. in coulombs per kilogram (C/kg)
B:Planet magnetic field vector
Bo:The strength of the field in weber-meters
:The masses of the primaries
:Mass of space craft (SC)
:The inertial spacecraft velocity vector
:The rotational velocity vector of the planet
:The small deviation from the original point
N:Mean motion of the primaries around their center of mass
:Dimensionless angular velocity of the smaller primary.

Data Availability

Data available on request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This article introduced by the students of Prof. Dr. Moustafa Kamal Ahmed group for space dynamics for his memory.