History
Historically, the apparent motions of the planets were described by European and Arabic philosophers using the idea of celestial spheres. This model posited the existence of perfect moving spheres or rings to which the stars and planets were attached. It assumed the heavens were fixed apart from the motion of the spheres and was developed without any understanding of gravity. After the planets' motions were more accurately measured, theoretical mechanisms such as deferent and epicycles were added. Although the model was capable of reasonably accurately predicting the planets' positions in the sky, more and more epicycles were required as the measurements became more accurate, hence the model became increasingly unwieldy. Originally geocentric, it was modified by Copernicus to place the Sun at the centre to help simplify the model. The model was further challenged during the 16th century, as comets were observed traversing the spheres.Planetary orbits
Within a planetary system, planets, dwarf planets, asteroids and other minor planets, comets, and space debris orbit the system's barycenter in elliptical orbits. A comet in a parabolic or hyperbolic orbit about a barycenter is not gravitationally bound to the star and therefore is not considered part of the star's planetary system. Bodies that are gravitationally bound to one of the planets in a planetary system, either natural or artificial satellites, follow orbits about a barycenter near or within that planet. Owing to mutual gravitational perturbations, the eccentricities of the planetary orbits vary over time. Mercury, the smallest planet in the Solar System, has the most eccentric orbit. At the present epoch, Mars has the next largest eccentricity while the smallest orbital eccentricities are seen with Venus and Neptune. As two objects orbit each other, the periapsis is that point at which the two objects are closest to each other and the apoapsis is that point at which they are the farthest. (More specific terms are used for specific bodies. For example, ''perigee'' and ''apogee'' are the lowest and highest parts of an orbit around Earth, while ''perihelion'' and ''aphelion'' are the closest and farthest points of an orbit around the Sun.) In the case of planets orbiting a star, the mass of the star and all its satellites are calculated to be at a single point called the barycenter. The paths of all the star's satellites are elliptical orbits about that barycenter. Each satellite in that system will have its own elliptical orbit with the barycenter at one focal point of that ellipse. At any point along its orbit, any satellite will have a certain value of kinetic and potential energy with respect to the barycenter, and the sum of those two energies is a constant value at every point along its orbit. As a result, as a planet approaches periapsis, the planet will increase in speed as its potential energy decreases; as a planet approaches apoapsis, its velocity will decrease as its potential energy increases.Principles
There are a few common ways of understanding orbits: * A force, such as gravity, pulls an object into a curved path as it attempts to fly off in a straight line. * As the object is pulled toward the massive body, it falls toward that body. However, if it has enough tangential velocity it will not fall into the body but will instead continue to follow the curved trajectory caused by that body indefinitely. The object is then said to be orbiting the body. The velocity relationship of two moving objects with mass can thus be considered in four practical classes, with subtypes: ; No orbit: ; Suborbital trajectories: Range of interrupted elliptical paths ; Orbital trajectories (or simply, orbits): ; Open (or escape) trajectories: Orbital rockets are launched vertically at first to lift the rocket above the atmosphere (which causes frictional drag), and then slowly pitch over and finish firing the rocket engine parallel to the atmosphere to achieve orbit speed. Once in orbit, their speed keeps them in orbit above the atmosphere. If e.g., an elliptical orbit dips into dense air, the object will lose speed and re-enter (i.e. fall). Occasionally a space craft will intentionally intercept the atmosphere, in an act commonly referred to as an aerobraking maneuver.Illustration
Newton's laws of motion
Newton's law of gravitation and laws of motion for two-body problems
In most situations, relativistic effects can be neglected, and Newton's laws give a sufficiently accurate description of motion. The acceleration of a body is equal to the sum of the forces acting on it, divided by its mass, and the gravitational force acting on a body is proportional to the product of the masses of the two attracting bodies and decreases inversely with the square of the distance between them. To this Newtonian approximation, for a system of two-point masses or spherical bodies, only influenced by their mutual gravitation (called a two-body problem), their trajectories can be exactly calculated. If the heavier body is much more massive than the smaller, as in the case of a satellite or small moon orbiting a planet or for the Earth orbiting the Sun, it is accurate enough and convenient to describe the motion in terms of a coordinate system that is centered on the heavier body, and we say that the lighter body is in orbit around the heavier. For the case where the masses of two bodies are comparable, an exact Newtonian solution is still sufficient and can be had by placing the coordinate system at the center of the mass of the system.Defining gravitational potential energy
Energy is associated with gravitational fields. A stationary body far from another can do external work if it is pulled towards it, and therefore has gravitational '' potential energy''. Since work is required to separate two bodies against the pull of gravity, their gravitational potential energy increases as they are separated, and decreases as they approach one another. For point masses, the gravitational energy decreases to zero as they approach zero separation. It is convenient and conventional to assign the potential energy as having zero value when they are an infinite distance apart, and hence it has a negative value (since it decreases from zero) for smaller finite distances.Orbital energies and orbit shapes
When only two gravitational bodies interact, their orbits follow a conic section. The orbit can be open (implying the object never returns) or closed (returning). Which it is depends on the totalKepler's laws
Bodies following closed orbits repeat their paths with a certain time called the period. This motion is described by the empirical laws of Kepler, which can be mathematically derived from Newton's laws. These can be formulated as follows: # The orbit of a planet around the Sun is an ellipse, with the Sun in one of the focal points of that ellipse. his focal point is actually the barycenter of the Solar System">Sun-planet system; for simplicity, this explanation assumes the Sun's mass is infinitely larger than that planet's.] The planet's orbit lies in a plane, called the orbital plane. The point on the orbit closest to the attracting body is the periapsis. The point farthest from the attracting body is called the apoapsis. There are also specific terms for orbits about particular bodies; things orbiting the Sun have a perihelion and aphelion, things orbiting the Earth have a perigee and apogee, and things orbiting the Moon have a perilune and apolune (or periselene and aposelene respectively). An orbit around any star, not just the Sun, has aLimitations of Newton's law of gravitation
Note that while bound orbits of a point mass or a spherical body with a Newtonian gravitational field are closed ellipses, which repeat the same path exactly and indefinitely, any non-spherical or non-Newtonian effects (such as caused by the slight oblateness of the Earth, or by relativistic effects, thereby changing the gravitational field's behavior with distance) will cause the orbit's shape to depart from the closed ellipses characteristic of Newtonian two-body motion. The two-body solutions were published by Newton in Principia in 1687. In 1912, Karl Fritiof Sundman developed a converging infinite series that solves the three-body problem; however, it converges too slowly to be of much use. Except for special cases like the Lagrangian points, no method is known to solve the equations of motion for a system with four or more bodies.Approaches to many-body problems
Rather than an exact closed form solution, orbits with many bodies can be approximated with arbitrarily high accuracy. These approximations take two forms: :One form takes the pure elliptic motion as a basis and adds perturbation terms to account for the gravitational influence of multiple bodies. This is convenient for calculating the positions of astronomical bodies. The equations of motion of the moons, planets, and other bodies are known with great accuracy, and are used to generate tables for celestial navigation. Still, there are secular phenomena that have to be dealt with by post-Newtonian methods. :The differential equation form is used for scientific or mission-planning purposes. According to Newton's laws, the sum of all the forces acting on a body will equal the mass of the body times its acceleration (''F = ma''). Therefore accelerations can be expressed in terms of positions. The perturbation terms are much easier to describe in this form. Predicting subsequent positions and velocities from initial values of position and velocity corresponds to solving an initial value problem. Numerical methods calculate the positions and velocities of the objects a short time in the future, then repeat the calculation ad nauseam. However, tiny arithmetic errors from the limited accuracy of a computer's math are cumulative, which limits the accuracy of this approach. Differential simulations with large numbers of objects perform the calculations in a hierarchical pairwise fashion between centers of mass. Using this scheme, galaxies, star clusters and other large assemblages of objects have been simulated.Formulation
Newtonian analysis of orbital motion
The following derivation applies to such an elliptical orbit. We start only with the Newtonian law of gravitation stating that the gravitational acceleration towards the central body is related to the inverse of the square of the distance between them, namely : where ''F''2 is the force acting on the mass ''m''2 caused by the gravitational attraction mass ''m''1 has for ''m''2, ''G'' is the universal gravitational constant, and ''r'' is the distance between the two masses centers. From Newton's second law, the summation of the forces acting on ''m''2 related to that body's acceleration: : where ''A''2 is the acceleration of ''m''2 caused by the force of gravitational attraction ''F''2 of ''m''1 acting on ''m''2. Combining Eq. 1 and 2: : Solving for the acceleration, ''A''2: : where is the standard gravitational parameter, in this case . It is understood that the system being described is ''m''2, hence the subscripts can be dropped. We assume that the central body is massive enough that it can be considered to be stationary and we ignore the more subtle effects of general relativity. When a pendulum or an object attached to a spring swings in an ellipse, the inward acceleration/force is proportional to the distance Due to the way vectors add, the component of the force in the or in the directions are also proportionate to the respective components of the distances, . Hence, the entire analysis can be done separately in these dimensions. This results in the harmonic parabolic equations and of the ellipse. The location of the orbiting object at the current time is located in the plane using vector calculus in polar coordinates both with the standard Euclidean basis and with the polar basis with the origin coinciding with the center of force. Let be the distance between the object and the center and be the angle it has rotated. Let and be the standard Euclidean bases and let and be the radial and transverse polar basis with the first being the unit vector pointing from the central body to the current location of the orbiting object and the second being the orthogonal unit vector pointing in the direction that the orbiting object would travel if orbiting in a counter clockwise circle. Then the vector to the orbiting object is : We use and to denote the standard derivatives of how this distance and angle change over time. We take the derivative of a vector to see how it changes over time by subtracting its location at time from that at time and dividing by . The result is also a vector. Because our basis vector moves as the object orbits, we start by differentiating it. From time to , the vector keeps its beginning at the origin and rotates from angle to which moves its head a distance in the perpendicular direction giving a derivative of . : We can now find the velocity and acceleration of our orbiting object. : The coefficients of and give the accelerations in the radial and transverse directions. As said, Newton gives this first due to gravity is and the second is zero. Equation (2) can be rearranged using integration by parts. : We can multiply through by because it is not zero unless the orbiting object crashes. Then having the derivative be zero gives that the function is a constant. which is actually the theoretical proof of Kepler's second law (A line joining a planet and the Sun sweeps out equal areas during equal intervals of time). The constant of integration, ''h'', is the angular momentum per unit mass. In order to get an equation for the orbit from equation (1), we need to eliminate time. (See also Binet equation.) In polar coordinates, this would express the distance of the orbiting object from the center as a function of its angle . However, it is easier to introduce the auxiliary variable and to express as a function of . Derivatives of with respect to time may be rewritten as derivatives of with respect to angle. : : (reworking (3)) : Plugging these into (1) gives : So for the gravitational force – or, more generally, for ''any'' inverse square force law – the right hand side of the equation becomes a constant and the equation is seen to be the harmonic equation (up to a shift of origin of the dependent variable). The solution is: : where ''A'' and ''θ''0 are arbitrary constants. This resulting equation of the orbit of the object is that of an ellipse in Polar form relative to one of the focal points. This is put into a more standard form by letting be the eccentricity, which when rearranged we see: : Note that by letting be the semi-major axis and letting so the long axis of the ellipse is along the positive ''x'' coordinate we yield: : When the two-body system is under the influence of torque, the angular momentum ''h'' is not a constant. After the following calculation: : we will get the Sturm-Liouville equation of two-body system.Relativistic orbital motion
The above classical ( Newtonian) analysis of orbital mechanics assumes that the more subtle effects of general relativity, such as frame dragging and gravitational time dilation are negligible. Relativistic effects cease to be negligible when near very massive bodies (as with the precession of Mercury's orbit about the Sun), or when extreme precision is needed (as with calculations of the orbital elements and time signal references for GPS satellites.).Specification
Six parameters are required to specify a Keplerian orbit about a body. For example, the three numbers that specify the body's initial position, and the three values that specify its velocity will define a unique orbit that can be calculated forwards (or backwards) in time. However, traditionally the parameters used are slightly different. The traditionally used set of orbital elements is called the set of Keplerian elements, after Johannes Kepler and his laws. The Keplerian elements are six: * Inclination (''i'') * Longitude of the ascending node (Ω) * Argument of periapsis (ω) * Eccentricity (''e'') * Semimajor axis (''a'') * Mean anomaly at epoch (''M''0). In principle, once the orbital elements are known for a body, its position can be calculated forward and backward indefinitely in time. However, in practice, orbits are affected or perturbed, by other forces than simple gravity from an assumed point source (see the next section), and thus the orbital elements change over time. Note that, unless the eccentricity is zero, ''a'' is not the average orbital radius. The time-averaged orbital distance is given by: :Orbital planes
The analysis so far has been two dimensional; it turns out that an unperturbed orbit is two-dimensional in a plane fixed in space, and thus the extension to three dimensions requires simply rotating the two-dimensional plane into the required angle relative to the poles of the planetary body involved. The rotation to do this in three dimensions requires three numbers to uniquely determine; traditionally these are expressed as three angles.Orbital period
The orbital period is simply how long an orbiting body takes to complete one orbit.Perturbations
An orbital perturbation is when a force or impulse which is much smaller than the overall force or average impulse of the main gravitating body and which is external to the two orbiting bodies causes an acceleration, which changes the parameters of the orbit over time.Radial, prograde and transverse perturbations
A small radial impulse given to a body in orbit changes the eccentricity, but not the orbital period (to first order). A prograde or retrograde impulse (i.e. an impulse applied along the orbital motion) changes both the eccentricity and the orbital period. Notably, a prograde impulse at periapsis raises the altitude at apoapsis, and vice versa and a retrograde impulse does the opposite. A transverse impulse (out of the orbital plane) causes rotation of the orbital plane without changing the period or eccentricity. In all instances, a closed orbit will still intersect the perturbation point.Orbital decay
If an orbit is about a planetary body with a significant atmosphere, its orbit can decay because of drag. Particularly at each periapsis, the object experiences atmospheric drag, losing energy. Each time, the orbit grows less eccentric (more circular) because the object loses kinetic energy precisely when that energy is at its maximum. This is similar to the effect of slowing a pendulum at its lowest point; the highest point of the pendulum's swing becomes lower. With each successive slowing more of the orbit's path is affected by the atmosphere and the effect becomes more pronounced. Eventually, the effect becomes so great that the maximum kinetic energy is not enough to return the orbit above the limits of the atmospheric drag effect. When this happens the body will rapidly spiral down and intersect the central body. The bounds of an atmosphere vary wildly. During a solar maximum, the Earth's atmosphere causes drag up to a hundred kilometres higher than during a solar minimum. Some satellites with long conductive tethers can also experience orbital decay because of electromagnetic drag from the Earth's magnetic field. As the wire cuts the magnetic field it acts as a generator, moving electrons from one end to the other. The orbital energy is converted to heat in the wire. Orbits can be artificially influenced through the use of rocket engines which change the kinetic energy of the body at some point in its path. This is the conversion of chemical or electrical energy to kinetic energy. In this way changes in the orbit shape or orientation can be facilitated. Another method of artificially influencing an orbit is through the use of solar sails or magnetic sails. These forms of propulsion require no propellant or energy input other than that of the Sun, and so can be used indefinitely. See statite for one such proposed use. Orbital decay can occur due to tidal forces for objects below the synchronous orbit for the body they're orbiting. The gravity of the orbiting object raises tidal bulges in the primary, and since below the synchronous orbit, the orbiting object is moving faster than the body's surface the bulges lag a short angle behind it. The gravity of the bulges is slightly off of the primary-satellite axis and thus has a component along with the satellite's motion. The near bulge slows the object more than the far bulge speeds it up, and as a result, the orbit decays. Conversely, the gravity of the satellite on the bulges applies torque on the primary and speeds up its rotation. Artificial satellites are too small to have an appreciable tidal effect on the planets they orbit, but several moons in the Solar System are undergoing orbital decay by this mechanism. Mars' innermost moon Phobos is a prime example and is expected to either impact Mars' surface or break up into a ring within 50 million years. Orbits can decay via the emission of gravitational waves. This mechanism is extremely weak for most stellar objects, only becoming significant in cases where there is a combination of extreme mass and extreme acceleration, such as with black holes or neutron stars that are orbiting each other closely.Oblateness
The standard analysis of orbiting bodies assumes that all bodies consist of uniform spheres, or more generally, concentric shells each of uniform density. It can be shown that such bodies are gravitationally equivalent to point sources. However, in the real world, many bodies rotate, and this introduces oblateness and distorts the gravity field, and gives a quadrupole moment to the gravitational field which is significant at distances comparable to the radius of the body. In the general case, the gravitational potential of a rotating body such as, e.g., a planet is usually expanded in multipoles accounting for the departures of it from spherical symmetry. From the point of view of satellite dynamics, of particular relevance are the so-called even zonal harmonic coefficients, or even zonals, since they induce secular orbital perturbations which are cumulative over time spans longer than the orbital period. They do depend on the orientation of the body's symmetry axis in the space, affecting, in general, the whole orbit, with the exception of the semimajor axis.Multiple gravitating bodies
The effects of other gravitating bodies can be significant. For example, the orbit of the Moon cannot be accurately described without allowing for the action of the Sun's gravity as well as the Earth's. One approximate result is that bodies will usually have reasonably stable orbits around a heavier planet or moon, in spite of these perturbations, provided they are orbiting well within the heavier body's Hill sphere. When there are more than two gravitating bodies it is referred to as an n-body problem. Most n-body problems have no closed form solution, although some special cases have been formulated.Light radiation and stellar wind
For smaller bodies particularly, light and stellar wind can cause significant perturbations to the attitude and direction of motion of the body, and over time can be significant. Of the planetary bodies, the motion of asteroids is particularly affected over large periods when the asteroids are rotating relative to the Sun.Strange orbits
Mathematicians have discovered that it is possible in principle to have multiple bodies in non-elliptical orbits that repeat periodically, although most such orbits are not stable regarding small perturbations in mass, position, or velocity. However, some special stable cases have been identified, including a planar figure-eight orbit occupied by three moving bodies. Further studies have discovered that nonplanar orbits are also possible, including one involving 12 masses moving in 4 roughly circular, interlocking orbits topologically equivalent to the edges of a cuboctahedron. Finding such orbits naturally occurring in the universe is thought to be extremely unlikely, because of the improbability of the required conditions occurring by chance.Astrodynamics
Orbital mechanics or astrodynamics is the application of ballistics and celestial mechanics to the practical problems concerning the motion of rockets and other spacecraft. The motion of these objects is usually calculated from Newton's laws of motion and Newton's law of universal gravitation. It is a core discipline within space mission design and control. Celestial mechanics treats more broadly the orbital dynamics of systems under the influence of gravity, including spacecraft and natural astronomical bodies such as star systems,Earth orbits
* Low Earth orbit (LEO): Geocentric orbits with altitudes up to 2,000 km (0–1,240 miles). * Medium Earth orbit (MEO): Geocentric orbits ranging in altitude from 2,000 km (1,240 miles) to just below geosynchronous orbit at . Also known as an intermediate circular orbit. These are "most commonly at , or , with an orbital period of 12 hours." * Both geosynchronous orbit (GSO) and geostationary orbit (GEO) are orbits around Earth matching Earth's sidereal rotation period. All geosynchronous and geostationary orbits have a semi-major axis of . All geostationary orbits are also geosynchronous, but not all geosynchronous orbits are geostationary. A geostationary orbit stays exactly above the equator, whereas a geosynchronous orbit may swing north and south to cover more of the Earth's surface. Both complete one full orbit of Earth per sidereal day (relative to the stars, not the Sun). * High Earth orbit: Geocentric orbits above the altitude of geosynchronous orbit 35,786 km (22,240 miles).Scaling in gravity
The gravitational constant ''G'' has been calculated as: * (6.6742 ± 0.001) × 10−11 (kg/m3)−1s−2. Thus the constant has dimension density−1 time−2. This corresponds to the following properties. Scaling of distances (including sizes of bodies, while keeping the densities the same) gives similar orbits without scaling the time: if for example distances are halved, masses are divided by 8, gravitational forces by 16 and gravitational accelerations by 2. Hence velocities are halved and orbital periods and other travel times related to gravity remain the same. For example, when an object is dropped from a tower, the time it takes to fall to the ground remains the same with a scale model of the tower on a scale model of the Earth. Scaling of distances while keeping the masses the same (in the case of point masses, or by adjusting the densities) gives similar orbits; if distances are multiplied by 4, gravitational forces and accelerations are divided by 16, velocities are halved and orbital periods are multiplied by 8. When all densities are multiplied by 4, orbits are the same; gravitational forces are multiplied by 16 and accelerations by 4, velocities are doubled and orbital periods are halved. When all densities are multiplied by 4, and all sizes are halved, orbits are similar; masses are divided by 2, gravitational forces are the same, gravitational accelerations are doubled. Hence velocities are the same and orbital periods are halved. In all these cases of scaling. if densities are multiplied by 4, times are halved; if velocities are doubled, forces are multiplied by 16. These properties are illustrated in the formula (derived from the formula for the orbital period) : for an elliptical orbit with semi-major axis ''a'', of a small body around a spherical body with radius ''r'' and average density ''ρ'', where ''T'' is the orbital period. See also Kepler’s third law.Tidal locking
Some bodies are tidally locked with other bodies, meaning that one side of the celestial body is permanently facing its host object. This is the case for Earth- Moon and Pluto-Charon system.See also
* Ephemeris is a compilation of positions of naturally occurring astronomical objects as well as artificial satellites in the sky at a given time or times. * Free drift * Klemperer rosette * List of orbits * Molniya orbit * Orbit determination * Orbital spaceflight * Perifocal coordinate system * Polar orbit * Radial trajectory * Rosetta orbit * VSOP modelReferences
Further reading
* * * Discusses new algorithms for determining the orbits of both natural and artificial celestial bodies. *External links