Journal of Astronomy and Space Sciences
The Korean Space Science Society
Technical Paper

Observational Arc-Length Effect on Orbit Determination for Korea Pathfinder Lunar Orbiter in the Earth-Moon Transfer Phase Using a Sequential Estimation

Young-Rok Kim†http://orcid.org/0000-0002-0862-0146, Young-Joo Song
Korea Aerospace Research Institute, Daejeon 34133, Korea
Corresponding Author Tel: +82-42-870-3523, E-mail: yrkim@kari.re.kr

© The Korean Space Science Society. All rights reserved. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Nov 15, 2019; Revised: Nov 28, 2019; Accepted: Nov 29, 2019

Abstract

In this study, the observational arc-length effect on orbit determination (OD) for the Korea Pathfinder Lunar Orbiter (KPLO) in the Earth-Moon Transfer phase was investigated. For the OD, we employed a sequential estimation using the extended Kalman filter and a fixed-point smoother. The mission periods, comprised between the perigee maneuvers (PM) and the lunar orbit insertion (LOI) maneuver in a 3.5 phasing loop of the KPLO, was the primary target. The total period was divided into three phases: launch–PM1, PM1–PM3, and PM3–LOI. The Doppler and range data obtained from three tracking stations [included in the deep space network (DSN) and Korea Deep Space Antenna (KDSA)] were utilized for the OD. Six arc-length cases (24 hrs, 48 hrs, 60 hrs, 3 days, 4 days, and 5 days) were considered for the arc-length effect investigation. In order to evaluate the OD accuracy, we analyzed the position uncertainties, the precision of orbit overlaps, and the position differences between true and estimated trajectories. The maximum performance of 3-day OD approach was observed in the case of stable flight dynamics operations and robust navigation capability. This study provides a guideline for the flight dynamics operations of the KPLO in the trans-lunar phase.

Keywords: Korea Pathfinder Lunar Orbiter (KPLO); orbit determination; arc length; sequential estimation; trans-lunar

1. INTRODUCTION

Korea’s first lunar mission, the Korea Pathfinder Lunar Orbiter (KPLO), is being progressed by the Korea Aerospace Research Institute (KARI) (Ju et al. 2013). The KPLO program has been recently modified based on development issues and the KPLO will be launched in the middle of 2022. The Earth-Moon Transfer phase of the KPLO consists of a 3.5 phasing loop, which includes two perigee (PM) and Lunar orbit insertion (LOI) maneuvers. Additional details on the KPLO mission and trajectory are reported in the literature (Ju et al. 2013;Choi et al. 2018;Song et al. 2018;Kim et al. 2019).

The flight dynamics system (FDS) of the ground system is essential in ensuring a stable and successful mission planning and operation. The most important elements for space navigation are orbit determination (OD) and state prediction capability. To ensure successful flight dynamics operations, KARI designed a dual-simulation engine, consisting of a commercial off-the-shelf (COTS) simulation engine and a KARI simulation engine (Song et al. 2018). Numerous studies (Song et al. 2014;Song et al. 2016;Bae et al. 2017;Song et al. 2017) have discussed development strategies, requirements, burn performances, and contingency designs for the development and implementation of the KPLO FDS. For the OD, the COTS simulation engine employs a sequential estimation based on an extended Kalman filter (EKF) and a backward smoother. The KARI simulation engine utilizes a batch estimation method that is based on a least-square filter. This engine uses in-house codes from previous flight dynamics heritage and developed prototypes (Lee et al. 2017, Kim et al. 2017). In order to guarantee a stable performance, the COTS simulation engine will be used as a primary tool for the flight dynamics operations of the KPLO. The KARI simulation engine can be used as a main tool for the space exploration after actual performance validation of the KPLO mission.

The navigation performance can vary depending on the OD and flight dynamics operation strategies. The performance of the OD is a critical factor for mission planning, operation, and payload data processing. The length of tracking arcs is an important aspect to consider in order to control the accuracy level of the OD; moreover, the arc length directly affects the OD accuracy and consistency. The analysis of the KPLO mission orbit periods was presented in a previous paper, in which the authors recommended a 48-hr tracking OD strategy (Kim et al. 2018). The choice of a long tracking arc can ensure accurate estimation results for the OD; however, limited data processing resources and validation time, an irregular maneuvers schedule, and operational constraints restrict the OD arc length to several days. The main issue related to the Earth-Moon Transfer phase is navigation performance, which is critical for accomplishing successful maneuvers. The chosen 3.5 phasing loop trajectory includes several PM and LOI burns; therefore, the performance of the OD during the maneuver time is particularly critical. As a result, to ensure the best possible navigation performance, it is necessary to investigate the arc-length effect on the OD during the Earth-Moon Transfer phase.

In previous lunar exploration missions, various OD arc lengths have been tested to achieve the target navigation performances. Kim et al. (2018) summarized the strategies used to select the OD arc lengths for various missions. Because each mission has its own requirements, mission constraints, maneuver plans, mission trajectory, and orbits, it is impossible determine “the best” arc length selection strategy for lunar missions. In order to the optimal OD arc length for a specific mission, it is important to analyze and determine an OD plan through a pre-launch OD analysis. For what concerns the Earth-Moon Transfer phase of most lunar missions, the OD performance is expected to be lower than that of the lunar mission phase; in fact, payload data processing and mission planning (e.g., measuring, imaging, and monitoring) are not commonly included in the objective of trans-lunar navigation. Therefore, it is indispensable to investigate the effect of the OD arc-length selection and determine the optimal arc length for the KPLO OD during the Earth-Moon Transfer phase.

In this study, focused on the KPLO trans-lunar orbits, the OD results were examined using a sequential estimation technique; moreover, the OD performance was estimated according to the length of the tracking arc for daily-OD. During the trans-lunar phase of real mission operations, various maneuver events (e.g., PM, LOI, momentum unloading, and test burns) affect the behavior of the OD; however, in order to examine in detail the arc-length effect, maneuvers and contingency events were ignored in this study.

Section 2 describes the true orbit and the measurement simulation step. The OD setting and strategy (e.g., dynamics and measurement modeling setting, estimation technique, and quality assessment philosophy) are summarized in Section 3. In the same section, various OD arc-length cases are introduced. Section 4 presents the results of the arclength effect investigation for the KPLO OD during the Earth-Moon Transfer phase (consisting of three parts: launch–PM1, PM1–PM3, and PM3–LOI). Finally, Section 5 summarizes the conclusions of this study.

2. MEASUREMENT SIMULATION

A simulated true orbit, presented in a previous paper (Kim et al. 2019), was used for the measurement generation of the KPLO in the trans-lunar phase. Fig. 1 shows the true orbit of the KPLO Earth-Moon Transfer phase in an Earth-centered inertial frame. The hypothetical Earth-Moon Transfer phase was set to include two maneuvers: PM1 and PM3. The translunar phase was programmed launch for December 26, 2020 at 07:22:17 (UTC), PM1 for January 3, 2021 at 05:25:10 (UTC) (by an impulsive burn), PM3 for January 19, 2021 at 17:10:15 (UTC), and the first LOI maneuver (LOI-1) for January 24, 2021 at 16:56:52 (UTC). Additional information about the maneuvers is shown in Table 1. Choi et al. (2018) and Kim et al. (2019) described a particularly detailed KPLO trans-lunar trajectory using a 3.5 phasing loop. The dynamic modeling information for the KPLO Earth-Moon Transfer true trajectory is summarized in Table 2. In order to account for the Earth gravity effect, the EGM2008 model was applied with a degree of 21 and an order of 21 (Pavlis et al. 2012). JPL DE430 was utilized for the determination of planetary ephemeris and perturbations caused by both solar radiation pressure and third bodies (Sun and Moon). The general relativity effect was not considered in the generation of the true trajectory. The Runge-Kutta 7–8 integration method (with a variable step) was employed for numerical integration.

jass-36-4-293_F1
Fig. 1. True trajectory of the KPLO Earth-Moon Transfer phase (Earth-centered inertial frame). KPLO, Korea Pathfinder Lunar Orbiter.
Download Original Figure
Table 1. Maneuver information of the simulated KPLO Earth-Moon Trajectory
jass-36-4-293_T1
Download Excel Table
Table 2. Dynamic model settings used for the true trajectory generation and orbit determination
jass-36-4-293_T2
Download Excel Table

The KPLO measurements were generated considering a three-ground-station configuration and simulating the true trajectory. The pseudo-noise (PN) range and the Doppler tracking data by two deep space network (DSN) antennas and the Korea Deep Space Antenna (KDSA) were simulated at 60-s intervals. The antennas located at Goldstone, Madrid (DSN), and Yeoju (KDSA) were regarded as communication and tracking stations. The continuous tracking duration was restricted to 2 hrs, followed by a 1-hr break. Multiple tracking by the two ground stations was not allowed. The simulated KPLO measurements for the trans-lunar phase are shown in Fig. 2.

jass-36-4-293_F2
Fig. 2. Tracking schedules of the ground stations used for the KPLO measurement simulation (trans-lunar phase). KPLO, Korea Pathfinder Lunar Orbiter.
Download Original Figure

The measurement error statistics were accounted for by using the specifications of the KDSA hardware and the values of the DNS service catalog for the tracking data simulation (JPL 2015). Additionally, we used the performance analysis results for the KPLO communication and tracking obtained through sequential and PN ranging techniques (Park & Moon 2018). In order to considering the transponder delay variation uncertainty and the tracking accuracy, the ranging noises of the DSN and KDSA were set to 13 m and 22 m, respectively. The sigma values of the Doppler measurement noises of the DSN and KDSA were set instead to 0.003 Hz and 0.15 Hz, respectively. The Doppler noise values, originally expressed as velocities (mm/s), were converted into cycles (Hz) to allow software implementation (Woodburn et al. 2015). Space environmental effects (e.g., tropospheric refraction and ionospheric delay), which can bias or distort the range and Doppler results, were considered in the measurement simulation. In particular, in order to determine measurement and time biases, we applied the Gauss-Markov and random walk models, respectively. Both of these models were chosen for their simplicity. The Gauss-Markov model is based on a scalar exponential Gauss-Markov sequence, which is the natural model for the measurement of bias modeling; meanwhile, the random walk model is able to determine the simple evolution characteristics of time biases.

3. ORBIT DETERMINATION

In this section, we demonstrate the strategy adopted for the study of the KPLO OD (sequential estimation for the Earth-Moon Transfer phase) and the configuration of the OD arc-length effect investigation. For the KPLO OD, we utilized a sequential filter based on the EKF and a backward smoother of the ODTK software (Vallado et al. 2010). The mission operation and analysis capability have been determined through sequential estimation techniques derived from the IBEX, LRO, and LADEE missions (Policastri et al. 2009;Slojkowski et al. 2015;Policastri et al. 2015a, b). The COTS simulation engine of the KPLO FDS employs the same sequential estimation methods of the ODTK; therefore, our analysis results can be used to examine the actual performance of the KPLO flight dynamics operation.

The dynamic modeling settings for the KPLO OD in the Earth-Moon Transfer phase are shown in Table 2. The same modeling configuration to establish the true trajectory was applied in the investigation of the arc-length effect. Media corrections (e.g., tropospheric reflection and ionosphere delay) were considered during measurement modeling, parallelly with the plate motion of the ground station and the antenna correction. The measurement and time bias estimations were accomplished during the OD process. The Runge-Kutta 7–8 integration method (with a variable step) was utilized for the numerical integration of the state propagation.

The OD for the KPLO was achieved considering a daily schedule and using a fixed tracking arc length [we used the same OD configuration applied to the trans-lunar phase in Kim et al. (2019)]. The fixed arc length-based daily OD was accomplished; afterward, in order to evaluate the orbit quality, we calculated the overlap period between the two arcs. The estimation parameters included: the state of the spacecraft (i.e., position and velocity), the solar radiation pressure coefficients (Cr), and the transponder’s and station’s biases. The initial orbit position uncertainties for the radial, along-track, and cross-track directions were assumed to be 300 m, 1,000 m, and 200 m, respectively [as in Kim et al. (2019)]. The OD calibrations by fine filter tuning were not included in our OD analysis. For the trans-lunar phase, the length of the OD arc can be shortened due to maneuver schedule and the operation; however, based on our arclength effect examination, we assumed a fixed arc length for OD. Six cases, corresponding to periods of 24 hrs, 48 hrs, 60 hrs, 3 days, 4 days, and 5 days, were considered for the OD arc-length investigation. Hence, we analyzed the OD results obtained for these different arc lengths during the launch– LOI phase. The orbit maneuvers (i.e., PM1 and PM3) were excluded in this study. Therefore, three main phases were finally considered: launch–PM1, PM1–PM3, and PM3–LOI.

In order to evaluate the quality of the KPLO OD orbit, we investigated the orbit uncertainty (by error covariance), the orbit overlap precision, and the external orbit comparison by considering the differences between the true and estimated trajectories. All the results (i.e., uncertainties and differences) were calculated in the radial, along-track, and cross-track directions, and their total 3D position values were used to determine the arc-length effect. Moreover, the mean and standard deviation values of the different phases were used to determine the OD performance, in accord with the arc-length selection strategy. Notably, the orbit uncertainties were calculated using 3 sigma values.

4. RESULTS

In this section, we describe the three phases (i.e., launch–PM1, PM1–PM3, and PM3–LOI) through which we performed the OD, as well as their accuracies with respect to the arc length. The investigation of the arc-length effect on the OD was conducted considering tracking arcs of 5 days, 4 days, 3 days, 60 hrs, 48 hrs, and 24 hrs, for a total of 29 days.

4.1 Launch–PM1 Phase (December 26, 2020 at 08:00:00– January 3, 2021 at 05:25:10)

The OD results relative to the launch–PM1 phase are presented in Figs. 35: the thick solid lines indicate the mean values of the position uncertainties for all cases, accompanied by the standard deviation error bars; the other lines indicate some results selected among all the OD dates. The OD results for all cases are shown in Tables 38. We obtained 14 OD results based on the 24-hr tracking arc for the launch–PM1 phase (Table 2), but only three results based on the 5-day tracking arc OD (Table 8). Since the results obtained for all cases could not be described in one figure, we created several figures; in each of them, we presented only the data corresponding to all the arclength cases having a specific period. A daily schedule was employed for the OD cases with 5-day, 4-day, 3-day, 60-hr, and 48-hr arcs, while a 12-hr schedule was employed for the 24-hr arc OD (for the organization of the orbit overlap period). Table 4, 5, 6, 7

jass-36-4-293_F3
Fig. 3. Total position uncertainty of the orbit determination based on various arc lengths (before PM1). OD, orbit determination; PM, perigee maneuvers.
Download Original Figure
jass-36-4-293_F4
Fig. 4. Orbit determination precision represented by orbit overlaps (before PM1). OD, orbit determination; PM, perigee maneuvers.
Download Original Figure
jass-36-4-293_F5
Fig. 5. Orbit determination accuracy represented by differences between the true and estimated orbits (before PM1). OD, orbit determination; PM, perigee maneuvers.
Download Original Figure
Table 3. Position uncertainties and differences based on the error covariance and true trajectory (24-hr arc length)
jass-36-4-293_T3
Download Excel Table
Table 4. Position uncertainties and differences based on the error covariance and true trajectory (48-hr arc length)
jass-36-4-293_T4
Download Excel Table
Table 5. Position uncertainties and differences based on the error covariance and true trajectory (60-hr arc length)
jass-36-4-293_T5
Download Excel Table
Table 6. Position uncertainties and differences based on the error covariance and true trajectory (3-day arc length)
jass-36-4-293_T6
Download Excel Table
Table 7. Position uncertainties and differences based on the error covariance and true trajectory (4-day arc length)
jass-36-4-293_T7
Download Excel Table
Table 8. Position uncertainties and differences based on the error covariance and true trajectory (5-day arc length)
jass-36-4-293_T8
Download Excel Table

Fig. 3 shows the 3D total position uncertainties of six arclength OD settings: the data show that longer arc lengths corresponded to more precise and stable OD results. The mean position uncertainty values were 615.9 m, 465.3 m, 441.5 m, 421.3 m, 391.5 m, and 357.9 m for the 24-hr, 48-hr, 60-hr, 3-day, 4-day, and 5-day arc-length cases, respectively (Table 9) ; meanwhile, the standard deviation values corresponded to 124.6 m, 61.0 m, 41.8 m, 52.8 m, 38.1 m, and 27.3 m, respectively. For the 24-hr arc OD, the mean and standard deviation values were larger, reflecting a large dispersion.

Table 9. Summary of the OD results (launch–PM1 phase)
jass-36-4-293_T9
Download Excel Table

Fig. 4 and Tables 38 show the orbit overlap differences before PM1. The orbit overlap differences were of 169.6 m, 160.1 m, 137.3 m, 110.6 m, 73.4 m, and 91.5 m for the 24-hr, 48-hr, 60-hr, 3-day, 4-day, and 5-day arc lengths, respectively; meanwhile, the standard deviation values corresponded to 127.5 m, 69.3 m, 49.5 m, 38.1 m, 18.4 m, and 60.0 m, respectively.

Fig. 5 shows the OD accuracy (i.e., the difference between the true and estimated orbits) before PM1. The orbit differences corresponded to 176.7 m, 162.4 m, 159.8 m, 155.9 m, 155.1 m, and 151.1 m in the case of the 24-hr, 48-hr, 60-hr, 3-day, 4-day, and 5-day arc lengths, respectively; meanwhile, the standard deviation values corresponded to 73.7 m, 36.7 m, 36.2 m, 28.0 m, 15.6 m, and 34.8 m, respectively. Figs. 35 and Table 9 show how OD precision and accuracy improved with increasing arc length; moreover, the OD performance did not improve drastically by increasing the arc length over 3 days.

Table 9 demonstrates that the precision and accuracy of the OD before PM1 were of 420 m (based on the 3 sigma values), 110 m, and 156 m (based on a > 3-days OD strategy).

4.2 PM1–PM3 Phase (January 3, 2021 at 05:25:10–January 19, 2021 at 17:10:15)

A lot of OD arcs were investigated for the PM1–PM3 phase. The OD results are shown in Figs. 68; moreover, the results of all arcs for the PM1–PM3 phase are contained in Tables 38. Fig. 6 shows the OD position uncertainties for various arc lengths in the PM1–PM3 phase. The mean values of the six arc-length cases were 594.9 m, 437.0 m, 404.2 m, 386.6 m, 347.7 m, and 312.6 m for the 24-hr, 48-hr, 60-hr, 3-day, 4-day, and 5-day arc-lengths, respectively; moreover, the correspondent standard deviation values were 113.0 m, 67.2 m, 60.8 m, 60.6 m, 55.7 m, and 46.5 m, respectively. Notably, the position precision was improved by increasing the tracking arc length from 24 hrs to 5 days; the standard deviation trend showed similar features. Fig. 7

jass-36-4-293_F6
Fig. 6. Total position uncertainty of the orbit determination based on various arc lengths (PM1–PM3). OD, orbit determination; PM, perigee maneuvers.
Download Original Figure
jass-36-4-293_F7
Fig. 7. Orbit determination precision represented by orbit overlaps (PM1–PM3). OD, orbit determination; PM, perigee maneuvers.
Download Original Figure
jass-36-4-293_F8
Fig. 8. Orbit determination accuracy represented by differences between the true and estimated orbits (PM1–PM3). OD, orbit determination; PM, perigee maneuvers.
Download Original Figure

Figs. 7 and 8 show the OD precision (obtained by considering the orbit overlap) and the OD accuracy (obtained by calculating the differences between the true and determined positions), respectively. The mean orbit overlap differences for the 24-hr, 48-hr, 60-hr, 3-day, 4-day, and 5-day arc lengths were of 205.3 m, 136.2 m, 122.8 m, 10.5 m, 87.7 m, and 90.5 m, respectively; moreover, the correspondent standard deviation values were of 223.0 m, 110.4 m, 64.5 m, 3.6 m, 57.1 m, and 59.8 m, respectively. The results of overlap comparison indicate that both the OD precision and stability increased proportionally with the tracking arc duration. Additionally, we found that the overlap precision did not drastically improve for arc length > 4 days. Similar trends were observed from the orbit difference assessments; however, the orbit accuracy remained almost unchanged when using arc lengths > 48 hrs. The orbit differences for the 24-hr, 48-hr, 60-hr, 3-day, 4-day, 5-day arc lengths were of 211.8 m, 176.5 m, 170.6 m, 175.8 m, 178.3 m, and 189.2 m, respectively; moreover, the correspondent standard deviation values were 172.1 m, 83.1 m, 59.8 m, 64.6 m, 63.3 m, and 60.5 m, respectively. All the results are summarized in Table 10. The OD accuracy was lower in the case of the 5-day tracking strategy than for shorter arc lengths; additionally, we noted that longer arc lengths corresponded always to a better performance. Since the 3.5 phasing loop trans-lunar trajectory has an elliptical orbit, the position of the spacecraft on the phasing loops can affect the OD accuracy. In this case, although the OD tracking data become longer, their accuracy may not be improved due to the tracking data around the perigee or apogee. Therefore, longer arcs are not expected to correspond always to better OD accuracies.

Table 10. Summary of the OD results (PM1–PM3 phase)
jass-36-4-293_T10
Download Excel Table
4.3 PM3–LOI Phase (January 19, 2021 at 17:10:15–January 24, 2021 at 16:56:52)

The number of OD arc cases investigated for the PM3– LOI phase is lower than that for the launch–PM1 and PM1–PM3 phases. As shown in Figs. 911, the OD results showed similar trends to those observed in other periods; however, fewer OD cases than those before PM3 revealed an indistinct tendency.

jass-36-4-293_F9
Fig. 9. Total position uncertainty of the orbit determination based on various arc lengths (PM3–LOI). OD, orbit determination; PM, perigee maneuvers; LOI, lunar orbit insertion.
Download Original Figure

Fig. 9 shows the position uncertainties by error covariance of the OD cases. The detailed results for all arcs are displayed in Tables 38. Only one result was obtained for the 5-day OD arc length. The mean values of the position uncertainties were 604.9 m, 415.5 m, 392.0 m, 376.8 m, 334.0 m, and 280.5 m for the 24-hr, 48-hr, 60-hr, 3-day, 4-day and 5-day arc lengths, respectively; moreover, the standard deviations were 174.2 m, 96.3 m, 89.9 m, 66.7 m, and 46.8 m for 24-hr, 48-hr, 60-hr, 3-day, and 4-day arc lengths, respectively. The standard deviation of the 5-day OD arc case was null, since only one value was available. The total position uncertainties obtained for the 5-day and 24-hr arc lengths indicated the highest and lowest precisions, respectively.

Fig. 10 shows the results of the orbit overlap comparison for five arc-length OD cases (no consecutive arcs were observed in the case of the 5-day arc length). The 24-hr arc OD delivered the worst performance, while the other arc settings delivered results with similar precisions. The mean overlap differences were of 373.8 m, 190.7 m, 171.6 m, 207.5 m, and 235.8 m for the 24-hr, 48-hr, 60-hr, 3-day, and 4-day arc lengths, respectively; moreover, the correspondent standard deviations were 262.1 m, 110.0 m, 115.3 m, 95.8 m, and 10.6 m, respectively.

jass-36-4-293_F10
Fig. 10. Orbit determination precision represented by orbit overlaps (PM3–LOI). OD, orbit determination; PM, perigee maneuvers; LOI, lunar orbit insertion.
Download Original Figure

Fig. 11 shows the orbit accuracy based on the differences between the true and estimated orbits. Similar trends were observed for the OD in the case of increasing arc lengths; however, the differences between the length cases were not significant, except in the case of the 24-hr arc length. The mean position differences were of 280.4 m, 185.9 m, 164.3 m, 187.1 m, 213.7 m, and 171.1 m for the 24-hr, 48-hr, 60- hr, 3-day, and 4-day arc settings, respectively; moreover, the correspondent standard deviations were 177.4 m, 48.0 m, 43.6 m, 74.2 m, and 53.4 m, respectively. The biggest difference was reported in the case of the 24-hr arc length, while the best accuracy was reported in the case of the 60 hr-tracking strategy, as shown in Table 11.

jass-36-4-293_F11
Fig. 11. Orbit determination accuracy represented by differences between the true and estimated orbits (PM3–LOI). OD, orbit determination; PM, perigee maneuvers; LOI, lunar orbit insertion.
Download Original Figure
Table 11. Summary of the OD results (PM3–LOI phase)
jass-36-4-293_T11
Download Excel Table

In this study, we described the KPLO OD results for the Earth-Moon Transfer phase and investigated the arc-length effect on the OD in six cases (arc-lengths of 24 hrs, 48 hrs, 60 hrs, 3 days, 4 days, and 5 days). We discovered that the OD orbit accuracy was higher in the case of longer arcs; however, unexpected trajectory correction maneuvers, abnormal spacecraft operations, and contingency situations cannot guarantee long arc lengths in the trans-lunar phase. Therefore, an OD analysis that considers various tracking arc lengths is useful for the preparation of actual KPLO mission operations in the Earth-Moon Transfer phase. Our analysis of the arc-length effect on the OD led to three main conclusions: 1) a 24-hr arc length is not suitable for the KPLO OD during the trans-lunar phase, since the orbit precision and accuracy values would be too large and unstable; 2) tracking arcs > 3 days can deliver precise and stable OD results for the KPLO mission; 3) the OD analysis should be conducted by applying various assessment methods. We used three different orbit quality checks, finding that their trends were all very similar; however, the OD performance was different in some cases. Therefore, it is necessary to check the OD precision and accuracy using various evaluation methods. Considering the quick processing of the sequential estimation and the mission operation effectiveness, we recommend a 3-day arc length for the KPLO OD during the Earth-Moon Transfer phase.

5. CONCLUSIONS

In this paper, we presented the results of an OD analysis conducted using a sequential estimation technique for the KPLO mission in the Earth-Moon transfer phase; moreover, we discussed the arc-length effect on the OD. The KPLO will be Korea’s first lunar exploration mission. To ensure a stable operation and the success of the mission during the Earth-Moon Transfer phase, it is critical to perform an OD analysis. To conduct such analysis, we utilized Doppler and range measurements tracked by three ground stations (one of the KDSA and two antennas of the DSN); additionally, the OD quality was assessed by analyzing the position uncertainty, the precision of the orbit overlaps, and the position differences between the true and determined orbits. Six OD arc-lengths cases were considered: of 24 hrs, 48 hrs, 60 hrs, 3 days, 4 days, and 5 days. Among those, the arcs > 3 days demonstrated a stable precision and highly accurate results. Considering the characteristics of the sequential estimation approach and effective mission operation, we concluded that a 3-day OD approach would deliver the most effective performance. This strategy is suitable in the case of stable flight dynamics operations and robust navigation capabilities during the trans-lunar phase of the KPLO. Overall, this study provides useful guidelines and new knowledge regarding the KPLO flight dynamics operations in the trans-lunar phase.

ACKNOWLEDGMENTS

This work was supported by the Korea Aerospace Research Institute (KARI); (under the Ministry of Science and ICT) through the “Development of Korea Pathfinder Lunar Orbiter and Key Technologies for the Second Stage Lunar Exploration” project (no. SR19060).

References

1.

Bae J, Song YJ, Kim YR, Kim B, Burn delay analysis of the lunar orbit insertion for Korea Pathfinder Lunar Orbiter, J. Astron. Space Sci. 34, 281-288 (2017).

2.

Choi SJ, Whitely R, Condon G, Loucks M, Park JI, et al., Trajectory design for the Korea Pathfinder Lunar Orbiter (KPLO), Proceedings of AAS/AIAA Astrodynamics Specialist Conference, Snowbird, UT, 19-23 Aug 2018 .

3.

JPL [Jet Propulsion Laboratory], Deep Space Network services catalog, NASA Jet Propulsion Laboratory, DSN No. 820-100 (2015).

4.

Ju G, Bae J, Choi SJ, Lee WB, Lee CJ, New Korean lunar exploration program (KLEP): an introduction to the objectives, approach, architecture, and analytical results, in 64th International Astronautical Congress, Beijing, China, 23-27 Sep 2013.

5.

Kim Y, Park SY, Lee E, Kim M, A deep space orbit determination software: overview and event prediction capability, J. Astron. Space Sci. 34, 139-151 (2017).

6.

Kim YR, Song YJ, Bae J, Choi SW, Observational arc-length effect on orbit determination for KPLO using a sequential estimation technique, J. Astron. Space Sci. 35, 295-308 (2018).

7.

Kim YR, Song YJ, Bae J, Choi SW, Orbit determination simulation for Korea Pathfinder Lunar Orbiter using sequential estimation approach, Proceedings of AAS/AIAA Space Flight Mechanics Meeting, Maui, HI, 13-17 Jan 2019.

8.

Lee E, Kim Y, Kim M, Park SY, Development, demonstration and validation of the deep space orbit determination software using lunar prospector tracking data, J. Astron. Space Sci. 34, 213-223 (2017).

9.

Park S, Moon S, Performance analysis of ranging techniques for the KPLO mission, J. Astron. Space Sci. 35, 39-46 (2018).

10.

Pavlis NK, Holmes SA, Kenyon SC, Factor JK, The development and evaluation of the Earth Gravitational Model 2008 (EGM2008), J. Geophys. Res. 117, B04406 (2012).

11.

Policastri L, Carrico J, Craychee T, Johnson T, Woodburn J, Orbit determination operations for the Interstellar Boundary Explorer, Proceedings of the 19th AAS/AIAA Space Flight Mechanics Meeting, Savannah, GA, 8-12 Feb 2009.

12.

Policastri L, Carrico J, Nickel C, Pre-launch orbit determination design and analysis for the LADEE mission, Proceedings of the 25th AAS/AIAA Space Flight Mechanics Meeting, Williamsburg, VA, 11-15 Jan 2015a.

13.

Policastri L, Carrico J, Nickel C, Kam A, Lebois R, et al., Orbit determination and acquisition for LADEE and LLCD mission operations, Proceedings of the 25th AAS/AIAA Space Flight Mechanics Meeting, Williamsburg, VA, 11-15 January, 2015b.

14.

Slojkowski S, Lowe J, Woodburn J, Orbit determination for the lunar reconnaissance orbiter using an extended Kalman filter, Proceedings of the 25th International Symposium on Space Flight Dynamics, Munich, Germany, 19-23 Oct 2015.

15.

Song YJ, Ahn SI, Sim ES, Development strategy of orbit determination system for Korea’s lunar mission: lessons from ESA, JAXA, ISRO and CNSA’s experiences, J. Astron. Space Sci. 31, 247-264 (2014).

16.

Song YJ, Bae J, Kim YR, Kim BY, Uncertainty requirement analysis for the orbit, attitude, and burn performance of the 1st lunar orbit insertion maneuver, J. Astron. Space Sci. 33, 323-333 (2016).

17.

Song YJ, Bae J, Kim YR, Kim BY, Early phase contingency trajectory design for the failure of the first lunar orbit insertion maneuver: direct recovery options, J. Astron. Space Sci. 34, 331- 342 (2017).

18.

Song YJ, Lee D, Bae J, Kim YR, Choi SJ, Flight dynamics and navigation for planetary missions in Korea: past efforts, recent status, and future preparations, J. Astron. Space Sci. 35, 119-131 (2018).

19.

Vallado DA, Hujsak RS, Johnson TM, Seago JH, Woodburn JW, Orbit determination using ODTK version 6, Proceedings of the 4th International Conference on Astrodynamics Tools and Techniques, Madrid, Spain, 3-6 May 2010.

20.

Woodburn J, Policastri L, Owens B, Generation of simulated tracking data for LADEE operational readiness testing, Proceedings of the 25th AAS/AIAA Space Flight Mechanics Meeting, Williamsburg, VA, 11-15 Jan 2015.