Can the Velocity Verlet Algorithm Prove Energy Conservation?

  • Context: Graduate 
  • Thread starter Thread starter whatisreality
  • Start date Start date
  • Tags Tags
    Algorithm Velocity
Click For Summary
SUMMARY

The Velocity Verlet Algorithm is confirmed to conserve energy in systems that inherently conserve energy, particularly in orbital mechanics applications. While it does not provide a straightforward proof, it is known to approximately conserve total energy and preserve total angular momentum. The properties of the propagation matrix of the discretized system, especially when analyzed through the Hamiltonian, offer significant insights into energy conservation. A key reference for further understanding is the paper by Hairer, Lubich, and Wanner, specifically chapter 5, which addresses energy conservation in detail.

PREREQUISITES
  • Understanding of the Velocity Verlet Algorithm
  • Familiarity with Hamiltonian mechanics
  • Knowledge of angular momentum conservation principles
  • Basic proficiency in numerical integration techniques
NEXT STEPS
  • Read the paper by Hairer, Lubich, and Wanner on energy conservation in numerical methods
  • Explore the mathematical foundations of the Hamiltonian in classical mechanics
  • Investigate the applications of the Velocity Verlet Algorithm in molecular dynamics simulations
  • Learn about the stability and accuracy of numerical integration methods in physics
USEFUL FOR

Researchers in computational physics, students studying classical mechanics, and professionals working on simulations in orbital mechanics or molecular dynamics will benefit from this discussion.

whatisreality
Messages
286
Reaction score
1
I've read that this algorithm conserves energy if the system it's applied to conserves energy. I can't find a proof, and it's not a particularly obvious statement, so how would you prove it?
 
Physics news on Phys.org
Actually, Verlet integration is approximately conserving the (total) energy, it does however preserve total angular momentum. This is the reason it is used for orbital mechanics problems. The properties of the propagation matrix of the discretized system can already give you a lot of information, especially if you do it for the hamiltonian. There is a very nice paper by Hairer, Lubich and Wanner on this topic, it is still on my 'have-to-really-read' list,

http://citeseerx.ist.psu.edu/viewdo...CA0E3FCC9?doi=10.1.1.7.7106&rep=rep1&type=pdf

(chapter 5 is on energy conservation)
 
  • Like
Likes   Reactions: DrClaude

Similar threads

  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 30 ·
2
Replies
30
Views
3K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 15 ·
Replies
15
Views
3K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 4 ·
Replies
4
Views
1K
  • · Replies 9 ·
Replies
9
Views
3K
  • · Replies 9 ·
Replies
9
Views
2K