Skip to main content

ORIGINAL RESEARCH article

Front. Space Technol., 08 December 2022
Sec. Space Robotics
This article is part of the Research Topic Advancing robotic exploration of Asteroids and Comets: a threat, an opportunity, or both? View all 6 articles

Redefining low Earth orbit as a parking orbit for flexible and economical Earth departure in deep space missions

  • Institute of Space and Astronautical Science (ISAS), Japan Aerospace Exploration Agency (JAXA), Sagamihara, Japan

There are a number of deep space probes that are currently in operation with diverse destinations and objectives. For example, the Japanese Hayabusa2 and the U.S. OSIRIS-REx missions are both sample returns, targeting different near-Earth asteroids. Europe’s ExoMars and the U.S. Perseverance are orbiting and roving Mars as precursors of future manned explorations. Conventionally, deep space missions require dedicated launch vehicles for each mission. The interplanetary Earth departure trajectory from the low Earth orbit (LEO) usually lacks flexibility and efficiency. Furthermore, innovative and reusable launch systems have been researched and developed by multiple organizations including private sector organizations such as SpaceX and Blue Origin. It is expected that the unit cost per launch weight to LEO be significantly reduced by rideshare mass transportation executed by using reusable mega launchers in the near future. This study aims to fill the transportation gap between LEO and deep space by realizing a flexible and economical interplanetary Earth departure without sacrificing the arbitrariness of LEO, target V-infinity vector, and target Earth departure epoch. Thus, the one-revolution Earth free-return orbit (1rEFRO) and the consequent Earth gravity assist (EGA) are introduced to separate the velocity increment and direction adjustment. The planetary free-return and EGAs are common in interplanetary missions; however, a comprehensive study on the flexibility, economic efficiency, and arbitrariness of the sequence (1rEFRO + EGA) originating from LEO was not explicitly found. After describing the necessary coordinate frames, LEO’s orbital elements, 1rEFRO, and the terms ‘flexibility’ and ‘economic efficiency’ are defined in Section 2. Then in Section 3, the two-body-based preliminary orbit design method is proposed and formulated. Section 4 aims to reveal LEO’s comprehensiveness as efficient parking orbits when adopting the 1rEFRO + EGA sequence, using the newly proposed “ΣVEt LEO i-Ω diagram”. Section 5 describes a detailed orbit design constructed based on multi-body propagation and optimization to confirm the feasibility, flexibility, and economics of the solution and the usefulness of the initial solution given by the preliminary design method formulated in Section 3.

1 Introduction

Currently, several spacecraft with deep space destinations are in operation or planned, with diverse objectives. For example, the asteroid sample-return missions, such as the Japanese Hayabusa2 and the U.S. OSIRIS-REx, aim to help elucidate the origin of the solar system, and Europe’s ExoMars Orbiter and the U.S. Perseverance Rover are significant as precursors for the future manned exploration of Mars. Their destinations are common in terms of deep space; however, the interplanetary departure conditions are likewise diverse. Thus, deep space missions conventionally require dedicated launch vehicles for each mission, and the interplanetary Earth departure from the low Earth orbit (LEO) usually lacks flexibility and efficiency.

Furthermore, innovative and reusable launch systems have been researched and developed by various organizations including private sector organizations such as SpaceX and Blue Origin. The unit cost per launch weight to LEO is expected to be significantly reduced by rideshare mass transportation by reusable mega launchers in the near future. This might accelerate deep space missions by multiple spacecrafts (i.e., asteroid mining and sample return, Lunar and Mars cargo transportation, reusable deep space transportation, deep space positioning system, etc.).

This study aims to fill the transportation gap between LEO and deep space by realizing flexible and economical interplanetary Earth departure without sacrificing the arbitrariness of LEO, target V-infinity vector, and target Earth departure epoch. Thus, we introduce the one-revolution Earth free-return orbit (1rEFRO) and the consequent Earth gravity-assist (EGA) to separate the velocity increment and direction adjustment. Planetary free-return and EGAs are common in interplanetary missions; however, a comprehensive study on the flexibility, economic efficiency, and arbitrariness of the sequence (1rEFRO + EGA) originated from LEO was not explicitly found.

A planetary free-return orbit (PFRO) is an interplanetary orbit in which departure from the vicinity of the subject planet is guaranteed to be followed by an engagement with the same planet in a certain period. Russel and Ocampo, (2005) formulated all possible PFROs that can be realized on any planet. A 1rEFRO stands for a PFRO having a one-revolution flight period from and back to the subject planet Earth. The 1rEFRO is also called a 1-year Earth synchronous orbit (ESO). The visualization of Earth departure conditions by the V-infinity globe (Strange et al., 2007) and V-infinity diagram (Kawakatsu, 2009) was proposed, and various related investigations have been conducted to date. For example, the flexible interplanetary parking method in combination with electrical propulsion (Ikenaga et al., 2015), the possibility of Earth departures from the geostationary transfer orbit (GTO) (Ikenaga et al., 2015; Ito et al., 2022), the utilization of free-return orbit to gain robustness in Mars orbit insertion (Takahashi et al., 2019), and the asteroid flyby cycler trajectory design (Ozaki et al., 2022) have been investigated. Furthermore, in actual deep space missions, Japanese Hayabusa (Kawaguchi et al., 2006), Hayabusa2 (Tsuda et al., 2016), and many other spacecrafts have incorporated 1rEFRO and the subsequent EGA into their trajectory design as the DVEGA or EDVEGA scheme amplifying the V-infinity norm in combination with the deep space maneuver following the fuel-free orbital maneuvering by EGA.

In this study, the possibilities of LEO, as a parking orbit, toward deep space are investigated by adopting the 1rEFRO + EGA sequence departing from LEO and assuming that all delta-V are impulsive. First, after describing the necessary coordinate frames, LEO’s orbital elements, 1rEFRO, and the terms ‘flexibility’ and ‘economic efficiency’ are defined in Section 2; their parameters are identified and the objectives are clarified. Then, in Section 3, the two-body-based preliminary orbit design method is proposed and formulated. Section 4 reveals LEO’s comprehensiveness as efficient parking orbits when adopting the 1rEFRO + EGA sequence using the newly proposed “ΣVEt LEO i-Ω Diagram”. The orbital elements i (inclination) and Ω (longitude of the ascending node) represent the vertical and horizontal axes on the proposed diagram, respectively. Contours, such as the ratio of achievable V-infinity directions among 4π steradians or the number of preliminary solutions that can achieve the specific target V-infinity vector after the 1rEFRO + EGA sequence, are expected to be drawn on the diagram. The proposed diagram is useful because the feasibility of the sequence departing from the specific LEO can be simply verified. Then, in Section 5, a detailed orbit design based on multi-body propagation and optimization is performed through the LEO-rideshare asteroid exploration mission tentatively set up for this study. Through the detailed design, the feasibility, flexibility, and economics of the solution, as well as the usefulness of the initial guess of the preliminary design method formulated in Section 3 are validated. Finally, Section 6 concludes this study.

2 Background

First, we define the coordinate system required for the discussion of the planetary free-return orbit and provide the definition of 1rPFRO in Section 2.2. Subsequently, in Section 2.3, flexibility and economics are defined as indicators of Earth departure to the interplanetary orbit, and the parameters and evaluation criteria for each indicator are summarized.

2.1 Coordinate system

An inertial coordinate system ΣVBt is introduced with respect to the orbital position and velocity vector of the planet at a certain time t. If the planet is the Earth, the suffix is changed to ΣVEt. If it is based on the Earth’s orbital velocity at time t=t1, it is denoted as ΣVEt1. Using the orbital position RBt and orbital velocity VBt of the planet at time t expressed in an arbitrary inertial coordinate system, each axis’ unit vector (XVBt, YVBt, ZVBt) in ΣVBt is defined as follows.

YVBt=VBtVBtZVBt=RBt×VBtRBt×VBtXVBt=YVBt×ZVBt(1)

The direction of the Earth departure V-infinity vector V is represented by azimuth α with respect to the X-axis direction of ΣVBt and elevation δ with respect to the XY plane. When the norm is denoted by V, the vector V on ΣVBt is expressed together with the domain of definition as follows.

VVBt=Vcosδcosαcosδsinαsinδ180°α180°90°δ90°(2)

Considering the LEO’s orbital element with respect to ΣVEt at time t, let iLEOVEt be the inclination, ΩLEOVEt be the longitude of ascending node, and aLEOVEt be the semi-major axis. When the orbital radius is rLEO, then aLEOVEt=rLEO. The relation of the elements to ΣVEt is described in the right side. iLEOVEt is the angle between LEO’s orbital plane and the XY plane of ΣVEt, ΩLEOVEt is the angle in the longitude of the ascending node with respect to the X-axis in ΣVEt, and the definition range of (iLEOVEt; ΩLEOVEt) is given as follows:

0degiLEOVEt180deg0degΩLEOVEt<360deg(3)

In addition, the J2000EQ coordinate frame is used in this study as a representative inertial coordinate system and is referred to as “J2000EQ”. The frame is defined based on the Earth’s Mean Equator and Mean Equinox (MEME) at the epoch J2000. The x-axis is aligned with the mean equinox, the z-axis is aligned with the Earth’s rotational axis at the epoch, and the y-axis is obtained by the outer product of the z- and x-axes.

2.2 Free-return trajectories

A planetary free-return orbit departs from the vicinity of a subject planet with the promise of meeting the same planet again at a certain period afterward. Russel and Ocampo, (2005) proposed a formulation of all possible free-return orbits for a given subject planet under the assumption that the planet’s own gravity is neglected. They proposed a method to systematically identify all possible free-return orbits, especially when the orbit of the subject planet is circular. In their method, free-return orbits are classified into three types: 1) full-revolution, 2) half-revolution, and 3) generic. In this study, full-revolution orbits, especially one-revolution planetary free-return orbits, are introduced and denoted as 1rPFRO. In particular, when the subject planet is the Earth, this study refers to 1rEFRO.

To derive conditions of 1rPFRO, a sphere whose radius is the desired V is placed at the center of the subject planet. Then, the inverse vector of the planet’s orbital velocity vector VB is placed from the planet center and a sphere of radius VB is drawn, centered on the obtained vector’s tip position. The circular ring at the intersection of those two spheres yields the directions of V satisfying the insertion to 1rPFRO. The aforementioned description is visualized at the bottom of Figure 1, taking the Earth as the target object as an example.

FIGURE 1
www.frontiersin.org

FIGURE 1. Inertial coordinate frame ΣVBt defined by position RBt and velocity VBt at epoch t of the planetary body and V-infinity vector VVBt expressed by azimuth α and elevation δ (top-left), the definition of LEO’s orbital elements (iLEOVEt, ΩLEOVEt,aLEOVEt) with respect to the inertial frame ΣVEt (top-right), and the example of the visualized circular ring representing the directions of V satisfying the insertion to 1rEFRO (bottom).

2.3 Definition of flexible and economical Earth departure

This study identifies flexibility and economics as desirable Earth departure indicators for deep space missions. The objectives for each indicator are defined in the subsequent sections.

2.3.1 Parameters for flexibility

The departure epoch determines the orbital phase of the Earth; V at Earth departure varies widely in magnitude and direction depending on the target planetary bodies or the interplanetary transfer trajectories. The parameters are summarized as follows.

• Departure epoch tdep

• Norm of the V-infinity vector V=V

• Direction of the V-infinity vector (α, δ)

All of these parameters take arbitrary values depending on the target object’s orbit. In this study, the term “flexibility” is used when an Earth departure has smaller constraints on the values of each of the aforementioned parameters. Efficient Earth departure conditions for rendezvous or flyby missions to various minor and major bodies in the solar system are widely distributed over the aforementioned four parameters.

2.3.2 Parameters for economics

We deal with orbital maneuvers to obtain the target V-infinity vector V starting from LEO. Therefore, the economic efficiency is evaluated by the addition of the amount of delta-V implemented in each orbital maneuver. Here, the leverage L is introduced, which is defined by the following equation.

L=Vi=1nVix+Viy+Viz2μerLEOvLEO(4)

where Vix Viy Viz is the delta-V for each axis in the spacecraft coordinate frame at the i-th maneuver out of n times, μE=GME is the Earth’s gravitational constant, and (rLEO, vLEO) is the distance and velocity of the spacecraft concerning the Earth’s center of gravity right before departing LEO. When each delta-V is implemented without division into each axis, the first term in the right-hand denominator of the aforementioned equation can be expressed as i=1nVi. The right-hand denominator is the sum of the delta-V to obtain the desired infinite velocity vector starting from LEO, minus the required minimum delta-V to achieve orbital energy C3=0 from LEO. The leverage L is the ratio of the obtained V to the sum of delta-V invested to achieve C3=V2 from the ideally reached C3=0 from LEO. As the value can assume L > 1, it is referred to as “leverage”.

For L, we consider the case where V is ideally obtained from LEO; if we let Lideal be the leverage when the desired Earth departure is established by a single impulsive tangential delta-V of Videal on LEO, the leverage is expressed as Eq. 4. Moreover, let Vesc=2μE/rLEO be the escape velocity on LEO.

Lideal=VVesc2+V2vLEOVescvLEO=VVesc2+V2Vesc=1VVesc2+V2+Vesc=VescV+VescV2+1(5)

As shown in the aforementioned equation, Lideal is expressed as a function of the LEO orbital radius rLEO and the desired V-infinity V. Furthermore, when expressing Vesc as the necessary amount to obtain the escape velocity above LEO and Vdep as the contribution to the net increase in V, each delta-V has the following relationship.

Videal=Vesc+VdepVescVescvLEO=2μErLEOvLEOVdepV2+Vesc2Vesc(6)

Generally, the smaller the V and lower the LEO’s altitude hLEO are, the larger the Lideal will be. However, the required Vdep increases in proportion to the desired V.

3 Designing flexible and economical Earth departure from LEO

Starting from a wide range of LEOs, to achieve a high level of both flexible and economical Earth departure as defined in the previous section, this study adopts an orbital sequence comprised of a transition to 1rEFRO and a subsequent EGA after 1 year and organizes a preliminary trajectory design method in this section. Although the sequence is commonly used in multiple deep space missions, no study has explicitly organized a LEO-based trajectory design method or conducted a comprehensive survey of compatible low Earth orbital elements.

Section 3.1 first describes the major challenges in terms of flexibility and economic efficiency when departing directly from the LEO to the desired interplanetary orbit. Then, the goal and overall picture of the introduction of the aforementioned orbital sequence are presented. Section 3.2 summarizes the preliminary design methodology for the desired Earth departure from the LEO via 1rEFRO and EGA under the assumption of the unperturbed two-body problem and impulsive delta-V.

3.1 Introduction of an orbital sequence departing LEO

This section describes the adopted orbital sequence for the flexible and economical Earth departure from the LEO and its necessity.

3.1.1 General difficulties in departing LEO toward arbitral deep space destinations

There are challenges in terms of the required ΔV for a direct Earth departure from an arbitrary LEO to a transiting orbit to smaller bodies or Mars. As an example, we assume that the LEO (altitude 500 km, orbital inclination i=31°, the longitude of ascending node Ω=0°) was obtained by the launch from the Japanese Tanegashima Space Center (TNSC). The target V-infinity vector Vtar (α=156.95°, δ=45.20°, V=4.70 km/s) is set from Hayabusa2’s first Earth swing-by on 3 December 2015. The transitions from the LEO to Vtar are considered to be directly made from each point on the LEO. Assuming the two-body problem of the Earth and the spacecraft and assuming that each position in the LEO is a perigee, hyperbolic trajectories can be obtained as shown in Figure 2.

FIGURE 2
www.frontiersin.org

FIGURE 2. Examples of hyperbolic departure trajectories that directly transition from a specific LEO to a specific Vtar (top-left, top-center). Corresponding delta-V quantities and leverage L (top-right) as well as the schematic of the adopted orbital sequence (1rEFRO-EGA) (bottom) are shown.

The in-plane and out-of-plane components of the ΔV quantity at each starting point and the leverage L were calculated based on Section 2.3.2 and are shown at the top-right side of Figure 2. The leverage L is significantly smaller than Lideal because the obtained hyperbolic trajectories’ orbital planes differ from the original LEO’s and the departure ΔV requires a large out-of-plane component. For the Vtar direction, it is necessary to be existing in the original LEO’s orbital plane to bring L close to the ideal value shown in the figure. However, this significantly restricts possible low Earth orbital elements and hence limits the flexibility. Furthermore, because the LEO must be designed in accordance with the departure epoch tdep to obtain an economical direct departure toward Vtar, the direct departure from the arbitrary LEOs cannot ensure a both economical and flexible Earth departure to the arbitrary Vtar.

3.1.2 Employed orbital sequence (1rEFRO + EGA)

Therefore, this study adopts the two consecutive orbital operations shown in Figure 2 that start from an arbitrary LEO and allow a transition to an arbitrary target V-infinity vector Vtar at an arbitrary time t=t2. 1) First, the spacecraft implements a tangential impulsive delta-V in LEO at t=t1 for the 1rEFRO insertion. 2) Second, the spacecraft transfers to Vtar through a consequent EGA at t=t2. The key feature of this sequence is the separation of velocity increment and direction adjustment. After obtaining the desired Vtar=Vtar in Eq. 1, the desired direction (αtar, δtar) is obtained in Eq. 2. The time required between 1) and 2), t2t1, is approximately 1 year, which is the Earth’s orbital period TE. The candidates of the departure time t1 in (Eq. 1) are designed referring to the default epoch of t1 by t1def=t2tarTE. The aim of this orbital sequence is to realize an economical Earth departure with high leverage L without compromising the arbitrariness of the orbital elements of the parking orbit, LEO, the target V-infinity vector Vtar, and the target departure time t2tar.

A preliminary design method based on the Sun–spacecraft and Earth–spacecraft two-body problem is proposed in this study for designing an Earth departure trajectory that uses the aforementioned sequence, which is described in detail in this section. In the preliminary design, according to the premise of 1rEFRO, an assumption that the output V-infinity vector at t1, VVEt11o, and the input V-infinity vector at t2, VVEt22i, coincide as follows.

VEt2V2iVEt1V1o(7)

3.2 Preliminary trajectory design method

This section details the preliminary design method for the abovementioned orbital sequence. The reason for the term “preliminary design” is that this method treats the interplanetary orbit as a two-body problem of the center body and the spacecraft as an infinitesimal object with no perturbation. The two bodies are distinguished at each design stage of the sequence: the transition from LEO to 1rEFRO is Earth + spacecraft, 1rEFRO itself is Sun + spacecraft, and the transition from 1rEFRO to EGA is Earth + spacecraft. The preliminary design results in this section are adopted as initial guesses for the orbit design using nonlinear numerical optimization accounting for N-body dynamics in Chapter 0, and the usefulness is evaluated. In the following part of this section, necessary conditions for the preliminary design are first formulated: 1) the insertion condition from LEO to 1rEFRO by a tangential impulsive ΔV on LEO, and 2) the transition condition from 1rEFRO to Vtar by EGA after 1rEFRO. The procedure of the specific preliminary design method is then explained by organizing the assumptions and input–output relations.

3.2.1 LEO to 1rEFRO transition (tangential impulsive ΔV)

Let V1o1rEFRO denote the V-infinity vector satisfying the transition condition to 1rEFRO with the desired norm of Vtar, and let V1o1rEFRO denote its set. As explained in Section 0, the endpoints of the direction vectors in V1o1rEFRO form a circular ring that coincides with the intersection of the two spheres. When expressed in terms of the azimuth and elevation (α, δ) on ΣVEt1, the vector V1o1rEFRO has a backward component (α < 0) relative to the Earth’s orbital velocity vector. Furthermore, the larger V1o1rEFROVEt1Vtar is the smaller the area inside the circular ring will be. These characteristics can be read from the equations, and are presented in Section 3.2.1.1, including the derivation process of V1o1rEFRO. Then in Section 3.2.1.2, the set V1o1rEFRO is further narrowed down to those capable of departing from a particular LEO to 1rEFRO.

3.2.1.1 Set of V to achieve 1rEFRO insertion

First, by V expressed using the azimuth and elevation angles (α, δ) on ΣVEt1 in Eq. 1, the orbital velocity vector of the spacecraft at time t1, VPt1, on ΣVEt1 , can be expressed as follows:

VPt1VEt1=VEt1VEt1+VVEt1=VEt1010VEt1+VcosδcosαcosδsinαsinδVEt1=VcosδcosαR˙e+VcosδsinαVsinδVEt1(8)

To establish a planetary return after one orbit as 1rEFRO, the orbital period of the spacecraft must be the same as that of the subject planet. Inputting (α1o, δ1o) as the azimuth and elevation angle that satisfy the condition of transition to 1rEFRO, the following equation holds.

VPt1=VEt1VPt12=VEt12V2cos2δ1ocos2α1o+VEt1+Vcosδ1osinα1o2+V2sin2δ1o=VEt12V2+2VEt1Vcosδ1osinα1o=0cosδ1osinα1o=V2VEt1(9)

Because cosδ0 always holds from the domain of definition and the value range y of cos1 is 0yπ, solving Eq. 9 for δ1o yields:

cosδ1o=V2VEt1sinα1oδ1o=cos1V2VEt1sinα1o,whenδ1o0δ1o=cos1V2VEt1sinα1o,whenδ1o<0(10)

Next, the possible range of α1o is considered. Here, assuming that V/2VEt11 holds as a prerequisite cosδ1o0, the inverse function of the cosine function is a real number, and the value range y of the sin1 function is π/2yπ/2. The possible range of α1o can be constrained by using Eq. 10 as follows.

0V2VEt1sinα1o12VEt1V1sinα1o0sinα1oV2VEt1πsin1V2VEt1α1sin1V2VEt1(11)

From Eq. 10 and Eq. 11, the possible range of δ1o can also be constrained as follows.

δ1ocos1V2VEt1,whenδ1o0δ1ocos1V2VEt1,whenδ1o<0(12)

Eq. 11 and Eq. 12 show that α1o<0, i.e., V1o1rEFROVEt1 is pointing backward with respect to the Earth’s orbital velocity vector, and the larger the V1o1rEFROVEt1 the smaller the range of possible values for (α1o, δ1o).

3.2.1.2 Set of V to achieve 1rEFRO from the specific LEO tangentially

A subset of V1o1rEFRO that is capable of being inserted into 1rEFRO from a particular LEO is denoted as V1o1rEFROLEO and is derived in the following part of this section. First, from this definition, the relation V1o1rEFROLEO=VtarVtar holds. Since, the end points of V1o1rEFRO form a circular ring, the problem of finding the set V1o1rEFROLEO is reduced to a geometrical analysis to find the intersections between the ring and the orbital plane of a given LEO. Hereafter, the equations will be organized in an orderly manner.

First, from this definition, the ring exists in the Y=Ly plane, and its center is on the Y-axis of ΣVEt and let the ring’s radius be Lr. Together with Ly, they can be expressed using the V-infinity vector norm Vtar and the Earth’s orbital velocity norm VEt at epoch t as follows.

Ly=Vtar22VEtLr=Vtar1Vtar2VEt2(13)

The existence of the intersections at time t1 is investigated for LEO’s orbital elements (iLEOVEt1; ΩLEOVEt1) and the number of intersections is denoted by n. Here, we focus on the intersection line between the Y=Ly plane where the ring exists, as well as LEO’s orbital plane. Based on the geometric relationship, a discriminant equation KiLEOVEt1,ΩLEOVEt1,Ly,Lr and the relationship between K and n were investigated; the latter is determined as Eq. 14. Because the problem is to find the intersection of a ring and a plane, the maximum n-value is two, and the minimum is zero. For simplicity, (iLEOVEt1,ΩLEOVEt1) are expressed as (i, Ω) in the equations in the rest of this section. We divided the problem into the following five cases.

n=0whenK>0n=1whenK=0n=2whenK<0(14)

• Case 1: (i=90°)∩{(Ω=0°)∪(Ω=180°)}

As the intersection line between Y=Ly and the LEO does not exist in this case, n=0 as well as V1o1rEFROLEOVEt1.

• Case 2: (i=90°)∩(Ω0°)∩(Ω180°)

The intersection line between Y=Ly and the LEO exists and is parallel to Z axis in this case. K and V1o1rEFROLEOVEt1 are derived as follows.

K=LytanΩ2Lr2V1o1rEFROLEOVEt1=LytanΩLy±K(15)

• Case 3: (i=0°)∪(i=180°)

The intersection line between Y=Ly and the LEO exists in this case. The line is parallel to the X-axis and the line’s Z component is zero. Therefore, it is always n=2, and V1o1rEFROLEOVEt1 is calculated as follows.

V1o1rEFROLEOVEt1=±LrLy0(16)

•Case 4: (i0°)∩(i90°)∩(i180°)∩{(Ω=0°)∪(Ω=180°)}

The intersection line between Y=Ly and the LEO exists in this case. The line is parallel to the X-axis, but the line’s Z component is NOT zero. K and V1oLEOVEt1 are expressed as:

K=Ly2tan2iLr2V1o1rEFROLEOVEt1=±KLyLycosΩtani(17)

•Case 5: (i0°)∩(i90°)∩(i180°)∩(Ω0°)∩(Ω180°)

The intersection between Y=Ly and the LEO exists in this case, but the line is not parallel to both the X and Z axes. K and V1o1rEFROLEOVEt1 are derived as in the following equation.

K=cos2Ωtan2i1+sin2Ωtan2iLy2Lr2V1o1rEFROLEOVEt1=sinΩcosΩtan2i1+sin2Ωtan2iLy±K1+sin2Ωtan2iLysinΩtanicosΩsinΩ1+sin2Ωtan2iLy±K1+sin2Ωtan2i(18)

Here, the set V1o1rEFROLEOVEt1 and the number n were derived. The V-infinity vectors in the set are classified in terms of their directions. There are a total of eight types of directions that satisfy 1rEFRO. For azimuth α, those satisfying 180°α<90° are referred to as “inbound”, and those satisfying 90°<α0° are referred to as “outbound”; the former flies inside the Earth’s orbit and then outside to make a return to Earth and the latter flies in the opposite order. In the case of α=90°, as the 1rEFRO orbit differs only in inclination from the Earth, the spacecraft returns after half a period. This is therefore technically identified as half-revolution EFRO (hrEFRO) rather than 1rEFRO.

For elevation δ, those with 90°δ<0° depart to the Southern Hemisphere against the ecliptic plane and return to Earth via the Northern Hemisphere, while those with 0°<δ90° fly in the inverse order. The following bullets show the total eight combinations of possible V1o1rEFROLEO directions by classifying α and δ. Notably, that there is no V-infinity vector that satisfies α=90°δ=0°.

• Type 1: Inbound (−180° < α < −90°), South first (−90°≦δ < 0°)

• Type 2: Inbound (−180° < α < −90°), in the ecliptic plane (δ = 0°)

• Type 3: Inbound (−180° < α < −90°), North first (0° < δ≦90°)

• Type 4: hrEFRO (α = −90°), South first (−90°≦δ < 0°)

• Type 5: hrEFRO (α = −90°), North first (0° < δ≦90°)

• Type 6: Outbound (−90 < α < 0), South first (−90°≦δ < 0°)

• Type 7: Outbound (−90 < α < 0), in the ecliptic plane (δ = 0°)

• Type 8: Outbound (−90 < α < 0), North first (0° < δ≦90°)

If two V1o1rEFROLEO are obtained, then, either or both signs of (α, δ) differ in between those two V-infinity vectors considering the symmetry of the intersection by the ring and the plane. In summary, the combination can be seen from one of the following seven types of pairs: (1, 3) (1, 6) (1, 8) (2, 7) (3, 6) (3, 8) (4, 5).

3.2.2 1rEFRO to Earth departure (unpowered EGA)

After departing the LEO to 1rEFRO, the spacecraft will meet Earth again at time t2, one orbital period later (approximately 1 year after t1), for an EGA with no ΔV. Here, the final Earth departure target condition is Vtar. The input and output V-infinity vectors of the gravity assist are denoted as V2i and V2o, respectively. Considering the important assumption of Eq. 7 according to the premise of 1rEFRO, their direction and magnitude are given by the following equations. Note that the display of the coordinate frames on the left shoulder of vectors is omitted in this section since the frame can be as long as they are an inertial coordinate frame centered on the Earth.

V2i=V1o1rEFROLEO,V2i=VtarV2o=Vtar,V2o=VtarV2o=V2i=Vtar(19)

As the input and output V-infinity vectors in the EGA are determined, the swing-by trajectory is uniquely assessed as well. In this section, the equations to determine the validity of the perigee radius and the position vector r2π when passing perigee at the epoch t2 will be developed.

The swing-by trajectory derived from V2i and V2o is symmetrical about r2π. When φ2 is the rotation angle between the input and output V-infinity vectors, S^2I and S^2O are the unit vectors of each V-infinity vector, r^2π is the unit vector of r2π, and they are derived together with r2π as follows.

S^2iV2iV2i=V2iVtarS^2oV2oV2o=V2oVtarφ2cos1V2iV2oV2iV2o=cos1V2iV2oVtar2=cos1S^2iS^2or2πμeV21sinφ221r^2πS^2iS^2o2sinφ22=V2iV2o2V2sinφ22r2πr2πr^2π=μe2V21sinφ221S^2iS^2osinφ22=μe2V41sinφ221V2iV2osinφ22(20)

The EGA at t2 is valid when the perigee radius r2π derived in Eq. 20 falls in between the thresholds (rπminr2πrπmax). As V2i=V1oV1o1rEFROLEO holds as explained previously, two vectors exist as V2i at maximum. The swing-by feasibility at time t2 is checked for each V2i to determine whether the corresponding V1o1rEFROLEO can be selected or not.

3.2.3 Proposal of preliminary trajectory design procedure

Using the relational equations derived up to the previous section based on the unperturbed two-body problem, this section proposes a preliminary orbit design procedure for an Earth departure starting from a specific LEO and going through 1rEFRO + EGA. First, if this procedure is assumed to be a function gpre, the input–output relationship is given as follows.

cpre,{t1},ΔV1},{t2=gpret0,rEP0,vEP0,t1mid,t2tar,V2tar(21)

As shown in the aforementioned equation, the inputs are the relative position velocity to Earth (rEP0,vEP0) in the LEO at the epoch t0, the center value t1mid of the search range to obtain the optimum LEO-departure epochs, the target Earth departure epoch t2tar, and the V-infinity vector V2tar. The output of the process is the epoch t1, the required delta-V amount of ΔV1, and the final Earth departure epoch t2. Assuming to obtain multiple sets of (t1,ΔV1,t2) as the result, their outputs of gpre are treated as the sets t1},ΔV1},{t2 of each variable and a number cpre of values exist in a set. This preliminary design procedure consists of a seven-block flow as shown in Figure 3, from (A) to (G). In block (E), the optimization routine calculates the optimum time t1opt at which the transition to 1rEFRO from a specific LEO is possible, and in block (F), the EGA feasibility to V2tar at t2 is evaluated. The branching in block (E) is intended to obtain t1opt corresponding to VV1o1rEFROLEO in each branch, where, at most, two V exist. For this reason, one of the two ends of the interval is given to each branch as the initial guess of the instance for their optimization routines. The following is a description of the process and the necessary equations for each block.

FIGURE 3
www.frontiersin.org

FIGURE 3. Flowchart of the preliminary design procedure of LEO-1rEFRO-EGA trajectories.

This method is referred to as “preliminary” because, as mentioned previously, it is mainly based on an unperturbed two-body problem. The specific assumptions for each phase are shown in the bullets below.

t0 to t1 propagation in LEO: Earth-centered, perturbation considered

• LEO to 1rEFRO transition at t1: Earth-centered, no perturbation

• 1rEFRO cruising: Sun-centered, no perturbation

• 1rEFRO to EGA at t2: Earth-centered, no perturbation

By simplifying the dynamics after LEO departure to an unperturbed two-body problem, this method can simply obtain the initial design of the Earth departure-orbit starting from LEO. In contrast, from t0 to t1, when the spacecraft remains in the LEO, there are non-negligible variations to the orbital elements due to Earth’s higher-order gravitational potential, with third bodies such as the Moon and Sun, atmospheric drag, etc. Because the variation of LEO significantly affects the entire trajectory after LEO departure, the spacecraft state at t1 is propagated with necessary perturbations.

3.2.3.1 (A) Read inputs

In this block, the conditions of the LEO (t0,rEP0,vEP0), the center value of search range t1mid, and the target Earth departure conditions (t2tar,V2tar) are read and can be referenced in the later flows.

3.2.3.2 (B) Load propagation, optimization, and swing-by conditions

In this block, the various conditions (LEO propagation, optimization, and swing-by) are read prior to specific calculations. For the swing-by condition, the minimum and maximum perigee distances rπmin and rπmax are read, which will be referred to in the swing-by validation in Block (F).

3.2.3.3 (C) Pre-processing

In this block, prior to the optimization calculation in Block (E), the necessary constants are prepared through the processes. The required V-infinity vector norm for the LEO departure to 1rEFRO is the same as V2tar. The main process of this block is to prepare the interval (search range) for the instance t1loop of the optimization routine in Block (E). The interval is both bounded and closed expressed as t1a,t1b. The midpoint is set to the input parameter t1mid, and the diameter is set to the specified LEO’s orbital period TLEO. If t1mid is not specified or null, the default value is calculated as t1mid=t1deft2tarTE. The left- and right-bound values are prepared as follows.

t1a=t1midTLEO/2t1b=t1mid+TLEO/2(22)

3.2.3.4 (D) Initialization and start of optimization routines

This block initializes and starts the optimization routines (a) and (b). The difference of the routines is the initial guess of the instance t1loop. As shown in Section 3.2.1.2, there are, at most, two directions of {V1o1rEFRO} for the same LEO orbital element, and they have different departure epochs and orbital phases in the LEO. Therefore, it is intended to obtain each of two different solutions in routines (a) and (b) by starting the optimization calculations from each end of the interval, t1loop=t1a for routine (a), and t1loop=t1b for routine (b), with a range of the LEO’s orbital periods TLEO and t1mid as the center value.

3.2.3.5 (E) Optimization routine

In this block, the processes are performed as an optimization calculation to determine the epoch optimum as t1. The main task is to calculate two independent V-infinity vectors. One is denoted as V1oloop, which is derived from the current position, and velocity vectors of the spacecraft at time t1loop, which is an instance of the routine, considering the desired tangential impulsive delta-V of ΔV1loop. The other is from the set V1o1rEFROLEOloop calculated from the osculating orbital elements (i1loop; Ω1loop) in LEO at t1loop (as described in Section 3.2.1). Then, the square of the angular difference between those vectors is calculated as the objective function value of fpre. When both V-infinity vectors are almost coincident with each other and satisfying the specified convergence conditions, the calculation is terminated, and the time t1loop of the latest loop is outputted as the optimal value t1opt.

When calculating the desired, required tangential delta-V norm V1oloop on LEO at t1loop, it is derived using the known V-infinity norm V1oloop=V2tar. Together with the resulting post-delta-V velocity vector VEP1loop, V1oloop is calculated as follows.

V1loop=V1oloop2+2μerLEOVEP1o=V2tar2+2μerLEOμerLEOVEP1loop=vEP1loop+V1loopvEP1loopvEP1loop(23)

The resulting V-infinity vector V1oloop can be obtained from the post-delta-V state of the spacecraft (rEP1loop,VEP1loop) using the angular momentum vector h and Laplace vector P along the definition of an unperturbed hyperbolic trajectory as the following equations.

h=rEP1loop×VEP1loopP=μerEP1looprEP1looph×VEP1loopQ=h×PP^=PP,Q^=QQe=PμEV1oloop=V1otar1eP^+11e2Q^(24)

3.2.3.6 (F) Post-processing

In this block, processes are performed to determine whether the target V-infinity vector V2tar can be obtained through the Earth swing-by 1 year after t1opt is obtained in the previous block. The major tasks are to calculate the EGA epoch t2 from t1opt by adding TE, preparing the input and output V-infinity vectors of the swing-by considering the assumption in Eq. 7, and calculating the perigee radius r2π as described in Section 3.2.2. Then, we determine whether rπminr2πrπmax is satisfied. If so, we assign one, and if not, we assign zero to the flag variable flagisEGAvalid and finish the processes in this block.

3.2.3.7 (G) Output variables

In this block, the results of blocks (E) and (F) in the two branches (a) and (b) are checked, respectively, and the output variables for the entire procedure are prepared. First, we initialize the integer variable that stores the number of preliminary solutions by cpre=0, then for each branch, if flagisEGAvalid=1 then increment cpre and store the variables to the set t1,V1,t2 (note that if the results are common among both branches, count as 1), and finally, the output variables are returned.

4 Comprehensiveness of LEO as a parking orbit toward various deep space missions

Based on the methods outlined previously, this section illustrates that the 1rEFRO + EGA orbital sequence can realize a flexible Earth departure for deep space economically from an extremely wide range of LEOs. First,” ΣVEt LEO i-Ω Diagram” is newly proposed as an indicator to explain LEO’s comprehensiveness toward deep space. The explanation of the diagram together with its usage is given in the following part of this chapter. Then, LEO is redefined as a parking orbit for various deep space missions, and its utilization is described.

4.1 Coverage of V2° directions achievable by the LEO-1rEFRO-EGA sequence

For the EGA at t2, using the norm of a V-infinity vector V2i=Vtar, the minimum perigee distance rπmin, the maximum perigee distance rπmax, and the Earth gravity constant μE, the maximum value φ2max, and minimum value φ2min of the rotating angle φ2 of the swing-by can be calculated as follows.

φ2max=2sin111+rπminVtar2μEφ2min=2sin111+rπmaxVtar2μEφ2minφ2φ2max(25)

This means that the set V2o of V that can be assumed after the EGA at t2 is the set of vectors obtained by rotating V2i around an arbitrary axis by minimum φ2min and maximum φ2max degrees, respectively. Let the fraction of directions that can be taken as {V2o occupying the unit sphere be defined as V2o coverage and denoted by ccov0,1. The rest of this section givesthe description to calculate the fraction from the given V2i.

First, for a given V2i, the directions that cannot be transitioned due to the perigee restriction are defined as “Steradians of V2o restriction” and are illustrated in Figure 4.

FIGURE 4
www.frontiersin.org

FIGURE 4. Schematic image of steradians of V2o restriction for given V2i (left) and surface areas Aβn,Aβf,Aβnf of the overlapping steradians of V2o restriction for given set of V2i (right).

For a given V2i, it is not possible to transition in the direction of half-peak angle θn (solid angle Θf) in the reverse direction and half-peak angle θf (solid angle Θf) in the forward direction. As there are, at most, two possible V from a particular LEO to 1rEFRO, the number of vectors cpre in V2i that can be taken for a given LEO is (0, 1, 2). When cpre=0, the EGA at t2 is likewise not satisfied, and ccov=0. When cpre=1, using the surface areas An and Af of the resulting steradians of V2o restriction by given V2i, ccov is calculated as follows.

ccov=1An+Af4π(26)

Then, when cpre=2, the V2o infeasible direction for one V2i might be the feasible direction of V2o for the other V2i. Therefore, in calculating ccov, the sum of the overlapped surface area of the steradians of V2o restriction from both given V2i is necessary. There are three possible types of overlaps, and the overlapped surface areas are denoted as Aβn,Aβf,Aβnf, as shown in Figure 4. In calculating each surface area, the angles βn,βf,βnf at which the circumferences of the cones overlap inside the cross-section cut by the plane stretched by both V2i are calculated, and the area of the overlapping spherical cap can be calculated through spherical trigonometry. Finally, ccov is given by using the overlapped surface areas Aβn,Aβf,Aβnf of the spherical caps as the following equation.

ccov=1Aβn+Aβf+2Aβnf4π(27)

4.2 Introduction of ΣVEt LEO i-Ω diagram

The variable ccov is an indicator that simultaneously expresses: 1) whether the transfer from LEO to 1rEFRO at t1 is capable or not, and 2) the percentage of directions achievable by EGA at t2. The variable takes the LEO’s orbital elements (iLEOVEt1,ΩLEOVEt1), the departure epoch t1 from LEO, and Vtar as its arguments. Since LEO is a circular orbit, the elements that determine the orientation of the orbit in inertial space can be limited to two: the inclination angle i and the longitude of the ascending node Ω. Therefore, in comprehensively understanding ccov for the possible LEOs, this section considers calculating the ccoviLEOVEt1,ΩLEOVEt1,Vtar,t1 against a given t1 and desired Vtar across the defined ranges of iLEOVEt1 and ΩLEOVEt1 (0°iLEOVEt1180°; 0°ΩLEOVEt1<360°) to investigate the sensitivity of ccov across the possible i and Ω. The obtained ccov array, plotting the horizontal axis ΩLEOVEt1 and vertical axis iLEOVEt1, is the newly proposed diagram in this study, referred to as” ΣVEt LEO i-Ω Diagram” and shown in the upper part of Figure 5. A ccov contour centered at (iLEO,ΩLEO) = (90°, 0°), (180°, 90°) is observed in the plane stretched by the LEO’s orbital elements (iLEO,ΩLEO). The perigee altitude constraint adopted in this calculation is 500hπ1e6km. The blue-filled area of ± about 4.5° from each center point implies ccov=0 having no V2o solution due to the infeasibility to enter 1rEFRO inside the LEO’s orbital plane. In the process of creating the diagram, the V diagrams at t1 and t2 for each of the (iLEO,ΩLEO) points can be drawn and are shown in the red boxes in Figure 5. The left diagram in each red box shows {V1o1rEFROLEO} calculated at t1. In the right diagram in each red box, {V2o} is identified among the 4π steradians centered at the Earth that satisfy the perigee altitude constraint at t2, and a contour map of its perigee altitude hπ is drawn. In contrast to the right diagram in the left red box, indicating ccov0.99, the right diagram in the right red box has a lower value of ccov0.79. This is because the gray-colored V2o infeasible region largely exists due to the minimum hπ constraint.

FIGURE 5
www.frontiersin.org

FIGURE 5. ΣVEt LEO i-Ω diagram (top left) and the diagram for specific V2tar (top-right) at t=t1 both overlayed with LEO propagated for ± 20 days around t1, and two examples of associated V diagrams at t1 and t2 for specific LEO (i,Ω) of (120 °, 60 °) [deg] (bottom left) and (85 °, 355 °)[deg] (bottom right).

The advantage of the ΣVEt LEO i-Ω diagram is that it can be used to determine whether to depart from the LEO orbit in a simplified manner or not. For the diagram in Figure 5, the time history of LEO’s orbital elements propagated along ±20 days from t1 with consideration of the perturbation is overlayed. The reference value of ccov can be checked when the departure time is shifted from t1.

To evaluate the feasibility of the 1rEFRO-EGA sequence in a more straightforward manner, the top right plot in Figure 5 shows the ΣVEt LEO i-Ω diagram for specific V2tar. In the figure, as an alternative to the aforementioned ccov values, it plots the number cval=0,1,2 of orbital solutions that is feasible for transferring to a desired V2tar and the time history of LEO’s propagated (i, Ω) is overlayed. In Figure 5, cval=1 holds for most of the propagation periods, which means that there is one orbital solution that is feasible up to the EGA even if the LEO departure time t1 is shifted. This is because the time variation of the diagram is small when the eccentricity of the subject planet’s orbit is close to zero.

4.3 Comprehensiveness of LEO when applying the 1rEFRO + EGA sequence

Figure 6 shows the results of the ΣVEt LEO i-Ω diagram for multiple V, while the conditions of t1 and the constraint of hπ are unchanged. As V increases, ccov gradually decreases; however, ccov>0.85 is satisfied in most of the areas even at V=6km/s. Hence, by going through 1rEFRO + EGA from most of the LEOs, the feasible V2o direction on the unit sphere is 100% at V=2.5km/s and >85% at V=6km/s. Although it depends on the hπ constraint, the range set of 500hπ1e6km is realistic in the actual mission by having an opportunity of deep space maneuver during the period of 1rEFRO.

FIGURE 6
www.frontiersin.org

FIGURE 6. Comparison of the ΣVEt LEO i-Ω diagram for various V2tar2.5,3.0,,6.0km/s.

These results demonstrate the comprehensiveness of LEOs by introducing the sequence 1rEFRO + EGA. It is possible to economically depart from the LEO to deep space and achieve an interplanetary orbit by a very wide V2o. In addition to the flexibility of V2o, the possible timing of departures can be as frequent as twice during the LEO’s orbital period, or twice in about 90 min for a 500 km altitude LEO. This provides a high level of time-flexibility compared to the launch window of deep space missions, which generally allows several days for the shift. These results imply that LEOs have great potential as parking orbits capable of hosting a variety of deep space missions economically.

4.4 Re-definition of LEO as a parking orbit for various interplanetary missions

As described in the previous section, the adoption of the 1rEFRO + EGA sequence used in this study has shown that economical and flexible Earth departures from a very wide range of LEOs to any V are possible. This means that LEO, which was thought to be highly constrained for interplanetary Earth departures as described in Section 3.1.1, can be used flexibly and economically as a parking orbit. The following applications are given as examples.

• Multiple deep space probes launched to the LEO by rideshare can independently depart from the LEO to their own interplanetary mission.

• Interplanetary Earth departure by a piggy-back deep space probe on an Earth observation satellite mission.

• Earth departure from a manned space station in the LEO.

• Extremely flexible launch window of the deep space probe to LEO can significantly relax the launch date and time constraints for a launch vehicle.

In utilizing LEO as a parking orbit, consideration must be given to the satellites’ congestion currently in progress in LEO due to the deployment of large-scale constellations. The congestion of LEO affects the probability of colliding with other spacecraft when staying in low earth parking orbit and during the EGA after 1rEFRO. There are several possible measures to this issue. The first is to perform a collision avoidance orbital maneuver for both during the LEO stay and the EGA. When the probability exceeds a certain level, collisions can be avoided by planning and performing a small ΔV in advance. The second is to avoid specific congested LEOs by adjusting the altitude of the parking orbit. A change in the LEO’s altitude affects the ideal leverage value Lideal of the entire sequence. For example, assuming an Earth departure V=4km, Lideal=5.64 at an altitude of 300 km (orbit radius 6,678 km), Lideal=5.56 at a 500 km altitude mainly discussed in this study, and Lideal=5.38 at 1,000 km. Lideal decreases when selecting the higher-altitude LEO but still provides significant leverage. Therefore, it can be said that it is possible to avoid congested LEO by changing the altitude without a significant sacrifice of economic efficiency.

5 Validation of preliminary design method

Based on the initial solution obtained by the preliminary orbit design method proposed in Section 3, we construct a detailed orbit design by numerical optimization using multi-body propagation to validate the usefulness of the initial solution and the method. First, the procedure of the detailed design method is described. Then, the LEO-rideshared asteroid exploration mission is tentatively set up, and the results of the detailed design and the discussions are presented.

5.1 Detailed trajectory design procedure

We describe the detailed orbit design method constructed to verify the initial solution obtained in Section 3. Nonlinear numerical optimization based on orbit propagation in a multi-body dynamical model is used to generate a reference orbit of 1rEFRO + EGA starting from a specific LEO from the departure epoch of t1 to the perigee passage epoch t2 during the EGA. The three phases that comprise this procedure and the instances of optimization and objective function at each phase are shown in Figure 7, together with an illustration.

FIGURE 7
www.frontiersin.org

FIGURE 7. Flowchart of the detailed trajectory design procedure composed of preliminary, multi-body first step, and the multi-body second step design phases (top) and the specifications of the optimization routine in each phase (bottom).

In phase (I), the initial solution of t1 is obtained through the preliminary design method. Then, from the viewpoint of convergence, the multi-body design is achieved in two steps of phases (II) and (III), respectively.

In phase (II) as Step 1, t1 and t2 are taken as instances, and the objective function f1 is defined as the square of the difference between the Earth–spacecraft distance and the target distance rEP2tar at t2. The perigee distance expected in the EGA from phase (I) is substituted in rEP2tar. The amount of tangential ΔV to be done at t1 is the constant necessary to obtain the desired V-infinity vector norm under the assumption of the Earth–spacecraft two-body problem.

In phase (III) as Step 2, the deep space maneuver VDSM, implemented at the intermediate epoch tdsm=t2t1/2, is set as an instance, with t1 and t2 obtained in the previous step and fixed as the orbit propagation period. VDSM has components in three orthogonal directions in inertial space. In this step, the objective function f2 is the sum of the squares of the difference between the infinite velocity vector V2o calculated from the spacecraft state at t2 by using the two-body solution and the target V-infinity vector V2otar. The slight difference (V2oV2otar) at the completion of phase (III) optimization is converted into a delta-V V2π at the perigee passage epoch t2π, also by a two-body solution.

The specifications of the optimization routine in each phase are shown in Figure 7 as well. In phases (I) and (III), no equality or inequality constraints are imposed, and only a range constraint is set. In phase (II), the inequality constraint is that the Earth’s distance at time t2 must be within a specified tolerance (about 10 km) about the target value of rEP2tar, and the equality constraint is that t2 must coincide with the time of the perigee passage. The range constraint values in the table in Figure 7 are examples.

5.2 Results and discussions

Reference orbits for multiple deep space probes starting from the same specific LEO are generated based on the detailed design method described in the previous section. The results are compared with those of the preliminary design, and the efficiency of paid ΔV to obtain the target V-infinity vector is evaluated in terms of leverage L.

5.2.1 LEO-rideshared asteroid exploration missions

This section provides an overview of the LEO-rideshare asteroid exploration mission group originating in the LEO. The interplanetary missions will be individual rendezvous or flyby missions by four spacecraft targeting a total of four bodies each (3361 Orpheus, 3352 McAuliffe, 2340 Hathor, and 433 Eros). Table 1 shows the target V-infinity vector V2otar at the target departure epoch t2tar for each asteroid. The default LEO departure time t1def is defined as the target departure time t2tar minus one sidereal year (365.25636 days) and is used during the design process and evaluation.

TABLE 1
www.frontiersin.org

TABLE 1. Target bodies and departure conditions of LEO-rideshare interplanetary missions.

This mission is a series of LEO-rideshare asteroid exploration missions where four spacecraft are launched to the LEO from the Tanegashima Space Center (TNSC) in Japan by using the same single launch vehicle. All four use the 1rEFRO + EGA sequence adopted in this study, but depart from the LEO at different epochs each. The launch date to the LEO was set to more than 1 year before t2tar of 3361 Orpheus, which has the earliest target departure epoch. Table 2 summarizes the initial conditions of the LEO. The initial epoch t0 was set considering the possible time in past launch achievements. The initial state (position and velocity vector) on the J2000EQ inertial frame was calculated considering the latitude and longitude of the TKSC launch pad 1 (LP1).

TABLE 2
www.frontiersin.org

TABLE 2. Initial conditions of the LEO parking orbit.

A ΣVEt LEO i-Ω diagram for specific V2tar for each target V2tar at the default LEO departure time t1def of each spacecraft was created, as shown in Figure 8. The propagated LEO for a t1def±10day period is overlayed, and a preliminary evaluation was made for each spacecraft to determine whether it would be feasible to depart. The perigee constraint at EGA was set to 500kmh2πInfkm, as the actual impact on the maximum side was limited. Consequently, all four probes were found to have at least one orbit solution at t1def±7days in terms of the preliminary design.

FIGURE 8
www.frontiersin.org

FIGURE 8. Earth departure conditions for near-Earth asteroids including four target bodies (left) and the ΣVEt LEO i-Ω diagram for specific V2tar at t=t1def of each spacecraft overlayed with LEO propagated for ± 10 days around t1def (right).

5.2.2 Detailed design and validation results

This section presents the results of the detailed design presented in Section 5.1 for the initial conditions and goals established in the previous section. The validity of the proposed method is evaluated by comparing the results of the preliminary design proposed in this study with those of the multi-body design. For the propagation of the LEO, the Earth is assumed as the central body, the 4 × 4 spherical harmonics model is adopted as the Earth’s gravity model, and the Sun and Moon are considered as third-body gravitational perturbations. Atmospheric drag and other perturbations were ignored. For the propagation of interplanetary orbits, the Sun (point mass) is assumed as the central body, and the Earth and the Moon are considered as third-body perturbations. Solar-radiation pressure and other perturbations were ignored.

First, the spacecraft positions and orbital elements at several candidate epochs of the LEO’s departure time t1 are shown in Figure 9 for each spacecraft with respect to the J2000EQ inertial frame. Each marker indicates the candidate’s departure time, and a filled marker is an EGA-assured candidate as a result of the preliminary design. The selected epochs as the output of the detailed design are indicated by the red line. As shown in the figure, there are two candidate departures per single orbital period of LEO, and one candidate per orbital period for SC1 to SC3 and the two candidates for SC4 are determined to be preliminarily EGA-assured. The basis of this decision is shown in the V diagram in the right-hand side of Figure 9. In addition to the spacecraft’s position at each t1 candidate and the direction (α, δ) of V1o resulting from the departure, the perigee altitude h2π expected for the transition to V2tar at time t2 is plotted along (α, δ). Because the V-infinity vector V1o, indicated by the circled markers, coincides with V2i by definition, the EGA is valid (filled-circle markers) if the circled marker is not within the grey shaded area.

FIGURE 9
www.frontiersin.org

FIGURE 9. States on LEO during the first to last t1 candidate epochs for each spacecraft (left) and V diagrams at t1 and t2 overlayed for each spacecraft as the results of the preliminary design (right).

As a result of the detailed design including the multi-body design phases, the resulting epochs are shown in Table 3. For the final selected orbit solution, t1, tdsm, and t2 obtained in each design phase are listed as UTC times. The resulting delta-V is summarized in Table 4. At t1, the amount of delta-V to achieve C3 = 0 is defined as V1esc, and the contribution to V is defined as V1dep. The total amount of delta-V at t1 is denoted as V1total=V1esc+V1dep. The delta-V of VDSM and V2π are shown by each component in the J2000EQ inertial frame, respectively; assuming that both delta-V are not divided into each axis, the norms VDSM and V2π are tabulated. The Vtotal representing the total amount of ∆V required by the spacecraft is defined by the following Eq. 28. According to Eq. 4 and Eq. 5, the actual LVtotal and the ideal LVtotal values of leverage were calculated.

VtotalV1total+VDSM+V2π=V1esc+V1dep+VDSM+V2π(28)

TABLE 3
www.frontiersin.org

TABLE 3. Detailed trajectory design results: Epochs.

TABLE 4
www.frontiersin.org

TABLE 4. Detailed trajectory design results: delta-V and leverage values.

First, it must be noted that the time t1 in Step 2 is within ±10 s in all four cases compared to the values obtained in the preliminary design. This indicates that the preliminary-designed t1 has a higher accuracy of less than 0.2% with respect to the LEO’s orbital period of about 90 min, and is useful as initial values for estimation in the detailed design.

Next, for the leverage value, the preliminary design assumes an ideal situation where the desired V2tar is obtained by a single tangential impulsive ∆V at time t1, in which case, the leverage value is Lideal. To adjust the EGA condition at t2 to obtain V2tar, the detailed design trajectory adds a total of six degrees of freedom, VDSM, and V2π, and the leverage value is expected to be below the ideal value. Consequently, the amount of delta-V required to adjust the EGA conditions VDSM and V2π totaled below 60 m/s, and as a ratio compared to Vtotal, the value was less than 1.5%. Its efficiency is good, and the ratio of the value LVtotal to the ideal leverage value Lideal was about 90%. The result is rather close to the ideal value, but there are still some possibilities to decrease the amount of delta-V such as optimizing the departure delta-V amount and direction at time t1, or by splitting VDSM into multiple opportunities.

Finally, for time t2, it must be noted that the solutions displayed in Table 3 are all within ±TLEO±90min, furthermore, they are within ±TLEO/2±45min. This result is expected, as there is at least one favorable transition to 1rEFRO during one orbital period of LEO. Therefore, as described in Section 4.3, the detailed orbit design confirms that the LEO has a high level of flexibility in terms of the target departure time constraint.

The orbit solutions obtained in the detailed design are shown in Figure 10. The left figure shows the orbit diagram expressed in the J2000EC inertial frame, and the right figure shows the diagram expressed in the Sun–Earth fixed frame. The arrow in the figure indicates the direction of DSM. As shown on the right figure, the four orbits are all inbound orbits. The left figure shows that the LEO departure position and the EGA position are almost the same, which means that the interplanetary orbit has almost the same orbital period with respect to the Earth’s. This means that the prescribed V is appropriately not allocated to the semi-major axis that affects the orbital period, but to the eccentricity and the inclination with respect to the ecliptic plane.

FIGURE 10
www.frontiersin.org

FIGURE 10. Designed trajectories in the J2000EC frame (left), Sun–Earth–Fixed frame (center), and the differences of optimized time instances with respect to the references (top-right and bottom-right).

Finally, the right-hand side of Figure 10 and Figure 11 show all the results of the t1 candidates. The vertical axes of the right-top and right-bottom plots in Figure 10 show the difference of t2 with respect to time t2tar and the difference of t1 in Step 2 with respect to t1 of the preliminary design, respectively. The horizontal axis shows the difference of t1 with respect to t1def [hours]. In Figure 11, each delta-V quantity (Vtotal; V1esc; V1dep; VDSM, V2π) is shown on the vertical axes, with the same horizontal axis. The cases selected in Table 3 and Table 4 are highlighted in both figures as either vertical or horizontal lines. First, the right-top plot in Figure 10 shows that the cases are selected by the absolute difference of t2 with respect to t2tar being the smallest. The right-bottom plot also shows that the differences between the t1 values from the preliminary design and the Step 2 design stayed within ±20 s for all candidates, which reiterates the usefulness of the preliminary design method. In addition, the manual adjustment of tmid in the preliminary design to bring t2 closer to t2tar in Step1 resulted in the plots that t1t1def seems to be inversely proportional to V. Qualitatively, it can be understood that a smaller V results in a larger Earth gravity effect as a third object perturbation during the interplanetary flight; however, we leave the quantitative discussion and prediction as a subject for future work.

FIGURE 11
www.frontiersin.org

FIGURE 11. Obtained delta-V solutions and selected solutions for each spacecraft.

6 Conclusion

To realize a flexible and economical Earth departure from LEO, this study adopts an orbital sequence 1rEFRO + EGA to separate the increment of velocity and the directional change under the assumption of impulsive delta-V. A hypothesis was made that the orbital sequence would enable the Earth departure to interplanetary space with a high level of flexibility and economic efficiency defined in Section 2, even departing from the initial parking orbit of the LEO, which was evaluated throughout this study.

In Section 3, a preliminary design method based on unperturbed two-body dynamics is proposed, and the necessary procedures and relations are summarized. The transition from a particular LEO to 1rEFRO is reduced to a geometric problem of finding the intersection of a ring and a plane, and a discriminant equation determining the number of intersections and an analytical expression for V that enables the transition to 1rEFRO is derived. The proposed method constructs a workflow to reliably identify the maximum two directions that can transition from LEO to 1rEFRO.

In Section 4, based on the preliminary design methodology developed in the previous section, the coverage ccov was defined as an indicator of the flexibility of the Earth departure from the LEO. The procedure to calculate the value from the spherical cap area of the steradians of V2o restriction is presented. The coverage is a function of the LEO’s (i, Ω), target V, and the departure epoch t from LEO. The coverage is plotted over the defined regions of (i, Ω) with given V and t. The plot is referred to as “ΣVEt LEO i-Ω Diagram” and is proposed for the first time in this study with the aid of the preliminary orbit design. In addition to the coverage ccov, plotting the number of feasible orbital solutions to a target V2tar allows for a simplified and preliminary investigation of the feasibility of the Earth departure. According to the ΣVEt LEO i-Ω diagram created with 2.5V6km/s and the perigee altitude constraint set to 500hπ1e6km, the coverage value for most of the LEOs, excluding some polar orbits, is higher than 85%. This indicates that an extremely wide range of LEOs can realize a flexible and economical Earth departure as parking orbits together with the 1rEFRO + EGA sequence.

In this section, using the preliminary orbit design results established in Section 3 as the initial solution, a detailed orbit design was designed in three phases by nonlinear numerical optimization based on a multi-body dynamics model to evaluate the usefulness of the preliminary design method. First, an LEO-rideshare mission of four asteroid explorers was set as an example, and it was confirmed that the preliminary orbit design was valid through the ΣVEt LEO i-Ω diagram. The detailed design results showed that the LEO departure epoch t1 of all four spacecraft remained within ±10 s between the preliminary and detailed designs. This is less than 0.2% of the LEO orbital period, indicating that the preliminary design solution is useful as an initial guess for the detailed design. Furthermore, the leverage LVtotal calculated from the total ∆V amount including DSM is all LVtotal>4.4, which is larger than about 90% of the ideal value Lideal. Finally, it was also confirmed that the EGA epoch t2 fell within the target Earth departure epoch t2tar±45 min, which possesses a high level of flexibility in terms of the Earth departure epoch.

In conclusion, this study revealed that LEO has both flexibility and economic efficiency as a parking orbit for deep space missions by adopting the 1rEFRO + EGA sequence. Based on these results, the study of a standardized deep space orbital transportation architecture starting from LEO, which is expected to significantly reduce the unit cost per launch weight in the future, can be accelerated. Also, possibilities of alternative sequences including other orbital manipulations such as powered EGA and Lunar gravity assist shall be investigated as future work.

Data availability statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Author contributions

The work was conceptualized and developed as a collaboration among all authors. The manuscript was prepared by YTa, and the authors TS and YTs reviewed the manuscript.

Acknowledgments

The authors thank Dr. Roger Gutierrez Ramon and Mr. Manel Caballero for sharing the dataset of Earth departure conditions for rendezvous missions to various near-Earth asteroids.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Ikenaga, T., Utashima, M., Ishii, N., Kawakatsu, Y., and Yoshikawa, M. (2015). Interplanetary parking method and its applications. Acta Astronaut. 116, 271–281. doi:10.1016/j.actaastro.2015.07.018

CrossRef Full Text | Google Scholar

Ito, D., Pushparaj, N., and Kawakatsu, Y. (2022). Expanding interplanetary transfer opportunities from geostationary transfer orbits via Earth synchronous orbits. J. Spacecr. Rockets, 1–16. doi:10.2514/1.a35270

CrossRef Full Text | Google Scholar

Kawaguchi, J., Kuninaka, H., Fujiwara, A., and Uesugi, T. (2006). MUSES-C, its launch and early orbit operations. Acta Astronaut. 59 (8–11), 669–678. doi:10.1016/j.actaastro.2005.07.002

CrossRef Full Text | Google Scholar

Kawakatsu, Y. (2009). “v∞ direction diagram and its application to swingby design,” in Proceedings of 21st International Symposium on Space Flight Dynamics, Toulouse, France, September 28-October 2, 2009.

Google Scholar

Ozaki, N., Pushparaj, N., Takeishi, N., Hyodo, R., and Hyodo, R. (2022). Asteroid flyby cycler trajectory design using deep neural networks. J. Guid. Control, Dyn. 45 (8), 1496–1511. doi:10.2514/1.G006487

CrossRef Full Text | Google Scholar

Russell, R. P., and Ocampo, C. A. (2005). Geometric analysis of free-return trajectories following a gravity-assisted flyby. J. Spacecr. Rockets 42 (1), 138–152. doi:10.2514/1.5571

CrossRef Full Text | Google Scholar

Strange, N., Russell, R., and Buffington, B. (2007). “Mapping the V-infinity globe,” in Proceedings of AAS Astrodynamics Specialist Conference (Michigan, USA, August 19-23, 2007, : Mackinac Island), 07–277.

Google Scholar

Takahashi, S., Ogawa, N., and Kawakatsu, Y. (2019). General characteristics of free-return ensured orbit insertion and trajectory design with MOI robustness in MMX mission. Trans. Jpn. Soc. Aeronautical Space Sci. Aerosp. Technol. Jpn. 17 (3), 404–411. doi:10.2322/tastj.17.404

CrossRef Full Text | Google Scholar

Tsuda, Y., Nakazawa, S., Kushiki, K., Yoshikawa, M., Kuninaka, H., and Watanabe, S. (2016). Flight status of robotic asteroid sample return mission Hayabusa2. Acta Astronaut. 127, 702–709. doi:10.1016/j.actaastro.2016.01.027

CrossRef Full Text | Google Scholar

Keywords: LEO (low Earth orbit), deep space missions, parking orbit, free-return, Earth gravity assist, V∞ diagram, LEO i-Ω diagram

Citation: Takei Y, Saiki T and Tsuda Y (2022) Redefining low Earth orbit as a parking orbit for flexible and economical Earth departure in deep space missions. Front. Space Technol. 3:1049262. doi: 10.3389/frspt.2022.1049262

Received: 20 September 2022; Accepted: 15 November 2022;
Published: 08 December 2022.

Edited by:

Xiaofeng Wu, The University of Sydney, Australia

Reviewed by:

Joseph N. Pelton, International Space University, United States
Yang Yang, University of New South Wales, Australia

Copyright © 2022 Takei, Saiki and Tsuda. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yuto Takei, takei.yuto@jaxa.jp

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.