Journal article
Preserving poisson structure and orthogonality in numerical integration of differential equations
Computers & mathematics with applications (1987), Vol.48(1), pp.237-255
2004
DOI: 10.1016/j.camwa.2003.02.013
Abstract
We consider the numerical integration of two types of systems of differential equations. We first consider Hamiltonian systems of differential equations with a Poisson structure. We show that symplectic Runge-Kutta methods preserve this structure when the Poisson tensor is constant. Using nonlinear changes of coordinates this structure can also be preserved for nonconstant Poisson tensors, as exemplified on the Euler equations for the free rigid body. We also consider orthogonal flows and the closely related class of isospectral flows. To numerically preserve the orthogonality property we take the approach of formulating an equivalent system of differential-algebraic equations (DAEs) and of integrating the system with a special combination of a particular class of Runge-Kutta methods. This approach requires only matrix-matrix products and can preserve geometric properties of the flow such as reversibility.
Details
- Title: Subtitle
- Preserving poisson structure and orthogonality in numerical integration of differential equations
- Creators
- L.O Jay - University of Iowa
- Resource Type
- Journal article
- Publication Details
- Computers & mathematics with applications (1987), Vol.48(1), pp.237-255
- DOI
- 10.1016/j.camwa.2003.02.013
- ISSN
- 0898-1221
- eISSN
- 1873-7668
- Publisher
- Elsevier Ltd
- Language
- English
- Date published
- 2004
- Academic Unit
- Mathematics
- Record Identifier
- 9984240765002771
Metrics
12 Record Views