Linearization-preserving self-adjoint and symplectic integrators

Springer Science and Business Media LLC - Tập 49 - Trang 177-197 - 2009
R. I. McLachlan1, G. R. W. Quispel2, P. S. P. Tse3
1IFS, Massey University, Palmerston North, New Zealand
2Department of Mathematics, Centre of Excellence for Mathematics and Statistics of Complex Systems, La Trobe University, Melbourne, Australia
3LSEC, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing, China

Tóm tắt

This article is concerned with geometric integrators which are linearization-preserving, i.e. numerical integrators which preserve the exact linearization at every fixed point of an arbitrary system of ODEs. For a canonical Hamiltonian system, we propose a new symplectic and self-adjoint B-series method which is also linearization-preserving. In a similar fashion, we show that it is possible to construct a self-adjoint and linearization-preserving B-series method for an arbitrary system of ODEs. Some numerical experiments on Hamiltonian ODEs are presented to test the behaviour of both proposed methods.

Tài liệu tham khảo

Arnold, V.I.: Mathematical Methods of Classical Mechanics. Springer, New York (1989) Barnett, S.: Matrices: Methods and Applications. Clarendon Press, Oxford (1990) Berland, H., Skaflestad, B., Wright, W.: EXPINT—A Matlab package for exponential integrators. Tech. Report 2005/4, Norwegian University of Science and Technology, Norway (2005) Budd, C.J., Piggott, M.D.: Geometric integration and its applications. In: Handbook of Numerical Analysis, vol. XI, pp. 35–139. North-Holland, Amsterdam (2003) Butcher, J.C.: Numerical Methods for Ordinary Differential Equations. Wiley, New York (2008) Channel, P.J., Scovel, C.: Symplectic integration of Hamiltonian systems. Nonlinearity 3, 231–259 (1990) Chartier, P., Faou, E., Murua, A.: An algebraic approach to invariant preserving integrators: the case of quadratic and Hamiltonian invariants. Numer. Math. 103, 575–590 (2006) Chartier, P., Hairer, E., Vilmart, G.: A substitution law for B-series vector fields. Tech. Report, 2005/5498, INRIA, France (2005) Chartier, P., Hairer, E., Vilmart, G.: Numerical integrators based on modified differential equations. Math. Comput. 76, 1941–1953 (2007) Dellnitz, M., Melbourne, I., Marsden, J.E.: Generic bifurcation of Hamiltonian vector fields with symmetry. Nonlinearity 5, 979–996 (1992) Feng, K.: Difference schemes for Hamiltonian formalism and symplectic geometry. J. Comput. Math. 4, 279–289 (1986) Hairer, E., Lubich, C., Wanner, G.: Geometric Numerical Integration. Structure-Preserving Algorithms for Ordinary Differential Equations. Springer, Berlin (2006) Hairer, E., Wanner, G.: Solving Ordinary Differential Equations. II Stiff and Differential-Algebraic Problems. Springer, Berlin (1996) Hénon, M., Heiles, C.: The applicability of the third integral of motion: Some numerical experiments. Astron. J. 69, 73–79 (1964) Hochbruck, M., Lubich, C., Selhofer, H.: Exponential integrators for large systems of differential equations. SIAM J. Sci. Comput. 19, 1552–1574 (1998) Leimkuhler, B., Reich, S.: Simulating Hamiltonian Dynamics. Cambridge University Press, Cambridge (2004) McLachlan, R.I., Quispel, G.R.W.: Geometric integrators for ODEs. J. Phys. A 39, 5251–5285 (2006) McLaren, D.I., Quispel, G.R.W., Wright, W.M.: Private communication Minchev, B.V., Wright, W.M.: A review of exponential integrators for first order semi-linear problems. Tech. Report, 2005/4, The Norwegian University of Science and Technology, Norway (2005) Sanz-Serna, J.M., Calvo, M.P.: Numerical Hamiltonian Problems. Chapman & Hall, London (1994) Tse, P.S.P.: Geometric Numerical Integration: On the Numerical Preservation of Multiple Geometric Properties for Ordinary Differential Equations. Ph.D. Thesis, Department of Mathematics, La Trobe University, Australia (2007) van der Meer, J.C.: The Hamiltonian Hopf Bifurcation. Springer, Berlin (1985)