<<

Mission Design for Deep Space Nano/Micro Utilizing Lunar Orbital Platform-Gateway Opportunities

Naoya Ozaki Institute of Space and Astronautical Science, Aerospace Exploration Agency

2020/10/10 1 2 Deep Space

©Morehead ©NASA State ©University of Tokyo University /JAXA /NASA

©NASA

©JAXA

©JPL/NASA

3 ©NASA

Lunar Orbital Platform-Gateway (In ) 4 LOP-G Related Launch Opportunities Starting from NASA’s -1, we can expect more than 10 CubeSats are launched to deep space every year. (Launch for LOP-G Construction, Resupply, etc…) ©NASA

NASA Updates Plans, NASA .com (Accessed on March 16, 2019) https://www.nasaspaceflight.com/2018/09/nasa-lunar-gateway-plans/ 5 A New World Opened by LOP-G and Nano/Micro Spacecraft LOP-G Innovation in (Low cost, short lead time) Similar innovation • Explosion in numbers • Frequent missions will happen in • Expansion of stakeholders deep space (Startups, universities, etc) missions!!

Number of Satellite Launched For A Year 6 7 http://unisec.jp/unisec/satellites.html(Info. in 2019,accessed on July 11, 2020) What is Lunar Orbital Platform-Gateway (LOP-G)?? Lunar Orbital Platform-Gateway (LOP-G) is a planned in lunar . NASA’s plays a major role to develop the Gateway in collaboration with commercial and international partners: ESA, JAXA, CSA, , etc.

8 Where is the Gateway??

Near Rectilinear

A type of Halo under the Earth and gravity.

9 Where is the Gateway??

Centrifugal force

Moon’s gravity

Earth’s gravity

The equilibrium point where the earth's gravity, the moon's gravity, and the centrifugal force balance each other is called the point. 10 Lagrange Points

Five types of Lagrange points exist in each three-body dynamical system. Ex) Earth-Moon L2 -Earth L1 Lagrange point 11 Geometry of Lagrange Points

Earth-Moon L3

Sun-Earth L1 Sun-Earth L2 Sun L4 L1 L5 L2

Sun-Earth line fixed rotational frame.

12 Geometry of Lagrange Points

Earth-Moon L3

Sun-Earth L1 Sun-Earth L2 Sun L4 L1 L5

L2 Gateway is here!!

Sun-Earth line fixed rotational frame.

13 Near Rectilinear Halo Orbit

14 the orbit, and over time the spacecraft drifts from the baseline. With sufficient drift, the spacecraft is at risk f lng eclie b he Eah had. The cen ineigaion explores phase control within the NRHO to maintain the eclipse-free characteristics of the virtual reference. The x-axis crossing control method is augmented to maintain periapse passage time, thus maintaining the eclipse-free phase achieved in the baseline NRHO. Then, a rendezvous strategy is developed to target back to the long-horizon NRHO in cases where the orbit evolves over time away from the virtual reference in the presence of large perturbations.

15-YEAR BASELINE NRHO: OSCULATING PARAMETERS AND ECLIP S ES

Historically, halo orbit missions including WIND7 and ARTEMIS8 have operated without a reference trajectory. Halo orbit stationkeeping is effective and inexpensive without targeting parameters from a pre- defined reference. However, as the L2 halo family approaches the Moon, the costs and computation time Nearassociated Rectilinear with orbit maintenance are decreased Halo by employing Orbit a baseline (NRHO) trajectory as a virtual reference, that is, as a catalog of targeting parameters.9 Adhering strictly to a reference orbit is unnecessarily expensive. CharacteristicsInstead, by targeting specific of NRHO parameters extracted from a virtual reference trajectory, a spacecraft can • maintain the orbit for low propellant costs while retaining characteristics of the reference. For the Gateway, Geometricremaining nea he relationship efeence i iman tof aidingthe Earth lng eclie is alwaysfm he Eah the had same a ell a f • Halofacilitating orbits mission are design unstable for spacecraft (easility visiting the reachable/escapableGateway, including , lunar ) and elements, spacecraftlogistics modules, cannot and others. stay in the orbits without station keeping maneuvers,The current 15 -butyear baselineNRHO orbit is for less the Gateway stable spacecraft than is designedgeneral by Lee Halo3 in an ephemerisorbits. model that includes n-body gravitational forces from the Sun, Earth, Moon, and barycenter. The Moon is • Gatewaymodeled with will an 8x8be gravityconstructed field, while thein other9:2 threesynodic bodie s areresonant considered pointNRHO masses. Non- (SRHO),gravitational whereforces, including the solarstation radiation never pressure (SRP), experiences are not included ineclipse the force modeling.. The • Fororbit 9:2extends SRHO, from January the 2020 perilune to February 2035altitude, and other is than 1458 small discontinuities-1820km, in velocityand theto maintain the almost-stable orbit (averaging less than 1.9 mm/s per revolution), the NRHO is a ballistic trajectory. The apolunefull 15-year ephemerisaltitude is plotted is 68267 in Figure -170112km.a in an Earth-Moon rotating view and in Figure 1b in a Sun-Earth rotating view.

Figure 1. 15-year reference NRHO in Earth-Mo on (a) and S un-Earth (b) rotating views 15 (D.C. Davis, F.S.Khoury, et al. , AAS, 2020) Over the 15-year span, the mean (time from one perilune passage to the next) of the Gateway NRHO ranges from 6.26 days to 6.76 days, with a mean value of 6.56 days. Similarly, the mean perilune radius is 3,366 km with a minimum value of 3,195 km and a maximum value of 3,557 km. The apolune radius can be as large as 71,849 km, or as small as 70,005 km, with an average value of 71,100 km. Osculating parameters are summarized in Table 1. Further details on the generation and characteristics of the baseline NRHO appear in a white paper.3

Table 1. Osculating Gateway orbital parameters over 15 years Orbital Perilune Perilune Apolune Apolune period (days) radius (km) altitude (km) radius (km) altitude (km) Minimum 6.26 3,195 1,458 70,005 68,267 Mean 6.56 3,366 1,629 71,100 69,363 Maximum 6.76 3,557 1,820 71,849 70,112

2

Launch Condition for Gateway Construction For gateway construction opportunities, the spacecraft is expected to be launched into lunar transfer trajectory (as was the case with ). Using lunar swing-by on this orbit, the spacecraft can fly to interplanetary space (to asteroid, , …), periodic orbits in Lagrange points, and so on. (a) L1 Single Rev (7 days). (b) L2 Short Stay (3-4 days). (c) L2 Long Stay (9.5-11 days). Earth

Moon

OPF RPF

NRHO (L2)

NRI NRD NRHO Stay (d) Full(J. Trajectory Williams, for L2 Short D.E. Stay. Lee, et al. , AAS, 2017) 16 Figure 10: Roundtrip Missions to Earth-Moon NRHOs(rp = 4500 km, South Families) Shown in the Earth-Moon Rotating Frame. For each mission type, Orion must perform four major trans- lational maneuvers (OPF, NRI, NRD, and RPF). The L2 cases include both “short stay” and “long stay” mission types. Stay times in the NRHO are shown in parentheses.

12 Mission Utilizing Gateway Construction Opportunity Launch & Separation Lunar transfer Lunar orbits, surface, trajectory Earth-Moon(EM) L1/L2, Sun-Earth(SE) L1/L2, ~1 week Interplanetary (Mars, Lunar swing-by asteroids,…) Trajectory (Vinf~1km/s) correction ΔV ~20m/s

Possible mission scenario for early Artemis opportunity (such as )

17 Mission Utilizing Gateway Construction Opportunity Launch & Separation Lunar transfer Lunar orbits, surface, trajectory Earth-Moon(EM) L1/L2, Sun-Earth(SE) L1/L2, ~1 week Interplanetary (Mars, Lunar swing-by asteroids,…) Trajectory (Vinf~1km/s) correction ΔV ~20m/s

1) To Sun-Earth L1/L2: Ballistic transfer without ΔV

2) To Lunar Surface:

3) To Lunar orbits, Earth-Moon L1/L2:

4) To Interplanetary (Mars, asteroids):

18 1) Transfer to Sun-Earth L1/L2 Points

There is a manifold structure in periodic orbits around Lagrange points. When a small ΔV disturb at each difference phase on the periodic orbit, the spacecraft leaves the periodic orbits on a group of orbits called the unstable manifold.

Because of the symmetry, the spacecraft can asymptoticaly approach a period orbit by riding on a group of orbits called the stable manifold.

(W.S., Koon, M.W., Lo, et al. , AAS, 2000)

19 1) Transfer to Sun-Earth L1/L2 Points

By selecting a single trajectory on a stable manifold, the spacecraft can ballistically tranfer to the periodic orbit around the LagrangeLOW-ENERGY MISSION DESIGN point. 97

(J.S.Parker and R.L.Anderson, JPL, p.97, 2013) 20

Figure 2-37 Anexampleunstablehaloorbit(green)abouttheSun–EarthL2 pointandits stableinvariantmanifold(blue).(Seeinsertforcolorrepresentationofthisfigure.)

barsindicatethelocationsinthemanifoldthathaveperigeealtitudesbelow500km. Forexample,onecanseethatthetrajectorywitha⌧ -valueof0.751encountersa closestapproachwiththeEarthwithaperigeealtitudeofapproximately185kmand aneclipticinclinationofapproximately34.8deg. Hence,aspacecraftinacircular lowEarthparkingorbitwithanaltitudeof185kmandaneclipticinclinationof 34.8degmayperformatangentialVtotransferontothemanifold; onceonthe manifold,thespacecraftballisticallyfollowsitandasymptoticallyarrivesonthehalo orbit. Therearetwostatementsintheseresultsthatneedtobeaddressed. Thefirstis thattheinclinationvaluesdisplayedinFig.2-38aretheinclinationvaluescomputed intheaxesoftheCRTBP:namely, inaplanethatisverysimilartotheecliptic. Theequatorialinclinationvaluesoftheseperigeepointsdependonwhichdatea spacecraft launches. Since the Earth’s rotational axis is tilted by approximately 23.45degwithrespecttotheecliptic[97],manyequatorialinclinationsmaybeused to inject onto a desired trajectory, depending on the date. The second statement thatshouldbeaddressedisthattheresultsshowninFig.2-38dependgreatlyon theperturbationmagnitude,✏,describedinSection2.6.10.Implementingadifferent perturbationmagnituderesultsinachangeinthe⌧-valuesrequiredtoobtainacertain trajectory. For example, if ✏ were reduced, the trajectories modeling the orbit’s manifold would spend more time asymptotically approaching/departing the orbit. Oncethetrajectoriesaresufficientlyfarfromtheorbit,theircharacteristicsarenearly Mission Utilizing Gateway Construction Opportunity Launch & Separation Lunar transfer Lunar orbits, surface, trajectory Earth-Moon(EM) L1/L2, Sun-Earth(SE) L1/L2, ~1 week Interplanetary (Mars, Lunar swing-by asteroids,…) Trajectory (Vinf~1km/s) correction ΔV ~20m/s

1) To Sun-Earth L1/L2: Ballistic transfer without ΔV

2) To Lunar Surface: Landing with 2.5km/s ΔV

3) To Lunar orbits, Earth-Moon L1/L2:

4) To Interplanetary (Mars, asteroids):

21 Mission Analysis for the EM-1 CubeSats EQUULEUS and OMOTENASHI

2) Landing on Lunar Surface OMOTENASHI Trajectory Figure 5. Launch & Lunar night portion and altitude for the science orbit. Separation Lunarenvironment transfer and soil mechanics to reduce the risks of Lunar orbits, surface, futuretrajectory human exploration. The spacecraft is composed by a small Surface Probe Earth-Moon(EM) L1/L2, with an airbag and a shock absorption mechanism; a Sun-Earth(SE) L1/L2, Retromotor Module with the solid to slow the Interplanetary (Mars, spacecraft~1 week at lunar approach; and the Orbit Module to guide Lunar landing asteroids,…) Trajectorythe spacecraft from SLS separation(Vinf~1km/s) to the solid rocket correctionignition. ΔV 20m/sOMOTENASHI introduces a unique approach for (S. Campagnola, et al., IEEE, 2019) ~lunar landing. Traditionally, lunar landing missions are Figure 6. characterized by an Earth–Moon transfer and OMOTENASHI trajectory: Approach and landing phases in the It is possibleinsertion, to estimate followed by the descent, the hovering,landing and landing ΔV byMoon assuming body-fixed frame (from a Hernando-Ayusotwo-body [17]). problem (patchedphases. This approach conics) allows for a flexible design, as the errors in all phases can be detected and corrected, but before impact, the solid motor executes a maneuver of In the vicinityrequires of athe full set moon, of sensors andvis a- largeviva restartable equation propul- (orbitalabout 2500-energy m/s to bring-invariance the spacecraft almost law) to a stop, gives sion system, both of which are not available to small followed by a free fall onto the lunar surface. The semi- . OMOTENASHI combines1 the maneuvers%& for 1 hard landing is enabled by a shock absorption mechanism ! " ! the lunar orbit insertion, descent, and landing# − into a single= # that allows up to 30 m/s impact vertical velocity, corre- 2 ' 2 # maneuver# = executed0.82 bykm a solid/s rocket motor, followed by a sponding to a few' hundred= ' meters of free fall. Suppose thatfree-fall # onto the lunar surface(example with impact speed of on Artemis the The1) and nominal trajectory"( islunar designed radius), to maximize te velocity isorder of 30 m/s. If proven to work, the OMOTENASHI robustness. With the current assumptions on knowledge approach will enable an entirely new# = class2.51 of lunar4 km/s and execution error, only 60% of the Monte Carlo runs exploration missions by small satellites, also exploiting successfully land with an impact velocity below 30 m/s. In order to landmore ride-shareon the opportunities moon, from we the need planned to Lunar cancelThe this major sourcevelocity, of failure i.e., are the errors on knowledge, ΔV~2.5km/s.Orbital Platform-Gateway [14]. thrust direction, and thrust duration. Following this analy- sis, the project is now considering solutions to improve the knowledge using Delta-DOR, and to increase the22 max- MISSION ANALYSIS imum impact velocity to 50 m/s. This paper provides a quick overview of the mission analysis and design approach. More details can be found DESIGN APPROACH in our cited papers [15], [16]. Figure 6 shows the current nominal trajectory and a close-up of the Lunar approach. OMOTENASHI trajectory design is unusual, as the nomi- One day after separation, Dv1 is executed by a cold-gas jet nal trajectory design is driven by knowledge and execution system to correct for launcher dispersions errors and to tar- errors. Therefore, multiple design iterations are performed get a lunar landing orbit, with shallow flight-path-angle at with different models. approach on a landing site visible from the Earth. Depend- The sensitivities to the errors at landing are first investi- ing on the launch geometry, Dv1 can be 5 to 16 m/s; a gated in a simplified, flat-Moon free-fall model, for different trajectory correction maneuver is also planned to correct approach angles (FPA) and post-Dv2 heights (h). Knowl- for Dv1 execution errors and knowledge errors. Shortly edge errors (along track and radial) and attitude pointing

42 IEEE A&E SYSTEMS MAGAZINE APRIL 2019 Mission Utilizing Gateway Construction Opportunity Launch & Separation Lunar transfer Lunar orbits, surface, trajectory Earth-Moon(EM) L1/L2, Sun-Earth(SE) L1/L2, ~1 week Interplanetary (Mars, Lunar swing-by asteroids,…) Trajectory (Vinf~1km/s) correction ΔV ~20m/s

1) To Sun-Earth L1/L2: Ballistic transfer without ΔV

2) To Lunar Surface: Landing with 2.5km/s ΔV

3) To Lunar orbits, Earth-Moon L1/L2: Direct insertion or Low-energy transfer/capture by reducing Vinf.

4) To Interplanetary (Mars, asteroids):

23 3) Transfer to Lunar Orbits or Earth-Moon L1/L2 points Launch & Separation Lunar transfer trajectory Low-energy transfer

~1 week Lunar swing-by 5 to 8 months Trajectory (Vinf~1km/s) Low-energy correction ΔV or capture ~20m/s Direct insertion

For insertion, it is possible to estimate the insertion ΔV by assuming a two-body problem (patched conics). The ΔV is about 0.5-1km/s. For low-energy transfer/capture, the tidal force can effectively reduce Vinf. For NRHO or other Halo orbits, only about 10m/s ΔV is required for the insertion.

24 3) Transfer to Lunar Orbits or Earth-Moon L1/L2 points Launch & Separation Lunar transfer trajectory Low-energy transfer

~1 week 5 to 8 months Trajectory Lunar swing-by correction ΔV (Vinf~1km/s) Low-energy ~20m/s capture For low-energy transfer/capture, the solar tidal force can effectively reduce Vinf without ΔV.

Figure 1. An illustration of GRAIL’s mission design, including a 21-day launch pe- riod and two deterministic maneuvers for both GRAIL-A and GRAIL-B, designed to separateGRAIL(NASA) their lunar orbit insertion trajectory times by 25 hours. (J.S. Graphic Parker, courtesy of Roncoli and Fujii.R.L.5 Anderson, AAS, 2011) EQUULEUS trajectory (example) 25 transfer time a mission may be designed that places multiple spacecraft into different orbits at the Moon using a single without requiring a large amount of fuel. The GRAIL mission is one practical example: GRAIL’s two spacecraft are launched aboard the same Delta II rocket and each performs two deterministic maneuvers during their trans-lunar cruises to separate their lunar arrival times by 25 hours.5, 6 This separation is a benefit to the spacecraft operations team. Finally, another compelling reason to use a low-energy transfer for a robotic mission to the Moon is that one can construct a realistic 21-day launch period using minimal fuel, as will be demonstrated in this paper. GRAIL’s mission involves at least 21 launch opportunities, such that any launch date sends the two spacecraft on two transfers that arrive at the Moon on the same two dates, again separated by 25 hours. Conventional lunar transfers can achieve the same result, though they typically require numerous Earth phasing orbits and/or lunar flybys that add complexity and radiation exposure to the mission. Numerous researchers have explored the trade space of low-energy lunar transfers since the 1960s using a variety of different techniques. In 1968, Charles Conley was among the first people to demonstrate that a trajectory may be designed to place a spacecraft in an orbit temporarily captured by the Moon without requiring an orbit insertion maneuver.13 His technique takes advantage of dynamical systems tools found in the planar circular restricted three-body problem. Later, in 1990, Belbruno and Miller developed a targeting technique to design a low-energy trajectory for the mission.1 Ivashkin is among many other people to employ similarly useful targeting techniques to generate low-energy lunar transfers.14, 15 Since 2000, several authors have continued to explore the dynamical systems methodology that Conley explored to generate low-energy transfers.16 No practical methods have ever been found to analytically generate a low-energy transfer; hence, all progress to date has involved some sort of numerical or iterative technique to build the transfers. Recent work has begun to systematically survey the trade space of low-energy lunar transfers by building entirely ballistic transfers between the Earth and (1) lunar libration orbits,9–11 (2) low

2 Mission Utilizing Gateway Construction Opportunity Launch & Separation Lunar transfer Lunar orbits, surface, trajectory Earth-Moon(EM) L1/L2, Sun-Earth(SE) L1/L2, ~1 week Interplanetary (Mars, Lunar swing-by asteroids,…) Trajectory (Vinf~1km/s) correction ΔV ~20m/s

1) To Sun-Earth L1/L2: Ballistic transfer without ΔV

2) To Lunar Surface: Landing with 2.5km/s ΔV

3) To Lunar orbits, Earth-Moon L1/L2: Direct insertion or Low-energy transfer/capture by reducing Vinf.

4) To Interplanetary (Mars, asteroids): Escaping by leveraging Vinf.

26 4) To Interplanetary Space (Mars, etc) Launch & Separation Lunar transfer Vinf leveraging by Escape from trajectory solar tidal force Earth

~1 week 5 to 8 months Trajectory Lunar swing-by Lunar swing-by correction ΔV (Vinf~1km/s) (Vinf~2.5km/s) ~20m/s Earth escape Vinf Exploiting solar gravity (tidal ~1.5km/s forces) can increase the Vinf with respect to the moon!! Sun-Earth L1 Lunar swingby #2 (Vinf=2.41km/s)

In this case, the maximum Vinf Sun w.r.t. the Earth is about 1.5km/s. Lunar swingby #1 The interplanetary trajectory can (Vinf=0.899km/s) be design under two-body problem (patched conics) with Vinf<1.5km/s.

DESTINY+ Trajectory (Ozaki, et al., 2019) 27 4) To Interplanetary Space (Mars, etc) Launch & Separation Lunar transfer Vinf leveraging by Escape trajectory solar tidal force from Earth

~1 week 5 to 8 months Trajectory Lunar swing-by Lunar swing-by correction ΔV (Vinf~1km/s) (Vinf~2.5km/s) Reachable Planet for each Earth departure Vinf (under Hohmann transfer assumption) Planet Earth departure Vinf, km/s Mercury 7.53 2.50 Asteroid Depending on the body Mars 2.94 Jupiter 8.49 10.29 We cannot reach of them with Vinf=1.5km/s 4) To Interplanetary Space (Mars, etc) Launch & Separation ΔVEGA/ Other EDVEGA planet

Lunar Lunar swing-by and Earth swing-by swing-by Earth escape m0=14kg V-infinity leveraging transfer (or Isp=2100s ΔVEGA/EDVEGA) can effectively T=1.1mN increase the Earth departure Vinf. Duty=70%

If we want to increase Vinf from 1.5km/s to 3.0km/s (reachable to Mars), we need 0.83km/s ΔV (about half of Vinf increment) ToF=745日

V-Infinity Leveraging Transfer to Mars 29 Gateway Metro Map

Earth Sun-Earth ΔV~10m/s L1/L2 ~0.5yr

Lunar Ballistic,0.5~1yr Transfer Orbit Interplanetary ΔV~500m/s <1week Earth Departure ΔV~10m/s Vinf<1.5km/s ΔV~10m/s ΔV~10m/s 0.5~1yr 0.5~1yr 0.5~1yr To Mars※: ΔV~10m/s ΔV~800m/s 0.5~1yr ~2yr Gateway ※Note: This number does ΔV~10m/s (L2 NRHO) not include Mars orbit L1 Several months insertion ΔV ΔV~500m/s Low- Lunar Orbit ΔV~1500m/s

Lunar Surface

30 Possible Small Sat Mission Utilizing LOP-G When we assume that we can deliver 6U CubeSat to LOP-G (or SLS/Artemis-1 like trajectory), the following missions are possible.

Possible Using Challenging, Target Current Technology Could Be Possible in the Future Moon Moon orbiter (Soft?) Landing Flyby to NEAs Rendezvous to NEAs, Asteroid Exploration to main belt asteroids Lagrange Points Earth-Moon halo, ー Sun-Earth halo Mars, Venus Flyby exploration, ー Orbiter? Lander? Outer Planet Dependent exploration ー (Stand alone mission could be possible if innovative technologies are developed) Bold: Possible missions by SLS, Artemis-2

31 Summary

What is the Lunar Gateway?

Which orbit can the spacecraft transfer from lunar transfer orbit (Gateway construction opportunities)

Which type of mission can the small spacecraft do by utilizing the Gateway opportunity.

32