1. Not finding help here? Sign up for a free 30min tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Hamiltonian formulation and the Kepler problem

  1. Oct 19, 2016 #1
    This isn't exactly homework, but something which you'd get as an assignment, so I'll still post it here in order to reach the right people..

    I'm attempting to freshen up my knowledge on Hamiltonian systems, so I've tried to formulate the Kepler problem in Hamiltonian dynamics. I[/PLAIN] [Broken] got results similar to the one from this thread:

    In polar coordinates (##\varphi##, r), the Lagrangian is ##L = \frac12 m(r^2\dot \varphi^2 + \dot r^2) + \frac{\mu}r##, where μ is the standard gravitational parameter, and m the mass of the object in orbit. The generalized momenta are ##p_\phi = \frac{\partial L}{\partial \dot \varphi} = mr^2 \dot \varphi## and ##p_r = m\dot r##. The associated Hamiltonian reads

    ##H = \frac12 \frac{p_\phi^2}{m r^2} + \frac12 \frac{p_r^2}{m} - \frac{\mu}r##

    The equations of motion should hence be

    ##\frac{dr}{dt} = \frac{\partial H}{\partial p_r} = \frac{p_r}m##
    ##\frac{d\varphi}{dt} = \frac{p_\phi}{m r^2}##
    ##\frac{dp_r}{dt} = -\frac{\partial H}{\partial r} = \frac{p_\phi^2}{m r^3} - \frac{\mu}{r^2}##
    ##p_\phi = const##

    At this point I thought I'd validate this by calculating the orbit of the ISS around earth. It's nearly circular, so I should be able to set ##\dot p_r = 0## and then check that the resulting equation

    ##\frac{p_\phi^2}{m r^3} - \frac{\mu}{r^2}##

    holds for the orbit parameters (which are available e.g. from Wolfram Alpha). Problem is: The result is off by a factor of the mass ##m##. This can be seen from the formula as well: I arrive at the expression

    ##m \dot \varphi^2 r^3 = \mu##

    for the orbit parameters, whereas according to Wikipedia on "Circular orbit", ##\dot \varphi^2 r^3 = \mu## would be correct.

    I am hence assuming that I made a mistake with the mass somewhere in the derivation, but I can't find it. I'd appreciate a hint! :-)
     
    Last edited by a moderator: May 8, 2017
  2. jcsd
  3. Oct 19, 2016 #2

    phyzguy

    User Avatar
    Science Advisor

    Shouldn't your Lagrangian be [itex] L = \frac12 m(r^2\dot \varphi^2 + \dot r^2) + \frac{\mu m}r [/itex] ?
     
  4. Oct 20, 2016 #3
    Right, thanks.. I was too worried about the details and forgot to double check the potential I guess.

    ##\int \frac{GMm}{r^2} dr = -\frac{GMm}{r} = -\frac{\mu m}{r}##

    of course. Thanks a lot, now everything works out :-)
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Have something to add?
Draft saved Draft deleted



Similar Discussions: Hamiltonian formulation and the Kepler problem
  1. Hamiltonian problem (Replies: 4)

  2. Kepler Problem (Replies: 1)

  3. Hamiltonian problem (Replies: 5)

Loading...