Modified Chebyshev-Picard Iteration Methods for Solution of Initial Value Problems

The Journal of the Astronautical Sciences - Tập 59 - Trang 327-351 - 2014
Xiaoli Bai1, John L. Junkins2
1Graduate Research Assistant, Department of Aerospace Engineering, 3141-TAMU
2Regents Professor, Distinguished Professor of Aerospace Engineering, Holder of the Royce E. Wisenbaker ’39 Chair in Engineering, Department of Aerospace Engineering, 3141-TAMU

Tóm tắt

The solution of initial value problems provide the state history for a given dynamic system, for prescribed initial conditions. Existing methods for solving these problems have not been very successful in exploiting parallel computation architectures, mainly because most of the integration methods implemented on parallel machines are only modified versions of forward integration approaches, which are typically poorly suited for parallel computation. This article proposes parallel-structured modified Chebyshev-Picard iteration (MCPI) methods, which iteratively refine estimates of the solutions until the iteration converges. Using Chebyshev polynomials as the orthogonal approximation basis, it is straightforward to distribute the computation of force functions and polynomial coefficients to different processors. A vector–matrix form is introduced that makes the methods computationally efficient. The power of the methods is illustrated through satellite motion propagation problems. Compared with a Runge–Kutta 4–5 forward integration method implemented in MATLAB, the proposed methods generate solutions with improved accuracy as well as several orders of magnitude speedup, even before parallel implementation. Allowing only to integrate position states or perturbation motion achieve further speedup. Parallel realization of the methods is implemented using a graphics processing unit to provide inexpensive parallel computation architecture. Significant further speedup is achieved from the parallel implementation.

Tài liệu tham khảo

PICARD, E. “Sur l’application des mbthodes d’approximations successives B l’8tude de certaines Bquatioiis differentielles ordinaires,” Journal de Mathématiques Pures et Appliquées, Vol. 9, 1893, pp. 217–271.

SINHA, S. and BUTCHER, E. “Symbolic Computation of Fundamental Solution Matrices for Linear Time Periodic Dynamic Systems,” Journal of Sound and Vibration, Vol. 206, No. 1, 1997, pp. 61–85.

NEAL, H. L., COFFEY, S. L., and KNOWLES, S. “Maintaining the Space Object Catalog with Special Perturbations,” presented as paper AAS 97–687 at the AAS/AIAA Spaceflight Mechanics Meeting, 1997.

COFFEY, S., NEAL, H., VISEL, C., and CONOLLY, P. “Demonstration of a Special-Perturbations-Based Catalog in the Naval Space Command System,” presented as paper AAS 98–113 at the AAS/AIAA Spaceflight Mechanics Meeting, 1998.