Abstract
This Letter proposes a new high precision energy-preserving integrator for system of oscillatory second-order differential equations q″(t) +Mq(t)=f(q(t)) with a symmetric and positive semi-definite matrix M and f(q)=-∇U(q). The system is equivalent to a separable Hamiltonian system with Hamiltonian H(p,q)=12pTp+12qTMq+U(q). The properties of the new energy-preserving integrator are analyzed. The well-known Fermi-Pasta-Ulam problem is performed numerically to show that the new integrator preserves the energy integral with higher accuracy than Average Vector Field (AVF) method and an energy-preserving collocation method.
| Original language | English |
|---|---|
| Pages (from-to) | 1185-1190 |
| Number of pages | 6 |
| Journal | Physics Letters, Section A: General, Atomic and Solid State Physics |
| Volume | 376 |
| Issue number | 14 |
| DOIs | |
| State | Published - 5 Mar 2012 |
| Externally published | Yes |
Keywords
- Energy-preserving integrators
- Fermi-Pasta-Ulam problem
- Hamiltonian systems
- Oscillatory differential equations
Fingerprint
Dive into the research topics of 'A new high precision energy-preserving integrator for system of oscillatory second-order differential equations'. Together they form a unique fingerprint.Cite this
- APA
- Author
- BIBTEX
- Harvard
- Standard
- RIS
- Vancouver