Monday, 2020-04-06Welcome, Guest! | RSS

1:29 AM
Dynamics of Motion in a Central Potential: Deriving Kepler's Laws

Dynamics of Motion in a Central Potential: Deriving Kepler's Laws

Conserved Quantities

The equation of motion is:

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image019.png.

Here we use the hat ^ to denote a unit vector, so http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image020.png gives the magnitude (and sign) of the force. For Kepler's problem, http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image021.png.

(Strictly speaking, we should be using the reduced mass for planetary motion, for our Solar System, that is a small correction. It can be put in at the end if needed.)

Let's see how using vector methods we can easily find constants of motion: first, angular momentum -- just act on the equation of motion with http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image022.png

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image023.png

Since http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image024.png, we have http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image025.png, which immediately integrates to

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image026.png,

a constant, the angular momentum, and note that http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image027.png so the motion will always stay in a plane, withhttp://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image028.pngperpendicular to the plane.

This establishes that motion in a purely central force obeys a conservation law: that of angular momentum.

(As we've discussed earlier in the course, conserved quantities in dynamical systems are always related to some underlying symmetry of the Hamiltonian. The conservation of angular momentum comes from the spherical symmetry of the system: the attraction depends only on distance, not angle. In quantum mechanics, the angular momentum operator is a rotation operator: the three components of the angular momentum vector are conserved, are constants of the motion, because the Hamiltonian is invariant under rotation. That is, the angular momentum operators commute with the Hamiltonian. The classical analogy is that they have zero Poisson brackets with the Hamiltonian.)

To get back to Kepler's statement of his Laws, notice that when the planet moves through an incremental distancehttp://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image029.png it "sweeps out" an area http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image030.png, so the rate of sweeping out area is http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image031.pngKepler's Second Law is just conservation of angular momentum!

Second, conservation of energy: this time, we act on the equation of motion with http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image032.png :

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image033.png

This immediately integrates to

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image034.png

Another conservation law coming from a simple integral: conservation of energy. What symmetry does thatcorrespond to? The answer is the invariance of the Hamiltonian under time: the central force is time invariant, and we're assuming there are time-dependent potential terms, (such as from another star passing close by).

Standard Calculus Derivation of Kepler's First Law

The first mathematical proof that an elliptic orbit about a focus meant an inverse-square attraction was given by Newton, using Euclidean geometry (even though he invented calculus!). The proof is notoriously difficult to follow. Bernoulli found a fairly straightforward calculus proof in polar coordinates by changing the variable tohttp://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image035.png

The first task is to express http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image036.png in polar, meaninghttp://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image037.png coordinates.

The simplest way to find the expression for acceleration is to parameterize the planar motion as a complex number: position http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image038.png, velocity http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image039.png, notice this means http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image040.pngsince thehttp://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image041.pngensures the http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image042.pngterm is in the positivehttp://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image043.pngdirection, and differentiating again gives

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image044.png

For a central force, the only acceleration is in the http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image045.pngdirection, so http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image046.pngwhich integrates to give

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image047.png

the constancy of angular momentum.

Equating the radial components,

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image048.png

which isn't ready to integrate yet, because http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image049.png varies too. But since the angular momentum http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image050.png is constant, we can eliminate http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image049.png from the equation, giving:

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image051.png

This doesn't look too promising, but Bernoulli came up with two clever tricks. The first was to change from the variable http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image045.pngto its inverse, http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image052.png. The other was to use the constancy of angular momentum to change the variablehttp://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image053.png to http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image054.png.

Putting these together:

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image055.png

so

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image056.png

Therefore

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image057.png

and similarly

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image058.png

Going from http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image045.pngto http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image059.png in the equation of motion

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image060.png

we get

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image061.png

or

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image062.png

This equation is easy to solve! The solution is

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image063.png

where http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image064.png is a constant of integration, determined by the initial conditions.

This proves that Kepler's First Law follows from the inverse-square nature of the force, because (see beginning of lecture) the equation above is exactly the standard http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image010.png equation of an ellipse of semi major axishttp://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image005.pngand eccentricity http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image001.png, with the origin at one focus:

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image065.png

Comparing the two equations, we can find the geometry of the ellipse in terms of the angular momentum, the gravitational attraction, and the initial conditions. The angular momentum is

http://galileoandeinstein.physics.virginia.edu/7010/CM_15_files/image066.png

Category: Education | Views: 318 | Added by: farrel | Tags: Dynamics of Motion in a Central Pot | Rating: 0.0/0
Total comments: 0
Name *:
Email *:
Code *: