A New Constellation Configuration Scheme for Communicating Architecture in Cislunar Space

A new constellation configuration is proposed to provide near-constant communications between any two points in the cislunar space. Four satellites were allocated on some natural orbits near LL 2 , LL 4 , LL 5 , and LL 3 points, the first three of which are expected to supply a global coverage on the Moon, and the last three afford a global coverage on the Earth. Considering the real dynamics and geometries, a full-ephemeris model is employed to refine the quasi-periodic trajectories near LL 2 , LL 3 , LL 4 , and LL 5 points by a multiple shooting corrector numerically. The revisiting intervals between the successive revisits on specified points by the constellation are demonstrated both on the geographical maps of the Earth and the Moon, with a consideration of equatorial inclinations with respect to the Earth-Moon orbital plane. Moreover, the coverage periodicity on the Moon is induced by Metonic cycle of 19 years. Cost analysis on stationkeeping and launching is implemented for this constellation in detail. Finally, some aspects about communications, such as satellite links, delays, ground stations, and communication devices, are surveyed preliminarily.


Introduction
With the increased interests in lunar exploration and development, autonomous and manned lunar surface missions will become more and more frequent in the future decades.In particular, with the construct of lunar laboratory and lunar base, a large number of data and images need to be transmitted from Moon to Earth (or reverse).Some of the missions, for example, the manned exploration on the lunar surface, require real-time communication.The communicating architecture in cislunar space, therefore, becomes a crucial problem.
Since 1960s many scholars have contributed a lot to the design of satellite configuration architecture.An important task of the constellation design is to make satellites own some desired coverage to the surface of the expected planet.Three satellites placed on the geostationary orbits equidistantly will almost cover the surfaces of the Earth [1].To achieve the multiple coverage, Walker [2] developed a common methodology to design the magical Walker- constellation, which has potential applications in the global position system (GPS).Ballard [3] designed rosette constellations with flowerlike orbital patterns to achieve invariant coverage on different deployment attitudes or orientations relative to Earth coordinates.Most of the significant works were employing Kepler orbit (or perturbed Kepler orbit) around the single planet to provide continuous communication between any points on the Earth.However, the constellation yielded in context of restricted two body problem is weak in building the link between the two separated planets, like Earth and Moon.
Because of the nature of the solutions, some useful architectures for lunar coverage may exist within the context of the restricted three-body problem (R3BP).Farquhar [4,5] used a multibody model to design trajectories for communications between the far side of the Moon and the Earth.As one family of natural solution of R3BP, libration point orbits have potential uses in lunar constellation architecture to ensure adequate coverage.Chow et al. [6] considered that the most cost for the low Earth orbit constellation was the launch costs and reduced cost by using a single large launch vehicle to carry an entire fleet of microsatellites to the  1 point.Grebow [7] created architectures by quasi-periodic orbits around  1 and  2 points for continuous coverage of the south pole of the Moon with only two spacecrafts.Bhasin and Hayden [8] adopted libration point orbits in deep space communicating architecture developed for NASA's future exploration and science missions.Kulkarni et al. [9] postulated a set of satellites placed in Halo orbits at Earth-Moon, Sun-Earth, and Sun-Mars libration points and a couple of relay satellites in a heliocentric earth orbit to support data and command communication for future space missions to Moon, Mars, and beyond.Lee et al. [10,11] proposed a conceptual constellation to substitute geostationary orbit by three satellites placed on  3 ,  4 , and  5 points; however, they did give the feasibility and coverage analysis in detail.Grebow et al. [12] investigated some natural orbits in full ephemeris model with periods ranging from 7 to 16 days and employed two spacecrafts in combination with Earth-Moon libration point orbits to achieve constant communications for lunar south pole exploration mission.
Quite different from the researches referred above, a new constellation configuration is proposed to provide nearconstant communications between any two points in the cislunar space in this paper, with four satellites allocated on some natural orbits near  2 ,  4 ,  5 , and  3 points, respectively.The first three satellites in this communicating architecture are expected to supply a global coverage on Moon, while the last three will afford a global coverage on Earth.The constellation is beneficial for the communications between different points on the grounds of Earth or Moon, even for the case with no direct links as one is located at the far side of the Moon and another is placed on the Earth.
Considering the real dynamics and geometries, a full ephemeris model is employed to refine the quasi-periodic trajectories near  2 ,  3 ,  4 , and  5 points by a multiple shooting corrector numerically.The revisiting intervals between the successive revisits on specified points by the constellation are demonstrated on the geographical maps of both the Earth and the Moon, with a consideration of equatorial inclinations with respect to the Earth-Moon orbital plane.Moreover, the coverage periodicity on the Moon is induced by Metonic cycle of 19 years.
Furthermore, cost analysis on stationkeeping and launching is implemented in detail in full ephemeris model.According to the stationkeeping and inserting strategies proposed in this paper, the total inserting cost of this constellation is 14.4731 km/s, and the total stationkeeping cost of this constellation is 261.8 m/s/year.
Compared with two geostationary orbit constellations that require six spacecrafts to cover both the Earth and Moon, this constellation architecture saves two spacecrafts to reduce the cost of the system a lot, which is considered as an alternative of limited geostationary orbit resources.

Constellation Configuration Scheme.
To describe the global coverage constellation in cislunar space, four different coordinates are referred to in this paper, as the syzygy frame in Earth-Moon system, the Earth-centered fixed frame, (i.e., WGS-84 frame), the Moon-centered fixed frame, and the Earth-centered inertial frame (i.e., J2000 frame), where the syzygy frame has the following features: the origin O is fixed at the barycenter of the Earth-Moon system; the z-axis is perpendicular to the lunar orbital plane; the x-axis points from Earth to Moon; and the y-axis is determined by the right-hand side rule.
Ignored other celestial bodies' gravitation and eccentricity of lunar orbit, the dynamics of system consisted of spacecraft, Earth, and Moon can be regarded as the circular restricted three-body problem (CR3BP).There are five libration points fixed in the syzygy frame, labeled as   ,  = 1, . . ., 5, where  2 and  3 locate at the far side of the Moon and Earth, respectively, and  4 and  5 form an equilateral triangle with Earth and Moon.
The combination of  2 ,  4 ,  5 , and  3 shapes a double-triangle configuration, where the first three points product one triangle focusing at the Moon, while the last three form the other triangle centering the Earth.By the inspiration of the global coverage constellation created by three satellites equispaced on the geostationary orbit, a new communicating architecture is proposed to provide some links (with time delay of 2-3 sec) between any two points in the cislunar space by  2 ,  4 ,  5 , and  3 points, shown in Figure 1.With the help of this constellation system, a direct or indirect link exists for any two points in cislunar space, even for the two points   and   located on the far side of Earth and Moon respectively (one of the feasible indirect links is as follows:   →  2 →  4 →  3 →   ).Moreover, the quasi-halo orbit near  2 inheriting from halo family is employed to avoid some links from  2 point to Earth blocked by the Moon.

Numerical Refinements by Multiple Shooting Corrector.
The eccentricity of the Moon's orbit and gravitational forces from other celestial bodies (like the Sun and the Jupiter) will make CR3BP different from the real dynamics, even the libration point cannot keep stationary due to this perturbations.
Therefore, a multiple shooting corrector is employed to create some natural orbits near the libration points in the full ephemeris model, which is useful for the computation of highly unstable periodic orbits with very long periods, as  2 invariant relative orbits [13,14] and solar sail's displaced orbits [15].
The general idea is to use a multiple shooting method similar to the one used for the numerical solution of boundary value problems.As in the iteration procedure, the total time span [ 0 ,   ] is split into a number of shorter subintervals by  equispaced points  0 ,  1 ,  2 , . . .,  −1 ,   .Denote the subinterval with equal length as Δ =  +1 −   ( = 0, 1, 2, . . .,  − 1), and denote the point on the quasi-periodic trajectory as the dynamics, some corrections to Q  are required to fulfill the matching conditions, as ) − ( In the th iteration of this procedure, the matching equation can be written as where .The differential function of F has the following expression: ) .
The iterations have been included in the software package OrbitalPredictor developed by Beihang University [16], which adapts JPL ephemeris model to generate some natural orbits near the libration points.
To accelerate the convergence of the algorithm, an initial iteration to the natural solution of the dynamics is required.For  3 ,  4 , and  5 , the libration points themselves can work as the initial iterations to create bounded trajectories near these points; however, for a quasi-halo orbit near  2 point, an effective initial iteration can employ the following Richardson expansion developed from the Lindstedt-Poincaré method [17,18]: where the detailed discussion on the coefficients , , , , and  can be found in [17].To keep the frequencies in -axis equal to and -axes, the amplitudes   and   are dependent on each other and also need to be no less than some critical values.Therefore, the amplitude   can be selected as the only parameter to identify different halo orbits, and then a halo orbit with the minimal amplitude   = 11500 km is created by the iteration algorithm in this paper.Following the initial iterations proposed for  2 ,  3 ,  4 , and  5 referred to above, some bounded trajectories can be refined by the OrbitalPredictor full-ephemeris model shown in Figure 2 and then are targeted in Satellite Tool Kit (STK) using Astrogator tool [12,19], which are shown in Figure 3 by STK.

Coverage Analysis of Constellation
The revisiting time is employed to parameterize the coverage of this constellation architecture, which is defined as the time intervals between the successive revisits on the specified ground point on the Earth or Moon.Therefore, this concept is measuring the interruption durations when communicating with any point.x (km) x (km) x (km) x (km) x ( k m )

Coverage to Moon by the Constellation.
As the preliminary investigation on the coverage to Moon, some hypotheses can be adopted to plot the covering region, including the following: the inclination of the equator plane with respect to the lunar orbital plane is neglected as 0 ∘ , instead of its true value of 6 ∘ 41  ; the eccentricity of lunar orbit and the gravitation from the Sun are ignored to place the spacecrafts on  3 ,  4 , and  5 points and the halo orbit near  2 without any stationkeeping maneuver; that is, the libration point orbits are yielded from CR3BP.Due to the osculating position on halo orbit, the maximal and common regions are achieved in Figure 4, where the former is defined as the region covered by the constellation at any time, while the latter is defined as the region covered all the time.Therefore, the constellation can cover the lunar surface by 99.45% (maximal region) and 98.95% (common region) except the south and north poles.
In a full investigation on the coverage to Moon, no hypothesis is employed, including the real geometries and the full dynamics.All the bounded trajectories in this constellation are refined by the OrbitalPredictor full-ephemeris model, and the revisiting intervals are counted statistically in real geometrical scenario by STK.
The Moon is in synchronous rotation: it rotates about its axis in about the same time it takes to orbit the Earth; however, due to the inclination between lunar orbital and equator planes and some others, there are several librations to change the angle from which the Moon is seen, allowing about 59% of its surface to be seen from the Earth.
Moreover, the 19-year Metonic cycle of the equatorial inclination, with respect to lunar orbital plane, makes the Moon go overhead at latitudes over the equator from 18.5 ∘ to 28.5 ∘ , so it is necessary to give results in term of time in a whole cycle.Figure 5 depicts the revisiting times on the lunar surface by the constellations created from different epoch times in a Metonic cycle.The maximal coverage for full geometries is investigated by the maximal revisiting times during this cycle on the specified point; however, the common coverage is yielded by the minimal revisiting times during this cycle, and then both the maximal and minimal revisiting times are revealed in Figure 6 on the global surface.Different from the preliminary results, some conclusions can be induced from Figures 5 and 6.Firstly, all the points on the lunar surface can be covered by this constellation, where most of them are accessed without any time interval (i.e., the blank regions in the top of Figure 6), and few of them are revisited after some intervals less than 5 days (i.e., the colored regions in the top of Figure 6).
Secondly, due to the lunar synchronous rotation, this constellation is weak in covering some longitude lines located near 150 ∘ W and 150 ∘ E shown in Figure 5.The discontinuous communications are interrupted by some intervals less than 2 hours.
Thirdly, only some sites keep discontinuous communications during a whole Metonic cycle, which are marked in the bottom of Figure 6 from 20 ∘ N140 ∘ W to 70 ∘ N140 ∘ W, from 30 ∘ S150 ∘ W to 2 ∘ N150 ∘ W, from 85 ∘ S55 ∘ E to 85 ∘ S100 ∘ E, and from 80 ∘ S170 ∘ E to 80 ∘ S130 ∘ W. Fourthly, the worst intervals occur at the high-latitude areas instead of the poles because the maximal revisiting time is 0.96 day at the north pole while 1.21 days at the south pole.Thus this constellation has a good coverage performance at the poles.

Coverage to Earth by the Constellation.
As the preliminary investigation on the coverage, the irrotational Earth is assumed in combination with no consideration of the eccentricity of lunar orbit and the gravitation from the Sun.Therefore, the constellation has the fixed configuration in the syzygy frame, which results in the maximal coverage region equal to the common coverage region shown in Figure 7. From the hypothesis referred above, the constellation can cover the Earth's surface by 99.1% except the south and north poles.
In a full investigation on the coverage to Earth, the hypothesis is unloaded with some considerations of the real geometries and full dynamics.All the quasi-periodic trajectories near  3 ,  4 , and  5 are refined by the OrbitalPredictor full-ephemeris model, and the revisiting intervals are counted statistically in real geometrical scenario by STK.
The 19-year Metonic cycle of the equatorial inclination, with respect to lunar orbital plane, has also an influence on the coverage by the constellation.However, because of rotating more quickly around its own axis, the Earth owns the constant revisiting time on the same latitude lines located near 70 ∘ S and 70 ∘ N shown in Figure 8, where the discontinuous communications are interrupted by some intervals less than 2.5 hours.It is quite different from the Moon where the discontinuous communications are just locating on the longitude lines near 150 ∘ W and 150 ∘ E.Moreover, the regions on other latitudes are accessed without any time interval (i.e., the blank regions in the top of Figure 8).
In addition, the conclusion that the revisiting time keeps invariant on the same latitude line will be true for a whole Metonic cycle.Therefore, the relationship between revisiting time and the latitudes will be investigated in the following section just on the longitude line of 0 ∘ E, shown in Figures 9  and 10.
Different from the preliminary results, some conclusions can be induced from Figures 9 and 10.Firstly, all the points on the Earth's surface can be covered by this constellation, and the discontinuous communications keep the invariant revisiting time on the same latitudes during a whole Metonic cycle.
Secondly, due to the rotation of the Earth, few of the regions located at the latitudes from 80 ∘ S to 52 ∘ S and from 55 ∘ N to 79 ∘ N are revisited after some intervals less than 2.5 hours, while the other regions are accessed without any time interval.
Thirdly, the worst intervals occur at the regions located from 71 ∘ S to 65 ∘ S and from 62 ∘ N to 70 ∘ N, which keep discontinuous communications during a whole Metonic cycle.Moreover, this constellation has a good coverage performance at the poles with no time interval.Therefore, two communication devices need to be equipped on the spacecrafts, with their antennas pointing to the Earth and Moon, respectively.However, due to the Metonic cycle of the equatorial inclination with respect to lunar orbital plane, the boundaries of the angle between the two antennas are illustrated points in Figure 11, where the angle between the two antennas near  4 point varies from 47 ∘ to 81 ∘ , and from 34 ∘ to 85 ∘ near  5 point.

Stationkeeping Strategies and Costs.
Stationkeeping correction is necessary because the true trajectory will diverge from the nominal one under the effects of the navigation error  1 , the initial error from the launch rocket  2 , and the error from inserting the nominal trajectory on velocity  3 [20,21].However, the presence of the trajectory correction will cause velocity errors  4 .These error resources are assumed to be Gauss distributed and their average values are zero. 1 and  2 have the same position error variances of 10 km and velocity error variances of 0.01 m/s (1);  3 and  4 have the velocity error variances of 3% and their velocity errors in the direction of radius are 2.06 ∘ (3) [22,23].The velocity error in the normal direction is uniformly distributed.Assume that there are  corrections during the flight on nominal trajectory, which is divided into  + 1 sections according to the amount of corrections, as it is shown in Figure 12.
The state errors between the actual and nominal trajectories at time   are R  and V  from the errors  1 ,  2 ,  3 , and  4 , respectively, and the Taylor expansion of the position vector at time   is obtained where the correction velocity ΔV  is used to correct the trajectory aiming at achieving the wanted position R( +1 , R  , V  ,   ) at time   .According to (5), ΔV  can be obtained The correction velocity ΔV  is a linear comparison of R  and V  , and this characteristic is definitely suitable for trajectory correction during the flight.
A correction schedule with an equi-interval correction time of 7 days is adopted in this paper to investigate the stationkeeping costs for this constellation in the full ephemeris model.The maneuver results are obtained in Table 1 for one complete year.Therefore, the total stationkeeping cost of this constellation is 261.8 m/s/year, which is accepted for the aerospace engineering.
Furthermore, some optimization technique to decrease the stationkeeping cost can be implemented in the future, including some selections for the correction schedule, and the application of the stochastic optimal control into the correction strategy.

Inserting Strategies and Costs into Nominal Trajectories.
After being launching into the parking orbit around the Earth, some maneuvers are required to guidance the spacecrafts to nominal trajectories located near  2 ,  3 ,  4 , and  5 , respectively, which are refined in Section 2.2.Some inserting strategies are included in this paper to analyze the inserting costs for this constellation.
The invariant manifold theory developed in CR3BP has potential applications in creating the inserting trajectories into a halo orbit near  1 point [24,25].However, According to the technique developed by Koon et al. [26], the insertion into the quasi-periodic orbit near  2 point needs some combinations between unstable manifolds stemming from  1 (or  2 ) point and stable manifolds stemming from  2 point.In fact, the similar results can be achieved by the weak stability boundary theory (WSB) proposed by Belbruno and Miller [27] to create some low-energy cislunar trajectory with the help of the Sun's gravitation.Therefore, following on the initial guess obtained in CR3BP, the iteration procedure defined in Section 2.2 can be used to product the real inserting trajectory to the quasi-periodic orbit near  2 point in the full ephemeris model, shown in Figure 13.Only one maneuver is required on the parking orbit with its altitude of 200 km to insert the nominal orbit, without other   maneuver at the end of the transfer trajectory.Thus, during this transfer, the total transfer duration is 91.3 days, and the total cost of maneuvers is 3.2018 km/s.For the bounded trajectories near  3 ,  4 , and  5 , all of the stable manifolds stemming from these libration points are quite far from the Earth, which is different from  2 point.The Hohmann transfer is employed in this paper to guide the spacecrafts from the parking orbits to their stable manifolds, and then to insert the nominal trajectory following its stable manifolds.Based on the initial guess obtained in CR3BP, the iteration procedure defined in Section 2.2 can be used to product the real inserting trajectory to the bounded orbit near  3 ,  4 , and  5 points in the full ephemeris model, shown in Figures 14, 15, and 16, respectively.Therefore, the whole transfer trajectory includes the Hohmann segment and coasting segment along the stable manifold, and two maneuvers are required on the parking orbit with its altitude of 200 km to insert the nominal orbit; the first one is imposed on the parking orbit, and the last one is imposed at the end of the Hohmann segment.
For  3 point, the total transfer duration is 45.1 days, and the total cost of maneuvers is 3.7518 km/s.For  4 point, the total transfer duration is 47.57days, and the total cost of maneuvers is 3.7585 km/s.For  5 point, the total transfer duration is 63.4 days, and the total cost of maneuvers is 3.761 km/s.Therefore, the total inserting cost of this constellation is 14.4731 km/s, which is accepted for the aerospace engineering.Furthermore, the inserting costs by Hohmann segments in the cases of  3 ,  4 , and  5 points can be optimized by some techniques developed in space manifold dynamics in the future.

Communications
There are several communicating links inside the constellation, which is shown in Table 2.The links make the constellation more functional than the single sum of the four satellites.Furthermore, the links make the constellation more robust, so the constellation will still work depending on the links at the event that one or two satellites happen to be broken.
The distance from Earth to Moon brings in a long time delay of 1.28 s.Due to this, direct voice communications may be impossible.Instead, data communications are preferred for this constellation architecture, including some links between ground sites on the Earth and the far side on the Moon or some data exchanges between the points on the Earth, and so on.Furthermore, some ground stations are required to track the four spacecrafts in the constellation, and also to deal with orbital determination.Due to the data rates and Bands, the diameter of antenna is not less than 50 m.According to the coverage region defined in Figure 8, the ground stations should be located at the latitudes from 80 ∘ S to 52 ∘ S and from 55 ∘ N to 79 ∘ N to track and manage this constellation.

Conclusion
To provide near-constant communications between any two points in the cislunar space, a new constellation configuration is proposed in this paper.Four satellites were allocated on some natural orbits near  2 ,  4 ,  5 , and  3 points, the first three of which are expected to supply a global coverage on the Moon, while the last three afford a global coverage on the Earth.Considering the real dynamics and geometries, a full ephemeris model is employed to refine the quasi-periodic trajectories near  2 ,  3 ,  4 , and  5 points by a multiple shooting corrector numerically.
The revisiting intervals between the successive revisits are demonstrated on the geographical maps of the Earth and the Moon, with a consideration of the Metonic cycle of the equatorial inclination with respect to lunar orbital plane.All the points on the lunar surface can be covered by this constellation, where most of them are accessed without any time interval, and few of them are revisited after some intervals less than 5 days.This constellation is weak in covering some longitude lines located near 150 ∘ W and 150 ∘ E. The discontinuous communications are interrupted by some intervals less than 2 hours.Only some sites keep discontinuous communications during a whole Metonic cycle, which are located from 20 ∘ N140 ∘ W to 70 ∘ N140 ∘ W, from 30 ∘ S150 ∘ W to 2 ∘ N150 ∘ W, from 85 ∘ S55 ∘ E to 85 ∘ S100 ∘ E, and from 80 ∘ S170 ∘ E to 80 ∘ S130 ∘ W. The worst intervals occur at the high-latitude areas instead of the poles because the maximal revisiting time is 0.96 day at the north pole while 1.21 days at the south pole.This constellation has a good coverage performance at the lunar poles.
All the points on the Earth's surface can be covered by this constellation, and the discontinuous communications keep the invariant revisiting time on the same latitudes during a whole Metonic cycle.Few of regions located at the latitudes from 80 ∘ S to 52 ∘ S and from 55 ∘ N to 79 ∘ N are revisited after some intervals less than 2.5 hours, while the other regions are accessed without any time interval.The worst intervals occur at the regions located from 71 ∘ S to 65 ∘ S and from 62 ∘ N to 70 ∘ N, which keep discontinuous communications during a whole Metonic cycle.This constellation has a good coverage performance at the Earth's poles with no time interval.
Moreover, cost analysis on stationkeeping and launching is implemented in detail in full ephemeris model.According to the stationkeeping and inserting strategies proposed in this paper, the total inserting cost of this constellation is 14.4731 km/s, and the total stationkeeping cost of this constellation is 261.8 m/s/year.Furthermore, some optimizations on this constellation can be implemented in the future, including optimizing the orbital elements in this constellation to promote the coverage on the Earth and Moon, selecting the correction schedule from viewpoint of the stochastic optimal control to decrease the stationkeeping cost, and optimizing the Hohmann segments for  3 ,  4 , and  5 points to decrease the launching cost based on space manifold dynamics.

Figure 3 :Figure 4 :
Figure 3: The quasi-periodic trajectories near  2 ,  3 ,  4 , and  5 created by STK with the epoch time of 12:00, July 1, 2015.The bounded trajectories are shown in (a) the syzygy frame in Earth-Moon system and (b) the Earth-centered inertial frame.The ground tracks of the trajectories are shown in (c) the Earth-centered fixed frames and (d) the Moon-centered fixed frames.

Figure 6 :Figure 7 :
Figure 6: The full analysis on the maximal and common coverages during a whole Metonic cycle from 2015 to 2034.The maximal revisiting times on the global surface during a whole Metonic cycle are implemented in (a), and the minimal revisiting times in (b).

Figure 8 :
Figure 8: The full analysis on revisiting time on the Earth's surface with the epoch time of 12:00, July 1, 2015.The conclusion that the revisiting time keeps invariant on the same latitude line will be true for a whole Metonic cycle.
4 and  5 .In this constellation, the  4 point takes simultaneously charge in the communications with Earth and Moon, as well as  5 .

Figure 10 :
Figure 10: The relationship between the latitudes and the maximal and common revisiting times during a whole Metonic cycle from 2015 to 2034.The maximal revisiting times on the latitudes during a whole Metonic cycle are implemented in (a), and the minimal revisiting times in (b).

Figure 11 :
Figure 11:  The angle between the two antennas equipped on the  4 and  5 spacecrafts.The maximal angle during a Metonic cycle is indicated by the upper boundary, while the minimal angle during a Metonic cycle is indicated by the lower boundary.

Figure 13 :Figure 14 :
Figure 13: The inserting trajectory into the nominal orbit near  2 point from the parking orbit with its altitude of 200 km in the syzygy frame.The quasi-periodic halo orbit near  2 point created in full ephemeris model with the epoch time of 12:00, July 1, 2015.Only one maneuver Δ = 3.2018 km/s on the parking orbit is required to insert the nominal orbit without other maneuvers at the end of the transfer trajectory.The transfer duration is 91.3 days.
,   ,   , ẋ  , ẏ  , ż  ]  ( = 0, 1, 2, . . .,  − 1).Let (Q  ) be the image of the point Q  under the flow associated to the full ephemeris dynamics after an amount of time Δ.To keep all the points Q  staying on the same trajectory solved from Figure 1: The communicating constellation architecture in cislunar space. 2 ,  4 , and  5 form one triangle to supply a global coverage on Moon;  3 ,  4 , and  5 form the other triangle points to supply a global coverage on Earth.