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.
On the numerical integration of the FPU-like systems
BENETTIN, GIANCARLO;PONNO, ANTONIO
2011
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.
Pubblicazioni consigliate
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.
