High Order and Degree Geopotential and Derivatives Computation Based on the Clenshaw Summation

Abstract:

An artificial satellite subject to the attraction of the Earth is disturbed due to nonsphericaldistribution and non-symmetrical Earth mass. This uneven distribution of mass is expressed by the socalledspherical harmonic coefficients of the Earth potential. For a faster computation, the accelerationderived from the potential is obtained by a series expansion in terms of these harmonics, the fullynormalized Legendre polynomials and their derivatives, and several recursions associated with thelongitude, geocentric latitude and altitude of the center of mass of the satellite. This paper analyzesthe detailed aspects of disturbances in artificial satellites, related with the modeling of the Earth'sgravitational potential as well as numerical implementation of a recursive algorithm to calculate theacceleration of the geopotential based on the Clenshaw summation. In general, one uses recursiveequations of high degree and order to calculate the Legendre polynomials in order to obtain fasterprocessing and numerical accuracy. However, the recursions can yield numerical errors at each stepof the recursion so that higher orders and degrees of harmonics, the accumulated error may be quitepronounced. The computational implementation of the algorithm is carried out by a PC computer.With the implementation of this algorithm it is possible to calculate the geopotential accelerationfor different orbits and different situations. Such approach aims at mitigating the numerical problemsarising from use of extended series expansion when computing recursively the Legendre polynomials.Once the favorable numerical properties are proven, the algorithm can be used in the solution ofpractical problems of orbital space mechanics and for the Brazilian Space Mission.

• Periodical: Applied Mechanics and Materials (Volume 706)
• Pages: 191-205
• DOI: https://doi.org/10.4028/www.scientific.net/AMM.706.191
• Online since: December 2014
• Authors: M.C. Zanardi*, H.K. Kuga, N.R. da Silveira, L. Morgan
• Keywords: Artificial Satellite, Clenshaw Summation, Geopotential, Legendre Functions, Recursive Calculation

References

[1] E. Fantino and S. Casotto, Methods of harmonic synthesis for global geopotential models and their first, second and third-order gradients. Journal of Geodesy, 83, 595-619 (2009).

DOI: 10.1007/s00190-008-0275-0

[2] H. K. Kuga, V. Carrara, Fortran and C codes for higher order and degree geopotential and derivatives computation, Proceedings of XVI SBSR Brazilian Symposium on Remote Sensing, Iguassu Falls, Brazil, Apr. 13-18, 2013. P. 2201-2210. Available from: <http: /urlib. net/3ERPFQRTRW34M/3E7GET9>. Access in: 2014, Jan. 10.

[3] J. R. Wertz, Spacecraft Attitude Determination and Control. D. Reidel. Dordrecht. Holanda. (1978).

[4] M. C. Zanardi, H. K. Kuga and N. R. Silveira, Using Clenshaw summation for recursive high order and degree geopotential for space applications. Proccedings of 11th International Conference on Vibration Problems, Z. Dimitrovová et al. (eds. ), 1-10, Lisboa, Portugal (2013).

DOI: 10.4028/www.scientific.net/amm.706.191

[5] N. R. Silveira, Development and implementation of recursive calculation of the highest order and degree geopotential for space applications, Graduation Work, UNESP, Guaratinguetá, Brazil (2012).

[6] S. A. Holmes and W. E. Featherstone, A unified approach to the Clenshaw summation and the recursive computation of very high degree and order normalized associated Legendre functions. Journal of Geodesy, 76, 279-299 (2002).

DOI: 10.1007/s00190-002-0216-2