Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Velocity verlet algorithm

  1. Nov 28, 2015 #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?
     
  2. jcsd
  3. Nov 28, 2015 #2
    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)
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook