Notices
Results 1 to 2 of 2

Thread: Symplectic Integration

  1. #1 Symplectic Integration 
    Forum Radioactive Isotope MagiMaster's Avatar
    Join Date
    Jul 2006
    Posts
    3,440
    Is anyone here familiar with symplectic numerical integration or Hamiltonian mechanics? I'm looking into numerical integration schemes for orbital mechanics and it seems that symplectic integrators are needed to keep the system from flying apart in long running simulations.

    So far, I know that the Verlet scheme is symplectic, but in looking for higher order schemes the only one I found (the Forest-Ruth algorithm) requires recomputing the forces at three different positions for each time step, which might be a bit much to do real time.

    I was looking a little closer at the Verlet scheme though, and I can cancel out the fourth order error term by using the accelerations from previous time steps. What I don't know how to check though is whether or not this modified scheme is still symplectic. (I suspect it isn't or I would think someone else would have mentioned it. That, or maybe there's something else wrong with it.)

    The original Verlet scheme is:
    The modified scheme is:


    Reply With Quote  
     

  2.  
     

  3. #2  
    Forum Radioactive Isotope MagiMaster's Avatar
    Join Date
    Jul 2006
    Posts
    3,440
    I guess no one's familiar enough with symplectic integration, but can anyone familiar with Hamiltonian mechanics and/or tensors help me work through this?

    According to Wikipedia, Hamilton's equations are:


    Where q is the position and p is the momentum.

    I've also got that , and that is the gradient of the potential energy. My first question is how to apply that definition to Newtonian gravity? (T is pretty straightforward, but I'm having trouble working out V.)

    Second, it says that the time evolution of Hamilton's equations is a symplectomorphism, meaning it conserves the two-form . As best I can tell, this is a tensor and I'm more or less completely unfamiliar with tensors. (A symplectic integrator is one that also conserves this two-form, which gives it certain good global properties.)


    Reply With Quote  
     

Similar Threads

  1. Integration
    By thyristor in forum Mathematics
    Replies: 5
    Last Post: February 14th, 2011, 10:33 AM
  2. Integration of x^x?
    By m84uily in forum Mathematics
    Replies: 6
    Last Post: November 20th, 2010, 12:42 PM
  3. More integration
    By Heinsbergrelatz in forum Mathematics
    Replies: 3
    Last Post: May 4th, 2010, 11:14 PM
  4. Integration
    By Heinsbergrelatz in forum Mathematics
    Replies: 22
    Last Post: April 3rd, 2010, 01:54 PM
  5. Need help with integration
    By 11rdc11 in forum Mathematics
    Replies: 8
    Last Post: August 1st, 2008, 03:46 PM
Bookmarks
Bookmarks
Posting Permissions
  • You may not post new threads
  • You may not post replies
  • You may not post attachments
  • You may not edit your posts
  •