|
|
| (One intermediate revision by one other user not shown) |
| Line 1: |
Line 1: |
| {{About|a special case of the [[two-body problem]] in [[classical mechanics]]|the problem of [[sphere packing|finding the densest packing of spheres]] in three-dimensional Euclidean space|Kepler conjecture}}
| | The author is called Irwin. Her family lives in Minnesota. My working day occupation is a librarian. One of the extremely best issues in the globe for him is to gather badges but he is struggling to discover time for it.<br><br>Review my webpage [http://Kpchurch.org/?document_srl=1359143 Kpchurch.org] |
| | |
| In [[classical mechanics]], the '''Kepler problem''' is a special case of the [[two-body problem]], in which the two bodies interact by a [[central force]] ''F'' that varies in strength as the [[inverse square law|inverse square]] of the distance ''r'' between them. The force may be either attractive or repulsive. The "problem" to be solved is to find the position or speed of the two bodies over time given their masses and initial positions and velocities. Using classical mechanics, the solution can be expressed as a [[Kepler orbit]] using six [[orbital elements]].
| |
| | |
| The Kepler problem is named after [[Johannes Kepler]], who proposed [[Kepler's laws of planetary motion]] (which are part of [[classical mechanics]] and solve the problem for the orbits of the planets) and investigated the types of forces that would result in orbits obeying those laws (called ''Kepler's inverse problem'').<ref name="goldstein_1980">{{cite book | last=Goldstein | first=H. | authorlink=Herbert Goldstein | year=1980 | title=Classical Mechanics | edition=2nd edition | publisher=Addison Wesley}}</ref> | |
| | |
| For a discussion of the Kepler problem specific to radial orbits, see: [[Radial trajectory]]. The [[Kepler problem in general relativity]] produces more accurate predictions, especially in strong gravitational fields.
| |
| | |
| ==Applications==
| |
| The Kepler problem arises in many contexts, some beyond the physics studied by Kepler himself. The Kepler problem is important in [[celestial mechanics]], since [[gravitation|Newtonian gravity]] obeys an [[inverse square law]]. Examples include a satellite moving about a planet, a planet about its sun, or two binary stars about each other. The Kepler problem is also important in the motion of two charged particles, since [[Coulomb’s law]] of [[electrostatics]] also obeys an [[inverse square law]]. Examples include the [[hydrogen]] atom, [[positronium]] and [[muonium]], which have all played important roles as model systems for testing physical theories and measuring constants of nature.{{Citation needed|date=August 2008}}
| |
| | |
| The Kepler problem and the [[simple harmonic oscillator]] problem are the two most fundamental problems in [[classical mechanics]]. They are the ''only'' two problems that have closed orbits for every possible set of initial conditions, i.e., return to their starting point with the same velocity ([[Bertrand's theorem]]). The Kepler problem has often been used to develop new methods in classical mechanics, such as [[Lagrangian mechanics]], [[Hamiltonian mechanics]], the [[Hamilton–Jacobi equation]], and [[action-angle coordinates]].{{Citation needed|date=August 2008}} The Kepler problem also conserves the [[Laplace–Runge–Lenz vector]], which has since been generalized to include other interactions. The solution of the Kepler problem allowed scientists to show that planetary motion could be explained entirely by classical mechanics and [[gravitation|Newton’s law of gravity]]; the scientific explanation of planetary motion played an important role in ushering in the [[Age of Enlightenment|Enlightenment]].
| |
| | |
| ==Mathematical definition==
| |
| | |
| The [[central force]] '''F''' that varies in strength as the [[inverse square law|inverse square]] of the distance ''r'' between them:
| |
| | |
| :<math>
| |
| \mathbf{F} = \frac{k}{r^{2}} \mathbf{\hat{r}}
| |
| </math>
| |
| | |
| where ''k'' is a constant and <math>\mathbf{\hat{r}}</math> represents the [[unit vector]] along the line between them.<ref>{{cite book | last = Arnold | first = VI | authorlink = Vladimir Arnold | year = 1989 | title = Mathematical Methods of Classical Mechanics, 2nd ed. | publisher = Springer-Verlag | location = New York | isbn = 0-387-96890-3 | page = 38}}</ref> The force may be either attractive (''k''<0) or repulsive (''k''>0). The corresponding [[scalar potential]] (the [[potential energy]] of the non-central body) is:
| |
| | |
| :<math>
| |
| V(r) = \frac{k}{r}
| |
| </math>
| |
| | |
| ==Solution of the Kepler problem==
| |
| | |
| The equation of motion for the radius <math>r</math> of a particle
| |
| of mass <math>m</math> moving in a [[central force|central potential]] <math>V(r)</math> is given by [[Euler–Lagrange equation|Lagrange's equations]]
| |
| | |
| :<math>
| |
| m\frac{d^{2}r}{dt^{2}} - mr \omega^{2} =
| |
| m\frac{d^{2}r}{dt^{2}} - \frac{L^{2}}{mr^{3}} = -\frac{dV}{dr}
| |
| </math>
| |
| | |
| :<math>\omega \equiv \frac{d\theta}{dt}</math> and the [[angular momentum]] <math>L = mr^{2}\omega</math> is conserved. For illustration, the first term on the left-hand side is zero for circular orbits, and the applied inwards force <math>\frac{dV}{dr}</math> equals the [[centripetal force|centripetal force requirement]] <math>mr \omega^{2}</math>, as expected.
| |
| | |
| If ''L'' is not zero the definition of [[angular momentum]] allows a change of independent variable from <math>t</math> to <math>\theta</math>
| |
| | |
| :<math>
| |
| \frac{d}{dt} = \frac{L}{mr^{2}} \frac{d}{d\theta}
| |
| </math>
| |
| | |
| giving the new equation of motion that is independent of time
| |
| | |
| :<math>
| |
| \frac{L}{r^{2}} \frac{d}{d\theta} \left( \frac{L}{mr^{2}} \frac{dr}{d\theta} \right)- \frac{L^{2}}{mr^{3}} = -\frac{dV}{dr}
| |
| </math>
| |
|
| |
| The expansion of the first term is
| |
| | |
| <math>\frac{L}{r^{2}} \frac{d}{d\theta} \left( \frac{L}{mr^{2}} \frac{dr}{d\theta} \right) = -\frac{{2}L^{2}}{mr^{5}} \left( \frac{dr}{d\theta} \right)^2 + \frac{L^{2}}{mr^{4}} \frac{d^{2}r}{d\theta^{2}}
| |
| </math>
| |
| | |
| This equation becomes quasilinear on making the change of variables <math>u \equiv \frac{1}{r}</math> and multiplying both sides by <math>\frac{mr^{2}}{L^{2}}</math>
| |
| | |
| :<math>
| |
| \frac{du}{d\theta} = \frac{-1}{r^{2}} \frac{dr}{d\theta}
| |
| </math>
| |
| | |
| :<math>
| |
| \frac{d^{2}u}{d^{2}\theta^{2}} = \frac{2}{r^{3}} \left( \frac{dr}{d\theta} \right)^{2} - \frac{1}{r^{2}} \frac{d^{2}r}{d\theta^{2}} </math>
| |
| | |
| After substitution and rearrangement:
| |
| | |
| :<math>
| |
| \frac{d^{2}u}{d\theta^{2}} + u = -\frac{m}{L^{2}} \frac{d}{du} V(1/u)
| |
| </math>
| |
| | |
| For an inverse-square force law such as the [[gravity|gravitational]] or [[electrostatics|electrostatic potential]], the [[potential]] can be written
| |
| | |
| :<math>
| |
| V(\mathbf{r}) = \frac{k}{r} = ku
| |
| </math>
| |
|
| |
| The orbit <math>u(\theta)</math> can be derived from the general equation
| |
| | |
| :<math>
| |
| \frac{d^{2}u}{d\theta^{2}} + u = -\frac{m}{L^{2}} \frac{d}{du} V(1/u) = -\frac{km}{L^{2}}
| |
| </math>
| |
| | |
| whose solution is the constant <math>-\frac{km}{L^{2}}</math> plus a simple sinusoid
| |
| | |
| :<math>
| |
| u \equiv \frac{1}{r} = -\frac{km}{L^{2}} \left[ 1 + e \cos \left( \theta - \theta_{0}\right) \right]
| |
| </math>
| |
| | |
| where <math>e</math> (the '''eccentricity''') and <math>\theta_{0}</math> (the '''phase offset''') are constants of integration.
| |
| | |
| This is the general formula for a [[conic section]] that has one focus at the origin; <math>e=0</math> corresponds to a [[circle]], <math>e<1</math> corresponds to an ellipse, <math>e=1</math> corresponds to a [[parabola]], and <math>e>1</math> corresponds to a [[hyperbola]]. The eccentricity <math>e</math> is related to the total [[energy]] <math>E</math> (cf. the [[Laplace–Runge–Lenz vector]])
| |
| | |
| :<math> | |
| e = \sqrt{1 + \frac{2EL^{2}}{k^{2}m^{3}}}
| |
| </math>
| |
| | |
| Comparing these formulae shows that <math>E<0</math> corresponds to an ellipse (all solutions which are [[orbit (dynamics)|closed orbits]] are ellipses), <math>E=0</math> corresponds to a [[parabola]], and <math>E>0</math> corresponds to a [[hyperbola]]. In particular, <math>E=-\frac{k^{2}m^{3}}{2L^{2}}</math> for perfectly [[circle|circular]] orbits (the central force exactly equals the [[centripetal force|centripetal force requirement]], which determines the required angular velocity for a given circular radius).
| |
| | |
| For a repulsive force (''k'' > 0) only ''e'' > 1 applies.
| |
| | |
| ==See also==
| |
| * [[Action-angle coordinates]]
| |
| * [[Bertrand's theorem]]
| |
| * [[Binet equation]]
| |
| * [[Hamilton–Jacobi equation]]
| |
| * [[Laplace–Runge–Lenz vector]]
| |
| * [[Kepler orbit]]
| |
| * [[Kepler problem in general relativity]]
| |
| * [[Kepler's equation]]
| |
| * [[Kepler's laws of planetary motion]]
| |
| | |
| ==References==
| |
| {{reflist|30em}}
| |
| | |
| [[Category:Classical mechanics]]
| |
| [[Category:Johannes Kepler]]
| |
The author is called Irwin. Her family lives in Minnesota. My working day occupation is a librarian. One of the extremely best issues in the globe for him is to gather badges but he is struggling to discover time for it.
Review my webpage Kpchurch.org