On the numerical integration of FPU-like systems |
| |
Authors: | G Benettin A Ponno |
| |
Institution: | Università di Padova, Dipartimento di Matematica Pura e Applicata, Via Trieste 63, 35131 Padova, Italy |
| |
Abstract: | This paper concerns the numerical integration of systems of harmonic oscillators coupled by nonlinear terms, like the common FPU models. We show that the most used integration algorithm, namely leap-frog, behaves very gently with such models, preserving in a beautiful way some peculiar features which are known to be very important in the dynamics, in particular the “selection rules” which regulate the interaction among normal modes. This explains why leap-frog, in spite of being a low order algorithm, behaves so well, as numerical experimentalists always observed. At the same time, we show how the algorithm can be improved by introducing, at a low cost, a “counterterm” which eliminates the dominant numerical error. |
| |
Keywords: | Fermi-Pasta-Ulam Leap-frog Geometric numerical integration Backward analysis |
本文献已被 ScienceDirect 等数据库收录! |
|