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

Development of Kinematic Ephemeris Generator for Korea Pathfinder Lunar Orbiter (KPLO)

Min-Sup Song1http://orcid.org/0000-0001-5705-1026, Sang-Young Park1,http://orcid.org/0000-0002-1962-4038, Youngkwang Kim1http://orcid.org/0000-0002-3062-7133, Jo Ryeong Yim2http://orcid.org/0000-0003-1516-4473
1Astrodynamics and Control Lab., Department of Astronomy, Yonsei University, Seoul 03722, Korea
2KPLO Mission and System Team, Korea Aerospace Research Institute, Daejeon 34133, Korea
Corresponding Author Tel: +82-2-2123-5687, E-mail: spark624@yonsei.ac.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: Jul 9, 2020; Revised: Aug 18, 2020; Accepted: Aug 20, 2020

Abstract

This paper presents a kinematic ephemeris generator for Korea Pathfinder Lunar Orbiter (KPLO) and its performance test results. The kinematic ephemeris generator consists of a ground ephemeris compressor and an onboard ephemeris calculator. The ground ephemeris compressor has to compress desired orbit propagation data by using an interpolation method in a ground system. The onboard ephemeris calculator can generate spacecraft ephemeris and the Sun/Moon ephemeris in onboard computer of the KPLO. Among many interpolation methods, polynomial interpolation with uniform node, Chebyshev interpolation, Hermite interpolation are tested for their performances. As a result of the test, it is shown that all the methods have some cases that meet requirements but there are some performance differences. It is also confirmed that, the Chebyshev interpolation shows better performance than other methods for spacecraft ephemeris generation, and the polynomial interpolation with uniform nodes yields good performance for the Sun/Moon ephemeris generation. Based on these results, a Kinematic ephemeris generator is developed for the KPLO mission. Then, the developed ephemeris generator can find an approximating function using interpolation method considering the size and accuracy of the data to be transmitted.

Keywords: lunar space mission; kinematic ephemeris generator; Korea Pathfinder Lunar Orbiter; KPLO

1. INTRODUCTION

Korea Pathfinder Lunar Orbiter (KPLO) is the first lunar spacecraft of Korea which will perform its missions over 100 km above the Moon’s surface in 2022 year (Ju 2017;Lee et al. 2017;Kim & Song 2019). The KPLO has several tasks, such as taking images of lunar surface, investigating lunar magnetic field distribution, radiation and surface compositions in polar orbits of the moon. For these purposes, the KPLO must be able to point to specific locations whenever needed. Furthermore, the solar panel must be oriented toward the Sun and the communication antenna must be directed to the Earth. To perform these orientations, the position of the spacecraft itself, the position of the Sun, and the position of the Moon should be calculated in the spacecraft’s onboard. The computation of the position of these space objects is called the ephemeris calculation. There have been many researches of onboard ephemeris generation for lunar exploration missions. For the Apollo missions from 1969 to 1971, coasting integration routine with Enke’s method was used to calculate the ephemeris of the spacecraft (Savely et al. 1972). In the case of the Sun and Moon ephemeris, an erasable memory was loaded showing 1° error for a year observed from the Earth (Savely et al. 1972). For Lunar Prospector mission in 1998, FDD (Flight Dynamics Division) in ground system sent data of the spacecraft and the Sun/ Moon/Planet ephemeris data (Andolz 1998). For GRAIL (Gravity Recovery And Interior Laboratory) mission in 2012, ephemeris data were sent to the spacecraft in Chebyshev polynomial coefficient form from the ground system (Antreasian et al. 2012). For Chang’E mission from 2007 to 2013, VLBI (Very Long Baseline Interferometry) system was used to get orbit determination data (Wei et al. 2015).

By referring these research, a KPLO ephemeris generator has been developed by using a kinematical method. The kinematic ephemeris generator is made by utilizing the ephemeris compression method of the NASA JPL planetary Ephemerides (Newhall 1988;Folkner et al. 2014). A kinematic ephemeris generator equipped in onboard computer of spacecraft generates spacecraft’s ephemeris at a desired time, by using uploaded information from a ground system in the Earth. Generally, the uploaded data are spacecraft’s orbit propagation data that contain time, position and velocity values in tabular form, which means it only contains information in discrete time. Saving all of this tabular data requires a lot of memory. Therefore, JPL planetary ephemerides provides only the coefficients of the orbit approximation function and the Lagrange multiplier method as an approximation function in order to calculate desired approximated ephemeris when needed. This tool can approximate orbit in the error of 0.5 mm (Newhall 1988). The ephemeris compression method of JPL planetary ephemerides has good performance in terms of accuracy but has a disadvantage of large data budget. In this method, Lagrange multiplier method is used to find the optimal approximate function and it uses information of both the position and velocity at a given node and boundary conditions. Therefore, the size of coefficients data be saved is large because it requires high order approximation function whose order is larger than the number of nodes given. Also, since only the position error is given for the requirements, using velocity information unnecessarily increases the accuracy by increasing the amount of the coefficient.

For the KPLO mission, an important constraint to be considered in developing a kinematic ephemeris generator is the limitation on the size of the ephemeris compression data because of limitation on onboard computer capability. Thus, a kinematic ephemeris generator for the KPLO has to utilizes a method of interpolating the position information only. The main contribution of this study is that a kinematic ephemeris generator algorithm based on very limited conditions for the KPLO is developed and validated. The algorithm consists of a ground ephemeris compressor and an onboard ephemeris calculator. The ground ephemeris compressor compresses spacecraft’s ephemeris and celestial bodies’ ephemeris in a ground system. This algorithm should be interfaced with a flight dynamic software of a ground system for a lunar mission. The onboard ephemeris calculator generates spacecraft’s ephemeris the Sun/Moon ephemeris at a desired time in an onboard computer of the KPLO, by using uploaded compressed information from ground system. The developed algorithm will be used for the KPLO mission.

In Section 2, an architecture of the kinematic ephemeris generator is presented for the KPLO. In section 3, the ground ephemeris compressor in a ground system is explained for spacecraft ephemeris and the Sun/Moon ephemeris. In section 4, the onboard ephemeris calculator for calculating the ephemeris with the compressed ephemeris data is introduced. In section 4, the kinematic ephemeris generator is verified by many numerical simulations. In the section 6, conclusion of this study is made.

2. A KINEMATIC EPHEMERIS GENERATOR FOR THE KPLO

For the KPLO mission, a kinematic ephemeris generator is developed by referring the planetary ephemerides by NASA JPL. The ephemeris compression method of JPL planetary ephemerides has good performance in terms of accuracy but has a disadvantage of large data budget because it requires both the position and velocity and boundary conditions to use high order approximation function. For the KPLO mission, a kinematic ephemeris generator should satisfy the limitation on the size of the ephemeris compression data and a desired accuracy requirement of ephemeris as well. To do this, a kinematic ephemeris generator adopts an interpolation method to find a proper approximation function. If an interpolation method is used to find a proper approximation function, only position information can be utilized. Consequently, lower order of approximation function can be allowed which has less number of coefficients hence, small data size. The spacecraft ephemeris can be generated by using Polynomial interpolation with Uniform node, Chebyshev interpolation and Hermite interpolation, respectively. In this study, after testing their performances of the three interpolation methods, an appropriate interpolation is chosen and validated for the KPLO mission.

Fig. 1 shows the structure of the kinematic ephemeris generator that will be operated at Flight Dynamics System( FDS) in a ground system of KPLO. The kinematic ephemeris generator is largely divided into two parts: a ground ephemeris compressor and an onboard ephemeris calculator. The ground ephemeris compressor in a ground system generates spacecraft ephemeris compression and the Sun/Moon ephemeris compression. The former compresses the ephemeris of orbit propagation data of spacecraft and the latter compresses the ephemeris of JPL planetary ephemerides. Compressed ephemeris data on the ground system are sent to the KPLO and then the onboard ephemeris calculator calculate ephemeris of the spacecraft and the Sun/Moon using these uploaded data.

jass-37-3-199_F1
Fig. 1. Structure of the kinematic ephemeris generator for the KPLO mission. KPLO, Korea Pathfinder Lunar Orbiter.
Download Original Figure
2.1 Ground Ephemeris Compressor

The ground ephemeris compressor compresses the information of spacecraft ephemeris and the Sun/Moon ephemeris. The compression algorithm for spacecraft ephemeris consists of four steps as shown in Fig. 2. The ephemeris compression algorithm requires information of spacecraft ephemeris (position and velocity values) that can be obtained by an orbit propagator. In the first step, it calculates the time corresponding to selected nodes and then calculates the ephemeris of spacecraft at desired time. In the second step, it interpolates given data point to find data value at the node point. In the third step, it normalizes time corresponding to the node. If it were not normalized, calculated coefficients of the approximating function would have very small value, so it would have round off error. In the last step, it computes polynomial coefficients of the interpolation function at the normalized time by using ephemeris calculated in the previous step. The polynomial coefficients of the interpolation function are calculated by using Vandermonde matrix.

jass-37-3-199_F2
Fig. 2. Ephemeris compression algorithm for spacecraft ephemeris.
Download Original Figure

Because spacecraft ephemeris compression algorithm does not have change of time system or coordinate axis, the time system and coordinate axis of the ephemeris compression file are the same as orbit propagation data. Therefore, if the time system and coordinate system of propagation data are same with those of the spacecraft, the process for adjusting is not needed.

However, for planetary ephemeris, JPL planetary ephemerides is used for ephemeris compression, so it is necessary to adjust the time system and coordinate system. For JPL planetary ephemerides, the ephemeris of the celestial body is represented in the ITRF (International Terrestrial Reference Frame), and only the ephemeris of the Moon is calculated in the ECI (Earth Centered Inertial) frame. Because the requirements are given based on the ECI frame, the coordinate system of the Sun ephemeris must be converted to the ECI frame. In addition, for time system, the JPL planetary ephemeris has TDB (Barycentric Dynamics Time) time system, so it should be converted to a spacecraft time system. The compression algorithm for the Sun/Moon ephemeris consists of four steps as shown in Fig. 3. In the first step, it selects nodes on the time axis of the spacecraft time system. Because JPL Planetary ephemerides provide ephemeris of celestials in TBD time system, the time system of the selected nodes should be changed into the TBD time system. In the second step, it calculates ephemeris of the nodes of TBD time system with JPL planetary ephemerides. In the third step, it normalizes selected nodes from step one. The normalized nodes have ephemeris value calculated in the step two. In the last step, it calculates coefficients of the polynomial interpolation function with the same method introduced in spacecraft ephemeris compression algorithm.

jass-37-3-199_F3
Fig. 3. Compression algorithm for the Sun/Moon ephemeris.
Download Original Figure
2.2 Onboard Ephemeris Calculator

Fig. 4 shows the algorithm of the onboard ephemeris calculator for calculating ephemeris with the compressed ephemeris data. The ephemeris compression data contains coefficient sets of the interpolation function and time information of the sets. In this algorithm, ephemeris compression file and onboard time is basically needed for calculating ephemeris at the desired time. In addition, previous time row number is needed to find coefficient set corresponding to the current onboard time without searching all of the input coefficient file. Finding coefficient set algorithm represented in right side of Fig. 4 contains finding previous time steps to cope with critical case of onboard time fault. After finding the coefficient set, the input onboard time should be normalized to evaluate the polynomials. Also velocity can be calculated by differentiate the polynomials.

jass-37-3-199_F4
Fig. 4. Onboard algorithm for calculating ephemeris.
Download Original Figure

When interpolation is used to find an approximating function, approximating function can be set as polynomial, logarithmic, or exponential, depending on the distribution of the given data. The orbit of the spacecraft or the orbit of the celestial body approximates a quadratic curve if there is no external force. However, it is better to use the polynomial function rather than the quadratic curve because the time boundary is not set wide. Orbit propagation result has position and velocity information. A type of interpolation method using only position information is polynomial interpolations. Polynomial interpolation methods can be classified by method of selecting node such as Uniform node, Chebyshev node. The interpolation method using both the position and velocity information is the Hermite interpolation.

An interpolation method is generally used for finding values of points that are located between the given data points. However, in the developed ephemeris generator, interpolation method is used for finding approximating functions by selecting a specific points (nodes) among many ephemeris data points calculated by orbit propagation. In this method, the interpolation error differs depending on the method of selecting nodes. But the increase of the bound values does not always mean the increase of interpolation values. In general, if number of nodes increases, the error of approximating function would be usually reduced. However, for a certain function like interpolation with uniform node, the interpolation error increases at the edge of interpolation interval as number of nodes increases (Atkinson 1989). This phenomenon is called Runge’s phenomenon (Heath 2018). Therefore, the interpolation uniform node works well with low order polynomial but it can have problem when high-order function is used. In order to relieve the Runge’s phenomenon, Chebyshev nodes can be used instead of uniform node because the maximum error bound is the smallest when Chebyshev interpolation is used (Deprit et al. 1979;Stewart 1996). However, the smallest maximum error bound means that it is likely to have the smallest error, not that it always shows the smallest maximum error. The Chebyshev interpolation has also the advantage that it does not show Runge’s phenomenon when high-order function is used for interpolation (Mathews & Fink 1999). The Hermite interpolation method is a method of finding approximate polynomial functions by using both the function value and its derivative value at the node (Atkinson 1989). Normally the orbit propagation results include position and velocity, so it is possible to apply Hermit interpolation because it can use both position and velocity information. Therefore, if the same node is used, the Hermite interpolation can show better performance than interpolation method using only position information. However, if same interpolation order is applied, the Hermite interpolation uses only half nodes. Therefore, this method would show worse performance if the interpolation interval is large.

3. PERFORMANCE ANALYSIS

Because the onboard ephemeris calculator utilizes interpolation with a small number of nodes to find orbit approximation functions, the performance depends on the shape of the orbit and the characteristics of interpolation method used. In this section, the onboard ephemeris calculator is performed with three interpolation methods such as Polynomial interpolation with Uniform node, Chebyshev interpolation and Hermite interpolation. For the KPLO mission, onboard computer must generate ephemeris data every second for the spacecraft and the Sun/Moon. Therefore, the tests were conducted by comparing the orbit propagation data (as true values) with the generated ephemeris every second. The test parameters contain maximum distance error, mean distance error, and standard deviation. The maximum distance error is the largest value of the error calculated in a specific period and mean distance error, standard deviation represents mean value and standard deviation of entire error calculation in a specific period. Data budget is also calculated according to interpolation order. More detailed anaysis can be found in Song (2019).

3.1 Performance Tests of Spacecraft Ephemeris Generation

To test the spacecraft ephemeris generation, ephemeris data which is similar to the mission orbit is required. The mission orbit is a circular polar orbit with 100km altitude. General Mission Analysis Tool (GMAT) is used to get the reference orbit by propagating orbit based on the mission orbit (NASA GSFC 2017). Table 1 shows the conditions of the orbit propagation. The time span of the propagation is set to 20 days because the ground track of the KPLO almost covers all of the lunar surface by this period. The KPLO spacecraft will has an orbit with 100 km altitude, so the spacecraft will move very fast in about 7 km/sec. Hence, the requirements of spacecraft ephemeris generation is 50 m maximum distance error and data budget of 10,000 byte/day. Table 2, Table 3, and Table 4 show the spacecraft ephemeris generations by using Polynomial interpolation with Uniform node, Chebyshev interpolation and Hermite interpolation, respectively. The tables are the results of the tests according to interpolation orders and interpolation intervals of the three interpolation methods. In these tables, the interpolation interval is significantly changed at the 8th order because the data size does not satisfy the requirement when the interpolation order is set to 8th and the interpolation interval is set to 30 minutes. In Table 2, when the interpolation order increases over the 10th order, the maximum error increases because of the Runge’s phenomenon. On the other hand, in Table 3, interpolation error decreases as both the interpolation order and interpolation interval are increased. In the case of Hermite interpolation of Table 4, the overall error is large compared to other methods, because the number of nodes is small and it does not reflect the rapidly changing orbit. Comparing Table 2 with Table 3, it is shown that polynomial interpolation with uniform nodes have less error when order is 6th, and Chebyshev interpolation has better results when order is over 7th. In the case of 7th order and 30 minute interval, uniform node shows better results in mean error but it has higher maximum error. Table 2, Table 3 and Table 4 have some cases that satisfy the requirements (50 m maximum distance error and data budget 10,000 byte/ day). For examples, when the interpolation interval is set to 30 minutes, the interpolation order of 6th and 7th satisfies the requirements. When the interpolation interval is set to 60 minutes, the interpolation order of 9th and more order satisfy the requirements. The orbit (position and velocity) of the KPLO mission changes rapidly, so it will be more stable to use higher order polynomials or very short interpolation intervals. When the same data budget is kept, higher order polynomials yield less mean distance errors. Therefore, the Chebyshev interpolation is suitable for the spacecraft ephemeris generation when the interpolation interval is set to 60 minutes with 10th or higher order interpolation as shown in Table 3. In the actual operation situation, it is best to select the case with the smallest maximum error as long as data budget requirement is satisfied.

Table 1. Orbit propagation conditions
Propagator Type RK89 Orbit element Semi major axis 1,850 km
Accuracy 1.00E-13 Eccentricity 0.01
Step size 1 sec Inclination 90°
Propagation span 20 days Right ascension of ascending node
Dynamic model Lunar model LP-165 165 ×165 Argument of perigee 120°
Coordinate system Lunar center inertial True anomaly
Download Excel Table
Table 2. Test result of Spacecraft Ephemeris Generation by Polynomial interpolation with Uniform node
Interpolation order Interpolation interval (min) Maximum distance error (m) Mean distance error (m) Standard deviation (m) Data budget (byte/day) Requirements satisfaction
6 30 96.95 15.46 14.39 7,296 no
60 3,690.54 919.04 825.95 3,648 no

7 30 39.784 2.21 2.86 8,448 yes
60 750.54 113.62 135.04 4,224 no

8 30 39.11 1.63 2.49 9,600 yes
60 210.83 18.21 25.10 4,800 no

9 30 38.00 1.49 2.61 10,752 no
60 135.19 8.67 13.23 5,376 no
90 972.41 80.95 136.47 3,854 no

10 60 163.40 8.05 14.34 5,952 no
90 573.76 30.25 53.34 3,968 no

11 60 233.82 8.18 16.70 6,528 no
90 489.13 23.16 44.51 4,352 no

12 60 322.72 8.89 21.05 7,104 no
90 799.65 25.81 56.40 4,736 no

13 60 501.86 10.38 28.44 7,680 no
90 1,090.88 31.51 77.75 5,120 no
Download Excel Table
Table 3. Test result of Spacecraft Ephemeris Generation by Chebyshev interpolation
Interpolation order Interpolation interval (min) Maximum distance error (m) Mean distance error (m) Standard deviation (m) Data budget (byte/day) Requirements satisfaction
6 30 153.88 33.75 25.29 7,296 no
60 6,442.82 1,989.35 1,442.55 3,648 no

7 30 32.82 2.99 3.09 8,448 yes
60 923.52 240.02 196.10 4,224 no

8 30 22.02 1.44 1.75 9,600 yes
60 172.02 31.20 28.15 4,800 no

9 30 15.04 1.00 1.24 10,752 no
60 58.26 7.16 6.70 5,376 no
90 642.14 128.07 120.11 3,854 no

10 60 37.66 4.28 4.55 5,952 yes
90 212.94 32.24 28.92 3,968 no

11 60 33.82 3.34 3.54 6,528 yes
90 93.04 11.02 10.39 4,352 no

12 60 29.28 2.57 2.76 7,104 yes
90 59.18 6.78 6.70 4,736 no

13 60 22.98 1.98 2.20 7,680 yes
90 54.01 5.30 5.37 5,120 no
Download Excel Table
Table 4. Test result of Spacecraft Ephemeris Generation by Hermite interpolation
Interpolation order Interpolation interval (min) Maximum distance error (m) Mean distance error (m) Standard deviation (m) Data budget (byte/day) Requirements satisfaction
6 30 140.61 47.98 35.20 7,296 no
60 6,744.78 2,838.89 1,998.27 3,648 no

8 30 15.28 2.47 2.63 9,600 yes
60 155.82 1.43 36.74 4,800 no

10 60 59.34 10.92 12.76 5,952 no
90 220.33 41.40 47.69 4,464 no

12 60 84.38 6.94 11.13 7,104 no
90 195.95 21.74 30.22 4,736 no
Download Excel Table
3.2 Performance Tests of the Sun/Moon Ephemeris Generation

The limitation on angle errors viewed from the Earth is required for the Sun/Moon ephemeris generation. For the Sun ephemeris, maximum angle error requirement is given as 0.01°, and for the Moon ephemeris it is given as 1 arcsec, which is 0.00028°. These requirements of angle errors can be converted to distance errors as 2,610.98 km for the Sun and 1.8785 km for the Moon. Data size requirements are not given because the size is much less than that of spacecraft ephemeris generation.

In the tests, JPL planetary ephemerides DE 430 is used. Table 5, Table 6, Table 7 are the test results of the Sun ephemeris generation by using Polynomial interpolation with Uniform node, Chebyshev interpolation and Hermite interpolation, respectively. The Hermite interpolation shows worse results because it uses only half number of nodes in the same order comparing to other interpolation methods. For all case, polynomial interpolation with uniform node yields the smallest mean error. Except for 8th order case it also has the smallest error in maximum error. Table 8, Table 9, and Table 10 show the errors of the Moon ephemeris generation by using Polynomial interpolation with Uniform node, Chebyshev interpolation and Hermite interpolation, respectively. For the Moon ephemeris generation, the Hermite interpolation has the greatest error, and the Polynomial interpolation with Uniform node in the Table 8 yields the smallest error for all cases. As a result of the performance tests of the Sun/Moon ephemeris generation, there are some cases that satisfy requirements although the interpolation intervals are set to days. In addition, the polynomial interpolation with uniform node shows better performance than the Chebyshev interpolation, because the Sun/Moon ephemeris generation meets requirements in low-order case. Therefore, the polynomial interpolation with uniform node is suitable for the Sun/Moon ephemeris generation.

Table 5. Test result of Sun Ephemeris Generation by Polynomial interpolation with Uniform node
Interpolation order Interpolation interval (days) Maximum distance error (km) Mean distance error (km) Standard deviation (km) Data budget (byte/day) Requirements satisfaction
4 8 72.19 19.62 14.98 390 yes
16 980.22 298.69 220.01 195 yes
32 11,026.53 4,032.19 2,681.76 97.5 no

5 8 7.64 1.46 1.41 480 yes
16 162.76 38.89 33.71 240 yes
32 2,359.21 831.34 561.10 120 yes

6 8 1.40 0.17 0.22 570 yes
16 53.67 7.77 9.80 285 yes
32 1,134.14 300.78 266.81 142.5 yes

7 8 0.30 0.02 0.04 660 yes
16 20.42 2.08 3.24 330 yes
32 637.00 119.58 137.59 165 yes

8 8 0.07 4.02E-3 8.28E-3 750 yes
16 8.48 0.63 1.21 375 yes
32 398.53 62.81 81.42 187.5 yes
Download Excel Table
Table 6. Test result of Sun Ephemeris Generation by Chebyshev interpolation
Interpolation order Interpolation interval (days) Maximum distance error (km) Mean distance error (km) Standard deviation (km) Data budget (byte/day) Requirements satisfaction
4 8 178.88 40.65 40.67 390 yes
16 2,462.02 612.83 597.09 195 yes
32 25,029.17 8,226.13 7,265.34 97.5 no

5 8 12.95 3.10 2.87 480 yes
16 267.41 79.80 65.44 240 yes
32 3,471.43 1,506.53 1,027.78 120 no

6 8 2.43 0.36 0.42 570 yes
16 88.03 15.68 17.11 285 yes
32 1,676.79 547.63 412.48 142.5 yes

7 8 0.37 0.05 0.06 660 yes
16 23.03 4.00 4.62 330 yes
32 612.10 188.20 162.04 165 yes

8 8 0.07 7.88E-3 0.01 750 yes
16 7.81 1.10 1.39 375 yes
32 298.24 80.53 69.93 187.5 yes
Download Excel Table
Table 7. Test result of Sun Ephemeris Generation by Hermite interpolation
Interpolation order Interpolation interval (days) Maximum distance error (km) Mean distance error (km) Standard deviation (km) Data budget (byte/day) Requirements satisfaction
4 8 348.83 95.62 82.33 390 yes
16 4,629.53 1,419.98 1,176.65 195 no
32 45,673.59 19,238.75 13,541.56 97.5 no

6 8 2.86 0.51 0.59 570 yes
16 101.43 22.35 23.95 285 yes
32 1,954.24 781.27 568.38 142.5 yes

8 8 0.09 9.08E-3 0.01 750 yes
16 10.08 1.31 1.88 375 yes
32 418.80 108.25 102.60 187.5 yes
Download Excel Table
Table 8. Test result of Moon Ephemeris Generation by Polynomial interpolation with Uniform node
Interpolation order Interpolation interval (days) Maximum distance error (km) Mean distance error (km) Standard deviation (km) Data budget (byte/day) Requirements satisfaction
4 2 15.80 5.07 3.35 390 no
4 241.20 80.47 52.67 195 no
8 3,614.65 1,242.55 793.25 97.5 no

5 2 0.69 0.13 0.12 480 yes
4 20.63 3.96 3.90 240 no
8 609.80 120.45 115.25 120 no

6 2 0.03 3.74E-3 5.19E-3 570 yes
4 2.00 0.23 0.32 285 no
8 115.35 13.79 18.50 142.5 no

7 2 1.88E-3 1.42E-4 2.46E-4 660 yes
4 0.22 0.02 0.03 330 yes
8 24.50 2.03 3.37 165 no

8 2 1.12E-4 6.02E-6 1.28E-5 750 yes
4 0.03 1.49E-3 3.13E-3 375 yes
8 5.64 0.33 0.68 187.5 no
Download Excel Table
Table 9. Test result of Moon Ephemeris Generation by Chebyshev interpolation
Interpolation order Interpolation interval (days) Maximum distance error (km) Mean distance error (km) Standard deviation (km) Data budget (byte/day) Requirements satisfaction
4 2 39.87 10.58 9.52 390 no
4 599.37 167.38 149.82 195 no
8 8,973.75 2,564.75 2,257.62 97.5 no

5 2 1.20 0.27 0.26 480 yes
4 35.48 8.50 8.12 240 no
8 1,034.67 255.56 235.32 120 no

6 2 0.06 8.27E-3 9.98E-3 570 yes
4 3.49 0.52 0.62 285 no
8 199.84 29.63 34.51 142.5 no

7 2 2.42E-3 3.10E-4 4.08E-4 660 yes
4 0.27 0.04 0.05 330 yes
8 30.21 4.26 5.35 165 no

8 2 1.22E-4 1.23E-5 1.78E-5 750 yes
4 0.03 3.02E-3 4.31E-3 375 yes
8 5.79 0.65 0.90 187.5 no
Download Excel Table
Table 10. Test result of Moon Ephemeris Generation by Hermite interpolation
Interpolation order Interpolation interval (days) Maximum distance error (km) Mean distance error (km) Standard deviation (km) Data budget (byte/day) Requirements satisfaction
4 2 79.50 25.12 19.04 390 no
4 1,185.71 396.15 297.62 195 no
8 17,196.97 5,988.12 4,375.10 97.5 no

6 2 0.07 0.01 0.01 570 yes
4 4.19 0.72 0.88 285 no
8 235.36 42.34 48.95 142.5 no

8 2 1.53E-4 1.46E-5 2.28E-5 750 yes
4 0.03 3.43E-3 5.54E-3 375 yes
8 7.35 0.75 1.17 187.5 no
Download Excel Table

4. CONCLUSIONS

In this study, a kinematic ephemeris generator for the KPLO is developed and its performance tests are conducted. The kinematic ephemeris generator consists of two parts: a ground ephemeris compressor and an onboard ephemeris calculator. To choose a proper method for the ephemeris generator, the polynomial interpolation with uniform node, Chebyshev interpolation and Hermite interpolation methods are investigated. The Chebyshev interpolation works well in high order interpolation so it shows good performance for spacecraft ephemeris generation because the mission orbit of the spacecraft changes rapidly so it needs high order interpolation to meet the requirements in ephemeris error and data budget. The Sun/Moon ephemeris generation does not need fine accuracy so the polynomial interpolation with uniform node is suitable. The developed kinematic ephemeris generator will be loaded as onboard S/W and will be used for actual mission.

ACKNOWLEDGEMENTS

This research was supported by Korea Aerospace Research Institute as a part of the Lunar exploration program (NRF-2017M1A3A 9014638).

References

1.

Andolz FJ, Lunar Prospector Mission Handbook (Lockheed Martin Missiles & Space, Sunnyvale, CA, 1998).

2.

Antreasian PG, Bhat RS, Broschart SB, Chung MK, Criddle KE, et al., The twin GRAIL spacecraft into science formation at the moon, in 23rd International Symposium on Space Flight Dynamics, Pasadena, CA, 29 Oct-2 Nov 2012.

3.

Atkinson KE, An Introduction To Numerical Analysis, 2nd ed. (John Wiley & Sons, New York, NY, 1989).

4.

Deprit A, Pickard H, Poplarchek W, Compression of ephemerides by discrete Chebyshev approximations, NAVIGATION. 26, 1-11 (1979).

5.

Folkner WM, Williams JG, Boggs DH, Park RS, Kuchynka P, The planetary and Lunar Ephemerides DE430 and DE431, IPN Progress Report 42-196 (2014).

6.

Heath MT, Scientific Computing: An Introductory Survey, 2nd ed. (SIAM, Philadelphia, PA, 2018).

7.

Ju G, Korean Pathfinder Lunar Orbiter (KPLO) status update, in Annual Meeting of the Lunar Exploration Analysis Group (LEAG) Conference, Columbia, MD, 10-12 Oct 2017.

8.

Kim YR, Song YJ, Observational arc-length effect on orbit determination for Korea Pathfinder Lunar Orbiter in the earth-moon transfer phase using a sequential estimation, J. Astron. Space Sci. 36, 293-306 (2019).

9.

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).

10.

Mathews JH, Fink KD, Numerical Methods using MATLAB, 3rd ed. (Pearson, New Delhi, India, 1999).

11.

NASA GSFC, General mission analysis tool (GMAT) (2017) [Internet], viewed 2020 Jul 6, available from: https://software. nasa.gov/software/GSC-17177-1

12.

Newhall XX, Numerical representation of planetary ephemeris, Celest. Mech. 45, 305-310 (1988).

13.

Savely RT, Cockrell BF, Pines S, Apollo experience report: onboard navigational and alignment software, NASA Technical Note D-6741 (1972).

14.

Song MS, Development and performance verification of ephemeris generator for Korea pathfinder lunar orbiter, Master Thesis, Yonsei University (2019).

15.

Stewart GW, Afternotes on Numerical Analysis (Society for Industrial and Applied Mathematics, Philadelphia, PA, 1996).

16.

Wei E, Yan W, Jin S, Wei J, Kutoglu H, et al., Contribution of simulated space VLBI to the Chang’E-1 orbit determination and EOPs estimation, Aerosp. Sci. Technol. 46, 256-263 (2015).