In astronomy, perturbation is the complex motion of a massive body subject to forces other than the gravitational attraction of a single other massive body. The other forces can include a third (fourth, fifth, etc.) body, resistance, as from an atmosphere, and the off-center attraction of an oblate or otherwise misshapen body.
The study of perturbations began with the first attempts to predict planetary motions in the sky. In ancient times the causes were a mystery. Newton, at the time he formulated his laws of motion and of gravitation, applied them to the first analysis of perturbations, recognizing the complex difficulties of their calculation. Many of the great mathematicians since then have given attention to the various problems involved; throughout the 18th and 19th centuries there was demand for accurate tables of the position of the Moon and planets for marine navigation.
The complex motions of gravitational perturbations can be broken down. The hypothetical motion that the body follows under the gravitational effect of one other body only is typically a conic section, and can be readily described with the methods of geometry. This is called a two-body problem, or an unperturbed Keplerian orbit. The differences between that and the actual motion of the body are perturbations due to the additional gravitational effects of the remaining body or bodies. If there is only one other significant body then the perturbed motion is a three-body problem; if there are multiple other bodies it is an n-body problem. A general analytical solution (a mathematical expression to predict the positions and motions at any future time) exists for the two-body problem; when more than two bodies are considered analytic solutions exist only for special cases. Even the two-body problem becomes insoluble if one of the bodies is irregular in shape.
Most systems that involve multiple gravitational attractions present one primary body which is dominant in its effects (for example, a star, in the case of the star and its planet, or a planet, in the case of the planet and its satellite). The gravitational effects of the other bodies can be treated as perturbations of the hypothetical unperturbed motion of the planet or satellite around its primary body.
In methods of general perturbations, general differential equations, either of motion or of change in the orbital elements, are solved analytically, usually by series expansions. The result is usually expressed in terms of algebraic and trigonometric functions of the orbital elements of the body in question and the perturbing bodies. This can be applied generally to many different sets of conditions, and is not specific to any particular set of gravitating objects. Historically, general perturbations were investigated first. The classical methods are known as variation of the elements, variation of parameters or variation of the constants of integration. In these methods, it is considered that the body is always moving in a conic section, however the conic section is constantly changing due to the perturbations. If all perturbations were to cease at any particular instant, the body would continue in this (now unchanging) conic section indefinitely; this conic is known as the osculating orbit and its orbital elements at any particular time are what are sought by the methods of general perturbations.
General perturbations takes advantage of the fact that in many problems of celestial mechanics, the two-body orbit changes rather slowly due to the perturbations; the two-body orbit is a good first approximation. General perturbations is applicable only if the perturbing forces are about one order of magnitude smaller, or less, than the gravitational force of the primary body. In the Solar System, this is usually the case; Jupiter, the second largest body, has a mass of about 1/1000 that of the Sun.
General perturbation methods are preferred for some types of problems, as the source of certain observed motions are readily found. This is not necessarily so for special perturbations; the motions would be predicted with similar accuracy, but no information on the configurations of the perturbing bodies (for instance, an orbital resonance) which caused them would be available.
In methods of special perturbations, numerical datasets, representing values for the positions, velocities and accelerative forces on the bodies of interest, are made the basis of numerical integration of the differential equations of motion. In effect, the positions and velocities are perturbed directly, and no attempt is made to calculate the curves of the orbits or the orbital elements.
Special perturbations can be applied to any problem in celestial mechanics, as it is not limited to cases where the perturbing forces are small. Once applied only to comets and minor planets, special perturbation methods are now the basis of the most accurate machine-generated planetary ephemerides of the great astronomical almanacs. Special perturbations are also used for modeling an orbit with computers.
Cowell's formulation (so named for Philip H. Cowell, who, with A.C.D. Cromellin, used a similar method to predict the return of Halley's comet) is perhaps the simplest of the special perturbation methods. In a system of mutually interacting bodies, this method mathematically solves for the Newtonian forces on body by summing the individual interactions from the other bodies:
where is the acceleration vector of body , is the gravitational constant, is the mass of body , and are the position vectors of objects and respectively, and is the distance from object to object . All vectors being referred to the barycenter of the system. This equation is resolved into components in , , and and these are integrated numerically to form the new velocity and position vectors. This process is repeated as many times as necessary. The advantage of Cowell's method is ease of application and programming. A disadvantage is that when perturbations become large in magnitude (as when an object makes a close approach to another) the errors of the method also become large. However, for many problems in celestial mechanics, this is never the case. Another disadvantage is that in systems with a dominant central body, such as the Sun, it is necessary to carry many significant digits in the arithmetic because of the large difference in the forces of the central body and the perturbing bodies, although with modern computers this is not nearly the limitation it once was.
Encke's method begins with the osculating orbit as a reference and integrates numerically to solve for the variation from the reference as a function of time. Its advantages are that perturbations are generally small in magnitude, so the integration can proceed in larger steps (with resulting lesser errors), and the method is much less affected by extreme perturbations. Its disadvantage is complexity; it cannot be used indefinitely without occasionally updating the osculating orbit and continuing from there, a process known as rectification. Encke's method is similar to the general perturbation method of variation of the elements, except the rectification is performed at discrete intervals rather than continuously.
, and the equation of motion of is simply
and are just the equations of motion of and
for the perturbed orbit and
for the unperturbed orbit,
which, in theory, could be integrated twice to find . Since the osculating orbit is easily calculated by two-body methods, and are accounted for and can be solved. In practice, the quantity in the brackets, , is the difference of two nearly equal vectors, and further manipulation is necessary to avoid the need for extra significant digits. Encke's method was more widely used before the advent of modern computers, when much orbit computation was performed on mechanical calculating machines.
In the Solar System, many of the disturbances of one planet by another are periodic, consisting of small impulses each time a planet passes another in its orbit. This causes the bodies to follow motions that are periodic or quasi-periodic – such as the Moon in its strongly perturbed orbit, which is the subject of lunar theory. This periodic nature led to the discovery of Neptune in 1846 as a result of its perturbations of the orbit of Uranus.
On-going mutual perturbations of the planets cause long-term quasi-periodic variations in their orbital elements, most apparent when two planets' orbital periods are nearly in sync. For instance, five orbits of Jupiter (59.31 years) is nearly equal to two of Saturn (58.91 years). This causes large perturbations of both, with a period of 918 years, the time required for the small difference in their positions at conjunction to make one complete circle, first discovered by Laplace. Venus currently has the orbit with the least eccentricity, i.e. it is the closest to circular, of all the planetary orbits. In 25,000 years' time, Earth will have a more circular (less eccentric) orbit than Venus. It has been shown that long-term periodic disturbances within the Solar System can become chaotic over very long time scales; under some circumstances one or more planets can cross the orbit of another, leading to collisions.
The orbits of many of the minor bodies of the Solar System, such as comets, are often heavily perturbed, particularly by the gravitational fields of the gas giants. While many of these perturbations are periodic, others are not, and these in particular may represent aspects of chaotic motion. For example, in April 1996, Jupiter's gravitational influence caused the period of Comet Hale–Bopp's orbit to decrease from 4,206 to 2,380 years, a change that will not revert on any periodic basis.
Orbit modeling is the process of creating mathematical models to simulate motion of a massive body as it moves in orbit around another massive body due to gravity. Other forces such as gravitational attraction from tertiary bodies, air resistance, solar pressure, or thrust from a propulsion system are typically modeled as secondary effects. Directly modeling an orbit can push the limits of machine precision due to the need to model small perturbations to very large orbits. Because of this, perturbation methods are often used to model the orbit in order to achieve better accuracy.Orbital perturbation
Orbital perturbation may refer to:
Perturbation (astronomy), the classical approach to the many-body problem of astronomy
Orbital perturbation analysis (spacecraft), the activity of determining why a satellite's orbit differs from the mathematical ideal orbitOrbital perturbation analysis
Orbital perturbation analysis is the activity of determining why a satellite's orbit differs from the mathematical ideal orbit. A satellite's orbit in an ideal two-body system describes a conic section, usually an ellipse. In reality, there are several factors that cause the conic section to continually change. These deviations from the ideal Kepler's orbit are called perturbations.Perturbation
Perturbation or perturb may refer to:
Perturbation theory, mathematical methods that give approximate solutions to problems that cannot be solved exactly
Perturbation (geology), changes in the nature of alluvial deposits over time
Perturbation (astronomy), alterations to an object's orbit (e.g., caused by gravitational interactions with other bodies)
Perturbation theory (quantum mechanics), a set of approximation schemes directly related to mathematical perturbation for describing a complicated quantum system in terms of a simpler one
Perturbation (biology), an alteration of the function of a biological system, induced by external or internal mechanisms
Perturbation function, mathematical function which relates the primal and dual problemsProper orbital elements
The proper orbital elements of an orbit are constants of motion of an object in space that remain practically unchanged over an astronomically long timescale. The term is usually used to describe the three quantities:
proper semimajor axis (ap),
proper eccentricity (ep), and
proper inclination (ip).The proper elements can be contrasted with the osculating Keplerian orbital elements observed at a particular time or epoch, such as the semi-major axis, eccentricity, and inclination. Those osculating elements change in a quasi-periodic and (in principle) predictable manner due to such effects as perturbations from planets or other bodies, and precession (e.g. perihelion precession). In the Solar System, such changes usually occur on timescales of thousands of years, while proper elements are meant to be practically constant over at least tens of millions of years.
For most bodies, the osculating elements are relatively close to the proper elements because precession and perturbation effects are relatively small (see diagram). For over 99% of asteroids in the asteroid belt, the differences are less than 0.02 AU (for semi-major axis a), 0.1 (for eccentricity e), and 2° (for inclination i).
Nevertheless, this difference is non-negligible for any purposes where precision is of importance. As an example, the asteroid Ceres has osculating orbital elements (at epoch November 26, 2005)
while its proper orbital elements (independent of epoch) are
A notable exception to this small-difference rule are asteroids lying in the Kirkwood gaps, which are in strong orbital resonance with Jupiter.
To obtain proper elements for an object, one usually conducts a detailed simulation of its motion over timespans of several millions of years. Such a simulation must take into account many details of celestial mechanics including perturbations by the planets. Subsequently, one extracts quantities from the simulation which remain unchanged over this long timespan; for example, the mean inclination, eccentricity, and semi-major axis. These are the proper orbital elements.
Historically, various approximate analytic calculations were made, starting with those of Kiyotsugu Hirayama in the early 20th century. Later analytic methods often included thousands of perturbing corrections for each particular object. Presently, the method of choice is to use a computer to numerically integrate the equations of celestial dynamics, and extract constants of motion directly from a numerical analysis of the predicted positions.
At present the most prominent use of proper orbital elements is in the study of asteroid families, following in the footsteps of the pioneering work of Hirayama. A Mars-crosser asteroid 132 Aethra is the lowest numbered asteroid to not have any proper orbital elements.