Registered users can edit this wiki.

Atmosphere

From Kerbal Space Program Wiki
Jump to: navigation, search
The pressures for all atmospheres
Planets Moons
TinyEve.png Eve TinyKerbin.png Kerbin TinyLaythe.png Laythe
TinyDuna.png Duna TinyJool.png Jool

The atmosphere of a celestial body slows the movement of any object passing through it, a force known as atmospheric drag (or simply drag). An atmosphere also allows for aerodynamic lift. The celestial bodies with atmospheres are the planets Eve, Kerbin, Duna and Jool, as well as Laythe, a moon of Jool. Only Kerbin and Laythe have atmospheres that contain oxygen.

Atmospheric pressure diminishes exponentially with increasing altitude. An atmosphere's scale height is the distance over which atmospheric pressure changes as a factor of e, or 2.718. For example, Kerbin's atmosphere has a scale height of 5000 m, meaning the atmospheric pressure at altitude n is 2.718 times greater than the pressure at altitude n + 5000.

Atmospheres vary in temperature, though this has no bearing on gameplay.

Atmospheres allow aerobraking and easier landing. When an atmosphere contains oxygen, it allows jet engines to work. However, an atmosphere makes taking off from a planet more difficult and increases the minimum stable orbit altitude.

Contents

Drag

A Mk1-2 pod with a Mk16-XL parachute being slowed by drag in Kerbin's atmosphere.

In the game, the force of atmospheric drag (FD) is modeled as follows:[1]

F_D = 0.5\, \rho\, v^2\, d\, A

where ρ is the atmospheric density (kg/m3), v is the ship's velocity (m/s), d is the coefficient of drag (dimensionless), and A is the cross-sectional area (m2).

Note that the cross-sectional area is not actually calculated in the game. It is instead assumed that it is directly proportional to the mass, which is an unrealistic simplification made by KSP. The parameter FlightGlobals.DragMultiplier indicates that the proportionality ratio is 0.008 m2/kg, so:

A = 0.008 \cdot m

where m is the ship's mass (kg).

The atmospheric density ρ is directly proportional to atmospheric pressure (p of unit atm), which is a function of altitude, the atmosphere's pressure at altitude 0 (p0), and scale height (H):

\begin{align}
  p &= p_0 \cdot e^\frac{-altitude}{H} \\
  \rho &= 1.2230948554874 \frac{\text{kg}}{\text{m}^3 \cdot \text{atm}} \cdot p
\end{align}

where p here is in units atm, and ρ in kg/m3. The conversion factor of 1.2230948554874 kg/(m3·atm) is given by FlightGlobals.getAtmDensity(1.0), which returns the density at 1 atmosphere (sea level on Kerbin) pressure.

The coefficient of drag (d) is calculated as the mass-weighted average of the max_drag values of all parts on the ship. For most ships without deployed parachutes, d will be very near 0.2, since this is the max_drag value of the vast majority of parts. Also a group of the same part have always the same drag coefficient.

As an example, the coefficient of drag for a craft consisting simply of a Mk1-2 Command Pod (mass 4, drag 0.2) and a deployed Mk16-XL Parachute (mass 0.3, drag 500) is:

\frac{4 \cdot 0.2 + 0.3 \cdot 500}{4 + 0.3} = 35.07

Terminal velocity

The terminal velocity of an object falling through an atmosphere is the velocity at which the force of gravity is equal to the force of drag. Terminal velocity changes as a function of altitude. Given enough time, an object falling into the atmosphere will slow to terminal velocity and then remain at terminal velocity for the rest of its fall.

Terminal velocity is important because:

  1. It describes the amount of velocity which a spacecraft must burn away when it is close to the ground.
  2. It represents the speed at which a ship should be traveling upward during a fuel-optimal ascent.

The force of gravity (FG) is:

F_G = m \cdot a = m \cdot \frac{GM}{r^2}

where m is still the ship's mass, G is the gravitational constant, M is the mass of the planet, and r is the distance from the center of the planet to the falling object.

To find terminal velocity, we set FG equal to FD:

\begin{align}
  m \cdot \frac{GM}{r^2} &= 0.5 \cdot \rho \cdot v^2 \cdot d \cdot \left(0.008 \frac{\text{m}^2}{\text{kg}} \cdot m\right) \\
  \frac{GM}{r^2} &= 0.004\frac{\text{m}^2}{\text{kg}} \cdot \rho \cdot v^2 \cdot d \\
  v  &= v_T = \sqrt{\frac{250 \frac{\text{kg}}{\text{m}^2} \cdot GM}{r^2 \cdot \rho \cdot d}}
\end{align}

Assuming d is 0.2 (which is a good approximation, provided parachutes are not in use), this simplifies to:

v_T = \sqrt{\frac{1250 \frac{\text{kg}}{\text{m}^2} \cdot GM}{r^2\, \rho}}

For the Mk16 pod and parachute example pictured above, the drag coefficient is 35.07, so its terminal velocity at sea level on Kerbin (which is 600 km from Kerbin's center) is:

v_T = \sqrt{\frac{250 \frac{\text{kg}}{\text{m}^2} \cdot GM}{r^2\, \rho \cdot 35.07}}
\rho = 1.2230948554874 \frac{\text{kg}}{\text{m}^3 \cdot \text{atm}} \cdot 1 \text{atm} \cdot e^\frac{-0 \text{m}}{5000 \text{m}}
v_T = \sqrt{\frac{250 \frac{\text{kg}}{\text{m}^2} \cdot 6.674 \cdot 10^{-11} \frac{\text{m}^3}{\text{kg} \cdot \text{s}^2} \cdot 5.2915793 \cdot 10^{22} \operatorname{kg}}{(600000 \operatorname{m})^2 \cdot 1.2230948554874 \frac{\text{kg}}{\text{m}^3} \cdot 35.07}} = 7.56 \frac{\text{m}}{\text{s}}

Examples

Altitude (m) vT (m/s)
Eve Kerbin Duna Jool Laythe
0 58.385 100.13 122.63 23.124 100.13
100 58.795 101.12 124.65 23.239 101.37
1000 62.618 110.48 144.42 24.305 113.23
10000 117.58 267.72 629.60 38.061 342.63

On-rails physics

If a ship is "on rails" (meaning it's further than 2.25 km from the actively-controlled ship) and its orbit passes through a planet's atmosphere, one of two things will happen based on atmospheric pressure at the ship's altitude:

  • below 0.01 atm: no atmospheric drag will occur — the ship will be completely unaffected
  • 0.01 atm or above: the ship will disappear

The following table gives the altitude of this 0.01 atm threshold for each celestial body with an atmosphere:

Body Altitude (m)
Eve 43 502
Kerbin 23 026
Duna 8 987
Jool 73 132
Laythe 17 528

Atmospheric height

The atmospheric height depend on the scale height of the celestial body and is where 0.000001th (0.0001%) of the surface pressure is remaining so the atmospheric pressure at the border isn't constant. Technically a craft in Jool's orbit can get lower into the atmosphere (or the atmosphere starts from a higher pressure).

alt_{\text{atmospheric height}} = -ln\left(10^{-6}\right) \cdot \text{scale height}
p_{\text{atmospheric height}} = p_0 \cdot 10^{-6}

Kerbin's atmosphere ends at 0.000001 atm and to calculate where the other celestial bodies should have the atmospheric height:

alt_{\text{atmospheric height (real)}} = -ln\left(\frac{10^{-6}}{p_0}\right) \cdot \text{scale height}

See also

Notes

  1. http://forum.kerbalspaceprogram.com/showthread.php/5235-Atmospheric-drag?p=88804&viewfull=1#post88804