Files in this item



application/pdf9314922.pdf (4MB)Restricted to U of Illinois
(no description provided)PDF


Title:Canonical integration methods for Hamiltonian dynamical systems
Author(s):Okunbor, Daniel Irowa
Doctoral Committee Chair(s):Skeel, Robert D.
Department / Program:Computer Science
Discipline:Computer Science
Degree Granting Institution:University of Illinois at Urbana-Champaign
Computer Science
Abstract:Hamiltonian systems possess dynamics (e.g., preservation of volume in phase space and symplectic structure) that call for special numerical integrators, namely canonical methods. Recent research in this aspect have shown that canonical numerical integrators may be needed for Hamiltonian systems. In this thesis, we explore possibilities for the construction of integrators that preserve the property of being canonical. We characterize explicit methods as compositions of 1-stage methods. We find coefficients for methods having minimum work per step. We study in more detail the relationship between canonical, symmetric and order conditions. We also analyze linearized stability. We perform numerical experiments to investigate the importance of being canonical. Examining the solution orbits we discover that numerical integrators possess the same characteristics as the flow of the original system. However, for a variety of test problems we demonstrate equally good results for integrations that merely have the right linearized stability property. We also perform numerical experiments to investigate the improvement that higher-order integrators have over lower-order integrators. This seems to be an interesting research problem, since lower-order methods are currently widely in use in real practical simulations. To direct interest to the application of higher-order methods, serious numerical experiments have to be made to find out their efficiency over lower-order methods. We show that for the two-body problem, the error growth in time for variable step size canonical methods is not in general quadratic. The error growth in time in the case of second-order integrators, canonical or non-canonical turns out to be approximately linear.
Issue Date:1993
Rights Information:Copyright 1993 Okunbor, Daniel Irowa
Date Available in IDEALS:2011-05-07
Identifier in Online Catalog:AAI9314922
OCLC Identifier:(UMI)AAI9314922

This item appears in the following Collection(s)

Item Statistics