A fast, accurate, and smooth planetary ephemeris retrieval system

Precise trajectory simulations typically require an ephemeris retrieval system, i.e. some mechanism to identify planetary body states and orientations at given times. However, the ephemeris systems most commonly used throughout industry and academia are, by design, general in their capabilities and application. Here, we introduce a new system called FIRE (Fast Interpolated Runtime Ephemeris) that is designed for custom trajectory applications that favor speed and smooth derivatives. The new system minimizes the overhead associated with ephemeris calls through the use of archived splines, a runtime ephemeris (stored in random access memory of the computer), and batch processing routines. Further, our approach naturally provides first and second time derivatives for a small additional computational cost. The derivative capability is particularly attractive for optimization and targeting where smooth and accurate derivatives are important. Relative performance comparisons with the Jet Propulsion Laboratory’s Spacecraft Planet Instrument C-matrix Events ephemeris system show typical speed improvements of approximately two orders of magnitude (250 times) for various state and orientation calls. Performance comparisons for high fidelity trajectory propagations are also considered and a factor of 70 in performance increase is achieved for typical cases. The new tool has potential value to any high precision application or software requiring fast, accurate, and smooth ephemeris data.

[1]  David G. Hull,et al.  Numerical Derivatives for Parameter Optimization , 1979 .

[2]  Ryan P. Russell,et al.  Primer Vector Theory Applied to Global Low-Thrust Trade Studies , 2006 .

[3]  Eric T. Campbell,et al.  PRELIMINARY DESIGN OF FEASIBLE ATHOS INTERCEPT TRAJECTORIES , 2004 .

[4]  Giacomo Giampieri,et al.  Rotation rate of Saturn's interior from magnetic field observations , 2004 .

[5]  Martin W. Lo,et al.  The role of invariant manifolds in lowthrust trajectory design (part III) , 2004 .

[6]  Robert A. Jacobson,et al.  Orbits and masses of Saturn's coorbital and F-ring shepherding satellites , 2004 .

[7]  J. Junkins,et al.  HIGHER-ORDER CAYLEY TRANSFORMS WITH APPLICATIONS TO ATTITUDE REPRESENTATIONS , 1997 .

[8]  D. Tholen,et al.  Report of the IAU/IAG Working Group on cartographic coordinates and rotational elements: 2006 , 2007 .

[9]  Ryan P. Russell,et al.  Designing Ephemeris Capture Trajectories at Europa Using Unstable Periodic Orbits , 2007 .

[10]  M. Shuster A survey of attitude representation , 1993 .

[11]  C. H. Acton,et al.  Ancillary data services of NASA's Navigation and Ancillary Information Facility , 1996 .

[12]  Eric Rignot,et al.  Performance analysis and application of grid interpolation techniques for fluid flows , 1993, Experiments in Fluids.

[13]  Anil N. Hirani,et al.  Characteristics of transfers to and captures at Europa , 2006 .

[14]  Carl de Boor,et al.  A Practical Guide to Splines , 1978, Applied Mathematical Sciences.

[15]  Michael Unser,et al.  Splines: a perfect fit for signal and image processing , 1999, IEEE Signal Process. Mag..

[16]  Carl G. Sauer MIDAS - Mission design and analysis software for the optimization of ballistic interplanetary trajectories , 1989 .

[17]  Kazuhito Yokoi,et al.  Planning walking patterns for a biped robot , 2001, IEEE Trans. Robotics Autom..

[18]  J. Dormand,et al.  High order embedded Runge-Kutta formulae , 1981 .

[19]  C. Lin,et al.  Formulation and optimization of cubic polynomial joint trajectories for industrial robots , 1983 .

[20]  C. Burrus,et al.  An in-place, in-order prime factor FFT algorithm , 1981 .

[21]  E. Brigham,et al.  The fast Fourier transform and its applications , 1988 .

[22]  J. Aiello,et al.  Numerical investigation of mapping orbits about Jupiter's icy moons , 2005 .

[23]  Ken Shoemake,et al.  Animating rotation with quaternion curves , 1985, SIGGRAPH.

[24]  Prasun N. Desai,et al.  Entry Dispersion Analysis for the Genesis Sample Return Capsule , 1999 .

[25]  Shane D. Ross,et al.  The Genesis Trajectory and Heteroclinic Connections , 1999 .

[26]  Rodney L. Anderson,et al.  Role of Invariant Manifolds in Low-Thrust Trajectory Design , 2004 .

[27]  Jan M. Rabaey,et al.  VLSI design and implementation fuels the signal-processing revolution , 1998 .

[28]  R. Russell Global search and optimization for free-return Earth-Mars cyclers , 2004 .

[29]  Ryan P. Russell,et al.  Long-Lifetime Lunar Repeat Ground Track Orbits , 2007 .

[30]  Peter D. Welch,et al.  The Fast Fourier Transform and Its Applications , 1969 .

[31]  Hsieh S. Hou,et al.  The Fast Hartley Transform Algorithm , 1987, IEEE Transactions on Computers.

[32]  Gregory J. Whiffen,et al.  The Jupiter Icy Moons Orbiter reference trajectory , 2006 .

[33]  ShoemakeKen Animating rotation with quaternion curves , 1985 .

[34]  D. Tholen,et al.  Report of the IAU Working Group on Cartographic Coordinates and Rotational Elements: 2009 , 2011 .

[35]  Ryan P. Russell,et al.  Designing capture trajectories to unstable periodic orbits around Europa , 2006 .