Lably
Skip to content

Orbital decay

Process that leads to gradual decrease of the distance between two orbiting bodies

  • 7 min read
  • 4 chapters

Related on Lably

Through Distance, Orbit and Planet

Orbital decay is caused by one or more mechanisms which absorb energy from the orbital motion, such as fluid friction, gravitational anomalies, or electromagnetic effects. For bodies in low Earth orbit, the most significant effect is atmospheric drag.

Due to atmospheric drag, the lowest altitude above the Earth at which an object in a circular orbit can complete at least one full revolution without propulsion is approximately 150 km (93 mi) while the lowest perigee of an elliptical revolution is approximately 90 km (56 mi).

Modeling

Simplified model

A simplified decay model for a near-circular two-body orbit about a central body (or planet) with an atmosphere, in terms of the rate of change of the orbital altitude, is given below.

dRdt=αo(R)⋅T(R)π

Where R is the distance of the spacecraft from the planet's origin, αo is the sum of all accelerations projected on the along-track direction of the spacecraft (or parallel to the spacecraft velocity vector), and T is the Keplerian period. Note that αo is often a function of R due to variations in atmospheric density in the altitude, and T is a function of R by virtue of Kepler's laws of planetary motion.

If only atmospheric drag is considered, one can approximate drag deceleration αo as a function of orbit radius R using the drag equation below:

αo=12ρ(R)v2cdAm
ρ(R) is the mass density of the atmosphere which is a function of the radius R from the origin,
v is the orbital velocity,
A is the drag reference area,
m is the mass of the satellite, and
cd is the dimensionless drag coefficient related to the satellite geometry, and accounting for skin friction and form drag (~2.2 for cube satellites).

Proof of simplified model

By the conservation of mechanical energy, the energy of the orbit is simply the sum of kinetic and gravitational potential energies, in an unperturbed two-body orbit. By substituting the vis-viva equation into the kinetic energy component, the orbital energy of a circular orbit is given by:

U=KE+GPE=−GMEm2R

Where G is the gravitational constant, ME is the mass of the central body and m is the mass of the orbiting satellite. We take the derivative of the orbital energy with respect to the radius.

dUdR=GMEm2R2

The total decelerating force, which is usually atmospheric drag for low Earth orbits, exerted on a satellite of constant mass m is given by some force F. The rate of loss of orbital energy is simply the rate at the external force does negative work on the satellite as the satellite traverses an infinitesimal circular arc-length ds, spanned by some infinitesimal angle dθ and angular rate ω.

dUdt=F⋅dsdt=F⋅R⋅dθdt=F⋅R⋅ω

The angular rate ω is also known as the Mean motion, where for a two-body circular orbit of radius R, it is expressed as:

ω=2πT=GMER3

and...

F=m⋅αo

Substituting ω into the rate of change of orbital energy above, and expressing the external drag or decay force in terms of the deceleration αo, the orbital energy rate of change with respect to time can be expressed as:

dUdt=m⋅αo⋅GMER

Having an equation for the rate of change of orbital energy with respect to both radial distance and time allows us to find the rate of change of the radial distance with respect to time as per below.

dRdt=((dUdR)−1⋅dUdt)
=2αo⋅R3GME
=αo⋅Tπ

The assumptions used in this derivation above are that the orbit stays very nearly circular throughout the decay process, so that the equations for orbital energy are more or less that of a circular orbit's case. This is often true for orbits that begin as circular, as drag forces are considered "re-circularizing", since drag magnitudes at the periapsis (lower altitude) is expectedly greater than that of the apoapsis, which has the effect of reducing the mean eccentricity.

Sign in to keep reading

3 more chapters are free with a Lably account.

More in Engineering