Abstract
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.
Similar content being viewed by others
References
GEAR, C. W. “Parallel Methods for Ordinary Differential Equations,” Calcolo, Vol. 25, Mar. 1988, pp. 1–20.
FRANKLIN, M. A. “Parallel Solution of Ordinary Differential Equations,” IEEE Transactions on Computers, Vol. 27, May 1978, pp. 413–420.
MIRANKER, W. L. and LINIGER, W. “Parallel Methods for the Numerical Integration of Ordinary Differential Equations,” Mathematics of Computation, Vol. 21, Jul. 1969, pp. 303-320.
INCE, E. L. Ordinary Differential Equations. New York, NY: Dover Publications, Inc., 1956, pp. 63, 66, 71.
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.
FOX, L. and PARKER, I. B. Chebyshev Polynomials in Numerical Analysis. London: Oxford University Press, 1972, pp. 31, 61, 68.
CLENSHAW, C. W. “The Numerical Solution of Linear Differential Equations in Chebyshev Series,” Mathematical Proceedings of the Cambridge Philosophical Society, Vol. 53:134–149, 1957, pp. 134–149.
SCRATON, R. E. “The Solution of Linear Differential Equations in Chebyshev Series,” The Computer Journal, Vol. 8, No. 1, 1965, pp. 57–61.
WRIGHT, K. “Chebyshev Collocation Methods for Ordinary Differential Equations,” The Computer Journal, Vol. 6, No. 4, 1964, pp. 358–365.
NORTON, H. J. “The Iterative Solution of Non-Linear Ordinary Differential Equations in Chebyshev Series,” The Computer Journal, Vol. 7, No. 2, 1964, pp. 76–85.
URABE, M. “Galerkin’s Procedure for Nonlinear Periodic Systems,” Archive for Rational Mechanics and Analysis, Vol. 20, Jan. 1965, pp. 120–152.
URABE, M. and REITER, A. “Numerical Computation of Nonlinear Forced Oscillations by Galerkin’s Procedure,” Journal of Mathematical Analysis and Application, Vol. 14, No. 1, 1966, pp. 107–140.
CHEN, B., GARCA-BOLSB, R., JDARB, L., and ROSELLB, M. “Chebyshev Polynomial Approximations for Nonlinear Differential Initial Value Problems,” Nonlinear Analysis, Vol. 63, No. 5–7, 2005, pp. 629–637.
CLENSHAW, C. W. and NORTON, H. J. “The Solution of Nonlinear Ordinary Differential Equations in Chebyshev Series,” The Computer Journal, Vol. 6, No. 1, 1963, pp. 88–92.
SHAVER, J. S. Formulation and Evaluation of Parallel Algorithms for the Orbit Determination Problem. Ph.D. Dissertation, Department of Aeronautics and Astronautics, Massachusetts Institute of Technology, Cambridge, MA, Mar. 1980.
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.
FEAGIN, T. and NACOZY, P. “Matrix Formulation of the Picard Method for Parallel Computation,” Celestial Mechanics and Dynamical Astronomy, Vol. 29, Feb. 1983, pp. 107–115.
FUKUSHIMA, T. “Vector Integration of Dynamical Motions by the Picard-Chebyshev Method,” The Astronomical Journal, Vol. 113, Jun. 1997, pp. 2325–2328.
BAI, X. and JUNKINS, J. L. “Solving Initial Value Problems by the Picard-Chebyshev Method with NVIDIA GPUS,” Proceedings of the 20th Spaceflight Mechanics Meeting, San Diego, CA, Feb. 2010.
BAI, X. Modified Chebyshev-Picard Iteration Methods for Solution of Initial Value and Boundary Value Problems. Ph.D. Dissertation, Texas A&M University, College Station, TX, 2010.
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.
COFFEY, S., HEALY, L., and NEAL, H. “Applications of Parallel Processing to Astro-dynamics,” Celestial Mechanics and Dynamical Astronomy, Vol. 66, Mar. 1996, pp. 61–70.
BERRY, M. M. and HEALY, L. M. “Implementation of Gauss-Jackson Integration for Orbit Propagation,” The Journal of the Astronautical Sciences, Vol. 52, Jul. Sep. 2004, pp. 331-357.
SCHAUB, H. and JUNKINS, J. L. Analytical Mechanics of Space Systems. Reston, Virginia: American Institute of Aeronautics and Astronautics, Inc., 1 ed., 2003.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Bai, X., Junkins, J.L. Modified Chebyshev-Picard Iteration Methods for Solution of Initial Value Problems. J of Astronaut Sci 59, 327–351 (2012). https://doi.org/10.1007/s40295-013-0021-6
Published:
Issue Date:
DOI: https://doi.org/10.1007/s40295-013-0021-6