Orbital mechanics

From Wikipedia, the free encyclopedia

https://en.wikipedia.org/wiki/Orbital_mechanics 

 

A satellite orbiting the earth has a tangential velocity and an inward acceleration.

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 law of universal gravitation. Orbital mechanics 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 both spacecraft and natural astronomical bodies such as star systems, planets, moons, and comets. Orbital mechanics focuses on spacecraft trajectories, including orbital maneuvers, orbital plane changes, and interplanetary transfers, and is used by mission planners to predict the results of propulsive maneuvers.

General relativity is a more exact theory than Newton's laws for calculating orbits, and is sometimes necessary for greater accuracy or in high-gravity situations (e.g. orbits near the Sun).

History

Until the rise of space travel in the twentieth century, there was little distinction between orbital and celestial mechanics. At the time of Sputnik, the field was termed 'space dynamics'. The fundamental techniques, such as those used to solve the Keplerian problem (determining position as a function of time), are therefore the same in both fields. Furthermore, the history of the fields are almost entirely shared.

Johannes Kepler was the first to successfully model planetary orbits to a high degree of accuracy, publishing his laws in 1605. Isaac Newton published more general laws of celestial motion in the first edition of Philosophiæ Naturalis Principia Mathematica (1687), which gave a method for finding the orbit of a body following a parabolic path from three observations. This was used by Edmund Halley to establish the orbits of various comets, including that which bears his name. Newton's method of successive approximation was formalised into an analytic method by Euler in 1744, whose work was in turn generalised to elliptical and hyperbolic orbits by Lambert in 1761–1777.

Another milestone in orbit determination was Carl Friedrich Gauss's assistance in the "recovery" of the dwarf planet Ceres in 1801. Gauss's method was able to use just three observations (in the form of pairs of right ascension and declination), to find the six orbital elements that completely describe an orbit. The theory of orbit determination has subsequently been developed to the point where today it is applied in GPS receivers as well as the tracking and cataloguing of newly observed minor planets. Modern orbit determination and prediction are used for operating all types of satellites and space probes, as it is necessary to know their future positions to a high degree of accuracy.

Astrodynamics was developed by astronomer Samuel Herrick beginning in the 1930s. He consulted the rocket scientist Robert Goddard and was encouraged to continue his work on space navigation techniques as Goddard believed they would be needed in the future. Astrodynamics numerical techniques were coupled with new powerful computers in the 1960s, and humans were ready to travel to the moon and return.

Practical techniques

Further information: List of orbits

Rules of thumb

The following rules of thumb are useful for situations approximated by classical mechanics under the standard assumptions of astrodynamics outlined below the rules. The specific example discussed is of a satellite orbiting a planet, but the rules of thumb could also apply to other situations, such as orbits of small bodies around a star such as the Sun.

• Kepler's laws of planetary motion:
  • Orbits are elliptical, with the heavier body at one focus of the ellipse. A special case of this is a circular orbit (a circle is a special case of ellipse) with the planet at the center.
  • A line drawn from the planet to the satellite sweeps out equal areas in equal times no matter which portion of the orbit is measured.
  • The square of a satellite's orbital period is proportional to the cube of its average distance from the planet.
• Without applying force (such as firing a rocket engine), the period and shape of the satellite's orbit won't change.
• A satellite in a low orbit (or low part of an elliptical orbit) moves more quickly with respect to the surface of the planet than a satellite in a higher orbit (or a high part of an elliptical orbit), due to the stronger gravitational attraction closer to the planet.
• If thrust is applied at only one point in the satellite's orbit, it will return to that same point on each subsequent orbit, though the rest of its path will change. Thus one cannot move from one circular orbit to another with only one brief application of thrust.
• From a circular orbit, thrust applied in a direction opposite to the satellite's motion changes orbit to elliptical; the satellite will descend and reach the lowest orbital point (the periapse) at 180 degrees away from the firing point; then it will ascend back. Thrust applied in the direction of the satellite's motion creates an elliptical orbit with its highest point (apoapse) 180 degrees away from the firing point.

The consequences of the rules of orbital mechanics are sometimes counter-intuitive. For example, if two spacecraft are in the same circular orbit and wish to dock, unless they are very close, the trailing craft cannot simply fire its engines to go faster. This will change the shape of its orbit, causing it to gain altitude and actually slow down relative to the leading craft, missing the target. The space rendezvous before docking normally takes multiple precisely calculated engine firings in multiple orbital periods requiring hours or even days to complete.

To the degree that the standard assumptions of astrodynamics do not hold, actual trajectories will vary from those calculated. For example, simple atmospheric drag is another complicating factor for objects in low Earth orbit. These rules of thumb are decidedly inaccurate when describing two or more bodies of similar mass, such as a binary star system (see n-body problem). Celestial mechanics uses more general rules applicable to a wider variety of situations. Kepler's laws of planetary motion, which can be mathematically derived from Newton's laws, hold strictly only in describing the motion of two gravitating bodies in the absence of non-gravitational forces; they also describe parabolic and hyperbolic trajectories. In the close proximity of large objects like stars the differences between classical mechanics and general relativity also become important.

Laws of astrodynamics

See also: Laplace–Runge–Lenz vector

The fundamental laws of astrodynamics are Newton's law of universal gravitation and Newton's laws of motion, while the fundamental mathematical tool is differential calculus.

Every orbit and trajectory outside atmospheres is in principle reversible, i.e., in the space-time function the time is reversed. The velocities are reversed and the accelerations are the same, including those due to rocket bursts. Thus if a rocket burst is in the direction of the velocity, in the reversed case it is opposite to the velocity. Of course in the case of rocket bursts there is no full reversal of events, both ways the same delta-v is used and the same mass ratio applies.

Standard assumptions in astrodynamics include non-interference from outside bodies, negligible mass for one of the bodies, and negligible other forces (such as from the solar wind, atmospheric drag, etc.). More accurate calculations can be made without these simplifying assumptions, but they are more complicated. The increased accuracy often does not make enough of a difference in the calculation to be worthwhile.

Kepler's laws of planetary motion may be derived from Newton's laws, when it is assumed that the orbiting body is subject only to the gravitational force of the central attractor. When an engine thrust or propulsive force is present, Newton's laws still apply, but Kepler's laws are invalidated. When the thrust stops, the resulting orbit will be different but will once again be described by Kepler's laws. The three laws are:

1. The orbit of every planet is an ellipse with the sun at one of the foci.
2. A line joining a planet and the sun sweeps out equal areas during equal intervals of time.
3. The squares of the orbital periods of planets are directly proportional to the cubes of the semi-major axis of the orbits.

Escape velocity

Main article: Escape velocity

The formula for an escape velocity is derived as follows. The specific energy (energy per unit mass) of any space vehicle is composed of two components, the specific potential energy and the specific kinetic energy. The specific potential energy associated with a planet of mass M is given by

${\displaystyle \epsilon _{p}=-{\frac {GM}{r}}\,}$

while the specific kinetic energy of an object is given by

${\displaystyle \epsilon _{k}={\frac {v^{2}}{2}}\,}$

and so the total specific orbital energy is

${\displaystyle \epsilon =\epsilon _{k}+\epsilon _{p}={\frac {v^{2}}{2}}-{\frac {GM}{r}}\,}$

Since energy is conserved, ${\displaystyle \epsilon }$ cannot depend on the distance, ${\displaystyle r}$, from the center of the central body to the space vehicle in question, i.e. v must vary with r to keep the specific orbital energy constant. Therefore, the object can reach infinite ${\displaystyle r}$ only if this quantity is nonnegative, which implies

${\displaystyle v\geq {\sqrt {\frac {2GM}{r}}}.}$

The escape velocity from the Earth's surface is about 11 km/s, but that is insufficient to send the body an infinite distance because of the gravitational pull of the Sun. To escape the Solar System from a location at a distance from the Sun equal to the distance Sun–Earth, but not close to the Earth, requires around 42 km/s velocity, but there will be "partial credit" for the Earth's orbital velocity for spacecraft launched from Earth, if their further acceleration (due to the propulsion system) carries them in the same direction as Earth travels in its orbit.

Formulae for free orbits

Orbits are conic sections, so the formula for the distance of a body for a given angle corresponds to the formula for that curve in polar coordinates, which is:

${\displaystyle r={\frac {p}{1+e\cos \theta }}}$

${\displaystyle \mu =G(m_{1}+m_{2})\,}$

${\displaystyle p=h^{2}/\mu \,}$

${\displaystyle \mu }$ is called the gravitational parameter. ${\displaystyle m_{1}}$ and ${\displaystyle m_{2}}$ are the masses of objects 1 and 2, and ${\displaystyle h}$ is the specific angular momentum of object 2 with respect to object 1. The parameter ${\displaystyle \theta }$ is known as the true anomaly, ${\displaystyle p}$ is the semi-latus rectum, while ${\displaystyle e}$ is the orbital eccentricity, all obtainable from the various forms of the six independent orbital elements.

Circular orbits

Main article: Circular orbit

All bounded orbits where the gravity of a central body dominates are elliptical in nature. A special case of this is the circular orbit, which is an ellipse of zero eccentricity. The formula for the velocity of a body in a circular orbit at distance r from the center of gravity of mass M can be derived as follows:

Centrifugal acceleration matches the acceleration due to gravity.

So,

$${\displaystyle {\frac {v^{2}}{r}}={\frac {GM}{r^{2}}}}$$

Therefore,

${\displaystyle \ v={\sqrt {{\frac {GM}{r}}\ }}}$

where ${\displaystyle G}$ is the gravitational constant, equal to

6.673 84 × 10⁻¹¹ m³/(kg·s²)

To properly use this formula, the units must be consistent; for example, ${\displaystyle M}$ must be in kilograms, and ${\displaystyle r}$ must be in meters. The answer will be in meters per second.

The quantity ${\displaystyle GM}$ is often termed the standard gravitational parameter, which has a different value for every planet or moon in the Solar System.

Once the circular orbital velocity is known, the escape velocity is easily found by multiplying by ${\displaystyle {\sqrt {2}}}$:

${\displaystyle \ v={\sqrt {2}}{\sqrt {{\frac {GM}{r}}\ }}={\sqrt {{\frac {2GM}{r}}\ }}.}$

To escape from gravity, the kinetic energy must at least match the negative potential energy. Therefore,

$${\displaystyle {\frac {1}{2}}mv^{2}={\frac {GMm}{r}}}$$

${\displaystyle v={\sqrt {{\frac {2GM}{r}}\ }}.}$

Elliptical orbits

If ${\displaystyle 0<e<1}$, then the denominator of the equation of free orbits varies with the true anomaly ${\displaystyle \theta }$, but remains positive, never becoming zero. Therefore, the relative position vector remains bounded, having its smallest magnitude at periapsis ${\displaystyle r_{p}}$, which is given by:

${\displaystyle r_{p}={\frac {p}{1+e}}}$

The maximum value ${\displaystyle r}$ is reached when ${\displaystyle \theta =180^{\circ }}$. This point is called the apoapsis, and its radial coordinate, denoted ${\displaystyle r_{a}}$, is

${\displaystyle r_{a}={\frac {p}{1-e}}}$

Let ${\displaystyle 2a}$ be the distance measured along the apse line from periapsis ${\displaystyle P}$ to apoapsis ${\displaystyle A}$, as illustrated in the equation below:

${\displaystyle 2a=r_{p}+r_{a}}$

Substituting the equations above, we get:

${\displaystyle a={\frac {p}{1-e^{2}}}}$

a is the semimajor axis of the ellipse. Solving for ${\displaystyle p}$, and substituting the result in the conic section curve formula above, we get:

${\displaystyle r={\frac {a(1-e^{2})}{1+e\cos \theta }}}$

Orbital period

Under standard assumptions the orbital period (${\displaystyle T\,\!}$) of a body traveling along an elliptic orbit can be computed as:

${\displaystyle T=2\pi {\sqrt {a^{3} \over {\mu }}}}$

where:

• ${\displaystyle \mu \,}$ is standard gravitational parameter,
• ${\displaystyle a\,\!}$ is length of semi-major axis.

Conclusions:

• The orbital period is equal to that for a circular orbit with the orbit radius equal to the semi-major axis (${\displaystyle a\,\!}$),
• For a given semi-major axis the orbital period does not depend on the eccentricity (See also: Kepler's third law).

Velocity

Under standard assumptions the orbital speed (${\displaystyle v\,}$) of a body traveling along an elliptic orbit can be computed from the Vis-viva equation as:

${\displaystyle v={\sqrt {\mu \left({2 \over {r}}-{1 \over {a}}\right)}}}$

where:

• ${\displaystyle \mu \,}$ is the standard gravitational parameter,
• ${\displaystyle r\,}$ is the distance between the orbiting bodies.
• ${\displaystyle a\,\!}$ is the length of the semi-major axis.

The velocity equation for a hyperbolic trajectory is ${\displaystyle v={\sqrt {\mu \left({2 \over {r}}+\left\vert {1 \over {a}}\right\vert \right)}}}$.

Energy

Under standard assumptions, specific orbital energy (${\displaystyle \epsilon \,}$) of elliptic orbit is negative and the orbital energy conservation equation (the Vis-viva equation) for this orbit can take the form:

${\displaystyle {v^{2} \over {2}}-{\mu \over {r}}=-{\mu \over {2a}}=\epsilon <0}$

where:

• ${\displaystyle v\,}$ is the speed of the orbiting body,
• ${\displaystyle r\,}$ is the distance of the orbiting body from the center of mass of the central body,
• ${\displaystyle a\,}$ is the semi-major axis,
• ${\displaystyle \mu \,}$ is the standard gravitational parameter.

Conclusions:

• For a given semi-major axis the specific orbital energy is independent of the eccentricity.

Using the virial theorem we find:

• the time-average of the specific potential energy is equal to ${\displaystyle 2\epsilon }$
• the time-average of ${\displaystyle r^{-1}}$ is ${\displaystyle a^{-1}}$
• the time-average of the specific kinetic energy is equal to ${\displaystyle -\epsilon }$

Parabolic orbits

If the eccentricity equals 1, then the orbit equation becomes:

${\displaystyle r={{h^{2}} \over {\mu }}{{1} \over {1+\cos \theta }}}$

where:

• ${\displaystyle r\,}$ is the radial distance of the orbiting body from the mass center of the central body,
• ${\displaystyle h\,}$ is specific angular momentum of the orbiting body,
• ${\displaystyle \theta \,}$ is the true anomaly of the orbiting body,
• ${\displaystyle \mu \,}$ is the standard gravitational parameter.

As the true anomaly θ approaches 180°, the denominator approaches zero, so that r tends towards infinity. Hence, the energy of the trajectory for which e=1 is zero, and is given by:

${\displaystyle \epsilon ={v^{2} \over 2}-{\mu \over {r}}=0}$

where:

• ${\displaystyle v\,}$ is the speed of the orbiting body.

In other words, the speed anywhere on a parabolic path is:

${\displaystyle v={\sqrt {2\mu \over {r}}}}$

Hyperbolic orbits

If ${\displaystyle e>1}$, the orbit formula,

${\displaystyle r={{h^{2}} \over {\mu }}{{1} \over {1+e\cos \theta }}}$

describes the geometry of the hyperbolic orbit. The system consists of two symmetric curves. The orbiting body occupies one of them; the other one is its empty mathematical image. Clearly, the denominator of the equation above goes to zero when ${\displaystyle \cos \theta =-1/e}$. we denote this value of true anomaly

${\displaystyle \theta _{\infty }=\cos ^{-1}\left(-{\frac {1}{e}}\right)}$

since the radial distance approaches infinity as the true anomaly approaches ${\displaystyle \theta _{\infty }}$, known as the true anomaly of the asymptote. Observe that ${\displaystyle \theta _{\infty }}$ lies between 90° and 180°. From the trigonometric identity ${\displaystyle \sin ^{2}\theta +\cos ^{2}\theta =1}$ it follows that:

${\displaystyle \sin \theta _{\infty }={\frac {1}{e}}{\sqrt {e^{2}-1}}}$

Energy

Under standard assumptions, specific orbital energy (${\displaystyle \epsilon \,}$) of a hyperbolic trajectory is greater than zero and the orbital energy conservation equation for this kind of trajectory takes form:

${\displaystyle \epsilon ={v^{2} \over 2}-{\mu \over {r}}={\mu \over {-2a}}}$

where:

• ${\displaystyle v\,}$ is the orbital velocity of orbiting body,
• ${\displaystyle r\,}$ is the radial distance of orbiting body from central body,
• ${\displaystyle a\,}$ is the negative semi-major axis of the orbit's hyperbola,
• ${\displaystyle \mu \,}$ is standard gravitational parameter.

Hyperbolic excess velocity

See also: Characteristic energy

Under standard assumptions the body traveling along a hyperbolic trajectory will attain at ${\displaystyle r=}$ infinity an orbital velocity called hyperbolic excess velocity (${\displaystyle v_{\infty }\,\!}$) that can be computed as:

${\displaystyle v_{\infty }={\sqrt {\mu \over {-a}}}\,\!}$

where:

• ${\displaystyle \mu \,\!}$ is standard gravitational parameter,
• ${\displaystyle a\,\!}$ is the negative semi-major axis of orbit's hyperbola.

The hyperbolic excess velocity is related to the specific orbital energy or characteristic energy by

${\displaystyle 2\epsilon =C_{3}=v_{\infty }^{2}\,\!}$

Calculating trajectories

Kepler's equation

One approach to calculating orbits (mainly used historically) is to use Kepler's equation:

${\displaystyle M=E-\epsilon \cdot \sin E}$.

where M is the mean anomaly, E is the eccentric anomaly, and ${\displaystyle \epsilon }$ is the eccentricity.

With Kepler's formula, finding the time-of-flight to reach an angle (true anomaly) of ${\displaystyle \theta }$ from periapsis is broken into two steps:

1. Compute the eccentric anomaly ${\displaystyle E}$ from true anomaly ${\displaystyle \theta }$
2. Compute the time-of-flight ${\displaystyle t}$ from the eccentric anomaly ${\displaystyle E}$

Finding the eccentric anomaly at a given time (the inverse problem) is more difficult. Kepler's equation is transcendental in ${\displaystyle E}$, meaning it cannot be solved for ${\displaystyle E}$ algebraically. Kepler's equation can be solved for ${\displaystyle E}$ analytically by inversion.

A solution of Kepler's equation, valid for all real values of ${\displaystyle \textstyle \epsilon }$ is:

$${\displaystyle E={\begin{cases}\displaystyle \sum _{n=1}^{\infty }{\frac {M^{\frac {n}{3}}}{n!}}\lim _{\theta \to 0}\left({\frac {\mathrm {d} ^{\,n-1}}{\mathrm {d} \theta ^{\,n-1}}}\left[\left({\frac {\theta }{\sqrt[{3}]{\theta -\sin(\theta )}}}\right)^{n}\right]\right),&\epsilon =1\\\displaystyle \sum _{n=1}^{\infty }{\frac {M^{n}}{n!}}\lim _{\theta \to 0}\left({\frac {\mathrm {d} ^{\,n-1}}{\mathrm {d} \theta ^{\,n-1}}}\left[\left({\frac {\theta }{\theta -\epsilon \cdot \sin(\theta )}}\right)^{n}\right]\right),&\epsilon \neq 1\end{cases}}}$$

Evaluating this yields:

$${\displaystyle E={\begin{cases}\displaystyle x+{\frac {1}{60}}x^{3}+{\frac {1}{1400}}x^{5}+{\frac {1}{25200}}x^{7}+{\frac {43}{17248000}}x^{9}+{\frac {1213}{7207200000}}x^{11}+{\frac {151439}{12713500800000}}x^{13}\cdots \ |\ x=(6M)^{\frac {1}{3}},&\epsilon =1\\\\\displaystyle {\frac {1}{1-\epsilon }}M-{\frac {\epsilon }{(1-\epsilon )^{4}}}{\frac {M^{3}}{3!}}+{\frac {(9\epsilon ^{2}+\epsilon )}{(1-\epsilon )^{7}}}{\frac {M^{5}}{5!}}-{\frac {(225\epsilon ^{3}+54\epsilon ^{2}+\epsilon )}{(1-\epsilon )^{10}}}{\frac {M^{7}}{7!}}+{\frac {(11025\epsilon ^{4}+4131\epsilon ^{3}+243\epsilon ^{2}+\epsilon )}{(1-\epsilon )^{13}}}{\frac {M^{9}}{9!}}\cdots ,&\epsilon \neq 1\end{cases}}}$$

Alternatively, Kepler's Equation can be solved numerically. First one must guess a value of ${\displaystyle E}$ and solve for time-of-flight; then adjust ${\displaystyle E}$ as necessary to bring the computed time-of-flight closer to the desired value until the required precision is achieved. Usually, Newton's method is used to achieve relatively fast convergence.

The main difficulty with this approach is that it can take prohibitively long to converge for the extreme elliptical orbits. For near-parabolic orbits, eccentricity ${\displaystyle \epsilon }$ is nearly 1, and substituting ${\displaystyle e=1}$ into the formula for mean anomaly, ${\displaystyle E-\sin E}$, we find ourselves subtracting two nearly-equal values, and accuracy suffers. For near-circular orbits, it is hard to find the periapsis in the first place (and truly circular orbits have no periapsis at all). Furthermore, the equation was derived on the assumption of an elliptical orbit, and so it does not hold for parabolic or hyperbolic orbits. These difficulties are what led to the development of the universal variable formulation, described below.

Conic orbits

For simple procedures, such as computing the delta-v for coplanar transfer ellipses, traditional approaches are fairly effective. Others, such as time-of-flight are far more complicated, especially for near-circular and hyperbolic orbits.

The patched conic approximation

Main article: Patched conic approximation

The Hohmann transfer orbit alone is a poor approximation for interplanetary trajectories because it neglects the planets' own gravity. Planetary gravity dominates the behavior of the spacecraft in the vicinity of a planet and in most cases Hohmann severely overestimates delta-v, and produces highly inaccurate prescriptions for burn timings. A relatively simple way to get a first-order approximation of delta-v is based on the 'Patched Conic Approximation' technique. One must choose the one dominant gravitating body in each region of space through which the trajectory will pass, and to model only that body's effects in that region. For instance, on a trajectory from the Earth to Mars, one would begin by considering only the Earth's gravity until the trajectory reaches a distance where the Earth's gravity no longer dominates that of the Sun. The spacecraft would be given escape velocity to send it on its way to interplanetary space. Next, one would consider only the Sun's gravity until the trajectory reaches the neighborhood of Mars. During this stage, the transfer orbit model is appropriate. Finally, only Mars's gravity is considered during the final portion of the trajectory where Mars's gravity dominates the spacecraft's behavior. The spacecraft would approach Mars on a hyperbolic orbit, and a final retrograde burn would slow the spacecraft enough to be captured by Mars. Friedrich Zander was one of the first to apply the patched-conics approach for astrodynamics purposes, when proposing the use of intermediary bodies' gravity for interplanetary travels, in what is known today as a gravity assist.

The size of the "neighborhoods" (or spheres of influence) vary with radius ${\displaystyle r_{SOI}}$:

${\displaystyle r_{SOI}=a_{p}\left({\frac {m_{p}}{m_{s}}}\right)^{2/5}}$

where ${\displaystyle a_{p}}$ is the semimajor axis of the planet's orbit relative to the Sun; ${\displaystyle m_{p}}$ and ${\displaystyle m_{s}}$ are the masses of the planet and Sun, respectively.

This simplification is sufficient to compute rough estimates of fuel requirements, and rough time-of-flight estimates, but it is not generally accurate enough to guide a spacecraft to its destination. For that, numerical methods are required.

The universal variable formulation

To address computational shortcomings of traditional approaches for solving the 2-body problem, the universal variable formulation was developed. It works equally well for the circular, elliptical, parabolic, and hyperbolic cases, the differential equations converging well when integrated for any orbit. It also generalizes well to problems incorporating perturbation theory.

Perturbations

The universal variable formulation works well with the variation of parameters technique, except now, instead of the six Keplerian orbital elements, we use a different set of orbital elements: namely, the satellite's initial position and velocity vectors ${\displaystyle x_{0}}$ and ${\displaystyle v_{0}}$ at a given epoch ${\displaystyle t=0}$. In a two-body simulation, these elements are sufficient to compute the satellite's position and velocity at any time in the future, using the universal variable formulation. Conversely, at any moment in the satellite's orbit, we can measure its position and velocity, and then use the universal variable approach to determine what its initial position and velocity would have been at the epoch. In perfect two-body motion, these orbital elements would be invariant (just like the Keplerian elements would be).

However, perturbations cause the orbital elements to change over time. Hence, we write the position element as ${\displaystyle x_{0}(t)}$ and the velocity element as ${\displaystyle v_{0}(t)}$, indicating that they vary with time. The technique to compute the effect of perturbations becomes one of finding expressions, either exact or approximate, for the functions ${\displaystyle x_{0}(t)}$ and ${\displaystyle v_{0}(t)}$.

The following are some effects which make real orbits differ from the simple models based on a spherical earth. Most of them can be handled on short timescales (perhaps less than a few thousand orbits) by perturbation theory because they are small relative to the corresponding two-body effects.

• Equatorial bulges cause precession of the node and the perigee
• Tesseral harmonics of the gravity field introduce additional perturbations
• Lunar and solar gravity perturbations alter the orbits
• Atmospheric drag reduces the semi-major axis unless make-up thrust is used

Over very long timescales (perhaps millions of orbits), even small perturbations can dominate, and the behavior can become chaotic. On the other hand, the various perturbations can be orchestrated by clever astrodynamicists to assist with orbit maintenance tasks, such as station-keeping, ground track maintenance or adjustment, or phasing of perigee to cover selected targets at low altitude.

Orbital maneuver

Main article: Orbital maneuver

In spaceflight, an orbital maneuver is the use of propulsion systems to change the orbit of a spacecraft. For spacecraft far from Earth—for example those in orbits around the Sun—an orbital maneuver is called a deep-space maneuver (DSM).

Orbital transfer

Transfer orbits are usually elliptical orbits that allow spacecraft to move from one (usually substantially circular) orbit to another. Usually they require a burn at the start, a burn at the end, and sometimes one or more burns in the middle.

• The Hohmann transfer orbit requires a minimal delta-v.
• A bi-elliptic transfer can require less energy than the Hohmann transfer, if the ratio of orbits is 11.94 or greater, but comes at the cost of increased trip time over the Hohmann transfer.
• Faster transfers may use any orbit that intersects both the original and destination orbits, at the cost of higher delta-v.
• Using low thrust engines (such as electrical propulsion), if the initial orbit is supersynchronous to the final desired circular orbit then the optimal transfer orbit is achieved by thrusting continuously in the direction of the velocity at apogee. This method however takes much longer due to the low thrust.

For the case of orbital transfer between non-coplanar orbits, the change-of-plane thrust must be made at the point where the orbital planes intersect (the "node"). As the objective is to change the direction of the velocity vector by an angle equal to the angle between the planes, almost all of this thrust should be made when the spacecraft is at the node near the apoapse, when the magnitude of the velocity vector is at its lowest. However, a small fraction of the orbital inclination change can be made at the node near the periapse, by slightly angling the transfer orbit injection thrust in the direction of the desired inclination change. This works because the cosine of a small angle is very nearly one, resulting in the small plane change being effectively "free" despite the high velocity of the spacecraft near periapse, as the Oberth Effect due to the increased, slightly angled thrust exceeds the cost of the thrust in the orbit-normal axis.

A Hohmann transfer from a low circular orbit to a higher circular orbit

 

A bi-elliptic transfer from a low circular starting orbit (dark blue), to a higher circular orbit (red)

 

Generic two-impulse elliptical transfer between two circular orbits

 

A general transfer from a low circular orbit to a higher circular orbit

 

An optimal sequence for transferring a satellite from a supersynchronous to a geosynchronous orbit using electric propulsion

Gravity assist and the Oberth effect

In a gravity assist, a spacecraft swings by a planet and leaves in a different direction, at a different speed. This is useful to speed or slow a spacecraft instead of carrying more fuel.

This maneuver can be approximated by an elastic collision at large distances, though the flyby does not involve any physical contact. Due to Newton's Third Law (equal and opposite reaction), any momentum gained by a spacecraft must be lost by the planet, or vice versa. However, because the planet is much, much more massive than the spacecraft, the effect on the planet's orbit is negligible.

The Oberth effect can be employed, particularly during a gravity assist operation. This effect is that use of a propulsion system works better at high speeds, and hence course changes are best done when close to a gravitating body; this can multiply the effective delta-v.

Interplanetary Transport Network and fuzzy orbits

Main article: Interplanetary Transport Network

See also: Low energy transfers

It is now possible to use computers to search for routes using the nonlinearities in the gravity of the planets and moons of the Solar System. For example, it is possible to plot an orbit from high earth orbit to Mars, passing close to one of the Earth's Trojan points. Collectively referred to as the Interplanetary Transport Network, these highly perturbative, even chaotic, orbital trajectories in principle need no fuel beyond that needed to reach the Lagrange point (in practice keeping to the trajectory requires some course corrections). The biggest problem with them is they can be exceedingly slow, taking many years. In addition launch windows can be very far apart.

They have, however, been employed on projects such as Genesis. This spacecraft visited the Earth-Sun L₁ point and returned using very little propellant.

Interplanetary Transport Network

From Wikipedia, the free encyclopedia

https://en.wikipedia.org/wiki/Interplanetary_Transport_Network 

This stylized depiction of the ITN is designed to show its (often convoluted) path through the Solar System. The green ribbon represents one path from among the many that are mathematically possible along the surface of the darker green bounding tube. Locations where the ribbon changes direction abruptly represent trajectory changes at Lagrange points, while constricted areas represent locations where objects linger in temporary orbit around a point before continuing on.

The Interplanetary Transport Network (ITN) is a collection of gravitationally determined pathways through the Solar System that require very little energy for an object to follow. The ITN makes particular use of Lagrange points as locations where trajectories through space can be redirected using little or no energy. These points have the peculiar property of allowing objects to orbit around them, despite lacking an object to orbit. While it would use little energy, transport along the network would take a long time.

History

Interplanetary transfer orbits are solutions to the gravitational three-body problem, which, for the general case, does not have analytical solutions, and is addressed by numerical analysis approximations. However, a small number of exact solutions exist, most notably the five orbits referred to as "Lagrange points", which are orbital solutions for circular orbits in the case when one body is significantly more massive.

The key to discovering the Interplanetary Transport Network was the investigation of the nature of the winding paths near the Earth-Sun and Earth-Moon Lagrange points. They were first investigated by Henri Poincaré in the 1890s. He noticed that the paths leading to and from any of those points would almost always settle, for a time, on an orbit about that point. There are in fact an infinite number of paths taking one to the point and away from it, and all of which require nearly zero change in energy to reach. When plotted, they form a tube with the orbit about the Lagrange point at one end.

The derivation of these paths traces back to mathematicians Charles C. Conley and Richard P. McGehee in 1968. Hiten, Japan's first lunar probe, was moved into lunar orbit using similar insight into the nature of paths between the Earth and the Moon. Beginning in 1997, Martin Lo, Shane D. Ross, and others wrote a series of papers identifying the mathematical basis that applied the technique to the Genesis solar wind sample return, and to lunar and Jovian missions. They referred to it as an Interplanetary Superhighway (IPS).

Paths

As it turns out, it is very easy to transit from a path leading to the point to one leading back out. This makes sense, since the orbit is unstable, which implies one will eventually end up on one of the outbound paths after spending no energy at all. Edward Belbruno coined the term "weak stability boundary"^[6] or "fuzzy boundary" for this effect.

With careful calculation, one can pick which outbound path one wants. This turned out to be useful, as many of these paths lead to some interesting points in space, such as the Earth's Moon or between the Galilean moons of Jupiter. As a result, for the cost of reaching the Earth–Sun L₂ point, which is rather low energy value, one can travel to a number of very interesting points for a little or no additional fuel cost. But the trip from Earth to Mars or other distant locations would likely take thousands of years.

The transfers are so low-energy that they make travel to almost any point in the Solar System possible. On the downside, these transfers are very slow. For trips from Earth to other planets, they are not useful for manned or unmanned probes, as the trip would take many generations. Nevertheless, they have already been used to transfer spacecraft to the Earth–Sun L₁ point, a useful point for studying the Sun that was employed in a number of recent missions, including the Genesis mission, the first to return solar wind samples to Earth. The network is also relevant to understanding Solar System dynamics; Comet Shoemaker–Levy 9 followed such a trajectory on its collision path with Jupiter.

Further explanation

The ITN is based around a series of orbital paths predicted by chaos theory and the restricted three-body problem leading to and from the orbits around the Lagrange points – points in space where the gravity between various bodies balances with the centrifugal force of an object there. For any two bodies in which one body orbits around the other, such as a star/planet or planet/moon system, there are five such points, denoted L₁ through L₅. For instance, the Earth–Moon L₁ point lies on a line between the two, where gravitational forces between them exactly balance with the centrifugal force of an object placed in orbit there. These five points have particularly low delta-v requirements, and appear to be the lowest-energy transfers possible, even lower than the common Hohmann transfer orbit that has dominated orbital navigation since the start of space travel.

Although the forces balance at these points, the first three points (the ones on the line between a certain large mass, e.g. a star, and a smaller, orbiting mass, e.g. a planet) are not stable equilibrium points. If a spacecraft placed at the Earth–Moon L₁ point is given even a slight nudge away from the equilibrium point, the spacecraft's trajectory will diverge away from the L₁ point. The entire system is in motion, so the spacecraft will not actually hit the Moon, but will travel in a winding path, off into space. There is, however, a semi-stable orbit around each of these points, called a halo orbit. The orbits for two of the points, L₄ and L₅, are stable, but the halo orbits for L₁ through L₃ are stable only on the order of months.

In addition to orbits around Lagrange points, the rich dynamics that arise from the gravitational pull of more than one mass yield interesting trajectories, also known as low energy transfers. For example, the gravity environment of the Sun–Earth–Moon system allows spacecraft to travel great distances on very little fuel, albeit on an often circuitous route.

Missions

Launched in 1978, the ISEE-3 spacecraft was sent on a mission to orbit around one of the Lagrange points. The spacecraft was able to maneuver around the Earth's neighborhood using little fuel by taking advantage of the unique gravity environment. After the primary mission was completed, ISEE-3 went on to accomplish other goals, including a flight through the geomagnetic tail and a comet flyby. The mission was subsequently renamed the International Cometary Explorer (ICE).

The first low energy transfer using what would later be called the ITN was the rescue of Japan's Hiten lunar mission in 1991.

Another example of the use of the ITN was NASA's 2001–2003 Genesis mission, which orbited the Sun–Earth L₁ point for over two years collecting material, before being redirected to the L₂ Lagrange point, and finally redirected from there back to Earth.

The 2003–2006 SMART-1 of the European Space Agency used another low energy transfer from the ITN.

In a more recent example, the Chinese spacecraft Chang'e 2 used the ITN to travel from lunar orbit to the Earth-Sun L₂ point, then on to fly by the asteroid 4179 Toutatis.

Asteroids

The asteroid 39P/Oterma's path from outside Jupiters orbit, to inside, and back to outside is said to use these low energy paths.